Rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

CMS_2011_S8978280

$K_s$, $\Lambda$, and Cascade$-$ transverse momentum and rapidity spectra at 900 and 7000 GeV.
Experiment: CMS (LHC)
Inspire ID: 890166
Status: VALIDATED
Authors:
  • Kevin Stenson
References: 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 $K_S$, $\Lambda$, and Cascade- particles were measured versus transverse-momentum ($p_\perp$) 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 $p_\perp$ range from 0 to 10 GeV ($K_S$ and $\Lambda$) and 0 to 6 GeV (Cascade-). Antiparticles are included in all measurements so only the sums of $\Lambda$ and $\bar{\Lambda}$, and Cascade$-$ and anti-Cascade$-$ are given. The rapidity distributions are shown versus $|y|$ but normalized to a unit of $y$. Ratios of $\Lambda$/$K_S$ and Cascade$-$/$\Lambda$ production versus $p_\perp$ 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_S8978280.cc
  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/UnstableFinalState.hh"

namespace Rivet {


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

    /// Constructor
    CMS_2011_S8978280()
      : Analysis("CMS_2011_S8978280")
    { }


    void init() {
      UnstableFinalState ufs(Cuts::absrap < 2);
      declare(ufs, "UFS");

      // Particle distributions versus rapidity and transverse momentum
      if (fuzzyEquals(sqrtS()/GeV, 900*GeV)){
        _h_dNKshort_dy  = bookHisto1D(1, 1, 1);
        _h_dNKshort_dpT = bookHisto1D(2, 1, 1);
        _h_dNLambda_dy  = bookHisto1D(3, 1, 1);
        _h_dNLambda_dpT = bookHisto1D(4, 1, 1);
        _h_dNXi_dy      = bookHisto1D(5, 1, 1);
        _h_dNXi_dpT     = bookHisto1D(6, 1, 1);
        //
        _h_LampT_KpT    = bookScatter2D(7, 1, 1);
        _h_XipT_LampT   = bookScatter2D(8, 1, 1);
        _h_Lamy_Ky      = bookScatter2D(9, 1, 1);
        _h_Xiy_Lamy     = bookScatter2D(10, 1, 1);

      } else if (fuzzyEquals(sqrtS()/GeV, 7000*GeV)){
        _h_dNKshort_dy  = bookHisto1D(1, 1, 2);
        _h_dNKshort_dpT = bookHisto1D(2, 1, 2);
        _h_dNLambda_dy  = bookHisto1D(3, 1, 2);
        _h_dNLambda_dpT = bookHisto1D(4, 1, 2);
        _h_dNXi_dy      = bookHisto1D(5, 1, 2);
        _h_dNXi_dpT     = bookHisto1D(6, 1, 2);
        //
        _h_LampT_KpT    = bookScatter2D(7, 1, 2);
        _h_XipT_LampT   = bookScatter2D(8, 1, 2);
        _h_Lamy_Ky      = bookScatter2D(9, 1, 2);
        _h_Xiy_Lamy     = bookScatter2D(10, 1, 2);
      }
    }


    void analyze(const Event& event) {
      const double weight = event.weight();

      const UnstableFinalState& parts = apply<UnstableFinalState>(event, "UFS");
      foreach (const Particle& p, parts.particles()) {
        switch (p.abspid()) {
        case PID::K0S:
          _h_dNKshort_dy->fill(p.absrap(), weight);
          _h_dNKshort_dpT->fill(p.pT(), weight);
          break;

        case PID::LAMBDA:
          // Lambda should not have Cascade or Omega ancestors since they should not decay. But just in case...
          if ( !( p.hasAncestor(3322) || p.hasAncestor(-3322) || p.hasAncestor(3312) || p.hasAncestor(-3312) || p.hasAncestor(3334) || p.hasAncestor(-3334) ) ) {
            _h_dNLambda_dy->fill(p.absrap(), weight);
            _h_dNLambda_dpT->fill(p.pT(), weight);
          }
          break;

        case PID::XIMINUS:
          // Cascade should not have Omega ancestors since it should not decay.  But just in case...
          if ( !( p.hasAncestor(3334) || p.hasAncestor(-3334) ) ) {
            _h_dNXi_dy->fill(p.absrap(), weight);
            _h_dNXi_dpT->fill(p.pT(), weight);
          }
          break;
        }

      }
    }


    void finalize() {
      divide(_h_dNLambda_dpT,_h_dNKshort_dpT, _h_LampT_KpT);
      divide(_h_dNXi_dpT,_h_dNLambda_dpT, _h_XipT_LampT);
      divide(_h_dNLambda_dy,_h_dNKshort_dy, _h_Lamy_Ky);
      divide(_h_dNXi_dy,_h_dNLambda_dy, _h_Xiy_Lamy);
      const double normpT = 1.0/sumOfWeights();
      const double normy = 0.5*normpT; // Accounts for using |y| instead of y
      scale(_h_dNKshort_dy, normy);
      scale(_h_dNKshort_dpT, normpT);
      scale(_h_dNLambda_dy, normy);
      scale(_h_dNLambda_dpT, normpT);
      scale(_h_dNXi_dy, normy);
      scale(_h_dNXi_dpT, normpT);
    }


  private:

    // Particle distributions versus rapidity and transverse momentum
    Histo1DPtr _h_dNKshort_dy, _h_dNKshort_dpT, _h_dNLambda_dy, _h_dNLambda_dpT, _h_dNXi_dy, _h_dNXi_dpT;
    Scatter2DPtr _h_LampT_KpT, _h_XipT_LampT, _h_Lamy_Ky, _h_Xiy_Lamy;

  };



  // The hook for the plugin system
  DECLARE_RIVET_PLUGIN(CMS_2011_S8978280);

}