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

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Oct 21 21:30:25 BST 2009


Author: fsiegert
Date: Wed Oct 21 21:30:24 2009
New Revision: 1949

Log:
Add CDF_1996_S3349578.
Have to test it with higher statistics, thus leaving UNVALIDATED for
now.

Added:
   trunk/data/anainfo/CDF_1996_S3349578.info
   trunk/data/plotinfo/CDF_1996_S3349578.plot
   trunk/data/refdata/CDF_1996_S3349578.aida
   trunk/src/Analyses/CDF_1996_S3349578.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/CDF_1996_S3349578.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/CDF_1996_S3349578.info	Wed Oct 21 21:30:24 2009	(r1949)
@@ -0,0 +1,20 @@
+Name: CDF_1996_S3349578
+Year: 1996
+Summary: Further properties of high-mass multijet events
+Experiment: CDF
+Collider: Tevatron Run 1
+SpiresID: 3349578
+Status: UNVALIDATED
+Authors:
+ - Frank Siegert <frank.siegert at durham.ac.uk>
+References:
+ - Phys.Rev.D54:4221-4233,1996
+ - doi:10.1103/PhysRevD.54.4221
+ - arXiv:hep-ex/9605004
+RunInfo:
+  Pure QCD events without underlying event.
+NumEvents: 10000000
+PtCuts: [20]
+Description:
+  Multijet distributions corresponding to (4N-4) variables that span the N-body
+  parameter space in inclusive N=3, 4, 5 jet events.

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Wed Oct 21 21:29:18 2009	(r1948)
+++ trunk/data/anainfo/Makefile.am	Wed Oct 21 21:30:24 2009	(r1949)
@@ -6,6 +6,7 @@
   CDF_1990_S2089246.info \
   CDF_1994_S2952106.info \
   CDF_1996_S3108457.info \
+  CDF_1996_S3349578.info \
   CDF_2000_S4155203.info \
   CDF_2001_S4751469.info \
   CDF_2002_S4796047.info \

Added: trunk/data/plotinfo/CDF_1996_S3349578.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/CDF_1996_S3349578.plot	Wed Oct 21 21:30:24 2009	(r1949)
@@ -0,0 +1,284 @@
+# BEGIN PLOT /CDF_1996_S3349578/d01-x01-y01
+Title=Multijet mass in inclusive 3-jet events
+XLabel=$m_{3J}$ [GeV]
+YLabel=$1/N \mathrm{d}N/\mathrm{d}m$ [1/GeV]
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d01-x01-y02
+Title=Multijet mass in inclusive 4-jet events
+XLabel=$m_{4J}$ [GeV]
+YLabel=$1/N \mathrm{d}N/\mathrm{d}m$ [1/GeV]
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d01-x01-y03
+Title=Multijet mass in inclusive 5-jet events
+XLabel=$m_{5J}$ [GeV]
+YLabel=$1/N \mathrm{d}N/\mathrm{d}m$ [1/GeV]
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d02-x01-y01
+Title=Dalitz distribution in inclusive 3-jet events
+XLabel=$X_3$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_3$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d03-x01-y01
+Title=Dalitz distribution in inclusive 3-jet events
+XLabel=$X_4$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_4$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d04-x01-y01
+Title=Dalitz distribution in inclusive 4-jet events
+XLabel=$X_3^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_3^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d05-x01-y01
+Title=Dalitz distribution in inclusive 4-jet events
+XLabel=$X_4^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_4^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d06-x01-y01
+Title=Dalitz distribution in inclusive 5-jet events
+XLabel=$X_3^{''}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_3^{''}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d07-x01-y01
+Title=Dalitz distribution in inclusive 5-jet events
+XLabel=$X_4^{''}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_4^{''}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d08-x01-y01
+Title=Leading jet angle in inclusive 3-jet events
+XLabel=$|\cos \theta_3|$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos\theta_3$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d09-x01-y01
+Title=Angular distribution in inclusive 3-jet events
+XLabel=$\Psi_3$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi_3$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d10-x01-y01
+Title=Leading jet angle in inclusive 4-jet events
+XLabel=$|\cos \theta_{3^'}|$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos\theta_{3^'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d11-x01-y01
+Title=Angular distribution in inclusive 4-jet events
+XLabel=$\Psi_{3^'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi_{3^'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d12-x01-y01
+Title=Leading jet angle in inclusive 5-jet events
+XLabel=$|\cos \theta_{3^{''}}|$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos\theta_{3^{''}}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d13-x01-y01
+Title=Angular distribution in inclusive 5-jet events
+XLabel=$\Psi_{3^{''}}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi_{3^{''}}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d14-x01-y01
+Title=Single-jet mass fraction in inclusive 3-jet events
+XLabel=$f_3$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_3$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d14-x01-y02
+Title=Single-jet mass fraction in inclusive 3-jet events
+XLabel=$f_4$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_4$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d14-x01-y03
+Title=Single-jet mass fraction in inclusive 3-jet events
+XLabel=$f_5$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_5$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d15-x01-y01
+Title=Single-jet mass fraction in inclusive 4-jet events
+XLabel=$f_3^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_3^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d15-x01-y02
+Title=Single-jet mass fraction in inclusive 4-jet events
+XLabel=$f_4^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_4^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d15-x01-y03
+Title=Single-jet mass fraction in inclusive 4-jet events
+XLabel=$f_5^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_5^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d16-x01-y01
+Title=Single-jet mass fraction in inclusive 5-jet events
+XLabel=$f_3^{''}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_3^{''}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d16-x01-y02
+Title=Single-jet mass fraction in inclusive 5-jet events
+XLabel=$f_4^{''}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_4^{''}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d16-x01-y03
+Title=Single-jet mass fraction in inclusive 5-jet events
+XLabel=$f_5^{''}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_5^{''}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d17-x01-y01
+Title=Two-body energy sharing in inclusive 4-jet events
+XLabel=$X_A$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_A$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d18-x01-y01
+Title=Two-body energy sharing in inclusive 5-jet events
+XLabel=$X_A^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_A^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d18-x01-y02
+Title=Two-body energy sharing in inclusive 5-jet events
+XLabel=$X_C$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}X_C$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d19-x01-y01
+Title=Two-body angular distribution in inclusive 4-jet events
+XLabel=$\Psi^{'}_{AB}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi^{'}_{AB}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d20-x01-y01
+Title=Two-body angular distribution in inclusive 5-jet events
+XLabel=$\Psi^{''}_{A^{'}B^{'}}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi^{''}_{A^{'}B^{'}}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d20-x01-y02
+Title=Two-body angular distribution in inclusive 5-jet events
+XLabel=$\Psi^{''}_{CD}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi^{''}_{CD}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d21-x01-y01
+Title=Single-body mass fraction in inclusive 4-jet events
+XLabel=$f_A$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_A$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d21-x01-y02
+Title=Single-body mass fraction in inclusive 4-jet events
+XLabel=$f_B$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_B$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d22-x01-y01
+Title=Single-body mass fraction in inclusive 5-jet events
+XLabel=$f_A^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_A^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d23-x01-y01
+Title=Single-body mass fraction in inclusive 5-jet events
+XLabel=$f_B^{'}$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_B^{'}$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d24-x01-y01
+Title=Single-body mass fraction in inclusive 5-jet events
+XLabel=$f_C$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_C$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT
+
+# BEGIN PLOT /CDF_1996_S3349578/d25-x01-y01
+Title=Single-body mass fraction in inclusive 5-jet events
+XLabel=$f_D$
+YLabel=$1/N \mathrm{d}N/\mathrm{d}f_D$
+LogY=0
+NormalizeToIntegral=1
+# END PLOT

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Wed Oct 21 21:29:18 2009	(r1948)
+++ trunk/data/plotinfo/Makefile.am	Wed Oct 21 21:30:24 2009	(r1949)
@@ -6,6 +6,7 @@
   CDF_1990_S2089246.plot \
   CDF_1994_S2952106.plot \
   CDF_1996_S3108457.plot \
+  CDF_1996_S3349578.plot \
   CDF_2000_S4155203.plot \
   CDF_2001_S4751469.plot \
   CDF_2002_S4796047.plot \

Added: trunk/data/refdata/CDF_1996_S3349578.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/refdata/CDF_1996_S3349578.aida	Wed Oct 21 21:30:24 2009	(r1949)
@@ -0,0 +1,1595 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd">
+<aida version="3.3">
+  <implementation version="1.0" package="HepData"/>
+  <dataPointSet name="d01-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DM IN PER 50 GEV" >
+    <dataPoint>
+      <measurement value="625.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.01" errorPlus="0.0010" errorMinus="0.0010"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="675.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0052" errorPlus="8.0E-4" errorMinus="8.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="725.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0025" errorPlus="3.0E-4" errorMinus="3.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="775.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0015" errorPlus="2.0E-4" errorMinus="2.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="825.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="8.4E-4" errorPlus="1.4E-4" errorMinus="1.4E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="875.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="2.8E-4" errorPlus="7.0E-5" errorMinus="7.0E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="925.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="1.8E-4" errorPlus="6.0E-5" errorMinus="6.0E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="975.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="6.0E-5" errorPlus="4.5E-5" errorMinus="4.5E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1025.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1075.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1125.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x01-y02" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DM IN PER 50 GEV" >
+    <dataPoint>
+      <measurement value="625.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="675.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0094" errorPlus="9.0E-4" errorMinus="9.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="725.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0051" errorPlus="5.0E-4" errorMinus="5.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="775.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0029" errorPlus="3.0E-4" errorMinus="3.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="825.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0014" errorPlus="2.0E-4" errorMinus="2.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="875.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="9.2E-4" errorPlus="1.1E-4" errorMinus="1.1E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="925.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="3.4E-4" errorPlus="8.0E-5" errorMinus="8.0E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="975.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="1.5E-4" errorPlus="6.0E-5" errorMinus="6.0E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1025.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1075.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1125.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x01-y03" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DM IN PER 50 GEV" >
+    <dataPoint>
+      <measurement value="625.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="675.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="725.0" errorPlus="25.0" errorMinus="25.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="775.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0067" errorPlus="6.0E-4" errorMinus="6.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="825.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0051" errorPlus="5.0E-4" errorMinus="5.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="875.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0031" errorPlus="4.0E-4" errorMinus="4.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="925.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0022" errorPlus="3.0E-4" errorMinus="3.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="975.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="0.0013" errorPlus="2.0E-4" errorMinus="2.0E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1025.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="6.8E-4" errorPlus="1.2E-4" errorMinus="1.2E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1075.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="5.3E-4" errorPlus="1.2E-4" errorMinus="1.2E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1125.0" errorPlus="25.0" errorMinus="25.0"/>
+      <measurement value="3.3E-4" errorPlus="1.0E-4" errorMinus="1.0E-4"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DX3" >
+    <dataPoint>
+      <measurement value="0.679" errorPlus="0.0129999999999999" errorMinus="0.013000000000000012"/>
+      <measurement value="0.5" errorPlus="0.18" errorMinus="0.18"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.705" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="1.77" errorPlus="0.27" errorMinus="0.27"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.731" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="2.64" errorPlus="0.27" errorMinus="0.27"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.757" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="3.68" errorPlus="0.39" errorMinus="0.39"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.783" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="4.36" errorPlus="0.36" errorMinus="0.36"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.8089999999999999" errorPlus="0.013000000000000012" errorMinus="0.0129999999999999"/>
+      <measurement value="5.68" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.835" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="6.05" errorPlus="0.55" errorMinus="0.55"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.861" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="7.05" errorPlus="0.55" errorMinus="0.55"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.887" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="8.05" errorPlus="0.55" errorMinus="0.55"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DX4" >
+    <dataPoint>
+      <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="0.84" errorPlus="0.14" errorMinus="0.14"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.16" errorPlus="0.25" errorMinus="0.25"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="5.34" errorPlus="0.33" errorMinus="0.33"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.75" errorPlus="0.27" errorMinus="0.27"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.5" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="2.56" errorPlus="0.25" errorMinus="0.25"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="0.66" errorPlus="0.12" errorMinus="0.12"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DX3" >
+    <dataPoint>
+      <measurement value="0.679" errorPlus="0.0129999999999999" errorMinus="0.013000000000000012"/>
+      <measurement value="0.88" errorPlus="0.18" errorMinus="0.18"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.705" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="1.88" errorPlus="0.24" errorMinus="0.24"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.731" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="3.0" errorPlus="0.32" errorMinus="0.32"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.757" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="4.03" errorPlus="0.36" errorMinus="0.36"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.783" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="4.27" errorPlus="0.38" errorMinus="0.38"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.8089999999999999" errorPlus="0.013000000000000012" errorMinus="0.0129999999999999"/>
+      <measurement value="5.7" errorPlus="0.41" errorMinus="0.41"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.835" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="5.52" errorPlus="0.42" errorMinus="0.42"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.861" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="7.42" errorPlus="0.53" errorMinus="0.53"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.887" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="7.15" errorPlus="0.45" errorMinus="0.45"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DX4" >
+    <dataPoint>
+      <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="1.17" errorPlus="0.11" errorMinus="0.11"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.82" errorPlus="0.31" errorMinus="0.31"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="5.74" errorPlus="0.35" errorMinus="0.35"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="4.45" errorPlus="0.25" errorMinus="0.25"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.87" errorPlus="0.21" errorMinus="0.21"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="1.71" errorPlus="0.16" errorMinus="0.16"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="0.34" errorPlus="0.15" errorMinus="0.15"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d06-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DX3" >
+    <dataPoint>
+      <measurement value="0.687" errorPlus="0.020999999999999908" errorMinus="0.02100000000000002"/>
+      <measurement value="0.32" errorPlus="0.11" errorMinus="0.11"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.7284999999999999" errorPlus="0.020500000000000074" errorMinus="0.020499999999999963"/>
+      <measurement value="1.38" errorPlus="0.19" errorMinus="0.19"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.77" errorPlus="0.02100000000000002" errorMinus="0.02100000000000002"/>
+      <measurement value="2.18" errorPlus="0.28" errorMinus="0.28"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.812" errorPlus="0.020999999999999908" errorMinus="0.02100000000000002"/>
+      <measurement value="2.77" errorPlus="0.24" errorMinus="0.24"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.854" errorPlus="0.02100000000000002" errorMinus="0.02100000000000002"/>
+      <measurement value="3.9" errorPlus="0.31" errorMinus="0.31"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.8955" errorPlus="0.020500000000000074" errorMinus="0.020499999999999963"/>
+      <measurement value="3.65" errorPlus="0.34" errorMinus="0.34"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.937" errorPlus="0.020999999999999908" errorMinus="0.02100000000000002"/>
+      <measurement value="3.19" errorPlus="0.31" errorMinus="0.31"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.979" errorPlus="0.02100000000000002" errorMinus="0.02100000000000002"/>
+      <measurement value="2.67" errorPlus="0.34" errorMinus="0.34"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d07-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DX4" >
+    <dataPoint>
+      <measurement value="0.531" errorPlus="0.031000000000000028" errorMinus="0.031000000000000028"/>
+      <measurement value="0.65" errorPlus="0.14" errorMinus="0.14"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5935" errorPlus="0.03149999999999997" errorMinus="0.03149999999999997"/>
+      <measurement value="2.62" errorPlus="0.28" errorMinus="0.28"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.6565" errorPlus="0.03149999999999997" errorMinus="0.03149999999999997"/>
+      <measurement value="4.98" errorPlus="0.39" errorMinus="0.39"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.719" errorPlus="0.031000000000000028" errorMinus="0.031000000000000028"/>
+      <measurement value="3.79" errorPlus="0.33" errorMinus="0.33"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.781" errorPlus="0.031000000000000028" errorMinus="0.031000000000000028"/>
+      <measurement value="2.39" errorPlus="0.24" errorMinus="0.24"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.8435" errorPlus="0.03149999999999997" errorMinus="0.03149999999999997"/>
+      <measurement value="1.44" errorPlus="0.14" errorMinus="0.14"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9065" errorPlus="0.03149999999999997" errorMinus="0.03149999999999997"/>
+      <measurement value="0.25" errorPlus="0.19" errorMinus="0.19"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d08-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCOS(THETA(P=3))" >
+    <dataPoint>
+      <measurement value="-0.55" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/>
+      <measurement value="1.39" errorPlus="0.12" errorMinus="0.12"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="0.99" errorPlus="0.11" errorMinus="0.11"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.35" errorPlus="0.04999999999999999" errorMinus="0.050000000000000044"/>
+      <measurement value="0.93" errorPlus="0.11" errorMinus="0.11"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="0.57" errorPlus="0.09" errorMinus="0.09"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.15000000000000002" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/>
+      <measurement value="0.52" errorPlus="0.07" errorMinus="0.07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.05" errorPlus="0.05" errorMinus="0.05"/>
+      <measurement value="0.6" errorPlus="0.08" errorMinus="0.08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+      <measurement value="0.5" errorPlus="0.08" errorMinus="0.08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15000000000000002" errorPlus="0.04999999999999999" errorMinus="0.05000000000000002"/>
+      <measurement value="0.71" errorPlus="0.08" errorMinus="0.08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="0.61" errorPlus="0.08" errorMinus="0.08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/>
+      <measurement value="0.74" errorPlus="0.08" errorMinus="0.08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="1.04" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="1.4" errorPlus="0.12" errorMinus="0.12"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d09-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCHI(P=3))" >
+    <dataPoint>
+      <measurement value="0.105" errorPlus="0.105" errorMinus="0.105"/>
+      <measurement value="0.99" errorPlus="0.07" errorMinus="0.07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.315" errorPlus="0.10499999999999998" errorMinus="0.10500000000000001"/>
+      <measurement value="0.52" errorPlus="0.05" errorMinus="0.05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.525" errorPlus="0.10499999999999998" errorMinus="0.10500000000000004"/>
+      <measurement value="0.31" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.735" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.21" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9450000000000001" errorPlus="0.10499999999999998" errorMinus="0.1050000000000001"/>
+      <measurement value="0.15" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.155" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.08" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.365" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.09" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.575" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.06" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.7799999999999998" errorPlus="0.10000000000000009" errorMinus="0.09999999999999987"/>
+      <measurement value="0.1" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.9849999999999999" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.09" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.195" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.13" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.405" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.18" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.615" errorPlus="0.10499999999999998" errorMinus="0.10500000000000043"/>
+      <measurement value="0.31" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.825" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.65" errorPlus="0.07" errorMinus="0.07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.035" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.88" errorPlus="0.07" errorMinus="0.07"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d10-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCOS(THETA(P=3))" >
+    <dataPoint>
+      <measurement value="-0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.53" errorPlus="0.11" errorMinus="0.11"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.6499999999999999" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="1.08" errorPlus="0.11" errorMinus="0.11"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.55" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/>
+      <measurement value="0.64" errorPlus="0.08" errorMinus="0.08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="0.45" errorPlus="0.07" errorMinus="0.07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.35" errorPlus="0.04999999999999999" errorMinus="0.050000000000000044"/>
+      <measurement value="0.46" errorPlus="0.06" errorMinus="0.06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="0.22" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.15000000000000002" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/>
+      <measurement value="0.41" errorPlus="0.06" errorMinus="0.06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.05" errorPlus="0.05" errorMinus="0.05"/>
+      <measurement value="0.19" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+      <measurement value="0.16" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15000000000000002" errorPlus="0.04999999999999999" errorMinus="0.05000000000000002"/>
+      <measurement value="0.22" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="0.29" errorPlus="0.05" errorMinus="0.05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/>
+      <measurement value="0.37" errorPlus="0.07" errorMinus="0.07"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="0.52" errorPlus="0.08" errorMinus="0.08"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="0.81" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.6499999999999999" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/>
+      <measurement value="0.97" errorPlus="0.13" errorMinus="0.13"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.66" errorPlus="0.11" errorMinus="0.11"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d11-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCHI(P=3))" >
+    <dataPoint>
+      <measurement value="0.105" errorPlus="0.105" errorMinus="0.105"/>
+      <measurement value="0.78" errorPlus="0.05" errorMinus="0.05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.315" errorPlus="0.10499999999999998" errorMinus="0.10500000000000001"/>
+      <measurement value="0.59" errorPlus="0.06" errorMinus="0.06"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.525" errorPlus="0.10499999999999998" errorMinus="0.10500000000000004"/>
+      <measurement value="0.36" errorPlus="0.05" errorMinus="0.05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.735" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.25" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.9450000000000001" errorPlus="0.10499999999999998" errorMinus="0.1050000000000001"/>
+      <measurement value="0.18" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.155" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.1" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.365" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.14" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.575" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.08" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.7799999999999998" errorPlus="0.10000000000000009" errorMinus="0.09999999999999987"/>
+      <measurement value="0.12" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.9849999999999999" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.13" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.195" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.17" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.405" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.23" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.615" errorPlus="0.10499999999999998" errorMinus="0.10500000000000043"/>
+      <measurement value="0.42" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.825" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.52" errorPlus="0.05" errorMinus="0.05"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.035" errorPlus="0.10499999999999998" errorMinus="0.10499999999999998"/>
+      <measurement value="0.69" errorPlus="0.06" errorMinus="0.06"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d12-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCOS(THETA(P=3))" >
+    <dataPoint>
+      <measurement value="-0.875" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="1.53" errorPlus="0.09" errorMinus="0.09"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.625" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="0.27" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.375" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="0.08" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.125" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="0.07" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.125" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="0.07" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.375" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="0.08" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.625" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="0.35" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.875" errorPlus="0.125" errorMinus="0.125"/>
+      <measurement value="1.55" errorPlus="0.09" errorMinus="0.09"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d13-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCHI(P=3))" >
+    <dataPoint>
+      <measurement value="0.195" errorPlus="0.195" errorMinus="0.195"/>
+      <measurement value="0.57" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.19999999999999996" errorMinus="0.20000000000000007"/>
+      <measurement value="0.31" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.985" errorPlus="0.19499999999999995" errorMinus="0.19499999999999995"/>
+      <measurement value="0.23" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.375" errorPlus="0.19500000000000006" errorMinus="0.19500000000000006"/>
+      <measurement value="0.19" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.7650000000000001" errorPlus="0.19499999999999984" errorMinus="0.19500000000000006"/>
+      <measurement value="0.19" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.16" errorPlus="0.19999999999999973" errorMinus="0.20000000000000018"/>
+      <measurement value="0.22" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.5549999999999997" errorPlus="0.19500000000000028" errorMinus="0.19499999999999984"/>
+      <measurement value="0.29" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.9450000000000003" errorPlus="0.19499999999999984" errorMinus="0.19500000000000028"/>
+      <measurement value="0.54" errorPlus="0.05" errorMinus="0.05"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d14-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(3)" >
+    <dataPoint>
+      <measurement value="0.0085" errorPlus="0.0085" errorMinus="0.0085"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0080" errorMinus="0.0080"/>
+      <measurement value="3.7" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0415" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="16.4" errorPlus="1.0" errorMinus="1.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0585" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="14.9" errorPlus="1.0" errorMinus="1.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.07500000000000001" errorPlus="0.007999999999999993" errorMinus="0.008000000000000007"/>
+      <measurement value="7.2" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0915" errorPlus="0.008500000000000008" errorMinus="0.008499999999999994"/>
+      <measurement value="6.1" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.10850000000000001" errorPlus="0.008499999999999994" errorMinus="0.008500000000000008"/>
+      <measurement value="4.5" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.125" errorPlus="0.008000000000000007" errorMinus="0.007999999999999993"/>
+      <measurement value="2.9" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14150000000000001" errorPlus="0.00849999999999998" errorMinus="0.008500000000000008"/>
+      <measurement value="1.4" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1585" errorPlus="0.008500000000000008" errorMinus="0.008500000000000008"/>
+      <measurement value="1.1" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.175" errorPlus="0.008000000000000007" errorMinus="0.00799999999999998"/>
+      <measurement value="0.6" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1915" errorPlus="0.008500000000000008" errorMinus="0.008500000000000008"/>
+      <measurement value="0.0" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d14-x01-y02" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(4)" >
+    <dataPoint>
+      <measurement value="0.0085" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="1.2" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0080" errorMinus="0.0080"/>
+      <measurement value="10.2" errorPlus="0.8" errorMinus="0.8"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0415" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="18.1" errorPlus="1.1" errorMinus="1.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0585" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="12.5" errorPlus="0.9" errorMinus="0.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.07500000000000001" errorPlus="0.007999999999999993" errorMinus="0.008000000000000007"/>
+      <measurement value="7.0" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0915" errorPlus="0.008500000000000008" errorMinus="0.008499999999999994"/>
+      <measurement value="4.0" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.10850000000000001" errorPlus="0.008499999999999994" errorMinus="0.008500000000000008"/>
+      <measurement value="3.4" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.125" errorPlus="0.008000000000000007" errorMinus="0.007999999999999993"/>
+      <measurement value="1.4" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14150000000000001" errorPlus="0.00849999999999998" errorMinus="0.008500000000000008"/>
+      <measurement value="1.1" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1585" errorPlus="0.008500000000000008" errorMinus="0.008500000000000008"/>
+      <measurement value="0.0" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.175" errorPlus="0.008000000000000007" errorMinus="0.00799999999999998"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1915" errorPlus="0.008500000000000008" errorMinus="0.008500000000000008"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d14-x01-y03" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(5)" >
+    <dataPoint>
+      <measurement value="0.0085" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="14.0" errorPlus="0.9" errorMinus="0.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0080" errorMinus="0.0080"/>
+      <measurement value="23.9" errorPlus="1.1" errorMinus="1.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0415" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="11.8" errorPlus="0.8" errorMinus="0.8"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0585" errorPlus="0.0085" errorMinus="0.0085"/>
+      <measurement value="5.6" errorPlus="0.9" errorMinus="0.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.07500000000000001" errorPlus="0.007999999999999993" errorMinus="0.008000000000000007"/>
+      <measurement value="2.9" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0915" errorPlus="0.008500000000000008" errorMinus="0.008499999999999994"/>
+      <measurement value="1.2" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.10850000000000001" errorPlus="0.008499999999999994" errorMinus="0.008500000000000008"/>
+      <measurement value="0.6" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.125" errorPlus="0.008000000000000007" errorMinus="0.007999999999999993"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14150000000000001" errorPlus="0.00849999999999998" errorMinus="0.008500000000000008"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1585" errorPlus="0.008500000000000008" errorMinus="0.008500000000000008"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.175" errorPlus="0.008000000000000007" errorMinus="0.00799999999999998"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1915" errorPlus="0.008500000000000008" errorMinus="0.008500000000000008"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d15-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(3)" >
+    <dataPoint>
+      <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/>
+      <measurement value="1.0" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/>
+      <measurement value="12.3" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/>
+      <measurement value="7.7" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/>
+      <measurement value="3.1" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/>
+      <measurement value="2.5" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/>
+      <measurement value="2.8" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="2.7" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/>
+      <measurement value="2.1" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="1.4" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="0.9" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/>
+      <measurement value="0.6" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="0.1" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d15-x01-y02" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(4)" >
+    <dataPoint>
+      <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/>
+      <measurement value="4.7" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/>
+      <measurement value="12.0" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/>
+      <measurement value="6.8" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/>
+      <measurement value="3.2" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/>
+      <measurement value="3.5" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/>
+      <measurement value="3.1" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="2.5" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/>
+      <measurement value="1.2" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="0.6" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d15-x01-y03" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(5)" >
+    <dataPoint>
+      <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/>
+      <measurement value="7.8" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/>
+      <measurement value="8.6" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/>
+      <measurement value="4.7" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/>
+      <measurement value="4.2" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/>
+      <measurement value="4.1" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/>
+      <measurement value="3.7" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="2.3" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/>
+      <measurement value="1.5" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="0.6" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="0.2" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d16-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(3)" >
+    <dataPoint>
+      <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/>
+      <measurement value="1.7" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/>
+      <measurement value="5.9" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/>
+      <measurement value="2.6" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/>
+      <measurement value="2.4" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/>
+      <measurement value="2.8" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/>
+      <measurement value="4.0" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="5.1" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/>
+      <measurement value="4.1" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="3.0" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="2.3" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/>
+      <measurement value="2.0" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="0.9" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3335" errorPlus="0.013499999999999956" errorMinus="0.013500000000000012"/>
+      <measurement value="0.5" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.36" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+      <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d16-x01-y02" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(4)" >
+    <dataPoint>
+      <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/>
+      <measurement value="6.5" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/>
+      <measurement value="10.8" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/>
+      <measurement value="5.8" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/>
+      <measurement value="3.6" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/>
+      <measurement value="3.7" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/>
+      <measurement value="2.9" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="1.5" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/>
+      <measurement value="1.8" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="0.8" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="0.3" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3335" errorPlus="0.013499999999999956" errorMinus="0.013500000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.36" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d16-x01-y03" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DF(5)" >
+    <dataPoint>
+      <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/>
+      <measurement value="5.9" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/>
+      <measurement value="9.4" errorPlus="0.7" errorMinus="0.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/>
+      <measurement value="6.4" errorPlus="0.6" errorMinus="0.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/>
+      <measurement value="5.0" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/>
+      <measurement value="4.2" errorPlus="0.5" errorMinus="0.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/>
+      <measurement value="2.8" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="1.6" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/>
+      <measurement value="1.2" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/>
+      <measurement value="0.6" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+      <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.3335" errorPlus="0.013499999999999956" errorMinus="0.013500000000000012"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.36" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d17-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DXA" >
+    <dataPoint>
+      <measurement value="0.525" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.7" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="3.1" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.4" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="2.6" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.5" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.0" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="2.3" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="1.4" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="0.3" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d18-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DXA" >
+    <dataPoint>
+      <measurement value="0.525" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.4" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="3.7" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.0" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="3.3" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.4" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="1.5" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="1.4" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="1.2" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="0.4" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d18-x01-y02" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DXC" >
+    <dataPoint>
+      <measurement value="0.525" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="3.4" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="2.9" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.7" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="2.5" errorPlus="0.4" errorMinus="0.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.6" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="2.5" errorPlus="0.3" errorMinus="0.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/>
+      <measurement value="1.8" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/>
+      <measurement value="1.4" errorPlus="0.2" errorMinus="0.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/>
+      <measurement value="0.2" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d19-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCHI(AB)" >
+    <dataPoint>
+      <measurement value="0.195" errorPlus="0.195" errorMinus="0.195"/>
+      <measurement value="0.35" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.19999999999999996" errorMinus="0.20000000000000007"/>
+      <measurement value="0.32" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.985" errorPlus="0.19499999999999995" errorMinus="0.19499999999999995"/>
+      <measurement value="0.31" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.375" errorPlus="0.19500000000000006" errorMinus="0.19500000000000006"/>
+      <measurement value="0.34" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.7650000000000001" errorPlus="0.19499999999999984" errorMinus="0.19500000000000006"/>
+      <measurement value="0.23" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.16" errorPlus="0.19999999999999973" errorMinus="0.20000000000000018"/>
+      <measurement value="0.32" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.5549999999999997" errorPlus="0.19500000000000028" errorMinus="0.19499999999999984"/>
+      <measurement value="0.36" errorPlus="0.02" errorMinus="0.02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.9450000000000003" errorPlus="0.19499999999999984" errorMinus="0.19500000000000028"/>
+      <measurement value="0.34" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d20-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCHI(AB)" >
+    <dataPoint>
+      <measurement value="0.195" errorPlus="0.195" errorMinus="0.195"/>
+      <measurement value="0.35" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.19999999999999996" errorMinus="0.20000000000000007"/>
+      <measurement value="0.36" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.985" errorPlus="0.19499999999999995" errorMinus="0.19499999999999995"/>
+      <measurement value="0.24" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.375" errorPlus="0.19500000000000006" errorMinus="0.19500000000000006"/>
+      <measurement value="0.31" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.7650000000000001" errorPlus="0.19499999999999984" errorMinus="0.19500000000000006"/>
+      <measurement value="0.3" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.16" errorPlus="0.19999999999999973" errorMinus="0.20000000000000018"/>
+      <measurement value="0.3" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.5549999999999997" errorPlus="0.19500000000000028" errorMinus="0.19499999999999984"/>
+      <measurement value="0.3" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.9450000000000003" errorPlus="0.19499999999999984" errorMinus="0.19500000000000028"/>
+      <measurement value="0.38" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d20-x01-y02" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DCHI(CD)" >
+    <dataPoint>
+      <measurement value="0.195" errorPlus="0.195" errorMinus="0.195"/>
+      <measurement value="0.26" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.5900000000000001" errorPlus="0.19999999999999996" errorMinus="0.20000000000000007"/>
+      <measurement value="0.3" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.985" errorPlus="0.19499999999999995" errorMinus="0.19499999999999995"/>
+      <measurement value="0.33" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.375" errorPlus="0.19500000000000006" errorMinus="0.19500000000000006"/>
+      <measurement value="0.35" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.7650000000000001" errorPlus="0.19499999999999984" errorMinus="0.19500000000000006"/>
+      <measurement value="0.41" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.16" errorPlus="0.19999999999999973" errorMinus="0.20000000000000018"/>
+      <measurement value="0.28" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.5549999999999997" errorPlus="0.19500000000000028" errorMinus="0.19499999999999984"/>
+      <measurement value="0.33" errorPlus="0.04" errorMinus="0.04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.9450000000000003" errorPlus="0.19499999999999984" errorMinus="0.19500000000000028"/>
+      <measurement value="0.28" errorPlus="0.03" errorMinus="0.03"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d21-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DFA" >
+    <dataPoint>
+      <measurement value="0.0025" errorPlus="0.0025" errorMinus="0.0025"/>
+      <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0075" errorPlus="0.0025000000000000005" errorMinus="0.0024999999999999996"/>
+      <measurement value="9.2" errorPlus="1.9" errorMinus="1.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0125" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000005"/>
+      <measurement value="21.0" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0175" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000022"/>
+      <measurement value="23.8" errorPlus="1.9" errorMinus="1.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0225" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999988"/>
+      <measurement value="27.6" errorPlus="2.4" errorMinus="2.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0275" errorPlus="0.0024999999999999988" errorMinus="0.0024999999999999988"/>
+      <measurement value="25.7" errorPlus="2.2" errorMinus="2.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0325" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="22.9" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.037500000000000006" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/>
+      <measurement value="20.3" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.042499999999999996" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999953"/>
+      <measurement value="15.9" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0475" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="9.8" errorPlus="1.6" errorMinus="1.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.052500000000000005" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/>
+      <measurement value="7.9" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.057499999999999996" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999953"/>
+      <measurement value="6.0" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0625" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="3.5" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0675" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="4.1" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.07250000000000001" errorPlus="0.0024999999999999883" errorMinus="0.0025000000000000022"/>
+      <measurement value="2.2" errorPlus="1.9" errorMinus="1.9"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d21-x01-y02" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DFB" >
+    <dataPoint>
+      <measurement value="0.0025" errorPlus="0.0025" errorMinus="0.0025"/>
+      <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0075" errorPlus="0.0025000000000000005" errorMinus="0.0024999999999999996"/>
+      <measurement value="24.1" errorPlus="2.1" errorMinus="2.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0125" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000005"/>
+      <measurement value="53.6" errorPlus="3.0" errorMinus="3.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0175" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000022"/>
+      <measurement value="45.6" errorPlus="2.9" errorMinus="2.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0225" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999988"/>
+      <measurement value="29.4" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0275" errorPlus="0.0024999999999999988" errorMinus="0.0024999999999999988"/>
+      <measurement value="17.6" errorPlus="1.7" errorMinus="1.7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0325" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="11.9" errorPlus="1.6" errorMinus="1.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.037500000000000006" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/>
+      <measurement value="3.9" errorPlus="1.6" errorMinus="1.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.042499999999999996" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999953"/>
+      <measurement value="3.6" errorPlus="2.2" errorMinus="2.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0475" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="2.6" errorPlus="1.6" errorMinus="1.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.052500000000000005" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/>
+      <measurement value="1.9" errorPlus="1.6" errorMinus="1.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.057499999999999996" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999953"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0625" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0675" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.07250000000000001" errorPlus="0.0024999999999999883" errorMinus="0.0025000000000000022"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d22-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DFA" >
+    <dataPoint>
+      <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/>
+      <measurement value="2.8" errorPlus="1.0" errorMinus="1.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/>
+      <measurement value="16.6" errorPlus="1.3" errorMinus="1.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/>
+      <measurement value="23.2" errorPlus="1.9" errorMinus="1.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/>
+      <measurement value="18.2" errorPlus="1.4" errorMinus="1.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/>
+      <measurement value="11.6" errorPlus="1.2" errorMinus="1.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/>
+      <measurement value="7.6" errorPlus="1.4" errorMinus="1.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.065" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="5.8" errorPlus="1.1" errorMinus="1.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.07500000000000001" errorPlus="0.0049999999999999906" errorMinus="0.0050000000000000044"/>
+      <measurement value="4.6" errorPlus="1.2" errorMinus="1.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.08499999999999999" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999906"/>
+      <measurement value="4.8" errorPlus="0.9" errorMinus="0.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.095" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="2.4" errorPlus="1.3" errorMinus="1.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.10500000000000001" errorPlus="0.0049999999999999906" errorMinus="0.0050000000000000044"/>
+      <measurement value="1.4" errorPlus="0.9" errorMinus="0.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.11499999999999999" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999906"/>
+      <measurement value="1.0" errorPlus="1.0" errorMinus="1.0"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d23-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DFB" >
+    <dataPoint>
+      <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/>
+      <measurement value="9.8" errorPlus="1.2" errorMinus="1.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/>
+      <measurement value="30.1" errorPlus="2.0" errorMinus="2.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/>
+      <measurement value="20.1" errorPlus="1.5" errorMinus="1.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/>
+      <measurement value="11.7" errorPlus="1.3" errorMinus="1.3"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/>
+      <measurement value="6.6" errorPlus="1.0" errorMinus="1.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/>
+      <measurement value="5.2" errorPlus="1.0" errorMinus="1.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.065" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="4.7" errorPlus="1.1" errorMinus="1.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.07500000000000001" errorPlus="0.0049999999999999906" errorMinus="0.0050000000000000044"/>
+      <measurement value="3.5" errorPlus="0.9" errorMinus="0.9"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.08499999999999999" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999906"/>
+      <measurement value="3.1" errorPlus="1.2" errorMinus="1.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.095" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/>
+      <measurement value="2.4" errorPlus="1.1" errorMinus="1.1"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d24-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DFC" >
+    <dataPoint>
+      <measurement value="0.0025" errorPlus="0.0025" errorMinus="0.0025"/>
+      <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0075" errorPlus="0.0025000000000000005" errorMinus="0.0024999999999999996"/>
+      <measurement value="29.4" errorPlus="2.6" errorMinus="2.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0125" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000005"/>
+      <measurement value="46.5" errorPlus="3.5" errorMinus="3.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0175" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000022"/>
+      <measurement value="43.4" errorPlus="3.1" errorMinus="3.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0225" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999988"/>
+      <measurement value="27.6" errorPlus="2.6" errorMinus="2.6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0275" errorPlus="0.0024999999999999988" errorMinus="0.0024999999999999988"/>
+      <measurement value="21.1" errorPlus="2.2" errorMinus="2.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0325" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="11.8" errorPlus="2.4" errorMinus="2.4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.037500000000000006" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/>
+      <measurement value="6.6" errorPlus="2.0" errorMinus="2.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.042499999999999996" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999953"/>
+      <measurement value="3.9" errorPlus="2.2" errorMinus="2.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0475" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/>
+      <measurement value="3.1" errorPlus="2.0" errorMinus="2.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.052500000000000005" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/>
+      <measurement value="2.2" errorPlus="2.0" errorMinus="2.0"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d25-x01-y01" dimension="2" path="/REF/CDF_1996_S3349578" title="(1/N)*D(N)/DFD" >
+    <dataPoint>
+      <measurement value="0.0025" errorPlus="0.0025" errorMinus="0.0025"/>
+      <measurement value="3.5" errorPlus="2.2" errorMinus="2.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0075" errorPlus="0.0025000000000000005" errorMinus="0.0024999999999999996"/>
+      <measurement value="56.1" errorPlus="3.5" errorMinus="3.5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0125" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000005"/>
+      <measurement value="73.7" errorPlus="4.2" errorMinus="4.2"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0175" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000022"/>
+      <measurement value="38.2" errorPlus="3.1" errorMinus="3.1"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0225" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999988"/>
+      <measurement value="18.0" errorPlus="2.0" errorMinus="2.0"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.0275" errorPlus="0.0024999999999999988" errorMinus="0.0024999999999999988"/>
+      <measurement value="7.9" errorPlus="2.6" errorMinus="2.6"/>
+    </dataPoint>
+  </dataPointSet>
+</aida>

Modified: trunk/data/refdata/Makefile.am
==============================================================================
--- trunk/data/refdata/Makefile.am	Wed Oct 21 21:29:18 2009	(r1948)
+++ trunk/data/refdata/Makefile.am	Wed Oct 21 21:30:24 2009	(r1949)
@@ -33,6 +33,7 @@
   CDF_1990_S2089246.aida         \
   CDF_1994_S2952106.aida         \
   CDF_1996_S3108457.aida         \
+  CDF_1996_S3349578.aida         \
   CDF_2000_S4155203.aida         \
   CDF_2001_S4751469.aida         \
   CDF_2002_S4796047.aida         \

Added: trunk/src/Analyses/CDF_1996_S3349578.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/CDF_1996_S3349578.cc	Wed Oct 21 21:30:24 2009	(r1949)
@@ -0,0 +1,438 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/FastJets.hh"
+
+namespace Rivet {
+
+
+  class CDF_1996_S3349578 : public Analysis {
+  public:
+
+    /// @name Constructors etc.
+    //@{
+
+    /// Constructor
+    CDF_1996_S3349578()
+      : Analysis("CDF_1996_S3349578") 
+    {
+      setBeams(PROTON, ANTIPROTON);
+      setNeedsCrossSection(true);
+    }
+
+    //@}
+
+
+  public:
+
+    /// @name Analysis methods
+    //@{
+
+    /// Book histograms and initialise projections before the run
+    void init() {
+
+      /// Initialise and register projections here
+      const FinalState fs(-4.2, 4.2);
+      addProjection(FastJets(fs, FastJets::CDFJETCLU, 0.7), "Jets");
+
+      /// Book histograms here, e.g.:
+      _h_3_mNJ = bookHistogram1D(1, 1, 1);
+      _h_3_X3 = bookHistogram1D(2, 1, 1);
+      _h_3_X4 = bookHistogram1D(3, 1, 1);
+      _h_3_costheta3 = bookHistogram1D(8, 1, 1);
+      _h_3_psi3 = bookHistogram1D(9, 1, 1);
+      _h_3_f3 = bookHistogram1D(14, 1, 1);
+      _h_3_f4 = bookHistogram1D(14, 1, 2);
+      _h_3_f5 = bookHistogram1D(14, 1, 3);
+      
+      _h_4_mNJ = bookHistogram1D(1, 1, 2);
+      _h_4_X3 = bookHistogram1D(4, 1, 1);
+      _h_4_X4 = bookHistogram1D(5, 1, 1);
+      _h_4_costheta3 = bookHistogram1D(10, 1, 1);
+      _h_4_psi3 = bookHistogram1D(11, 1, 1);
+      _h_4_f3 = bookHistogram1D(15, 1, 1);
+      _h_4_f4 = bookHistogram1D(15, 1, 2);
+      _h_4_f5 = bookHistogram1D(15, 1, 3);
+      _h_4_XA = bookHistogram1D(17, 1, 1);
+      _h_4_psiAB = bookHistogram1D(19, 1, 1);
+      _h_4_fA = bookHistogram1D(21, 1, 1);
+      _h_4_fB = bookHistogram1D(21, 1, 2);
+      
+      _h_5_mNJ = bookHistogram1D(1, 1, 3);
+      _h_5_X3 = bookHistogram1D(6, 1, 1);
+      _h_5_X4 = bookHistogram1D(7, 1, 1);
+      _h_5_costheta3 = bookHistogram1D(12, 1, 1);
+      _h_5_psi3 = bookHistogram1D(13, 1, 1);
+      _h_5_f3 = bookHistogram1D(16, 1, 1);
+      _h_5_f4 = bookHistogram1D(16, 1, 2);
+      _h_5_f5 = bookHistogram1D(16, 1, 3);
+      _h_5_XA = bookHistogram1D(18, 1, 1);
+      _h_5_XC = bookHistogram1D(18, 1, 2);
+      _h_5_psiAB = bookHistogram1D(20, 1, 1);
+      _h_5_psiCD = bookHistogram1D(20, 1, 2);
+      _h_5_fA = bookHistogram1D(22, 1, 1);
+      _h_5_fB = bookHistogram1D(23, 1, 1);
+      _h_5_fC = bookHistogram1D(24, 1, 1);
+      _h_5_fD = bookHistogram1D(25, 1, 1);
+      
+    }
+
+
+    /// Perform the per-event analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+
+      /// Do the event by event analysis here
+      Jets jets;
+      double sumEt = 0.0;
+      FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0);
+      foreach (const Jet& jet, applyProjection<FastJets>(event, "Jets").jetsByEt()) {
+        double Et = jet.momentum().Et();
+        if (Et > 20.0*GeV) {
+          jets.push_back(jet);
+          sumEt += Et;
+          jetsystem += jet.momentum();
+        }
+      }
+      /// @todo include gaussian jet energy resolution smearing?
+      
+      if (jets.size() < 3) {
+        vetoEvent;
+      }
+      
+      if (sumEt < 420.0*GeV) {
+        vetoEvent;
+      }
+      
+      if (jets.size() > 2) _threeJetAnalysis(jets, weight);
+      if (jets.size() > 3) _fourJetAnalysis(jets, weight);
+      if (jets.size() > 4) _fiveJetAnalysis(jets, weight);
+    }
+    
+    void _threeJetAnalysis(const Jets& jets, const double& weight) {
+      getLog() << Log::DEBUG << "3 jet analysis" << std::endl;
+      FourMomentum jjj(jets[0].momentum()+jets[1].momentum()+jets[2].momentum());
+      const double m3J = jjj.mass();
+      if (m3J<600*GeV) {
+        return;
+      }
+    
+      LorentzTransform cms_boost(-jjj.boostVector());
+      vector<FourMomentum> jets_boosted;
+      foreach (Jet jet, jets) {
+        jets_boosted.push_back(cms_boost.transform(jet.momentum()));
+      }
+      std::sort(jets_boosted.begin(), jets_boosted.end(), FourMomentum::byEDescending());
+      FourMomentum p3(jets_boosted[0]);
+      FourMomentum p4(jets_boosted[1]);
+      FourMomentum p5(jets_boosted[2]);
+      
+      double costheta3 = cos(p3.theta());
+      if (fabs(costheta3)>0.6) {
+        return;
+      }
+      
+      double X3 = 2.0*p3.E()/m3J;
+      if (X3>0.9) {
+        return;
+      }
+      
+      
+      // fill histograms
+      const double X4 = 2.0*p4.E()/m3J;
+      Vector3 beam1(0.0, 0.0, 1.0);
+      Vector3 p1xp3 = beam1.cross(p3.vector3());
+      Vector3 p4xp5 = p4.vector3().cross(p5.vector3());
+      const double cospsi3 = p1xp3.dot(p4xp5)/p1xp3.mod()/p4xp5.mod();
+      const double f3 = p3.mass()/m3J;
+      const double f4 = p4.mass()/m3J;
+      const double f5 = p5.mass()/m3J;
+      
+      _h_3_mNJ->fill(m3J, weight);
+      _h_3_X3->fill(X3, weight);
+      _h_3_X4->fill(X4, weight);
+      _h_3_costheta3->fill(costheta3, weight);
+      _h_3_psi3->fill(mapAngle0ToPi(acos(cospsi3)), weight);
+      _h_3_f3->fill(f3, weight);
+      _h_3_f4->fill(f4, weight);
+      _h_3_f5->fill(f5, weight);
+      
+    }
+
+    void _fourJetAnalysis(const Jets& jets, const double& weight) {
+      getLog() << Log::DEBUG << "4 jet analysis" << std::endl;
+      FourMomentum jjjj(0.0, 0.0, 0.0, 0.0);
+      vector<FourMomentum> jetmoms;
+      for (size_t i=0; i<4; ++i) {
+        jetmoms.push_back(jets[i].momentum());
+        jjjj += jets[i].momentum();
+      }
+      const double m4J = jjjj.mass();
+      if (m4J < 650*GeV) return;
+      
+      FourMomentum pA, pB;
+      vector<FourMomentum> jetmoms3(_reduce(jetmoms, pA, pB));
+      LorentzTransform cms_boost(-jjjj.boostVector());
+      vector<FourMomentum> jetmoms3_boosted;
+      foreach (FourMomentum mom, jetmoms3) {
+        jetmoms3_boosted.push_back(cms_boost.transform(mom));
+      }
+      pA = cms_boost.transform(pA);
+      pB = cms_boost.transform(pB);
+      
+      sort(jetmoms3_boosted.begin(), jetmoms3_boosted.end(), FourMomentum::byEDescending());
+      if (pB.E()>pA.E()) std::swap(pA, pB);
+      FourMomentum p3(jetmoms3_boosted[0]);
+      FourMomentum p4(jetmoms3_boosted[1]);
+      FourMomentum p5(jetmoms3_boosted[2]);
+      
+      const double costheta3 = cos(p3.theta());
+      if (fabs(costheta3)>0.8) {
+        return;
+      }
+      
+      const double X3 = 2.0*p3.E()/m4J;
+      if (X3>0.9) {
+        return;
+      }
+      
+      // fill histograms
+      const double X4 = 2.0*p4.E()/m4J;
+      Vector3 beam1(0.0, 0.0, 1.0);
+      Vector3 p1xp3 = beam1.cross(p3.vector3());
+      Vector3 p4xp5 = p4.vector3().cross(p5.vector3());
+      const double cospsi3 = p1xp3.dot(p4xp5)/p1xp3.mod()/p4xp5.mod();
+      const double f3 = p3.mass()/m4J;
+      const double f4 = p4.mass()/m4J;
+      const double f5 = p5.mass()/m4J;
+      const double fA = pA.mass()/m4J;
+      const double fB = pB.mass()/m4J;
+      const double XA = pA.E()/(pA.E()+pB.E());
+      FourMomentum pAB = pA+pB;
+      Vector3 pABxp1 = pAB.vector3().cross(beam1);
+      Vector3 pAxpB = pA.vector3().cross(pB.vector3());
+      const double cospsiAB = pAxpB.dot(pABxp1)/pAxpB.mod()/pABxp1.mod();
+      
+      _h_4_mNJ->fill(m4J, weight);
+      _h_4_X3->fill(X3, weight);
+      _h_4_X4->fill(X4, weight);
+      _h_4_costheta3->fill(costheta3, weight);
+      _h_4_psi3->fill(mapAngle0ToPi(acos(cospsi3)), weight);
+      _h_4_f3->fill(f3, weight);
+      _h_4_f4->fill(f4, weight);
+      _h_4_f5->fill(f5, weight);
+      _h_4_XA->fill(XA, weight);
+      _h_4_psiAB->fill(mapAngle0ToPi(acos(cospsiAB)), weight);
+      _h_4_fA->fill(fA, weight);
+      _h_4_fB->fill(fB, weight);
+    }
+      
+      
+    void _fiveJetAnalysis(const Jets& jets, const double& weight) {
+      getLog() << Log::DEBUG << "5 jet analysis" << std::endl;
+      FourMomentum jjjjj(0.0, 0.0, 0.0, 0.0);
+      vector<FourMomentum> jetmoms;
+      for (size_t i=0; i<5; ++i) {
+        jetmoms.push_back(jets[i].momentum());
+        jjjjj += jets[i].momentum();
+      }
+      const double m5J = jjjjj.mass();
+      if (m5J < 750*GeV) return;
+      
+      FourMomentum pA, pB, pC, pD;
+      vector<FourMomentum> jetmoms4(_reduce(jetmoms, pC, pD));
+      vector<FourMomentum> jetmoms3(_reduce(jetmoms4, pA, pB));
+      
+      LorentzTransform cms_boost(-jjjjj.boostVector());
+      vector<FourMomentum> jetmoms3_boosted;
+      foreach (FourMomentum mom, jetmoms3) {
+        jetmoms3_boosted.push_back(cms_boost.transform(mom));
+      }
+      pA = cms_boost.transform(pA);
+      pB = cms_boost.transform(pB);
+      pC = cms_boost.transform(pC);
+      pD = cms_boost.transform(pD);
+      
+      sort(jetmoms3_boosted.begin(), jetmoms3_boosted.end(), FourMomentum::byEDescending());
+      if (pB.E()>pA.E()) std::swap(pA, pB);
+      if (pD.E()>pC.E()) std::swap(pD, pC);
+      FourMomentum p3(jetmoms3_boosted[0]);
+      FourMomentum p4(jetmoms3_boosted[1]);
+      FourMomentum p5(jetmoms3_boosted[2]);
+      
+      // fill histograms
+      const double costheta3 = cos(p3.theta());
+      const double X3 = 2.0*p3.E()/m5J;
+      const double X4 = 2.0*p4.E()/m5J;
+      Vector3 beam1(0.0, 0.0, 1.0);
+      Vector3 p1xp3 = beam1.cross(p3.vector3());
+      Vector3 p4xp5 = p4.vector3().cross(p5.vector3());
+      const double cospsi3 = p1xp3.dot(p4xp5)/p1xp3.mod()/p4xp5.mod();
+      const double f3 = p3.mass()/m5J;
+      const double f4 = p4.mass()/m5J;
+      const double f5 = p5.mass()/m5J;
+      const double fA = pA.mass()/m5J;
+      const double fB = pB.mass()/m5J;
+      const double XA = pA.E()/(pA.E()+pB.E());
+      FourMomentum pAB = pA+pB;
+      Vector3 pABxp1 = pAB.vector3().cross(beam1);
+      Vector3 pAxpB = pA.vector3().cross(pB.vector3());
+      const double cospsiAB = pAxpB.dot(pABxp1)/pAxpB.mod()/pABxp1.mod();
+      const double fC = pC.mass()/m5J;
+      const double fD = pD.mass()/m5J;
+      const double XC = pC.E()/(pC.E()+pD.E());
+      FourMomentum pCD = pC+pD;
+      Vector3 pCDxp1 = pCD.vector3().cross(beam1);
+      Vector3 pCxpD = pC.vector3().cross(pD.vector3());
+      const double cospsiCD = pCxpD.dot(pCDxp1)/pCxpD.mod()/pCDxp1.mod();
+      
+      _h_5_mNJ->fill(m5J, weight);
+      _h_5_X3->fill(X3, weight);
+      _h_5_X4->fill(X4, weight);
+      _h_5_costheta3->fill(costheta3, weight);
+      _h_5_psi3->fill(mapAngle0ToPi(acos(cospsi3)), weight);
+      _h_5_f3->fill(f3, weight);
+      _h_5_f4->fill(f4, weight);
+      _h_5_f5->fill(f5, weight);
+      _h_5_XA->fill(XA, weight);
+      _h_5_psiAB->fill(mapAngle0ToPi(acos(cospsiAB)), weight);
+      _h_5_fA->fill(fA, weight);
+      _h_5_fB->fill(fB, weight);
+      _h_5_XC->fill(XC, weight);
+      _h_5_psiCD->fill(mapAngle0ToPi(acos(cospsiCD)), weight);
+      _h_5_fC->fill(fC, weight);
+      _h_5_fD->fill(fD, weight);
+    }
+      
+      
+    /// Normalise histograms etc., after the run
+    void finalize() {
+      
+      /// Normalise, scale and otherwise manipulate histograms here
+      scale(_h_3_mNJ, crossSection()/sumOfWeights());
+      scale(_h_3_X3, crossSection()/sumOfWeights());
+      scale(_h_3_X4, crossSection()/sumOfWeights());
+      scale(_h_3_costheta3, crossSection()/sumOfWeights());
+      scale(_h_3_psi3, crossSection()/sumOfWeights());
+      scale(_h_3_f3, crossSection()/sumOfWeights());
+      scale(_h_3_f4, crossSection()/sumOfWeights());
+      scale(_h_3_f5, crossSection()/sumOfWeights());
+      
+      scale(_h_4_mNJ, crossSection()/sumOfWeights());
+      scale(_h_4_X3, crossSection()/sumOfWeights());
+      scale(_h_4_X4, crossSection()/sumOfWeights());
+      scale(_h_4_costheta3, crossSection()/sumOfWeights());
+      scale(_h_4_psi3, crossSection()/sumOfWeights());
+      scale(_h_4_f3, crossSection()/sumOfWeights());
+      scale(_h_4_f4, crossSection()/sumOfWeights());
+      scale(_h_4_f5, crossSection()/sumOfWeights());
+      scale(_h_4_XA, crossSection()/sumOfWeights());
+      scale(_h_4_psiAB, crossSection()/sumOfWeights());
+      scale(_h_4_fA, crossSection()/sumOfWeights());
+      scale(_h_4_fB, crossSection()/sumOfWeights());
+      
+      scale(_h_5_mNJ, crossSection()/sumOfWeights());
+      scale(_h_5_X3, crossSection()/sumOfWeights());
+      scale(_h_5_X4, crossSection()/sumOfWeights());
+      scale(_h_5_costheta3, crossSection()/sumOfWeights());
+      scale(_h_5_psi3, crossSection()/sumOfWeights());
+      scale(_h_5_f3, crossSection()/sumOfWeights());
+      scale(_h_5_f4, crossSection()/sumOfWeights());
+      scale(_h_5_f5, crossSection()/sumOfWeights());
+      scale(_h_5_XA, crossSection()/sumOfWeights());
+      scale(_h_5_XC, crossSection()/sumOfWeights());
+      scale(_h_5_psiAB, crossSection()/sumOfWeights());
+      scale(_h_5_psiCD, crossSection()/sumOfWeights());
+      scale(_h_5_fA, crossSection()/sumOfWeights());
+      scale(_h_5_fB, crossSection()/sumOfWeights());
+      scale(_h_5_fC, crossSection()/sumOfWeights());
+      scale(_h_5_fD, crossSection()/sumOfWeights());
+      
+    }
+
+    //@}
+
+
+  private:
+    vector<FourMomentum> _reduce(const vector<FourMomentum>& jets,
+                                 FourMomentum& combined1,
+                                 FourMomentum& combined2) {
+      double minMass2 = 1e9;
+      FourMomentum combined;
+      size_t idx1(jets.size()), idx2(jets.size());
+      for (size_t i=0; i<jets.size(); ++i) {
+        for (size_t j=i+1; j<jets.size(); ++j) {
+          double mass2 = FourMomentum(jets[i]+jets[j]).mass2();
+          if (mass2<minMass2) {
+            combined = jets[i]+jets[j];
+            idx1=i;
+            idx2=j;
+          }
+        }
+      }
+      vector<FourMomentum> newjets;
+      for (size_t i=0; i<jets.size(); ++i) {
+        if (i!=idx1 && i!=idx2) newjets.push_back(jets[i]);
+      }
+      newjets.push_back(combined);
+      combined1 = jets[idx1];
+      combined2 = jets[idx2];
+      return newjets;
+    }
+    
+
+  private:
+
+    /// @name Histograms
+    //@{
+    AIDA::IHistogram1D *_h_3_mNJ;
+    AIDA::IHistogram1D *_h_3_X3;
+    AIDA::IHistogram1D *_h_3_X4;
+    AIDA::IHistogram1D *_h_3_costheta3;
+    AIDA::IHistogram1D *_h_3_psi3;
+    AIDA::IHistogram1D *_h_3_f3;
+    AIDA::IHistogram1D *_h_3_f4;
+    AIDA::IHistogram1D *_h_3_f5;
+    
+    AIDA::IHistogram1D *_h_4_mNJ;
+    AIDA::IHistogram1D *_h_4_X3;
+    AIDA::IHistogram1D *_h_4_X4;
+    AIDA::IHistogram1D *_h_4_costheta3;
+    AIDA::IHistogram1D *_h_4_psi3;
+    AIDA::IHistogram1D *_h_4_f3;
+    AIDA::IHistogram1D *_h_4_f4;
+    AIDA::IHistogram1D *_h_4_f5;
+    AIDA::IHistogram1D *_h_4_XA;
+    AIDA::IHistogram1D *_h_4_psiAB;
+    AIDA::IHistogram1D *_h_4_fA;
+    AIDA::IHistogram1D *_h_4_fB;
+    
+    AIDA::IHistogram1D *_h_5_mNJ;
+    AIDA::IHistogram1D *_h_5_X3;
+    AIDA::IHistogram1D *_h_5_X4;
+    AIDA::IHistogram1D *_h_5_costheta3;
+    AIDA::IHistogram1D *_h_5_psi3;
+    AIDA::IHistogram1D *_h_5_f3;
+    AIDA::IHistogram1D *_h_5_f4;
+    AIDA::IHistogram1D *_h_5_f5;
+    AIDA::IHistogram1D *_h_5_XA;
+    AIDA::IHistogram1D *_h_5_XC;
+    AIDA::IHistogram1D *_h_5_psiAB;
+    AIDA::IHistogram1D *_h_5_psiCD;
+    AIDA::IHistogram1D *_h_5_fA;
+    AIDA::IHistogram1D *_h_5_fB;
+    AIDA::IHistogram1D *_h_5_fC;
+    AIDA::IHistogram1D *_h_5_fD;
+    //@}
+
+  };
+
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<CDF_1996_S3349578> plugin_CDF_1996_S3349578;
+
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Wed Oct 21 21:29:18 2009	(r1948)
+++ trunk/src/Analyses/Makefile.am	Wed Oct 21 21:30:24 2009	(r1949)
@@ -15,6 +15,7 @@
     CDF_1990_S2089246.cc \
     CDF_1994_S2952106.cc \
     CDF_1996_S3108457.cc \
+    CDF_1996_S3349578.cc \
     CDF_2000_S4155203.cc \
     CDF_2001_S4751469.cc \
     CDF_2002_S4796047.cc \


More information about the Rivet-svn mailing list