rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

CMS_2013_I1265659

Probing color coherence effects in $pp$ collisions at $\sqrt{s} = 7$ TeV
Experiment: CMS (LHC)
Inspire ID: 1265659
Status: VALIDATED
Authors:
• Maxime Gouzevitch
• Chawon Park
• Inkyu Park
References:
• CMS-SMP-12-010
• CERN-PH-EP-2013-200
• arXiv: 1311.5815
Beams: p+ p+
Beam energies: (3500.0, 3500.0) GeV
Run details:
• $pp$ QCD interactions at $\sqrt{s} = 7$ TeV. Data collected by CMS during the year 2010.

A study of color coherence effects in $pp$ collisions at a center-of-mass energy of 7 TeV is presented. The data used in the analysis were collected in 2010 with the CMS detector at the LHC and correspond to an integrated luminosity of 36/pb. Events are selected that contain at least three jets and where the two jets with the largest transverse momentum exhibit a back-to-back topology. The measured angular correlation between the second- and third-leading jet is shown to be sensitive to color coherence effects, and is compared to the predictions of Monte Carlo models with various implementations of color coherence. None of the models describe the data satisfactorily.

Source code: CMS_2013_I1265659.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 // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/FastJets.hh" namespace Rivet { class CMS_2013_I1265659 : public Analysis { public: /// Constructor CMS_2013_I1265659() : Analysis("CMS_2013_I1265659") { } /// Book histograms and initialise projections before the run void init() { const FastJets jets(FinalState((Cuts::etaIn(-10, 10))), FastJets::ANTIKT, 0.5); declare(jets, "Jets"); book(_h_hTotD ,1, 1, 1); book(_h_hTotDF ,1, 1, 2); } /// Perform the per-event analysis void analyze(const Event& event) { const Jets& jets = apply(event, "Jets").jetsByPt(30.0*GeV); if (jets.size() < 3) vetoEvent; const FourMomentum jet1 = jets[0].momentum(); const FourMomentum jet2 = jets[1].momentum(); const FourMomentum jet3 = jets[2].momentum(); // Cut on lead jet pT and lead/sublead jet centrality if (jet1.pT() < 100*GeV) vetoEvent; if (jet1.abseta() > 2.5 || jet2.abseta() > 2.5) vetoEvent; // Construct eta & phi distances between 2nd and 3rd jets double dEta23 = jet3.eta() - jet2.eta(); ///< Note not abs double dPhi23 = jet3.phi() - jet2.phi(); ///< Note not abs if (dPhi23 > M_PI) dPhi23 -= 2*M_PI; ///< @todo Use mapTo... functions? if (dPhi23 < -M_PI) dPhi23 += 2*M_PI; ///< @todo Use mapTo... functions? // Cut on distance between 2nd and 3rd jets const double R23 = add_quad(dPhi23, dEta23); if (!inRange(R23, 0.5, 1.5)) vetoEvent; // Cut on dijet mass const FourMomentum diJet = jet1 + jet2; if (diJet.mass() < 220*GeV) vetoEvent; // Calc beta and fill histogram (choose central or fwd histo inline) double beta = fabs(atan2(dPhi23, sign(jet2.eta())*dEta23)); ((jet2.abseta() < 0.8) ? _h_hTotD : _h_hTotDF)->fill(beta, 1.0); } /// Normalise histograms etc., after the run void finalize() { const double width = _h_hTotD->bin(0).xWidth(); normalize(_h_hTotD, width); normalize(_h_hTotDF, width); } private: /// @name Histograms Histo1DPtr _h_hTotD; Histo1DPtr _h_hTotDF; //@} }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(CMS_2013_I1265659); }