Rivet Analyses Reference
MC_JETS
Monte Carlo validation observables for jet production
Experiment: ()
Status: VALIDATED
Authors:No references listed
Beams: * *
Beam energies: ANY
Run details:- Pure QCD jet production events at an arbitrary collider.
Jets with $p_\perp>20$ GeV are constructed with an anti-$k_\perp$ jet finder with $R=0.4$ and projected onto many different observables.
Source code:
MC_JETS.cc 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
| // -*- C++ -*-
#include "Rivet/Analyses/MC_JetAnalysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FastJets.hh"
#include "fastjet/contrib/SoftDrop.hh"
namespace Rivet {
/// @brief MC validation analysis for jet events
class MC_JETS : public MC_JetAnalysis {
public:
MC_JETS()
: MC_JetAnalysis("MC_JETS", 4, "Jets")
{ }
void init() {
FinalState fs;
FastJets jetpro(fs, FastJets::ANTIKT, 0.4);
const string groomopt = getOption("GROOM", "");
if (groomopt == "SD") {
jetpro.addTrf(new fastjet::contrib::SoftDrop(0.0, 0.1));
} else if (groomopt == "TRIM") {
jetpro.addTrf(new fastjet::Filter(fastjet::JetDefinition(fastjet::kt_algorithm, 0.2), fastjet::SelectorPtFractionMin(0.05)));
} else if (groomopt != "") {
MSG_WARNING("Unknown GROOM=" + groomopt + " option. Not applying jet grooming");
}
declare(jetpro, "Jets");
MC_JetAnalysis::init();
}
void analyze(const Event& event) {
MC_JetAnalysis::analyze(event);
}
void finalize() {
MC_JetAnalysis::finalize();
}
};
// The hook for the plugin system
RIVET_DECLARE_PLUGIN(MC_JETS);
}
|