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

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Nov 4 12:13:54 GMT 2009


Author: hoeth
Date: Wed Nov  4 12:13:53 2009
New Revision: 1997

Log:
New analysis STAR_2006_S6500200 (pion and proton pT spectra in pp
collisions at 200 GeV). It is still unclear if they used a cut in
rapidity or pseudorapidity, thus the analysis is declared "UNDER
DEVELOPMENT" and "DO NOT USE".

Added:
   trunk/data/anainfo/STAR_2006_S6500200.info
   trunk/data/plotinfo/STAR_2006_S6500200.plot
   trunk/data/refdata/STAR_2006_S6500200.aida
   trunk/data/refdata/STAR_2006_S6500200.dat
   trunk/src/Analyses/STAR_2006_S6500200.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 Nov  4 12:13:17 2009	(r1996)
+++ trunk/ChangeLog	Wed Nov  4 12:13:53 2009	(r1997)
@@ -1,3 +1,10 @@
+2009-11-04  Hendrik Hoeth <hendrik.hoeth at cern.ch>
+
+	* New analysis STAR_2006_S6500200 (pion and proton pT spectra
+	in pp collisions at 200 GeV). It is still unclear if they used
+	a cut in rapidity or pseudorapidity, thus the analysis is declared
+	"UNDER DEVELOPMENT" and "DO NOT USE".
+
 2009-11-04  Andy Buckley  <andy at insectnation.org>
 
 	* Adding consistence checking on beam ID and sqrt(s) vs. those

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Wed Nov  4 12:13:17 2009	(r1996)
+++ trunk/data/anainfo/Makefile.am	Wed Nov  4 12:13:53 2009	(r1997)
@@ -67,6 +67,7 @@
   PDG_HADRON_MULTIPLICITIES.info \
   PDG_HADRON_MULTIPLICITIES_RATIOS.info \
   SFM_1984_S1178091.info \
+  STAR_2006_S6500200.info \
   STAR_2006_S6870392.info \
   STAR_2008_S7993412.info \
   STAR_2009_UE_HELEN.info \

Added: trunk/data/anainfo/STAR_2006_S6500200.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/STAR_2006_S6500200.info	Wed Nov  4 12:13:53 2009	(r1997)
@@ -0,0 +1,20 @@
+Name: STAR_2006_S6500200
+Year: 2006
+Summary: Identified hadron spectra in pp at 200 GeV
+Experiment: STAR
+Collider: RHIC pp 200 GeV
+SpiresID: 6500200
+Status: UNDER DEVELOPMENT, UNVALIDATED, DO NOT USE.
+Authors:
+ - Hendrik Hoeth <hendrik.hoeth at cern.ch>
+References:
+ - Phys. Lett. B637, 161
+ - nucl-ex/0601033
+RunInfo:
+  pp at 200 GeV
+NumEvents: 2500000
+PtCuts: [0]
+Description:
+  pT distributions of charged pions and (anti)protons in pp collisions 
+  at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC in
+  non-single-diffractive minbias events.

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Wed Nov  4 12:13:17 2009	(r1996)
+++ trunk/data/plotinfo/Makefile.am	Wed Nov  4 12:13:53 2009	(r1997)
@@ -66,6 +66,7 @@
   PDG_HADRON_MULTIPLICITIES.plot \
   PDG_HADRON_MULTIPLICITIES_RATIOS.plot \
   SFM_1984_S1178091.plot \
+  STAR_2006_S6500200.plot \
   STAR_2006_S6870392.plot \
   STAR_2008_S7993412.plot \
   STAR_2009_UE_HELEN.plot \

Added: trunk/data/plotinfo/STAR_2006_S6500200.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/STAR_2006_S6500200.plot	Wed Nov  4 12:13:53 2009	(r1997)
@@ -0,0 +1,28 @@
+# BEGIN PLOT /STAR_2006_S6500200/d01-x01-y01
+Title=$\pi^+$ transverse momentum
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/(2\pi) \, d^2N/(p_\perp dy dp_\perp)$ [GeV$^{-2}$]
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /STAR_2006_S6500200/d01-x02-y01
+Title=$\pi^-$ transverse momentum
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/(2\pi) \, d^2N/(p_\perp dy dp_\perp)$ [GeV$^{-2}$]
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /STAR_2006_S6500200/d01-x03-y01
+Title=Proton transverse momentum
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/(2\pi) \, d^2N/(p_\perp dy dp_\perp)$ [GeV$^{-2}$]
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /STAR_2006_S6500200/d01-x04-y01
+Title=Anti-proton transverse momentum
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/(2\pi) \, d^2N/(p_\perp dy dp_\perp)$ [GeV$^{-2}$]
+FullRange=1
+# END PLOT
+

Modified: trunk/data/refdata/Makefile.am
==============================================================================
--- trunk/data/refdata/Makefile.am	Wed Nov  4 12:13:17 2009	(r1996)
+++ trunk/data/refdata/Makefile.am	Wed Nov  4 12:13:53 2009	(r1997)
@@ -20,6 +20,7 @@
   H1_2007_S7223935.aida          \
   ZEUS_2001_S4815815.aida        \
   PHENIX_2003_S5538505.aida      \
+  STAR_2006_S6500200.aida        \
   STAR_2006_S6870392.aida        \
   STAR_2008_S7993412.aida        \
   STAR_2009_UE_HELEN.aida        \

Added: trunk/data/refdata/STAR_2006_S6500200.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/refdata/STAR_2006_S6500200.aida	Wed Nov  4 12:13:53 2009	(r1997)
@@ -0,0 +1,381 @@
+<?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="d01-x01-y01" dimension="2"
+    path="/REF/STAR_2006_S6500200" title="Transverse momentum spectrum for $\pi^+$">
+    <annotation>
+      <item key="Title" value="Transverse momentum spectrum for $\pi^+$" sticky="true"/>
+      <item key="AidaPath" value="/REF/STAR_2006_S6500200/d01-x01-y01" sticky="true"/>
+      <item key="FullPath" value="/STAR_2006_S6500200.aida/REF/STAR_2006_S6500200/d01-x01-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="3.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="7.850000e-02" value="9.710000e-01" errorMinus="7.850000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="4.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="4.330000e-02" value="5.320000e-01" errorMinus="4.330000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="5.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="2.560000e-02" value="3.140000e-01" errorMinus="2.560000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="6.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.440000e-02" value="1.740000e-01" errorMinus="1.440000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="7.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="9.030000e-03" value="1.080000e-01" errorMinus="9.030000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="8.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="5.470000e-03" value="6.420000e-02" errorMinus="5.470000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="9.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="3.510000e-03" value="4.030000e-02" errorMinus="3.510000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.050000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="2.160000e-03" value="2.400000e-02" errorMinus="2.160000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.150000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.460000e-03" value="1.550000e-02" errorMinus="1.460000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.300000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="7.600000e-04" value="8.190000e-03" errorMinus="7.600000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.500000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="3.890000e-04" value="3.770000e-03" errorMinus="3.890000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.700000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="2.450000e-04" value="1.840000e-03" errorMinus="2.450000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.900000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="1.780000e-04" value="8.320000e-04" errorMinus="1.780000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.070000e-04" value="3.650000e-04" errorMinus="1.070000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.330000e-05" value="1.020000e-04" errorMinus="1.330000e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="3.800000e-06" value="2.900000e-05" errorMinus="3.800000e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.330000e-06" value="1.000000e-05" errorMinus="1.330000e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="5.370000e-07" value="3.930000e-06" errorMinus="5.370000e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="2.070000e-07" value="1.400000e-06" errorMinus="2.070000e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="5.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="7.120000e-08" value="4.530000e-07" errorMinus="7.120000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="6.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="2.520000e-08" value="1.070000e-07" errorMinus="2.520000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="7.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.320000e-08" value="3.270000e-08" errorMinus="1.320000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e+00" value="9.000000e+00" errorMinus="1.000000e+00"/>
+      <measurement errorPlus="4.690000e-09" value="7.830000e-09" errorMinus="4.690000e-09"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x02-y01" dimension="2"
+    path="/REF/STAR_2006_S6500200" title="Transverse momentum spectrum for $\pi^-$">
+    <annotation>
+      <item key="Title" value="Transverse momentum spectrum for $\pi^-$" sticky="true"/>
+      <item key="AidaPath" value="/REF/STAR_2006_S6500200/d01-x02-y01" sticky="true"/>
+      <item key="FullPath" value="/STAR_2006_S6500200.aida/REF/STAR_2006_S6500200/d01-x02-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="3.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="7.850000e-02" value="9.710000e-01" errorMinus="7.850000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="4.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="4.440000e-02" value="5.470000e-01" errorMinus="4.440000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="5.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="2.520000e-02" value="3.090000e-01" errorMinus="2.520000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="6.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.520000e-02" value="1.840000e-01" errorMinus="1.520000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="7.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="8.390000e-03" value="1.000000e-01" errorMinus="8.390000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="8.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="5.420000e-03" value="6.360000e-02" errorMinus="5.420000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="9.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="3.330000e-03" value="3.800000e-02" errorMinus="3.330000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.050000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="2.200000e-03" value="2.440000e-02" errorMinus="2.200000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.150000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.470000e-03" value="1.570000e-02" errorMinus="1.470000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.300000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="8.030000e-04" value="8.700000e-03" errorMinus="8.030000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.500000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="3.790000e-04" value="3.620000e-03" errorMinus="3.790000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.700000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="2.210000e-04" value="1.690000e-03" errorMinus="2.210000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.900000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="2.060000e-04" value="1.100000e-03" errorMinus="2.060000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.010000e-04" value="3.500000e-04" errorMinus="1.010000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.320000e-05" value="1.000000e-04" errorMinus="1.320000e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="3.720000e-06" value="2.830000e-05" errorMinus="3.720000e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.200000e-06" value="9.000000e-06" errorMinus="1.200000e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="5.270000e-07" value="3.860000e-06" errorMinus="5.270000e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="2.190000e-07" value="1.500000e-06" errorMinus="2.190000e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="5.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="7.000000e-08" value="4.400000e-07" errorMinus="7.000000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="6.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="2.510000e-08" value="1.070000e-07" errorMinus="2.510000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="7.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="1.320000e-08" value="3.270000e-08" errorMinus="1.320000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e+00" value="9.000000e+00" errorMinus="1.000000e+00"/>
+      <measurement errorPlus="4.640000e-09" value="7.830000e-09" errorMinus="4.640000e-09"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x03-y01" dimension="2"
+    path="/REF/STAR_2006_S6500200" title="Transverse momentum spectrum for $p$">
+    <annotation>
+      <item key="Title" value="Transverse momentum spectrum for $p$" sticky="true"/>
+      <item key="AidaPath" value="/REF/STAR_2006_S6500200/d01-x03-y01" sticky="true"/>
+      <item key="FullPath" value="/STAR_2006_S6500200.aida/REF/STAR_2006_S6500200/d01-x03-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="4.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="7.860000e-03" value="4.510000e-02" errorMinus="7.860000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="5.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="5.430000e-03" value="3.330000e-02" errorMinus="5.430000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="6.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="4.240000e-03" value="2.600000e-02" errorMinus="4.240000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="7.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="3.380000e-03" value="2.030000e-02" errorMinus="3.380000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="8.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.970000e-03" value="1.140000e-02" errorMinus="1.970000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="9.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.660000e-03" value="8.930000e-03" errorMinus="1.660000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.050000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.140000e-03" value="6.980000e-03" errorMinus="1.140000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.150000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="7.030000e-04" value="4.680000e-03" errorMinus="7.030000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.300000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="3.900000e-04" value="2.900000e-03" errorMinus="3.900000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.500000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="2.300000e-04" value="1.340000e-03" errorMinus="2.300000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.700000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="8.330000e-05" value="6.610000e-04" errorMinus="8.330000e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.900000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="6.730000e-05" value="4.730000e-04" errorMinus="6.730000e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="2.670000e-05" value="1.840000e-04" errorMinus="2.670000e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="4.290000e-06" value="3.500000e-05" errorMinus="4.290000e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.550000e-06" value="1.210000e-05" errorMinus="1.550000e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="5.140000e-07" value="3.620000e-06" errorMinus="5.140000e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="2.090000e-07" value="1.290000e-06" errorMinus="2.090000e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="8.480000e-08" value="4.340000e-07" errorMinus="8.480000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="5.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="3.120000e-08" value="1.430000e-07" errorMinus="3.120000e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="6.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="9.190000e-09" value="2.380000e-08" errorMinus="9.190000e-09"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x04-y01" dimension="2"
+    path="/REF/STAR_2006_S6500200" title="Transverse momentum spectrum for $\bar p$">
+    <annotation>
+      <item key="Title" value="Transverse momentum spectrum for $\bar p$" sticky="true"/>
+      <item key="AidaPath" value="/REF/STAR_2006_S6500200/d01-x04-y01" sticky="true"/>
+      <item key="FullPath" value="/STAR_2006_S6500200.aida/REF/STAR_2006_S6500200/d01-x04-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="4.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="5.511769e-03" value="3.740000e-02" errorMinus="5.511769e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="5.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="4.101573e-03" value="2.840000e-02" errorMinus="4.101573e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="6.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="3.265884e-03" value="2.270000e-02" errorMinus="3.265884e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="7.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.916331e-03" value="1.320000e-02" errorMinus="1.916331e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="8.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.547676e-03" value="1.050000e-02" errorMinus="1.547676e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="9.500000e-01" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="1.077137e-03" value="7.170000e-03" errorMinus="1.077137e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.050000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="5.855194e-04" value="5.600000e-03" errorMinus="5.855194e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-02" value="1.150000e+00" errorMinus="5.000000e-02"/>
+      <measurement errorPlus="4.137596e-04" value="3.740000e-03" errorMinus="4.137596e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.300000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="2.367636e-04" value="2.310000e-03" errorMinus="2.367636e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.500000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="1.175619e-04" value="9.690000e-04" errorMinus="1.175619e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.700000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="8.178025e-05" value="5.950000e-04" errorMinus="8.178025e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.000000e-01" value="1.900000e+00" errorMinus="1.000000e-01"/>
+      <measurement errorPlus="5.578862e-05" value="3.570000e-04" errorMinus="5.578862e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="2.099393e-05" value="1.120000e-04" errorMinus="2.099393e-05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="2.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="3.554665e-06" value="2.890000e-05" errorMinus="3.554665e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.245080e-06" value="9.640000e-06" errorMinus="1.245080e-06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="3.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="6.766284e-07" value="3.060000e-06" errorMinus="6.766284e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.250000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="1.720809e-07" value="1.050000e-06" errorMinus="1.720809e-07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="2.500000e-01" value="4.750000e+00" errorMinus="2.500000e-01"/>
+      <measurement errorPlus="8.075172e-08" value="4.030000e-07" errorMinus="8.075172e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="5.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="2.520714e-08" value="1.040000e-07" errorMinus="2.520714e-08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.000000e-01" value="6.500000e+00" errorMinus="5.000000e-01"/>
+      <measurement errorPlus="7.840312e-09" value="1.910000e-08" errorMinus="7.840312e-09"/>
+    </dataPoint>
+  </dataPointSet>
+</aida>

Added: trunk/data/refdata/STAR_2006_S6500200.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/refdata/STAR_2006_S6500200.dat	Wed Nov  4 12:13:53 2009	(r1997)
@@ -0,0 +1,230 @@
+# Identified hadron spectra at large transverse momentum
+# in p+p and d+Au collisions at sqrt(s[NN])=200 GeV
+# 
+# J. Adams et al., Phys. Lett. B 637 (2006) 161.
+
+### Figure 2: Invariant yield as function of transverse momentum for |y|<0.5.
+### Errors are statistical plus point-to-point systematic.
+
+# BEGIN HISTOGRAM pT_piplus
+AidaPath=/REF/STAR_2006_S6500200/d01-x01-y01
+Title=Transverse momentum spectrum for $\pi^+$
+0.3	0.4	9.71e-01	7.85e-02
+0.4	0.5	5.32e-01	4.33e-02
+0.5	0.6	3.14e-01	2.56e-02
+0.6	0.7	1.74e-01	1.44e-02
+0.7	0.8	1.08e-01	9.03e-03
+0.8	0.9	6.42e-02	5.47e-03
+0.9	1.0	4.03e-02	3.51e-03
+1.0	1.1	2.40e-02	2.16e-03
+1.1	1.2	1.55e-02	1.46e-03
+1.2	1.4	8.19e-03	7.60e-04
+1.4	1.6	3.77e-03	3.89e-04
+1.6	1.8	1.84e-03	2.45e-04
+1.8	2.0	8.32e-04	1.78e-04
+2.0	2.5	3.65e-04	1.07e-04
+2.5	3.0	1.02e-04	1.33e-05
+3.0	3.5	2.90e-05	3.80e-06
+3.5	4.0	1.00e-05	1.33e-06
+4.0	4.5	3.93e-06	5.37e-07
+4.5	5.0	1.40e-06	2.07e-07
+5.0	6.0	4.53e-07	7.12e-08
+6.0	7.0	1.07e-07	2.52e-08
+7.0	8.0	3.27e-08	1.32e-08
+8.0	10.0	7.83e-09	4.69e-09
+# END HISTOGRAM
+
+# BEGIN HISTOGRAM pT_piminus
+AidaPath=/REF/STAR_2006_S6500200/d01-x02-y01
+Title=Transverse momentum spectrum for $\pi^-$
+0.3	0.4	9.71e-01	7.85e-02
+0.4	0.5	5.47e-01	4.44e-02
+0.5	0.6	3.09e-01	2.52e-02
+0.6	0.7	1.84e-01	1.52e-02
+0.7	0.8	1.00e-01	8.39e-03
+0.8	0.9	6.36e-02	5.42e-03
+0.9	1.0	3.80e-02	3.33e-03
+1.0	1.1	2.44e-02	2.20e-03
+1.1	1.2	1.57e-02	1.47e-03
+1.2	1.4	8.70e-03	8.03e-04
+1.4	1.6	3.62e-03	3.79e-04
+1.6	1.8	1.69e-03	2.21e-04
+1.8	2.0	1.10e-03	2.06e-04
+2.0	2.5	3.50e-04	1.01e-04
+2.5	3.0	1.00e-04	1.32e-05
+3.0	3.5	2.83e-05	3.72e-06
+3.5	4.0	9.00e-06	1.20e-06
+4.0	4.5	3.86e-06	5.27e-07
+4.5	5.0	1.50e-06	2.19e-07
+5.0	6.0	4.40e-07	7.00e-08
+6.0	7.0	1.07e-07	2.51e-08
+7.0	8.0	3.27e-08	1.32e-08
+8.0	10.0	7.83e-09	4.64e-09
+# END HISTOGRAM
+
+# BEGIN HISTOGRAM pT_proton
+AidaPath=/REF/STAR_2006_S6500200/d01-x03-y01
+Title=Transverse momentum spectrum for $p$
+0.4	0.5	4.51e-02	7.86e-03
+0.5	0.6	3.33e-02	5.43e-03
+0.6	0.7	2.60e-02	4.24e-03
+0.7	0.8	2.03e-02	3.38e-03
+0.8	0.9	1.14e-02	1.97e-03
+0.9	1.0	8.93e-03	1.66e-03
+1.0	1.1	6.98e-03	1.14e-03
+1.1	1.2	4.68e-03	7.03e-04
+1.2	1.4	2.90e-03	3.90e-04
+1.4	1.6	1.34e-03	2.30e-04
+1.6	1.8	6.61e-04	8.33e-05
+1.8	2.0	4.73e-04	6.73e-05
+2.0	2.5	1.84e-04	2.67e-05
+2.5	3.0	3.50e-05	4.29e-06
+3.0	3.5	1.21e-05	1.55e-06
+3.5	4.0	3.62e-06	5.14e-07
+4.0	4.5	1.29e-06	2.09e-07
+4.5	5.0	4.34e-07	8.48e-08
+5.0	6.0	1.43e-07	3.12e-08
+6.0	7.0	2.38e-08	9.19e-09
+# END HISTOGRAM
+
+# BEGIN HISTOGRAM pT_antiproton
+AidaPath=/REF/STAR_2006_S6500200/d01-x04-y01
+Title=Transverse momentum spectrum for $\bar p$
+0.4	0.5	3.74e-02	5.511769e-03
+0.5	0.6	2.84e-02	4.101573e-03
+0.6	0.7	2.27e-02	3.265884e-03
+0.7	0.8	1.32e-02	1.916331e-03
+0.8	0.9	1.05e-02	1.547676e-03
+0.9	1.0	7.17e-03	1.077137e-03
+1.0	1.1	5.60e-03	5.855194e-04
+1.1	1.2	3.74e-03	4.137596e-04
+1.2	1.4	2.31e-03	2.367636e-04
+1.4	1.6	9.69e-04	1.175619e-04
+1.6	1.8	5.95e-04	8.178025e-05
+1.8	2.0	3.57e-04	5.578862e-05
+2.0	2.5	1.12e-04	2.099393e-05
+2.5	3.0	2.89e-05	3.554665e-06
+3.0	3.5	9.64e-06	1.245080e-06
+3.5	4.0	3.06e-06	6.766284e-07
+4.0	4.5	1.05e-06	1.720809e-07
+4.5	5.0	4.03e-07	8.075172e-08
+5.0	6.0	1.04e-07	2.520714e-08
+6.0	7.0	1.91e-08	7.840312e-09
+# END HISTOGRAM
+
+# 
+# Figure 4: Ratio in minbias dAu and pp collisions at |y|<0.5.
+# 
+# Errors are statistical and systematic. In addition to these, there about 7% correlated error for the ratios in p+p
+# collisions and about 10% in d+Au collisions.
+# +----------------------------------------------+
+# |                       pi-/pi+                |
+# |----------------------------------------------|
+# |p[T] [GeV/c]|                      pp NSD     |
+# |------------+---------------------------------|
+# |3.50e-01    |1.000+/-0.017(stat)+/-0.114(syst)|
+# |4.50e-01    |1.028+/-0.022(stat)+/-0.117(syst)|
+# |5.50e-01    |0.984+/-0.024(stat)+/-0.111(syst)|
+# |6.50e-01    |1.057+/-0.024(stat)+/-0.111(syst)|
+# |7.50e-01    |0.926+/-0.033(stat)+/-0.104(syst)|
+# |8.50e-01    |0.991+/-0.042(stat)+/-0.111(syst)|
+# |9.5e-01     |0.943+/-0.048(stat)+/-0.105(syst)|
+# |1.05e+00    |1.017+/-0.060(stat)+/-0.116(syst)|
+# |1.15e+00    |1.013+/-0.072(stat)+/-0.115(syst)|
+# |1.30e+00    |1.063+/-0.070(stat)+/-0.121(syst)|
+# |1.50e+00    |0.960+/-0.092(stat)+/-0.108(syst)|
+# |1.70e+00    |0.918+/-0.140(stat)+/-0.106(syst)|
+# |1.90e+00    |1.322+/-0.306(stat)+/-0.223(syst)|
+# |2.25e+00    |0.959+/-0.205(stat)+/-0.339(syst)|
+# |2.75e+00    |0.976+/-0.026(stat)+/-0.108(syst)|
+# |3.25e+00    |0.976+/-0.026(stat)+/-0.107(syst)|
+# |3.75e+00    |0.900+/-0.040(stat)+/-0.099(syst)|
+# |4.25e+00    |0.981+/-0.057(stat)+/-0.108(syst)|
+# |4.75e+00    |1.071+/-0.105(stat)+/-0.118(syst)|
+# |5.50e+00    |0.970+/-0.123(stat)+/-0.107(syst)|
+# |6.50e+00    |1.00+/-0.265(stat)+/-0.110(syst) |
+# |7.50e+00    |1.00+/-0.420(stat)+/-0.110(syst) |
+# |9.00e+00    |1.00+/-0.560(stat)+/-0.110(syst) |
+# +----------------------------------------------+
+# +----------------------------------------------+
+# |                   pbar/p                     |
+# |----------------------------------------------|
+# |p[T] [GeV/c]|             pp NSD              |
+# |------------+---------------------------------|
+# |4.68e-01    |0.829+/-0.112(stat)+/-0.153(syst)|
+# |5.63e-01    |0.853+/-0.100(stat)+/-0.156(syst)|
+# |6.61e-01    |0.873+/-0.102(stat)+/-0.162(syst)|
+# |7.59e-01    |0.650+/-0.079(stat)+/-0.118(syst)|
+# |8.58e-01    |0.921+/-0.122(stat)+/-0.172(syst)|
+# |9.57e-01    |0.803+/-0.124(stat)+/-0.149(syst)|
+# |1.06e+00    |0.802+/-0.127(stat)+/-0.091(syst)|
+# |1.16e+00    |0.806+/-0.121(stat)+/-0.092(syst)|
+# |1.31e+00    |0.796+/-0.099(stat)+/-0.088(syst)|
+# |1.50e+00    |0.723+/-0.126(stat)+/-0.083(syst)|
+# |1.70e+00    |0.900+/-0.126(stat)+/-0.102(syst)|
+# |1.90e+00    |0.755+/-0.135(stat)+/-0.086(syst)|
+# |2.25e+00    |0.609+/-0.126(stat)+/-0.069(syst)|
+# |2.75e+00    |0.826+/-0.024(stat)+/-0.090(syst)|
+# |3.25e+00    |0.795+/-0.036(stat)+/-0.087(syst)|
+# |3.75e+00    |0.846+/-0.156(stat)+/-0.093(syst)|
+# |4.25e+00    |0.808+/-0.101(stat)+/-0.088(syst)|
+# |4.75e+00    |0.929+/-0.180(stat)+/-0.102(syst)|
+# |5.50e+00    |0.727+/-0.173(stat)+/-0.080(syst)|
+# |6.50e+00    |0.804+/-0.413(stat)+/-0.088(syst)|
+# +----------------------------------------------+
+# +----------------------------------------------+
+# |                   p/pi+                      |
+# |----------------------------------------------|
+# |p[T] [GeV/c]|             pp NSD              |
+# |------------+---------------------------------|
+# |4.68e-01    |0.085+/-0.009(stat)+/-0.013(syst)|
+# |5.63e-01    |0.106+/-0.011(stat)+/-0.016(syst)|
+# |6.61e-01    |0.149+/-0.015(stat)+/-0.023(syst)|
+# |7.59e-01    |0.188+/-0.020(stat)+/-0.028(syst)|
+# |8.58e-01    |0.177+/-0.021(stat)+/-0.027(syst)|
+# |9.57e-01    |0.221+/-0.031(stat)+/-0.034(syst)|
+# |1.06e+00    |0.290+/-0.043(stat)+/-0.033(syst)|
+# |1.16e+00    |0.299+/-0.041(stat)+/-0.034(syst)|
+# |1.31e+00    |0.354+/-0.041(stat)+/-0.040(syst)|
+# |1.50e+00    |0.355+/-0.058(stat)+/-0.040(syst)|
+# |1.70e+00    |0.359+/-0.052(stat)+/-0.041(syst)|
+# |1.90e+00    |0.568+/-0.122(stat)+/-0.082(syst)|
+# |2.25e+00    |0.504+/-0.099(stat)+/-0.132(syst)|
+# |2.75e+00    |0.343+/-0.008(stat)+/-0.048(syst)|
+# |3.25e+00    |0.418+/-0.015(stat)+/-0.059(syst)|
+# |3.75e+00    |0.362+/-0.023(stat)+/-0.052(syst)|
+# |4.25e+00    |0.329+/-0.031(stat)+/-0.048(syst)|
+# |4.75e+00    |0.310+/-0.047(stat)+/-0.045(syst)|
+# |5.50e+00    |0.316+/-0.056(stat)+/-0.047(syst)|
+# |6.50e+00    |0.223+/-0.088(stat)+/-0.034(syst)|
+# +----------------------------------------------+
+# +----------------------------------------------+
+# |                  pbar/pi-                    |
+# |----------------------------------------------|
+# |p[T] [GeV/c]|             pp NSD              |
+# |------------+---------------------------------|
+# |4.68e-01    |0.068+/-0.005(stat)+/-0.010(syst)|
+# |5.63e-01    |0.091+/-0.006(stat)+/-0.014(syst)|
+# |6.61e-01    |0.123+/-0.008(stat)+/-0.019(syst)|
+# |7.59e-01    |0.132+/-0.009(stat)+/-0.020(syst)|
+# |8.58e-01    |0.165+/-0.012(stat)+/-0.026(syst)|
+# |9.57e-01    |0.188+/-0.016(stat)+/-0.029(syst)|
+# |1.06e+00    |0.229+/-0.018(stat)+/-0.026(syst)|
+# |1.16e+00    |0.238+/-0.022(stat)+/-0.027(syst)|
+# |1.31e+00    |0.265+/-0.021(stat)+/-0.030(syst)|
+# |1.50e+00    |0.268+/-0.030(stat)+/-0.030(syst)|
+# |1.70e+00    |0.352+/-0.054(stat)+/-0.040(syst)|
+# |1.90e+00    |0.324+/-0.064(stat)+/-0.046(syst)|
+# |2.25e+00    |0.320+/-0.071(stat)+/-0.084(syst)|
+# |2.75e+00    |0.289+/-0.007(stat)+/-0.040(syst)|
+# |3.25e+00    |0.340+/-0.013(stat)+/-0.048(syst)|
+# |3.75e+00    |0.340+/-0.061(stat)+/-(0.049syst)|
+# |4.25e+00    |0.271+/-0.027(stat)+/-0.039(syst)|
+# |4.75e+00    |0.269+/-0.042(stat)+/-0.039(syst)|
+# |5.50e+00    |0.236+/-0.048(stat)+/-(0.035syst)|
+# |6.50e+00    |0.179+/-0.075(stat)+/-0.027(syst)|
+# |------------+---------------------------------|
+# -----------------------------------------------------------------------------------------------------------------------
+# 
+# If you have question or find any problems with this web page please send mail to James Dunlop.
+# This page was last updated on 

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Wed Nov  4 12:13:17 2009	(r1996)
+++ trunk/src/Analyses/Makefile.am	Wed Nov  4 12:13:53 2009	(r1997)
@@ -79,6 +79,7 @@
 
 lib_LTLIBRARIES += RivetRHICAnalyses.la
 RivetRHICAnalyses_la_SOURCES = \
+    STAR_2006_S6500200.cc \
     STAR_2006_S6870392.cc \
     STAR_2008_S7993412.cc \
     STAR_2009_UE_HELEN.cc

Added: trunk/src/Analyses/STAR_2006_S6500200.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/STAR_2006_S6500200.cc	Wed Nov  4 12:13:53 2009	(r1997)
@@ -0,0 +1,89 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/IdentifiedFinalState.hh"
+#include "Rivet/RivetAIDA.hh"
+
+namespace Rivet {
+
+  /// @brief identified hadron spectra in pp at 200 GeV
+  class STAR_2006_S6500200 : public Analysis {
+  public:
+
+    /// Constructor
+    STAR_2006_S6500200()
+      : Analysis("STAR_2006_S6500200")
+    {
+      setBeams(PROTON, PROTON);
+    }
+
+    /// Book projections and histograms
+    void init() {
+      IdentifiedFinalState pionfs(-3.0, 3.0, 0.3*GeV);
+      IdentifiedFinalState protonfs(-3.0, 3.0, 0.4*GeV);
+      pionfs.acceptIdPair(PIPLUS);
+      protonfs.acceptIdPair(PROTON);
+      addProjection(pionfs, "PIONFS");
+      addProjection(protonfs, "PROTONFS");
+
+      _h_pT_piplus     = bookHistogram1D(1, 1, 1);
+      _h_pT_piminus    = bookHistogram1D(1, 2, 1);
+      _h_pT_proton     = bookHistogram1D(1, 3, 1);
+      _h_pT_antiproton = bookHistogram1D(1, 4, 1);
+    }
+
+
+    /// Do the analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+
+      const IdentifiedFinalState& pionfs = applyProjection<IdentifiedFinalState>(event, "PIONFS");
+      const IdentifiedFinalState& protonfs = applyProjection<IdentifiedFinalState>(event, "PROTONFS");
+      foreach (const Particle& p, pionfs.particles()) {
+        if (fabs(p.momentum().rapidity()) < 0.5) {
+          const double pT = p.momentum().pT() / GeV;
+          if (p.pdgId()>0) {
+            _h_pT_piplus->fill(pT, weight/pT);
+          }
+          else {
+            _h_pT_piminus->fill(pT, weight/pT);
+          }
+        }
+      }
+      foreach (const Particle& p, protonfs.particles()) {
+        if (fabs(p.momentum().rapidity()) < 0.5) {
+          const double pT = p.momentum().pT() / GeV;
+          if (p.pdgId()>0) {
+            _h_pT_proton->fill(pT, weight/pT);
+          }
+          else {
+            _h_pT_antiproton->fill(pT, weight/pT);
+          }
+        }
+      }
+    }
+
+
+    /// Finalize
+    void finalize() {
+      scale(_h_pT_piplus,     1./(2*M_PI*sumOfWeights()));
+      scale(_h_pT_piminus,    1./(2*M_PI*sumOfWeights()));
+      scale(_h_pT_proton,     1./(2*M_PI*sumOfWeights()));
+      scale(_h_pT_antiproton, 1./(2*M_PI*sumOfWeights()));
+    }
+
+  private:
+
+    AIDA::IHistogram1D * _h_pT_piplus;
+    AIDA::IHistogram1D * _h_pT_piminus;
+    AIDA::IHistogram1D * _h_pT_proton;
+    AIDA::IHistogram1D * _h_pT_antiproton;
+
+  };
+
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<STAR_2006_S6500200> plugin_STAR_2006_S6500200;
+
+}


More information about the Rivet-svn mailing list