Rivet analyses referenceCMS_2012_I1090423Dijet angular distributions in $pp$ collisions at $\sqrt{s} = 7$ TeVExperiment: CMS (LHC) Inspire ID: 1090423 Status: VALIDATED Authors:
Beam energies: (3500.0, 3500.0) GeV Run details:
Measurement of dijet angular distributions in proton-proton collisions at a center-of-mass energy of 7 TeV. The data sample has a total integrated luminosity of 2.2 inverse femtobarns, recorded by the CMS experiment at the LHC. Normalized dijet angular distributions have been measured for dijet invariant masses from 0.4 TeV to above 3 TeV. Source code: CMS_2012_I1090423.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/FastJets.hh"
4
5namespace Rivet {
6
7
8 class CMS_2012_I1090423 : public Analysis {
9 public:
10
11 RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2012_I1090423);
12
13 void init() {
14
15 FinalState fs;
16 FastJets antikt(fs, JetAlg::ANTIKT, 0.5);
17 declare(antikt, "ANTIKT");
18
19 book(_h_chi_dijet, {400., 600., 800., 1000., 1200., 1500., 1900., 2400., 3000., 7000.});
20 for (auto& b : _h_chi_dijet->bins()) {
21 book(b, _h_chi_dijet->numBins() - b.index() + 1, 1, 1);
22 }
23 }
24
25
26 void analyze(const Event& event) {
27 const Jets& jets = apply<JetFinder>(event, "ANTIKT").jetsByPt();
28 if (jets.size() < 2) vetoEvent;
29
30 const double y0 = jets[0].rapidity();
31 const double y1 = jets[1].rapidity();
32 if (fabs(y0+y1)/2 > 1.11) vetoEvent;
33
34 const double chi = exp(fabs(y0-y1));
35 if (chi > 16) vetoEvent;
36
37 const FourMomentum jj = jets[0].momentum() + jets[1].momentum();
38 _h_chi_dijet->fill(jj.mass(), chi);
39 }
40
41
42 void finalize() {
43 normalize(_h_chi_dijet);
44 }
45
46
47 private:
48
49 Histo1DGroupPtr _h_chi_dijet;
50
51 };
52
53
54
55 RIVET_DECLARE_PLUGIN(CMS_2012_I1090423);
56
57}
|