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
| //-*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/UnstableParticles.hh"
#include "Rivet/Tools/ParticleUtils.hh"
namespace Rivet {
class ALICE_2017_I1620477 : public Analysis {
public:
/// Constructor
ALICE_2017_I1620477()
: Analysis("ALICE_2017_I1620477"),
_rapmax(0.8)
{ }
void init() {
const UnstableParticles ufs(Cuts::absrap < _rapmax);
declare(ufs, "UFS");
book(_h_pi0,1,1,1);
book(_h_eta,2,1,1);
book(_h_etaToPion,8,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(8,1,1));
book(_temp_h_eta, "TMP/h_eta", refData(8,1,1));
}
void analyze(const Event& event) {
const UnstableParticles& ufs = applyProjection<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. /(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. /(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_I1620477);
}
|