[Rivet-svn] r2836 - in trunk: . data/anainfo data/plotinfo data/refdata src/Analyses

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Dec 8 09:54:35 GMT 2010


Author: hoeth
Date: Wed Dec  8 09:54:35 2010
New Revision: 2836

Log:
Added ATLAS_2010_CONF_2010_046 (only two plots, not the whole analysis).
Data is read off the plots, as usual.

Added:
   trunk/data/anainfo/ATLAS_2010_CONF_2010_046.info
   trunk/data/plotinfo/ATLAS_2010_CONF_2010_046.plot
   trunk/data/refdata/ATLAS_2010_CONF_2010_046.aida
   trunk/src/Analyses/ATLAS_2010_CONF_2010_046.cc
Modified:
   trunk/ChangeLog
   trunk/data/anainfo/Makefile.am
   trunk/data/plotinfo/Makefile.am
   trunk/data/refdata/Makefile.am
   trunk/src/Analyses/Makefile.am

Modified: trunk/ChangeLog
==============================================================================
--- trunk/ChangeLog	Wed Dec  8 08:29:25 2010	(r2835)
+++ trunk/ChangeLog	Wed Dec  8 09:54:35 2010	(r2836)
@@ -1,3 +1,11 @@
+2010-12-08  Hendrik Hoeth <hendrik.hoeth at cern.ch>
+
+	* Fixes in ATLAS_2010_CONF_083, declaring it validated
+
+	* Added ATLAS_2010_CONF_046, only two plots are implemented.
+	The paper will be out soon, and we don't need the other plots
+	right now. Data is read off the plots in the note.
+
 2010-12-07  Andy Buckley  <andy at insectnation.org>
 
 	* Update the documentation to explain this latest bump to path

Added: trunk/data/anainfo/ATLAS_2010_CONF_2010_046.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/ATLAS_2010_CONF_2010_046.info	Wed Dec  8 09:54:35 2010	(r2836)
@@ -0,0 +1,23 @@
+Name: ATLAS_2010_CONF_2010_046
+Year: 2010
+Summary: Charged particles at 900 and 7000 GeV in ATLAS
+Experiment: ATLAS
+Collider: LHC 900GeV, LHC 7000GeV
+SpiresID: None
+Status: VALIDATED -- but PRELIMINARY DATA!
+Authors:
+ - Hendrik Hoeth <hendrik.hoeth at cern.ch>
+References:
+ - ATLAS-CONF-2010-046
+RunInfo:
+  pp QCD interactions at 900 or 7000 GeV including diffractive events.
+NumEvents: 1000000
+Beams: [p+, p+]
+Energies: [900, 7000]
+PtCuts: [0]
+Description:
+  Minimum bias data measured by ATLAS at 900 and 7000 GeV. Same as
+  ATLAS_2010_S8591806, but with an additional cut on the charged particle
+  multiplicity of $N_\mathrm{ch} \ge 2$. $|\eta| < 2.5$ and $p_\perp > 100$ MeV.
+  All data is corrected to the particle level. ATTENTION - Data read from plots!
+  Only two of the plots are included!

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Wed Dec  8 08:29:25 2010	(r2835)
+++ trunk/data/anainfo/Makefile.am	Wed Dec  8 09:54:35 2010	(r2836)
@@ -9,6 +9,7 @@
   ATLAS_2010_S8591806.info \
   ATLAS_2010_S8817804.info \
   ATLAS_2010_CONF_2010_031.info \
+  ATLAS_2010_CONF_2010_046.info \
   ATLAS_2010_CONF_2010_049.info \
   ATLAS_2010_CONF_2010_081.info \
   ATLAS_2010_CONF_2010_083.info \

Added: trunk/data/plotinfo/ATLAS_2010_CONF_2010_046.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/ATLAS_2010_CONF_2010_046.plot	Wed Dec  8 09:54:35 2010	(r2836)
@@ -0,0 +1,73 @@
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d02-x01-y01
+Title=Charged multiplicity vs. $\eta$, $N_\mathrm{ch}\ge 2$, $\sqrt{s}={}$900\,GeV
+XLabel=$\eta$
+YLabel=$1/N_\mathrm{ev} \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta$
+LogY=0
+YMin=1.00
+CustomLegend={\large ATLAS preliminary}
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d02-x02-y01
+Title=Charged particle density, $\sqrt{s}={}$900\,GeV
+XLabel=$N_\mathrm{ch}$
+YLabel=$1/N_\mathrm{ev} \mathrm{d}N_\mathrm{ev}/\mathrm{d}N_\mathrm{ch}$
+FullRange=1
+CustomLegend={\large ATLAS preliminary}
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d02-x03-y01
+Title=Charged multiplicity as function of $p_\perp$, $N_\mathrm{ch}\ge 2$, $\sqrt{s}={}$900\,GeV
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/N_\mathrm{ev} \; 1/(2\pi p_\perp) \; \mathrm{d}^2N_\mathrm{ch}/\mathrm{d}\eta\mathrm{d}p_\perp$ [GeV$^{-2}$]
+LogX=1
+FullRange=1
+CustomLegend={\large ATLAS preliminary}
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d02-x04-y01
+Title=Average transverse momentum as function of $N_\mathrm{ch}$, $\sqrt{s}={}$900\,GeV
+XLabel=$N_\mathrm{ch}$
+YLabel=$\langle p_\perp \rangle$ [GeV]
+YMin=0.65
+YMax=1.4
+LogY=0
+CustomLegend={\large ATLAS preliminary}
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d03-x01-y01
+Title=Charged multiplicity vs. $\eta$, $N_\mathrm{ch}\ge 2$, $\sqrt{s}={}$7000\,GeV
+XLabel=$\eta$
+YLabel=$1/N_\mathrm{ev} \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta$
+LogY=0
+YMin=1.00
+CustomLegend={\large ATLAS preliminary}
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d03-x02-y01
+Title=Charged particle density, $\sqrt{s}={}$7000\,GeV
+XLabel=$N_\mathrm{ch}$
+YLabel=$1/N_\mathrm{ev} \mathrm{d}N_\mathrm{ev}/\mathrm{d}N_\mathrm{ch}$
+FullRange=1
+LegendXPos=0.37
+CustomLegend={\large ATLAS preliminary}
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d03-x03-y01
+Title=Charged multiplicity as function of $p_\perp$, $N_\mathrm{ch}\ge 2$, $\sqrt{s}={}$7000\,GeV
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/N_\mathrm{ev} \; 1/(2\pi p_\perp) \; \mathrm{d}^2N_\mathrm{ch}/\mathrm{d}\eta\mathrm{d}p_\perp$ [GeV$^{-2}$]
+LogX=1
+FullRange=1
+LegendXPos=0.37
+CustomLegend={\large ATLAS preliminary}
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2010_CONF_2010_046/d03-x04-y01
+Title=Average transverse momentum as function of $N_\mathrm{ch}$, $\sqrt{s}={}$7000\,GeV
+XLabel=$N_\mathrm{ch}$
+YLabel=$\langle p_\perp \rangle$ [GeV]
+YMin=0.2
+YMax=0.8
+LogY=0
+CustomLegend={\large ATLAS preliminary}
+# END PLOT

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Wed Dec  8 08:29:25 2010	(r2835)
+++ trunk/data/plotinfo/Makefile.am	Wed Dec  8 09:54:35 2010	(r2836)
@@ -9,6 +9,7 @@
   ATLAS_2010_S8591806.plot \
   ATLAS_2010_S8817804.plot \
   ATLAS_2010_CONF_2010_031.plot \
+  ATLAS_2010_CONF_2010_046.plot \
   ATLAS_2010_CONF_2010_049.plot \
   ATLAS_2010_CONF_2010_081.plot \
   ATLAS_2010_CONF_2010_083.plot \

Added: trunk/data/refdata/ATLAS_2010_CONF_2010_046.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/refdata/ATLAS_2010_CONF_2010_046.aida	Wed Dec  8 09:54:35 2010	(r2836)
@@ -0,0 +1,317 @@
+<?xml version="1.0" encoding="ISO-8859-1" ?>
+<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd">
+<aida version="3.3">
+  <implementation version="1.1" package="FreeHEP"/>
+  <dataPointSet name="d03-x02-y01" dimension="2"
+    path="/REF/ATLAS_2010_CONF_2010_046" title="">
+    <annotation>
+      <item key="Title" value="" sticky="true"/>
+      <item key="AidaPath" value="/REF/ATLAS_2010_CONF_2010_046/d03-x02-y01" sticky="true"/>
+      <item key="FullPath" value="/ATLAS_2010_CONF_2010_046.aida/REF/ATLAS_2010_CONF_2010_046/d03-x02-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="7.802000e-03" value="2.141000e-02" errorMinus="7.890000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="3.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.543000e-03" value="2.206000e-02" errorMinus="4.529000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="4.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="3.580000e-03" value="2.454000e-02" errorMinus="3.486000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="5.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="2.999000e-03" value="2.672000e-02" errorMinus="2.900000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="6.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="2.552000e-03" value="2.872000e-02" errorMinus="2.567000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="7.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="2.396000e-03" value="3.010000e-02" errorMinus="2.338000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="8.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="2.194000e-03" value="3.114000e-02" errorMinus="2.050000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="9.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.810000e-03" value="3.181000e-02" errorMinus="1.967000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.000000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.667000e-03" value="3.181000e-02" errorMinus="1.584000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.100000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.360000e-03" value="3.127000e-02" errorMinus="1.303000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.200000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="9.140000e-04" value="3.023000e-02" errorMinus="1.012000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.300000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="7.460000e-04" value="2.884000e-02" errorMinus="7.272000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.400000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.666000e-04" value="2.718000e-02" errorMinus="5.722000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.500000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="3.262000e-04" value="2.539000e-02" errorMinus="4.285000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.600000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="3.047000e-04" value="2.372000e-02" errorMinus="3.009000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.700000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.902000e-04" value="2.225000e-02" errorMinus="1.886000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.800000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.784000e-04" value="2.087000e-02" errorMinus="1.769000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.900000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.688000e-04" value="1.975000e-02" errorMinus="1.674000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.000000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.604000e-04" value="1.877000e-02" errorMinus="7.970000e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.100000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.531000e-04" value="1.791000e-02" errorMinus="1.518000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.200000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="7.319000e-05" value="1.716000e-02" errorMinus="1.455000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.300000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.406000e-04" value="1.645000e-02" errorMinus="1.394000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.400000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.359000e-04" value="1.590000e-02" errorMinus="1.347000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.500000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.308000e-04" value="1.530000e-02" errorMinus="1.297000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.600000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.259000e-04" value="1.472000e-02" errorMinus="6.253000e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.700000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.216000e-04" value="1.423000e-02" errorMinus="1.206000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.800000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.176000e-04" value="1.376000e-02" errorMinus="1.166000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.900000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.141000e-04" value="1.335000e-02" errorMinus="1.132000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="3.000000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.103000e-04" value="1.290000e-02" errorMinus="1.094000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="3.100000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.596000e-04" value="1.242000e-02" errorMinus="1.053000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e+00" value="3.300000e+01" errorMinus="1.500000e+00"/>
+      <measurement errorPlus="1.504000e-04" value="1.170000e-02" errorMinus="1.484000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e+00" value="3.700000e+01" errorMinus="2.500000e+00"/>
+      <measurement errorPlus="1.768000e-04" value="1.030000e-02" errorMinus="1.738000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="4.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="2.464000e-04" value="8.149000e-03" errorMinus="2.728000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="5.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="3.340000e-04" value="5.872000e-03" errorMinus="3.397000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="6.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="3.618000e-04" value="4.072000e-03" errorMinus="3.640000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="7.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="3.338000e-04" value="2.741000e-03" errorMinus="3.388000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e+01" value="8.950000e+01" errorMinus="1.000000e+01"/>
+      <measurement errorPlus="2.603000e-04" value="1.442000e-03" errorMinus="2.563000e-04"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d03-x04-y01" dimension="2"
+    path="/REF/ATLAS_2010_CONF_2010_046" title="">
+    <annotation>
+      <item key="Title" value="" sticky="true"/>
+      <item key="AidaPath" value="/REF/ATLAS_2010_CONF_2010_046/d03-x04-y01" sticky="true"/>
+      <item key="FullPath" value="/ATLAS_2010_CONF_2010_046.aida/REF/ATLAS_2010_CONF_2010_046/d03-x04-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.426000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="3.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.420000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="4.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.415000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="5.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.415000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="6.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.426000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="7.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.437000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="8.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.459000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="9.000000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.487000e-01" errorMinus="4.989000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.000000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.434000e-03" value="4.520000e-01" errorMinus="4.989000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.100000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.989000e-03" value="4.553000e-01" errorMinus="4.434000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.200000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.989000e-03" value="4.597000e-01" errorMinus="4.989000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.300000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.989000e-03" value="4.642000e-01" errorMinus="4.989000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.400000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="4.989000e-03" value="4.692000e-01" errorMinus="5.543000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.500000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="5.543000e-03" value="4.742000e-01" errorMinus="5.543000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.600000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="4.791000e-01" errorMinus="5.543000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.700000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="5.543000e-03" value="4.847000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.800000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="4.897000e-01" errorMinus="5.543000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="1.900000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="4.947000e-01" errorMinus="5.543000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.000000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="4.997000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.100000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.046000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.200000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.091000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.300000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.135000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.400000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.174000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.500000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.652000e-03" value="5.213000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.600000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.252000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.700000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.290000e-01" errorMinus="6.652000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.800000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.324000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="2.900000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.357000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="3.000000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.390000e-01" errorMinus="6.097000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="3.100000e+01" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="6.097000e-03" value="5.423000e-01" errorMinus="6.652000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e+00" value="3.300000e+01" errorMinus="1.500000e+00"/>
+      <measurement errorPlus="6.097000e-03" value="5.484000e-01" errorMinus="6.652000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e+00" value="3.700000e+01" errorMinus="2.500000e+00"/>
+      <measurement errorPlus="6.652000e-03" value="5.595000e-01" errorMinus="6.652000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="4.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="6.652000e-03" value="5.789000e-01" errorMinus="7.206000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="5.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="6.652000e-03" value="6.017000e-01" errorMinus="7.206000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="6.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="7.206000e-03" value="6.211000e-01" errorMinus="7.206000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e+00" value="7.450000e+01" errorMinus="5.000000e+00"/>
+      <measurement errorPlus="7.760000e-03" value="6.382000e-01" errorMinus="7.206000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e+01" value="8.950000e+01" errorMinus="1.000000e+01"/>
+      <measurement errorPlus="3.104000e-02" value="6.576000e-01" errorMinus="3.104000e-02"/>
+    </dataPoint>
+  </dataPointSet>
+</aida>

Modified: trunk/data/refdata/Makefile.am
==============================================================================
--- trunk/data/refdata/Makefile.am	Wed Dec  8 08:29:25 2010	(r2835)
+++ trunk/data/refdata/Makefile.am	Wed Dec  8 09:54:35 2010	(r2836)
@@ -11,6 +11,7 @@
   ATLAS_2010_S8591806.aida       \
   ATLAS_2010_S8817804.aida       \
   ATLAS_2010_CONF_2010_031.aida  \
+  ATLAS_2010_CONF_2010_046.aida  \
   ATLAS_2010_CONF_2010_049.aida  \
   ATLAS_2010_CONF_2010_081.aida  \
   ATLAS_2010_CONF_2010_083.aida  \

Added: trunk/src/Analyses/ATLAS_2010_CONF_2010_046.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/ATLAS_2010_CONF_2010_046.cc	Wed Dec  8 09:54:35 2010	(r2836)
@@ -0,0 +1,85 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Tools/Logging.hh"
+
+namespace Rivet {
+
+
+  /// @brief ATLAS minimum bias analysis at 900/7000 GeV with nCh>=2
+  class ATLAS_2010_CONF_2010_046 : public Analysis {
+  public:
+
+    ATLAS_2010_CONF_2010_046()
+      : Analysis("ATLAS_2010_CONF_2010_046"),
+        _Nevt_after_cuts(0.0)
+    {
+      setNeedsCrossSection(false);
+    }
+
+
+    void init() {
+      ChargedFinalState cfs(-2.5, 2.5, 0.1*GeV);
+      addProjection(cfs, "CFS");
+
+      //if (sqrtS()/GeV >= 890 && sqrtS()/GeV <= 910) {
+      //  _h_dNch_deta = bookHistogram1D(2, 1, 1);
+      //  _h_dNevt_dNch = bookHistogram1D(2, 2, 1);
+      //  _h_dNch_dpT = bookHistogram1D(2, 3, 1);
+      //  _p_meanpT_Nch = bookProfile1D(2, 4, 1);
+      //}
+      if (sqrtS()/GeV >= 6990 && sqrtS()/GeV <= 7010) {
+      //  _h_dNch_deta = bookHistogram1D(3, 1, 1);
+        _h_dNevt_dNch = bookHistogram1D(3, 2, 1);
+      //  _h_dNch_dpT = bookHistogram1D(3, 3, 1);
+        _p_meanpT_Nch = bookProfile1D(3, 4, 1);
+      }
+    }
+
+
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+
+      const ChargedFinalState& charged = applyProjection<ChargedFinalState>(event, "CFS");
+      if (charged.size() < 2) {
+        vetoEvent;
+      }
+      _Nevt_after_cuts += weight;
+
+      _h_dNevt_dNch->fill(charged.size(), weight);
+      foreach (const Particle& p, charged.particles()) {
+        double pT = p.momentum().pT()/GeV;
+        //_h_dNch_deta->fill(p.momentum().eta(), weight);
+        //_h_dNch_dpT->fill(pT, weight/pT);
+        _p_meanpT_Nch->fill(charged.size(), pT, weight);
+      }
+    }
+
+
+    void finalize() {
+      //double deta = 5.0;
+      //scale(_h_dNch_deta, 1.0/_Nevt_after_cuts);
+      //scale(_h_dNch_dpT, 1.0/_Nevt_after_cuts/TWOPI/deta);
+      scale(_h_dNevt_dNch, 1.0/_Nevt_after_cuts);
+    }
+
+
+  private:
+
+    //AIDA::IHistogram1D* _h_dNch_deta;
+    //AIDA::IHistogram1D* _h_dNch_dpT;
+    AIDA::IHistogram1D* _h_dNevt_dNch;
+    AIDA::IProfile1D*   _p_meanpT_Nch;
+
+    double _Nevt_after_cuts;
+
+  };
+
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<ATLAS_2010_CONF_2010_046> plugin_ATLAS_2010_CONF_2010_046;
+
+}
+

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Wed Dec  8 08:29:25 2010	(r2835)
+++ trunk/src/Analyses/Makefile.am	Wed Dec  8 09:54:35 2010	(r2836)
@@ -54,6 +54,7 @@
     ATLAS_2010_S8591806.cc \
     ATLAS_2010_S8817804.cc \
     ATLAS_2010_CONF_2010_031.cc \
+    ATLAS_2010_CONF_2010_046.cc \
     ATLAS_2010_CONF_2010_049.cc \
     ATLAS_2010_CONF_2010_081.cc \
     ATLAS_2010_CONF_2010_083.cc


More information about the Rivet-svn mailing list