[Rivet-svn] r2914 - in trunk: data/anainfo data/plotinfo data/refdata src/Analyses

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Thu Feb 3 13:50:28 GMT 2011


Author: fsiegert
Date: Thu Feb  3 13:50:27 2011
New Revision: 2914

Log:
Add D0_2010_S8821313. Sherpa events look reasonable (within 20% of the data), so marking as VALIDATED.

Added:
   trunk/data/anainfo/D0_2010_S8821313.info
   trunk/data/plotinfo/D0_2010_S8821313.plot
   trunk/data/refdata/D0_2010_S8821313.aida
   trunk/src/Analyses/D0_2010_S8821313.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_S8821313.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/D0_2010_S8821313.info	Thu Feb  3 13:50:27 2011	(r2914)
@@ -0,0 +1,36 @@
+Name: D0_2010_S8821313
+Year: 2010
+Summary: Precise study of Z pT using novel technique
+Experiment: D0
+Collider: Tevatron Run 2
+SpiresID: 8821313
+Status: VALIDATED
+Authors:
+ - Frank Siegert <frank.siegert at cern.ch>
+References:
+ - arXiv:1010.0262
+RunInfo:
+  Inclusive $Z/\gamma^*$ production in both electron and muon channels.
+  Cut on invariant lepton mass should be wider than $70<m_{\ell\ell}<110$ GeV.
+NumEvents: 1000000
+Beams: [p-, p+]
+Energies: [1960]
+Description:
+  Using 7.3 pb-1 the distribution of the variable $\phi^*$ is measured, which
+  probes the same physical effects as the $Z/\gamma^*$ boson transverse
+  momentum, but is less susceptible to the effects of experimental resolution
+  and efficiency.
+BibKey: Abazov:2010mk
+BibTeX: '@Article{Abazov:2010mk,
+     author    = "Abazov, Victor Mukhamedovich and others",
+ collaboration = "D0",
+     title     = "{Precise study of the Z/gamma* boson transverse momentum
+                  distribution in ppbar collisions using a novel technique}",
+     year      = "2010",
+     eprint    = "1010.0262",
+     archivePrefix = "arXiv",
+     primaryClass  =  "hep-ex",
+     SLACcitation  = "%%CITATION = 1010.0262;%%"
+}'
+ToDo:
+

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Wed Feb  2 21:34:24 2011	(r2913)
+++ trunk/data/anainfo/Makefile.am	Thu Feb  3 13:50:27 2011	(r2914)
@@ -71,6 +71,7 @@
   D0_2010_S8566488.info \
   D0_2010_S8570965.info \
   D0_2010_S8671338.info \
+  D0_2010_S8821313.info \
   E735_1998_S3905616.info \
   DELPHI_1995_S3137023.info \
   DELPHI_1996_S3430090.info \

Added: trunk/data/plotinfo/D0_2010_S8821313.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/D0_2010_S8821313.plot	Thu Feb  3 13:50:27 2011	(r2914)
@@ -0,0 +1,39 @@
+# BEGIN PLOT /D0_2010_S8821313/d01-x01-y01
+Title=Electron channel ($|y_Z|<1$)
+XLabel=$\phi^*_\eta$
+YLabel=$1/\sigma\;\mathrm{d}\sigma/\mathrm{d}\phi^*_\eta$
+LogX=1
+XMin=5e-3
+# END PLOT
+
+# BEGIN PLOT /D0_2010_S8821313/d01-x01-y02
+Title=Electron channel ($1<|y_Z|<2$)
+XLabel=$\phi^*_\eta$
+YLabel=$1/\sigma\;\mathrm{d}\sigma/\mathrm{d}\phi^*_\eta$
+LogX=1
+XMin=5e-3
+# END PLOT
+
+# BEGIN PLOT /D0_2010_S8821313/d01-x01-y03
+Title=Electron channel ($|y_Z|>2$)
+XLabel=$\phi^*_\eta$
+YLabel=$1/\sigma\;\mathrm{d}\sigma/\mathrm{d}\phi^*_\eta$
+LogX=1
+XMin=5e-3
+# END PLOT
+
+# BEGIN PLOT /D0_2010_S8821313/d02-x01-y01
+Title=Muon channel ($|y_Z|<1$)
+XLabel=$\phi^*_\eta$
+YLabel=$1/\sigma\;\mathrm{d}\sigma/\mathrm{d}\phi^*_\eta$
+LogX=1
+XMin=5e-3
+# END PLOT
+
+# BEGIN PLOT /D0_2010_S8821313/d02-x01-y02
+Title=Muon channel ($1<|y_Z|<2$)
+XLabel=$\phi^*_\eta$
+YLabel=$1/\sigma\;\mathrm{d}\sigma/\mathrm{d}\phi^*_\eta$
+LogX=1
+XMin=5e-3
+# END PLOT

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Wed Feb  2 21:34:24 2011	(r2913)
+++ trunk/data/plotinfo/Makefile.am	Thu Feb  3 13:50:27 2011	(r2914)
@@ -70,6 +70,7 @@
   D0_2010_S8566488.plot \
   D0_2010_S8570965.plot \
   D0_2010_S8671338.plot \
+  D0_2010_S8821313.plot \
   E735_1998_S3905616.plot \
   DELPHI_1995_S3137023.plot \
   DELPHI_1996_S3430090.plot \

Added: trunk/data/refdata/D0_2010_S8821313.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/refdata/D0_2010_S8821313.aida	Thu Feb  3 13:50:27 2011	(r2914)
@@ -0,0 +1,579 @@
+<?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_S8821313" title="(1/SIG)*D(SIG)/DPHI" >
+    <dataPoint>
+      <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/>
+      <measurement value="13.242" errorPlus="0.06800735254367722" errorMinus="0.06800735254367722"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/>
+      <measurement value="12.006" errorPlus="0.06315061361538778" errorMinus="0.06315061361538778"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/>
+      <measurement value="10.429" errorPlus="0.05922837157984339" errorMinus="0.05922837157984339"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/>
+      <measurement value="8.756" errorPlus="0.05360037313302959" errorMinus="0.05360037313302959"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/>
+      <measurement value="7.183" errorPlus="0.04850773134253961" errorMinus="0.04850773134253961"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/>
+      <measurement value="5.911" errorPlus="0.043289721643826724" errorMinus="0.043289721643826724"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0655" errorPlus="0.005499999999999991" errorMinus="0.005500000000000005"/>
+      <measurement value="4.762" errorPlus="0.039319206502675" errorMinus="0.039319206502675"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.076" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="4.07" errorPlus="0.03512833614050059" errorMinus="0.03512833614050059"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.087" errorPlus="0.006000000000000005" errorMinus="0.0059999999999999915"/>
+      <measurement value="3.387" errorPlus="0.031400636936215164" errorMinus="0.031400636936215164"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0995" errorPlus="0.006499999999999992" errorMinus="0.006500000000000006"/>
+      <measurement value="2.806" errorPlus="0.02630589287593181" errorMinus="0.02630589287593181"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.11349999999999999" errorPlus="0.007500000000000007" errorMinus="0.007499999999999993"/>
+      <measurement value="2.279" errorPlus="0.022203603311174516" errorMinus="0.022203603311174516"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.13" errorPlus="0.009000000000000008" errorMinus="0.009000000000000008"/>
+      <measurement value="1.83" errorPlus="0.018248287590894658" errorMinus="0.018248287590894658"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15050000000000002" errorPlus="0.011499999999999982" errorMinus="0.01150000000000001"/>
+      <measurement value="1.414" errorPlus="0.014317821063276354" errorMinus="0.014317821063276354"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.176" errorPlus="0.014000000000000012" errorMinus="0.013999999999999985"/>
+      <measurement value="1.084" errorPlus="0.011180339887498947" errorMinus="0.011180339887498947"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.20850000000000002" errorPlus="0.01849999999999999" errorMinus="0.018500000000000016"/>
+      <measurement value="0.75" errorPlus="0.00824621125123532" errorMinus="0.00824621125123532"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.251" errorPlus="0.02400000000000002" errorMinus="0.023999999999999994"/>
+      <measurement value="0.513" errorPlus="0.00608276253029822" errorMinus="0.00608276253029822"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.30600000000000005" errorPlus="0.030999999999999972" errorMinus="0.031000000000000028"/>
+      <measurement value="0.333" errorPlus="0.00412310562561766" errorMinus="0.00412310562561766"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3775" errorPlus="0.04049999999999998" errorMinus="0.04049999999999998"/>
+      <measurement value="0.197" errorPlus="0.0030" errorMinus="0.0030"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.47050000000000003" errorPlus="0.05249999999999999" errorMinus="0.05250000000000005"/>
+      <measurement value="0.115" errorPlus="0.0020" errorMinus="0.0020"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.06699999999999995" errorMinus="0.06700000000000006"/>
+      <measurement value="0.0617" errorPlus="0.0012041594578792293" errorMinus="0.0012041594578792293"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.742" errorPlus="0.08499999999999996" errorMinus="0.08499999999999996"/>
+      <measurement value="0.0321" errorPlus="8.062257748298549E-4" errorMinus="8.062257748298549E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9339999999999999" errorPlus="0.10699999999999998" errorMinus="0.10699999999999998"/>
+      <measurement value="0.0165" errorPlus="5.099019513592785E-4" errorMinus="5.099019513592785E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.1749999999999998" errorPlus="0.13400000000000012" errorMinus="0.1339999999999999"/>
+      <measurement value="0.00796" errorPlus="3.623534186398688E-4" errorMinus="3.623534186398688E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.4745" errorPlus="0.16549999999999998" errorMinus="0.16549999999999998"/>
+      <measurement value="0.00388" errorPlus="2.193171219946131E-4" errorMinus="2.193171219946131E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.8445" errorPlus="0.2044999999999999" errorMinus="0.20450000000000013"/>
+      <measurement value="0.00201" errorPlus="1.4317821063276353E-4" errorMinus="1.4317821063276353E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.298" errorPlus="0.2490000000000001" errorMinus="0.2490000000000001"/>
+      <measurement value="0.00107" errorPlus="9.486832980505139E-5" errorMinus="9.486832980505139E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.849" errorPlus="0.3019999999999996" errorMinus="0.30200000000000005"/>
+      <measurement value="7.02E-4" errorPlus="7.261542535852834E-5" errorMinus="7.261542535852834E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.5145" errorPlus="0.36350000000000016" errorMinus="0.36350000000000016"/>
+      <measurement value="3.89E-4" errorPlus="4.7434164902525696E-5" errorMinus="4.7434164902525696E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.3134999999999994" errorPlus="0.4355000000000002" errorMinus="0.43549999999999933"/>
+      <measurement value="2.84E-4" errorPlus="3.8275318418009274E-5" errorMinus="3.8275318418009274E-5"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x01-y02" dimension="2" path="/REF/D0_2010_S8821313" title="(1/SIG)*D(SIG)/DPHI" >
+    <dataPoint>
+      <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/>
+      <measurement value="14.235" errorPlus="0.11672617529928751" errorMinus="0.11672617529928751"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/>
+      <measurement value="12.782" errorPlus="0.10916501271011697" errorMinus="0.10916501271011697"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/>
+      <measurement value="11.035" errorPlus="0.09676776322722355" errorMinus="0.09676776322722355"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/>
+      <measurement value="9.023" errorPlus="0.08514693182963201" errorMinus="0.08514693182963201"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/>
+      <measurement value="7.268" errorPlus="0.07481310045707235" errorMinus="0.07481310045707235"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/>
+      <measurement value="5.911" errorPlus="0.06736467917239718" errorMinus="0.06736467917239718"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0655" errorPlus="0.005499999999999991" errorMinus="0.005500000000000005"/>
+      <measurement value="4.988" errorPlus="0.061294371682887816" errorMinus="0.061294371682887816"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.076" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="4.029" errorPlus="0.05346026561849464" errorMinus="0.05346026561849464"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.087" errorPlus="0.006000000000000005" errorMinus="0.0059999999999999915"/>
+      <measurement value="3.346" errorPlus="0.04738143096192853" errorMinus="0.04738143096192853"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0995" errorPlus="0.006499999999999992" errorMinus="0.006500000000000006"/>
+      <measurement value="2.696" errorPlus="0.04079215610874228" errorMinus="0.04079215610874228"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.11349999999999999" errorPlus="0.007500000000000007" errorMinus="0.007499999999999993"/>
+      <measurement value="2.267" errorPlus="0.0354400902933387" errorMinus="0.0354400902933387"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.13" errorPlus="0.009000000000000008" errorMinus="0.009000000000000008"/>
+      <measurement value="1.78" errorPlus="0.028792360097775936" errorMinus="0.028792360097775936"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15050000000000002" errorPlus="0.011499999999999982" errorMinus="0.01150000000000001"/>
+      <measurement value="1.308" errorPlus="0.02213594362117866" errorMinus="0.02213594362117866"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.176" errorPlus="0.014000000000000012" errorMinus="0.013999999999999985"/>
+      <measurement value="1.015" errorPlus="0.01676305461424021" errorMinus="0.01676305461424021"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.20850000000000002" errorPlus="0.01849999999999999" errorMinus="0.018500000000000016"/>
+      <measurement value="0.683" errorPlus="0.013000000000000001" errorMinus="0.013000000000000001"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.251" errorPlus="0.02400000000000002" errorMinus="0.023999999999999994"/>
+      <measurement value="0.444" errorPlus="0.008544003745317531" errorMinus="0.008544003745317531"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.30600000000000005" errorPlus="0.030999999999999972" errorMinus="0.031000000000000028"/>
+      <measurement value="0.282" errorPlus="0.00670820393249937" errorMinus="0.00670820393249937"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3775" errorPlus="0.04049999999999998" errorMinus="0.04049999999999998"/>
+      <measurement value="0.157" errorPlus="0.004472135954999579" errorMinus="0.004472135954999579"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.47050000000000003" errorPlus="0.05249999999999999" errorMinus="0.05250000000000005"/>
+      <measurement value="0.089" errorPlus="0.0031622776601683794" errorMinus="0.0031622776601683794"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.06699999999999995" errorMinus="0.06700000000000006"/>
+      <measurement value="0.0453" errorPlus="0.0016492422502470644" errorMinus="0.0016492422502470644"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.742" errorPlus="0.08499999999999996" errorMinus="0.08499999999999996"/>
+      <measurement value="0.0229" errorPlus="0.001004987562112089" errorMinus="0.001004987562112089"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9339999999999999" errorPlus="0.10699999999999998" errorMinus="0.10699999999999998"/>
+      <measurement value="0.0109" errorPlus="6.08276253029822E-4" errorMinus="6.08276253029822E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.1749999999999998" errorPlus="0.13400000000000012" errorMinus="0.1339999999999999"/>
+      <measurement value="0.00491" errorPlus="3.956008088970497E-4" errorMinus="3.956008088970497E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.4745" errorPlus="0.16549999999999998" errorMinus="0.16549999999999998"/>
+      <measurement value="0.00285" errorPlus="2.6925824035672517E-4" errorMinus="2.6925824035672517E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.8445" errorPlus="0.2044999999999999" errorMinus="0.20450000000000013"/>
+      <measurement value="0.00133" errorPlus="1.6492422502470644E-4" errorMinus="1.6492422502470644E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.298" errorPlus="0.2490000000000001" errorMinus="0.2490000000000001"/>
+      <measurement value="9.21E-4" errorPlus="1.2781236246936368E-4" errorMinus="1.2781236246936368E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.849" errorPlus="0.3019999999999996" errorMinus="0.30200000000000005"/>
+      <measurement value="3.63E-4" errorPlus="7.156814934033155E-5" errorMinus="7.156814934033155E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.5145" errorPlus="0.36350000000000016" errorMinus="0.36350000000000016"/>
+      <measurement value="2.26E-4" errorPlus="5.1971145840745135E-5" errorMinus="5.1971145840745135E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.3134999999999994" errorPlus="0.4355000000000002" errorMinus="0.43549999999999933"/>
+      <measurement value="2.14E-4" errorPlus="4.972926703662543E-5" errorMinus="4.972926703662543E-5"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x01-y03" dimension="2" path="/REF/D0_2010_S8821313" title="(1/SIG)*D(SIG)/DPHI" >
+    <dataPoint>
+      <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/>
+      <measurement value="15.625" errorPlus="0.3623285801589491" errorMinus="0.3623285801589491"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/>
+      <measurement value="14.288" errorPlus="0.3447680379617577" errorMinus="0.3447680379617577"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/>
+      <measurement value="12.13" errorPlus="0.32040755296965145" errorMinus="0.32040755296965145"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/>
+      <measurement value="9.514" errorPlus="0.28229417280560365" errorMinus="0.28229417280560365"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/>
+      <measurement value="7.572" errorPlus="0.2501619475459847" errorMinus="0.2501619475459847"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/>
+      <measurement value="6.311" errorPlus="0.22656566377101364" errorMinus="0.22656566377101364"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0655" errorPlus="0.005499999999999991" errorMinus="0.005500000000000005"/>
+      <measurement value="5.052" errorPlus="0.20241788458533008" errorMinus="0.20241788458533008"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.076" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="3.991" errorPlus="0.17518276170902203" errorMinus="0.17518276170902203"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.087" errorPlus="0.006000000000000005" errorMinus="0.0059999999999999915"/>
+      <measurement value="3.206" errorPlus="0.15221038072352358" errorMinus="0.15221038072352358"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0995" errorPlus="0.006499999999999992" errorMinus="0.006500000000000006"/>
+      <measurement value="2.533" errorPlus="0.12614277624977183" errorMinus="0.12614277624977183"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.11349999999999999" errorPlus="0.007500000000000007" errorMinus="0.007499999999999993"/>
+      <measurement value="1.796" errorPlus="0.09918165152889924" errorMinus="0.09918165152889924"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.13" errorPlus="0.009000000000000008" errorMinus="0.009000000000000008"/>
+      <measurement value="1.658" errorPlus="0.08720665112249179" errorMinus="0.08720665112249179"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15050000000000002" errorPlus="0.011499999999999982" errorMinus="0.01150000000000001"/>
+      <measurement value="1.223" errorPlus="0.06702984409947557" errorMinus="0.06702984409947557"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.176" errorPlus="0.014000000000000012" errorMinus="0.013999999999999985"/>
+      <measurement value="0.767" errorPlus="0.047265209192385894" errorMinus="0.047265209192385894"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.20850000000000002" errorPlus="0.01849999999999999" errorMinus="0.018500000000000016"/>
+      <measurement value="0.605" errorPlus="0.03712142238654117" errorMinus="0.03712142238654117"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.251" errorPlus="0.02400000000000002" errorMinus="0.023999999999999994"/>
+      <measurement value="0.378" errorPlus="0.025179356624028346" errorMinus="0.025179356624028346"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.30600000000000005" errorPlus="0.030999999999999972" errorMinus="0.031000000000000028"/>
+      <measurement value="0.195" errorPlus="0.0161245154965971" errorMinus="0.0161245154965971"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3775" errorPlus="0.04049999999999998" errorMinus="0.04049999999999998"/>
+      <measurement value="0.114" errorPlus="0.01104536101718726" errorMinus="0.01104536101718726"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.47050000000000003" errorPlus="0.05249999999999999" errorMinus="0.05250000000000005"/>
+      <measurement value="0.05" errorPlus="0.00608276253029822" errorMinus="0.00608276253029822"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.06699999999999995" errorMinus="0.06700000000000006"/>
+      <measurement value="0.0245" errorPlus="0.003847076812334269" errorMinus="0.003847076812334269"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.742" errorPlus="0.08499999999999996" errorMinus="0.08499999999999996"/>
+      <measurement value="0.00414" errorPlus="0.0013852075656738234" errorMinus="0.0013852075656738234"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9339999999999999" errorPlus="0.10699999999999998" errorMinus="0.10699999999999998"/>
+      <measurement value="0.00291" errorPlus="0.001110180165558726" errorMinus="0.001110180165558726"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.1749999999999998" errorPlus="0.13400000000000012" errorMinus="0.1339999999999999"/>
+      <measurement value="5.03E-4" errorPlus="3.568080716575789E-4" errorMinus="3.568080716575789E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.4745" errorPlus="0.16549999999999998" errorMinus="0.16549999999999998"/>
+      <measurement value="1.54E-4" errorPlus="1.5463505424062165E-4" errorMinus="1.5463505424062165E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.8445" errorPlus="0.2044999999999999" errorMinus="0.20450000000000013"/>
+      <measurement value="1.01E-4" errorPlus="1.0140019723846694E-4" errorMinus="1.0140019723846694E-4"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/D0_2010_S8821313" title="(1/SIG)*D(SIG)/DPHI" >
+    <dataPoint>
+      <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/>
+      <measurement value="12.992" errorPlus="0.06397655820689326" errorMinus="0.06397655820689326"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/>
+      <measurement value="11.958" errorPlus="0.059236812878479546" errorMinus="0.059236812878479546"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/>
+      <measurement value="10.263" errorPlus="0.05316013544000805" errorMinus="0.05316013544000805"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/>
+      <measurement value="8.62" errorPlus="0.04876474136094643" errorMinus="0.04876474136094643"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/>
+      <measurement value="7.115" errorPlus="0.04414748010928823" errorMinus="0.04414748010928823"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/>
+      <measurement value="5.863" errorPlus="0.0396232255123179" errorMinus="0.0396232255123179"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0655" errorPlus="0.005499999999999991" errorMinus="0.005500000000000005"/>
+      <measurement value="4.868" errorPlus="0.03551056180912941" errorMinus="0.03551056180912941"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.076" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="4.084" errorPlus="0.03157530680769389" errorMinus="0.03157530680769389"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.087" errorPlus="0.006000000000000005" errorMinus="0.0059999999999999915"/>
+      <measurement value="3.4" errorPlus="0.02745906043549196" errorMinus="0.02745906043549196"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0995" errorPlus="0.006499999999999992" errorMinus="0.006500000000000006"/>
+      <measurement value="2.831" errorPlus="0.02433105012119288" errorMinus="0.02433105012119288"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.11349999999999999" errorPlus="0.007500000000000007" errorMinus="0.007499999999999993"/>
+      <measurement value="2.32" errorPlus="0.020223748416156685" errorMinus="0.020223748416156685"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.13" errorPlus="0.009000000000000008" errorMinus="0.009000000000000008"/>
+      <measurement value="1.85" errorPlus="0.016278820596099707" errorMinus="0.016278820596099707"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15050000000000002" errorPlus="0.011499999999999982" errorMinus="0.01150000000000001"/>
+      <measurement value="1.439" errorPlus="0.013341664064126334" errorMinus="0.013341664064126334"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.176" errorPlus="0.014000000000000012" errorMinus="0.013999999999999985"/>
+      <measurement value="1.061" errorPlus="0.01019803902718557" errorMinus="0.01019803902718557"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.20850000000000002" errorPlus="0.01849999999999999" errorMinus="0.018500000000000016"/>
+      <measurement value="0.779" errorPlus="0.007280109889280519" errorMinus="0.007280109889280519"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.251" errorPlus="0.02400000000000002" errorMinus="0.023999999999999994"/>
+      <measurement value="0.526" errorPlus="0.005385164807134504" errorMinus="0.005385164807134504"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.30600000000000005" errorPlus="0.030999999999999972" errorMinus="0.031000000000000028"/>
+      <measurement value="0.331" errorPlus="0.00412310562561766" errorMinus="0.00412310562561766"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3775" errorPlus="0.04049999999999998" errorMinus="0.04049999999999998"/>
+      <measurement value="0.207" errorPlus="0.0031622776601683794" errorMinus="0.0031622776601683794"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.47050000000000003" errorPlus="0.05249999999999999" errorMinus="0.05250000000000005"/>
+      <measurement value="0.117" errorPlus="0.0022360679774997894" errorMinus="0.0022360679774997894"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.06699999999999995" errorMinus="0.06700000000000006"/>
+      <measurement value="0.066" errorPlus="0.0013" errorMinus="0.0013"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.742" errorPlus="0.08499999999999996" errorMinus="0.08499999999999996"/>
+      <measurement value="0.0342" errorPlus="8.544003745317531E-4" errorMinus="8.544003745317531E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9339999999999999" errorPlus="0.10699999999999998" errorMinus="0.10699999999999998"/>
+      <measurement value="0.0171" errorPlus="5.8309518948453E-4" errorMinus="5.8309518948453E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.1749999999999998" errorPlus="0.13400000000000012" errorMinus="0.1339999999999999"/>
+      <measurement value="0.00875" errorPlus="4.0607881008493905E-4" errorMinus="4.0607881008493905E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.4745" errorPlus="0.16549999999999998" errorMinus="0.16549999999999998"/>
+      <measurement value="0.00477" errorPlus="2.4698178070456936E-4" errorMinus="2.4698178070456936E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.8445" errorPlus="0.2044999999999999" errorMinus="0.20450000000000013"/>
+      <measurement value="0.00249" errorPlus="1.5652475842498526E-4" errorMinus="1.5652475842498526E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.298" errorPlus="0.2490000000000001" errorMinus="0.2490000000000001"/>
+      <measurement value="0.00139" errorPlus="1.0770329614269009E-4" errorMinus="1.0770329614269009E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.849" errorPlus="0.3019999999999996" errorMinus="0.30200000000000005"/>
+      <measurement value="8.02E-4" errorPlus="7.057619995437555E-5" errorMinus="7.057619995437555E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.5145" errorPlus="0.36350000000000016" errorMinus="0.36350000000000016"/>
+      <measurement value="5.36E-4" errorPlus="5.679788728465171E-5" errorMinus="5.679788728465171E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.3134999999999994" errorPlus="0.4355000000000002" errorMinus="0.43549999999999933"/>
+      <measurement value="3.07E-4" errorPlus="3.6055512754639895E-5" errorMinus="3.6055512754639895E-5"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y02" dimension="2" path="/REF/D0_2010_S8821313" title="(1/SIG)*D(SIG)/DPHI" >
+    <dataPoint>
+      <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/>
+      <measurement value="13.404" errorPlus="0.1167304587500623" errorMinus="0.1167304587500623"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/>
+      <measurement value="12.008" errorPlus="0.10332473082471592" errorMinus="0.10332473082471592"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/>
+      <measurement value="10.647" errorPlus="0.09535197952848173" errorMinus="0.09535197952848173"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/>
+      <measurement value="8.755" errorPlus="0.08697700845625814" errorMinus="0.08697700845625814"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/>
+      <measurement value="7.23" errorPlus="0.07864477096412704" errorMinus="0.07864477096412704"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/>
+      <measurement value="5.804" errorPlus="0.07156814934033157" errorMinus="0.07156814934033157"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0655" errorPlus="0.005499999999999991" errorMinus="0.005500000000000005"/>
+      <measurement value="4.972" errorPlus="0.06395310782127793" errorMinus="0.06395310782127793"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.076" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="4.045" errorPlus="0.0565685424949238" errorMinus="0.0565685424949238"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.087" errorPlus="0.006000000000000005" errorMinus="0.0059999999999999915"/>
+      <measurement value="3.441" errorPlus="0.05063595560468866" errorMinus="0.05063595560468866"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0995" errorPlus="0.006499999999999992" errorMinus="0.006500000000000006"/>
+      <measurement value="2.82" errorPlus="0.04393176527297759" errorMinus="0.04393176527297759"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.11349999999999999" errorPlus="0.007500000000000007" errorMinus="0.007499999999999993"/>
+      <measurement value="2.33" errorPlus="0.036496575181789316" errorMinus="0.036496575181789316"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.13" errorPlus="0.009000000000000008" errorMinus="0.009000000000000008"/>
+      <measurement value="1.824" errorPlus="0.029274562336608894" errorMinus="0.029274562336608894"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15050000000000002" errorPlus="0.011499999999999982" errorMinus="0.01150000000000001"/>
+      <measurement value="1.414" errorPlus="0.023194827009486402" errorMinus="0.023194827009486402"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.176" errorPlus="0.014000000000000012" errorMinus="0.013999999999999985"/>
+      <measurement value="1.066" errorPlus="0.01811077027627483" errorMinus="0.01811077027627483"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.20850000000000002" errorPlus="0.01849999999999999" errorMinus="0.018500000000000016"/>
+      <measurement value="0.756" errorPlus="0.013038404810405297" errorMinus="0.013038404810405297"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.251" errorPlus="0.02400000000000002" errorMinus="0.023999999999999994"/>
+      <measurement value="0.514" errorPlus="0.009219544457292887" errorMinus="0.009219544457292887"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.30600000000000005" errorPlus="0.030999999999999972" errorMinus="0.031000000000000028"/>
+      <measurement value="0.326" errorPlus="0.007280109889280519" errorMinus="0.007280109889280519"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3775" errorPlus="0.04049999999999998" errorMinus="0.04049999999999998"/>
+      <measurement value="0.2" errorPlus="0.005099019513592785" errorMinus="0.005099019513592785"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.47050000000000003" errorPlus="0.05249999999999999" errorMinus="0.05250000000000005"/>
+      <measurement value="0.107" errorPlus="0.0030" errorMinus="0.0030"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.06699999999999995" errorMinus="0.06700000000000006"/>
+      <measurement value="0.0543" errorPlus="0.0019924858845171275" errorMinus="0.0019924858845171275"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.742" errorPlus="0.08499999999999996" errorMinus="0.08499999999999996"/>
+      <measurement value="0.0259" errorPlus="0.0012369316876852982" errorMinus="0.0012369316876852982"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9339999999999999" errorPlus="0.10699999999999998" errorMinus="0.10699999999999998"/>
+      <measurement value="0.0123" errorPlus="7.615773105863909E-4" errorMinus="7.615773105863909E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.1749999999999998" errorPlus="0.13400000000000012" errorMinus="0.1339999999999999"/>
+      <measurement value="0.0052" errorPlus="4.716990566028302E-4" errorMinus="4.716990566028302E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.4745" errorPlus="0.16549999999999998" errorMinus="0.16549999999999998"/>
+      <measurement value="0.00254" errorPlus="2.941088233970548E-4" errorMinus="2.941088233970548E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.8445" errorPlus="0.2044999999999999" errorMinus="0.20450000000000013"/>
+      <measurement value="0.00126" errorPlus="2.2847319317591724E-4" errorMinus="2.2847319317591724E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.298" errorPlus="0.2490000000000001" errorMinus="0.2490000000000001"/>
+      <measurement value="5.03E-4" errorPlus="1.1921828718783037E-4" errorMinus="1.1921828718783037E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.849" errorPlus="0.3019999999999996" errorMinus="0.30200000000000005"/>
+      <measurement value="3.5E-4" errorPlus="1.0630145812734648E-4" errorMinus="1.0630145812734648E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.5145" errorPlus="0.36350000000000016" errorMinus="0.36350000000000016"/>
+      <measurement value="2.09E-4" errorPlus="6.453681120105021E-5" errorMinus="6.453681120105021E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.3134999999999994" errorPlus="0.4355000000000002" errorMinus="0.43549999999999933"/>
+      <measurement value="6.9E-5" errorPlus="4.368065933568311E-5" errorMinus="4.368065933568311E-5"/>
+    </dataPoint>
+  </dataPointSet>
+</aida>

Modified: trunk/data/refdata/Makefile.am
==============================================================================
--- trunk/data/refdata/Makefile.am	Wed Feb  2 21:34:24 2011	(r2913)
+++ trunk/data/refdata/Makefile.am	Thu Feb  3 13:50:27 2011	(r2914)
@@ -105,6 +105,7 @@
   D0_2010_S8566488.aida          \
   D0_2010_S8570965.aida          \
   D0_2010_S8671338.aida          \
+  D0_2010_S8821313.aida          \
   E735_1992_S2485869.aida        \
   E735_1993_S2896508.aida        \
   E735_1998_S3905616.aida        \

Added: trunk/src/Analyses/D0_2010_S8821313.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/D0_2010_S8821313.cc	Thu Feb  3 13:50:27 2011	(r2914)
@@ -0,0 +1,131 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Tools/BinnedHistogram.hh"
+#include "Rivet/Tools/ParticleIdUtils.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ZFinder.hh"
+
+
+namespace Rivet {
+
+
+  class D0_2010_S8821313 : public Analysis {
+  public:
+
+    /// @name Constructors etc.
+    //@{
+
+    /// Constructor
+    D0_2010_S8821313()
+      : Analysis("D0_2010_S8821313")
+    {
+      setNeedsCrossSection(false);
+    }
+
+    //@}
+
+
+  public:
+
+    /// @name Analysis methods
+    //@{
+
+    /// Book histograms and initialise projections before the run
+    void init() {
+
+      /// Initialise and register projections
+      vector<pair<double, double> > etaRanges_ee;
+      etaRanges_ee.push_back(make_pair(-3.0, -1.5));
+      etaRanges_ee.push_back(make_pair(-1.1, 1.1));
+      etaRanges_ee.push_back(make_pair(1.5, 3.0));
+      ZFinder zfinder_ee(etaRanges_ee, 20.0*GeV, ELECTRON, 70.0*GeV, 110.0*GeV, 0.2);
+      addProjection(zfinder_ee, "zfinder_ee");
+
+      ZFinder zfinder_mm(-2.0, 2.0, 15.0*GeV, MUON, 70.0*GeV, 110.0*GeV, 0.0);
+      addProjection(zfinder_mm, "zfinder_mm");
+
+      /// Book histograms here
+      _h_phistar_ee.addHistogram(0.0, 1.0, bookHistogram1D(1, 1, 1));
+      _h_phistar_ee.addHistogram(1.0, 2.0, bookHistogram1D(1, 1, 2));
+      _h_phistar_ee.addHistogram(2.0, 10.0, bookHistogram1D(1, 1, 3));
+
+      _h_phistar_mm.addHistogram(0.0, 1.0, bookHistogram1D(2, 1, 1));
+      _h_phistar_mm.addHistogram(1.0, 2.0, bookHistogram1D(2, 1, 2));
+    }
+
+
+    /// Perform the per-event analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+
+      const ZFinder& zfinder_ee = applyProjection<ZFinder>(event, "zfinder_ee");
+      if (zfinder_ee.particles().size()==1) {
+        ParticleVector ee=zfinder_ee.constituentsFinalState().particlesByPt();
+        FourMomentum eminus=PID::threeCharge(ee[0].pdgId())<0.0?ee[0].momentum():ee[1].momentum();
+        FourMomentum eplus=PID::threeCharge(ee[0].pdgId())<0.0?ee[1].momentum():ee[0].momentum();
+        double phi_acop=M_PI-mapAngle0ToPi(eminus.phi()-eplus.phi());
+        double costhetastar=tanh((eminus.eta()-eplus.eta())/2.0);
+        double sin2thetastar=1.0-sqr(costhetastar);
+        if (sin2thetastar<0.0) sin2thetastar=0.0;
+        double phistar=tan(phi_acop/2.0)*sqrt(sin2thetastar);
+
+        FourMomentum Zmom=zfinder_ee.particles()[0].momentum();
+        _h_phistar_ee.fill(Zmom.rapidity(), phistar, weight);
+      }
+
+      const ZFinder& zfinder_mm = applyProjection<ZFinder>(event, "zfinder_mm");
+      if (zfinder_mm.particles().size()==1) {
+        ParticleVector mm=zfinder_mm.constituentsFinalState().particlesByPt();
+        FourMomentum mminus=PID::threeCharge(mm[0].pdgId())<0.0?mm[0].momentum():mm[1].momentum();
+        FourMomentum mplus=PID::threeCharge(mm[0].pdgId())<0.0?mm[1].momentum():mm[0].momentum();
+        double phi_acop=M_PI-mapAngle0ToPi(mminus.phi()-mplus.phi());
+        double costhetastar=tanh((mminus.eta()-mplus.eta())/2.0);
+        double sin2thetastar=1.0-sqr(costhetastar);
+        if (sin2thetastar<0.0) sin2thetastar=0.0;
+        double phistar=tan(phi_acop/2.0)*sqrt(sin2thetastar);
+
+        FourMomentum Zmom=zfinder_mm.particles()[0].momentum();
+        _h_phistar_mm.fill(Zmom.rapidity(), phistar, weight);
+      }
+    }
+
+
+    /// Normalise histograms etc., after the run
+    void finalize() {
+      foreach (AIDA::IHistogram1D* hist, _h_phistar_ee.getHistograms()) {
+        normalize(hist, 1.0);
+      }
+      foreach (AIDA::IHistogram1D* hist, _h_phistar_mm.getHistograms()) {
+        normalize(hist, 1.0);
+      }
+    }
+
+    //@}
+
+
+  private:
+
+    // Data members like post-cuts event weight counters go here
+
+
+  private:
+
+    /// @name Histograms
+    //@{
+
+    BinnedHistogram<double> _h_phistar_ee;
+    BinnedHistogram<double> _h_phistar_mm;
+    //@}
+
+
+  };
+
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<D0_2010_S8821313> plugin_D0_2010_S8821313;
+
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Wed Feb  2 21:34:24 2011	(r2913)
+++ trunk/src/Analyses/Makefile.am	Thu Feb  3 13:50:27 2011	(r2914)
@@ -111,7 +111,8 @@
     D0_2009_S8349509.cc \
     D0_2010_S8566488.cc \
     D0_2010_S8570965.cc \
-    D0_2010_S8671338.cc
+    D0_2010_S8671338.cc \
+    D0_2010_S8821313.cc
 
 
 lib_LTLIBRARIES += RivetHERAAnalyses.la


More information about the Rivet-svn mailing list