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_I723509 : public Analysis {
11  public:
12
13    RIVET_DEFAULT_ANALYSIS_CTOR(STAR_2006_I723509);
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, JetAlg::CDFMIDPOINT, 0.4,
24                             JetMuons::ALL, JetInvisibles::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_I723509, STAR_2006_S6870392);
80
81}