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 dijet mass spectrum
10  class CDF_2008_I805902 : public Analysis {
11  public:
12
13    RIVET_DEFAULT_ANALYSIS_CTOR(CDF_2008_I805902);
14
15
16    /// @name Analysis methods
17    /// @{
18
19    /// Book histograms
20    void init() {
21      FinalState fs;
22      FastJets conefinder(fs, JetAlg::CDFMIDPOINT, 0.7);
23      declare(conefinder, "ConeFinder");
24
25      book(_h_m_dijet ,1, 1, 1);
26    }
27
28
29    /// Do the analysis
30    void analyze(const Event & e) {
31      const JetFinder& jetpro = apply<JetFinder>(e, "ConeFinder");
32      const Jets& jets = jetpro.jetsByPt();
33
34      if (jets.size() < 2) vetoEvent;
35
36      const FourMomentum j0(jets[0].momentum());
37      const FourMomentum j1(jets[1].momentum());
38      if (j1.absrap() > 1.0 || j0.absrap() > 1.0) {
39        vetoEvent;
40      }
41
42      double mjj = FourMomentum(j0+j1).mass();
43      _h_m_dijet->fill(mjj);
44    }
45
46
47    /// Finalize
48    void finalize() {
49      scale(_h_m_dijet, crossSection()/picobarn/sumOfWeights());
50    }
51
52    /// @}
53
54
55  private:
56
57    /// Histogram
58    Histo1DPtr _h_m_dijet;
59
60  };
61
62
63
64  RIVET_DECLARE_ALIASED_PLUGIN(CDF_2008_I805902, CDF_2008_S8093652);
65
66}