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      FJPluginPtr cdfMidPointPlugin = FastJets::mkPlugin<JetAlg::CDFMIDPOINT>(0.4, 0.5, 0.5);
24      declare(FastJets(fs, cdfMidPointPlugin), "MidpointJets");
25
26      book(_h_jet_pT_MB ,1, 1, 1);
27      book(_h_jet_pT_HT ,2, 1, 1);
28    }
29
30
31    /// Do the analysis
32    void analyze(const Event& event) {
33      // Skip if the event is empty
34      const FinalState& fs = apply<FinalState>(event, "FS");
35      if (fs.empty()) {
36        MSG_DEBUG("Skipping event " << numEvents() << " because no final state found ");
37        vetoEvent;
38      }
39
40      // Find jets
41      const FastJets& jetpro = apply<FastJets>(event, "MidpointJets");
42      const Jets& jets = jetpro.jetsByPt();
43      if (!jets.empty()) {
44        const Jet& j1 = jets.front();
45        if (inRange(fabs(j1.eta()), 0.2, 0.8)) {
46          for (const Jet& j : jets) {
47            const FourMomentum pj = j.momentum();
48            _h_jet_pT_MB->fill(pj.pT());
49            _h_jet_pT_HT->fill(pj.pT());
50          }
51        }
52      }
53    }
54
55
56    /// Finalize
57    void finalize() {
58      double normalisation = crossSection()/picobarn/sumOfWeights()/(2*0.6*2*M_PI);
59      scale(_h_jet_pT_MB, normalisation);
60      scale(_h_jet_pT_HT, normalisation);
61    }
62
63    /// @}
64
65
66  private:
67
68    /// @name Histograms
69    /// @{
70    Histo1DPtr _h_jet_pT_MB;
71    Histo1DPtr _h_jet_pT_HT;
72    /// @}
73
74  };
75
76
77
78  RIVET_DECLARE_ALIASED_PLUGIN(STAR_2006_I723509, STAR_2006_S6870392);
79
80}