1// -*- C++ -*-
 2#include "Rivet/Analyses/MC_KTSPLITTINGS_BASE.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/IdentifiedFinalState.hh"
 5#include "Rivet/Projections/FastJets.hh"
 6
 7namespace Rivet {
 8
 9
10
11
12  /// Generic analysis looking at kt splitting scales of partons
13  class MC_QCD_PARTONS : public MC_KTSPLITTINGS_BASE {
14  public:
15
16    /// Constructor
17    MC_QCD_PARTONS()
18      : MC_KTSPLITTINGS_BASE("MC_QCD_PARTONS", 4, "Jets")
19    {    }
20
21
22  public:
23
24    /// @name Analysis methods
25    /// @{
26
27    /// Book histograms and initialise projections before the run
28    void init() {
29
30      // Projections
31      IdentifiedFinalState partonfs;
32      for (int i=1; i<6; ++i) partonfs.acceptIdPair(i);
33      partonfs.acceptId(PID::GLUON);
34      declare(FastJets(partonfs, JetAlg::KT, 0.6), "Jets");
35
36      MC_KTSPLITTINGS_BASE::init();
37    }
38
39
40
41    /// Perform the per-event analysis
42    void analyze(const Event& event) {
43      MC_KTSPLITTINGS_BASE::analyze(event);
44    }
45
46
47
48    /// Finalize
49    void finalize() {
50      MC_KTSPLITTINGS_BASE::finalize();
51    }
52
53    /// @}
54
55
56  private:
57
58    /// @name Histograms
59    /// @{
60    /// @}
61
62  };
63
64
65  RIVET_DECLARE_PLUGIN(MC_QCD_PARTONS);
66
67}