Rivet analyses

title: CMS_2022_I1870319

Prompt J/$\psi$ production in jets at 5.02 TeV

Experiment: CMS (LHC)

Inspire ID: 1870319

Status: VALIDATED

Authors: - Peter Richardson

References: - Phys.Lett.B 825 (2022) 136842 - Expt page: CMS-HIN-19-007 - arXiv: 2106.13235

Beams: p+ p+

Beam energies: (2510.0, 2510.0)GeV

Run details: - prompt J/psi production

Measurement of the fragmentation function for the production of J/$\psi$ in jets at 5.02 TeV by CMS

Source code:CMS_2022_I1870319.cc

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

include "Rivet/Analysis.hh"

include "Rivet/Projections/FinalState.hh"

include "Rivet/Projections/UnstableParticles.hh"

include "fastjet/ClusterSequence.hh"

include "fastjet/JetDefinition.hh"

namespace Rivet {

/// @brief J/psi in jets at 5.02 TeV class CMS_2022_I1870319 : public Analysis { public:

/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2022_I1870319);


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

/// Book histograms and initialise projections before the run
void init() {
  declare(UnstableParticles(), "UFS");
  declare(FinalState(), "FS");
  for (unsigned int ix = 0; ix < 2; ++ix) {
    book(_h_frag[ix], ix + 1, 1, 1);
  }
}


/// Perform the per-event analysis
void analyze(const Event& event) {
  // first see if we have any prompt J/psi in the region
  Particles Jpsi;
  for (const Particle& p :
       apply<UnstableParticles>(event, "UFS").particles(Cuts::pid == 443 && Cuts::pT > 6.5)) {
    if (!p.fromBottom()) Jpsi.push_back(p);
  }
  // no jpsi veto
  if (Jpsi.empty()) vetoEvent;
  // now get the final-state particles for the jets
  const Particles& fsParticles = apply<FinalState>(event, "FS").particles();
  vector<PseudoJet> particles;
  particles.reserve(fsParticles.size());
  // fs for fastjet omitting any J/psi decay products
  for (const Particle& p : fsParticles) {
    if (p.abspid() == 12 || p.abspid() == 14 or p.abspid() == 16) continue;
    // skip anything coming from the decay of one of the jpsis
    Particle parent = p;
    while (!parent.parents().empty()) {
      if (parent.pid() == 443) break;
      parent = parent.parents()[0];
    }
    bool match = parent.pid() == 443;
    if (match) {
      match = false;
      for (const Particle& psi : Jpsi) {
        match = fuzzyEquals(parent.momentum(), psi.momentum());
        if (match) break;
      }
    }
    if (!match) {
      PseudoJet j = p.pseudojet();
      j.set_user_index(0);
      particles.push_back(j);
    }
  }
  // add the jpsis to the particles for fastjet
  for (const Particle& p : Jpsi) {
    PseudoJet j = p.pseudojet();
    j.set_user_index(1);
    particles.push_back(j);
  }
  JetDefinition jet_def(fastjet::antikt_algorithm, 0.3);
  fastjet::ClusterSequence clu = ClusterSequence(particles, jet_def);
  vector<PseudoJet> jets = clu.inclusive_jets();
  for (const PseudoJet& jet : jets) {
    // pt and eta cut
    if (jet.perp() < 30. || jet.perp() > 40. || abs(jet.eta()) > 2.) continue;
    // loop over constituents and find jpsi
    for (const PseudoJet& sub : jet.constituents()) {
      if (sub.user_index() == 0) continue;
      const double z = sub.perp() / jet.perp();
      _h_frag[0]->fill(z);
      _h_frag[1]->fill(z);
    }
  }
}


/// Normalise histograms etc., after the run
void finalize() {
  normalize(_h_frag[0]);
  scale(_h_frag[1], crossSection() / nanobarn / sumOfWeights());
}

/// @}


/// @name Histograms
/// @{
Histo1DPtr _h_frag[2];
/// @}

};

RIVET_DECLARE_PLUGIN(CMS_2022_I1870319);

} ```