Rivet analyses referenceCMS_2015_I1327224Search for contact interactions \& extra dimensions using dijet angular distributions in proton-proton collisions at $\sqrt{s} = 8$ TeVExperiment: CMS (LHC) Inspire ID: 1327224 Status: VALIDATED Authors:
Beams: p+ p+ Beam energies: (4000.0, 4000.0) GeV Run details:
A search is presented for quark contact interactions and extra spatial dimensions in proton-proton collisions at sqrt(s) = 8 TeV using dijet angular distributions. The search is based on a data set corresponding to an integrated luminosity of 19.7 fb-1 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. Limits on the contact interaction scale from a variety of models at next-to-leading order in QCD corrections are obtained. A benchmark model in which only left-handed quarks participate is excluded up to a scale of 9.0 (11.7) TeV for destructive (constructive) interference at 95% confidence level. Lower limits between 6.0 and 8.4 TeV on the scale of virtual graviton exchange are extracted for the Arkani-Hamed--Dimopoulos--Dvali model of extra spatial dimensions. Source code: CMS_2015_I1327224.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/FastJets.hh"
4
5namespace Rivet {
6
7
8 class CMS_2015_I1327224 : public Analysis {
9 public:
10
11 RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2015_I1327224);
12
13
14 void init() {
15
16 FinalState fs;
17 FastJets antikt(fs, JetAlg::ANTIKT, 0.5);
18 declare(antikt, "ANTIKT");
19
20 book(_h_chi_dijet, {1900., 2400., 3000., 3600., 4200., 8000.});
21 for (auto& b : _h_chi_dijet->bins()) {
22 book(b, _h_chi_dijet->numBins() - b.index() + 1, 1, 1);
23 }
24 }
25
26
27 void analyze(const Event& event) {
28
29 const Jets& jets = apply<JetFinder>(event, "ANTIKT").jetsByPt();
30 if (jets.size() < 2) vetoEvent;
31
32 FourMomentum j0(jets[0].momentum());
33 FourMomentum j1(jets[1].momentum());
34 const double y0 = j0.rapidity();
35 const double y1 = j1.rapidity();
36 if (fabs(y0 + y1) / 2. > 1.11) vetoEvent;
37
38 const double mjj = FourMomentum(j0 + j1).mass();
39 if (mjj/GeV <1900) vetoEvent;
40
41 const double chi = exp(fabs(y0 - y1));
42 if (chi >= 16.) vetoEvent;
43
44 // Fill the histogram
45 _h_chi_dijet->fill(mjj/GeV, chi);
46 }
47
48 void finalize() {
49 normalize(_h_chi_dijet);
50 }
51
52 private:
53
54 Histo1DGroupPtr _h_chi_dijet;
55
56 };
57
58 RIVET_DECLARE_PLUGIN(CMS_2015_I1327224);
59}
|