rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

STAR_2006_S6870392

Inclusive jet cross-section in pp at 200 GeV
Experiment: STAR (RHIC pp 200 GeV)
Inspire ID: 723509
Status: VALIDATED
Authors:
  • Hendrik Hoeth
References:
  • Phys. Rev. Lett. 97, 252001
  • hep-ex/0608030
Beams: p+ p+
Beam energies: (100.0, 100.0) GeV
Run details:
  • pp at 200 GeV

Inclusive jet cross section as a function of pT in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC.

Source code: STAR_2006_S6870392.cc
 1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// STAR inclusive jet cross-section in pp at 200 GeV
10  class STAR_2006_S6870392 : public Analysis {
11  public:
12
13    RIVET_DEFAULT_ANALYSIS_CTOR(STAR_2006_S6870392);
14
15
16    /// @name Analysis methods
17    /// @{
18
19    /// Book projections and histograms
20    void init() {
21      FinalState fs((Cuts::etaIn(-2.0, 2.0)));
22      declare(fs, "FS");
23      declare(FastJets(fs, FastJets::CDFMIDPOINT, 0.4,
24                             JetAlg::Muons::ALL, JetAlg::Invisibles::NONE,
25                             nullptr, 0.5), "MidpointJets");
26
27      book(_h_jet_pT_MB ,1, 1, 1);
28      book(_h_jet_pT_HT ,2, 1, 1);
29    }
30
31
32    /// Do the analysis
33    void analyze(const Event& event) {
34      // Skip if the event is empty
35      const FinalState& fs = apply<FinalState>(event, "FS");
36      if (fs.empty()) {
37        MSG_DEBUG("Skipping event " << numEvents() << " because no final state found ");
38        vetoEvent;
39      }
40
41      // Find jets
42      const FastJets& jetpro = apply<FastJets>(event, "MidpointJets");
43      const Jets& jets = jetpro.jetsByPt();
44      if (!jets.empty()) {
45        const Jet& j1 = jets.front();
46        if (inRange(fabs(j1.eta()), 0.2, 0.8)) {
47          for (const Jet& j : jets) {
48            const FourMomentum pj = j.momentum();
49            _h_jet_pT_MB->fill(pj.pT());
50            _h_jet_pT_HT->fill(pj.pT());
51          }
52        }
53      }
54    }
55
56
57    /// Finalize
58    void finalize() {
59      double normalisation = crossSection()/picobarn/sumOfWeights()/(2*0.6*2*M_PI);
60      scale(_h_jet_pT_MB, normalisation);
61      scale(_h_jet_pT_HT, normalisation);
62    }
63
64    /// @}
65
66
67  private:
68
69    /// @name Histograms
70    /// @{
71    Histo1DPtr _h_jet_pT_MB;
72    Histo1DPtr _h_jet_pT_HT;
73    /// @}
74
75  };
76
77
78
79  RIVET_DECLARE_ALIASED_PLUGIN(STAR_2006_S6870392, STAR_2006_I723509);
80
81}