1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5#include "Rivet/Projections/Thrust.hh"
 6
 7namespace Rivet {
 8
 9
10  /// Event shapes at 7 TeV
11  class CMS_2011_I886332 : public Analysis {
12  public:
13
14    RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2011_I886332);
15
16
17    /// Initialization, called once before running
18    void init() {
19      // Projections
20      const FastJets jets(FinalState((Cuts::etaIn(-5.0, 5.0))), JetAlg::ANTIKT, 0.5);
21      declare(jets, "Jets");
22
23      // Book histograms
24      book(_hist_T_90  ,1, 1, 1);
25      book(_hist_m_90  ,2, 1, 1);
26      book(_hist_T_125 ,3, 1, 1);
27      book(_hist_m_125 ,4, 1, 1);
28      book(_hist_T_200 ,5, 1, 1);
29      book(_hist_m_200 ,6, 1, 1);
30    }
31
32
33    void analyze(const Event& event) {
34      const Jets& jets = apply<FastJets>(event, "Jets").jetsByPt(Cuts::pT > 30.0*GeV);
35      if (jets.size() < 2 ||
36          fabs(jets[0].eta()) >= 1.3 ||
37          fabs(jets[1].eta()) >= 1.3 ||
38          jets[0].pT() < 90*GeV) {
39        vetoEvent;
40      }
41      std::vector<Vector3> momenta;
42      for (const Jet& j : jets) {
43        if (j.abseta() < 1.3) {
44          Vector3 mom = j.p3();
45          mom.setZ(0.0);
46          momenta.push_back(mom);
47        }
48      }
49      if (momenta.size() == 2) {
50        // We need to use a ghost so that Thrust.calc() doesn't return 1.
51        momenta.push_back(Vector3(1e-10*MeV, 0., 0.));
52      }
53      Thrust thrust;
54      thrust.calc(momenta);
55
56      // The lowest bin also includes the underflow:
57      const double T = max(log(1-thrust.thrust()), -12.0);
58      const double M = max(log(thrust.thrustMajor()), -6.0);
59      if (jets[0].pT()/GeV > 200) {
60        _hist_T_200->fill(T);
61        _hist_m_200->fill(M);
62      } else if (jets[0].pT()/GeV > 125) {
63        _hist_T_125->fill(T);
64        _hist_m_125->fill(M);
65      } else if (jets[0].pT()/GeV > 90) {
66        _hist_T_90->fill(T);
67        _hist_m_90->fill(M);
68      }
69    }
70
71
72    void finalize() {
73      normalize(_hist_T_90);
74      normalize(_hist_m_90);
75      normalize(_hist_T_125);
76      normalize(_hist_m_125);
77      normalize(_hist_T_200);
78      normalize(_hist_m_200);
79    }
80
81
82  private:
83
84    /// @{
85    Histo1DPtr _hist_T_90;
86    Histo1DPtr _hist_m_90;
87    Histo1DPtr _hist_T_125;
88    Histo1DPtr _hist_m_125;
89    Histo1DPtr _hist_T_200;
90    Histo1DPtr _hist_m_200;
91    /// @}
92
93  };
94
95
96
97  RIVET_DECLARE_ALIASED_PLUGIN(CMS_2011_I886332, CMS_2011_S8957746);
98
99}