|
[Rivet-svn] r1836 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Sep 24 16:04:52 BST 2009
Author: fsiegert Date: Thu Sep 24 16:04:52 2009 New Revision: 1836 Log: Add ALEPH_1996_S3196992 (photon fragmentation function measurement). Added: trunk/data/anainfo/ALEPH_1996_S3196992.info trunk/data/plotinfo/ALEPH_1996_S3196992.plot trunk/data/refdata/ALEPH_1996_S3196992.aida trunk/src/Analyses/ALEPH_1996_S3196992.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/ALEPH_1996_S3196992.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/ALEPH_1996_S3196992.info Thu Sep 24 16:04:52 2009 (r1836) @@ -0,0 +1,26 @@ +Name: ALEPH_1996_S3196992 +Year: 1996 +Summary: Measurement of the quark to photon fragmentation function +Experiment: ALEPH +Collider: LEP Run 1 +SpiresID: 3196992 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - Z.Phys.C69:365-378,1996 + - doi:10.1007/s002880050037 +RunInfo: + e+ e- -> jets at shower level. The experimental results have been corrected + for hadronisation +NumEvents: 1000000 +PtCuts: [0] +Description: + Earlier measurements at LEP of isolated hard photons in hadronic Z decays, + attributed to radiation from primary quark pairs, have been extended in the + ALEPH experiment to include hard photon production inside hadron jets. Events + are selected where all particles combine democratically to form hadron jets, + one of which contains a photon with a fractional energy z > 0.7. After + statistical subtraction of non-prompt photons, the quark-to-photon + fragmentation function, D(z), is extracted directly from the measured 2-jet + rate. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Sep 24 12:51:39 2009 (r1835) +++ trunk/data/anainfo/Makefile.am Thu Sep 24 16:04:52 2009 (r1836) @@ -1,6 +1,7 @@ dist_pkgdata_DATA = \ ALEPH_1991_S2435284.info \ ALEPH_1996_S3486095.info \ + ALEPH_1996_S3196992.info \ CDF_1988_S1865951.info \ CDF_1990_S2089246.info \ CDF_1994_S2952106.info \ Added: trunk/data/plotinfo/ALEPH_1996_S3196992.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/ALEPH_1996_S3196992.plot Thu Sep 24 16:04:52 2009 (r1836) @@ -0,0 +1,47 @@ +# BEGIN PLOT /ALEPH_1996_S3196992/d01-x01-y01 +Title=Photon Fragmentation in 2-jet events with $y_\mathrm{cut}=0.01$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT + +# BEGIN PLOT /ALEPH_1996_S3196992/d02-x01-y01 +Title=Photon Fragmentation in 2-jet events with $y_\mathrm{cut}=0.06$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT + +# BEGIN PLOT /ALEPH_1996_S3196992/d03-x01-y01 +Title=Photon Fragmentation in 2-jet events with $y_\mathrm{cut}=0.1$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT + +# BEGIN PLOT /ALEPH_1996_S3196992/d04-x01-y01 +Title=Photon Fragmentation in 2-jet events with $y_\mathrm{cut}=0.33$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT + +# BEGIN PLOT /ALEPH_1996_S3196992/d05-x01-y01 +Title=Photon Fragmentation in 3-jet events with $y_\mathrm{cut}=0.01$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT + +# BEGIN PLOT /ALEPH_1996_S3196992/d06-x01-y01 +Title=Photon Fragmentation in 3-jet events with $y_\mathrm{cut}=0.06$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT + +# BEGIN PLOT /ALEPH_1996_S3196992/d07-x01-y01 +Title=Photon Fragmentation in 3-jet events with $y_\mathrm{cut}=0.1$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT + +# BEGIN PLOT /ALEPH_1996_S3196992/d08-x01-y01 +Title=Photon Fragmentation in 4-jet events with $y_\mathrm{cut}=0.01$ +XLabel=$z_\gamma$ +YLabel=$1/\sigma \mathrm{d}\sigma(\mathrm{2-jet})/\mathrm{d}z \times 10^3$ +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Sep 24 12:51:39 2009 (r1835) +++ trunk/data/plotinfo/Makefile.am Thu Sep 24 16:04:52 2009 (r1836) @@ -1,6 +1,7 @@ dist_pkgdata_DATA = \ ALEPH_1991_S2435284.plot \ ALEPH_1996_S3486095.plot \ + ALEPH_1996_S3196992.plot \ CDF_1988_S1865951.plot \ CDF_1990_S2089246.plot \ CDF_1994_S2952106.plot \ Added: trunk/data/refdata/ALEPH_1996_S3196992.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ALEPH_1996_S3196992.aida Thu Sep 24 16:04:52 2009 (r1836) @@ -0,0 +1,485 @@ +<?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/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.01" errorPlus="1.0280077820717117" errorMinus="1.0280077820717117"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.77" errorPlus="0.7060453243241541" errorMinus="0.7060453243241541"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.96" errorPlus="0.6657326790837296" errorMinus="0.6657326790837296"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.85" errorPlus="0.5660388679233962" errorMinus="0.5660388679233962"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.64" errorPlus="0.32449961479175904" errorMinus="0.32449961479175904"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="0.27" errorPlus="0.223606797749979" errorMinus="0.223606797749979"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="5.84" errorPlus="1.6297545827516486" errorMinus="1.6297545827516486"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.82" errorPlus="1.0857716150277643" errorMinus="1.0857716150277643"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.63" errorPlus="0.8373768566183329" errorMinus="0.8373768566183329"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.43" errorPlus="0.6603029607687672" errorMinus="0.6603029607687672"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.83" errorPlus="0.5103920062069938" errorMinus="0.5103920062069938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.53" errorPlus="0.3124099870362662" errorMinus="0.3124099870362662"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="0.92" errorPlus="0.27730849247724093" errorMinus="0.27730849247724093"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="12.74" errorPlus="2.409336008115099" errorMinus="2.409336008115099"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.95" errorPlus="1.0499523798725352" errorMinus="1.0499523798725352"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.68" errorPlus="0.8091971329657563" errorMinus="0.8091971329657563"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.25" errorPlus="0.6580273550544841" errorMinus="0.6580273550544841"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.81" errorPlus="0.5103920062069938" errorMinus="0.5103920062069938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.49" errorPlus="0.29832867780352595" errorMinus="0.29832867780352595"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="1.1" errorPlus="0.30463092423455634" errorMinus="0.30463092423455634"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="16.06" errorPlus="2.606338427756457" errorMinus="2.606338427756457"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.09" errorPlus="1.0499523798725352" errorMinus="1.0499523798725352"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.62" errorPlus="0.8015609770940699" errorMinus="0.8015609770940699"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.31" errorPlus="0.6640030120413613" errorMinus="0.6640030120413613"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.04" errorPlus="0.5103920062069938" errorMinus="0.5103920062069938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.53" errorPlus="0.29832867780352595" errorMinus="0.29832867780352595"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="1.77" errorPlus="0.340147027033899" errorMinus="0.340147027033899"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="23.68" errorPlus="3.0492294108512072" errorMinus="3.0492294108512072"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.32" errorPlus="0.8989994438263019" errorMinus="0.8989994438263019"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.08" errorPlus="0.73824115301167" errorMinus="0.73824115301167"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="2.15" errorPlus="0.7300684899377592" errorMinus="0.7300684899377592"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.49" errorPlus="0.6425729530566938" errorMinus="0.6425729530566938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="4.72" errorPlus="0.567978872846517" errorMinus="0.567978872846517"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="15.57" errorPlus="1.3285330255586423" errorMinus="1.3285330255586423"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="97.32" errorPlus="7.408994533673243" errorMinus="7.408994533673243"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.62" errorPlus="0.3080584360149873" errorMinus="0.3080584360149873"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.01" errorPlus="0.30083217912982646" errorMinus="0.30083217912982646"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.08" errorPlus="0.29120439557122074" errorMinus="0.29120439557122074"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.7" errorPlus="0.24515301344262524" errorMinus="0.24515301344262524"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="1.16" errorPlus="0.25" errorMinus="0.25"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="6.1" errorPlus="0.7038465741907111" errorMinus="0.7038465741907111"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="47.87" errorPlus="4.8582404222104945" errorMinus="4.8582404222104945"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d07-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.25" errorPlus="0.20223748416156687" errorMinus="0.20223748416156687"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.0" errorPlus="0.044721359549995794" errorMinus="0.044721359549995794"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.09" errorPlus="0.194164878389476" errorMinus="0.194164878389476"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.5" errorPlus="0.22360679774997896" errorMinus="0.22360679774997896"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.56" errorPlus="0.18439088914585774" errorMinus="0.18439088914585774"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="3.28" errorPlus="0.5536244214266564" errorMinus="0.5536244214266564"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="24.48" errorPlus="3.311223338888514" errorMinus="3.311223338888514"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d08-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="10E3*1/SIG(Q=HAD)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.84" errorPlus="0.5197114584074514" errorMinus="0.5197114584074514"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.36" errorPlus="0.35510561809129404" errorMinus="0.35510561809129404"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.91" errorPlus="0.5557877292636102" errorMinus="0.5557877292636102"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.78" errorPlus="0.31400636936215165" errorMinus="0.31400636936215165"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="1.73" errorPlus="0.383275357934736" errorMinus="0.383275357934736"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="3.5" errorPlus="0.6676076692189807" errorMinus="0.6676076692189807"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="21.54" errorPlus="3.3261539351028238" errorMinus="3.3261539351028238"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d09-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.01" errorPlus="1.0280077820717117" errorMinus="1.0280077820717117"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.77" errorPlus="0.7060453243241541" errorMinus="0.7060453243241541"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.96" errorPlus="0.6657326790837296" errorMinus="0.6657326790837296"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.85" errorPlus="0.5660388679233962" errorMinus="0.5660388679233962"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.64" errorPlus="0.32449961479175904" errorMinus="0.32449961479175904"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="0.27" errorPlus="0.223606797749979" errorMinus="0.223606797749979"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="5.84" errorPlus="1.6297545827516486" errorMinus="1.6297545827516486"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d10-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.82" errorPlus="1.0857716150277643" errorMinus="1.0857716150277643"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.63" errorPlus="0.8373768566183329" errorMinus="0.8373768566183329"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.43" errorPlus="0.6603029607687672" errorMinus="0.6603029607687672"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.83" errorPlus="0.5103920062069938" errorMinus="0.5103920062069938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.53" errorPlus="0.3124099870362662" errorMinus="0.3124099870362662"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="0.92" errorPlus="0.27730849247724093" errorMinus="0.27730849247724093"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="12.74" errorPlus="2.409336008115099" errorMinus="2.409336008115099"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d11-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.95" errorPlus="1.0499523798725352" errorMinus="1.0499523798725352"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.68" errorPlus="0.8091971329657563" errorMinus="0.8091971329657563"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.25" errorPlus="0.6580273550544841" errorMinus="0.6580273550544841"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.81" errorPlus="0.5103920062069938" errorMinus="0.5103920062069938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.49" errorPlus="0.29832867780352595" errorMinus="0.29832867780352595"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="1.1" errorPlus="0.30463092423455634" errorMinus="0.30463092423455634"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="16.06" errorPlus="2.606338427756457" errorMinus="2.606338427756457"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d12-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.09" errorPlus="1.0499523798725352" errorMinus="1.0499523798725352"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.62" errorPlus="0.8015609770940699" errorMinus="0.8015609770940699"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="1.31" errorPlus="0.6640030120413613" errorMinus="0.6640030120413613"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="1.04" errorPlus="0.5103920062069938" errorMinus="0.5103920062069938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.53" errorPlus="0.29832867780352595" errorMinus="0.29832867780352595"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="1.77" errorPlus="0.340147027033899" errorMinus="0.340147027033899"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="23.68" errorPlus="3.0492294108512072" errorMinus="3.0492294108512072"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d13-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.32" errorPlus="0.8989994438263019" errorMinus="0.8989994438263019"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="2.08" errorPlus="0.73824115301167" errorMinus="0.73824115301167"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="2.15" errorPlus="0.7300684899377592" errorMinus="0.7300684899377592"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="3.49" errorPlus="0.6425729530566938" errorMinus="0.6425729530566938"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="4.72" errorPlus="0.567978872846517" errorMinus="0.567978872846517"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="15.57" errorPlus="1.3285330255586423" errorMinus="1.3285330255586423"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="97.32" errorPlus="7.408994533673243" errorMinus="7.408994533673243"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d14-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.62" errorPlus="0.3080584360149873" errorMinus="0.3080584360149873"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.01" errorPlus="0.30083217912982646" errorMinus="0.30083217912982646"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.08" errorPlus="0.29120439557122074" errorMinus="0.29120439557122074"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.7" errorPlus="0.24515301344262524" errorMinus="0.24515301344262524"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="1.16" errorPlus="0.25" errorMinus="0.25"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="6.1" errorPlus="0.7038465741907111" errorMinus="0.7038465741907111"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="47.87" errorPlus="4.8582404222104945" errorMinus="4.8582404222104945"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d15-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.25" errorPlus="0.20223748416156687" errorMinus="0.20223748416156687"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.0" errorPlus="0.044721359549995794" errorMinus="0.044721359549995794"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.09" errorPlus="0.194164878389476" errorMinus="0.194164878389476"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.5" errorPlus="0.22360679774997896" errorMinus="0.22360679774997896"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="0.56" errorPlus="0.18439088914585774" errorMinus="0.18439088914585774"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="3.28" errorPlus="0.5536244214266564" errorMinus="0.5536244214266564"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="24.48" errorPlus="3.311223338888514" errorMinus="3.311223338888514"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d16-x01-y01" dimension="2" path="/REF/ALEPH_1996_S3196992" title="(1000/SIG)*D(SIG)/D(Z)" > + <dataPoint> + <measurement value="0.725" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.84" errorPlus="0.5197114584074514" errorMinus="0.5197114584074514"/> + </dataPoint> + <dataPoint> + <measurement value="0.775" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.36" errorPlus="0.35510561809129404" errorMinus="0.35510561809129404"/> + </dataPoint> + <dataPoint> + <measurement value="0.825" errorPlus="0.025000000000000022" errorMinus="0.02499999999999991"/> + <measurement value="0.91" errorPlus="0.5557877292636102" errorMinus="0.5557877292636102"/> + </dataPoint> + <dataPoint> + <measurement value="0.875" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.78" errorPlus="0.31400636936215165" errorMinus="0.31400636936215165"/> + </dataPoint> + <dataPoint> + <measurement value="0.925" errorPlus="0.02499999999999991" errorMinus="0.025000000000000022"/> + <measurement value="1.73" errorPlus="0.383275357934736" errorMinus="0.383275357934736"/> + </dataPoint> + <dataPoint> + <measurement value="0.97" errorPlus="0.020000000000000018" errorMinus="0.020000000000000018"/> + <measurement value="3.5" errorPlus="0.6676076692189807" errorMinus="0.6676076692189807"/> + </dataPoint> + <dataPoint> + <measurement value="0.995" errorPlus="0.0051" errorMinus="0.0050000000000000044"/> + <measurement value="21.54" errorPlus="3.3261539351028238" errorMinus="3.3261539351028238"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Sep 24 12:51:39 2009 (r1835) +++ trunk/data/refdata/Makefile.am Thu Sep 24 16:04:52 2009 (r1836) @@ -2,6 +2,7 @@ JADE_OPAL_2000_S4300807.aida \ ALEPH_1991_S2435284.aida \ ALEPH_1996_S3486095.aida \ + ALEPH_1996_S3196992.aida \ DELPHI_1994_S3021912.aida \ DELPHI_1995_S3137023.aida \ DELPHI_1996_S3430090.aida \ Added: trunk/src/Analyses/ALEPH_1996_S3196992.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/ALEPH_1996_S3196992.cc Thu Sep 24 16:04:52 2009 (r1836) @@ -0,0 +1,160 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" +#include "Rivet/Projections/FastJets.hh" + +namespace Rivet { + + + class ALEPH_1996_S3196992 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + ALEPH_1996_S3196992() : Analysis("ALEPH_1996_S3196992") + { + setBeams(ELECTRON, POSITRON); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms before the run + void init() { + + _h_z_2jet_001 = bookHistogram1D(1, 1, 1); + _h_z_2jet_006 = bookHistogram1D(2, 1, 1); + _h_z_2jet_01 = bookHistogram1D(3, 1, 1); + _h_z_2jet_033 = bookHistogram1D(4, 1, 1); + + _h_z_3jet_001 = bookHistogram1D(5, 1, 1); + _h_z_3jet_006 = bookHistogram1D(6, 1, 1); + _h_z_3jet_01 = bookHistogram1D(7, 1, 1); + + _h_z_4jet_001 = bookHistogram1D(8, 1, 1); + + + FinalState fs; + addProjection(FastJets(fs, FastJets::DURHAM, 0.7), "DurhamJets"); + + IdentifiedFinalState ifs(-MAXRAPIDITY, +MAXRAPIDITY, 5.0*GeV); + ifs.acceptId(PHOTON); + addProjection(ifs, "Photons"); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + const ParticleVector photons = applyProjection<IdentifiedFinalState>(event, "Photons").particles(); + if (photons.size()<1) { + vetoEvent; + } + + const FastJets& durjet = applyProjection<FastJets>(event, "DurhamJets"); + + foreach (const Particle& photon, photons) { + + PseudoJets jets_001 = durjet.clusterSeq()->exclusive_jets_ycut(0.01); + foreach (const fastjet::PseudoJet& jet, jets_001) { + if (particleInJet(photon, jet, durjet.clusterSeq())) { + double zgamma = photon.momentum().E()/jet.E(); + if (jets_001.size()==2) _h_z_2jet_001->fill(zgamma, weight); + else if (jets_001.size()==3) _h_z_3jet_001->fill(zgamma, weight); + else if (jets_001.size()==4) _h_z_4jet_001->fill(zgamma, weight); + break; + } + } + + PseudoJets jets_006 = durjet.clusterSeq()->exclusive_jets_ycut(0.06); + foreach (const fastjet::PseudoJet& jet, jets_006) { + if (particleInJet(photon, jet, durjet.clusterSeq())) { + double zgamma = photon.momentum().E()/jet.E(); + if (jets_006.size()==2) _h_z_2jet_006->fill(zgamma, weight); + else if (jets_006.size()==3) _h_z_3jet_006->fill(zgamma, weight); + break; + } + } + + PseudoJets jets_01 = durjet.clusterSeq()->exclusive_jets_ycut(0.1); + foreach (const fastjet::PseudoJet& jet, jets_01) { + if (particleInJet(photon, jet, durjet.clusterSeq())) { + double zgamma = photon.momentum().E()/jet.E(); + if (jets_01.size()==2) _h_z_2jet_01->fill(zgamma, weight); + else if (jets_01.size()==3) _h_z_3jet_01->fill(zgamma, weight); + break; + } + } + + PseudoJets jets_033 = durjet.clusterSeq()->exclusive_jets_ycut(0.33); + foreach (const fastjet::PseudoJet& jet, jets_033) { + if (particleInJet(photon, jet, durjet.clusterSeq())) { + double zgamma = photon.momentum().E()/jet.E(); + if (jets_033.size()==2) _h_z_2jet_033->fill(zgamma, weight); + break; + } + } + + } + } + + + bool particleInJet(const Particle& p, const fastjet::PseudoJet& jet, + const fastjet::ClusterSequence* cseq ) { + foreach (const fastjet::PseudoJet& jetpart, cseq->constituents(jet)) { + if (fuzzyEquals(jetpart.E(), p.momentum().E()) && + fuzzyEquals(jetpart.px(), p.momentum().x()) && + fuzzyEquals(jetpart.py(), p.momentum().y()) && + fuzzyEquals(jetpart.pz(), p.momentum().z())) { + return true; + } + } + return false; + } + + + /// Normalise histograms etc., after the run + void finalize() { + scale(_h_z_2jet_001, 1000.0/sumOfWeights()); + scale(_h_z_2jet_006, 1000.0/sumOfWeights()); + scale(_h_z_2jet_01, 1000.0/sumOfWeights()); + scale(_h_z_2jet_033, 1000.0/sumOfWeights()); + scale(_h_z_3jet_001, 1000.0/sumOfWeights()); + scale(_h_z_3jet_006, 1000.0/sumOfWeights()); + scale(_h_z_3jet_01, 1000.0/sumOfWeights()); + scale(_h_z_4jet_001, 1000.0/sumOfWeights()); + } + + //@} + + + private: + + /// @name Histograms + //@{ + + AIDA::IHistogram1D *_h_z_2jet_001, *_h_z_2jet_006, *_h_z_2jet_01, *_h_z_2jet_033; + AIDA::IHistogram1D *_h_z_3jet_001, *_h_z_3jet_006, *_h_z_3jet_01; + AIDA::IHistogram1D *_h_z_4jet_001; + //@} + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<ALEPH_1996_S3196992> plugin_ALEPH_1996_S3196992; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Sep 24 12:51:39 2009 (r1835) +++ trunk/src/Analyses/Makefile.am Thu Sep 24 16:04:52 2009 (r1836) @@ -66,6 +66,7 @@ RivetLEPAnalyses_la_SOURCES = \ ALEPH_1991_S2435284.cc \ ALEPH_1996_S3486095.cc \ + ALEPH_1996_S3196992.cc \ DELPHI_1995_S3137023.cc \ DELPHI_1996_S3430090.cc \ DELPHI_2002_069_CONF_603.cc \
More information about the Rivet-svn mailing list |