|
|
Rivet analyses reference
ALICE_2017_I1512110
Invariant cross section of neutral pions measured at mid-rapidity in pp collisions at $2.76\,$TeV.
Experiment: ALICE (LHC)
Inspire ID: 1512110
Status: VALIDATED
No authors listed
No references listed
Beams: p+ p+
Beam energies: (1380.0, 1380.0) GeV
Run details:
- Transverse momentum spectra of neutral pions and $\eta$ mesons measured at 0.4<p_T<40 GeV/c and 0.6<p_T<20 GeV/c, respectively, and the ratio $\pi^{0}/\eta$, obtained at mid-rapidity in pp collisions at $\sqrt(s) = 2.76\,$TeV with ALICE at the LHC.
Transverse momentum spectra of neutral pions and $\eta$ mesons measured at 0.4<p_T<40 GeV/c and 0.6<p_T<20 GeV/c, respectively, and the ratio $\pi^{0}/\eta$, obtained at mid-rapidity in pp collisions at $\sqrt(s) = 2.76\,$TeV with ALICE at the LHC.
Source code:
ALICE_2017_I1512110.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 | //-*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/UnstableParticles.hh"
namespace Rivet {
class ALICE_2017_I1512110 : public Analysis {
public:
/// Constructor
ALICE_2017_I1512110()
: Analysis("ALICE_2017_I1512110"),
_rapmax(0.8)
{ }
void init() {
const UnstableParticles ufs(Cuts::absrap < _rapmax);
declare(ufs, "UFS");
book(_h_pi0, 3,1,1);
book(_h_eta, 4,1,1);
book(_h_etaToPion, 5,1,1);
// temporary plots with the binning of _h_etaToPion
// to construct the eta/pi0 ratio in the end
book(_temp_h_pion, "TMP/h_pion",refData(5,1,1));
book(_temp_h_eta , "TMP/h_eta", refData(5,1,1));
}
void analyze(const Event& event) {
const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");
for (const Particle& p : ufs.particles()) {
if (p.pid() == 111) {
// neutral pion; ALICE corrects for pi0 feed-down
if ( !(p.hasAncestor(310) || p.hasAncestor(130) || // K0_s, K0_l
p.hasAncestor(321) || p.hasAncestor(-321) || // K+,K-
p.hasAncestor(3122) || p.hasAncestor(-3122) || // Lambda, Anti-Lambda
p.hasAncestor(3212) || p.hasAncestor(-3212) || // Sigma0
p.hasAncestor(3222) || p.hasAncestor(-3222) || // Sigmas
p.hasAncestor(3112) || p.hasAncestor(-3112) || // Sigmas
p.hasAncestor(3322) || p.hasAncestor(-3322) || // Cascades
p.hasAncestor(3312) || p.hasAncestor(-3312) )) // Cascades
{
_h_pi0->fill(p.pT()/GeV, 1.0 /(TWOPI*p.pT()/GeV*2*_rapmax));
_temp_h_pion->fill(p.pT()/GeV);
}
}
else if (p.pid() == 221){
// eta meson
_h_eta->fill(p.pT()/GeV, 1.0 /(TWOPI*p.pT()/GeV*2*_rapmax));
_temp_h_eta->fill(p.pT()/GeV);
}
}
}
void finalize() {
scale(_h_pi0, crossSection()/picobarn/sumOfWeights());
scale(_h_eta, crossSection()/picobarn/sumOfWeights());
divide(_temp_h_eta, _temp_h_pion, _h_etaToPion);
}
private:
double _rapmax;
Histo1DPtr _h_pi0;
Histo1DPtr _h_eta;
Histo1DPtr _temp_h_pion;
Histo1DPtr _temp_h_eta;
Scatter2DPtr _h_etaToPion;
};
RIVET_DECLARE_PLUGIN(ALICE_2017_I1512110);
}
|
|