|
[Rivet-svn] r2024 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Nov 5 18:03:30 GMT 2009
Author: fsiegert Date: Thu Nov 5 18:03:30 2009 New Revision: 2024 Log: Add CDF_2000_S4266730 (not tested yet). Added: trunk/data/anainfo/CDF_2000_S4266730.info trunk/data/plotinfo/CDF_2000_S4266730.plot trunk/data/refdata/CDF_2000_S4266730.aida trunk/src/Analyses/CDF_2000_S4266730.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_2000_S4266730.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_2000_S4266730.info Thu Nov 5 18:03:30 2009 (r2024) @@ -0,0 +1,21 @@ +Name: CDF_2000_S4266730 +Year: 2000 +Summary: Differential Dijet Mass Cross Section +Experiment: CDF +Collider: Tevatron Run 1 +SpiresID: 4266730 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - Phys.Rev.D61:091101,2000 + - 10.1103/PhysRevD.61.091101 + - arXiv:hep-ex/9912022 +RunInfo: + Dijet events at Tevatron with $\sqrt{s}=1.8 TeV$ +NumEvents: 1000000 +PtCuts: [0] +Description: + Measurement of the cross section for production of two or more jets as a + function of dijet mass in the range 180 to 1000 GeV. + It is based on an integrated luminosity of $86 \mathrm{pb}^-1$. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Nov 5 17:45:35 2009 (r2023) +++ trunk/data/anainfo/Makefile.am Thu Nov 5 18:03:30 2009 (r2024) @@ -8,6 +8,7 @@ CDF_1996_S3108457.info \ CDF_1996_S3349578.info \ CDF_2000_S4155203.info \ + CDF_2000_S4266730.info \ CDF_2001_S4517016.info \ CDF_2001_S4563131.info \ CDF_2001_S4751469.info \ Added: trunk/data/plotinfo/CDF_2000_S4266730.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_2000_S4266730.plot Thu Nov 5 18:03:30 2009 (r2024) @@ -0,0 +1,7 @@ +# BEGIN PLOT /CDF_2000_S4266730/d01-x01-y01 +Title=Dijet mass +XLabel=$m_{jj}$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}m$ [pb/GeV] +LogX=1 +# END PLOT + Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Nov 5 17:45:35 2009 (r2023) +++ trunk/data/plotinfo/Makefile.am Thu Nov 5 18:03:30 2009 (r2024) @@ -8,6 +8,7 @@ CDF_1996_S3108457.plot \ CDF_1996_S3349578.plot \ CDF_2000_S4155203.plot \ + CDF_2000_S4266730.plot \ CDF_2001_S4517016.plot \ CDF_2001_S4563131.plot \ CDF_2001_S4751469.plot \ Added: trunk/data/refdata/CDF_2000_S4266730.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_2000_S4266730.aida Thu Nov 5 18:03:30 2009 (r2024) @@ -0,0 +1,79 @@ +<?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_2000_S4266730" title="D(SIG)/DM IN PB/GEV" > + <dataPoint> + <measurement value="188.0" errorPlus="10.0" errorMinus="8.0"/> + <measurement value="607.0" errorPlus="104.7377677822093" errorMinus="104.7377677822093"/> + </dataPoint> + <dataPoint> + <measurement value="207.0" errorPlus="10.0" errorMinus="9.0"/> + <measurement value="342.0" errorPlus="59.665735560705194" errorMinus="59.665735560705194"/> + </dataPoint> + <dataPoint> + <measurement value="228.0" errorPlus="13.0" errorMinus="11.0"/> + <measurement value="181.0" errorPlus="29.068883707497267" errorMinus="29.068883707497267"/> + </dataPoint> + <dataPoint> + <measurement value="252.0" errorPlus="13.0" errorMinus="11.0"/> + <measurement value="98.1" errorPlus="15.762296786953353" errorMinus="15.762296786953353"/> + </dataPoint> + <dataPoint> + <measurement value="277.0" errorPlus="15.0" errorMinus="12.0"/> + <measurement value="49.8" errorPlus="8.54751425854324" errorMinus="8.54751425854324"/> + </dataPoint> + <dataPoint> + <measurement value="305.0" errorPlus="16.0" errorMinus="13.0"/> + <measurement value="27.8" errorPlus="4.709564735726647" errorMinus="4.709564735726647"/> + </dataPoint> + <dataPoint> + <measurement value="335.0" errorPlus="18.0" errorMinus="14.0"/> + <measurement value="14.3" errorPlus="2.4083189157584592" errorMinus="2.4083189157584592"/> + </dataPoint> + <dataPoint> + <measurement value="368.0" errorPlus="20.0" errorMinus="15.0"/> + <measurement value="7.41" errorPlus="1.3373481222179961" errorMinus="1.3373481222179961"/> + </dataPoint> + <dataPoint> + <measurement value="405.0" errorPlus="22.0" errorMinus="17.0"/> + <measurement value="3.83" errorPlus="0.6906518659932802" errorMinus="0.6906518659932802"/> + </dataPoint> + <dataPoint> + <measurement value="446.0" errorPlus="24.0" errorMinus="19.0"/> + <measurement value="1.89" errorPlus="0.36055512754639896" errorMinus="0.36055512754639896"/> + </dataPoint> + <dataPoint> + <measurement value="491.0" errorPlus="26.0" errorMinus="21.0"/> + <measurement value="0.907" errorPlus="0.1726528308484978" errorMinus="0.1726528308484978"/> + </dataPoint> + <dataPoint> + <measurement value="539.0" errorPlus="29.0" errorMinus="22.0"/> + <measurement value="0.45" errorPlus="0.09055385138137416" errorMinus="0.09055385138137416"/> + </dataPoint> + <dataPoint> + <measurement value="592.0" errorPlus="33.0" errorMinus="24.0"/> + <measurement value="0.19" errorPlus="0.04044749683231337" errorMinus="0.04044749683231337"/> + </dataPoint> + <dataPoint> + <measurement value="652.0" errorPlus="36.0" errorMinus="27.0"/> + <measurement value="0.0742" errorPlus="0.016737084572887834" errorMinus="0.016737084572887834"/> + </dataPoint> + <dataPoint> + <measurement value="716.0" errorPlus="40.0" errorMinus="28.0"/> + <measurement value="0.0292" errorPlus="0.007051950084905593" errorMinus="0.007051950084905593"/> + </dataPoint> + <dataPoint> + <measurement value="784.0" errorPlus="48.0" errorMinus="28.0"/> + <measurement value="0.0118" errorPlus="0.0032695565448543633" errorMinus="0.0032695565448543633"/> + </dataPoint> + <dataPoint> + <measurement value="865.0" errorPlus="50.0" errorMinus="33.0"/> + <measurement value="0.00357" errorPlus="0.0011700427342623003" errorMinus="0.0011700427342623003"/> + </dataPoint> + <dataPoint> + <measurement value="968.0" errorPlus="32.0" errorMinus="53.0"/> + <measurement value="9.03E-4" errorPlus="3.909130338067535E-4" errorMinus="3.909130338067535E-4"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Nov 5 17:45:35 2009 (r2023) +++ trunk/data/refdata/Makefile.am Thu Nov 5 18:03:30 2009 (r2024) @@ -36,6 +36,7 @@ CDF_1996_S3108457.aida \ CDF_1996_S3349578.aida \ CDF_2000_S4155203.aida \ + CDF_2000_S4266730.aida \ CDF_2001_S4517016.aida \ CDF_2001_S4563131.aida \ CDF_2001_S4751469.aida \ Added: trunk/src/Analyses/CDF_2000_S4266730.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_2000_S4266730.cc Thu Nov 5 18:03:30 2009 (r2024) @@ -0,0 +1,91 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/FastJets.hh" + +namespace Rivet { + + + class CDF_2000_S4266730 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + CDF_2000_S4266730() + : Analysis("CDF_2000_S4266730") + { + 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, 0.0*GeV), "Jets"); + + _h_mjj = bookHistogram1D(1, 1, 1); + + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + Jets jets = applyProjection<FastJets>(event, "Jets").jetsByEt(0.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()); + if (eta1>2.0 || eta2>2.0) { + vetoEvent; + } + if (fabs(tanh((eta1-eta2)/2))>2.0/3.0) { + vetoEvent; + } + double mjj = FourMomentum(jet1+jet2).mass()/GeV; + if (mjj<180.0) { + vetoEvent; + } + _h_mjj->fill(mjj, weight); + + } + + + /// Normalise histograms etc., after the run + void finalize() { + scale(_h_mjj, crossSection()/picobarn/sumOfWeights()); + } + + //@} + private: + + /// @name Histograms + //@{ + AIDA::IHistogram1D *_h_mjj; + //@} + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<CDF_2000_S4266730> plugin_CDF_2000_S4266730; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Nov 5 17:45:35 2009 (r2023) +++ trunk/src/Analyses/Makefile.am Thu Nov 5 18:03:30 2009 (r2024) @@ -17,6 +17,7 @@ CDF_1996_S3108457.cc \ CDF_1996_S3349578.cc \ CDF_2000_S4155203.cc \ + CDF_2000_S4266730.cc \ CDF_2001_S4517016.cc \ CDF_2001_S4563131.cc \ CDF_2001_S4751469.cc \
More information about the Rivet-svn mailing list |