Rivet is hosted by Hepforge, IPPP Durham

## Rivet analyses reference

### CMS_2011_S8957746

Event shapes
Experiment: CMS (LHC)
Inspire ID: 886332
Status: VALIDATED
Authors:
• Hendrik Hoeth
References:
Beams: p+ p+
Beam energies: (3500.0, 3500.0) GeV
Run details:
• pp QCD interactions at 7000 GeV. Particles with c*tau>10mm are stable.

Central transverse Thrust and Minor have been measured in proton-proton collisions at $\sqrt{s} = 7$ TeV, with a data sample collected with the CMS detector at the LHC. The sample corresponds to an integrated luminosity of 3.2 inverse picobarns. Input for the variables are anti-$k_t$ jets with $R = 0.5$.

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