1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FastJets.hh"
 4#include "Rivet/Tools/BinnedHistogram.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// Search for new physics with dijet angular distributions at 13 TeV
10  class CMS_2017_I1519995 : public Analysis {
11  public:
12
13    RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2017_I1519995);
14
15
16    /// Book projections and histograms
17    void init() {
18      FastJets antikt(FinalState(), FastJets::ANTIKT, 0.4);
19      declare(antikt, "ANTIKT");
20      /// @todo We need a better way!!
21      {Histo1DPtr tmp; _h_chi_dijet.add(4800., 8000., book(tmp, 1, 1, 1));}
22      {Histo1DPtr tmp; _h_chi_dijet.add(4200., 4800., book(tmp, 2, 1, 1));}
23      {Histo1DPtr tmp; _h_chi_dijet.add(3600., 4200., book(tmp, 3, 1, 1));}
24      {Histo1DPtr tmp; _h_chi_dijet.add(3000., 3600., book(tmp, 4, 1, 1));}
25      {Histo1DPtr tmp; _h_chi_dijet.add(2400., 3000., book(tmp, 5, 1, 1));}
26      {Histo1DPtr tmp; _h_chi_dijet.add(1900., 2400., book(tmp, 6, 1, 1));}
27    }
28
29
30    /// Per-event analysis
31    void analyze(const Event& event) {
32      const Jets& jets = apply<JetAlg>(event, "ANTIKT").jetsByPt();
33      if (jets.size() < 2) vetoEvent;
34
35      const FourMomentum j0(jets[0].mom()), j1(jets[1].mom());
36      if (fabs(j0.rap()+j1.rap())/2 > 1.11) vetoEvent;
37
38      const double mjj = (j0+j1).mass();
39      const double chi = exp(fabs(j0.rap()-j1.rap()));
40      if (chi < 16) _h_chi_dijet.fill(mjj/GeV, chi, 1.0);
41    }
42
43
44    /// Normalize histograms
45    void finalize() {
46      for (Histo1DPtr hist : _h_chi_dijet.histos()) normalize(hist);
47    }
48
49
50    BinnedHistogram _h_chi_dijet;
51
52  };
53
54
55  // The hook for the plugin system
56  RIVET_DECLARE_PLUGIN(CMS_2017_I1519995);
57
58}