rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

CDF_2007_S7057202

CDF Run II inclusive jet cross-section using the kT algorithm
Experiment: CDF (Tevatron Run 2)
Inspire ID: 743342
Status: VALIDATED
Authors:
  • David Voong
  • Frank Siegert
References:
  • Phys.Rev.D75:092006,2007
  • Erratum-ibid.D75:119901,2007
  • FNAL-PUB 07/026-E
  • hep-ex/0701051
Beams: p- p+
Beam energies: (980.0, 980.0) GeV
Run details:
  • p-pbar collisions at 1960 GeV. Jet pT bins from 54 GeV to 700 GeV. Jet rapidity $< |2.1|$.

CDF Run II measurement of inclusive jet cross sections at a p-pbar collision energy of 1.96 TeV. Jets are reconstructed in the central part of the detector ($|y|<2.1$) using the kT algorithm with an $R$ parameter of 0.7. The minimum jet pT considered is 54 GeV, with a maximum around 700 GeV. The inclusive jet pT is plotted in bins of rapidity $|y|<0.1$, $0.1<|y|<0.7$, $0.7<|y|<1.1$, $1.1<|y|<1.6$ and $1.6<|y|<2.1$.

Source code: CDF_2007_S7057202.cc
 1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Tools/BinnedHistogram.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// @brief CDF inclusive jet cross-section using the \f$ k_\perp \f$ algorithm
10  class CDF_2007_S7057202 : public Analysis {
11  public:
12
13    RIVET_DEFAULT_ANALYSIS_CTOR(CDF_2007_S7057202);
14
15
16    void init() {
17      // Set up projections
18      const FinalState fs;
19      declare(FastJets(fs, FastJets::KT, 0.5), "JetsD05");
20      declare(FastJets(fs, FastJets::KT, 0.7), "JetsD07");
21      declare(FastJets(fs, FastJets::KT, 1.0), "JetsD10");
22
23      // Book histos
24      {Histo1DPtr tmp; _binnedHistosD07.add(  0, 0.1, book(tmp, 1, 1, 1));}
25      {Histo1DPtr tmp; _binnedHistosD07.add(0.1, 0.7, book(tmp, 2, 1, 1));}
26      {Histo1DPtr tmp; _binnedHistosD07.add(0.7, 1.1, book(tmp, 3, 1, 1));}
27      {Histo1DPtr tmp; _binnedHistosD07.add(1.1, 1.6, book(tmp, 4, 1, 1));}
28      {Histo1DPtr tmp; _binnedHistosD07.add(1.6, 2.1, book(tmp, 5, 1, 1));}
29      book(_histoD05 ,6, 1, 1);
30      book(_histoD10 ,7, 1, 1);
31    }
32
33
34    void analyze(const Event& event) {
35      const double weight = 1.0;
36
37      for (const Jet& jet : apply<JetAlg>(event, "JetsD07").jets(Cuts::pT > 54*GeV))
38        _binnedHistosD07.fill(jet.absrap(), jet.pT(), weight);
39
40      for (const Jet& jet : apply<JetAlg>(event, "JetsD05").jets(Cuts::pT > 54*GeV))
41        if (inRange(jet.absrap(), 0.1, 0.7)) _histoD05->fill(jet.pT(), weight);
42
43      for (const Jet& jet : apply<JetAlg>(event, "JetsD10").jets(Cuts::pT > 54*GeV))
44        if (inRange(jet.absrap(), 0.1, 0.7)) _histoD10->fill(jet.pT(), weight);
45    }
46
47
48    // Normalise histograms to cross-section
49    void finalize() {
50      const double xSec = crossSectionPerEvent()/nanobarn;
51
52      scale(_histoD05, xSec);
53      scale(_histoD10, xSec);
54      // scale to xSec/yBinWidth and take into account the double yBinWidth due
55      // to the absolute value of y
56      _binnedHistosD07.scale(xSec/2.0, this);
57    }
58
59
60  private:
61
62    BinnedHistogram _binnedHistosD07;
63
64    /// Single histogram for the \f$R=0.5\f$ \f$k_\perp\f$ jets
65    Histo1DPtr _histoD05;
66
67    /// Single histogram for the \f$R=1.0\f$ \f$k_\perp\f$ jets
68    Histo1DPtr _histoD10;
69
70  };
71
72
73
74  RIVET_DECLARE_ALIASED_PLUGIN(CDF_2007_S7057202, CDF_2007_I743342);
75
76}