Rivet analyses

Ks, Λ, and Cascade− transverse momentum and rapidity spectra at 900 and 7000 GeV.

Experiment: CMS (LHC)

Inspire ID: 890166

Status: VALIDATED

Authors: - Kevin Stenson

References: - JHEP 05 (2011) 064 - DOI: 10.1007/JHEP05(2011)064 - arXiv: 1102.4282

Beams: p+ p+

Beam energies: (450.0, 450.0); (3500.0, 3500.0)GeV

Run details: - Non-single-diffractive (NSD) events only. Should include double-diffractive (DD) events and non-diffractive (ND) events but NOT single-diffractive (SD) events. For example, in Pythia6 the SD processes to be turned off are 92 and 93, and in Pythia8 the SD processes are 103 and 104 (also called SoftQCD:singleDiffractive).

The spectra of KS, Λ, and Cascade- particles were measured versus transverse-momentum ($\pT$) and rapidity (y) in proton-proton collisions at center-of-mass energies 900 and 7000 GeV. The production is normalized to all non-single-diffractive (NSD) events using corrections for trigger and selection efficiency, acceptance, and branching ratios. The results cover a rapidity range of |y| < 2 and a $\pT$ range from 0 to 10 GeV (KS and Λ) and 0 to 6 GeV (Cascade-). Antiparticles are included in all measurements so only the sums of Λ and Λ̄, and Cascade− and anti-Cascade− are given. The rapidity distributions are shown versus |y| but normalized to a unit of y. Ratios of Λ/KS and Cascade−/Λ production versus $\pT$ and |y| are also given, with somewhat smaller systematic uncertainties than obtained from taking the ratio of the individual distributions.’ The data were corrected according to the SD/DD/ND content of the CMS trigger, as predicted by PYTHIA6. The uncertainties connected with correct or uncorrect modelling of diffraction were included in the systematic errors.

Source code:CMS_2011_I890166.cc

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/UnstableParticles.hh"

namespace Rivet {


  /// @brief CMS strange particle spectra (Ks, Lambda, Cascade) in pp at 900 and 7000 GeV
  ///
  /// @author Kevin Stenson
  class CMS_2011_I890166 : public Analysis {
  public:

    RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2011_I890166);


    void init() {

      UnstableParticles ufs(Cuts::absrap < 2);
      declare(ufs, "UFS");

      // Particle distributions versus rapidity and transverse momentum
      for (double eVal : allowedEnergies()) {
        const string en = toString(round(eVal));
        if (isCompatibleWithSqrtS(eVal))  _sqs = en;
        bool offset(en == "7000"s);

        book(_h[en+"dNKshort_dy"],  1, 1, 1+offset);
        book(_h[en+"dNKshort_dpT"], 2, 1, 1+offset);
        book(_h[en+"dNLambda_dy"],  3, 1, 1+offset);
        book(_h[en+"dNLambda_dpT"], 4, 1, 1+offset);
        book(_h[en+"dNXi_dy"],      5, 1, 1+offset);
        book(_h[en+"dNXi_dpT"],     6, 1, 1+offset);
        //
        book(_e[en+"LampT_KpT"],   7, 1, 1+offset);
        book(_e[en+"XipT_LampT"],  8, 1, 1+offset);
        book(_e[en+"Lamy_Ky"],     9, 1, 1+offset);
        book(_e[en+"Xiy_Lamy"],   10, 1, 1+offset);

      }
      raiseBeamErrorIf(_sqs.empty());
    }


    void analyze(const Event& event) {

      const UnstableParticles& parts = apply<UnstableParticles>(event, "UFS");
      for (const Particle& p : parts.particles()) {
        switch (p.abspid()) {
        case PID::K0S:
          _h[_sqs+"dNKshort_dy"]->fill(p.absrap());
          _h[_sqs+"dNKshort_dpT"]->fill(p.pT()/GeV);
          break;

        case PID::LAMBDA:
          // Lambda should not have Cascade or Omega ancestors since they should not decay. But just in case...
          if ( !( p.hasAncestorWith(Cuts::pid ==  3322) ||
                  p.hasAncestorWith(Cuts::pid == -3322) ||
                  p.hasAncestorWith(Cuts::pid ==  3312) ||
                  p.hasAncestorWith(Cuts::pid == -3312) ||
                  p.hasAncestorWith(Cuts::pid ==  3334) ||
                  p.hasAncestorWith(Cuts::pid == -3334) ) ) {
            _h[_sqs+"dNLambda_dy"]->fill(p.absrap());
            _h[_sqs+"dNLambda_dpT"]->fill(p.pT()/GeV);
          }
          break;

        case PID::XIMINUS:
          // Cascade should not have Omega ancestors since it should not decay.  But just in case...
          if ( !( p.hasAncestorWith(Cuts::pid ==  3334) ||
                  p.hasAncestorWith(Cuts::pid == -3334) ) ) {
            _h[_sqs+"dNXi_dy"]->fill(p.absrap());
            _h[_sqs+"dNXi_dpT"]->fill(p.pT()/GeV);
          }
          break;
        }
      }
    }


    void finalize() {

      for (double eVal : allowedEnergies()) {

        const string en = toString(round(eVal));

        divide(_h[en+"dNLambda_dpT"], _h[en+"dNKshort_dpT"], _e[en+"LampT_KpT"]);

        YODA::Histo1D denom = _h[en+"dNLambda_dpT"]->clone();
        denom.rebinXTo(_h[en+"dNXi_dpT"]->xEdges());
        divide(*_h[en+"dNXi_dpT"], denom, _e[en+"XipT_LampT"]);

        divide(_h[en+"dNLambda_dy"], _h[en+"dNKshort_dy"], _e[en+"Lamy_Ky"]);
        divide(_h[en+"dNXi_dy"], _h[en+"dNLambda_dy"], _e[en+"Xiy_Lamy"]);
      }

      const double normpT = 1.0/sumOfWeights();
      const double normy = 0.5*normpT; // Accounts for using |y| instead of y
      for (auto& item : _h) {
        if (item.first.find("_dy") != string::npos) {
          scale(item.second, normy);
        }
        else {
          scale(item.second, normpT);
        }
      }
    }


  private:

    /// @name Particle distributions versus rapidity and transverse momentum
    /// @{
    map<string,Histo1DPtr> _h;
    map<string,Estimate1DPtr> _e;

    string _sqs = "";
    /// @}

  };



  RIVET_DECLARE_ALIASED_PLUGIN(CMS_2011_I890166, CMS_2011_S8978280);

}

Aliases: - CMS_2011_S8978280