|
[Rivet-svn] r2475 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgWed Jun 9 22:21:05 BST 2010
Author: fsiegert Date: Wed Jun 9 22:21:03 2010 New Revision: 2475 Log: Add OPAL_2001_S4553896. Has been validated using Sherpa events. Closes ticket #405. Added: trunk/data/anainfo/OPAL_2001_S4553896.info trunk/data/plotinfo/OPAL_2001_S4553896.plot trunk/data/refdata/OPAL_2001_S4553896.aida trunk/src/Analyses/OPAL_2001_S4553896.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Tue Jun 8 23:44:05 2010 (r2474) +++ trunk/data/anainfo/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) @@ -78,6 +78,7 @@ MC_WJETS.info \ MC_ZJETS.info \ OPAL_1998_S3780481.info \ + OPAL_2001_S4553896.info \ OPAL_2004_S6132243.info \ PDG_HADRON_MULTIPLICITIES.info \ PDG_HADRON_MULTIPLICITIES_RATIOS.info \ Added: trunk/data/anainfo/OPAL_2001_S4553896.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/OPAL_2001_S4553896.info Wed Jun 9 22:21:03 2010 (r2475) @@ -0,0 +1,44 @@ +Name: OPAL_2001_S4553896 +Year: 2001 +Summary: Four-jet angles using Durham algorithm +Experiment: OPAL +Collider: LEP Run 1 +SpiresID: 4553896 +Status: VALIDATED +Authors: + - Frank Siegert <frank.siegert at cern.ch> +References: + - Eur.Phys.J.C20:601-615,2001 + - doi:10.1007/s100520100699 + - arXiv:hep-ex/0101044 +RunInfo: + Hadronic Z decay events generated on the Z pole (sqrt(s) = 91.2 GeV) + Hadronisation should be turned off because the data is corrected back to the + parton level. +NumEvents: 10000000 +Beams: [e+, e-] +Energies: [91.2] +PtCuts: [0] +Description: + 'Angles between the leading (in energy) four jets defined using the Durham + algorithm with $y_\mathrm{cut}=0.008$. The data is presented at the parton + level and includes the Bengtsson-Zerwas, Korner-Schierholz-Willrodt and + Nachtmann-Reiter angles as well as the angle between the two softest jets.' +BibKey: Abbiendi:2001qn +BibTeX: '@Article{Abbiendi:2001qn, + author = "Abbiendi, G. and others", + collaboration = "OPAL", + title = "{A simultaneous measurement of the QCD colour factors and + the strong coupling}", + journal = "Eur. Phys. J.", + volume = "C20", + year = "2001", + pages = "601-615", + eprint = "hep-ex/0101044", + archivePrefix = "arXiv", + doi = "10.1007/s100520100699", + SLACcitation = "%%CITATION = HEP-EX/0101044;%%" +}' +ToDo: + - Add y23 and R4(ycut) histograms + - Get data from HepData when available (so far just read off the plots) Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Tue Jun 8 23:44:05 2010 (r2474) +++ trunk/data/plotinfo/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) @@ -76,6 +76,7 @@ MC_WJETS.plot \ MC_ZJETS.plot \ OPAL_1998_S3780481.plot \ + OPAL_2001_S4553896.plot \ OPAL_2004_S6132243.plot \ PDG_HADRON_MULTIPLICITIES.plot \ PDG_HADRON_MULTIPLICITIES_RATIOS.plot \ Added: trunk/data/plotinfo/OPAL_2001_S4553896.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/OPAL_2001_S4553896.plot Wed Jun 9 22:21:03 2010 (r2475) @@ -0,0 +1,38 @@ +# BEGIN PLOT /OPAL_2001_S4553896/d03-x01-y01 +Title=Bengtsson-Zerwas angle (parton level) +XLabel=$|\cos(\chi_\text{BZ})|$ +YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}|\cos(\chi_\text{BZ})|$ +LogY=0 +RatioPlotYMin=0.9 +RatioPlotYMax=1.1 +# END PLOT + +# BEGIN PLOT /OPAL_2001_S4553896/d04-x01-y01 +Title=K\"orner-Schierholz-Willrodt angle (parton level) +XLabel=$\cos(\phi_\text{KSW})$ +YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}\,\cos(\phi_\text{KSW})$ +LogY=0 +RatioPlotYMin=0.9 +RatioPlotYMax=1.1 +# END PLOT + +# BEGIN PLOT /OPAL_2001_S4553896/d05-x01-y01 +Title=Modified Nachtmann-Reiter angle (parton level) +XLabel=$|\cos(\theta^*_\text{NR})|$ +YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}|\cos(\theta^*_\text{NR})|$ +LogY=0 +LegendXPos=0.5 +RatioPlotYMin=0.9 +RatioPlotYMax=1.1 +# END PLOT + +# BEGIN PLOT /OPAL_2001_S4553896/d06-x01-y01 +Title=Angle between the two softest jets (parton level) +XLabel=$\cos(\alpha_{34})$ +YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}\,\cos(\alpha_{34})$ +LogY=0 +LegendXPos=0.5 +LegendYPos=0.3 +RatioPlotYMin=0.9 +RatioPlotYMax=1.1 +# END PLOT Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Tue Jun 8 23:44:05 2010 (r2474) +++ trunk/data/refdata/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) @@ -11,6 +11,7 @@ DELPHI_2002_069_CONF_603.aida \ DELPHI_2003_WUD_03_11.aida \ OPAL_2004_S6132243.aida \ + OPAL_2001_S4553896.aida \ OPAL_1998_S3780481.aida \ CLEO_1997_S3486664.aida \ CLEO_2001_S4557530.aida \ Added: trunk/data/refdata/OPAL_2001_S4553896.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/OPAL_2001_S4553896.aida Wed Jun 9 22:21:03 2010 (r2475) @@ -0,0 +1,365 @@ +<?xml version="1.0" encoding="ISO-8859-1" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.1" package="FreeHEP"/> +<dataPointSet name="d03-x01-y01" dimension="2" + path="/REF/OPAL_2001_S4553896" title="BZ"> + <dimension dim="0" title="BZ" /> + <dimension dim="1" title="DSIGMADBZ" /> + <annotation> + <item key="Title" value="BZ" sticky="true"/> + <item key="AidaPath" value="/REF/OPAL_2001_S4553896/d03-x01-y01" sticky="true"/> + <item key="FullPath" value="/OPAL_2001_S4553896.aida/REF/OPAL_2001_S4553896/d03-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.500000e-02" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.305936e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.500000e-02" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.351598e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.465753e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.511416e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.579909e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.808219e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.785388e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.694064e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.922374e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="8.127854e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="8.630137e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="8.584475e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="9.360731e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="9.406393e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.027397e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.036530e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.125571e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.230594e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.420091e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="2.947489e+00" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" + path="/REF/OPAL_2001_S4553896" title="KSW"> + <dimension dim="0" title="KSW" /> + <dimension dim="1" title="DSIGMADKSW" /> + <annotation> + <item key="Title" value="KSW" sticky="true"/> + <item key="AidaPath" value="/REF/OPAL_2001_S4553896/d04-x01-y01" sticky="true"/> + <item key="FullPath" value="/OPAL_2001_S4553896.aida/REF/OPAL_2001_S4553896/d04-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-9.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.258797e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-8.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.499673e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-7.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.383293e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-6.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.500589e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-5.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.680842e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-4.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.689954e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-3.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.816260e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-2.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.960522e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-1.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.916005e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-5.000000e-02" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.916473e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="5.000000e-02" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.341645e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="1.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.189577e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="2.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.144808e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="3.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="3.974736e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="4.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="3.956923e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="5.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="3.930718e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="6.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.182532e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="7.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.488576e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="8.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.486797e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="9.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.197716e+00" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" dimension="2" + path="/REF/OPAL_2001_S4553896" title="NR"> + <dimension dim="0" title="NR" /> + <dimension dim="1" title="DSIGMADNR" /> + <annotation> + <item key="Title" value="NR" sticky="true"/> + <item key="AidaPath" value="/REF/OPAL_2001_S4553896/d05-x01-y01" sticky="true"/> + <item key="FullPath" value="/OPAL_2001_S4553896.aida/REF/OPAL_2001_S4553896/d05-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.500000e-02" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.856180e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.500000e-02" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.829213e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.955056e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="1.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="8.134831e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.865169e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="2.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="8.332584e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="8.386517e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="3.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="8.647191e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="9.078652e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="4.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="9.258427e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="9.770787e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="5.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="9.986517e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.058876e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="6.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.084944e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.159551e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="7.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.173933e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.190112e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="8.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.204494e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.250000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.244944e+00" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.500000e-02" value="9.750000e-01" errorMinus="2.500000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.550562e+00" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y01" dimension="2" + path="/REF/OPAL_2001_S4553896" title="ALPHA34"> + <dimension dim="0" title="ALPHA34" /> + <dimension dim="1" title="DSIGMADALPHA34" /> + <annotation> + <item key="Title" value="ALPHA34" sticky="true"/> + <item key="AidaPath" value="/REF/OPAL_2001_S4553896/d06-x01-y01" sticky="true"/> + <item key="FullPath" value="/OPAL_2001_S4553896.aida/REF/OPAL_2001_S4553896/d06-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-9.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="6.894629e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-8.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="6.140462e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-7.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="6.022773e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-6.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.731969e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-5.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.766183e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-4.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.558160e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-3.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.592592e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-2.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.329236e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-1.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.156239e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="-5.000000e-02" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.349310e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="5.000000e-02" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.715740e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="1.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.756854e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="2.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="6.053874e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="3.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.776920e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="4.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.506878e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="5.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="5.167419e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="6.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="4.468826e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="7.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="3.230592e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="8.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="1.114061e-01" errorMinus="0.000000e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-02" value="9.500000e-01" errorMinus="5.000000e-02"/> + <measurement errorPlus="0.000000e+00" value="7.176745e-04" errorMinus="0.000000e+00"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Tue Jun 8 23:44:05 2010 (r2474) +++ trunk/src/Analyses/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) @@ -105,6 +105,7 @@ DELPHI_1996_S3430090.cc \ DELPHI_2002_069_CONF_603.cc \ OPAL_1998_S3780481.cc \ + OPAL_2001_S4553896.cc \ OPAL_2004_S6132243.cc Added: trunk/src/Analyses/OPAL_2001_S4553896.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/OPAL_2001_S4553896.cc Wed Jun 9 22:21:03 2010 (r2475) @@ -0,0 +1,159 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FastJets.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + namespace { + + /// @name Jet angle calculator functions + //@{ + /// @todo Move to utils? (taken from DELPHI_2003) + + /// @todo Use Jet or FourMomentum interface rather than PseudoJet + /// @todo Move to utils? + double calc_BZ(const vector<fastjet::PseudoJet>& jets) { + assert(jets.size() == 4); + Vector3 p12 = cross( momentum3(jets[0]), momentum3(jets[1])); + Vector3 p34 = cross( momentum3(jets[2]), momentum3(jets[3])); + return dot(p12,p34) / (p12.mod()*p34.mod()); + } + + + /// @todo Use Jet or FourMomentum interface rather than PseudoJet + /// @todo Move to utils? + double calc_KSW(const vector<fastjet::PseudoJet>& jets) { + assert(jets.size() == 4); + Vector3 p13 = cross( momentum3(jets[0]), momentum3(jets[2])); + Vector3 p24 = cross( momentum3(jets[1]), momentum3(jets[3])); + Vector3 p14 = cross( momentum3(jets[0]), momentum3(jets[3])); + Vector3 p23 = cross( momentum3(jets[1]), momentum3(jets[2])); + return cos (0.5*( acos (dot(p14,p23) / (p14.mod()*p23.mod())) + + acos (dot(p13,p24) / (p13.mod()*p24.mod())) )); + } + + + /// @todo Use Jet or FourMomentum interface rather than PseudoJet + /// @todo Move to utils? + double calc_NR(const vector<fastjet::PseudoJet>& jets) { + assert(jets.size() == 4); + Vector3 p12 = momentum3(jets[0]) - momentum3(jets[1]); + Vector3 p34 = momentum3(jets[2]) - momentum3(jets[3]); + return dot(p12,p34) / (p12.mod()*p34.mod()); + } + + /// @todo Use Jet or FourMomentum interface rather than PseudoJet + /// @todo Move to utils? + double calc_ALPHA34(const vector<fastjet::PseudoJet>& jets) { + assert(jets.size() == 4); + Vector3 p3 = momentum3(jets[2]); + Vector3 p4 = momentum3(jets[3]); + return dot(p3,p4) / (p3.mod()*p4.mod()); + } + + //@} + + } + + + class OPAL_2001_S4553896 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + OPAL_2001_S4553896() + : Analysis("OPAL_2001_S4553896") + { + setBeams(ELECTRON, POSITRON); + setNeedsCrossSection(false); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + /// Initialise and register projections here + const FinalState fs; + addProjection(fs, "FS"); + addProjection(FastJets(fs, FastJets::DURHAM, 0.7), "Jets"); + + + /// @todo Book histograms here, e.g.: + _h_BZ = bookHistogram1D(3, 1, 1); + _h_KSW = bookHistogram1D(4, 1, 1); + _h_NR = bookHistogram1D(5, 1, 1); + _h_ALPHA34 = bookHistogram1D(6, 1, 1); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + // Even if we only generate hadronic events, we still need a cut on numCharged >= 2. + if (applyProjection<FinalState>(event, "FS").particles().size() < 2) { + vetoEvent; + } + + const FastJets& fastjets = applyProjection<FastJets>(event, "Jets"); + vector<fastjet::PseudoJet> jets; + if (fastjets.clusterSeq()) { + jets = fastjet::sorted_by_E(fastjets.clusterSeq()->exclusive_jets_ycut(0.008)); + if (jets.size() == 4) { + _h_BZ->fill(fabs(calc_BZ(jets)), weight); + _h_KSW->fill(calc_KSW(jets), weight); + _h_NR->fill(fabs(calc_NR(jets)), weight); + _h_ALPHA34->fill(calc_ALPHA34(jets), weight); + } + } + + + } + + + /// Normalise histograms etc., after the run + void finalize() { + + /// Normalise, scale and otherwise manipulate histograms here + normalize(_h_BZ); + normalize(_h_KSW); + normalize(_h_NR); + normalize(_h_ALPHA34); + + } + + //@} + + + private: + + /// @name Histograms + //@{ + + AIDA::IHistogram1D *_h_BZ; + AIDA::IHistogram1D *_h_KSW; + AIDA::IHistogram1D *_h_NR; + AIDA::IHistogram1D *_h_ALPHA34; + //@} + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<OPAL_2001_S4553896> plugin_OPAL_2001_S4553896; + + +}
More information about the Rivet-svn mailing list |