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_{\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// -*- 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.hasAncestor(310) || p.hasAncestor(-310) || // K0s
36 p.hasAncestor(130) || p.hasAncestor(-130) || // K0l
37 p.hasAncestor(3322) || p.hasAncestor(-3322) || // Xi0
38 p.hasAncestor(3122) || p.hasAncestor(-3122) || // Lambda
39 p.hasAncestor(3222) || p.hasAncestor(-3222) || // Sigma+/-
40 p.hasAncestor(3312) || p.hasAncestor(-3312) || // Xi-/+
41 p.hasAncestor(3334) || p.hasAncestor(-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 // The hook for the plugin system
118 RIVET_DECLARE_PLUGIN(ALICE_2014_I1300380);
119
120}
|