1// -*- C++ -*-
 2#include "Rivet/Analyses/MC_KTSPLITTINGS_BASE.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5
 6namespace Rivet {
 7
 8
 9
10
11  /// @brief MC validation analysis for jet events
12  class MC_KTSPLITTINGS : public MC_KTSPLITTINGS_BASE {
13  public:
14
15    MC_KTSPLITTINGS()
16      : MC_KTSPLITTINGS_BASE("MC_KTSPLITTINGS", 4, "Jets")
17    {    }
18
19
20  public:
21
22    void init() {
23      // set clustering radius from input option
24      const double R = getOption<double>("R", 0.6);
25
26      FastJets jetpro(FinalState(), JetAlg::KT, R);
27      declare(jetpro, "Jets");
28      MC_KTSPLITTINGS_BASE::init();
29    }
30
31
32    void analyze(const Event& event) {
33      MC_KTSPLITTINGS_BASE::analyze(event);
34    }
35
36
37    void finalize() {
38      MC_KTSPLITTINGS_BASE::finalize();
39    }
40
41  };
42
43
44  RIVET_DECLARE_PLUGIN(MC_KTSPLITTINGS);
45
46}