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 enrico.fragiacomo@ts.infn.it
  • Andreas Morsch andreas.morsch@cern.ch
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_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm 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
  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
115
116
117
118
119
120
// -*- 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);
      declare(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 = apply<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);
}