|
|
Rivet analyses reference
MC_WWJETS
Monte Carlo validation observables for $W^+[e^+ \, \nu]W^-[\mu^- \, \nu]$ + jets production
Experiment: ()
Status: VALIDATED
Authors:
No references listed
Beams: * *
Beam energies: ANY
Run details:
In addition to the typical jet observables this analysis contains observables related to properties of the WW-pair momentum, correlations between the WW, properties of the W bosons, properties of the leptons, correlations between the opposite charge leptons and correlations with jets.
Source code:
MC_WWJETS.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
163
164
165
166
167 | // -*- C++ -*-
#include "Rivet/Analyses/MC_JetAnalysis.hh"
#include "Rivet/Projections/WFinder.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
namespace Rivet {
/// @brief MC validation analysis for W^+[enu]W^-[munu] + jets events
class MC_WWJETS : public MC_JetAnalysis {
public:
/// Default constructor
MC_WWJETS()
: MC_JetAnalysis("MC_WWJETS", 4, "Jets")
{ }
/// @name Analysis methods
//@{
/// Book histograms
void init() {
FinalState fs;
// set FS cuts from input options
const double etaecut = getOption<double>("ABSETAEMAX", 3.5);
const double ptecut = getOption<double>("PTEMIN", 25.);
Cut cute = Cuts::abseta < etaecut && Cuts::pT > ptecut*GeV;
WFinder wenufinder(fs, cute, PID::ELECTRON, 60.0*GeV, 100.0*GeV, 25.0*GeV, 0.2);
declare(wenufinder, "WenuFinder");
VetoedFinalState wmnuinput;
wmnuinput.addVetoOnThisFinalState(wenufinder);
// set FS cuts from input options
const double etamucut = getOption<double>("ABSETAMUMAX", 3.5);
const double ptmucut = getOption<double>("PTMUMIN", 25.);
Cut cutmu = Cuts::abseta < etamucut && Cuts::pT > ptmucut*GeV;
WFinder wmnufinder(wmnuinput, cutmu, PID::MUON, 60.0*GeV, 100.0*GeV, 25.0*GeV, 0.2);
declare(wmnufinder, "WmnuFinder");
VetoedFinalState jetinput;
jetinput
.addVetoOnThisFinalState(wenufinder)
.addVetoOnThisFinalState(wmnufinder);
// 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(jetinput, clusterAlgo, R);
declare(jetpro, "Jets");
// correlations with jets
book(_h_WW_jet1_deta ,"WW_jet1_deta", 70, -7.0, 7.0);
book(_h_WW_jet1_dR ,"WW_jet1_dR", 25, 1.5, 7.0);
book(_h_We_jet1_dR ,"We_jet1_dR", 25, 0.0, 7.0);
// global stuff
book(_h_HT ,"HT", logspace(100, 100.0, 0.5*(sqrtS()>0.?sqrtS():14000.)));
book(_h_jets_m_12 ,"jets_m_12", logspace(100, 1.0, 0.25*(sqrtS()>0.?sqrtS():14000.)));
MC_JetAnalysis::init();
}
/// Do the analysis
void analyze(const Event& e) {
const double weight = 1.0;
const WFinder& wenufinder = apply<WFinder>(e, "WenuFinder");
if (wenufinder.bosons().size() !=1 ) vetoEvent;
const WFinder& wmnufinder = apply<WFinder>(e, "WmnuFinder");
if (wmnufinder.bosons().size() !=1 ) vetoEvent;
FourMomentum wenu = wenufinder.bosons()[0].momentum();
FourMomentum wmnu = wmnufinder.bosons()[0].momentum();
FourMomentum ww = wenu + wmnu;
// find leptons
FourMomentum ep = wenufinder.constituentLeptons()[0].momentum();
FourMomentum enu = wenufinder.constituentNeutrinos()[0].momentum();
FourMomentum mm = wmnufinder.constituentLeptons()[0].momentum();
FourMomentum mnu = wmnufinder.constituentNeutrinos()[0].momentum();
const Jets& jets = apply<FastJets>(e, "Jets").jetsByPt(_jetptcut);
if (jets.size() > 0) {
_h_WW_jet1_deta->fill(ww.eta()-jets[0].eta(), weight);
_h_WW_jet1_dR->fill(deltaR(ww, jets[0].momentum()), weight);
_h_We_jet1_dR->fill(deltaR(ep, jets[0].momentum()), weight);
}
double HT = ep.pT() + mm.pT() + FourMomentum(enu+mnu).pT();
for (const Jet& jet : jets) HT += jet.pT();
if (HT > 0.0) _h_HT->fill(HT/GeV, weight);
if (jets.size() > 1) {
const FourMomentum jet1 = jets[0].momentum();
const FourMomentum jet2 = jets[1].momentum();
_h_jets_m_12->fill((jet1+jet2).mass()/GeV, weight);
}
MC_JetAnalysis::analyze(e);
}
/// Finalize
void finalize() {
const double norm = crossSection()/picobarn/sumOfWeights();
scale(_h_WW_jet1_deta, norm);
scale(_h_WW_jet1_dR, norm);
scale(_h_We_jet1_dR, norm);
scale(_h_jets_m_12, norm);
scale(_h_HT, norm);
MC_JetAnalysis::finalize();
}
//@}
private:
/// @name Histograms
//@{
Histo1DPtr _h_WW_jet1_deta;
Histo1DPtr _h_WW_jet1_dR;
Histo1DPtr _h_We_jet1_dR;
Histo1DPtr _h_jets_m_12;
Histo1DPtr _h_HT;
//@}
};
// The hook for the plugin system
RIVET_DECLARE_PLUGIN(MC_WWJETS);
}
|
|