|
[Rivet-svn] r2696 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgSun Sep 19 23:35:03 BST 2010
Author: holsch Date: Sun Sep 19 23:35:01 2010 New Revision: 2696 Log: Add STAR analysis which makes use of the LossyFinalState projection to account for track-reco- and vertex-efficiencies. Data have been read from plot. There is something going on with the 0 bin in the multiplicity distribution, looks fine otherwise. Not sure about the normalisation to data. Added: trunk/data/anainfo/STAR_2008_S7869363.info trunk/data/plotinfo/STAR_2008_S7869363.plot trunk/data/refdata/STAR_2008_S7869363.aida trunk/src/Analyses/STAR_2008_S7869363.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Sep 16 09:07:06 2010 (r2695) +++ trunk/data/anainfo/Makefile.am Sun Sep 19 23:35:01 2010 (r2696) @@ -94,6 +94,7 @@ STAR_2006_S6500200.info \ STAR_2006_S6860818.info \ STAR_2006_S6870392.info \ + STAR_2008_S7869363.info \ STAR_2008_S7993412.info \ STAR_2009_UE_HELEN.info \ ZEUS_2001_S4815815.info \ Added: trunk/data/anainfo/STAR_2008_S7869363.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/STAR_2008_S7869363.info Sun Sep 19 23:35:01 2010 (r2696) @@ -0,0 +1,39 @@ +Name: STAR_2008_S7869363 +Year: 2008 +Summary: Multiplicities and pT spectra from STAR for pp at 200 GeV +Experiment: STAR +Collider: RHIC +SpiresID: 7869363 +Status: UNVALIDATED +Authors: + - Holger Schulz <holger.schulz at physik.hu-berlin.de> +References: + - arXiv:0808.2041 + - http://drupal.star.bnl.gov/STAR/files/starpublications/124/data.html +RunInfo: + QCD (pp) events at 200 GeV +NumEvents: 1.000.000 +Beams: [p+, p+] +Energies: [200] +PtCuts: +Description: + 'Charged Mulitplicity and identified charged particle spectra' +BibKey: :2008ez +BibTeX: '@Article{:2008ez, + author = "Abelev, B. I. and others", + collaboration = "STAR", + title = "{Systematic Measurements of Identified Particle Spectra in + $p p, d^+$ Au and Au+Au Collisions from STAR}", + journal = "Phys. Rev.", + volume = "C79", + year = "2009", + pages = "034909", + eprint = "0808.2041", + archivePrefix = "arXiv", + primaryClass = "nucl-ex", + doi = "10.1103/PhysRevC.79.034909", + SLACcitation = "%%CITATION = 0808.2041;%%" +}' +ToDo: + - Understand first bin in multiplicity distribution + Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Sep 16 09:07:06 2010 (r2695) +++ trunk/data/plotinfo/Makefile.am Sun Sep 19 23:35:01 2010 (r2696) @@ -92,6 +92,7 @@ STAR_2006_S6500200.plot \ STAR_2006_S6860818.plot \ STAR_2006_S6870392.plot \ + STAR_2008_S7869363.plot \ STAR_2008_S7993412.plot \ STAR_2009_UE_HELEN.plot \ ZEUS_2001_S4815815.plot \ Added: trunk/data/plotinfo/STAR_2008_S7869363.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/STAR_2008_S7869363.plot Sun Sep 19 23:35:01 2010 (r2696) @@ -0,0 +1,43 @@ +# BEGIN PLOT /STAR_2008_S7869363/d01-x01-y01 +Title=Raw charged multiplicity ($|\eta|<0.5 \quad p_\perp>0.2$ [GeV]) +XLabel=$N_\text{ch}^\text{raw}$ +YLabel=$P(N_\text{ch}^\text{raw})$ +# END PLOT + +# BEGIN PLOT /STAR_2008_S7869363/d02-x01-y01 +Title=$\pi^-\quad p_\perp$ spectrum +XLabel=$p_\perp$ [\unit{GeV}] +YLabel=$\frac{1}{2\pi} \frac{\text{d}^2N}{p_\perp \text{d}p_\perp \text{d}y} [\text{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /STAR_2008_S7869363/d02-x01-y02 +Title=$\pi^+\quad p_\perp$ spectrum +XLabel=$p_\perp$ [\unit{GeV}] +YLabel=$\frac{1}{2\pi} \frac{\text{d}^2N}{p_\perp \text{d}p_\perp \text{d}y} [\text{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /STAR_2008_S7869363/d02-x01-y03 +Title=$K^-\quad p_\perp$ spectrum +XLabel=$p_\perp$ [\unit{GeV}] +YLabel=$\frac{1}{2\pi} \frac{\text{d}^2N}{p_\perp \text{d}p_\perp \text{d}y} [\text{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /STAR_2008_S7869363/d02-x01-y04 +Title=$K^+\quad p_\perp$ spectrum +XLabel=$p_\perp$ [\unit{GeV}] +YLabel=$\frac{1}{2\pi} \frac{\text{d}^2N}{p_\perp \text{d}p_\perp \text{d}y} [\text{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /STAR_2008_S7869363/d02-x01-y05 +Title=Antiproton $p_\perp$ spectrum +XLabel=$p_\perp$ [\unit{GeV}] +YLabel=$\frac{1}{2\pi} \frac{\text{d}^2N}{p_\perp \text{d}p_\perp \text{d}y} [\text{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /STAR_2008_S7869363/d02-x01-y06 +Title=Proton $p_\perp$ spectrum +XLabel=$p_\perp$ [\unit{GeV}] +YLabel=$\frac{1}{2\pi} \frac{\text{d}^2N}{p_\perp \text{d}p_\perp \text{d}y} [\text{GeV}^{-2}]$ +# END PLOT + +# ... add more histograms as you need them ... Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Sep 16 09:07:06 2010 (r2695) +++ trunk/data/refdata/Makefile.am Sun Sep 19 23:35:01 2010 (r2696) @@ -96,6 +96,7 @@ E735_1993_S2896508.aida \ E735_1998_S3905616.aida \ SFM_1984_S1178091.aida \ + STAR_2008_S7869363.aida \ PDG_HADRON_MULTIPLICITIES.aida \ PDG_HADRON_MULTIPLICITIES_RATIOS.aida \ BELLE_2006_S6265367.aida Added: trunk/data/refdata/STAR_2008_S7869363.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/STAR_2008_S7869363.aida Sun Sep 19 23:35:01 2010 (r2696) @@ -0,0 +1,479 @@ +<?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_2008_S7869363" title="Charged particle multiplicity distribution"> + <dimension dim="0" title="$N_\text{ch}$" /> + <dimension dim="1" title="$P(N_\text{ch})$" /> + <annotation> + <item key="Title" value="Charged particle multiplicity distribution" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2008_S7869363/d01-x01-y01" sticky="true"/> + <item key="FullPath" value="/STAR_2008_S7869363.aida/REF/STAR_2008_S7869363/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="0.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="2.070000e-04" value="1.563670e-01" errorMinus="2.070000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="2.630000e-04" value="2.598680e-01" errorMinus="2.630000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="2.230000e-04" value="2.003350e-01" errorMinus="2.230000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="3.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.890000e-04" value="1.483870e-01" errorMinus="1.890000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="4.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.570000e-04" value="1.046980e-01" errorMinus="1.570000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="5.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.300000e-04" value="7.193700e-02" errorMinus="1.300000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="6.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.060000e-04" value="4.816900e-02" errorMinus="1.060000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="7.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="8.700000e-05" value="3.216000e-02" errorMinus="8.700000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="8.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="7.000000e-05" value="2.089900e-02" errorMinus="7.000000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="9.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="5.600000e-05" value="1.340700e-02" errorMinus="5.600000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.000000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="4.500000e-05" value="8.675000e-03" errorMinus="4.500000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.100000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="3.500000e-05" value="5.349000e-03" errorMinus="3.500000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.200000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="2.800000e-05" value="3.346000e-03" errorMinus="2.800000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.300000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="2.200000e-05" value="2.057000e-03" errorMinus="2.200000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.400000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.700000e-05" value="1.234000e-03" errorMinus="1.700000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.500000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.300000e-05" value="7.370000e-04" errorMinus="1.300000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.600000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.000000e-05" value="4.500000e-04" errorMinus="1.000000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.700000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="8.000000e-06" value="2.670000e-04" errorMinus="8.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.800000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="6.000000e-06" value="1.500000e-04" errorMinus="6.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.900000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="5.000000e-06" value="9.000000e-05" errorMinus="5.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.000000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="4.000000e-06" value="5.300000e-05" errorMinus="4.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.100000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="3.000000e-06" value="3.200000e-05" errorMinus="3.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.200000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="2.000000e-06" value="2.000000e-05" errorMinus="2.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.300000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.000000e-06" value="7.000000e-06" errorMinus="1.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.400000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.000000e-06" value="4.000000e-06" errorMinus="1.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.500000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.000000e-06" value="2.000000e-06" errorMinus="1.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.600000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.000000e-06" value="2.000000e-06" errorMinus="1.000000e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.700000e+01" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.000000e-06" value="1.000000e-06" errorMinus="1.000000e-06"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" + path="/REF/STAR_2008_S7869363" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2008_S7869363/d02-x01-y01" sticky="true"/> + <item key="FullPath" value="/STAR_2008_S7869363.aida/REF/STAR_2008_S7869363/d02-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.000000e-02" value="2.020000e+00" errorMinus="6.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-02" value="1.520000e+00" errorMinus="3.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.000000e-02" value="1.130000e+00" errorMinus="2.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="9.000000e-03" value="8.440000e-01" errorMinus="9.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="7.000000e-03" value="6.350000e-01" errorMinus="7.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-03" value="4.690000e-01" errorMinus="5.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-03" value="3.540000e-01" errorMinus="4.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-03" value="2.670000e-01" errorMinus="3.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-03" value="2.020000e-01" errorMinus="4.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-03" value="1.530000e-01" errorMinus="3.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-03" value="1.160000e-01" errorMinus="3.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-03" value="8.650000e-02" errorMinus="3.000000e-03"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y02" dimension="2" + path="/REF/STAR_2008_S7869363" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2008_S7869363/d02-x01-y02" sticky="true"/> + <item key="FullPath" value="/STAR_2008_S7869363.aida/REF/STAR_2008_S7869363/d02-x01-y02" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.000000e-02" value="2.070000e+00" errorMinus="6.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-02" value="1.540000e+00" errorMinus="3.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.000000e-02" value="1.140000e+00" errorMinus="2.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="9.000000e-03" value="8.570000e-01" errorMinus="9.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="7.000000e-03" value="6.380000e-01" errorMinus="7.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-03" value="4.760000e-01" errorMinus="5.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-03" value="3.590000e-01" errorMinus="4.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-03" value="2.730000e-01" errorMinus="3.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-03" value="2.070000e-01" errorMinus="4.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-03" value="1.550000e-01" errorMinus="3.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.000000e-03" value="1.160000e-01" errorMinus="3.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.100000e-03" value="8.950000e-02" errorMinus="3.100000e-03"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y03" dimension="2" + path="/REF/STAR_2008_S7869363" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2008_S7869363/d02-x01-y03" sticky="true"/> + <item key="FullPath" value="/STAR_2008_S7869363.aida/REF/STAR_2008_S7869363/d02-x01-y03" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.100000e-02" value="1.430000e-01" errorMinus="1.100000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-03" value="1.260000e-01" errorMinus="5.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.000000e-03" value="1.080000e-01" errorMinus="2.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.400000e-03" value="8.770000e-02" errorMinus="2.400000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.600000e-03" value="7.340000e-02" errorMinus="2.600000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.300000e-03" value="6.170000e-02" errorMinus="6.300000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-03" value="4.870000e-02" errorMinus="5.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.200000e-03" value="4.110000e-02" errorMinus="4.200000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.800000e-03" value="3.700000e-02" errorMinus="2.800000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.200000e-03" value="2.860000e-02" errorMinus="2.200000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.700000e-03" value="2.420000e-02" errorMinus="2.700000e-03"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y04" dimension="2" + path="/REF/STAR_2008_S7869363" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2008_S7869363/d02-x01-y04" sticky="true"/> + <item key="FullPath" value="/STAR_2008_S7869363.aida/REF/STAR_2008_S7869363/d02-x01-y04" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.100000e-02" value="1.520000e-01" errorMinus="1.100000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-03" value="1.300000e-01" errorMinus="5.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.000000e-03" value="1.080000e-01" errorMinus="2.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.500000e-03" value="9.160000e-02" errorMinus="2.500000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.700000e-03" value="7.470000e-02" errorMinus="2.700000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.400000e-03" value="6.260000e-02" errorMinus="6.400000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.400000e-03" value="5.260000e-02" errorMinus="5.400000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.500000e-03" value="4.410000e-02" errorMinus="4.500000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.800000e-03" value="3.810000e-02" errorMinus="2.800000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.400000e-03" value="3.060000e-02" errorMinus="2.400000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.800000e-03" value="2.490000e-02" errorMinus="2.800000e-03"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y05" dimension="2" + path="/REF/STAR_2008_S7869363" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2008_S7869363/d02-x01-y05" sticky="true"/> + <item key="FullPath" value="/STAR_2008_S7869363.aida/REF/STAR_2008_S7869363/d02-x01-y05" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.300000e-03" value="5.540000e-02" errorMinus="1.300000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.100000e-03" value="4.810000e-02" errorMinus="1.100000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.000000e-03" value="4.250000e-02" errorMinus="1.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="9.000000e-04" value="3.770000e-02" errorMinus="9.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="8.000000e-04" value="3.380000e-02" errorMinus="8.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="7.000000e-04" value="2.780000e-02" errorMinus="7.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.000000e-04" value="2.490000e-02" errorMinus="6.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.000000e-04" value="2.030000e-02" errorMinus="6.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.000000e-04" value="1.750000e-02" errorMinus="6.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-04" value="1.520000e-02" errorMinus="5.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-04" value="1.270000e-02" errorMinus="4.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-04" value="1.050000e-02" errorMinus="4.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.400000e-04" value="8.950000e-03" errorMinus="3.400000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.025000e+00" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.400000e-04" value="7.360000e-03" errorMinus="3.400000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.075000e+00" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.200000e-04" value="6.590000e-03" errorMinus="3.200000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.125000e+00" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.900000e-04" value="5.210000e-03" errorMinus="2.900000e-04"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y06" dimension="2" + path="/REF/STAR_2008_S7869363" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2008_S7869363/d02-x01-y06" sticky="true"/> + <item key="FullPath" value="/STAR_2008_S7869363.aida/REF/STAR_2008_S7869363/d02-x01-y06" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="2.500000e-03" value="5.070000e-02" errorMinus="2.500000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.900000e-03" value="4.690000e-02" errorMinus="1.900000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.400000e-03" value="4.080000e-02" errorMinus="1.400000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="1.000000e-03" value="3.420000e-02" errorMinus="1.000000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="7.000000e-04" value="2.870000e-02" errorMinus="7.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="7.000000e-04" value="2.410000e-02" errorMinus="7.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="6.000000e-04" value="2.130000e-02" errorMinus="6.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-04" value="1.820000e-02" errorMinus="5.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="5.000000e-04" value="1.540000e-02" errorMinus="5.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-04" value="1.310000e-02" errorMinus="4.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-04" value="1.110000e-02" errorMinus="4.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.025000e+00" errorMinus="2.500000e-02"/> + <measurement errorPlus="4.000000e-04" value="9.780000e-03" errorMinus="4.000000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.075000e+00" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.700000e-04" value="8.560000e-03" errorMinus="3.700000e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.125000e+00" errorMinus="2.500000e-02"/> + <measurement errorPlus="3.800000e-04" value="7.380000e-03" errorMinus="3.800000e-04"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Sep 16 09:07:06 2010 (r2695) +++ trunk/src/Analyses/Makefile.am Sun Sep 19 23:35:01 2010 (r2696) @@ -37,6 +37,7 @@ MC_TTBAR.cc \ OPAL_1993_S2692198.cc \ SFM_1984_S1178091.cc \ + STAR_2008_S7869363.cc \ STAR_2008_S7993412.cc \ STAR_2009_UE_HELEN.cc \ ZEUS_2001_S4815815.cc Added: trunk/src/Analyses/STAR_2008_S7869363.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/STAR_2008_S7869363.cc Sun Sep 19 23:35:01 2010 (r2696) @@ -0,0 +1,195 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/ChargedFinalState.hh" +#include "Rivet/Projections/LossyFinalState.hh" +#include "Rivet/Tools/ParticleIdUtils.hh" +/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder... + +namespace Rivet { + + class STARRandomFilter { + public: + + STARRandomFilter() + { + // Tracking efficiency, binned in pT of track between 0 and 600 MeV, spacing = 50 MeV + // above pT, the efficiency is constant at 88pct. + //_trkeff = {0,0,0.38,0.72,0.78,0.81,0.82,0.84,0.85,0.86,0.87,0.88}; + } + + // Return true to throw away a particle + bool operator()(const Particle& p) { + /// @todo Use a better RNG + size_t idx = floor(p.momentum().pT()/MeV/50); + if (idx > 11) idx = 11; + return (rand()/static_cast<double>(RAND_MAX) > _trkeff[idx]); + } + + + int compare(const STARRandomFilter& other) const { + return true; + } + + private: + + const static double _trkeff[12]; + + }; + + + // Here we have the track reconstruction efficiencies for tracks with pT from 0 to 600 MeV + // in steps of 50 MeV. The efficiency is assumed to be 0.88 for pT >= 600 MeV + const double STARRandomFilter::_trkeff[12] = {0,0,0.38,0.72,0.78,0.81,0.82,0.84,0.85,0.86,0.87,0.88}; + + + + class STAR_2008_S7869363 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + STAR_2008_S7869363() + : Analysis("STAR_2008_S7869363"), + nCutsPassed(0), nPi(0), nPiPlus(0), nKaon(0), nKaonPlus(0), nProton(0), nAntiProton(0) + { + /// @todo Set approriate for your analysis + setBeams(PROTON, PROTON); + /// @todo Set whether your finalize method needs the generator cross section + setNeedsCrossSection(false); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + const ChargedFinalState cfs(-0.5, 0.5, 0.2*GeV); + const LossyFinalState<STARRandomFilter> lfs(cfs, STARRandomFilter()); + addProjection(lfs, "FS"); + + _h_dNch = bookHistogram1D(1, 1, 1); + _h_dpT_Pi = bookHistogram1D(2, 1, 1); + _h_dpT_Piplus = bookHistogram1D(2, 1, 2); + _h_dpT_Kaon = bookHistogram1D(2, 1, 3); + _h_dpT_Kaonplus = bookHistogram1D(2, 1, 4); + _h_dpT_AntiProton = bookHistogram1D(2, 1, 5); + _h_dpT_Proton = bookHistogram1D(2, 1, 6); + + + + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& charged = applyProjection<FinalState>(event, "FS"); + + //if (charged.particles().size() == 0) { + //vetoEvent; + //} + + // Vertex reconstrucion efficiency, for events with more than 23 reconstructed tracks, + // the efficiency is constant at 100pct. + double vtxeffs[24] = {0.,0.512667,0.739365,0.847131,0.906946,0.940922,0.959328, + 0.96997,0.975838,0.984432,0.988311,0.990327,0.990758,0.995767,0.99412,0.992271, + 0.996631,0.994802,0.99635,0.997384,0.998986,0.996441,0.994513,1.}; + + double vtxeff = 1.0; + if (charged.particles().size() < 24) vtxeff = vtxeffs[charged.particles().size()]; + + const double weight = vtxeff * event.weight(); + + foreach (const Particle& p, charged.particles()) { + double pT = p.momentum().pT()/GeV; + double y = p.momentum().rapidity(); + + if ( fabs(y) < 0.1 ) { + nCutsPassed+=weight; + const PdgId id = p.pdgId(); + switch (id) { + case -211: + _h_dpT_Pi->fill(pT, weight/(TWOPI*pT*0.2)); + nPi+=weight; + break; + case 211: + _h_dpT_Piplus->fill(pT, weight/(TWOPI*pT*0.2)); + nPiPlus+=weight; + break; + case -321: + _h_dpT_Kaon->fill(pT, weight/(TWOPI*pT*0.2)); + nKaon+=weight; + break; + case 321: + _h_dpT_Kaonplus->fill(pT, weight/(TWOPI*pT*0.2)); + nKaonPlus+=weight; + break; + case -2212: + _h_dpT_AntiProton->fill(pT, weight/(TWOPI*pT*0.2)); + nAntiProton+=weight; + break; + case 2212: + _h_dpT_Proton->fill(pT, weight/(TWOPI*pT*0.2)); + nProton+=weight; + break; + } + } + else { + continue; + } + } + _h_dNch->fill(charged.particles().size(), weight); + } + + + /// Normalise histograms etc., after the run + void finalize() { + double nTot = nPi + nPiPlus + nKaon + nKaonPlus + nProton + nAntiProton; + + + normalize(_h_dNch); + + // Norm to data, I know, I know + + normalize(_h_dpT_Pi , 0.389825 ); + normalize(_h_dpT_Piplus , 0.396025 ); + normalize(_h_dpT_Kaon , 0.03897 ); + normalize(_h_dpT_Kaonplus , 0.04046 ); + normalize(_h_dpT_AntiProton, 0.0187255); + normalize(_h_dpT_Proton , 0.016511 ); + + + } + + + + private: + + + AIDA::IHistogram1D *_h_dNch; + + AIDA::IHistogram1D *_h_dpT_Pi, *_h_dpT_Piplus; + AIDA::IHistogram1D *_h_dpT_Kaon, *_h_dpT_Kaonplus; + AIDA::IHistogram1D *_h_dpT_AntiProton, *_h_dpT_Proton; + + AIDA::IProfile1D *_h_pT_vs_Nch; + double nCutsPassed, nPi, nPiPlus, nKaon, nKaonPlus, nProton, nAntiProton; + + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<STAR_2008_S7869363> plugin_STAR_2008_S7869363; + + +}
More information about the Rivet-svn mailing list |