Rivet analyses referenceALICE_2014_I1300380Production of $\Sigma(1385)^{\pm}$ and $\Xi(1530)^{0}$ in proton-proton collisions at $\sqrt{s} =$ 7 TeVExperiment: ALICE (LHC) Inspire ID: 1300380 Status: VALIDATED Authors:
Beams: p+ p+ Beam energies: (3500.0, 3500.0) GeV Run details:
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 if(_edges.empty()) _edges=_histAveragePt->xEdges();
33 const UnstableParticles& cfs = apply<UnstableParticles>(event, "CFS");
34 for (const Particle& p : cfs.particles()) {
35 // protections against mc generators decaying long-lived particles
36 if ( !(p.hasAncestorWith(Cuts::pid == 310) || p.hasAncestorWith(Cuts::pid == -310) || // K0s
37 p.hasAncestorWith(Cuts::pid == 130) || p.hasAncestorWith(Cuts::pid == -130) || // K0l
38 p.hasAncestorWith(Cuts::pid == 3322) || p.hasAncestorWith(Cuts::pid == -3322) || // Xi0
39 p.hasAncestorWith(Cuts::pid == 3122) || p.hasAncestorWith(Cuts::pid == -3122) || // Lambda
40 p.hasAncestorWith(Cuts::pid == 3222) || p.hasAncestorWith(Cuts::pid == -3222) || // Sigma+/-
41 p.hasAncestorWith(Cuts::pid == 3312) || p.hasAncestorWith(Cuts::pid == -3312) || // Xi-/+
42 p.hasAncestorWith(Cuts::pid == 3334) || p.hasAncestorWith(Cuts::pid == -3334) )) // Omega-/+
43 {
44 int aid = p.abspid();
45 if (aid == 211 ) _histAveragePt->fill(_edges[0], p.pT()/GeV); // pi+
46 else if (aid == 321 ) _histAveragePt->fill(_edges[1], p.pT()/GeV); // K+
47 else if (aid == 313 ) _histAveragePt->fill(_edges[2], p.pT()/GeV); // K*(892)0
48 else if (aid == 2212) _histAveragePt->fill(_edges[3], p.pT()/GeV); // proton
49 else if (aid == 333 ) _histAveragePt->fill(_edges[4], p.pT()/GeV); // phi(1020)
50 } // end if "rejection of long-lived particles"
51 switch (p.pid()) {
52 case 3224:
53 _histPtSigmaStarPlus->fill(p.pT()/GeV);
54 _histAveragePt->fill(_edges[6], p.pT()/GeV);
55 break;
56 case -3224:
57 _histPtSigmaStarPlusAnti->fill(p.pT()/GeV);
58 _histAveragePt->fill(_edges[6], p.pT()/GeV);
59 break;
60 case 3114:
61 _histPtSigmaStarMinus->fill(p.pT()/GeV);
62 _histAveragePt->fill(_edges[6], p.pT()/GeV);
63 break;
64 case -3114:
65 _histPtSigmaStarMinusAnti->fill(p.pT()/GeV);
66 _histAveragePt->fill(_edges[6], p.pT()/GeV);
67 break;
68 case 3324:
69 _histPtXiStar->fill(p.pT()/GeV);
70 _histAveragePt->fill(_edges[7], p.pT()/GeV);
71 break;
72 case -3324:
73 _histPtXiStar->fill(p.pT()/GeV);
74 _histAveragePt->fill(_edges[7], p.pT()/GeV);
75 break;
76 case 3312:
77 _histAveragePt->fill(_edges[5], p.pT()/GeV);
78 break;
79 case -3312:
80 _histAveragePt->fill(_edges[5], p.pT()/GeV);
81 break;
82 case 3334:
83 _histAveragePt->fill(_edges[8], p.pT()/GeV);
84 break;
85 case -3334:
86 _histAveragePt->fill(_edges[8], p.pT()/GeV);
87 break;
88 }
89 }
90 }
91
92
93 void finalize() {
94 scale(_histPtSigmaStarPlus, 1./sumOfWeights());
95 scale(_histPtSigmaStarPlusAnti, 1./sumOfWeights());
96 scale(_histPtSigmaStarMinus, 1./sumOfWeights());
97 scale(_histPtSigmaStarMinusAnti, 1./sumOfWeights());
98 scale(_histPtXiStar, 1./sumOfWeights()/ 2.);
99 }
100
101
102 private:
103 // plots from the paper
104 Histo1DPtr _histPtSigmaStarPlus;
105 Histo1DPtr _histPtSigmaStarPlusAnti;
106 Histo1DPtr _histPtSigmaStarMinus;
107 Histo1DPtr _histPtSigmaStarMinusAnti;
108 Histo1DPtr _histPtXiStar;
109 BinnedProfilePtr<string> _histAveragePt;
110 vector<string> _edges;
111 };
112
113
114 RIVET_DECLARE_PLUGIN(ALICE_2014_I1300380);
115
116}
|