Rivet analyses

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 |η| < 2.2 and 5.3 < |η| < 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++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/Beam.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);

}

Rivet analyses

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

Source code:CMSTOTEM_2014_I1294140.cc

Source code:`CMSTOTEM_2014_I1294140.cc`