rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

ALICE_2014_I1300380

Production of $\Sigma(1385)^{\pm}$ and $\Xi(1530)^{0}$ in proton-proton collisions at $\sqrt{s} =$ 7 TeV
Experiment: ALICE (LHC)
Inspire ID: 1300380
Status: VALIDATED
Authors:
  • Enrico Fragiacomo
  • Andreas Morsch
No references listed
Beams: p+ p+
Beam energies: (3500.0, 3500.0) GeV
Run details:
  • inelastic collisions pp at sqrt(s) = 7 TeV

Transverse momentum spectra $\frac{1}{N_\mathrm{inel}} \frac{\mathrm{d}^2N}{\mathrm{d}p_\mathrm{T}\mathrm{d}y}$ of $\Sigma(1385)^{\pm}$, $\overline{\Sigma}(1385)^{\pm}$ and $(\Xi(1530)^{0}+\overline{\Xi}(1530)^{0})/2$ at mid-rapidity ($|y| < 0.5$) in inelastic proton-proton collisions at $\sqrt{s} =$ 7 TeV

Source code: ALICE_2014_I1300380.cc
  1// -*- C++ -*-
  2#include "Rivet/Analysis.hh"
  3#include "Rivet/Projections/UnstableParticles.hh"
  4namespace Rivet {
  5
  6
  7  class ALICE_2014_I1300380 : public Analysis {
  8  public:
  9
 10    ALICE_2014_I1300380()
 11      : Analysis("ALICE_2014_I1300380")
 12    {}
 13
 14
 15  public:
 16
 17    void init() {
 18      const UnstableParticles cfs(Cuts::absrap<0.5);
 19      declare(cfs, "CFS");
 20
 21      // Plots from the paper
 22      book(_histPtSigmaStarPlus        ,"d01-x01-y01");    // Sigma*+
 23      book(_histPtSigmaStarMinus       ,"d01-x01-y02");    // Sigma*-
 24      book(_histPtSigmaStarPlusAnti    ,"d01-x01-y03");    // anti Sigma*-
 25      book(_histPtSigmaStarMinusAnti   ,"d01-x01-y04");    // anti Sigma*+
 26      book(_histPtXiStar               ,"d02-x01-y01");    // 0.5 * (xi star + anti xi star)
 27      book(_histAveragePt              ,"d03-x01-y01");  // <pT> profile
 28    }
 29
 30
 31    void analyze(const Event& event) {
 32      const UnstableParticles& cfs = apply<UnstableParticles>(event, "CFS");
 33      for (const Particle& p : cfs.particles()) {
 34	// protections against mc generators decaying long-lived particles
 35	if ( !(p.hasAncestorWith(Cuts::pid == 310)  || p.hasAncestorWith(Cuts::pid == -310)   || // K0s
 36	       p.hasAncestorWith(Cuts::pid == 130)  || p.hasAncestorWith(Cuts::pid == -130)   ||     // K0l
 37	       p.hasAncestorWith(Cuts::pid == 3322) || p.hasAncestorWith(Cuts::pid == -3322)  ||     // Xi0
 38	       p.hasAncestorWith(Cuts::pid == 3122) || p.hasAncestorWith(Cuts::pid == -3122)  ||     // Lambda
 39	       p.hasAncestorWith(Cuts::pid == 3222) || p.hasAncestorWith(Cuts::pid == -3222)  ||     // Sigma+/-
 40	       p.hasAncestorWith(Cuts::pid == 3312) || p.hasAncestorWith(Cuts::pid == -3312)  ||     // Xi-/+
 41	       p.hasAncestorWith(Cuts::pid == 3334) || p.hasAncestorWith(Cuts::pid == -3334)  ))     // Omega-/+
 42	{
 43	  int aid = p.abspid();
 44	  if (aid == 211  || // pi+
 45            aid == 321  || // K+
 46            aid == 313  || // K*(892)0
 47            aid == 2212 || // proton
 48            aid == 333  )  // phi(1020)
 49          {
 50	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 51	  }
 52	} // end if "rejection of long-lived particles"
 53
 54
 55        switch (p.pid()) {
 56	  case 3224:
 57	    _histPtSigmaStarPlus->fill(p.pT()/GeV);
 58	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 59	    break;
 60	  case -3224:
 61	    _histPtSigmaStarPlusAnti->fill(p.pT()/GeV);
 62	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 63	    break;
 64	  case 3114:
 65	    _histPtSigmaStarMinus->fill(p.pT()/GeV);
 66	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 67	    break;
 68	  case -3114:
 69	    _histPtSigmaStarMinusAnti->fill(p.pT()/GeV);
 70	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 71	    break;
 72	  case 3324:
 73	    _histPtXiStar->fill(p.pT()/GeV);
 74	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 75	    break;
 76	  case -3324:
 77	    _histPtXiStar->fill(p.pT()/GeV);
 78	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 79	    break;
 80	  case 3312:
 81	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 82	    break;
 83	  case -3312:
 84	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 85	    break;
 86	  case 3334:
 87	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 88	    break;
 89	  case -3334:
 90	    _histAveragePt->fill(p.mass()/GeV, p.pT()/GeV);
 91	    break;
 92        }
 93      }
 94    }
 95
 96
 97    void finalize() {
 98      scale(_histPtSigmaStarPlus,       1./sumOfWeights());
 99      scale(_histPtSigmaStarPlusAnti,   1./sumOfWeights());
100      scale(_histPtSigmaStarMinus,      1./sumOfWeights());
101      scale(_histPtSigmaStarMinusAnti,  1./sumOfWeights());
102      scale(_histPtXiStar,              1./sumOfWeights()/ 2.);
103    }
104
105
106  private:
107    // plots from the paper
108    Histo1DPtr   _histPtSigmaStarPlus;
109    Histo1DPtr   _histPtSigmaStarPlusAnti;
110    Histo1DPtr   _histPtSigmaStarMinus;
111    Histo1DPtr   _histPtSigmaStarMinusAnti;
112    Histo1DPtr   _histPtXiStar;
113    Profile1DPtr _histAveragePt;
114  };
115
116
117  RIVET_DECLARE_PLUGIN(ALICE_2014_I1300380);
118
119}