ALICE_2014_I1300380.cc

Go to the documentation of this file.

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/UnstableFinalState.hh"
namespace Rivet {


  class ALICE_2014_I1300380 : public Analysis {
  public:

    ALICE_2014_I1300380()
      : Analysis("ALICE_2014_I1300380")
    {}


  public:

    void init() {
      const UnstableFinalState cfs(Cuts::absrap<0.5);
      addProjection(cfs, "CFS");

      // Plots from the paper
      _histPtSigmaStarPlus        = bookHisto1D("d01-x01-y01");    // Sigma*+
      _histPtSigmaStarMinus       = bookHisto1D("d01-x01-y02");    // Sigma*- 
      _histPtSigmaStarPlusAnti    = bookHisto1D("d01-x01-y03");    // anti Sigma*-
      _histPtSigmaStarMinusAnti   = bookHisto1D("d01-x01-y04");    // anti Sigma*+
      _histPtXiStar               = bookHisto1D("d02-x01-y01");    // 0.5 * (xi star + anti xi star)
      _histAveragePt              = bookProfile1D("d03-x01-y01");  // <pT> profile
    }


    void analyze(const Event& event) {
      const double weight = event.weight();
      const UnstableFinalState& cfs = applyProjection<UnstableFinalState>(event, "CFS");
      foreach (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-/+     
    {   
      int aid = abs(p.pdgId());
      if (aid == 211  || // pi+ 
          aid == 321  || // K+
          aid == 313  || // K*(892)0
          aid == 2212 || // proton
          aid == 333  ) {  // phi(1020)
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
      }
    } // end if "rejection of long-lived particles"
      
      
        switch (p.pdgId()) {
      case 3224:
        _histPtSigmaStarPlus->fill(p.pT()/GeV, weight);
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case -3224:
        _histPtSigmaStarPlusAnti->fill(p.pT()/GeV, weight);
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case 3114:
        _histPtSigmaStarMinus->fill(p.pT()/GeV, weight);
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case -3114:
        _histPtSigmaStarMinusAnti->fill(p.pT()/GeV, weight);
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case 3324:
        _histPtXiStar->fill(p.pT()/GeV, weight);
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case -3324:
        _histPtXiStar->fill(p.pT()/GeV, weight);
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case 3312:
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case -3312:
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case 3334:
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
      case -3334:
        _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV, weight);
        break;
        }
      }
    }


    void finalize() {
      scale(_histPtSigmaStarPlus,       1./sumOfWeights());
      scale(_histPtSigmaStarPlusAnti,   1./sumOfWeights());
      scale(_histPtSigmaStarMinus,      1./sumOfWeights());
      scale(_histPtSigmaStarMinusAnti,  1./sumOfWeights());
      scale(_histPtXiStar,              1./sumOfWeights()/ 2.); 
    }


  private:
    // plots from the paper
    Histo1DPtr   _histPtSigmaStarPlus;
    Histo1DPtr   _histPtSigmaStarPlusAnti;
    Histo1DPtr   _histPtSigmaStarMinus;
    Histo1DPtr   _histPtSigmaStarMinusAnti;
    Histo1DPtr   _histPtXiStar;
    Profile1DPtr _histAveragePt;
  };



  // The hook for the plugin system
  DECLARE_RIVET_PLUGIN(ALICE_2014_I1300380);
}