|
[Rivet-svn] r1898 - in trunk: data/anainfo data/plotinfo data/refdata src src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgFri Oct 9 17:23:24 BST 2009
Author: fsiegert Date: Fri Oct 9 17:23:23 2009 New Revision: 1898 Log: Finish CDF_2009_S8383952 Added: trunk/data/refdata/CDF_2009_S8383952.aida Modified: trunk/data/anainfo/CDF_2009_S8383952.info trunk/data/plotinfo/CDF_2009_S8383952.plot trunk/data/refdata/Makefile.am trunk/src/Analyses/CDF_2009_S8383952.cc trunk/src/Makefile.am Modified: trunk/data/anainfo/CDF_2009_S8383952.info ============================================================================== --- trunk/data/anainfo/CDF_2009_S8383952.info Fri Oct 9 09:53:07 2009 (r1897) +++ trunk/data/anainfo/CDF_2009_S8383952.info Fri Oct 9 17:23:23 2009 (r1898) @@ -4,7 +4,7 @@ Experiment: CDF Collider: Tevatron Run 2 SpiresID: 8383952 -Status: UNFINISHED +Status: VALIDATED Authors: - Frank Siegert <frank.siegert at durham.ac.uk> References: Modified: trunk/data/plotinfo/CDF_2009_S8383952.plot ============================================================================== --- trunk/data/plotinfo/CDF_2009_S8383952.plot Fri Oct 9 09:53:07 2009 (r1897) +++ trunk/data/plotinfo/CDF_2009_S8383952.plot Fri Oct 9 17:23:23 2009 (r1898) @@ -1,5 +1,19 @@ # BEGIN PLOT /CDF_2009_S8383952/d01-x01-y01 +Title=Total XS for $66 < M_{ee}/\mathrm{GeV} < 116$ +XLabel=$\sqrt{s}$ +YLabel=$\sigma$ [pb] +LogY=0 +XMajorTickMarks=20 +XMinorTickMarks=0 +YMin=150 +YMax=350 +XMin=1955 +XMax=1965 +# END PLOT + +# BEGIN PLOT /CDF_2009_S8383952/d02-x01-y01 Title=$e^+ e^-$ pair rapidity XLabel=$|y_{e^+ e^-}|$ YLabel=$\mathrm{d}\sigma/\mathrm{d}y$ [pb] +LogY=0 # END PLOT Added: trunk/data/refdata/CDF_2009_S8383952.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_2009_S8383952.aida Fri Oct 9 17:23:23 2009 (r1898) @@ -0,0 +1,129 @@ +<?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_2009_S8383952" title="SIG IN PB" > + <dataPoint> + <measurement value="1960.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="256.0" errorPlus="2.118962010041709" errorMinus="2.118962010041709"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/CDF_2009_S8383952" title="D(SIG)/YRAP IN PB" > + <dataPoint> + <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/> + <measurement value="69.41" errorPlus="0.8792041856133306" errorMinus="0.8792041856133306"/> + </dataPoint> + <dataPoint> + <measurement value="0.15" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/> + <measurement value="70.94" errorPlus="0.8875246475450695" errorMinus="0.8875246475450695"/> + </dataPoint> + <dataPoint> + <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="70.86" errorPlus="0.8792041856133306" errorMinus="0.8792041856133306"/> + </dataPoint> + <dataPoint> + <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/> + <measurement value="69.68" errorPlus="0.8653323061113574" errorMinus="0.8653323061113574"/> + </dataPoint> + <dataPoint> + <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="67.62" errorPlus="0.8431488599292535" errorMinus="0.8431488599292535"/> + </dataPoint> + <dataPoint> + <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="67.79" errorPlus="0.8431488599292535" errorMinus="0.8431488599292535"/> + </dataPoint> + <dataPoint> + <measurement value="0.65" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="66.07" errorPlus="0.834865258589672" errorMinus="0.834865258589672"/> + </dataPoint> + <dataPoint> + <measurement value="0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="66.19" errorPlus="0.840535543567314" errorMinus="0.840535543567314"/> + </dataPoint> + <dataPoint> + <measurement value="0.85" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/> + <measurement value="64.31" errorPlus="0.8381527307120106" errorMinus="0.8381527307120106"/> + </dataPoint> + <dataPoint> + <measurement value="0.95" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/> + <measurement value="63.88" errorPlus="0.8381527307120106" errorMinus="0.8381527307120106"/> + </dataPoint> + <dataPoint> + <measurement value="1.05" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="61.91" errorPlus="0.8280096617793781" errorMinus="0.8280096617793781"/> + </dataPoint> + <dataPoint> + <measurement value="1.15" errorPlus="0.050000000000000044" errorMinus="0.04999999999999982"/> + <measurement value="61.27" errorPlus="0.814002456998749" errorMinus="0.814002456998749"/> + </dataPoint> + <dataPoint> + <measurement value="1.25" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="58.21" errorPlus="0.802122185206219" errorMinus="0.802122185206219"/> + </dataPoint> + <dataPoint> + <measurement value="1.35" errorPlus="0.04999999999999982" errorMinus="0.050000000000000044"/> + <measurement value="55.59" errorPlus="0.7655063683601855" errorMinus="0.7655063683601855"/> + </dataPoint> + <dataPoint> + <measurement value="1.45" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="53.11" errorPlus="0.7462573282722255" errorMinus="0.7462573282722255"/> + </dataPoint> + <dataPoint> + <measurement value="1.55" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="49.91" errorPlus="0.7220110802473879" errorMinus="0.7220110802473879"/> + </dataPoint> + <dataPoint> + <measurement value="1.65" errorPlus="0.050000000000000044" errorMinus="0.04999999999999982"/> + <measurement value="46.48" errorPlus="0.6946221994724903" errorMinus="0.6946221994724903"/> + </dataPoint> + <dataPoint> + <measurement value="1.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="40.83" errorPlus="0.6673080248281149" errorMinus="0.6673080248281149"/> + </dataPoint> + <dataPoint> + <measurement value="1.85" errorPlus="0.04999999999999982" errorMinus="0.050000000000000044"/> + <measurement value="37.03" errorPlus="0.644980619863884" errorMinus="0.644980619863884"/> + </dataPoint> + <dataPoint> + <measurement value="1.95" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="33.21" errorPlus="0.6313477647065839" errorMinus="0.6313477647065839"/> + </dataPoint> + <dataPoint> + <measurement value="2.05" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/> + <measurement value="28.05" errorPlus="0.5860034129593445" errorMinus="0.5860034129593445"/> + </dataPoint> + <dataPoint> + <measurement value="2.15" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/> + <measurement value="22.43" errorPlus="0.563648826841678" errorMinus="0.563648826841678"/> + </dataPoint> + <dataPoint> + <measurement value="2.25" errorPlus="0.04999999999999982" errorMinus="0.04999999999999982"/> + <measurement value="19.07" errorPlus="0.551543289325507" errorMinus="0.551543289325507"/> + </dataPoint> + <dataPoint> + <measurement value="2.35" errorPlus="0.04999999999999982" errorMinus="0.050000000000000266"/> + <measurement value="14.85" errorPlus="0.5502726596879042" errorMinus="0.5502726596879042"/> + </dataPoint> + <dataPoint> + <measurement value="2.45" errorPlus="0.04999999999999982" errorMinus="0.050000000000000266"/> + <measurement value="9.4" errorPlus="0.49244289008980524" errorMinus="0.49244289008980524"/> + </dataPoint> + <dataPoint> + <measurement value="2.55" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/> + <measurement value="6.16" errorPlus="0.4903060268852505" errorMinus="0.4903060268852505"/> + </dataPoint> + <dataPoint> + <measurement value="2.65" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/> + <measurement value="3.27" errorPlus="0.49040799340956914" errorMinus="0.49040799340956914"/> + </dataPoint> + <dataPoint> + <measurement value="2.75" errorPlus="0.04999999999999982" errorMinus="0.04999999999999982"/> + <measurement value="1.48" errorPlus="0.41593268686170837" errorMinus="0.41593268686170837"/> + </dataPoint> + <dataPoint> + <measurement value="2.85" errorPlus="0.04999999999999982" errorMinus="0.050000000000000266"/> + <measurement value="0.96" errorPlus="0.5770615218501403" errorMinus="0.5770615218501403"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Fri Oct 9 09:53:07 2009 (r1897) +++ trunk/data/refdata/Makefile.am Fri Oct 9 17:23:23 2009 (r1898) @@ -52,6 +52,7 @@ CDF_2008_NOTE_9351.aida \ CDF_2008_LEADINGJETS.aida \ CDF_2009_S8233977.aida \ + CDF_2009_S8383952.aida \ D0_1996_S3214044.aida \ D0_1996_S3324664.aida \ D0_2001_S4674421.aida \ Modified: trunk/src/Analyses/CDF_2009_S8383952.cc ============================================================================== --- trunk/src/Analyses/CDF_2009_S8383952.cc Fri Oct 9 09:53:07 2009 (r1897) +++ trunk/src/Analyses/CDF_2009_S8383952.cc Fri Oct 9 17:23:23 2009 (r1898) @@ -34,13 +34,16 @@ void init() { /// Initialise and register projections here - ZFinder zfinder(-2.8, 2.8, 20.0*GeV, ELECTRON, + // this seems to have been corrected completely for all selection cuts, + // i.e. eta cuts and pT cuts on leptons. + ZFinder zfinder(-MAXRAPIDITY, MAXRAPIDITY, 0.0*GeV, ELECTRON, 66.0*GeV, 116.0*GeV, 0.2); addProjection(zfinder, "ZFinder"); - /// Book histograms here, e.g.: - _h_yZ = bookHistogram1D(1, 1, 1); + /// Book histograms here + _h_xs = bookHistogram1D(1, 1, 1); + _h_yZ = bookHistogram1D(2, 1, 1); } @@ -53,6 +56,7 @@ if (zfinder.particles().size() == 1) { double yZ = fabs(zfinder.particles()[0].momentum().rapidity()); _h_yZ->fill(yZ, weight); + _h_xs->fill(1960.0, weight); } else { getLog() << Log::DEBUG << "no unique lepton pair found." << endl; @@ -63,7 +67,10 @@ /// Normalise histograms etc., after the run void finalize() { - scale(_h_yZ, crossSection()/sumOfWeights()); + scale(_h_xs, crossSection()/sumOfWeights()); + // Data seems to have been normalized for the avg of the two sides + // (+ve & -ve rapidity) rather than the sum, hence the 0.5: + scale(_h_yZ, 0.5*crossSection()/sumOfWeights()); } //@} @@ -73,8 +80,8 @@ /// @name Histograms //@{ - AIDA::IHistogram1D *_h_yZ; + AIDA::IHistogram1D *_h_xs; //@} }; Modified: trunk/src/Makefile.am ============================================================================== --- trunk/src/Makefile.am Fri Oct 9 09:53:07 2009 (r1897) +++ trunk/src/Makefile.am Fri Oct 9 17:23:23 2009 (r1898) @@ -1,4 +1,4 @@ -SUBDIRS = Projections +SUBDIRS = . Projections if ENABLE_ANALYSES SUBDIRS += Analyses endif
More information about the Rivet-svn mailing list |