|
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
|
#ifndef __Cent_Bin_h__
#define __Cent_Bin_h__
#include <TNamed.h>
#include <TFile.h>
#include <vector>
#include <map>
class CBin : public TObject {
public:
CBin() { ; }
~CBin() override { ; }
float bin_edge;
float n_part_mean;
float n_part_var;
float n_coll_mean;
float n_coll_var;
float n_hard_mean;
float n_hard_var;
float b_mean;
float b_var;
float eccRP_mean;
float eccRP_var;
float ecc2_mean;
float ecc2_var;
float ecc3_mean;
float ecc3_var;
float s_mean;
float s_var;
ClassDefOverride(CBin, 1)
};
class CentralityBins : public TNamed {
public:
typedef std::map<int, const CentralityBins*> RunMap;
CentralityBins() { ; }
CentralityBins(const char* name, const char* title, int nbins) : TNamed(name, title) {
table_.reserve(nbins);
for (int j = 0; j < nbins; ++j) {
CBin b;
table_.push_back(b);
}
}
~CentralityBins() override { ; }
int getBin(double value) const;
int getNbins() const { return table_.size(); }
float lowEdge(double value) const { return lowEdgeOfBin(getBin(value)); }
float lowEdgeOfBin(int bin) const { return table_[bin].bin_edge; }
float NpartMean(double value) const { return NpartMeanOfBin(getBin(value)); }
float NpartMeanOfBin(int bin) const { return table_[bin].n_part_mean; }
float NpartSigma(double value) const { return NpartSigmaOfBin(getBin(value)); }
float NpartSigmaOfBin(int bin) const { return table_[bin].n_part_var; }
float NcollMean(double value) const { return NcollMeanOfBin(getBin(value)); }
float NcollMeanOfBin(int bin) const { return table_[bin].n_coll_mean; }
float NcollSigma(double value) const { return NcollSigmaOfBin(getBin(value)); }
float NcollSigmaOfBin(int bin) const { return table_[bin].n_coll_var; }
float NhardMean(double value) const { return NhardMeanOfBin(getBin(value)); }
float NhardMeanOfBin(int bin) const { return table_[bin].n_hard_mean; }
float NhardSigma(double value) const { return NhardSigmaOfBin(getBin(value)); }
float NhardSigmaOfBin(int bin) const { return table_[bin].n_hard_var; }
float bMean(double value) const { return bMeanOfBin(getBin(value)); }
float bMeanOfBin(int bin) const { return table_[bin].b_mean; }
float bSigma(double value) const { return bSigmaOfBin(getBin(value)); }
float bSigmaOfBin(int bin) const { return table_[bin].b_var; }
float eccentricityRPMean(double value) const { return eccentricityRPMeanOfBin(getBin(value)); }
float eccentricityRPMeanOfBin(int bin) const { return table_[bin].eccRP_mean; }
float eccentricityRPSigma(double value) const { return eccentricityRPSigmaOfBin(getBin(value)); }
float eccentricityRPSigmaOfBin(int bin) const { return table_[bin].eccRP_var; }
float eccentricityMean(double value) const { return eccentricityMeanOfBin(getBin(value)); }
float eccentricityMeanOfBin(int bin) const { return table_[bin].ecc2_mean; }
float eccentricitySigma(double value) const { return eccentricitySigmaOfBin(getBin(value)); }
float eccentricitySigmaOfBin(int bin) const { return table_[bin].ecc2_var; }
float triangularityMean(double value) const { return triangularityMeanOfBin(getBin(value)); }
float triangularityMeanOfBin(int bin) const { return table_[bin].ecc3_mean; }
float triangularitySigma(double value) const { return triangularitySigmaOfBin(getBin(value)); }
float triangularitySigmaOfBin(int bin) const { return table_[bin].ecc3_var; }
float areaMean(double value) const { return areaMeanOfBin(getBin(value)); }
float areaMeanOfBin(int bin) const { return table_[bin].s_mean; }
float areaSigma(double value) const { return areaSigmaOfBin(getBin(value)); }
float areaSigmaOfBin(int bin) const { return table_[bin].s_var; }
// private:
std::vector<CBin> table_;
ClassDefOverride(CentralityBins, 1)
};
CentralityBins::RunMap getCentralityFromFile(
TDirectoryFile*, const char* dir, const char* tag, int firstRun = 0, int lastRun = 10);
CentralityBins::RunMap getCentralityFromFile(TDirectoryFile*, const char* tag, int firstRun = 0, int lastRun = 10);
#endif
|