Rivet Analyses Reference

  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

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Tools/ParticleName.hh"

namespace Rivet {


  /// Measurement of charged pion, kaon, and proton production in 13 TeV pp collisions
  class CMS_2017_I1608166 : public Analysis {
  public:

    /// Constructor
    RIVET_DEFAULT_ANALYSIS_CTOR(CMS_2017_I1608166);


    /// @name Analysis methods
    //@{

    /// Book histograms and initialise projections before the run
    void init() {
      const ChargedFinalState cfs(Cuts::absrap < 1.);
      declare(cfs, "CFS");
      //
      // pt spectra
      book(_h[PID::PIPLUS],  "d01-x01-y01");
      book(_h[PID::KPLUS],   "d01-x01-y02");
      book(_h[PID::PROTON],  "d01-x01-y03");
      book(_h[PID::PIMINUS], "d02-x01-y01");
      book(_h[PID::KMINUS],  "d02-x01-y02");
      book(_h[PID::PBAR],    "d02-x01-y03");

      // negative/positive ratios
      book(_s["pi-/pi+"], "d43-x01-y01");
      book(_s["k-/k+"],   "d44-x01-y01");
      book(_s["p~/p"],    "d45-x01-y01");

      // k/pi and p/pi ratios
      book(_hkpi[PID::PIPLUS], "TMP/hkpi/pi", refData(46, 1, 1));
      book(_hkpi[PID::KPLUS],  "TMP/hkpi/k",  refData(46, 1, 1));
      book(_hppi[PID::PIPLUS], "TMP/hppi/pi", refData(47, 1, 1));
      book(_hppi[PID::PROTON], "TMP/hppi/p",  refData(47, 1, 1));
      book(_s["k/pi"],    "d46-x01-y01");
      book(_s["p/pi"],    "d47-x01-y01");
    }


    void analyze(const Event& event) {

      const ChargedFinalState& cfs = apply<ChargedFinalState>(event, "CFS");
      for (const Particle& p : cfs.particles()) {

        // protections against MC generators decaying long-lived particles
        if (p.hasAncestor(310)  || p.hasAncestor(-310)  ||  // K0s
            p.hasAncestor(130)  || p.hasAncestor(-130)  ||  // K0l
            p.hasAncestor(3322) || p.hasAncestor(-3322) ||  // Xi0
            p.hasAncestor(3122) || p.hasAncestor(-3122) ||  // Lambda
            p.hasAncestor(3222) || p.hasAncestor(-3222) ||  // Sigma+/-
            p.hasAncestor(3312) || p.hasAncestor(-3312) ||  // Xi-/+
            p.hasAncestor(3334) || p.hasAncestor(-3334))    // Omega-/+
          continue;

        if (theParticles.find(p.pid()) != theParticles.end()) {
          // fill pt spectra
          _h[p.pid()]->fill(p.pt() / GeV);
          // fill tmp histos for ratios
          if (p.abspid() != PID::PROTON)
            _hkpi[p.abspid()]->fill(p.pt() / GeV);
          if (p.abspid() != PID::KPLUS)
            _hppi[p.abspid()]->fill(p.pt() / GeV);
        }

      }

    }


    void finalize() {

      divide(_h[PID::PIMINUS], _h[PID::PIPLUS], _s["pi-/pi+"]);
      divide(_h[PID::KMINUS],  _h[PID::KPLUS],  _s["k-/k+"]);
      divide(_h[PID::PBAR],    _h[PID::PROTON], _s["p~/p"]);

      divide(_hkpi[PID::KPLUS],  _hkpi[PID::PIPLUS], _s["k/pi"]);
      divide(_hppi[PID::PROTON], _hppi[PID::PIPLUS], _s["p/pi"]);

      scale(_h, 1./2./sumOfWeights());

    }


    set<int> theParticles = {PID::PIPLUS, PID::KPLUS, PID::PROTON, PID::PIMINUS, PID::KMINUS, PID::PBAR};

    map<int, Histo1DPtr> _h;
    map<int, Histo1DPtr> _hkpi, _hppi;
    map<string, Scatter2DPtr> _s;

  };


  RIVET_DECLARE_PLUGIN(CMS_2017_I1608166);

}