Rivet analyses referenceCDF_2001_S4563131Inclusive jet cross-sectionExperiment: CDF (Tevatron Run 1) Inspire ID: 552797 Status: VALIDATED Authors:
Beam energies: (900.0, 900.0) GeV Run details:
Measurement of the inclusive jet cross-section for jet transverse energies from 40 to 465 GeV in the pseudorapidity range $0.1<|\eta|<0.7$. The results are based on 87 $\mathrm{pb}^{-1}$ of data. Source code: CDF_2001_S4563131.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 /// @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}
|