|
Rivet analyses reference
ATLAS_2011_S8971293
Dijet azimuthal decorrelations
Experiment: ATLAS (LHC)
Inspire ID: 889546
Status: VALIDATED
Authors:
References:
Beams: p+ p+
Beam energies: (3500.0, 3500.0) GeV
Run details:
- pp QCD interactions at 7000 GeV. The distributions are binned in leading pT starting at 110 GeV with the last bin starting at 800 GeV.
Dijet azimuthal decorrelation measured by ATLAS at 7 TeV. Jets are anti-$k_t$ with $R = 0.6$, $p_\perp > 100$ GeV, $|\eta| < 0.8$. The analysis is binned in leading jet $p_\perp$ bins. All data is fully corrected for detector effects.
Source code:
ATLAS_2011_S8971293.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 | // -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/BinnedHistogram.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FastJets.hh"
namespace Rivet {
class ATLAS_2011_S8971293 : public Analysis {
public:
/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(ATLAS_2011_S8971293);
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
/// Initialise and register projections
declare(FastJets(FinalState(), FastJets::ANTIKT, 0.6), "AntiKtJets06");
/// Book histograms
/// @todo Provide a nicer way!
{Histo1DPtr tmp; _h_deltaPhi.add(110., 160., book(tmp, 1, 1, 1));}
{Histo1DPtr tmp; _h_deltaPhi.add(160., 210., book(tmp, 1, 1, 2));}
{Histo1DPtr tmp; _h_deltaPhi.add(210., 260., book(tmp, 1, 1, 3));}
{Histo1DPtr tmp; _h_deltaPhi.add(260., 310., book(tmp, 1, 1, 4));}
{Histo1DPtr tmp; _h_deltaPhi.add(310., 400., book(tmp, 1, 1, 5));}
{Histo1DPtr tmp; _h_deltaPhi.add(400., 500., book(tmp, 1, 1, 6));}
{Histo1DPtr tmp; _h_deltaPhi.add(500., 600., book(tmp, 1, 1, 7));}
{Histo1DPtr tmp; _h_deltaPhi.add(600., 800., book(tmp, 1, 1, 8));}
{Histo1DPtr tmp; _h_deltaPhi.add(800.,10000.,book(tmp, 1, 1, 9));}
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = 1.0;
Jets jets06;
for (const Jet& jet : apply<FastJets>(event, "AntiKtJets06").jetsByPt(100.0*GeV)) {
if (jet.absrap() < 2.8) {
jets06.push_back(jet);
}
}
if (jets06.size()>1){
if (fabs(jets06[0].rapidity())<0.8 && fabs(jets06[1].rapidity())<0.8) {
double observable = mapAngle0ToPi(jets06[0].phi()-jets06[1].phi()) / M_PI;
_h_deltaPhi.fill(jets06[0].pT(), observable, weight);
}
}
}
/// Normalise histograms etc., after the run
void finalize() {
for (Histo1DPtr hist : _h_deltaPhi.histos()) {
normalize(hist, 1/M_PI);
}
}
/// @}
private:
/// Histograms
BinnedHistogram _h_deltaPhi;
};
RIVET_DECLARE_ALIASED_PLUGIN(ATLAS_2011_S8971293, ATLAS_2011_I889546);
}
|
|