|
[Rivet-svn] r2026 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Nov 5 22:05:32 GMT 2009
Author: fsiegert Date: Thu Nov 5 22:05:32 2009 New Revision: 2026 Log: Add CDF_1996_S3418421 (not yet tested). Added: trunk/data/anainfo/CDF_1996_S3418421.info trunk/data/plotinfo/CDF_1996_S3418421.plot trunk/data/refdata/CDF_1996_S3418421.aida trunk/src/Analyses/CDF_1996_S3418421.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_1996_S3418421.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_1996_S3418421.info Thu Nov 5 22:05:32 2009 (r2026) @@ -0,0 +1,20 @@ +Name: CDF_1996_S3418421 +Year: 1996 +Summary: Dijet angular distributions +Experiment: CDF +Collider: Tevatron Run 1 +SpiresID: 3418421 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - Phys.Rev.Lett.77:5336-5341,1996 + - doi:10.1103/PhysRevLett.77.5336 + - arXiv:hep-ex/9609011 +RunInfo: + QCD dijet events at Tevatron $\sqrt{s}=1.8 TeV$ without MPI. +NumEvents: 1000000 +PtCuts: [0] +Description: + Measurement of jet angular distributions in events with two jets in the final + state in 5 bins of dijet invariant mass. Based on $106 pb^{-1}$ Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Nov 5 20:54:09 2009 (r2025) +++ trunk/data/anainfo/Makefile.am Thu Nov 5 22:05:32 2009 (r2026) @@ -7,6 +7,7 @@ CDF_1994_S2952106.info \ CDF_1996_S3108457.info \ CDF_1996_S3349578.info \ + CDF_1996_S3418421.info \ CDF_1998_S3618439.info \ CDF_2000_S4155203.info \ CDF_2000_S4266730.info \ Added: trunk/data/plotinfo/CDF_1996_S3418421.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_1996_S3418421.plot Thu Nov 5 22:05:32 2009 (r2026) @@ -0,0 +1,34 @@ +# BEGIN PLOT /CDF_1996_S3418421/d01-x01-y01 +Title=Dijet events with $241 < m_\mathrm{dijet}/\mathrm{GeV} < 300$ +XLabel=$\chi = \exp(|\eta_1-\eta_2|)$ +YLabel=1/N dN/d$\chi$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1996_S3418421/d01-x01-y02 +Title=Dijet events with $300 < m_\mathrm{dijet}/\mathrm{GeV} < 400$ +XLabel=$\chi = \exp(|\eta_1-\eta_2|)$ +YLabel=1/N dN/d$\chi$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1996_S3418421/d01-x01-y03 +Title=Dijet events with $400 < m_\mathrm{dijet}/\mathrm{GeV} < 517$ +XLabel=$\chi = \exp(|\eta_1-\eta_2|)$ +YLabel=1/N dN/d$\chi$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1996_S3418421/d01-x01-y04 +Title=Dijet events with $517 < m_\mathrm{dijet}/\mathrm{GeV} < 625$ +XLabel=$\chi = \exp(|\eta_1-\eta_2|)$ +YLabel=1/N dN/d$\chi$ +LogY=0 +# END PLOT + +# BEGIN PLOT /CDF_1996_S3418421/d01-x01-y05 +Title=Dijet events with $625 < m_\mathrm{dijet}/\mathrm{GeV}$ +XLabel=$\chi = \exp(|\eta_1-\eta_2|)$ +YLabel=1/N dN/d$\chi$ +LogY=0 +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Nov 5 20:54:09 2009 (r2025) +++ trunk/data/plotinfo/Makefile.am Thu Nov 5 22:05:32 2009 (r2026) @@ -7,6 +7,7 @@ CDF_1994_S2952106.plot \ CDF_1996_S3108457.plot \ CDF_1996_S3349578.plot \ + CDF_1996_S3418421.plot \ CDF_1998_S3618439.plot \ CDF_2000_S4155203.plot \ CDF_2000_S4266730.plot \ Added: trunk/data/refdata/CDF_1996_S3418421.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_1996_S3418421.aida Thu Nov 5 22:05:32 2009 (r2026) @@ -0,0 +1,197 @@ +<?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_1996_S3418421" title="(100/N)*D(N)/DCHI" > + <dataPoint> + <measurement value="1.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="31.1" errorPlus="0.7" errorMinus="0.7"/> + </dataPoint> + <dataPoint> + <measurement value="1.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="26.8" errorPlus="0.6" errorMinus="0.6"/> + </dataPoint> + <dataPoint> + <measurement value="2.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.0" errorPlus="0.6" errorMinus="0.6"/> + </dataPoint> + <dataPoint> + <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.4" errorPlus="0.6" errorMinus="0.6"/> + </dataPoint> + <dataPoint> + <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.3" errorPlus="0.7" errorMinus="0.7"/> + </dataPoint> + <dataPoint> + <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="22.5" errorPlus="0.6" errorMinus="0.6"/> + </dataPoint> + <dataPoint> + <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.6" errorPlus="0.8" errorMinus="0.8"/> + </dataPoint> + <dataPoint> + <measurement value="4.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.4" errorPlus="0.8" errorMinus="0.8"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y02" dimension="2" path="/REF/CDF_1996_S3418421" title="(100/N)*D(N)/DCHI" > + <dataPoint> + <measurement value="1.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="31.7" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="1.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="26.5" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="2.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.8" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.2" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.8" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.0" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.3" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="4.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.7" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y03" dimension="2" path="/REF/CDF_1996_S3418421" title="(100/N)*D(N)/DCHI" > + <dataPoint> + <measurement value="1.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="31.9" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="1.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="26.3" errorPlus="0.4" errorMinus="0.4"/> + </dataPoint> + <dataPoint> + <measurement value="2.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.3" errorPlus="0.4" errorMinus="0.4"/> + </dataPoint> + <dataPoint> + <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.9" errorPlus="0.4" errorMinus="0.4"/> + </dataPoint> + <dataPoint> + <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.5" errorPlus="0.4" errorMinus="0.4"/> + </dataPoint> + <dataPoint> + <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.3" errorPlus="0.4" errorMinus="0.4"/> + </dataPoint> + <dataPoint> + <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.6" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + <dataPoint> + <measurement value="4.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.1" errorPlus="0.5" errorMinus="0.5"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y04" dimension="2" path="/REF/CDF_1996_S3418421" title="(100/N)*D(N)/DCHI" > + <dataPoint> + <measurement value="1.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="32.6" errorPlus="1.2" errorMinus="1.2"/> + </dataPoint> + <dataPoint> + <measurement value="1.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="27.2" errorPlus="1.1" errorMinus="1.1"/> + </dataPoint> + <dataPoint> + <measurement value="2.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="25.1" errorPlus="1.1" errorMinus="1.1"/> + </dataPoint> + <dataPoint> + <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="23.0" errorPlus="1.0" errorMinus="1.0"/> + </dataPoint> + <dataPoint> + <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="21.2" errorPlus="1.0" errorMinus="1.0"/> + </dataPoint> + <dataPoint> + <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.1" errorPlus="1.1" errorMinus="1.1"/> + </dataPoint> + <dataPoint> + <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="22.8" errorPlus="1.1" errorMinus="1.1"/> + </dataPoint> + <dataPoint> + <measurement value="4.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.1" errorPlus="1.2" errorMinus="1.2"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y05" dimension="2" path="/REF/CDF_1996_S3418421" title="(100/N)*D(N)/DCHI" > + <dataPoint> + <measurement value="1.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="31.7" errorPlus="2.4" errorMinus="2.4"/> + </dataPoint> + <dataPoint> + <measurement value="1.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="26.5" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="2.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="26.3" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="25.3" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="21.6" errorPlus="2.0" errorMinus="2.0"/> + </dataPoint> + <dataPoint> + <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="22.2" errorPlus="2.1" errorMinus="2.1"/> + </dataPoint> + <dataPoint> + <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="22.4" errorPlus="2.1" errorMinus="2.1"/> + </dataPoint> + <dataPoint> + <measurement value="4.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="24.0" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/CDF_1996_S3418421" title="N(Q=LT)/N(Q=GT)" > + <dataPoint> + <measurement value="263.0" errorPlus="37.0" errorMinus="22.0"/> + <measurement value="0.678" errorPlus="0.021633307652783935" errorMinus="0.021633307652783935"/> + </dataPoint> + <dataPoint> + <measurement value="334.0" errorPlus="66.0" errorMinus="34.0"/> + <measurement value="0.695" errorPlus="0.026925824035672525" errorMinus="0.026925824035672525"/> + </dataPoint> + <dataPoint> + <measurement value="440.0" errorPlus="77.0" errorMinus="40.0"/> + <measurement value="0.703" errorPlus="0.03420526275297414" errorMinus="0.03420526275297414"/> + </dataPoint> + <dataPoint> + <measurement value="557.0" errorPlus="68.0" errorMinus="50.0"/> + <measurement value="0.738" errorPlus="0.05869412236331675" errorMinus="0.05869412236331675"/> + </dataPoint> + <dataPoint> + <measurement value="698.0" errorPlus="1102.0" errorMinus="73.0"/> + <measurement value="0.732" errorPlus="0.11280514172678477" errorMinus="0.11280514172678477"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Nov 5 20:54:09 2009 (r2025) +++ trunk/data/refdata/Makefile.am Thu Nov 5 22:05:32 2009 (r2026) @@ -35,6 +35,7 @@ CDF_1994_S2952106.aida \ CDF_1996_S3108457.aida \ CDF_1996_S3349578.aida \ + CDF_1996_S3418421.aida \ CDF_1998_S3618439.aida \ CDF_2000_S4155203.aida \ CDF_2000_S4266730.aida \ Added: trunk/src/Analyses/CDF_1996_S3418421.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_1996_S3418421.cc Thu Nov 5 22:05:32 2009 (r2026) @@ -0,0 +1,130 @@ +// -*- 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_1996_S3418421 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + CDF_1996_S3418421() + : Analysis("CDF_1996_S3418421"), _nevt(0) + { + setBeams(PROTON, ANTIPROTON); + } + + //@} + + + 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_chi.addHistogram(241.0, 300.0, bookHistogram1D(1, 1, 1)); + _h_chi.addHistogram(300.0, 400.0, bookHistogram1D(1, 1, 2)); + _h_chi.addHistogram(400.0, 517.0, bookHistogram1D(1, 1, 3)); + _h_chi.addHistogram(517.0, 625.0, bookHistogram1D(1, 1, 4)); + _h_chi.addHistogram(625.0, 1800.0, bookHistogram1D(1, 1, 5)); + + _h_ratio = bookHistogram1D(2,1,1); + _chi_above_25.resize(_h_ratio->axis().bins()); + _chi_below_25.resize(_h_ratio->axis().bins()); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + Jets jets = applyProjection<FastJets>(event, "Jets").jetsByEt(50.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 chi = exp(fabs(eta1-eta2)); + if (eta2>2.0 || eta1>2.0 || chi>5.0) { + vetoEvent; + } + + double m = FourMomentum(jet1+jet2).mass(); + _h_chi.fill(m, chi, weight); + + // fill ratio counter + int bin = _h_ratio->coordToIndex(m); + _nevt++; + if (chi>2.5) { + _chi_above_25[bin] += weight; + } + else { + _chi_below_25[bin] += weight; + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + + foreach (AIDA::IHistogram1D* hist, _h_chi.getHistograms()) { + // because HepData contains 100/N instead of 1/N this is 100 in the aida + normalize(hist, 100.0); + } + + // now mimic filling of the ratio histogram + std::vector<double> ratios(_h_ratio->axis().bins()); + for (size_t i=0; i<ratios.size(); ++i) { + ratios[i] = _chi_below_25[i]/_chi_above_25[i]/double(_nevt); + } + for (size_t bin=0; bin<ratios.size(); ++bin) { + double coord = _h_ratio->binMean(bin); + for (size_t n=0; n<_nevt; ++n) { + _h_ratio->fill(coord, ratios[bin]); + } + } + + } + + //@} + + + private: + + // Data members like post-cuts event weight counters go here + std::vector<double> _chi_above_25; + std::vector<double> _chi_below_25; + size_t _nevt; + + private: + + /// @name Histograms + //@{ + BinnedHistogram<double> _h_chi; + AIDA::IHistogram1D* _h_ratio; + //@} + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<CDF_1996_S3418421> plugin_CDF_1996_S3418421; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Nov 5 20:54:09 2009 (r2025) +++ trunk/src/Analyses/Makefile.am Thu Nov 5 22:05:32 2009 (r2026) @@ -16,6 +16,7 @@ CDF_1994_S2952106.cc \ CDF_1996_S3108457.cc \ CDF_1996_S3349578.cc \ + CDF_1996_S3418421.cc \ CDF_1998_S3618439.cc \ CDF_2000_S4155203.cc \ CDF_2000_S4266730.cc \
More information about the Rivet-svn mailing list |