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  /// @brief CDF Run I inclusive jet cross-section
10  class CDF_2001_S4563131 : public Analysis {
11  public:
12
13    RIVET_DEFAULT_ANALYSIS_CTOR(CDF_2001_S4563131);
14
15
16    /// @name Analysis methods
17    //@{
18
19    /// Book histograms and initialise projections before the run
20    void init() {
21      FinalState fs(Cuts::abseta < 4.2);
22      declare(FastJets(fs, FastJets::CDFJETCLU, 0.7), "Jets");
23      book(_h_ET ,1, 1, 1);
24    }
25
26
27    /// Perform the per-event analysis
28    void analyze(const Event& event) {
29      Jets jets = apply<FastJets>(event, "Jets").jets(Cuts::Et > 40*GeV && Cuts::abseta >= 0.1 && Cuts::abseta <= 0.7, cmpMomByEt);
30      for (const Jet& jet : jets) {
31        _h_ET->fill(jet.Et());
32      }
33    }
34
35
36    /// Normalise histograms etc., after the run
37    void finalize() {
38      const double deta = 1.2;
39      scale(_h_ET, crossSection()/sumOfWeights()/deta/nanobarn);
40    }
41
42    //@}
43
44
45  private:
46
47    /// Histogram
48    Histo1DPtr _h_ET;
49
50  };
51
52
53
54  RIVET_DECLARE_ALIASED_PLUGIN(CDF_2001_S4563131, CDF_2001_I552797);
55
56}