|
[Rivet-svn] r2294 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgTue Mar 2 16:06:02 GMT 2010
Author: fsiegert Date: Tue Mar 2 16:06:02 2010 New Revision: 2294 Log: Add CDF_1997_S3541940 analysis and fix the similar CDF_1996_S3541940 one. Both are still unvalidated though. Added: trunk/data/anainfo/CDF_1997_S3541940.info trunk/data/plotinfo/CDF_1997_S3541940.plot trunk/data/refdata/CDF_1997_S3541940.aida trunk/src/Analyses/CDF_1997_S3541940.cc Modified: trunk/src/Analyses/CDF_1996_S3349578.cc trunk/src/Analyses/Makefile.am Added: trunk/data/anainfo/CDF_1997_S3541940.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_1997_S3541940.info Tue Mar 2 16:06:02 2010 (r2294) @@ -0,0 +1,20 @@ +Name: CDF_1997_S3541940 +Year: 1997 +Summary: Properties of six jet events with large six jet mass +Experiment: CDF +Collider: Tevatron Run 1 +SpiresID: 3541940 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - Phys.Rev.D56:2532-2543,1997 + - doi:10.1103/PhysRevD.56.2532 + - http://lss.fnal.gov/archive/1997/pub/Pub-97-093-E.pdf +RunInfo: + Pure QCD events without underlying event. +NumEvents: 1000000 +PtCuts: [0] +Description: + Multijet distributions corresponding to 20 variables that span the 6-body + parameter space in inclusive 6-jet events. Added: trunk/data/plotinfo/CDF_1997_S3541940.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_1997_S3541940.plot Tue Mar 2 16:06:02 2010 (r2294) @@ -0,0 +1,138 @@ +# BEGIN PLOT /CDF_1997_S3541940/d01-x01-y01 +Title=Multijet mass +XLabel=$m_{6J}$ [GeV] +YLabel=$1/N \mathrm{d}N/\mathrm{d}m$ [1/GeV] +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d02-x01-y01 +Title=Dalitz distribution +XLabel=$X_3^{'''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}X_3^{'''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d03-x01-y01 +Title=Dalitz distribution +XLabel=$X_4^{'''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}X_4^{'''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d04-x01-y01 +Title=Leading jet angle +XLabel=$|\cos \theta_{3^{'''}}|$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos\theta_{3^{'''}}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d05-x01-y01 +Title=Angular distribution +XLabel=$\Psi_{3^{'''}}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi_{3^{'''}}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d06-x01-y01 +Title=Single-jet mass fraction +XLabel=$f_3^{'''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_3^{'''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d06-x01-y02 +Title=Single-jet mass fraction +XLabel=$f_4^{'''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_4^{'''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d06-x01-y03 +Title=Single-jet mass fraction +XLabel=$f_5^{'''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_5^{'''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d07-x01-y01 +Title=Two-body energy sharing +XLabel=$X_A^{''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}X_A^{''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d08-x01-y01 +Title=Two-body energy sharing +XLabel=$X_C^{'}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}X_C^{'}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d09-x01-y01 +Title=Two-body energy sharing +XLabel=$X_E$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}X_E$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d10-x01-y01 +Title=Two-body angular distribution +XLabel=$\Psi_{A^{''}B^{''}}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi_{A^{''}B^{''}}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d11-x01-y01 +Title=Two-body angular distribution +XLabel=$\Psi_{C^{'}D^{'}}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi_{C^{'}D^{'}}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d12-x01-y01 +Title=Two-body angular distribution +XLabel=$\Psi_{EF}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}\Psi_{EF}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d13-x01-y01 +Title=Single-body mass fraction +XLabel=$f_A^{''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_A^{''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d14-x01-y01 +Title=Single-body mass fraction +XLabel=$f_B^{''}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_B^{''}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d15-x01-y01 +Title=Single-body mass fraction +XLabel=$f_C^{'}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_C^{'}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d16-x01-y01 +Title=Single-body mass fraction +XLabel=$f_D^{'}$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_D^{'}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d17-x01-y01 +Title=Single-body mass fraction +XLabel=$f_E$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_E$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1997_S3541940/d18-x01-y01 +Title=Single-body mass fraction +XLabel=$f_F$ +YLabel=$1/N \mathrm{d}N/\mathrm{d}f_F$ +LogY=0 +# END PLOT Added: trunk/data/refdata/CDF_1997_S3541940.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_1997_S3541940.aida Tue Mar 2 16:06:02 2010 (r2294) @@ -0,0 +1,889 @@ +<?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_1997_S3541940" title="(1/N)*D(N)/DM IN PER 40 GEV" > + <dataPoint> + <measurement value="540.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.00692" errorPlus="9.0E-4" errorMinus="9.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="580.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.00549" errorPlus="6.8E-4" errorMinus="6.8E-4"/> + </dataPoint> + <dataPoint> + <measurement value="620.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.00418" errorPlus="5.4E-4" errorMinus="5.4E-4"/> + </dataPoint> + <dataPoint> + <measurement value="660.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.00267" errorPlus="3.4E-4" errorMinus="3.4E-4"/> + </dataPoint> + <dataPoint> + <measurement value="700.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.00187" errorPlus="2.3E-4" errorMinus="2.3E-4"/> + </dataPoint> + <dataPoint> + <measurement value="740.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="0.00143" errorPlus="1.9E-4" errorMinus="1.9E-4"/> + </dataPoint> + <dataPoint> + <measurement value="780.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="9.6E-4" errorPlus="1.24E-4" errorMinus="1.24E-4"/> + </dataPoint> + <dataPoint> + <measurement value="820.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="4.48E-4" errorPlus="8.0E-5" errorMinus="8.0E-5"/> + </dataPoint> + <dataPoint> + <measurement value="860.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="2.53E-4" errorPlus="6.7E-5" errorMinus="6.7E-5"/> + </dataPoint> + <dataPoint> + <measurement value="900.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="1.35E-4" errorPlus="5.0E-5" errorMinus="5.0E-5"/> + </dataPoint> + <dataPoint> + <measurement value="940.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="1.16E-4" errorPlus="4.7E-5" errorMinus="4.7E-5"/> + </dataPoint> + <dataPoint> + <measurement value="980.0" errorPlus="20.0" errorMinus="20.0"/> + <measurement value="7.61E-5" errorPlus="3.37E-5" errorMinus="3.37E-5"/> + </dataPoint> + <dataPoint> + <measurement value="1020.0" errorPlus="60.0" errorMinus="20.0"/> + <measurement value="2.86E-5" errorPlus="1.63E-5" errorMinus="1.63E-5"/> + </dataPoint> + <dataPoint> + <measurement value="1171.0" errorPlus="69.0" errorMinus="91.0"/> + <measurement value="1.44E-5" errorPlus="8.1E-6" errorMinus="8.1E-6"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DX(P=3,Q=THREE-JET)" > + <dataPoint> + <measurement value="0.666" errorPlus="0.015000000000000013" errorMinus="0.016000000000000014"/> + <measurement value="0.32" errorPlus="0.15" errorMinus="0.15"/> + </dataPoint> + <dataPoint> + <measurement value="0.697" errorPlus="0.016000000000000014" errorMinus="0.015999999999999903"/> + <measurement value="1.82" errorPlus="0.22" errorMinus="0.22"/> + </dataPoint> + <dataPoint> + <measurement value="0.728" errorPlus="0.016000000000000014" errorMinus="0.015000000000000013"/> + <measurement value="3.15" errorPlus="0.3" errorMinus="0.3"/> + </dataPoint> + <dataPoint> + <measurement value="0.76" errorPlus="0.015000000000000013" errorMinus="0.016000000000000014"/> + <measurement value="4.19" errorPlus="0.3" errorMinus="0.3"/> + </dataPoint> + <dataPoint> + <measurement value="0.791" errorPlus="0.015000000000000013" errorMinus="0.016000000000000014"/> + <measurement value="5.21" errorPlus="0.39" errorMinus="0.39"/> + </dataPoint> + <dataPoint> + <measurement value="0.822" errorPlus="0.016000000000000014" errorMinus="0.015999999999999903"/> + <measurement value="5.52" errorPlus="0.38" errorMinus="0.38"/> + </dataPoint> + <dataPoint> + <measurement value="0.853" errorPlus="0.016000000000000014" errorMinus="0.015000000000000013"/> + <measurement value="6.37" errorPlus="0.38" errorMinus="0.38"/> + </dataPoint> + <dataPoint> + <measurement value="0.884" errorPlus="0.016000000000000014" errorMinus="0.015000000000000013"/> + <measurement value="5.43" errorPlus="0.36" errorMinus="0.36"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DX(P=4,Q=THREE-JET)" > + <dataPoint> + <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.36" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="4.36" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="6.3" errorPlus="0.33" errorMinus="0.33"/> + </dataPoint> + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="4.34" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.33" errorPlus="0.18" errorMinus="0.18"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.07" errorPlus="0.18" errorMinus="0.18"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.18" errorPlus="0.19" errorMinus="0.19"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DCOS(THETA(P=3,Q=THREE-JET))" > + <dataPoint> + <measurement value="-0.8500000000000001" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/> + <measurement value="2.201" errorPlus="0.059" errorMinus="0.059"/> + </dataPoint> + <dataPoint> + <measurement value="-0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="1.126" errorPlus="0.106" errorMinus="0.106"/> + </dataPoint> + <dataPoint> + <measurement value="-0.6499999999999999" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="0.653" errorPlus="0.078" errorMinus="0.078"/> + </dataPoint> + <dataPoint> + <measurement value="-0.55" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/> + <measurement value="0.342" errorPlus="0.053" errorMinus="0.053"/> + </dataPoint> + <dataPoint> + <measurement value="-0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="0.249" errorPlus="0.047" errorMinus="0.047"/> + </dataPoint> + <dataPoint> + <measurement value="-0.35" errorPlus="0.04999999999999999" errorMinus="0.050000000000000044"/> + <measurement value="0.199" errorPlus="0.037" errorMinus="0.037"/> + </dataPoint> + <dataPoint> + <measurement value="-0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="0.18" errorPlus="0.04" errorMinus="0.04"/> + </dataPoint> + <dataPoint> + <measurement value="-0.15000000000000002" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/> + <measurement value="0.112" errorPlus="0.04" errorMinus="0.04"/> + </dataPoint> + <dataPoint> + <measurement value="-0.05" errorPlus="0.05" errorMinus="0.05"/> + <measurement value="0.137" errorPlus="0.04" errorMinus="0.04"/> + </dataPoint> + <dataPoint> + <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/> + <measurement value="0.162" errorPlus="0.037" errorMinus="0.037"/> + </dataPoint> + <dataPoint> + <measurement value="0.15000000000000002" errorPlus="0.04999999999999999" errorMinus="0.05000000000000002"/> + <measurement value="0.149" errorPlus="0.037" errorMinus="0.037"/> + </dataPoint> + <dataPoint> + <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="0.143" errorPlus="0.037" errorMinus="0.037"/> + </dataPoint> + <dataPoint> + <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/> + <measurement value="0.193" errorPlus="0.037" errorMinus="0.037"/> + </dataPoint> + <dataPoint> + <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="0.199" errorPlus="0.031" errorMinus="0.031"/> + </dataPoint> + <dataPoint> + <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="0.323" errorPlus="0.04" errorMinus="0.04"/> + </dataPoint> + <dataPoint> + <measurement value="0.6499999999999999" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/> + <measurement value="0.529" errorPlus="0.065" errorMinus="0.065"/> + </dataPoint> + <dataPoint> + <measurement value="0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="0.92" errorPlus="0.093" errorMinus="0.093"/> + </dataPoint> + <dataPoint> + <measurement value="0.8500000000000001" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="2.32" errorPlus="0.177" errorMinus="0.177"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DPSI(P=3,Q=THREE-JET)" > + <dataPoint> + <measurement value="0.157" errorPlus="0.157" errorMinus="0.157"/> + <measurement value="0.737" errorPlus="0.042" errorMinus="0.042"/> + </dataPoint> + <dataPoint> + <measurement value="0.471" errorPlus="0.15700000000000003" errorMinus="0.15699999999999997"/> + <measurement value="0.383" errorPlus="0.031" errorMinus="0.031"/> + </dataPoint> + <dataPoint> + <measurement value="0.7855" errorPlus="0.15749999999999997" errorMinus="0.15749999999999997"/> + <measurement value="0.225" errorPlus="0.023" errorMinus="0.023"/> + </dataPoint> + <dataPoint> + <measurement value="1.0999999999999999" errorPlus="0.15700000000000003" errorMinus="0.15699999999999992"/> + <measurement value="0.189" errorPlus="0.021" errorMinus="0.021"/> + </dataPoint> + <dataPoint> + <measurement value="1.414" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.122" errorPlus="0.015" errorMinus="0.015"/> + </dataPoint> + <dataPoint> + <measurement value="1.728" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.139" errorPlus="0.02" errorMinus="0.02"/> + </dataPoint> + <dataPoint> + <measurement value="2.042" errorPlus="0.15700000000000003" errorMinus="0.1569999999999998"/> + <measurement value="0.143" errorPlus="0.02" errorMinus="0.02"/> + </dataPoint> + <dataPoint> + <measurement value="2.3564999999999996" errorPlus="0.1575000000000002" errorMinus="0.15749999999999975"/> + <measurement value="0.236" errorPlus="0.022" errorMinus="0.022"/> + </dataPoint> + <dataPoint> + <measurement value="2.671" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.392" errorPlus="0.026" errorMinus="0.026"/> + </dataPoint> + <dataPoint> + <measurement value="2.985" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.629" errorPlus="0.04" errorMinus="0.04"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(P=3,Q=THREE-BODY)/DM(C=6JET))" > + <dataPoint> + <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/> + <measurement value="0.23" errorPlus="0.26" errorMinus="0.26"/> + </dataPoint> + <dataPoint> + <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/> + <measurement value="3.32" errorPlus="0.34" errorMinus="0.34"/> + </dataPoint> + <dataPoint> + <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/> + <measurement value="1.43" errorPlus="0.28" errorMinus="0.28"/> + </dataPoint> + <dataPoint> + <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/> + <measurement value="1.09" errorPlus="0.28" errorMinus="0.28"/> + </dataPoint> + <dataPoint> + <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/> + <measurement value="2.53" errorPlus="0.3" errorMinus="0.3"/> + </dataPoint> + <dataPoint> + <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/> + <measurement value="3.66" errorPlus="0.34" errorMinus="0.34"/> + </dataPoint> + <dataPoint> + <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/> + <measurement value="3.96" errorPlus="0.32" errorMinus="0.32"/> + </dataPoint> + <dataPoint> + <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/> + <measurement value="4.72" errorPlus="0.36" errorMinus="0.36"/> + </dataPoint> + <dataPoint> + <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/> + <measurement value="5.25" errorPlus="0.38" errorMinus="0.38"/> + </dataPoint> + <dataPoint> + <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/> + <measurement value="4.42" errorPlus="0.32" errorMinus="0.32"/> + </dataPoint> + <dataPoint> + <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/> + <measurement value="4.19" errorPlus="0.34" errorMinus="0.34"/> + </dataPoint> + <dataPoint> + <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/> + <measurement value="2.15" errorPlus="0.32" errorMinus="0.32"/> + </dataPoint> + <dataPoint> + <measurement value="0.3335" errorPlus="0.013499999999999956" errorMinus="0.013500000000000012"/> + <measurement value="0.75" errorPlus="0.26" errorMinus="0.26"/> + </dataPoint> + <dataPoint> + <measurement value="0.36" errorPlus="0.013000000000000012" errorMinus="0.013000000000000012"/> + <measurement value="0.19" errorPlus="0.3" errorMinus="0.3"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y02" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(P=4,Q=THREE-BODY)/DM(C=6JET))" > + <dataPoint> + <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/> + <measurement value="3.04" errorPlus="0.35" errorMinus="0.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/> + <measurement value="5.69" errorPlus="0.4" errorMinus="0.4"/> + </dataPoint> + <dataPoint> + <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/> + <measurement value="3.88" errorPlus="0.33" errorMinus="0.33"/> + </dataPoint> + <dataPoint> + <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/> + <measurement value="3.88" errorPlus="0.35" errorMinus="0.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/> + <measurement value="4.01" errorPlus="0.33" errorMinus="0.33"/> + </dataPoint> + <dataPoint> + <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/> + <measurement value="4.63" errorPlus="0.38" errorMinus="0.38"/> + </dataPoint> + <dataPoint> + <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/> + <measurement value="3.93" errorPlus="0.38" errorMinus="0.38"/> + </dataPoint> + <dataPoint> + <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/> + <measurement value="3.62" errorPlus="0.31" errorMinus="0.31"/> + </dataPoint> + <dataPoint> + <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/> + <measurement value="2.74" errorPlus="0.33" errorMinus="0.33"/> + </dataPoint> + <dataPoint> + <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/> + <measurement value="1.28" errorPlus="0.35" errorMinus="0.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/> + <measurement value="0.57" errorPlus="0.31" errorMinus="0.31"/> + </dataPoint> + <dataPoint> + <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/> + <measurement value="0.26" errorPlus="0.31" errorMinus="0.31"/> + </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="d06-x01-y03" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(P=5,Q=THREE-BODY)/DM(C=6JET))" > + <dataPoint> + <measurement value="0.0135" errorPlus="0.0135" errorMinus="0.0135"/> + <measurement value="3.22" errorPlus="0.29" errorMinus="0.29"/> + </dataPoint> + <dataPoint> + <measurement value="0.04" errorPlus="0.012999999999999998" errorMinus="0.013000000000000001"/> + <measurement value="7.01" errorPlus="0.44" errorMinus="0.44"/> + </dataPoint> + <dataPoint> + <measurement value="0.0665" errorPlus="0.013499999999999998" errorMinus="0.013500000000000005"/> + <measurement value="4.5" errorPlus="0.35" errorMinus="0.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.0935" errorPlus="0.013499999999999998" errorMinus="0.013499999999999998"/> + <measurement value="4.81" errorPlus="0.4" errorMinus="0.4"/> + </dataPoint> + <dataPoint> + <measurement value="0.12" errorPlus="0.013000000000000012" errorMinus="0.012999999999999998"/> + <measurement value="5.6" errorPlus="0.33" errorMinus="0.33"/> + </dataPoint> + <dataPoint> + <measurement value="0.14650000000000002" errorPlus="0.013499999999999984" errorMinus="0.013500000000000012"/> + <measurement value="4.41" errorPlus="0.38" errorMinus="0.38"/> + </dataPoint> + <dataPoint> + <measurement value="0.1735" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/> + <measurement value="3.22" errorPlus="0.33" errorMinus="0.33"/> + </dataPoint> + <dataPoint> + <measurement value="0.2" errorPlus="0.012999999999999984" errorMinus="0.013000000000000012"/> + <measurement value="2.25" errorPlus="0.35" errorMinus="0.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.22649999999999998" errorPlus="0.013500000000000012" errorMinus="0.013499999999999984"/> + <measurement value="1.76" errorPlus="0.33" errorMinus="0.33"/> + </dataPoint> + <dataPoint> + <measurement value="0.2535" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/> + <measurement value="0.62" errorPlus="0.37" errorMinus="0.37"/> + </dataPoint> + <dataPoint> + <measurement value="0.28" errorPlus="0.012999999999999956" errorMinus="0.013000000000000012"/> + <measurement value="0.18" errorPlus="0.38" errorMinus="0.38"/> + </dataPoint> + <dataPoint> + <measurement value="0.3065" errorPlus="0.013500000000000012" errorMinus="0.013500000000000012"/> + <measurement value="0.09" errorPlus="0.35" errorMinus="0.35"/> + </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="d07-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DX(C=A)" > + <dataPoint> + <measurement value="0.525" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="5.26" errorPlus="0.25" errorMinus="0.25"/> + </dataPoint> + <dataPoint> + <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="5.09" errorPlus="0.27" errorMinus="0.27"/> + </dataPoint> + <dataPoint> + <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.59" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="2.95" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.01" errorPlus="0.22" errorMinus="0.22"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.87" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.3" errorPlus="0.22" errorMinus="0.22"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.03" errorPlus="0.2" errorMinus="0.2"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d08-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DX(C=C)" > + <dataPoint> + <measurement value="0.525" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.57" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="3.41" errorPlus="0.24" errorMinus="0.24"/> + </dataPoint> + <dataPoint> + <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.84" errorPlus="0.25" errorMinus="0.25"/> + </dataPoint> + <dataPoint> + <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="3.15" errorPlus="0.21" errorMinus="0.21"/> + </dataPoint> + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.96" errorPlus="0.24" errorMinus="0.24"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.63" errorPlus="0.16" errorMinus="0.16"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.99" errorPlus="0.17" errorMinus="0.17"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.45" errorPlus="0.16" errorMinus="0.16"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d09-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DX(C=E)" > + <dataPoint> + <measurement value="0.525" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="4.51" errorPlus="0.27" errorMinus="0.27"/> + </dataPoint> + <dataPoint> + <measurement value="0.575" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="4.56" errorPlus="0.24" errorMinus="0.24"/> + </dataPoint> + <dataPoint> + <measurement value="0.625" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.54" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.675" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="2.98" errorPlus="0.23" errorMinus="0.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.47" errorPlus="0.19" errorMinus="0.19"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.4" errorPlus="0.16" errorMinus="0.16"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.59" errorPlus="0.19" errorMinus="0.19"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.13" errorPlus="0.17" errorMinus="0.17"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d10-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DPSI(C=AB)" > + <dataPoint> + <measurement value="0.157" errorPlus="0.157" errorMinus="0.157"/> + <measurement value="0.384" errorPlus="0.043" errorMinus="0.043"/> + </dataPoint> + <dataPoint> + <measurement value="0.471" errorPlus="0.15700000000000003" errorMinus="0.15699999999999997"/> + <measurement value="0.363" errorPlus="0.031" errorMinus="0.031"/> + </dataPoint> + <dataPoint> + <measurement value="0.7855" errorPlus="0.15749999999999997" errorMinus="0.15749999999999997"/> + <measurement value="0.323" errorPlus="0.031" errorMinus="0.031"/> + </dataPoint> + <dataPoint> + <measurement value="1.0999999999999999" errorPlus="0.15700000000000003" errorMinus="0.15699999999999992"/> + <measurement value="0.295" errorPlus="0.032" errorMinus="0.032"/> + </dataPoint> + <dataPoint> + <measurement value="1.414" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.262" errorPlus="0.022" errorMinus="0.022"/> + </dataPoint> + <dataPoint> + <measurement value="1.728" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.259" errorPlus="0.027" errorMinus="0.027"/> + </dataPoint> + <dataPoint> + <measurement value="2.042" errorPlus="0.15700000000000003" errorMinus="0.1569999999999998"/> + <measurement value="0.255" errorPlus="0.027" errorMinus="0.027"/> + </dataPoint> + <dataPoint> + <measurement value="2.3564999999999996" errorPlus="0.1575000000000002" errorMinus="0.15749999999999975"/> + <measurement value="0.266" errorPlus="0.025" errorMinus="0.025"/> + </dataPoint> + <dataPoint> + <measurement value="2.671" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.399" errorPlus="0.032" errorMinus="0.032"/> + </dataPoint> + <dataPoint> + <measurement value="2.985" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.395" errorPlus="0.031" errorMinus="0.031"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d11-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DPSI(C=CD)" > + <dataPoint> + <measurement value="0.157" errorPlus="0.157" errorMinus="0.157"/> + <measurement value="0.28" errorPlus="0.025" errorMinus="0.025"/> + </dataPoint> + <dataPoint> + <measurement value="0.471" errorPlus="0.15700000000000003" errorMinus="0.15699999999999997"/> + <measurement value="0.347" errorPlus="0.029" errorMinus="0.029"/> + </dataPoint> + <dataPoint> + <measurement value="0.7855" errorPlus="0.15749999999999997" errorMinus="0.15749999999999997"/> + <measurement value="0.349" errorPlus="0.031" errorMinus="0.031"/> + </dataPoint> + <dataPoint> + <measurement value="1.0999999999999999" errorPlus="0.15700000000000003" errorMinus="0.15699999999999992"/> + <measurement value="0.32" errorPlus="0.029" errorMinus="0.029"/> + </dataPoint> + <dataPoint> + <measurement value="1.414" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.309" errorPlus="0.028" errorMinus="0.028"/> + </dataPoint> + <dataPoint> + <measurement value="1.728" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.328" errorPlus="0.029" errorMinus="0.029"/> + </dataPoint> + <dataPoint> + <measurement value="2.042" errorPlus="0.15700000000000003" errorMinus="0.1569999999999998"/> + <measurement value="0.301" errorPlus="0.024" errorMinus="0.024"/> + </dataPoint> + <dataPoint> + <measurement value="2.3564999999999996" errorPlus="0.1575000000000002" errorMinus="0.15749999999999975"/> + <measurement value="0.293" errorPlus="0.027" errorMinus="0.027"/> + </dataPoint> + <dataPoint> + <measurement value="2.671" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.333" errorPlus="0.033" errorMinus="0.033"/> + </dataPoint> + <dataPoint> + <measurement value="2.985" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.312" errorPlus="0.04" errorMinus="0.04"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d12-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/DPSI(C=EF)" > + <dataPoint> + <measurement value="0.157" errorPlus="0.157" errorMinus="0.157"/> + <measurement value="0.257" errorPlus="0.023" errorMinus="0.023"/> + </dataPoint> + <dataPoint> + <measurement value="0.471" errorPlus="0.15700000000000003" errorMinus="0.15699999999999997"/> + <measurement value="0.235" errorPlus="0.023" errorMinus="0.023"/> + </dataPoint> + <dataPoint> + <measurement value="0.7855" errorPlus="0.15749999999999997" errorMinus="0.15749999999999997"/> + <measurement value="0.327" errorPlus="0.028" errorMinus="0.028"/> + </dataPoint> + <dataPoint> + <measurement value="1.0999999999999999" errorPlus="0.15700000000000003" errorMinus="0.15699999999999992"/> + <measurement value="0.4" errorPlus="0.035" errorMinus="0.035"/> + </dataPoint> + <dataPoint> + <measurement value="1.414" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.422" errorPlus="0.031" errorMinus="0.031"/> + </dataPoint> + <dataPoint> + <measurement value="1.728" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.435" errorPlus="0.036" errorMinus="0.036"/> + </dataPoint> + <dataPoint> + <measurement value="2.042" errorPlus="0.15700000000000003" errorMinus="0.1569999999999998"/> + <measurement value="0.378" errorPlus="0.03" errorMinus="0.03"/> + </dataPoint> + <dataPoint> + <measurement value="2.3564999999999996" errorPlus="0.1575000000000002" errorMinus="0.15749999999999975"/> + <measurement value="0.286" errorPlus="0.03" errorMinus="0.03"/> + </dataPoint> + <dataPoint> + <measurement value="2.671" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.249" errorPlus="0.028" errorMinus="0.028"/> + </dataPoint> + <dataPoint> + <measurement value="2.985" errorPlus="0.15700000000000003" errorMinus="0.15700000000000003"/> + <measurement value="0.219" errorPlus="0.023" errorMinus="0.023"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d13-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(C=A)/M(C=6JET))" > + <dataPoint> + <measurement value="0.0075" errorPlus="0.0075" errorMinus="0.0075"/> + <measurement value="2.88" errorPlus="1.26" errorMinus="1.26"/> + </dataPoint> + <dataPoint> + <measurement value="0.022" errorPlus="0.007000000000000003" errorMinus="0.006999999999999999"/> + <measurement value="13.69" errorPlus="1.35" errorMinus="1.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.0365" errorPlus="0.0075" errorMinus="0.007499999999999996"/> + <measurement value="17.12" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.051000000000000004" errorPlus="0.006999999999999999" errorMinus="0.007000000000000006"/> + <measurement value="9.91" errorPlus="1.17" errorMinus="1.17"/> + </dataPoint> + <dataPoint> + <measurement value="0.0655" errorPlus="0.007499999999999993" errorMinus="0.0075"/> + <measurement value="7.39" errorPlus="1.17" errorMinus="1.17"/> + </dataPoint> + <dataPoint> + <measurement value="0.07999999999999999" errorPlus="0.007000000000000006" errorMinus="0.006999999999999992"/> + <measurement value="6.13" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.0945" errorPlus="0.007499999999999993" errorMinus="0.007500000000000007"/> + <measurement value="4.68" errorPlus="1.17" errorMinus="1.17"/> + </dataPoint> + <dataPoint> + <measurement value="0.109" errorPlus="0.007000000000000006" errorMinus="0.007000000000000006"/> + <measurement value="3.24" errorPlus="1.35" errorMinus="1.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.1235" errorPlus="0.007500000000000007" errorMinus="0.007499999999999993"/> + <measurement value="1.98" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.138" errorPlus="0.0069999999999999785" errorMinus="0.007000000000000006"/> + <measurement value="1.26" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.1525" errorPlus="0.007500000000000007" errorMinus="0.007500000000000007"/> + <measurement value="0.72" errorPlus="0.9" errorMinus="0.9"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d14-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(C=B)/M(C=6JET))" > + <dataPoint> + <measurement value="0.0075" errorPlus="0.0075" errorMinus="0.0075"/> + <measurement value="2.88" errorPlus="1.26" errorMinus="1.26"/> + </dataPoint> + <dataPoint> + <measurement value="0.022" errorPlus="0.007000000000000003" errorMinus="0.006999999999999999"/> + <measurement value="13.69" errorPlus="1.35" errorMinus="1.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.0365" errorPlus="0.0075" errorMinus="0.007499999999999996"/> + <measurement value="17.12" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.051000000000000004" errorPlus="0.006999999999999999" errorMinus="0.007000000000000006"/> + <measurement value="9.91" errorPlus="1.17" errorMinus="1.17"/> + </dataPoint> + <dataPoint> + <measurement value="0.0655" errorPlus="0.007499999999999993" errorMinus="0.0075"/> + <measurement value="7.39" errorPlus="1.17" errorMinus="1.17"/> + </dataPoint> + <dataPoint> + <measurement value="0.07999999999999999" errorPlus="0.007000000000000006" errorMinus="0.006999999999999992"/> + <measurement value="6.13" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.0945" errorPlus="0.007499999999999993" errorMinus="0.007500000000000007"/> + <measurement value="4.68" errorPlus="1.17" errorMinus="1.17"/> + </dataPoint> + <dataPoint> + <measurement value="0.109" errorPlus="0.007000000000000006" errorMinus="0.007000000000000006"/> + <measurement value="3.24" errorPlus="1.35" errorMinus="1.35"/> + </dataPoint> + <dataPoint> + <measurement value="0.1235" errorPlus="0.007500000000000007" errorMinus="0.007499999999999993"/> + <measurement value="1.98" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.138" errorPlus="0.0069999999999999785" errorMinus="0.007000000000000006"/> + <measurement value="1.26" errorPlus="1.08" errorMinus="1.08"/> + </dataPoint> + <dataPoint> + <measurement value="0.1525" errorPlus="0.007500000000000007" errorMinus="0.007500000000000007"/> + <measurement value="0.72" errorPlus="0.9" errorMinus="0.9"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d15-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(C=C)/M(C=6JET))" > + <dataPoint> + <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/> + <measurement value="2.38" errorPlus="2.23" errorMinus="2.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/> + <measurement value="19.6" errorPlus="1.78" errorMinus="1.78"/> + </dataPoint> + <dataPoint> + <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/> + <measurement value="27.92" errorPlus="2.23" errorMinus="2.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/> + <measurement value="19.01" errorPlus="2.23" errorMinus="2.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/> + <measurement value="13.07" errorPlus="2.23" errorMinus="2.23"/> + </dataPoint> + <dataPoint> + <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/> + <measurement value="8.02" errorPlus="1.78" errorMinus="1.78"/> + </dataPoint> + <dataPoint> + <measurement value="0.065" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/> + <measurement value="4.46" errorPlus="1.78" errorMinus="1.78"/> + </dataPoint> + <dataPoint> + <measurement value="0.07500000000000001" errorPlus="0.0049999999999999906" errorMinus="0.0050000000000000044"/> + <measurement value="2.67" errorPlus="1.63" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="0.08499999999999999" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999906"/> + <measurement value="2.08" errorPlus="1.78" errorMinus="1.78"/> + </dataPoint> + <dataPoint> + <measurement value="0.095" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/> + <measurement value="0.59" errorPlus="1.93" errorMinus="1.93"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d16-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(C=D)/M(C=6JET))" > + <dataPoint> + <measurement value="0.0050" errorPlus="0.0050" errorMinus="0.0050"/> + <measurement value="5.09" errorPlus="1.79" errorMinus="1.79"/> + </dataPoint> + <dataPoint> + <measurement value="0.015" errorPlus="0.005000000000000001" errorMinus="0.004999999999999999"/> + <measurement value="40.04" errorPlus="2.09" errorMinus="2.09"/> + </dataPoint> + <dataPoint> + <measurement value="0.025" errorPlus="0.0049999999999999975" errorMinus="0.005000000000000001"/> + <measurement value="29.91" errorPlus="1.94" errorMinus="1.94"/> + </dataPoint> + <dataPoint> + <measurement value="0.035" errorPlus="0.0049999999999999975" errorMinus="0.0050000000000000044"/> + <measurement value="11.73" errorPlus="1.94" errorMinus="1.94"/> + </dataPoint> + <dataPoint> + <measurement value="0.045" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999975"/> + <measurement value="4.59" errorPlus="1.94" errorMinus="1.94"/> + </dataPoint> + <dataPoint> + <measurement value="0.055" errorPlus="0.0049999999999999975" errorMinus="0.0049999999999999975"/> + <measurement value="2.82" errorPlus="1.79" errorMinus="1.79"/> + </dataPoint> + <dataPoint> + <measurement value="0.065" errorPlus="0.0050000000000000044" errorMinus="0.0050000000000000044"/> + <measurement value="2.85" errorPlus="1.94" errorMinus="1.94"/> + </dataPoint> + <dataPoint> + <measurement value="0.07500000000000001" errorPlus="0.0049999999999999906" errorMinus="0.0050000000000000044"/> + <measurement value="2.28" errorPlus="1.94" errorMinus="1.94"/> + </dataPoint> + <dataPoint> + <measurement value="0.08499999999999999" errorPlus="0.0050000000000000044" errorMinus="0.0049999999999999906"/> + <measurement value="1.7" errorPlus="1.64" errorMinus="1.64"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d17-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(C=E)/M(C=6JET))" > + <dataPoint> + <measurement value="0.0025" errorPlus="0.0025" errorMinus="0.0025"/> + <measurement value="0.53" errorPlus="3.45" errorMinus="3.45"/> + </dataPoint> + <dataPoint> + <measurement value="0.0075" errorPlus="0.0025000000000000005" errorMinus="0.0024999999999999996"/> + <measurement value="14.87" errorPlus="3.19" errorMinus="3.19"/> + </dataPoint> + <dataPoint> + <measurement value="0.0125" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000005"/> + <measurement value="51.5" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.0175" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000022"/> + <measurement value="49.38" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.0225" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999988"/> + <measurement value="31.33" errorPlus="3.98" errorMinus="3.98"/> + </dataPoint> + <dataPoint> + <measurement value="0.0275" errorPlus="0.0024999999999999988" errorMinus="0.0024999999999999988"/> + <measurement value="20.18" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.0325" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/> + <measurement value="9.03" errorPlus="3.45" errorMinus="3.45"/> + </dataPoint> + <dataPoint> + <measurement value="0.037500000000000006" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/> + <measurement value="5.31" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.042499999999999996" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999953"/> + <measurement value="3.72" errorPlus="2.92" errorMinus="2.92"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d18-x01-y01" dimension="2" path="/REF/CDF_1997_S3541940" title="(1/N)*D(N)/D(M(C=F)/M(C=6JET))" > + <dataPoint> + <measurement value="0.0025" errorPlus="0.0025" errorMinus="0.0025"/> + <measurement value="0.53" errorPlus="2.92" errorMinus="2.92"/> + </dataPoint> + <dataPoint> + <measurement value="0.0075" errorPlus="0.0025000000000000005" errorMinus="0.0024999999999999996"/> + <measurement value="27.61" errorPlus="3.45" errorMinus="3.45"/> + </dataPoint> + <dataPoint> + <measurement value="0.0125" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000005"/> + <measurement value="75.93" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.0175" errorPlus="0.0024999999999999988" errorMinus="0.0025000000000000022"/> + <measurement value="52.04" errorPlus="3.45" errorMinus="3.45"/> + </dataPoint> + <dataPoint> + <measurement value="0.0225" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999988"/> + <measurement value="19.65" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.0275" errorPlus="0.0024999999999999988" errorMinus="0.0024999999999999988"/> + <measurement value="6.9" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.0325" errorPlus="0.0025000000000000022" errorMinus="0.0025000000000000022"/> + <measurement value="2.65" errorPlus="3.72" errorMinus="3.72"/> + </dataPoint> + <dataPoint> + <measurement value="0.037500000000000006" errorPlus="0.0024999999999999953" errorMinus="0.0025000000000000022"/> + <measurement value="2.12" errorPlus="4.78" errorMinus="4.78"/> + </dataPoint> + <dataPoint> + <measurement value="0.042499999999999996" errorPlus="0.0025000000000000022" errorMinus="0.0024999999999999953"/> + <measurement value="0.0" errorPlus="3.43" errorMinus="3.43"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/src/Analyses/CDF_1996_S3349578.cc ============================================================================== --- trunk/src/Analyses/CDF_1996_S3349578.cc Tue Mar 2 13:53:07 2010 (r2293) +++ trunk/src/Analyses/CDF_1996_S3349578.cc Tue Mar 2 16:06:02 2010 (r2294) @@ -19,7 +19,6 @@ : Analysis("CDF_1996_S3349578") { setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); } //@} @@ -80,56 +79,71 @@ } - /// 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) { + bool separated=true; + foreach (const Jet& ref, jets) { + if (deltaR(jet.momentum(), ref.momentum())<0.9) { + separated=false; + break; + } + } + if (!separated) continue; jets.push_back(jet); - sumEt += Et; jetsystem += jet.momentum(); } + if (jets.size()>=5) break; } /// @todo include gaussian jet energy resolution smearing? - if (jets.size() < 3) { - vetoEvent; - } - - if (sumEt < 420.0*GeV) { - vetoEvent; + if (jets.size() > 4) { + _fiveJetAnalysis(jets, weight); + jets.resize(4); + } + if (jets.size() > 3) { + _fourJetAnalysis(jets, weight); + jets.resize(3); } - 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(); + + double sumEt=0.0; + FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0); + foreach (const Jet& jet, jets) { + sumEt+=jet.momentum().Et(); + jetsystem+=jet.momentum(); + } + if (sumEt < 420.0*GeV) return; + + const double m3J = jetsystem.mass(); if (m3J<600*GeV) { return; } - LorentzTransform cms_boost(-jjj.boostVector()); - vector<FourMomentum> jets_boosted; + LorentzTransform cms_boost(-jetsystem.boostVector()); + vector<FourMomentum> jets3; foreach (Jet jet, jets) { - jets_boosted.push_back(cms_boost.transform(jet.momentum())); + jets3.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]); + std::sort(jets3.begin(), jets3.end(), FourMomentum::byEDescending()); + FourMomentum p3(jets3[0]); + FourMomentum p4(jets3[1]); + FourMomentum p5(jets3[2]); - double costheta3 = cos(p3.theta()); + FourMomentum pAV = cms_boost.transform(_avg_beam_in_lab(m3J, jetsystem.y())); + double costheta3=pAV.vector3().unit().dot(p3.vector3().unit()); if (fabs(costheta3)>0.6) { return; } @@ -138,14 +152,9 @@ 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 psi3 = _psi(p3, pAV, p4, p5); const double f3 = p3.mass()/m3J; const double f4 = p4.mass()/m3J; const double f5 = p5.mass()/m3J; @@ -154,41 +163,46 @@ _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_psi3->fill(psi3, 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(); + + double sumEt=0.0; + FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0); + foreach (const Jet& jet, jets) { + sumEt+=jet.momentum().Et(); + jetsystem+=jet.momentum(); + } + if (sumEt < 420.0*GeV) return; + + const double m4J = jetsystem.mass(); if (m4J < 650*GeV) return; + LorentzTransform cms_boost(-jetsystem.boostVector()); + vector<FourMomentum> jets4; + foreach (Jet jet, jets) { + jets4.push_back(cms_boost.transform(jet.momentum())); + } + std::sort(jets4.begin(), jets4.end(), FourMomentum::byEDescending()); + 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]); + vector<FourMomentum> jets3(_reduce(jets4, pA, pB)); + std::sort(jets3.begin(), jets3.end(), FourMomentum::byEDescending()); + FourMomentum p3(jets3[0]); + FourMomentum p4(jets3[1]); + FourMomentum p5(jets3[2]); - const double costheta3 = cos(p3.theta()); + FourMomentum pAV = cms_boost.transform(_avg_beam_in_lab(m4J, jetsystem.y())); + double costheta3=pAV.vector3().unit().dot(p3.vector3().unit()); if (fabs(costheta3)>0.8) { return; } @@ -200,108 +214,96 @@ // 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 psi3 = _psi(p3, pAV, p4, p5); 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(); + const double psiAB = _psi(pA, pB, pA+pB, pAV); _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_psi3->fill(psi3, 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_psiAB->fill(psiAB, 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(); + + double sumEt=0.0; + FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0); + foreach (const Jet& jet, jets) { + sumEt+=jet.momentum().Et(); + jetsystem+=jet.momentum(); + } + if (sumEt < 420.0*GeV) return; + + const double m5J = jetsystem.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]); + + LorentzTransform cms_boost(-jetsystem.boostVector()); + vector<FourMomentum> jets5; + foreach (Jet jet, jets) { + jets5.push_back(cms_boost.transform(jet.momentum())); + } + std::sort(jets5.begin(), jets5.end(), FourMomentum::byEDescending()); + + FourMomentum pC, pD; + vector<FourMomentum> jets4(_reduce(jets5, pC, pD)); + std::sort(jets4.begin(), jets4.end(), FourMomentum::byEDescending()); + + FourMomentum pA, pB; + vector<FourMomentum> jets3(_reduce(jets4, pA, pB)); + std::sort(jets3.begin(), jets3.end(), FourMomentum::byEDescending()); + FourMomentum p3(jets3[0]); + FourMomentum p4(jets3[1]); + FourMomentum p5(jets3[2]); // fill histograms - const double costheta3 = cos(p3.theta()); + FourMomentum pAV = cms_boost.transform(_avg_beam_in_lab(m5J, jetsystem.y())); + const double costheta3 = pAV.vector3().unit().dot(p3.vector3().unit()); 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 psi3 = _psi(p3, pAV, p4, p5); 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 psiAB = _psi(pA, pB, pA+pB, pAV); 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(); + const double psiCD = _psi(pC, pD, pC+pD, pAV); _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_psi3->fill(psi3, 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_psiAB->fill(psiAB, 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_psiCD->fill(psiCD, weight); _h_5_fC->fill(fC, weight); _h_5_fD->fill(fD, weight); } @@ -360,13 +362,11 @@ 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; } @@ -376,11 +376,35 @@ for (size_t i=0; i<jets.size(); ++i) { if (i!=idx1 && i!=idx2) newjets.push_back(jets[i]); } - newjets.push_back(combined); + newjets.push_back(jets[idx1]+jets[idx2]); combined1 = jets[idx1]; combined2 = jets[idx2]; return newjets; } + + FourMomentum _avg_beam_in_lab(const double& m, const double& y) { + FourMomentum boostvec(cosh(y), 0.0, 0.0, sinh(y)); + LorentzTransform cms_boost(-boostvec.boostVector()); + cms_boost = cms_boost.inverse(); + const double mt = m/2.0; + FourMomentum beam1(mt, 0, 0, mt); + FourMomentum beam2(mt, 0, 0, -mt); + beam1=cms_boost.transform(beam1); + beam2=cms_boost.transform(beam2); + if (beam1.E()>beam2.E()) { + return beam1-beam2; + } + else { + return beam2-beam1; + } + } + + double _psi(const FourMomentum& p1, const FourMomentum& p2, + const FourMomentum& p3, const FourMomentum& p4) { + Vector3 p1xp2 = p1.vector3().cross(p2.vector3()); + Vector3 p3xp4 = p3.vector3().cross(p4.vector3()); + return mapAngle0ToPi(acos(p1xp2.unit().dot(p3xp4.unit()))); + } private: Added: trunk/src/Analyses/CDF_1997_S3541940.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_1997_S3541940.cc Tue Mar 2 16:06:02 2010 (r2294) @@ -0,0 +1,265 @@ +// -*- 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_1997_S3541940 : public Analysis { + public: + + CDF_1997_S3541940() + : Analysis("CDF_1997_S3541940") + { + setBeams(PROTON, ANTIPROTON); + } + + + public: + + void init() { + + const FinalState fs(-4.2, 4.2); + addProjection(FastJets(fs, FastJets::CDFJETCLU, 0.7), "Jets"); + + _h_m6J = bookHistogram1D(1, 1, 1); + _h_X3ppp = bookHistogram1D(2, 1, 1); + _h_X4ppp = bookHistogram1D(3, 1, 1); + _h_costheta3ppp = bookHistogram1D(4, 1, 1); + _h_psi3ppp = bookHistogram1D(5, 1, 1); + _h_f3ppp = bookHistogram1D(6, 1, 1); + _h_f4ppp = bookHistogram1D(6, 1, 2); + _h_f5ppp = bookHistogram1D(6, 1, 3); + _h_XApp = bookHistogram1D(7, 1, 1); + _h_XCp = bookHistogram1D(8, 1, 1); + _h_XE = bookHistogram1D(9, 1, 1); + _h_psiAppBpp = bookHistogram1D(10, 1, 1); + _h_psiCpDp = bookHistogram1D(11, 1, 1); + _h_psiEF = bookHistogram1D(12, 1, 1); + _h_fApp = bookHistogram1D(13, 1, 1); + _h_fBpp = bookHistogram1D(14, 1, 1); + _h_fCp = bookHistogram1D(15, 1, 1); + _h_fDp = bookHistogram1D(16, 1, 1); + _h_fE = bookHistogram1D(17, 1, 1); + _h_fF = bookHistogram1D(18, 1, 1); + } + + + void analyze(const Event& event) { + const double weight = event.weight(); + + 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(); + double eta = fabs(jet.momentum().eta()); + if (Et > 20.0*GeV && eta<3.0) { + bool separated=true; + foreach (const Jet& ref, jets) { + if (deltaR(jet.momentum(), ref.momentum())<0.9) { + separated=false; + break; + } + } + if (!separated) continue; + jets.push_back(jet); + sumEt += Et; + jetsystem += jet.momentum(); + } + if (jets.size()>=6) break; + } + + if (jets.size()<6) { + vetoEvent; + } + + if (sumEt<320.0*GeV) { + vetoEvent; + } + + double m6J=jetsystem.mass(); + if (m6J<520.0*GeV) { + vetoEvent; + } + + LorentzTransform cms_boost(-jetsystem.boostVector()); + vector<FourMomentum> jets6; + foreach (Jet jet, jets) { + jets6.push_back(cms_boost.transform(jet.momentum())); + } + std::sort(jets6.begin(), jets6.end(), FourMomentum::byEDescending()); + + FourMomentum pE, pF; + vector<FourMomentum> jets5(_reduce(jets6, pE, pF)); + std::sort(jets5.begin(), jets5.end(), FourMomentum::byEDescending()); + + FourMomentum pCp, pDp; + vector<FourMomentum> jets4(_reduce(jets5, pCp, pDp)); + std::sort(jets4.begin(), jets4.end(), FourMomentum::byEDescending()); + + FourMomentum pApp, pBpp; + vector<FourMomentum> jets3(_reduce(jets4, pApp, pBpp)); + std::sort(jets3.begin(), jets3.end(), FourMomentum::byEDescending()); + FourMomentum p3ppp(jets3[0]); + FourMomentum p4ppp(jets3[1]); + FourMomentum p5ppp(jets3[2]); + + double X3ppp = 2.0*p3ppp.E()/m6J; + if (X3ppp>0.9) { + vetoEvent; + } + + FourMomentum pAV = cms_boost.transform(_avg_beam_in_lab(m6J, jetsystem.y())); + double costheta3ppp=pAV.vector3().unit().dot(p3ppp.vector3().unit()); + if (fabs(costheta3ppp)>0.9) { + vetoEvent; + } + + // 3-jet-system variables + _h_m6J->fill(m6J, weight); + _h_X3ppp->fill(X3ppp, weight); + _h_X4ppp->fill(2.0*p4ppp.E()/m6J, weight); + _h_costheta3ppp->fill(costheta3ppp, weight); + double psi3ppp=_psi(p3ppp, pAV, p4ppp, p5ppp); + _h_psi3ppp->fill(psi3ppp, weight); + _h_f3ppp->fill(p3ppp.mass()/m6J, weight); + _h_f4ppp->fill(p4ppp.mass()/m6J, weight); + _h_f5ppp->fill(p5ppp.mass()/m6J, weight); + + // 4 -> 3 jet variables + _h_fApp->fill(pApp.mass()/m6J, weight); + _h_fBpp->fill(pApp.mass()/m6J, weight); + _h_XApp->fill(pApp.E()/(pApp.E()+pBpp.E()), weight); + double psiAppBpp=_psi(pApp, pBpp, pApp+pBpp, pAV); + _h_psiAppBpp->fill(psiAppBpp, weight); + + // 5 -> 4 jet variables + _h_fCp->fill(pCp.mass()/m6J, weight); + _h_fDp->fill(pDp.mass()/m6J, weight); + _h_XCp->fill(pCp.E()/(pCp.E()+pDp.E()), weight); + double psiCpDp=_psi(pCp, pDp, pCp+pDp, pAV); + _h_psiCpDp->fill(psiCpDp, weight); + + // 6 -> 5 jet variables + _h_fE->fill(pE.mass()/m6J, weight); + _h_fF->fill(pF.mass()/m6J, weight); + _h_XE->fill(pE.E()/(pE.E()+pF.E()), weight); + double psiEF=_psi(pE, pF, pE+pF, pAV); + _h_psiEF->fill(psiEF, weight); + } + + + void finalize() { + + normalize(_h_m6J); + normalize(_h_X3ppp); + normalize(_h_X4ppp); + normalize(_h_costheta3ppp); + normalize(_h_psi3ppp); + normalize(_h_f3ppp); + normalize(_h_f4ppp); + normalize(_h_f5ppp); + normalize(_h_XApp); + normalize(_h_XCp); + normalize(_h_XE); + normalize(_h_psiAppBpp); + normalize(_h_psiCpDp); + normalize(_h_psiEF); + normalize(_h_fApp); + normalize(_h_fBpp); + normalize(_h_fCp); + normalize(_h_fDp); + normalize(_h_fE); + normalize(_h_fF); + + } + + + + private: + + vector<FourMomentum> _reduce(const vector<FourMomentum>& jets, + FourMomentum& combined1, + FourMomentum& combined2) { + double minMass2 = 1e9; + 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) { + 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(jets[idx1]+jets[idx2]); + combined1 = jets[idx1]; + combined2 = jets[idx2]; + return newjets; + } + + FourMomentum _avg_beam_in_lab(const double& m, const double& y) { + FourMomentum boostvec(cosh(y), 0.0, 0.0, sinh(y)); + LorentzTransform cms_boost(-boostvec.boostVector()); + cms_boost = cms_boost.inverse(); + const double mt = m/2.0; + FourMomentum beam1(mt, 0, 0, mt); + FourMomentum beam2(mt, 0, 0, -mt); + beam1=cms_boost.transform(beam1); + beam2=cms_boost.transform(beam2); + if (beam1.E()>beam2.E()) { + return beam1-beam2; + } + else { + return beam2-beam1; + } + } + + double _psi(const FourMomentum& p1, const FourMomentum& p2, + const FourMomentum& p3, const FourMomentum& p4) { + Vector3 p1xp2 = p1.vector3().cross(p2.vector3()); + Vector3 p3xp4 = p3.vector3().cross(p4.vector3()); + return mapAngle0ToPi(acos(p1xp2.unit().dot(p3xp4.unit()))); + } + + + private: + + AIDA::IHistogram1D *_h_m6J; + AIDA::IHistogram1D *_h_X3ppp; + AIDA::IHistogram1D *_h_X4ppp; + AIDA::IHistogram1D *_h_costheta3ppp; + AIDA::IHistogram1D *_h_psi3ppp; + AIDA::IHistogram1D *_h_f3ppp; + AIDA::IHistogram1D *_h_f4ppp; + AIDA::IHistogram1D *_h_f5ppp; + AIDA::IHistogram1D *_h_XApp; + AIDA::IHistogram1D *_h_XCp; + AIDA::IHistogram1D *_h_XE; + AIDA::IHistogram1D *_h_psiAppBpp; + AIDA::IHistogram1D *_h_psiCpDp; + AIDA::IHistogram1D *_h_psiEF; + AIDA::IHistogram1D *_h_fApp; + AIDA::IHistogram1D *_h_fBpp; + AIDA::IHistogram1D *_h_fCp; + AIDA::IHistogram1D *_h_fDp; + AIDA::IHistogram1D *_h_fE; + AIDA::IHistogram1D *_h_fF; + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<CDF_1997_S3541940> plugin_CDF_1997_S3541940; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Tue Mar 2 13:53:07 2010 (r2293) +++ trunk/src/Analyses/Makefile.am Tue Mar 2 16:06:02 2010 (r2294) @@ -16,6 +16,7 @@ CDF_1996_S3349578.cc \ CDF_1996_S3418421.cc \ CDF_1994_S2952106.cc \ + CDF_1997_S3541940.cc \ CDF_2006_S6653332.cc \ CDF_2007_S7057202.cc \ CDF_2008_S7541902.cc \
More information about the Rivet-svn mailing list |