|
|
[Rivet-svn] r2402 - in trunk: data/anainfo data/plotinfo data/refdata src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Mon Apr 19 12:28:07 BST 2010
Author: fsiegert
Date: Mon Apr 19 12:28:06 2010
New Revision: 2402
Log:
Update to D0_2010 diphoton analysis, because version 2 of the preprint fixes
a problem with the dM binwidth in the normalisation of the double differential
cross sections.
Also I have implemented exactly the cuts that they apply to the theoretical
predictions (see caption of Fig. 4), since it seems that's what the data has
been corrected to. Doesn't make a whole lot of a difference, but is neater
and more correct.
Modified:
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/D0_2010_S8570965.info
==============================================================================
--- trunk/data/anainfo/D0_2010_S8570965.info Tue Apr 13 21:04:54 2010 (r2401)
+++ trunk/data/anainfo/D0_2010_S8570965.info Mon Apr 19 12:28:06 2010 (r2402)
@@ -25,6 +25,10 @@
scattering angle of the photons. Also available are double differential cross
sections considering the last three kinematic variables in three diphoton
mass bins.
+ Note: The numbers in version 1 of the arXiv preprint were missing the dM
+ normalisation in the double differential cross sections. This has been
+ reported to and fixed by the authors in v2 and the journal submission.
+ HepData as well as the Rivet analysis have also been updated.
BibKey: Abazov:2010ah
BibTeX: '@Article{Abazov:2010ah,
author = "Abazov, V. M. and others",
Modified: trunk/data/plotinfo/D0_2010_S8570965.plot
==============================================================================
--- trunk/data/plotinfo/D0_2010_S8570965.plot Tue Apr 13 21:04:54 2010 (r2401)
+++ trunk/data/plotinfo/D0_2010_S8570965.plot Mon Apr 19 12:28:06 2010 (r2402)
@@ -28,7 +28,7 @@
# BEGIN PLOT /D0_2010_S8570965/d05-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}p_\perp^{\gamma\gamma}/\mathrm{d}M_{\gamma\gamma}$ [pb/GeV$^2$]
LogX=1
XMin=0.1
# END PLOT
@@ -36,19 +36,19 @@
# BEGIN PLOT /D0_2010_S8570965/d06-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}/\mathrm{d}M_{\gamma\gamma}$ [pb/rad/GeV]
# END PLOT
# BEGIN PLOT /D0_2010_S8570965/d07-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}|\cos\theta^*|/\mathrm{d}M_{\gamma\gamma}$ [pb/GeV]
# END PLOT
# BEGIN PLOT /D0_2010_S8570965/d08-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}p_\perp^{\gamma\gamma}/\mathrm{d}M_{\gamma\gamma}$ [pb/GeV$^2$]
LogX=1
XMin=0.1
# END PLOT
@@ -56,19 +56,19 @@
# BEGIN PLOT /D0_2010_S8570965/d09-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}/\mathrm{d}M_{\gamma\gamma}$ [pb/rad/GeV]
# END PLOT
# BEGIN PLOT /D0_2010_S8570965/d10-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}|\cos\theta^*|/\mathrm{d}M_{\gamma\gamma}$ [pb/GeV]
# END PLOT
# BEGIN PLOT /D0_2010_S8570965/d11-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}p_\perp^{\gamma\gamma}/\mathrm{d}M_{\gamma\gamma}$ [pb/GeV$^2$]
LogX=1
XMin=0.1
# END PLOT
@@ -76,12 +76,12 @@
# BEGIN PLOT /D0_2010_S8570965/d12-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}/\mathrm{d}M_{\gamma\gamma}$ [pb/rad/GeV]
# END PLOT
# BEGIN PLOT /D0_2010_S8570965/d13-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]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}|\cos\theta^*|/\mathrm{d}M_{\gamma\gamma}$ [pb/GeV]
# END PLOT
Modified: trunk/data/refdata/D0_2010_S8570965.aida
==============================================================================
--- trunk/data/refdata/D0_2010_S8570965.aida Tue Apr 13 21:04:54 2010 (r2401)
+++ trunk/data/refdata/D0_2010_S8570965.aida Mon Apr 19 12:28:06 2010 (r2402)
@@ -153,179 +153,179 @@
<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"/>
+ <measurement value="0.00511" errorPlus="0.0010484843394157112" errorMinus="0.0010484843394157112"/>
</dataPoint>
<dataPoint>
<measurement value="7.0" errorPlus="3.0" errorMinus="2.0"/>
- <measurement value="0.0731" errorPlus="0.016669364714949397" errorMinus="0.016669364714949397"/>
+ <measurement value="0.00365" errorPlus="8.323280603223707E-4" errorMinus="8.323280603223707E-4"/>
</dataPoint>
<dataPoint>
<measurement value="12.2" errorPlus="2.8000000000000007" errorMinus="2.1999999999999993"/>
- <measurement value="0.0433" errorPlus="0.01021916694256435" errorMinus="0.01021916694256435"/>
+ <measurement value="0.00217" errorPlus="5.121383895003381E-4" errorMinus="5.121383895003381E-4"/>
</dataPoint>
<dataPoint>
<measurement value="23.4" errorPlus="26.6" errorMinus="8.399999999999999"/>
- <measurement value="0.00716" errorPlus="0.0016898206768766915" errorMinus="0.0016898206768766915"/>
+ <measurement value="3.58E-4" errorPlus="8.449103384383457E-5" errorMinus="8.449103384383457E-5"/>
</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"/>
+ <measurement value="0.0148" errorPlus="0.003374919258293449" errorMinus="0.003374919258293449"/>
</dataPoint>
<dataPoint>
<measurement value="2.7" errorPlus="0.1299999999999999" errorMinus="0.1900000000000004"/>
- <measurement value="0.907" errorPlus="0.19974602999809535" errorMinus="0.19974602999809535"/>
+ <measurement value="0.0454" errorPlus="0.00999831285767754" errorMinus="0.00999831285767754"/>
</dataPoint>
<dataPoint>
<measurement value="2.92" errorPlus="0.06000000000000005" errorMinus="0.08999999999999986"/>
- <measurement value="1.89" errorPlus="0.4460560166615848" errorMinus="0.4460560166615848"/>
+ <measurement value="0.0945" errorPlus="0.02230280083307924" errorMinus="0.02230280083307924"/>
</dataPoint>
<dataPoint>
<measurement value="3.08" errorPlus="0.06000000000000005" errorMinus="0.10000000000000009"/>
- <measurement value="3.13" errorPlus="0.6654471278771891" errorMinus="0.6654471278771891"/>
+ <measurement value="0.157" errorPlus="0.0333786578519868" errorMinus="0.0333786578519868"/>
</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"/>
+ <measurement value="0.255" errorPlus="0.05232162554814214" errorMinus="0.05232162554814214"/>
</dataPoint>
<dataPoint>
<measurement value="0.15" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/>
- <measurement value="4.18" errorPlus="0.8867119036079306" errorMinus="0.8867119036079306"/>
+ <measurement value="0.209" errorPlus="0.04433559518039652" errorMinus="0.04433559518039652"/>
</dataPoint>
<dataPoint>
<measurement value="0.28" errorPlus="0.12" errorMinus="0.08000000000000002"/>
- <measurement value="1.77" errorPlus="0.43965887913244744" errorMinus="0.43965887913244744"/>
+ <measurement value="0.0884" errorPlus="0.021958104471925626" errorMinus="0.021958104471925626"/>
</dataPoint>
<dataPoint>
<measurement value="0.44" errorPlus="0.25999999999999995" errorMinus="0.03999999999999998"/>
- <measurement value="0.36" errorPlus="0.13708391590555036" errorMinus="0.13708391590555036"/>
+ <measurement value="0.018" errorPlus="0.006854195795277516" errorMinus="0.006854195795277516"/>
</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"/>
+ <measurement value="0.00368" errorPlus="7.550728706555413E-4" errorMinus="7.550728706555413E-4"/>
</dataPoint>
<dataPoint>
<measurement value="7.3" errorPlus="2.7" errorMinus="2.3"/>
- <measurement value="0.148" errorPlus="0.027289851593586947" errorMinus="0.027289851593586947"/>
+ <measurement value="0.00492" errorPlus="9.072031745976201E-4" errorMinus="9.072031745976201E-4"/>
</dataPoint>
<dataPoint>
<measurement value="12.3" errorPlus="2.6999999999999993" errorMinus="2.3000000000000007"/>
- <measurement value="0.088" errorPlus="0.01742311108843653" errorMinus="0.01742311108843653"/>
+ <measurement value="0.00293" errorPlus="5.801104032854436E-4" errorMinus="5.801104032854436E-4"/>
</dataPoint>
<dataPoint>
<measurement value="17.3" errorPlus="2.6999999999999993" errorMinus="2.3000000000000007"/>
- <measurement value="0.0557" errorPlus="0.012701554235604397" errorMinus="0.012701554235604397"/>
+ <measurement value="0.00186" errorPlus="4.2414525813687934E-4" errorMinus="4.2414525813687934E-4"/>
</dataPoint>
<dataPoint>
<measurement value="24.1" errorPlus="5.899999999999999" errorMinus="4.100000000000001"/>
- <measurement value="0.0247" errorPlus="0.005632466599989741" errorMinus="0.005632466599989741"/>
+ <measurement value="8.22E-4" errorPlus="1.874448398862983E-4" errorMinus="1.874448398862983E-4"/>
</dataPoint>
<dataPoint>
<measurement value="39.8" errorPlus="40.2" errorMinus="9.799999999999997"/>
- <measurement value="0.00403" errorPlus="8.87515436485473E-4" errorMinus="8.87515436485473E-4"/>
+ <measurement value="1.34E-4" errorPlus="2.9510438831030623E-5" errorMinus="2.9510438831030623E-5"/>
</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"/>
+ <measurement value="0.00619" errorPlus="0.0017736267391985272" errorMinus="0.0017736267391985272"/>
</dataPoint>
<dataPoint>
<measurement value="2.38" errorPlus="0.1299999999999999" errorMinus="0.17999999999999972"/>
- <measurement value="0.582" errorPlus="0.14208430736714028" errorMinus="0.14208430736714028"/>
+ <measurement value="0.0194" errorPlus="0.004736143578904677" errorMinus="0.004736143578904677"/>
</dataPoint>
<dataPoint>
<measurement value="2.6" errorPlus="0.06999999999999984" errorMinus="0.0900000000000003"/>
- <measurement value="1.35" errorPlus="0.31861144047256057" errorMinus="0.31861144047256057"/>
+ <measurement value="0.0449" errorPlus="0.010596780501642939" errorMinus="0.010596780501642939"/>
</dataPoint>
<dataPoint>
<measurement value="2.76" errorPlus="0.07000000000000028" errorMinus="0.08999999999999986"/>
- <measurement value="1.99" errorPlus="0.4230798033468391" errorMinus="0.4230798033468391"/>
+ <measurement value="0.0664" errorPlus="0.014116833639311615" errorMinus="0.014116833639311615"/>
</dataPoint>
<dataPoint>
<measurement value="2.92" errorPlus="0.06000000000000005" errorMinus="0.08999999999999986"/>
- <measurement value="3.55" errorPlus="0.7028641404994281" errorMinus="0.7028641404994281"/>
+ <measurement value="0.118" errorPlus="0.023362808050403524" errorMinus="0.023362808050403524"/>
</dataPoint>
<dataPoint>
<measurement value="3.07" errorPlus="0.07000000000000028" errorMinus="0.08999999999999986"/>
- <measurement value="6.89" errorPlus="1.185400421798474" errorMinus="1.185400421798474"/>
+ <measurement value="0.23" errorPlus="0.03957069622839608" errorMinus="0.03957069622839608"/>
</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"/>
+ <measurement value="0.177" errorPlus="0.03381580251894076" errorMinus="0.03381580251894076"/>
</dataPoint>
<dataPoint>
<measurement value="0.15" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/>
- <measurement value="4.49" errorPlus="0.8889746453077276" errorMinus="0.8889746453077276"/>
+ <measurement value="0.15" errorPlus="0.029698484809835" errorMinus="0.029698484809835"/>
</dataPoint>
<dataPoint>
<measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
- <measurement value="4.58" errorPlus="0.8750077713940602" errorMinus="0.8750077713940602"/>
+ <measurement value="0.153" errorPlus="0.029230608957050482" errorMinus="0.029230608957050482"/>
</dataPoint>
<dataPoint>
<measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/>
- <measurement value="3.44" errorPlus="0.7313540319161439" errorMinus="0.7313540319161439"/>
+ <measurement value="0.115" errorPlus="0.024449335369289695" errorMinus="0.024449335369289695"/>
</dataPoint>
<dataPoint>
<measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
- <measurement value="3.17" errorPlus="0.6981200827937841" errorMinus="0.6981200827937841"/>
+ <measurement value="0.106" errorPlus="0.023344078478277953" errorMinus="0.023344078478277953"/>
</dataPoint>
<dataPoint>
<measurement value="0.58" errorPlus="0.12" errorMinus="0.07999999999999996"/>
- <measurement value="1.52" errorPlus="0.37107929071830453" errorMinus="0.37107929071830453"/>
+ <measurement value="0.0508" errorPlus="0.01240186050558544" errorMinus="0.01240186050558544"/>
</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"/>
+ <measurement value="1.64E-4" errorPlus="4.8233847037117E-5" errorMinus="4.8233847037117E-5"/>
</dataPoint>
<dataPoint>
<measurement value="9.3" errorPlus="5.699999999999999" errorMinus="4.300000000000001"/>
- <measurement value="0.0275" errorPlus="0.0056425282453878775" errorMinus="0.0056425282453878775"/>
+ <measurement value="1.02E-4" errorPlus="2.0928650219256852E-5" errorMinus="2.0928650219256852E-5"/>
</dataPoint>
<dataPoint>
<measurement value="24.3" errorPlus="15.7" errorMinus="9.3"/>
- <measurement value="0.012" errorPlus="0.0024738633753705963" errorMinus="0.0024738633753705963"/>
+ <measurement value="4.46E-5" errorPlus="9.194525545127382E-6" errorMinus="9.194525545127382E-6"/>
</dataPoint>
<dataPoint>
<measurement value="58.1" errorPlus="41.9" errorMinus="18.1"/>
- <measurement value="0.0018" errorPlus="4.5429946070846267E-4" errorMinus="4.5429946070846267E-4"/>
+ <measurement value="6.67E-6" errorPlus="1.6834318905141366E-6" errorMinus="1.6834318905141366E-6"/>
</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"/>
+ <measurement value="3.63E-4" errorPlus="9.874667589341931E-5" errorMinus="9.874667589341931E-5"/>
</dataPoint>
<dataPoint>
<measurement value="2.87" errorPlus="0.10999999999999988" errorMinus="0.20000000000000018"/>
- <measurement value="0.928" errorPlus="0.1837346260235125" errorMinus="0.1837346260235125"/>
+ <measurement value="0.00344" errorPlus="6.810852516388825E-4" errorMinus="6.810852516388825E-4"/>
</dataPoint>
<dataPoint>
<measurement value="3.08" errorPlus="0.06000000000000005" errorMinus="0.10000000000000009"/>
- <measurement value="3.21" errorPlus="0.5715242514539518" errorMinus="0.5715242514539518"/>
+ <measurement value="0.0119" errorPlus="0.002118734763957018" errorMinus="0.002118734763957018"/>
</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"/>
+ <measurement value="0.00758" errorPlus="0.0013976829397256017" errorMinus="0.0013976829397256017"/>
</dataPoint>
<dataPoint>
<measurement value="0.3" errorPlus="0.10000000000000003" errorMinus="0.09999999999999998"/>
- <measurement value="1.38" errorPlus="0.26364862980869064" errorMinus="0.26364862980869064"/>
+ <measurement value="0.00511" errorPlus="9.762641292191371E-4" errorMinus="9.762641292191371E-4"/>
</dataPoint>
<dataPoint>
<measurement value="0.53" errorPlus="0.16999999999999993" errorMinus="0.13"/>
- <measurement value="0.762" errorPlus="0.17983845751117863" errorMinus="0.17983845751117863"/>
+ <measurement value="0.00282" errorPlus="6.655438978760154E-4" errorMinus="6.655438978760154E-4"/>
</dataPoint>
</dataPointSet>
</aida>
Modified: trunk/src/Analyses/D0_2010_S8570965.cc
==============================================================================
--- trunk/src/Analyses/D0_2010_S8570965.cc Tue Apr 13 21:04:54 2010 (r2401)
+++ trunk/src/Analyses/D0_2010_S8570965.cc Mon Apr 19 12:28:06 2010 (r2402)
@@ -4,7 +4,6 @@
#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 {
@@ -63,17 +62,16 @@
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 (p.genParticle().barcode()!=photon.genParticle().barcode() &&
+ deltaR(eta_P, phi_P, p.momentum().eta(), p.momentum().phi()) < 0.4) {
+ Etsum += p.momentum().Et();
}
}
- if (mom_in_cone.E()/photon.momentum().E() < 1.1 && Etsum<1.5*GeV) {
+ if (Etsum<2.5*GeV) {
isolated_photons.push_back(photon);
}
}
@@ -91,12 +89,20 @@
FourMomentum yy=y1+y2;
double Myy = yy.mass()/GeV;
+ if (Myy<30.0 || Myy>350.0) {
+ vetoEvent;
+ }
+
double pTyy = yy.pT()/GeV;
if (Myy<pTyy) {
vetoEvent;
}
double dPhiyy = mapAngle0ToPi(y1.phi()-y2.phi());
+ if (dPhiyy<0.5*M_PI) {
+ vetoEvent;
+ }
+
double costhetayy = fabs(tanh((y1.eta()-y2.eta())/2.0));
_h_M->fill(Myy, weight);
@@ -116,13 +122,10 @@
scale(_h_pT, crossSection()/sumOfWeights());
scale(_h_dPhi, crossSection()/sumOfWeights());
scale(_h_costheta, crossSection()/sumOfWeights());
- // despite what it says in the paper, the mass bin width doesn't seem
- // to have been taken into account for these:
- for (size_t i=0; i<3; ++i) {
- scale(_h_pT_M.getHistograms()[i], crossSection()/sumOfWeights());
- scale(_h_dPhi_M.getHistograms()[i], crossSection()/sumOfWeights());
- scale(_h_costheta_M.getHistograms()[i], crossSection()/sumOfWeights());
- }
+
+ _h_pT_M.scale(crossSection()/sumOfWeights(), this);
+ _h_dPhi_M.scale(crossSection()/sumOfWeights(), this);
+ _h_costheta_M.scale(crossSection()/sumOfWeights(), this);
}
More information about the Rivet-svn
mailing list
|