Rivet Analyses Reference
CDF_2008_S8093652
Dijet mass spectrum
Experiment: CDF (Tevatron Run 2)
Inspire ID: 805902
Status: VALIDATED
Authors:References:Beams: p- p+
Beam energies: (980.0, 980.0) GeV
Run details:- $p \bar{p} \to$ jets at 1960 GeV
Dijet mass spectrum from 0.2 TeV to 1.4 TeV in $p \bar{p}$ collisions at $\sqrt{s} = 1.96$ TeV, based on an integrated luminosity of 1.13 fb$^{-1}$.
Source code:
CDF_2008_S8093652.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
| // -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FastJets.hh"
namespace Rivet {
/// @brief CDF dijet mass spectrum
class CDF_2008_S8093652 : public Analysis {
public:
RIVET_DEFAULT_ANALYSIS_CTOR(CDF_2008_S8093652);
/// @name Analysis methods
//@{
/// Book histograms
void init() {
FinalState fs;
FastJets conefinder(fs, FastJets::CDFMIDPOINT, 0.7);
declare(conefinder, "ConeFinder");
book(_h_m_dijet ,1, 1, 1);
}
/// Do the analysis
void analyze(const Event & e) {
const JetAlg& jetpro = apply<JetAlg>(e, "ConeFinder");
const Jets& jets = jetpro.jetsByPt();
if (jets.size() < 2) vetoEvent;
const FourMomentum j0(jets[0].momentum());
const FourMomentum j1(jets[1].momentum());
if (j1.absrap() > 1.0 || j0.absrap() > 1.0) {
vetoEvent;
}
double mjj = FourMomentum(j0+j1).mass();
_h_m_dijet->fill(mjj);
}
/// Finalize
void finalize() {
scale(_h_m_dijet, crossSection()/sumOfWeights());
}
//@}
private:
/// Histogram
Histo1DPtr _h_m_dijet;
};
RIVET_DECLARE_ALIASED_PLUGIN(CDF_2008_S8093652, CDF_2008_I805902);
}
|