|
[Rivet-svn] r2020 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Nov 5 17:06:39 GMT 2009
Author: fsiegert Date: Thu Nov 5 17:06:38 2009 New Revision: 2020 Log: Add CDF_2001_S4517016 (still untested). Added: trunk/data/anainfo/CDF_2001_S4517016.info trunk/data/plotinfo/CDF_2001_S4517016.plot trunk/data/refdata/CDF_2001_S4517016.aida trunk/src/Analyses/CDF_2001_S4517016.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Added: trunk/data/anainfo/CDF_2001_S4517016.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_2001_S4517016.info Thu Nov 5 17:06:38 2009 (r2020) @@ -0,0 +1,25 @@ +Name: CDF_2001_S4517016 +Year: 2001 +Summary: Two jet triply-differential cross-section +Experiment: CDF +Collider: Tevatron Run 1 +SpiresID: 4517016 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - Phys.Rev.D64:012001,2001 + - doi:10.1103/PhysRevD.64.012001 + - arXiv:hep-ex/0012013 +RunInfo: + Dijet events at Tevatron with $\sqrt{s}=1.8 TeV$ +NumEvents: 1000000 +PtCuts: [0] +Description: + A measurement of the two-jet differential cross section, + $\mathrm{d}^3\sigma/\mathrm{d}E_T \mathrm{d}\eta_1 \mathrm{d}\eta_2$, + based on an integrated luminosity of $86 \mathrm{pb}^-1$. + The differential cross section is measured as a function of the transverse + energy, E_T, of a jet in the pseudorapidity region $0.1 < |\eta_1| < 0.7$ for + four different pseudorapidity bins of a second jet restricted to + $0.1 < |\eta_2| < 3.0$. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Nov 5 17:06:16 2009 (r2019) +++ trunk/data/anainfo/Makefile.am Thu Nov 5 17:06:38 2009 (r2020) @@ -8,6 +8,7 @@ CDF_1996_S3108457.info \ CDF_1996_S3349578.info \ CDF_2000_S4155203.info \ + CDF_2001_S4517016.info \ CDF_2001_S4751469.info \ CDF_2002_S4796047.info \ CDF_2004_S5839831.info \ Added: trunk/data/plotinfo/CDF_2001_S4517016.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_2001_S4517016.plot Thu Nov 5 17:06:38 2009 (r2020) @@ -0,0 +1,27 @@ +# BEGIN PLOT /CDF_2001_S4517016/d01-x01-y01 +Title=$E_\perp$ of leading jet for events with $0.1 < |\eta_2| < 0.7$ +XLabel=$E_\perp$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp$ [nb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /CDF_2001_S4517016/d02-x01-y01 +Title=$E_\perp$ of leading jet for events with $0.7 < |\eta_2| < 1.4$ +XLabel=$E_\perp$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp$ [nb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /CDF_2001_S4517016/d03-x01-y01 +Title=$E_\perp$ of leading jet for events with $1.4 < |\eta_2| < 2.1$ +XLabel=$E_\perp$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp$ [nb/GeV] +LogX=1 +# END PLOT + +# BEGIN PLOT /CDF_2001_S4517016/d04-x01-y01 +Title=$E_\perp$ of leading jet for events with $2.1 < |\eta_2| < 3.0$ +XLabel=$E_\perp$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp$ [nb/GeV] +LogX=1 +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Nov 5 17:06:16 2009 (r2019) +++ trunk/data/plotinfo/Makefile.am Thu Nov 5 17:06:38 2009 (r2020) @@ -8,6 +8,7 @@ CDF_1996_S3108457.plot \ CDF_1996_S3349578.plot \ CDF_2000_S4155203.plot \ + CDF_2001_S4517016.plot \ CDF_2001_S4751469.plot \ CDF_2002_S4796047.plot \ CDF_2004_S5839831.plot \ Added: trunk/data/refdata/CDF_2001_S4517016.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_2001_S4517016.aida Thu Nov 5 17:06:38 2009 (r2020) @@ -0,0 +1,217 @@ +<?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_2001_S4517016" title="D(SIG)/DET IN NB/GEV" > + <dataPoint> + <measurement value="44.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="12.3" errorPlus="2.4083189157584592" errorMinus="2.4083189157584592"/> + </dataPoint> + <dataPoint> + <measurement value="50.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="6.48" errorPlus="1.195073219514185" errorMinus="1.195073219514185"/> + </dataPoint> + <dataPoint> + <measurement value="58.9" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="2.78" errorPlus="0.4816637831516918" errorMinus="0.4816637831516918"/> + </dataPoint> + <dataPoint> + <measurement value="75.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.754" errorPlus="0.1213424904969401" errorMinus="0.1213424904969401"/> + </dataPoint> + <dataPoint> + <measurement value="94.3" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.222" errorPlus="0.03413209633175202" errorMinus="0.03413209633175202"/> + </dataPoint> + <dataPoint> + <measurement value="106.6" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.11" errorPlus="0.0161245154965971" errorMinus="0.0161245154965971"/> + </dataPoint> + <dataPoint> + <measurement value="119.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0583" errorPlus="0.008521150157109076" errorMinus="0.008521150157109076"/> + </dataPoint> + <dataPoint> + <measurement value="132.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0313" errorPlus="0.004517742799230606" errorMinus="0.004517742799230606"/> + </dataPoint> + <dataPoint> + <measurement value="150.8" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0142" errorPlus="0.0021023796041628635" errorMinus="0.0021023796041628635"/> + </dataPoint> + <dataPoint> + <measurement value="174.4" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.00553" errorPlus="8.409518416651456E-4" errorMinus="8.409518416651456E-4"/> + </dataPoint> + <dataPoint> + <measurement value="209.4" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.00167" errorPlus="2.7073972741361766E-4" errorMinus="2.7073972741361766E-4"/> + </dataPoint> + <dataPoint> + <measurement value="264.1" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="3.1E-4" errorPlus="5.755866572463264E-5" errorMinus="5.755866572463264E-5"/> + </dataPoint> + <dataPoint> + <measurement value="318.2" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="6.06E-5" errorPlus="1.3008074415531301E-5" errorMinus="1.3008074415531301E-5"/> + </dataPoint> + <dataPoint> + <measurement value="382.4" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="1.14E-5" errorPlus="3.046309242345563E-6" errorMinus="3.046309242345563E-6"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/CDF_2001_S4517016" title="D(SIG)/DET IN NB/GEV" > + <dataPoint> + <measurement value="43.1" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="12.9" errorPlus="2.80713376952364" errorMinus="2.80713376952364"/> + </dataPoint> + <dataPoint> + <measurement value="49.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="6.41" errorPlus="1.264792473095883" errorMinus="1.264792473095883"/> + </dataPoint> + <dataPoint> + <measurement value="58.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="2.65" errorPlus="0.4716990566028302" errorMinus="0.4716990566028302"/> + </dataPoint> + <dataPoint> + <measurement value="75.2" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.673" errorPlus="0.11031772296417289" errorMinus="0.11031772296417289"/> + </dataPoint> + <dataPoint> + <measurement value="94.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.187" errorPlus="0.029068883707497266" errorMinus="0.029068883707497266"/> + </dataPoint> + <dataPoint> + <measurement value="106.2" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0936" errorPlus="0.014179562757715768" errorMinus="0.014179562757715768"/> + </dataPoint> + <dataPoint> + <measurement value="119.1" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0454" errorPlus="0.0069180922225711905" errorMinus="0.0069180922225711905"/> + </dataPoint> + <dataPoint> + <measurement value="132.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0241" errorPlus="0.003712142238654117" errorMinus="0.003712142238654117"/> + </dataPoint> + <dataPoint> + <measurement value="150.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0103" errorPlus="0.0016031219541881396" errorMinus="0.0016031219541881396"/> + </dataPoint> + <dataPoint> + <measurement value="173.2" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.00385" errorPlus="6.506919393998976E-4" errorMinus="6.506919393998976E-4"/> + </dataPoint> + <dataPoint> + <measurement value="206.9" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="9.92E-4" errorPlus="1.875233318816621E-4" errorMinus="1.875233318816621E-4"/> + </dataPoint> + <dataPoint> + <measurement value="260.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="1.33E-4" errorPlus="3.0413812651491098E-5" errorMinus="3.0413812651491098E-5"/> + </dataPoint> + <dataPoint> + <measurement value="313.7" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="1.98E-5" errorPlus="5.8872744118140105E-6" errorMinus="5.8872744118140105E-6"/> + </dataPoint> + <dataPoint> + <measurement value="373.9" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="3.37E-6" errorPlus="1.3737903770226372E-6" errorMinus="1.3737903770226372E-6"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/CDF_2001_S4517016" title="D(SIG)/DET IN NB/GEV" > + <dataPoint> + <measurement value="42.1" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="11.4" errorPlus="2.6076809620810595" errorMinus="2.6076809620810595"/> + </dataPoint> + <dataPoint> + <measurement value="48.9" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="5.04" errorPlus="1.0547511554864493" errorMinus="1.0547511554864493"/> + </dataPoint> + <dataPoint> + <measurement value="58.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="1.99" errorPlus="0.3920459156782532" errorMinus="0.3920459156782532"/> + </dataPoint> + <dataPoint> + <measurement value="74.3" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.472" errorPlus="0.08811923740024082" errorMinus="0.08811923740024082"/> + </dataPoint> + <dataPoint> + <measurement value="93.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.108" errorPlus="0.02009975124224178" errorMinus="0.02009975124224178"/> + </dataPoint> + <dataPoint> + <measurement value="104.9" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0483" errorPlus="0.008967719888578144" errorMinus="0.008967719888578144"/> + </dataPoint> + <dataPoint> + <measurement value="117.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0221" errorPlus="0.004219004621945797" errorMinus="0.004219004621945797"/> + </dataPoint> + <dataPoint> + <measurement value="130.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0103" errorPlus="0.002009975124224178" errorMinus="0.002009975124224178"/> + </dataPoint> + <dataPoint> + <measurement value="147.4" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.00348" errorPlus="7.310950690573695E-4" errorMinus="7.310950690573695E-4"/> + </dataPoint> + <dataPoint> + <measurement value="169.8" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="9.59E-4" errorPlus="2.2081893034792103E-4" errorMinus="2.2081893034792103E-4"/> + </dataPoint> + <dataPoint> + <measurement value="200.8" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="1.88E-4" errorPlus="4.9497474683058326E-5" errorMinus="4.9497474683058326E-5"/> + </dataPoint> + <dataPoint> + <measurement value="252.7" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="1.06E-5" errorPlus="3.891015291668744E-6" errorMinus="3.891015291668744E-6"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/CDF_2001_S4517016" title="D(SIG)/DET IN NB/GEV" > + <dataPoint> + <measurement value="40.9" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="5.81" errorPlus="1.5945218719101975" errorMinus="1.5945218719101975"/> + </dataPoint> + <dataPoint> + <measurement value="47.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="2.39" errorPlus="0.6029925372672534" errorMinus="0.6029925372672534"/> + </dataPoint> + <dataPoint> + <measurement value="56.2" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.785" errorPlus="0.1875419952970534" errorMinus="0.1875419952970534"/> + </dataPoint> + <dataPoint> + <measurement value="71.7" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.14" errorPlus="0.033955853692699295" errorMinus="0.033955853692699295"/> + </dataPoint> + <dataPoint> + <measurement value="90.4" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.0193" errorPlus="0.004866210024238576" errorMinus="0.004866210024238576"/> + </dataPoint> + <dataPoint> + <measurement value="101.8" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.00747" errorPlus="0.002003521899056758" errorMinus="0.002003521899056758"/> + </dataPoint> + <dataPoint> + <measurement value="114.2" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="0.00228" errorPlus="6.530696746902279E-4" errorMinus="6.530696746902279E-4"/> + </dataPoint> + <dataPoint> + <measurement value="126.0" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="8.13E-4" errorPlus="2.511035642917081E-4" errorMinus="2.511035642917081E-4"/> + </dataPoint> + <dataPoint> + <measurement value="142.5" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="1.89E-4" errorPlus="6.280127387243033E-5" errorMinus="6.280127387243033E-5"/> + </dataPoint> + <dataPoint> + <measurement value="163.7" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="2.39E-5" errorPlus="9.377632963600143E-6" errorMinus="9.377632963600143E-6"/> + </dataPoint> + <dataPoint> + <measurement value="191.4" errorPlus="0.0" errorMinus="0.0"/> + <measurement value="3.64E-6" errorPlus="1.9039432764659771E-6" errorMinus="1.9039432764659771E-6"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Nov 5 17:06:16 2009 (r2019) +++ trunk/data/refdata/Makefile.am Thu Nov 5 17:06:38 2009 (r2020) @@ -36,6 +36,7 @@ CDF_1996_S3108457.aida \ CDF_1996_S3349578.aida \ CDF_2000_S4155203.aida \ + CDF_2001_S4517016.aida \ CDF_2001_S4751469.aida \ CDF_2002_S4796047.aida \ CDF_2004_S5839831.aida \ Added: trunk/src/Analyses/CDF_2001_S4517016.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_2001_S4517016.cc Thu Nov 5 17:06:38 2009 (r2020) @@ -0,0 +1,96 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Tools/BinnedHistogram.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/FastJets.hh" + +namespace Rivet { + + + class CDF_2001_S4517016 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + CDF_2001_S4517016() + : Analysis("CDF_2001_S4517016") + { + setBeams(PROTON, ANTIPROTON); + setNeedsCrossSection(true); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + FinalState fs(-4.2, 4.2); + addProjection(FastJets(fs, FastJets::CDFJETCLU, 0.7, 10.0*GeV), "Jets"); + + _h_ET.addHistogram(0.1, 0.7, bookHistogram1D(1, 1, 1)); + _h_ET.addHistogram(0.7, 1.4, bookHistogram1D(2, 1, 1)); + _h_ET.addHistogram(1.4, 2.1, bookHistogram1D(3, 1, 1)); + _h_ET.addHistogram(2.1, 3.0, bookHistogram1D(4, 1, 1)); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + Jets jets = applyProjection<FastJets>(event, "Jets").jetsByEt(10.0*GeV); + if (jets.size()<2) { + vetoEvent; + } + FourMomentum jet1 = jets[0].momentum(); + FourMomentum jet2 = jets[1].momentum(); + double eta1 = fabs(jet1.eta()); + double eta2 = fabs(jet2.eta()); + double ET1 = jet1.Et(); + double ET2 = jet2.Et(); + if (eta1<0.1 || eta1>0.7 || ET1<40.0*GeV) { + vetoEvent; + } + if (eta2<0.1 || eta2>3.0) { + vetoEvent; + } + _h_ET.fill(eta2, ET1, weight); + if (eta2<0.7 && ET2>40.0*GeV) _h_ET.fill(eta1, ET2, weight); + } + + + /// Normalise histograms etc., after the run + void finalize() { + foreach (AIDA::IHistogram1D* hist, _h_ET.getHistograms()) { + scale(hist, crossSection()/nanobarn/sumOfWeights()); + } + } + + //@} + + + private: + + /// @name Histograms + //@{ + BinnedHistogram<double> _h_ET; + //@} + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<CDF_2001_S4517016> plugin_CDF_2001_S4517016; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Nov 5 17:06:16 2009 (r2019) +++ trunk/src/Analyses/Makefile.am Thu Nov 5 17:06:38 2009 (r2020) @@ -17,6 +17,7 @@ CDF_1996_S3108457.cc \ CDF_1996_S3349578.cc \ CDF_2000_S4155203.cc \ + CDF_2001_S4517016.cc \ CDF_2001_S4751469.cc \ CDF_2002_S4796047.cc \ CDF_2004_S5839831.cc \
More information about the Rivet-svn mailing list |