Rivet analyses referenceCMS_2017_I1519995Search for new physics with dijet angular distributions in proton-proton collisions at $\sqrt{s} = 13 \text{TeV}$Experiment: CMS (CERN LHC) Inspire ID: 1519995 Status: VALIDATED Authors:
Beam energies: (6500.0, 6500.0) GeV Run details:
Dijet angular distributions are measured in proton-proton collisions at $\sqrt{s}=13 \text{TeV}$, using a data set corresponding to an integrated luminosity of 2.6/fb collected by the CMS detector at the CERN LHC. Dijet angular distributions are found to be in agreement with the perturbative QCD predictions that include electroweak corrections. Source code: CMS_2017_I1519995.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/FastJets.hh"
4
5namespace Rivet {
6
7
8 /// Search for new physics with dijet angular distributions at 13 TeV
9 class CMS_2017_I1519995 : public Analysis {
10 public:
11
12 RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2017_I1519995);
13
14
15 /// Book projections and histograms
16 void init() {
17
18 FastJets antikt(FinalState(), JetAlg::ANTIKT, 0.4);
19 declare(antikt, "ANTIKT");
20
21 book(_h_chi_dijet, {1900., 2400., 3600., 4200., 4800., 8000.});
22 for (auto& b : _h_chi_dijet->bins()) {
23 book(b, _h_chi_dijet->numBins() - b.index() + 1, 1, 1);
24 }
25 }
26
27
28 /// Per-event analysis
29 void analyze(const Event& event) {
30 const Jets& jets = apply<JetFinder>(event, "ANTIKT").jetsByPt();
31 if (jets.size() < 2) vetoEvent;
32
33 const FourMomentum j0(jets[0].mom()), j1(jets[1].mom());
34 if (fabs(j0.rap()+j1.rap())/2 > 1.11) vetoEvent;
35
36 const double mjj = (j0+j1).mass();
37 const double chi = exp(fabs(j0.rap()-j1.rap()));
38 if (chi < 16) _h_chi_dijet->fill(mjj/GeV, chi);
39 }
40
41
42 /// Normalize histograms
43 void finalize() {
44 normalize(_h_chi_dijet);
45 }
46
47
48 Histo1DGroupPtr _h_chi_dijet;
49
50 };
51
52
53 RIVET_DECLARE_PLUGIN(CMS_2017_I1519995);
54
55}
|