|
[Rivet-svn] r3502 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Dec 5 15:36:00 GMT 2011
Author: fsiegert Date: Mon Dec 5 15:36:00 2011 New Revision: 3502 Log: Add LHCB_2011_I919315 (thanks to Michael Kaballo): Inclusive differential $\Phi$ production cross-section. Still UNVALIDATED since the data has yet to be submitted to HepData. Michael validated it by running Perugia0 through it and comparing with the plots in the paper. I have additionally run some Sherpa events, plots are at http://fsiegert.web.cern.ch/fsiegert/tmp/lhcb-phi/ Added: trunk/data/anainfo/LHCB_2011_I919315.info trunk/data/plotinfo/LHCB_2011_I919315.plot trunk/data/refdata/LHCB_2011_I919315.aida trunk/src/Analyses/LHCB_2011_I919315.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Added: trunk/data/anainfo/LHCB_2011_I919315.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/LHCB_2011_I919315.info Mon Dec 5 15:36:00 2011 (r3502) @@ -0,0 +1,42 @@ +Name: LHCB_2011_I919315 +Year: 2011 +Summary: Inclusive differential $\Phi$ production cross-section as a function of $p_\text{T}$ and $y$ +Experiment: LHCB +Collider: LHC +InspireID: 919315 +Status: UNVALIDATED +Authors: + - Friederike Blatt friederike.blatt at tu-dortmund.de + - Michael Kaballo michael.kaballo at tu-dortmund.de + - Till Moritz Karbach moritz.karbach at tu-dortmund.de +References: + - Phys.Lett.B703:267-273,2011 + - arXiv:1107.3935 [hep-ex] +RunInfo: + $pp$ collisions, QCD-Events, $\sqrt{s}=7$TeV +Beams: [p+, p+] +Energies: [7000] +Description: + 'Measurement of the inclusive differential $\Phi$ cross-section in $pp$ + collisions at $\sqrt {s}=7$TeV in the rapidity range of + $ 2.44 < y < 4.06$ and the $p_\text{T}$ range of + 0.6 GeV/c < $p_\text{T}$ < 5.0 GeV/c.' +NeedCrossSection: yes +BibKey: Aaij:2011uk +BibTeX: '@article{Aaij:2011uk, + author = "Aaij, R. and others", + title = "{Measurement of the inclusive phi cross-section in pp + collisions at sqrt(s) = 7 TeV}", + collaboration = "LHCb Collaboration", + journal = "Phys.Lett.", + volume = "B703", + pages = "267-273", + doi = "10.1016/j.physletb.2011.08.017", + year = "2011", + note = "* Temporary entry *", + eprint = "1107.3935", + archivePrefix = "arXiv", + primaryClass = "Unknown", +}' +ToDo: + - Get reference data into HepData \ No newline at end of file Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Mon Dec 5 14:48:53 2011 (r3501) +++ trunk/data/anainfo/Makefile.am Mon Dec 5 15:36:00 2011 (r3502) @@ -112,6 +112,7 @@ JADE_OPAL_2000_S4300807.info \ JADE_1998_S3612880.info \ LHCB_2010_S8758301.info \ + LHCB_2011_I919315.info \ MC_GENERIC.info \ MC_LEADINGJETS.info \ MC_DIJET.info \ Added: trunk/data/plotinfo/LHCB_2011_I919315.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/LHCB_2011_I919315.plot Mon Dec 5 15:36:00 2011 (r3502) @@ -0,0 +1,56 @@ +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y.. +Title=Transverse momentum of $\Phi$-mesons +XLabel=$p_T^\Phi$ [MeV] +YLabel=$\text{d}\sigma /\text{d}p_\perp$ $[\mu \text{b}/\text{MeV}]$ +LogY=1 +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x02-y01 +Title=Rapidity of $\Phi$-mesons +XLabel=$y^\Phi$ +YLabel=$\text{d}\sigma /\text{d}y$ $[\mu \text{b}]$ +LogY=0 +CustomLegend=$(0.6 < p_\perp^\Phi/\text{GeV} < 5.0)$ +LegendXPos=0.05 +LegendYPos=0.4 +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y10 +CustomLegend=$(2.44<y^\Phi<4.06)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y01 +CustomLegend=$(2.44<y^\Phi<2.62)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y02 +CustomLegend=$(2.62<y^\Phi<2.80)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y03 +CustomLegend=$(2.80<y^\Phi<2.98)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y04 +CustomLegend=$(2.98<y^\Phi<3.16)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y05 +CustomLegend=$(3.16<y^\Phi<3.34)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y06 +CustomLegend=$(3.34<y^\Phi<3.52)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y07 +CustomLegend=$(3.52<y^\Phi<3.70)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y08 +CustomLegend=$(3.70<y^\Phi<3.88)$ +# END PLOT + +# BEGIN PLOT /LHCB_2011_I919315/d01-x01-y09 +CustomLegend=$(3.88<y^\Phi<4.06)$ +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Mon Dec 5 14:48:53 2011 (r3501) +++ trunk/data/plotinfo/Makefile.am Mon Dec 5 15:36:00 2011 (r3502) @@ -103,6 +103,7 @@ JADE_OPAL_2000_S4300807.plot \ JADE_1998_S3612880.plot \ LHCB_2010_S8758301.plot \ + LHCB_2011_I919315.plot \ MC_GENERIC.plot \ MC_LEADINGJETS.plot \ MC_DIJET.plot \ Added: trunk/data/refdata/LHCB_2011_I919315.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/LHCB_2011_I919315.aida Mon Dec 5 15:36:00 2011 (r3502) @@ -0,0 +1,609 @@ +<?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/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y01" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.001000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="9.590000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.580000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.480000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.690000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.220000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.340000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.090000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.270000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.800000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.000000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.400000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y02" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y02" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y02" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.530000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.970000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.760000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.270000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.110000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.810000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.230000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.920000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.120000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.900000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.800000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.400000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y03" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y03" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y03" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.069000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.190000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.950000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.040000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.210000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.090000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.760000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.010000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.110000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.300000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.400000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.100000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y04" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y04" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y04" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.171000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.032000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.180000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.480000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.840000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.080000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.200000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.060000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.090000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.500000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.100000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.000000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y05" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y05" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y05" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.060000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.620000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.510000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.930000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.990000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.820000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.080000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.940000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.060000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.700000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.900000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.000000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y06" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y06" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y06" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.131000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.170000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.810000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.610000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.700000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.480000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.550000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.690000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.060000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.300000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.500000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="9.000000e-03" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y07" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y07" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y07" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.341000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.160000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.850000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.090000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.840000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.360000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.310000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.640000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.200000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.900000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.200000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="8.000000e-03" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y08" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y08" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y08" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.164000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.065000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="6.900000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.610000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.330000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.150000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.280000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.400000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.800000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.900000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.900000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.000000e-03" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y09" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y09" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y09" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.341000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="9.750000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.600000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.310000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.090000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.790000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.130000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.310000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.000000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="3.900000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.200000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.000000e-03" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y10" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn/MeV/c, Projection"> + <annotation> + <item key="Title" value="Data, SIG in mubarn/MeV/c, Projection" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x01-y10" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x01-y10" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="7.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.823500e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="9.000000e+02" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.529000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.100000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.262300e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.300000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.004000e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.500000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="7.704000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.700000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="5.901000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+02" value="1.900000e+03" errorMinus="1.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.477000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.200000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="2.894000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="2.600000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.623000e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.000000e+02" value="3.000000e+03" errorMinus="2.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="9.600000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.000000e+02" value="3.600000e+03" errorMinus="4.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="4.670000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+02" value="4.500000e+03" errorMinus="5.000000e+02"/> + <measurement errorPlus="0.000000e+00" value="1.610000e-02" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x02-y01" dimension="2" + path="/REF/LHCB_2011_I919315" title="Data, SIG in mubarn, Projection"> + <annotation> + <item key="Title" value="Data, SIG in mubarn, Projection" sticky="true"/> + <item key="AidaPath" value="/REF/LHCB_2011_I919315/d01-x02-y01" sticky="true"/> + <item key="FullPath" value="/LHCB_2011_I919315.aida/REF/LHCB_2011_I919315/d01-x02-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="2.530000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.130000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="2.710000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.047000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="2.890000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.087000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="3.070000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.163000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="3.250000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.107000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="3.430000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.123000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="3.610000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.068000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="3.790000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.021000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="9.000000e-02" value="3.970000e+00" errorMinus="9.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.022000e+03" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Mon Dec 5 14:48:53 2011 (r3501) +++ trunk/data/refdata/Makefile.am Mon Dec 5 15:36:00 2011 (r3502) @@ -1,5 +1,6 @@ dist_pkgdata_DATA = \ LHCB_2010_S8758301.aida \ + LHCB_2011_I919315.aida \ JADE_OPAL_2000_S4300807.aida \ JADE_1998_S3612880.aida \ TASSO_1990_S2148048.aida \ Added: trunk/src/Analyses/LHCB_2011_I919315.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/LHCB_2011_I919315.cc Mon Dec 5 15:36:00 2011 (r3502) @@ -0,0 +1,98 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/BinnedHistogram.hh" +#include "Rivet/Projections/UnstableFinalState.hh" +#include "Rivet/Tools/ParticleIdUtils.hh" + +namespace Rivet { + + + class LHCB_2011_I919315 : public Analysis { + public: + /// @name Constructors etc. + //@{ + + /// Constructor + LHCB_2011_I919315() + : Analysis("LHCB_2011_I919315") + { + } + + //@} + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + addProjection(UnstableFinalState(), "UFS"); + + _h_Phi_pT_y.addHistogram( 2.44, 2.62, bookHistogram1D(1, 1, 1)); + _h_Phi_pT_y.addHistogram( 2.62, 2.80, bookHistogram1D(1, 1, 2)); + _h_Phi_pT_y.addHistogram( 2.80, 2.98, bookHistogram1D(1, 1, 3)); + _h_Phi_pT_y.addHistogram( 2.98, 3.16, bookHistogram1D(1, 1, 4)); + _h_Phi_pT_y.addHistogram( 3.16, 3.34, bookHistogram1D(1, 1, 5)); + _h_Phi_pT_y.addHistogram( 3.34, 3.52, bookHistogram1D(1, 1, 6)); + _h_Phi_pT_y.addHistogram( 3.52, 3.70, bookHistogram1D(1, 1, 7)); + _h_Phi_pT_y.addHistogram( 3.70, 3.88, bookHistogram1D(1, 1, 8)); + _h_Phi_pT_y.addHistogram( 3.88, 4.06, bookHistogram1D(1, 1, 9)); + _h_Phi_pT = bookHistogram1D(1, 1, 10); + _h_Phi_y = bookHistogram1D(1, 2, 1); + } + + + /// Perform the per-event analysis + void analyze (const Event& event) { + const double weight = event.weight(); + const UnstableFinalState& ufs = applyProjection<UnstableFinalState> (event, "UFS"); + + foreach (const Particle& p, ufs.particles()) { + const PdgId id = abs(p.pdgId()); + + if (id == 333) { // id 333 = phi-meson + double y = p.momentum().rapidity(); + double pT = p.momentum().perp(); + + if (pT < 0.6*GeV || pT > 5.0*GeV || y < 2.44 || y > 4.06) { + continue; + } + + _h_Phi_y->fill (y, weight); + _h_Phi_pT->fill (pT/MeV, weight); + _h_Phi_pT_y.fill(y, pT/MeV, weight); + } + } + } + + /// Normalise histograms etc., after the run + void finalize() { + + double scale_factor = crossSectionPerEvent()/microbarn; + scale (_h_Phi_y, scale_factor); + scale (_h_Phi_pT, scale_factor); + _h_Phi_pT_y.scale(scale_factor, this); + } + + //@} + + private: + + /// @name Histograms + //@{ + AIDA::IHistogram1D * _h_Phi_y; + AIDA::IHistogram1D * _h_Phi_pT; + BinnedHistogram<double> _h_Phi_pT_y; + //@} + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(LHCB_2011_I919315); + +} + +//@} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Mon Dec 5 14:48:53 2011 (r3501) +++ trunk/src/Analyses/Makefile.am Mon Dec 5 15:36:00 2011 (r3502) @@ -93,7 +93,8 @@ RivetLHCbAnalyses_la_SOURCES = if ENABLE_UNVALIDATED RivetLHCbAnalyses_la_SOURCES += \ - LHCB_2010_S8758301.cc + LHCB_2010_S8758301.cc \ + LHCB_2011_I919315.cc endif
More information about the Rivet-svn mailing list |