|
[Rivet-svn] r1750 - in trunk: data/anainfo data/plotinfo data/refdata include/Rivet include/Rivet/Analyses src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgFri Aug 7 17:29:42 BST 2009
Author: fsiegert Date: Fri Aug 7 17:29:42 2009 New Revision: 1750 Log: Add analysis CDF_2005_S6080774 for measurements in prompt diphoton production. Plots look good compared to a merged Sherpa sample, so I set status to VALIDATED. Added: trunk/data/anainfo/CDF_2005_S6080774.info trunk/data/plotinfo/CDF_2005_S6080774.plot trunk/data/refdata/CDF_2005_S6080774.aida trunk/include/Rivet/Analyses/CDF_2005_S6080774.hh trunk/src/Analyses/CDF_2005_S6080774.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/include/Rivet/Makefile.am trunk/src/Analyses/Makefile.am trunk/src/Analyses/StdAnalyses.cc Added: trunk/data/anainfo/CDF_2005_S6080774.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_2005_S6080774.info Fri Aug 7 17:29:42 2009 (r1750) @@ -0,0 +1,25 @@ +Name: CDF_2005_S6080774 +Year: 2005 +Summary: Differential cross sections for Prompt Diphoton Production +Experiment: CDF +Collider: Tevatron Run 2 +SpiresID: 6080774 +Status: VALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - Phys. Rev. Lett. 95, 022003 + - doi:10.1103/PhysRevLett.95.022003 + - arXiv:hep-ex/0412050 +RunInfo: + p pbar -> photon photon [+ jets] at 1960 GeV. The analysis uses photons + with pT larger then 13 GeV. To allow for shifts in the shower, the ME + cut on the transverse photon momentum shouldn't be too hard, e.g. 5 GeV. +NumEvents: 1000000 +PtCuts: [5] +Description: + Measurement of the cross section of prompt diphoton production in p-pbar + collisions at $\sqrt{s} = 1.96$ TeV using a data sample of + $207 \mathrm{pb}^{-1}$ as a function of the diphoton mass, the transverse + momentum of the diphoton system, and the azimuthal angle between the two + photons. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Aug 6 22:05:44 2009 (r1749) +++ trunk/data/anainfo/Makefile.am Fri Aug 7 17:29:42 2009 (r1750) @@ -8,6 +8,7 @@ CDF_2001_S4751469.info \ CDF_2002_S4796047.info \ CDF_2004_S5839831.info \ + CDF_2005_S6080774.info \ CDF_2005_S6217184.info \ CDF_2006_S6450792.info \ CDF_2006_S6653332.info \ Added: trunk/data/plotinfo/CDF_2005_S6080774.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_2005_S6080774.plot Fri Aug 7 17:29:42 2009 (r1750) @@ -0,0 +1,71 @@ +# BEGIN PLOT /CDF_2005_S6080774/d01-x01-y01 +Title=Invariant mass of diphoton pair +XLabel=$M_{\gamma \gamma}$ (GeV/$c^2$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}M_{\gamma \gamma}$ (pb/(GeV/$c^2$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d01-x01-y02 +Title=Invariant mass of diphoton pair (compared to DIPHOX) +XLabel=$M_{\gamma \gamma}$ (GeV/$c^2$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}M_{\gamma \gamma}$ (pb/(GeV/$c^2$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d01-x01-y03 +Title=Invariant mass of diphoton pair (compared to RESBOS) +XLabel=$M_{\gamma \gamma}$ (GeV/$c^2$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}M_{\gamma \gamma}$ (pb/(GeV/$c^2$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d01-x01-y04 +Title=Invariant mass of diphoton pair (compared to PYTHIA) +XLabel=$M_{\gamma \gamma}$ (GeV/$c^2$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}M_{\gamma \gamma}$ (pb/(GeV/$c^2$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d02-x01-y01 +Title=Transverse momentum of diphoton pair +XLabel=$p_{\perp}^{\gamma \gamma}$ (GeV/$c$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_{\perp}^{\gamma \gamma}$ (pb/(GeV/$c$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d02-x01-y02 +Title=Transverse momentum of diphoton pair (compared to DIPHOX) +XLabel=$p_{\perp}^{\gamma \gamma}$ (GeV/$c$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_{\perp}^{\gamma \gamma}$ (pb/(GeV/$c$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d02-x01-y03 +Title=Transverse momentum of diphoton pair (compared to RESBOS) +XLabel=$p_{\perp}^{\gamma \gamma}$ (GeV/$c$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_{\perp}^{\gamma \gamma}$ (pb/(GeV/$c$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d02-x01-y04 +Title=Transverse momentum of diphoton pair (compared to PYTHIA) +XLabel=$p_{\perp}^{\gamma \gamma}$ (GeV/$c$) +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_{\perp}^{\gamma \gamma}$ (pb/(GeV/$c$)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d03-x01-y01 +Title=Azimuthal angle between photons +XLabel=$\Delta\Phi_{\gamma \gamma}$ ($\pi\cdot$rad) +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\Phi_{\gamma \gamma}$ (pb/($\pi\cdot$rad)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d03-x01-y02 +Title=Azimuthal angle between photons (compared to DIPHOX) +XLabel=$\Delta\Phi_{\gamma \gamma}$ ($\pi\cdot$rad) +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\Phi_{\gamma \gamma}$ (pb/($\pi\cdot$rad)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d03-x01-y03 +Title=Azimuthal angle between photons (compared to RESBOS) +XLabel=$\Delta\Phi_{\gamma \gamma}$ ($\pi\cdot$rad) +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\Phi_{\gamma \gamma}$ (pb/($\pi\cdot$rad)) +# END PLOT + +# BEGIN PLOT /CDF_2005_S6080774/d03-x01-y04 +Title=Azimuthal angle between photons (compared to PYTHIA) +XLabel=$\Delta\Phi_{\gamma \gamma}$ ($\pi\cdot$rad) +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\Phi_{\gamma \gamma}$ (pb/($\pi\cdot$rad)) +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Aug 6 22:05:44 2009 (r1749) +++ trunk/data/plotinfo/Makefile.am Fri Aug 7 17:29:42 2009 (r1750) @@ -8,6 +8,7 @@ CDF_2001_S4751469.plot \ CDF_2002_S4796047.plot \ CDF_2004_S5839831.plot \ + CDF_2005_S6080774.plot \ CDF_2005_S6217184.plot \ CDF_2006_S6450792.plot \ CDF_2006_S6653332.plot \ Added: trunk/data/refdata/CDF_2005_S6080774.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_2005_S6080774.aida Fri Aug 7 17:29:42 2009 (r1750) @@ -0,0 +1,349 @@ +<?xml version="1.0" encoding="utf-8" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.0" package="HepData"/> + <dataPointSet name="d01-x01-y01" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG)/DM IN PB/GEV" > + <dataPoint> + <measurement value="17.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.03" errorPlus="0.03162277660168379" errorMinus="0.03162277660168379"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.44" errorPlus="0.17691806012954134" errorMinus="0.17691806012954134"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.61" errorPlus="0.23345235059857505" errorMinus="0.23345235059857505"/> + </dataPoint> + <dataPoint> + <measurement value="40.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.46" errorPlus="0.17204650534085256" errorMinus="0.17204650534085256"/> + </dataPoint> + <dataPoint> + <measurement value="52.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.16" errorPlus="0.06403124237432849" errorMinus="0.06403124237432849"/> + </dataPoint> + <dataPoint> + <measurement value="80.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.01" errorPlus="0.022360679774997897" errorMinus="0.022360679774997897"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y02" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=DIPHOX))/DM IN PB/GEV" > + <dataPoint> + <measurement value="17.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.04" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.41" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.7" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="40.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.46" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="52.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.19" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="80.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.04" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y03" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=RESBOS))/DM IN PB/GEV" > + <dataPoint> + <measurement value="17.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.01" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.31" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.65" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="40.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.43" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="52.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.16" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="80.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.04" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y04" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=PYTHIA))/DM IN PB/GEV" > + <dataPoint> + <measurement value="17.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.01" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.18" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.38" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="40.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.24" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="52.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="0.09" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="80.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.02" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG)/DPT IN PB/GEV" > + <dataPoint> + <measurement value="0.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.7" errorPlus="0.33105890714493696" errorMinus="0.33105890714493696"/> + </dataPoint> + <dataPoint> + <measurement value="1.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="1.18" errorPlus="0.5131276644266999" errorMinus="0.5131276644266999"/> + </dataPoint> + <dataPoint> + <measurement value="3.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.92" errorPlus="0.4482186966202994" errorMinus="0.4482186966202994"/> + </dataPoint> + <dataPoint> + <measurement value="6.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.96" errorPlus="0.3940812099047606" errorMinus="0.3940812099047606"/> + </dataPoint> + <dataPoint> + <measurement value="10.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.29" errorPlus="0.2469817807045694" errorMinus="0.2469817807045694"/> + </dataPoint> + <dataPoint> + <measurement value="14.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.42" errorPlus="0.18439088914585774" errorMinus="0.18439088914585774"/> + </dataPoint> + <dataPoint> + <measurement value="20.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.19" errorPlus="0.10295630140987" errorMinus="0.10295630140987"/> + </dataPoint> + <dataPoint> + <measurement value="28.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.12" errorPlus="0.0670820393249937" errorMinus="0.0670820393249937"/> + </dataPoint> + <dataPoint> + <measurement value="36.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.1" errorPlus="0.07071067811865477" errorMinus="0.07071067811865477"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y02" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=DIPHOX))/DPT IN PB/GEV" > + <dataPoint> + <measurement value="0.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="-2.45" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="1.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="5.59" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="3.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="2.06" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="6.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="1.17" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="10.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.44" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="14.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.24" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="20.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.13" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="28.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.09" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="36.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.06" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y03" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=RESBOS))/DPT IN PB/GEV" > + <dataPoint> + <measurement value="0.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.34" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="1.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.95" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="3.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="1.03" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="6.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.94" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="10.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.59" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="14.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.36" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="20.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.19" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="28.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.07" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="36.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.03" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y04" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=PYTHIA))/DPT IN PB/GEV" > + <dataPoint> + <measurement value="0.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.53" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="1.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="1.15" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="3.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.94" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="6.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.46" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="10.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.21" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="14.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="0.12" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="20.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.07" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="28.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.03" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="36.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="0.01" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG)/D(PHI(P=3)-PHI(P=4)) IN PB/RAD" > + <dataPoint> + <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/> + <measurement value="1.06" errorPlus="0.6212889826803627" errorMinus="0.6212889826803627"/> + </dataPoint> + <dataPoint> + <measurement value="0.30000000000000004" errorPlus="0.09999999999999998" errorMinus="0.10000000000000003"/> + <measurement value="0.89" errorPlus="0.6105735008989499" errorMinus="0.6105735008989499"/> + </dataPoint> + <dataPoint> + <measurement value="0.5" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="0.51" errorPlus="0.6580273550544841" errorMinus="0.6580273550544841"/> + </dataPoint> + <dataPoint> + <measurement value="0.7" errorPlus="0.10000000000000009" errorMinus="0.09999999999999998"/> + <measurement value="3.34" errorPlus="1.5138031576133009" errorMinus="1.5138031576133009"/> + </dataPoint> + <dataPoint> + <measurement value="0.9" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="15.56" errorPlus="5.366386121031546" errorMinus="5.366386121031546"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y02" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=DIPHOX))/D(PHI(P=3)-PHI(P=4)) IN PB/RAD" > + <dataPoint> + <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/> + <measurement value="0.69" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.30000000000000004" errorPlus="0.09999999999999998" errorMinus="0.10000000000000003"/> + <measurement value="0.56" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.5" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="0.71" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.7" errorPlus="0.10000000000000009" errorMinus="0.09999999999999998"/> + <measurement value="1.83" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.9" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="23.37" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y03" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=RESBOS))/D(PHI(P=3)-PHI(P=4)) IN PB/RAD" > + <dataPoint> + <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/> + <measurement value="0.01" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.30000000000000004" errorPlus="0.09999999999999998" errorMinus="0.10000000000000003"/> + <measurement value="0.23" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.5" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="0.73" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.7" errorPlus="0.10000000000000009" errorMinus="0.09999999999999998"/> + <measurement value="3.08" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.9" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="17.52" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y04" dimension="2" path="/REF/CDF_2005_S6080774" title="D(SIG(C=PYTHIA))/D(PHI(P=3)-PHI(P=4)) IN PB/RAD" > + <dataPoint> + <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/> + <measurement value="0.02" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.30000000000000004" errorPlus="0.09999999999999998" errorMinus="0.10000000000000003"/> + <measurement value="0.09" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.5" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="0.44" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.7" errorPlus="0.10000000000000009" errorMinus="0.09999999999999998"/> + <measurement value="1.09" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.9" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/> + <measurement value="10.68" errorPlus="0.0" errorMinus="0.0"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Aug 6 22:05:44 2009 (r1749) +++ trunk/data/refdata/Makefile.am Fri Aug 7 17:29:42 2009 (r1750) @@ -36,6 +36,7 @@ CDF_2001_S4751469.aida \ CDF_2002_S4796047.aida \ CDF_2004_S5839831.aida \ + CDF_2005_S6080774.aida \ CDF_2005_S6217184.aida \ CDF_2006_S6450792.aida \ CDF_2006_S6653332.aida \ Added: trunk/include/Rivet/Analyses/CDF_2005_S6080774.hh ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/include/Rivet/Analyses/CDF_2005_S6080774.hh Fri Aug 7 17:29:42 2009 (r1750) @@ -0,0 +1,51 @@ +// -*- C++ -*- +#ifndef RIVET_CDF_2005_S6080774_HH +#define RIVET_CDF_2005_S6080774_HH + +#include "Rivet/Analysis.hh" + +namespace Rivet { + + + class CDF_2005_S6080774 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + CDF_2005_S6080774(); + + /// Factory method + static Analysis* create() { + return new CDF_2005_S6080774(); + } + //@} + + + public: + + /// @name Analysis methods + //@{ + void init(); + void analyze(const Event& event); + void finalize(); + //@} + + private: + + /// @name Histograms + //@{ + + std::vector<AIDA::IHistogram1D*> _h_m_PP; + std::vector<AIDA::IHistogram1D*> _h_pT_PP; + std::vector<AIDA::IHistogram1D*> _h_dphi_PP; + //@} + + }; + + +} + +#endif + Modified: trunk/include/Rivet/Makefile.am ============================================================================== --- trunk/include/Rivet/Makefile.am Thu Aug 6 22:05:44 2009 (r1749) +++ trunk/include/Rivet/Makefile.am Fri Aug 7 17:29:42 2009 (r1750) @@ -44,6 +44,7 @@ Analyses/CDF_2001_S4751469.hh \ Analyses/CDF_2002_S4796047.hh \ Analyses/CDF_2004_S5839831.hh \ + Analyses/CDF_2005_S6080774.hh \ Analyses/CDF_2005_S6217184.hh \ Analyses/CDF_2006_S6653332.hh \ Analyses/CDF_2007_S7057202.hh \ Added: trunk/src/Analyses/CDF_2005_S6080774.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_2005_S6080774.cc Fri Aug 7 17:29:42 2009 (r1750) @@ -0,0 +1,92 @@ +// -*- C++ -*- +#include "Rivet/Rivet.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Analyses/CDF_2005_S6080774.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" + +namespace Rivet { + + + CDF_2005_S6080774::CDF_2005_S6080774() : Analysis("CDF_2005_S6080774") { + setBeams(PROTON, ANTIPROTON); + setNeedsCrossSection(true); + + FinalState fs; + addProjection(fs, "FS"); + + IdentifiedFinalState ifs(-0.9, 0.9, 13.0*GeV); + ifs.acceptId(PHOTON); + addProjection(ifs, "IFS"); + } + + + void CDF_2005_S6080774::init() { + + for (size_t yAxisId=1; yAxisId<5; ++yAxisId) { + _h_m_PP.push_back(bookHistogram1D(1, 1, yAxisId)); + _h_pT_PP.push_back(bookHistogram1D(2, 1, yAxisId)); + _h_dphi_PP.push_back(bookHistogram1D(3, 1, yAxisId)); + } + + } + + + void CDF_2005_S6080774::analyze(const Event& event) { + + const double weight = event.weight(); + + ParticleVector photons = applyProjection<IdentifiedFinalState>(event, "IFS").particles(); + + if (photons.size() < 2 || + (photons[0].momentum().pT() < 14.0*GeV && photons[1].momentum().pT() < 14.0*GeV)) { + vetoEvent; + } + + // isolate photons with ET_sum in cone + ParticleVector isolated_photons; + ParticleVector fs = applyProjection<FinalState>(event, "FS").particles(); + foreach (const Particle& photon, photons) { + FourMomentum mom_in_cone; + double eta_P = photon.momentum().eta(); + double phi_P = photon.momentum().phi(); + foreach (const Particle& p, fs) { + if (deltaR(eta_P, phi_P, p.momentum().eta(), p.momentum().phi()) < 0.4) { + mom_in_cone += p.momentum(); + } + } + if (mom_in_cone.Et()-photon.momentum().Et() < 1.0*GeV) { + isolated_photons.push_back(photon); + } + } + + if (isolated_photons.size()!=2) { + vetoEvent; + } + + FourMomentum mom_PP = isolated_photons[0].momentum() + isolated_photons[1].momentum(); + + for (size_t i=0; i<4; ++i) { + _h_m_PP[i]->fill(mom_PP.mass(), weight); + _h_pT_PP[i]->fill(mom_PP.pT(), weight); + _h_dphi_PP[i]->fill(mapAngle0ToPi(isolated_photons[0].momentum().phi()- + isolated_photons[1].momentum().phi())/M_PI, + weight); + } + + } + + + void CDF_2005_S6080774::finalize() { + + for (size_t i=0; i<4; ++i) { + scale(_h_m_PP[i], crossSection()/sumOfWeights()); + scale(_h_pT_PP[i], crossSection()/sumOfWeights()); + scale(_h_dphi_PP[i], crossSection()/sumOfWeights()); + } + + } + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Aug 6 22:05:44 2009 (r1749) +++ trunk/src/Analyses/Makefile.am Fri Aug 7 17:29:42 2009 (r1750) @@ -12,6 +12,7 @@ CDF_2001_S4751469.cc \ CDF_2002_S4796047.cc \ CDF_2004_S5839831.cc \ + CDF_2005_S6080774.cc \ CDF_2005_S6217184.cc \ CDF_2006_S6450792.cc \ CDF_2006_S6653332.cc \ Modified: trunk/src/Analyses/StdAnalyses.cc ============================================================================== --- trunk/src/Analyses/StdAnalyses.cc Thu Aug 6 22:05:44 2009 (r1749) +++ trunk/src/Analyses/StdAnalyses.cc Fri Aug 7 17:29:42 2009 (r1750) @@ -33,6 +33,7 @@ #include "Rivet/Analyses/CDF_2001_S4751469.hh" #include "Rivet/Analyses/CDF_2002_S4796047.hh" #include "Rivet/Analyses/CDF_2004_S5839831.hh" +#include "Rivet/Analyses/CDF_2005_S6080774.hh" #include "Rivet/Analyses/CDF_2005_S6217184.hh" #include "Rivet/Analyses/CDF_2006_S6450792.hh" #include "Rivet/Analyses/CDF_2006_S6653332.hh" @@ -124,6 +125,7 @@ fns["CDF_2001_S4751469"] = Rivet::CDF_2001_S4751469::create; fns["CDF_2002_S4796047"] = Rivet::CDF_2002_S4796047::create; fns["CDF_2004_S5839831"] = Rivet::CDF_2004_S5839831::create; + fns["CDF_2005_S6080774"] = Rivet::CDF_2005_S6080774::create; fns["CDF_2005_S6217184"] = Rivet::CDF_2005_S6217184::create; fns["CDF_2006_S6450792"] = Rivet::CDF_2006_S6450792::create; fns["CDF_2006_S6653332"] = Rivet::CDF_2006_S6653332::create;
More information about the Rivet-svn mailing list |