|
|
Rivet analyses reference
MC_HHJETS
Monte Carlo validation observables for $HH$ (stable) + jets production
Experiment: ()
Status: VALIDATED
Authors:
No references listed
Beams: * *
Beam energies: ANY
Run details:
- $HH$ production with stable Higgses
The available observables are the Higgs boson pair invariant mass, the separation between the two Higgs bosons, pT of the di-Higgs (any, hardest, second hardest), pT and pseudorapidity of Higgs bosons (any, hardest, second hardest), pT of jets 1--4, jet multiplicity, $\Delta\eta(h, \text{jet1})$, $\Delta R(\text{jet2}, \text{jet3})$, differential jet rates 0->1, 1->2, 2->3, 3->4, and integrated 0--4 jet rates.
Source code:
MC_HHJETS.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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162 | // -*- C++ -*-
#include "Rivet/Analyses/MC_JetAnalysis.hh"
#include "Rivet/Projections/ZFinder.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
#include "Rivet/Projections/IdentifiedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"
namespace Rivet {
/// @brief MC validation analysis for higgs pairs events (stable Higgses)
class MC_HHJETS : public MC_JetAnalysis {
public:
/// Default constructor
MC_HHJETS()
: MC_JetAnalysis("MC_HHJETS", 4, "Jets")
{ }
/// @name Analysis methods
//@{
/// Book histograms
void init() {
IdentifiedFinalState ifs(Cuts::abseta < 10.0 && Cuts::pT > 0*GeV);
ifs.acceptId(25);
declare(ifs,"IFS");
VetoedFinalState vfs;
vfs.addVetoPairId(25);
// set ptcut from input option
const double jetptcut = getOption<double>("PTJMIN", 20.0);
_jetptcut = jetptcut * GeV;
// set clustering radius from input option
const double R = getOption<double>("R", 0.4);
// set clustering algorithm from input option
FastJets::Algo clusterAlgo;
const string algoopt = getOption("ALGO", "ANTIKT");
if ( algoopt == "KT" ) {
clusterAlgo = FastJets::KT;
} else if ( algoopt == "CA" ) {
clusterAlgo = FastJets::CA;
} else if ( algoopt == "ANTIKT" ) {
clusterAlgo = FastJets::ANTIKT;
} else {
MSG_WARNING("Unknown jet clustering algorithm option " + algoopt + ". "
"Defaulting to anti-kT");
clusterAlgo = FastJets::ANTIKT;
}
FastJets jetpro(vfs, clusterAlgo, R);
declare(jetpro, "Jets");
book(_h_HH_mass ,"HH_mass", 250, 240, 4000.0);
book(_h_HH_dR ,"HH_dR", 25, 0.5, 10.0);
book(_h_HH_dPhi ,"HH_dPhi", 64, 0, 3.2);
book(_h_HH_deta,"HH_deta", 50, -5, 5);
book(_h_H_pT ,"H_pT", 50, 0, 2000.0);
book(_h_HH_pT ,"HH_pT", 200, 0, 2000.0);
book(_h_H_pT1 ,"H_pT1", 200, 0, 2000.0);
book(_h_H_pT2 ,"H_pT2", 200, 0, 2000.0);
book(_h_H_eta ,"H_eta", 50, -5.0, 5.0);
book(_h_H_eta1 ,"H_eta1", 50, -5.0, 5.0);
book(_h_H_eta2 ,"H_eta2", 50, -5.0, 5.0);
book(_h_H_phi ,"H_phi", 25, 0.0, TWOPI);
book(_h_H_jet1_deta ,"H_jet1_deta", 50, -5.0, 5.0);
book(_h_H_jet1_dR ,"H_jet1_dR", 25, 0.5, 7.0);
MC_JetAnalysis::init();
}
/// Do the analysis
void analyze(const Event & e) {
const IdentifiedFinalState& ifs = apply<IdentifiedFinalState>(e, "IFS");
Particles allp = ifs.particlesByPt();
if (allp.empty()) vetoEvent;
const double weight = 1.0;
FourMomentum hmom = allp[0].momentum();
if (allp.size() > 1) {
FourMomentum hmom2(allp[1].momentum());
_h_HH_dR->fill(deltaR(hmom, hmom2), weight);
_h_HH_dPhi->fill(deltaPhi(hmom, hmom2), weight);
_h_HH_deta->fill(hmom.eta()-hmom2.eta(), weight);
_h_HH_pT->fill((hmom+hmom2).pT(), weight);
_h_HH_mass->fill((hmom+hmom2).mass(), weight);
if (hmom.pT() > hmom2.pT()) {
_h_H_pT1->fill(hmom.pT(), weight);
_h_H_eta1->fill(hmom.eta(), weight);
_h_H_pT2->fill(hmom2.pT(), weight);
_h_H_eta2->fill(hmom2.eta(), weight);
} else {
_h_H_pT1->fill(hmom2.pT(), weight);
_h_H_eta1->fill(hmom2.eta(), weight);
_h_H_pT2->fill(hmom.pT(), weight);
_h_H_eta2->fill(hmom.eta(), weight);
}
}
_h_H_pT->fill(hmom.pT(), weight);
_h_H_eta->fill(hmom.eta(), weight);
_h_H_phi->fill(hmom.azimuthalAngle(), weight);
// Get the jet candidates
Jets jets = apply<FastJets>(e, "Jets").jetsByPt(20.0*GeV);
if (!jets.empty()) {
_h_H_jet1_deta->fill(deltaEta(hmom, jets[0]), weight);
_h_H_jet1_dR->fill(deltaR(hmom, jets[0]), weight);
}
MC_JetAnalysis::analyze(e);
}
/// Finalize
void finalize() {
scale(_h_HH_mass, crossSection()/sumOfWeights());
scale(_h_HH_dR, crossSection()/sumOfWeights());
scale(_h_HH_deta, crossSection()/sumOfWeights());
scale(_h_HH_dPhi, crossSection()/sumOfWeights());
scale(_h_H_pT, crossSection()/sumOfWeights());
scale(_h_H_pT1, crossSection()/sumOfWeights());
scale(_h_H_pT2, crossSection()/sumOfWeights());
scale(_h_HH_pT, crossSection()/sumOfWeights());
scale(_h_H_eta, crossSection()/sumOfWeights());
scale(_h_H_eta1, crossSection()/sumOfWeights());
scale(_h_H_eta2, crossSection()/sumOfWeights());
scale(_h_H_phi, crossSection()/sumOfWeights());
scale(_h_H_jet1_deta, crossSection()/sumOfWeights());
scale(_h_H_jet1_dR, crossSection()/sumOfWeights());
MC_JetAnalysis::finalize();
}
//@}
private:
/// @name Histograms
//@{
Histo1DPtr _h_HH_mass, _h_HH_pT, _h_HH_dR, _h_HH_deta, _h_HH_dPhi;
Histo1DPtr _h_H_pT, _h_H_pT1, _h_H_pT2, _h_H_eta, _h_H_eta1, _h_H_eta2, _h_H_phi;
Histo1DPtr _h_H_jet1_deta, _h_H_jet1_dR;
//@}
};
// The hook for the plugin system
RIVET_DECLARE_PLUGIN(MC_HHJETS);
}
|
|