Rivet analyses

title: CMS_2019_I1672011

Prompt $\psi(2S)$ production at 5.02 TeV

Experiment: CMS (LHC)

Inspire ID: 1672011

Status: VALIDATED

Authors: - Peter Richardson

References: - Phys.Lett.B 790 (2019) 509-532, 2019 - arXiv: 1805.02248 - Expt page: CMS-HIN-16-015

Beams: p+ p+

Beam energies: (2510.0, 2510.0)GeV

Run details: - hadronic events with psi(2S) production

Measurement of the double differential cross section for prompt $\psi(2S)$ production at 5.02 TeV by the CMS collaboration. Only the proton-proton results are implemented.

Source code:CMS_2019_I1672011.cc

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

include "Rivet/Analysis.hh"

include "Rivet/Projections/UnstableParticles.hh"

namespace Rivet {

/// @brief psi(2S) at 5.02 class CMS_2019_I1672011 : public Analysis { public:

/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2019_I1672011);


/// @name Analysis methods
/// @{

/// Book histograms and initialise projections before the run
void init() {
  declare(UnstableParticles(), "UFS");
  // histograms
  book(_h_psi2S, {0.0, 0.9, 1.5, 1.93, 2.4});
  for (int iy = 0; iy < 4; ++iy) {
    book(_h_psi2S->bin(iy + 1), 3, 1, 4 - iy);
  }
}


/// Perform the per-event analysis
void analyze(const Event& event) {

  // Final state of unstable particles to get particle spectra
  const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");

  for (const Particle& p : ufs.particles(Cuts::pid == 100443)) {
    // prompt
    if (p.fromBottom()) continue;
    const double absrap = p.absrap();
    const double xp = p.perp();
    _h_psi2S->fill(absrap, xp);
  }
}


/// Normalise histograms etc., after the run
void finalize() {
  // br to muons PDG 2021 (e+e- due large errors on mu+mu-)
  const double br = 0.00793;
  // factor of 1/2 as +/- y
  const double factor = 0.5 * br * crossSection() / nanobarn / sumOfWeights();
  scale(_h_psi2S, factor);
  divByGroupWidth(_h_psi2S);
}

/// @}


/// @name Histograms
/// @{
Histo1DGroupPtr _h_psi2S;
/// @}

};

RIVET_DECLARE_PLUGIN(CMS_2019_I1672011);

} ```