-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathpset.h
53 lines (46 loc) · 6.2 KB
/
pset.h
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
#pragma once
#define PSET_H
#include <vector>
//std::vector<Double_t> pt_bins = {0., 0.5, 1., 2., 3., 4., 5., 7., 9., 12., 15., 18., 21., 24., 30., 35., 40., 45., 50.};
std::vector<Double_t> phi_bins = {-3.2, -2.9, -2.6, -2.4, -2.1, -1.8, -1.5, -1.3, -1.2, -1.1, -1.0, -0.9, -0.7, -0.4, -0.2, 0., 0.2, 0.4, 0.7, 0.9, 1.0, 1.1, 1.2, 1.3, 1.5, 1.8, 2.1, 2.4, 2.6, 2.9, 3.2 };
std::vector<Double_t> pt_bins = {
0., 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3,
1.4,
1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4,
2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4,
3.5, 3.6, 3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.3, 4.4,
4.5, 4.6, 4.7, 4.8, 4.9, 5.0, 5.1, 5.2, 5.3, 5.4,
5.5, 5.6, 5.7, 5.8, 5.9, 6.0, 6.1, 6.2, 6.3, 6.4,
6.5, 6.6, 6.7, 6.8, 6.9, 7.0, 7.1, 7.2, 7.3, 7.4,
7.5, 7.6, 7.7, 7.8, 7.9, 8.0, 8.1, 8.2, 8.3, 8.4,
8.5, 8.6, 8.7, 8.8, 8.9, 9., 10.,
11., 12.,13.,14., 15., 18., 21.,
24., 30., 35., 40., 45., 50.
};
std::vector<Double_t> pt_bins_fine = {0.,0.1, 0.3, 0.5, 0.75, 1., 1.5, 2., 2.5, 3., 3.5, 4., 4.5, 5., 6., 7., 8, 9., 10., 11., 12., 13., 14., 15., 16., 17., 18., 19., 20., 21., 22, 23, 24., 25, 27., 29., 32., 38, 42., 45., 50.};
std::vector<Double_t> rap_bins = {-2.4, -2.1, -1.9, -1.75, -1.6, -1.5, -1.35, -1.2, -1, -0.8, -0.6,-0.3, 0, 0.3, 0.6, 0.8, 1., 1.2, 1.35, 1.5, 1.6, 1.75, 1.9, 2.1, 2.4};
//std::vector<Double_t> cent_bins = { 0., 2., 4., 8., 14., 20. ,40, 60, 80, 100, 120, 140, 160, 180, 200};
//std::vector<Double_t> cent_bins = { 0., 2., 4., 8., 14., 16., 18., 20., 22., 24., 26., 28., 30., 32., 34., 36., 38.,40., 45., 50., 55., 60, 80, 100, 120, 140, 160, 180, 200};
//std::vector<Double_t> cent_bins = { 0., 2., 4., 8., 14., 16., 18., 20., 22., 24., 26., 28., 30., 32., 34., 36., 38.,40., 45., 50., 55., 60, 80, 100, 120, 140, 160, 180, 200};
std::vector<Double_t> cent_bins = { 0., 4., 8., 12., 16., 20., 24., 28., 32., 36., 40., 44., 48., 52., 56., 60., 64., 68., 72., 76., 80., 84., 88., 92., 96., 100., 104., 108., 112., 116., 120., 124., 128., 132., 136., 140., 144., 148., 152., 156., 160., 164., 168., 172., 176., 180., 184., 188., 192., 196., 200.};
std::vector<Double_t> derivedPtCuts = { 0.5, 1., 1.5, 2., 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3., 3.1, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7, 3.8, 3.9, 4., 4.5, 5., 5.5, 6., 6.5, 7.};
std::vector<Double_t> derivedMassCuts = { 0.5, 1.0 };
float findNcoll(int hiBin){
const int nbins = 200;
const float Ncoll[nbins] = {1919.78, 1885.19, 1856.99, 1825.04, 1784.78, 1746.71, 1703.82, 1662.5, 1616.36, 1579.02, 1534.86, 1501.37, 1458.44, 1423.87, 1389.69, 1353.63, 1319.72, 1287.51, 1258.84, 1225.49, 1195.78, 1166.55, 1140.71, 1110.7, 1083.73, 1058.86, 1030.54, 1005.22, 978.904, 956.635, 933.093, 911.224, 887.728, 867.682, 845.888, 823.128, 802.103, 780.984, 762.997, 743.821, 724.782, 705.008, 684.918, 669.106, 650.536, 635.356, 617.694, 602.554, 586.821, 569.964, 556.196, 540.463, 525.911, 510.929, 498.874, 484.433, 472.776, 458.498, 444.539, 432.389, 420.628, 407.479, 396.521, 385.629, 374.165, 363.421, 352.297, 344.028, 333.112, 321.054, 312.064, 301.815, 293.326, 283.747, 275.412, 265.697, 258.378, 248.661, 242.796, 234.84, 226.734, 219.617, 212.156, 205.325, 198.439, 191.398, 184.682, 177.485, 171.602, 165.937, 160.209, 154.521, 149.497, 143.61, 137.743, 133.381, 128.111, 124.338, 118.93, 114.865, 110.556, 106.003, 102.698, 97.8558, 93.4104, 90.1664, 86.6085, 83.0393, 79.8933, 76.3645, 73.2382, 70.0091, 67.5678, 64.8457, 61.7102, 59.1422, 56.4246, 53.8929, 51.7162, 49.3104, 47.3088, 44.6109, 42.6843, 40.8715, 39.0845, 37.3144, 35.425, 33.7925, 32.1438, 30.715, 29.0353, 27.5403, 26.4516, 25.3269, 24.0873, 22.8437, 21.6271, 20.4867, 19.5775, 18.638, 17.7247, 16.7578, 16.1765, 14.9389, 14.2802, 13.6783, 12.9599, 12.2249, 11.6932, 11.022, 10.3871, 9.86386, 9.32372, 8.85135, 8.46894, 8.07039, 7.62505, 7.19262, 6.82553, 6.45724, 6.0716, 5.80577, 5.45966, 5.17205, 4.88161, 4.62344, 4.36325, 4.11154, 3.89572, 3.65409, 3.39633, 3.26805, 3.04014, 2.83118, 2.67305, 2.47943, 2.34147, 2.18737, 1.99697, 1.88241, 1.76099, 1.64482, 1.53348, 1.42174, 1.27572, 1.19665, 1.07725, 0.934853, 0.859419, 0.754135, 0.683542, 0.604478, 0.557281, 0.541549, 0.52138, 0.507664, 0.503227, 0.50121, 0.5, 0.5};
return Ncoll[hiBin];
};
const int nBins = 200;
const Double_t binTable[nBins+1] = {0, 12.0726, 12.9775, 13.7646, 14.6009, 15.4167, 16.2718, 17.1496, 18.0357, 18.933, 19.8565, 20.8485, 21.8614, 22.9305, 23.9966, 25.0743, 26.2168, 27.3867, 28.5948, 29.892, 31.1614, 32.5169, 33.96, 35.4747, 37.0106, 38.5893, 40.2316, 41.9823, 43.8428, 45.8105, 47.8386, 50.0019, 52.2309, 54.4736, 56.8414, 59.323, 61.9717, 64.7627, 67.6805, 70.6497, 73.7457, 76.9925, 80.1679, 83.6124, 87.1841, 90.8269, 94.685, 98.676, 102.906, 107.15, 111.736, 116.464, 121.338, 126.239, 131.41, 136.801, 142.322, 148.098, 154.006, 160.177, 166.526, 172.901, 179.823, 186.767, 194.153, 201.624, 209.203, 217.071, 225.281, 233.668, 242.615, 251.737, 261.417, 270.915, 280.69, 290.442, 301.228, 312.217, 322.842, 334.313, 346.203, 358.49, 370.721, 383.106, 396.469, 409.893, 423.605, 437.406, 451.505, 465.856, 481.029, 496.263, 512.126, 528.386, 544.806, 561.834, 579.189, 596.839, 614.239, 632.817, 651.443, 670.052, 689.384, 709.959, 729.154, 749.497, 770.253, 792.303, 813.733, 835.81, 858.891, 881.69, 906.947, 931.019, 955.307, 981.046, 1007.22, 1033.61, 1059.87, 1086.71, 1114.7, 1142.3, 1169.1, 1198.04, 1227.5, 1258.53, 1289.54, 1320.27, 1351.44, 1383.87, 1417.43, 1451.81, 1486.47, 1519.98, 1555.82, 1591.49, 1628.13, 1664.54, 1701.71, 1739.11, 1778.24, 1816.52, 1855.62, 1895.81, 1938.61, 1981.62, 2023.59, 2067.99, 2111.18, 2155.79, 2201.01, 2247.7, 2295.59, 2343.33, 2392.58, 2443.16, 2494.46, 2545.55, 2599.57, 2653.17, 2706.87, 2762.46, 2819.43, 2876.81, 2934.32, 2993.26, 3054.02, 3114.36, 3177.81, 3241.03, 3303.68, 3369.72, 3436.8, 3504.28, 3576.17, 3649.11, 3719.79, 3791.38, 3866.09, 3941.29, 4017.09, 4096.05, 4177.37, 4259.28, 4344.66, 4431.01, 4522.42, 4614.38, 4708.27, 4805.19, 4901.19, 5001.78, 5105.13, 5206.38, 5312.61, 5421.24, 5536.04, 5659.64, 5797.53, 5975.82, 7029.37};
Int_t getHiBinFromhiHF(const Double_t hiHF)
{
Int_t binPos=-1;
for(int i=0;i<nBins;++i){
if(hiHF>=binTable[i]&&hiHF<binTable[i+1]){
binPos=i;
break;
}
}
binPos=nBins-1-binPos;
return(Int_t)(200*((Double_t)binPos)/((Double_t)nBins));
}