Rivet analyses

title: CMSTOTEM_2014_I1294140

Charged particle pseudorapidity distribution at $\sqrt{s}$=8 TeV

Experiment: CMS, TOTEM (LHC)

Inspire ID: 1294140

Status: VALIDATED

Authors: - Panos Katsas

References: - CERN-PH-EP-2014-063 - arXiv: 1405.0722 - Submitted to EPJ C

Beams: p+ p+

Beam energies: (4000.0, 4000.0)GeV

Run details: - Inclusive, NSD-enhanced and SD-enhanced pp events.

The pseudorapidity distribution of charged particles produced in proton-proton collisions at a centre-of-mass energy of 8 TeV are measured in the ranges $|\eta| < 2.2$ and $5.3 < |\eta| < 6.4$, with the CMS and TOTEM detectors, respectively. The measurement is performed with a one-side TOTEM trigger, which is sensitive to 99\% of non-diffractive interactions and diffractive interactions with masses above 3.6 GeV, for three different event selections. An inclusive sample with the least selection bias, a sample enhanced in non-single diffractive events, and a sample enhanced in single-diffractive events were selected.

Source code:CMSTOTEM_2014_I1294140.cc

```c++ // -- C++ --

include "Rivet/Analysis.hh"

include "Rivet/Projections/Beam.hh"

include "Rivet/Projections/ChargedFinalState.hh"

namespace Rivet {

class CMSTOTEM_2014_I1294140 : public Analysis { public:

RIVET_DEFAULT_ANALYSIS_CTOR(CMSTOTEM_2014_I1294140);

void init() {
  ChargedFinalState cfs(Cuts::abseta < 7.0);
  declare(cfs, "CFS");

  book(_Nevt_after_cuts_or, "Nevt_or");
  book(_Nevt_after_cuts_and, "Nevt_and");
  book(_Nevt_after_cuts_xor, "Nevt_xor");
  book(_h_dNch_dEta_OR, 1, 1, 1);
  book(_h_dNch_dEta_AND, 2, 1, 1);
  book(_h_dNch_dEta_XOR, 3, 1, 1);
}


void analyze(const Event& event) {
  // Count forward and backward charged particles
  const ChargedFinalState& charged = apply<ChargedFinalState>(event, "CFS");
  int count_plus = 0, count_minus = 0;
  for (const Particle& p : charged.particles()) {
    if (inRange(p.eta(), 5.3, 6.5)) ++count_plus;
    if (inRange(p.eta(), -6.5, -5.3)) ++count_minus;
  }

  // Cut combinations
  const bool cutsor = (count_plus || count_minus);
  const bool cutsand = (count_plus && count_minus);
  const bool cutsxor = (count_plus && !count_minus) || (!count_plus && count_minus);

  // Increment counters and fill histos
  if (cutsor) _Nevt_after_cuts_or->fill();
  if (cutsand) _Nevt_after_cuts_and->fill();
  if (cutsxor) _Nevt_after_cuts_xor->fill();
  for (const Particle& p : charged.particles()) {
    if (cutsor) _h_dNch_dEta_OR->fill(p.abseta());
    if (cutsand) _h_dNch_dEta_AND->fill(p.abseta());
    if (cutsxor) _h_dNch_dEta_XOR->fill(p.abseta());
  }
}


void finalize() {
  scale(_h_dNch_dEta_OR, 0.5 / *_Nevt_after_cuts_or);
  scale(_h_dNch_dEta_AND, 0.5 / *_Nevt_after_cuts_and);
  scale(_h_dNch_dEta_XOR, 0.5 / *_Nevt_after_cuts_xor);
}

private:

Histo1DPtr _h_dNch_dEta_OR, _h_dNch_dEta_AND, _h_dNch_dEta_XOR;
CounterPtr _Nevt_after_cuts_or, _Nevt_after_cuts_and, _Nevt_after_cuts_xor;

};

RIVET_DECLARE_PLUGIN(CMSTOTEM_2014_I1294140);

} ```