|
[Rivet-svn] r2312 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Mar 4 17:02:03 GMT 2010
Author: fsiegert Date: Thu Mar 4 17:02:02 2010 New Revision: 2312 Log: Add untested first version of D0_2010_S8570965. Added: trunk/data/anainfo/D0_2010_S8570965.info trunk/data/plotinfo/D0_2010_S8570965.plot trunk/data/refdata/D0_2010_S8570965.aida trunk/src/Analyses/D0_2010_S8570965.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/D0_2010_S8570965.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/D0_2010_S8570965.info Thu Mar 4 17:02:02 2010 (r2312) @@ -0,0 +1,25 @@ +Name: D0_2010_S8570965 +Year: 2010 +Summary: Direct photon pair production +Experiment: CDF +Collider: Tevatron Run 2 +SpiresID: 8570965 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - arXiv:1002.4917 +RunInfo: + All processes that can produce prompt photon pairs, e.g. $jj \to jj$, + $jj \to j\gamma$ and $jj \to \gamma \gamma$. + Non-prompt photons from hadron decays like $\pi$ and $\eta$ have been + corrected for. +NumEvents: 10000000 +PtCuts: [0] +Description: + Direct photon pair production cross sections are measured using 4.2 fb$^{-1}$ + of data. They are binned in diphoton mass, the transverse momentum of the + diphoton system, the azimuthal angle between the photons, and the polar + scattering angle of the photons. Also available are double differential cross + sections considering the last three kinematic variables in three diphoton + mass bins. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Mar 4 14:36:14 2010 (r2311) +++ trunk/data/anainfo/Makefile.am Thu Mar 4 17:02:02 2010 (r2312) @@ -52,6 +52,7 @@ D0_2009_S8320160.info \ D0_2009_S8349509.info \ D0_2010_S8566488.info \ + D0_2010_S8570965.info \ E735_1998_S3905616.info \ DELPHI_1995_S3137023.info \ DELPHI_1996_S3430090.info \ Added: trunk/data/plotinfo/D0_2010_S8570965.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/D0_2010_S8570965.plot Thu Mar 4 17:02:02 2010 (r2312) @@ -0,0 +1,83 @@ +# BEGIN PLOT /D0_2010_S8570965/d01-x01-y01 +Title=Diphoton mass +XLabel=$M_{\gamma\gamma}$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}M_{\gamma\gamma}$ [pb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d02-x01-y01 +Title=$p_\perp$ of the diphoton system +XLabel=$p_\perp^{\gamma\gamma}$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_\perp^{\gamma\gamma}$ [pb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d03-x01-y01 +Title=Azimuthal angle between the photons +XLabel=$\Delta\phi_{\gamma\gamma}$ [rad] +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}$ [pb/rad] +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d04-x01-y01 +Title=Polar scattering angle of the photons +XLabel=$|\cos\theta^*|$ +YLabel=$\mathrm{d}\sigma/\mathrm{d}|\cos\theta^*|$ [pb] +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d02-x01-y01 +Title=$p_\perp$ of the diphoton system (30 GeV$<M_{\gamma\gamma}<$50 GeV) +XLabel=$p_\perp^{\gamma\gamma}$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_\perp^{\gamma\gamma}$ [pb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d03-x01-y01 +Title=Azimuthal angle between the photons (30 GeV$<M_{\gamma\gamma}<$50 GeV) +XLabel=$\Delta\phi_{\gamma\gamma}$ [rad] +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}$ [pb/rad] +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d04-x01-y01 +Title=Polar scattering angle of the photons (30 GeV$<M_{\gamma\gamma}<$50 GeV) +XLabel=$|\cos\theta^*|$ +YLabel=$\mathrm{d}\sigma/\mathrm{d}|\cos\theta^*|$ [pb] +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d02-x01-y01 +Title=$p_\perp$ of the diphoton system (50 GeV$<M_{\gamma\gamma}<$80 GeV) +XLabel=$p_\perp^{\gamma\gamma}$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_\perp^{\gamma\gamma}$ [pb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d03-x01-y01 +Title=Azimuthal angle between the photons (50 GeV$<M_{\gamma\gamma}<$80 GeV) +XLabel=$\Delta\phi_{\gamma\gamma}$ [rad] +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}$ [pb/rad] +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d04-x01-y01 +Title=Polar scattering angle of the photons (50 GeV$<M_{\gamma\gamma}<$80 GeV) +XLabel=$|\cos\theta^*|$ +YLabel=$\mathrm{d}\sigma/\mathrm{d}|\cos\theta^*|$ [pb] +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d02-x01-y01 +Title=$p_\perp$ of the diphoton system (80 GeV$<M_{\gamma\gamma}<$350 GeV) +XLabel=$p_\perp^{\gamma\gamma}$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}p_\perp^{\gamma\gamma}$ [pb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d03-x01-y01 +Title=Azimuthal angle between the photons (80 GeV$<M_{\gamma\gamma}<$350 GeV) +XLabel=$\Delta\phi_{\gamma\gamma}$ [rad] +YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}$ [pb/rad] +# END PLOT + +# BEGIN PLOT /D0_2010_S8570965/d04-x01-y01 +Title=Polar scattering angle of the photons (80 GeV$<M_{\gamma\gamma}<$350 GeV) +XLabel=$|\cos\theta^*|$ +YLabel=$\mathrm{d}\sigma/\mathrm{d}|\cos\theta^*|$ [pb] +# END PLOT + Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Mar 4 14:36:14 2010 (r2311) +++ trunk/data/plotinfo/Makefile.am Thu Mar 4 17:02:02 2010 (r2312) @@ -51,6 +51,7 @@ D0_2009_S8320160.plot \ D0_2009_S8349509.plot \ D0_2010_S8566488.plot \ + D0_2010_S8570965.plot \ E735_1998_S3905616.plot \ DELPHI_1995_S3137023.plot \ DELPHI_1996_S3430090.plot \ Added: trunk/data/refdata/D0_2010_S8570965.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/D0_2010_S8570965.aida Thu Mar 4 17:02:02 2010 (r2312) @@ -0,0 +1,331 @@ +<?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/D0_2010_S8570965" title="D(SIG)/DM(P=3_4) IN PB/GEV" > + <dataPoint> + <measurement value="43.0" errorPlus="2.0" errorMinus="13.0"/> + <measurement value="0.0311" errorPlus="0.010154042840169623" errorMinus="0.010154042840169623"/> + </dataPoint> + <dataPoint> + <measurement value="47.6" errorPlus="2.3999999999999986" errorMinus="2.6000000000000014"/> + <measurement value="0.174" errorPlus="0.03820082721617426" errorMinus="0.03820082721617426"/> + </dataPoint> + <dataPoint> + <measurement value="54.7" errorPlus="5.299999999999997" errorMinus="4.700000000000003"/> + <measurement value="0.119" errorPlus="0.02347046867874606" errorMinus="0.02347046867874606"/> + </dataPoint> + <dataPoint> + <measurement value="64.6" errorPlus="5.400000000000006" errorMinus="4.599999999999994"/> + <measurement value="0.0789" errorPlus="0.015319608904929654" errorMinus="0.015319608904929654"/> + </dataPoint> + <dataPoint> + <measurement value="74.6" errorPlus="5.400000000000006" errorMinus="4.599999999999994"/> + <measurement value="0.0561" errorPlus="0.010113571327676488" errorMinus="0.010113571327676488"/> + </dataPoint> + <dataPoint> + <measurement value="88.6" errorPlus="11.400000000000006" errorMinus="8.599999999999994"/> + <measurement value="0.0239" errorPlus="0.004591040078239353" errorMinus="0.004591040078239353"/> + </dataPoint> + <dataPoint> + <measurement value="108.9" errorPlus="11.099999999999994" errorMinus="8.900000000000006"/> + <measurement value="0.0112" errorPlus="0.002298047867212517" errorMinus="0.002298047867212517"/> + </dataPoint> + <dataPoint> + <measurement value="132.9" errorPlus="17.099999999999994" errorMinus="12.900000000000006"/> + <measurement value="0.00365" errorPlus="9.82792577302047E-4" errorMinus="9.82792577302047E-4"/> + </dataPoint> + <dataPoint> + <measurement value="170.7" errorPlus="29.30000000000001" errorMinus="20.69999999999999"/> + <measurement value="0.00167" errorPlus="4.076989575655057E-4" errorMinus="4.076989575655057E-4"/> + </dataPoint> + <dataPoint> + <measurement value="248.8" errorPlus="101.19999999999999" errorMinus="48.80000000000001"/> + <measurement value="3.3E-4" errorPlus="9.744783219754044E-5" errorMinus="9.744783219754044E-5"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D(SIG)/DPT(P=3_4) IN PB/GEV" > + <dataPoint> + <measurement value="1.5" errorPlus="1.0" errorMinus="1.5"/> + <measurement value="0.192" errorPlus="0.04647827879773518" errorMinus="0.04647827879773518"/> + </dataPoint> + <dataPoint> + <measurement value="3.7" errorPlus="1.2999999999999998" errorMinus="1.2000000000000002"/> + <measurement value="0.334" errorPlus="0.0676298454825974" errorMinus="0.0676298454825974"/> + </dataPoint> + <dataPoint> + <measurement value="6.2" errorPlus="1.2999999999999998" errorMinus="1.2000000000000002"/> + <measurement value="0.306" errorPlus="0.05941445278717965" errorMinus="0.05941445278717965"/> + </dataPoint> + <dataPoint> + <measurement value="8.7" errorPlus="1.3000000000000007" errorMinus="1.1999999999999993"/> + <measurement value="0.238" errorPlus="0.049524591871109844" errorMinus="0.049524591871109844"/> + </dataPoint> + <dataPoint> + <measurement value="11.2" errorPlus="1.3000000000000007" errorMinus="1.1999999999999993"/> + <measurement value="0.166" errorPlus="0.03529208409827903" errorMinus="0.03529208409827903"/> + </dataPoint> + <dataPoint> + <measurement value="13.7" errorPlus="1.3000000000000007" errorMinus="1.1999999999999993"/> + <measurement value="0.11" errorPlus="0.02804460732476032" errorMinus="0.02804460732476032"/> + </dataPoint> + <dataPoint> + <measurement value="17.3" errorPlus="2.6999999999999993" errorMinus="2.3000000000000007"/> + <measurement value="0.088" errorPlus="0.019950979925808153" errorMinus="0.019950979925808153"/> + </dataPoint> + <dataPoint> + <measurement value="22.3" errorPlus="2.6999999999999993" errorMinus="2.3000000000000007"/> + <measurement value="0.063" errorPlus="0.015172409169278294" errorMinus="0.015172409169278294"/> + </dataPoint> + <dataPoint> + <measurement value="27.3" errorPlus="2.6999999999999993" errorMinus="2.3000000000000007"/> + <measurement value="0.042" errorPlus="0.010992451955774018" errorMinus="0.010992451955774018"/> + </dataPoint> + <dataPoint> + <measurement value="34.3" errorPlus="5.700000000000003" errorMinus="4.299999999999997"/> + <measurement value="0.0299" errorPlus="0.006398879433150776" errorMinus="0.006398879433150776"/> + </dataPoint> + <dataPoint> + <measurement value="47.8" errorPlus="12.200000000000003" errorMinus="7.799999999999997"/> + <measurement value="0.00758" errorPlus="0.0019414272687896398" errorMinus="0.0019414272687896398"/> + </dataPoint> + <dataPoint> + <measurement value="73.4" errorPlus="26.599999999999994" errorMinus="13.400000000000006"/> + <measurement value="9.92E-4" errorPlus="4.134391331260262E-4" errorMinus="4.134391331260262E-4"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D(SIG)/DABS(PHI(P=3)-PHI(P=4))) IN PB/RAD" > + <dataPoint> + <measurement value="1.75" errorPlus="0.1299999999999999" errorMinus="0.17999999999999994"/> + <measurement value="0.432" errorPlus="0.12528" errorMinus="0.12528"/> + </dataPoint> + <dataPoint> + <measurement value="2.06" errorPlus="0.14000000000000012" errorMinus="0.18000000000000016"/> + <measurement value="0.53" errorPlus="0.15287537407967314" errorMinus="0.15287537407967314"/> + </dataPoint> + <dataPoint> + <measurement value="2.38" errorPlus="0.1299999999999999" errorMinus="0.17999999999999972"/> + <measurement value="1.15" errorPlus="0.26021529547664946" errorMinus="0.26021529547664946"/> + </dataPoint> + <dataPoint> + <measurement value="2.6" errorPlus="0.06999999999999984" errorMinus="0.0900000000000003"/> + <measurement value="2.43" errorPlus="0.5735005928506091" errorMinus="0.5735005928506091"/> + </dataPoint> + <dataPoint> + <measurement value="2.76" errorPlus="0.07000000000000028" errorMinus="0.08999999999999986"/> + <measurement value="3.99" errorPlus="0.7747178647740094" errorMinus="0.7747178647740094"/> + </dataPoint> + <dataPoint> + <measurement value="2.92" errorPlus="0.06000000000000005" errorMinus="0.08999999999999986"/> + <measurement value="6.7" errorPlus="1.2641534716955851" errorMinus="1.2641534716955851"/> + </dataPoint> + <dataPoint> + <measurement value="3.08" errorPlus="0.06000000000000005" errorMinus="0.10000000000000009"/> + <measurement value="13.4" errorPlus="2.3402093923407796" errorMinus="2.3402093923407796"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D(SIG)/DABS(COS(THETA)) IN PB" > + <dataPoint> + <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/> + <measurement value="13.8" errorPlus="2.5927846034717192" errorMinus="2.5927846034717192"/> + </dataPoint> + <dataPoint> + <measurement value="0.15" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/> + <measurement value="10.0" errorPlus="1.8357559750685821" errorMinus="1.8357559750685821"/> + </dataPoint> + <dataPoint> + <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="7.78" errorPlus="1.4679274641480076" errorMinus="1.4679274641480076"/> + </dataPoint> + <dataPoint> + <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/> + <measurement value="6.38" errorPlus="1.276" errorMinus="1.276"/> + </dataPoint> + <dataPoint> + <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="4.77" errorPlus="1.0141159105348854" errorMinus="1.0141159105348854"/> + </dataPoint> + <dataPoint> + <measurement value="0.57" errorPlus="0.13" errorMinus="0.06999999999999995"/> + <measurement value="2.35" errorPlus="0.5153952366873408" errorMinus="0.5153952366873408"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D2(SIG)/DM(P=3_4)/PT(P=3_4) IN PB/GEV**2" > + <dataPoint> + <measurement value="2.4" errorPlus="2.6" errorMinus="2.4"/> + <measurement value="0.102" errorPlus="0.020928650219256852" errorMinus="0.020928650219256852"/> + </dataPoint> + <dataPoint> + <measurement value="7.0" errorPlus="3.0" errorMinus="2.0"/> + <measurement value="0.0731" errorPlus="0.016669364714949397" errorMinus="0.016669364714949397"/> + </dataPoint> + <dataPoint> + <measurement value="12.2" errorPlus="2.8000000000000007" errorMinus="2.1999999999999993"/> + <measurement value="0.0433" errorPlus="0.01021916694256435" errorMinus="0.01021916694256435"/> + </dataPoint> + <dataPoint> + <measurement value="23.4" errorPlus="26.6" errorMinus="8.399999999999999"/> + <measurement value="0.00716" errorPlus="0.0016898206768766915" errorMinus="0.0016898206768766915"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D(2SIG)/DM(P=3_4)/DABS(PHI(P=3)-PHI(P=4))) IN PB/GEV/RAD" > + <dataPoint> + <measurement value="2.16" errorPlus="0.34999999999999964" errorMinus="0.5900000000000001"/> + <measurement value="0.295" errorPlus="0.06727035008084914" errorMinus="0.06727035008084914"/> + </dataPoint> + <dataPoint> + <measurement value="2.7" errorPlus="0.1299999999999999" errorMinus="0.1900000000000004"/> + <measurement value="0.907" errorPlus="0.19974602999809535" errorMinus="0.19974602999809535"/> + </dataPoint> + <dataPoint> + <measurement value="2.92" errorPlus="0.06000000000000005" errorMinus="0.08999999999999986"/> + <measurement value="1.89" errorPlus="0.4460560166615848" errorMinus="0.4460560166615848"/> + </dataPoint> + <dataPoint> + <measurement value="3.08" errorPlus="0.06000000000000005" errorMinus="0.10000000000000009"/> + <measurement value="3.13" errorPlus="0.6654471278771891" errorMinus="0.6654471278771891"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d07-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D2(SIG)/DM(P=3_4)/DABS(COS(THETA)) IN PB/GEV" > + <dataPoint> + <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/> + <measurement value="5.1" errorPlus="1.0464325109628427" errorMinus="1.0464325109628427"/> + </dataPoint> + <dataPoint> + <measurement value="0.15" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/> + <measurement value="4.18" errorPlus="0.8867119036079306" errorMinus="0.8867119036079306"/> + </dataPoint> + <dataPoint> + <measurement value="0.28" errorPlus="0.12" errorMinus="0.08000000000000002"/> + <measurement value="1.77" errorPlus="0.43965887913244744" errorMinus="0.43965887913244744"/> + </dataPoint> + <dataPoint> + <measurement value="0.44" errorPlus="0.25999999999999995" errorMinus="0.03999999999999998"/> + <measurement value="0.36" errorPlus="0.13708391590555036" errorMinus="0.13708391590555036"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d08-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D2(SIG)/DM(P=3_4)/PT(P=3_4) IN PB/GEV**2" > + <dataPoint> + <measurement value="2.8" errorPlus="2.2" errorMinus="2.8"/> + <measurement value="0.11" errorPlus="0.02257011298155151" errorMinus="0.02257011298155151"/> + </dataPoint> + <dataPoint> + <measurement value="7.3" errorPlus="2.7" errorMinus="2.3"/> + <measurement value="0.148" errorPlus="0.027289851593586947" errorMinus="0.027289851593586947"/> + </dataPoint> + <dataPoint> + <measurement value="12.3" errorPlus="2.6999999999999993" errorMinus="2.3000000000000007"/> + <measurement value="0.088" errorPlus="0.01742311108843653" errorMinus="0.01742311108843653"/> + </dataPoint> + <dataPoint> + <measurement value="17.3" errorPlus="2.6999999999999993" errorMinus="2.3000000000000007"/> + <measurement value="0.0557" errorPlus="0.012701554235604397" errorMinus="0.012701554235604397"/> + </dataPoint> + <dataPoint> + <measurement value="24.1" errorPlus="5.899999999999999" errorMinus="4.100000000000001"/> + <measurement value="0.0247" errorPlus="0.005632466599989741" errorMinus="0.005632466599989741"/> + </dataPoint> + <dataPoint> + <measurement value="39.8" errorPlus="40.2" errorMinus="9.799999999999997"/> + <measurement value="0.00403" errorPlus="8.87515436485473E-4" errorMinus="8.87515436485473E-4"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d09-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D(2SIG)/DM(P=3_4)/DABS(PHI(P=3)-PHI(P=4))) IN PB/GEV/RAD" > + <dataPoint> + <measurement value="1.98" errorPlus="0.2200000000000002" errorMinus="0.4099999999999999"/> + <measurement value="0.186" errorPlus="0.05329476146864718" errorMinus="0.05329476146864718"/> + </dataPoint> + <dataPoint> + <measurement value="2.38" errorPlus="0.1299999999999999" errorMinus="0.17999999999999972"/> + <measurement value="0.582" errorPlus="0.14208430736714028" errorMinus="0.14208430736714028"/> + </dataPoint> + <dataPoint> + <measurement value="2.6" errorPlus="0.06999999999999984" errorMinus="0.0900000000000003"/> + <measurement value="1.35" errorPlus="0.31861144047256057" errorMinus="0.31861144047256057"/> + </dataPoint> + <dataPoint> + <measurement value="2.76" errorPlus="0.07000000000000028" errorMinus="0.08999999999999986"/> + <measurement value="1.99" errorPlus="0.4230798033468391" errorMinus="0.4230798033468391"/> + </dataPoint> + <dataPoint> + <measurement value="2.92" errorPlus="0.06000000000000005" errorMinus="0.08999999999999986"/> + <measurement value="3.55" errorPlus="0.7028641404994281" errorMinus="0.7028641404994281"/> + </dataPoint> + <dataPoint> + <measurement value="3.07" errorPlus="0.07000000000000028" errorMinus="0.08999999999999986"/> + <measurement value="6.89" errorPlus="1.185400421798474" errorMinus="1.185400421798474"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d10-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D2(SIG)/DM(P=3_4)/DABS(COS(THETA)) IN PB/GEV" > + <dataPoint> + <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/> + <measurement value="5.31" errorPlus="1.0144740755682224" errorMinus="1.0144740755682224"/> + </dataPoint> + <dataPoint> + <measurement value="0.15" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/> + <measurement value="4.49" errorPlus="0.8889746453077276" errorMinus="0.8889746453077276"/> + </dataPoint> + <dataPoint> + <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="4.58" errorPlus="0.8750077713940602" errorMinus="0.8750077713940602"/> + </dataPoint> + <dataPoint> + <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/> + <measurement value="3.44" errorPlus="0.7313540319161439" errorMinus="0.7313540319161439"/> + </dataPoint> + <dataPoint> + <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="3.17" errorPlus="0.6981200827937841" errorMinus="0.6981200827937841"/> + </dataPoint> + <dataPoint> + <measurement value="0.58" errorPlus="0.12" errorMinus="0.07999999999999996"/> + <measurement value="1.52" errorPlus="0.37107929071830453" errorMinus="0.37107929071830453"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d11-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D2(SIG)/DM(P=3_4)/PT(P=3_4) IN PB/GEV**2" > + <dataPoint> + <measurement value="2.8" errorPlus="2.2" errorMinus="2.8"/> + <measurement value="0.0443" errorPlus="0.013029020876489529" errorMinus="0.013029020876489529"/> + </dataPoint> + <dataPoint> + <measurement value="9.3" errorPlus="5.699999999999999" errorMinus="4.300000000000001"/> + <measurement value="0.0275" errorPlus="0.0056425282453878775" errorMinus="0.0056425282453878775"/> + </dataPoint> + <dataPoint> + <measurement value="24.3" errorPlus="15.7" errorMinus="9.3"/> + <measurement value="0.012" errorPlus="0.0024738633753705963" errorMinus="0.0024738633753705963"/> + </dataPoint> + <dataPoint> + <measurement value="58.1" errorPlus="41.9" errorMinus="18.1"/> + <measurement value="0.0018" errorPlus="4.5429946070846267E-4" errorMinus="4.5429946070846267E-4"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d12-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D(2SIG)/DM(P=3_4)/DABS(PHI(P=3)-PHI(P=4))) IN PB/GEV/RAD" > + <dataPoint> + <measurement value="2.42" errorPlus="0.25" errorMinus="0.8499999999999999"/> + <measurement value="0.0979" errorPlus="0.026631679256104" errorMinus="0.026631679256104"/> + </dataPoint> + <dataPoint> + <measurement value="2.87" errorPlus="0.10999999999999988" errorMinus="0.20000000000000018"/> + <measurement value="0.928" errorPlus="0.1837346260235125" errorMinus="0.1837346260235125"/> + </dataPoint> + <dataPoint> + <measurement value="3.08" errorPlus="0.06000000000000005" errorMinus="0.10000000000000009"/> + <measurement value="3.21" errorPlus="0.5715242514539518" errorMinus="0.5715242514539518"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d13-x01-y01" dimension="2" path="/REF/D0_2010_S8570965" title="D2(SIG)/DM(P=3_4)/DABS(COS(THETA)) IN PB/GEV" > + <dataPoint> + <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/> + <measurement value="2.05" errorPlus="0.3780013227490084" errorMinus="0.3780013227490084"/> + </dataPoint> + <dataPoint> + <measurement value="0.3" errorPlus="0.10000000000000003" errorMinus="0.09999999999999998"/> + <measurement value="1.38" errorPlus="0.26364862980869064" errorMinus="0.26364862980869064"/> + </dataPoint> + <dataPoint> + <measurement value="0.53" errorPlus="0.16999999999999993" errorMinus="0.13"/> + <measurement value="0.762" errorPlus="0.17983845751117863" errorMinus="0.17983845751117863"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Mar 4 14:36:14 2010 (r2311) +++ trunk/data/refdata/Makefile.am Thu Mar 4 17:02:02 2010 (r2312) @@ -82,6 +82,7 @@ D0_2009_S8320160.aida \ D0_2009_S8349509.aida \ D0_2010_S8566488.aida \ + D0_2010_S8570965.aida \ E735_1992_S2485869.aida \ E735_1993_S2896508.aida \ E735_1998_S3905616.aida \ Added: trunk/src/Analyses/D0_2010_S8570965.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/D0_2010_S8570965.cc Thu Mar 4 17:02:02 2010 (r2312) @@ -0,0 +1,150 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" +#include "Rivet/Tools/ParticleIdUtils.hh" +#include "Rivet/Tools/BinnedHistogram.hh" + +namespace Rivet { + + typedef std::pair<double, double> doublepair; + + class D0_2010_S8570965 : public Analysis { + public: + + D0_2010_S8570965() + : Analysis("D0_2010_S8570965") + { + setBeams(PROTON, ANTIPROTON); + setNeedsCrossSection(true); + } + + + public: + + void init() { + FinalState fs; + addProjection(fs, "FS"); + + IdentifiedFinalState ifs(-0.9, 0.9, 20.0*GeV); + ifs.acceptId(PHOTON); + addProjection(ifs, "IFS"); + + _h_M = bookHistogram1D(1, 1, 1); + _h_pT = bookHistogram1D(2, 1, 1); + _h_dPhi = bookHistogram1D(3, 1, 1); + _h_costheta = bookHistogram1D(4, 1, 1); + + std::vector<doublepair> M_ranges; + M_ranges.push_back(std::make_pair(30.0, 50.0)); + M_ranges.push_back(std::make_pair(50.0, 80.0)); + M_ranges.push_back(std::make_pair(80.0, 350.0)); + int i=0; + + foreach (const doublepair& M_range, M_ranges) { + _h_pT_M.addHistogram(M_range.first, M_range.second, bookHistogram1D(5+i, 1, 1)); + _h_dPhi_M.addHistogram(M_range.first, M_range.second, bookHistogram1D(6+i, 1, 1)); + _h_costheta_M.addHistogram(M_range.first, M_range.second, bookHistogram1D(7+i, 1, 1)); + ++i; + } + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + ParticleVector photons = applyProjection<IdentifiedFinalState>(event, "IFS").particlesByPt(); + if (photons.size() < 2 || + (photons[0].momentum().pT() < 21.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(); + double Etsum=0.0; + 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 (PID::threeCharge(p.pdgId())!=0) Etsum += p.momentum().Et(); + } + } + if (mom_in_cone.E()/photon.momentum().E() < 1.1 && Etsum<1.5*GeV) { + isolated_photons.push_back(photon); + } + } + + if (isolated_photons.size() != 2) { + vetoEvent; + } + + FourMomentum y1=isolated_photons[0].momentum(); + FourMomentum y2=isolated_photons[0].momentum(); + if (deltaR(y1, y2)<0.4) { + vetoEvent; + } + + FourMomentum yy=y1+y2; + double Myy = yy.mass()/GeV; + double pTyy = yy.pT()/GeV; + if (Myy<pTyy) { + vetoEvent; + } + + double dPhiyy = mapAngle0ToPi(y1.phi()-y2.phi()); + double costhetayy = fabs(tanh(y1.eta()-y2.eta())/2.0); + + _h_M->fill(Myy, weight); + _h_pT->fill(pTyy, weight); + _h_dPhi->fill(dPhiyy, weight); + _h_costheta->fill(costhetayy, weight); + + _h_pT_M.fill(Myy, pTyy, weight); + _h_dPhi_M.fill(Myy, dPhiyy, weight); + _h_costheta_M.fill(Myy, costhetayy, weight); + } + + + void finalize() { + + scale(_h_M, crossSection()/sumOfWeights()); + scale(_h_pT, crossSection()/sumOfWeights()); + scale(_h_dPhi, crossSection()/sumOfWeights()); + scale(_h_costheta, crossSection()/sumOfWeights()); + std::vector<double> dMyy; + dMyy += 20.0, 30.0, 270.0; + for (size_t i=0; i<3; ++i) { + scale(_h_pT_M.getHistograms()[i], crossSection()/sumOfWeights()/dMyy[i]); + scale(_h_dPhi_M.getHistograms()[i], crossSection()/sumOfWeights()/dMyy[i]); + scale(_h_costheta_M.getHistograms()[i], crossSection()/sumOfWeights()/dMyy[i]); + } + + } + + + private: + + AIDA::IHistogram1D *_h_M; + AIDA::IHistogram1D *_h_pT; + AIDA::IHistogram1D *_h_dPhi; + AIDA::IHistogram1D *_h_costheta; + BinnedHistogram<double> _h_pT_M; + BinnedHistogram<double> _h_dPhi_M; + BinnedHistogram<double> _h_costheta_M; + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<D0_2010_S8570965> plugin_D0_2010_S8570965; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Mar 4 14:36:14 2010 (r2311) +++ trunk/src/Analyses/Makefile.am Thu Mar 4 17:02:02 2010 (r2312) @@ -25,6 +25,7 @@ D0_1996_S3214044.cc \ D0_1996_S3324664.cc \ D0_1998_S3711838.cc \ + D0_2010_S8570965.cc \ DELPHI_2003_WUD_03_11.cc \ E735_1998_S3905616.cc \ H1_1995_S3167097.cc \
More information about the Rivet-svn mailing list |