Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
pAu_bins.h
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pAu_bins.h
1 #ifndef pAu_bins__h
2 #define pAu_bins__h
3 
4 #include "noiBinVec.h"
5 
6 const int ITER_BASE = 6;
7 const int ITER_LO = 4;
8 const int ITER_HI = 8;
9 
10 static const int _trigEt { 0 };
11 static const int _geantid { 0 };
12 static const int _bbc { 1 };
13 static const int _zdcx { 2 };
14 static const int _pt { 3 };
15 static const int _dphi { 4 }; // used for jets relative to trigger
16 static const int _eta { 4 };
17 static const int _dca { 5 };
18 static const int _ptmatch { 5 }; // used for Aj jet pairs
19 static const int _nhitfit { 6 };
20 static const int _Aj { 6 }; // used for Aj jet pairs
21 static const int _nhitrat { 7 };
22 static const int _ptreco { 8 };
23 
24 static const int _rec_absdphi_event { 3 };
25 static const int _rec_eta_event { 4 };
26 static const int _lead_eta_event { 5 };
27 static const int _lead_deta_event { 6 };
28 
29 static const int _rec_absdphi_track { 8 };
30 static const int _rec_eta_track { 9 };
31 static const int _lead_eta_track { 10 };
32 static const int _lead_deta_track { 11 };
33 
34 // used for Aj and acoplanarity
35 /* static const int _dphiAj { 6 }; */
36 
37 static const int _etaAu { 0 }; // use like (*data)[_etaAu][_PUfit_m]
38 static const int _etaMid { 1 };
39 static const int _etaPP { 2 };
40 static const int _etaAll { 3 };
41 
42 static const int _phi_trig { 1 };
43 static const int _phi_trans { 2 };
44 static const int _phi_recoil { 3 };
45 
46 static const int _PUfit_b { 0 };
47 static const int _PUfit_m { 1 };
48 
49 // bin boundaries for 6 bins of ZDCx (the 7th is for all ZDCx)
50 static const array<int, 9> _0_zdcx_3kHz { 1, 4, 7, 10, 12, 15, 18, 21, 1 }; // for binning in 3kHz increments 4-28 kHz
51 static const array<int, 9> _1_zdcx_3kHz { 3, 6, 9, 12, 14, 17, 20, 23, 23 };
52 static const array<int,9> _i0_zdcx { 6, 9, 12, 15, 18, 21, 24, 27, 6 }; // binning by 3 for zdcx
53 static const array<int,9> _i1_zdcx { 8, 11, 14, 17, 20, 23, 26, 29, 29};
54 
55 static array<int,5> _50_zdcx { 6, 9, 12, 15, 18 }; // binning by 3 for zdcx
56 static array<int,5> _51_zdcx { 8, 11, 14, 17, 20 };
57 
58 static array<string,3> name_eta_range { "#eta_{Au}#in[-1.,-0.3]", "#eta_{mid}#in[-0.3,0.3]", "#eta_{#it{p}}#in[0.3,1]" };
59 // bin kHz for the various bins (bin 1 at 4 kHz, etc...)
60 static const double kHz_ZDCx (int bin) { return (float)bin+3; };
61 
62 
63 #ifdef ioClass__h
64  static ioBinVec bin_id {{ -0.5, -0.5, 6, 5.5 }};
65  static ioBinVec bin_bbc {{ 0., 2802.38, 5460.32, 8420.66, 11722, 15410.1, 19584.9, 24371, 30105.3, 37666.3, 64000}};
66  static ioBinVec bin_bbc3070 {{ 0., 8420.66, 24371, 64000}};
67  static ioBinVec bin_zdcx {{ 0., 0., 30, 30000. }};
68 
69  static ioBinVec bin_pt_v0 {{ // track pT bins
70  0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4,
71  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,
72  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,
73  4.5, 4.7, 4.9, 5.1, 5.3, 5.5, 5.7, 5.9, 6.1, 6.3, 6.5, 6.7, 6.9, 7.2, 7.5,
74  7.8, 8.1, 8.4, 8.7, 9.1, 9.5, 9.9, 10.4, 10.9, 11.4, 12.0, 12.6, 13.3, 15.0
75  }};
76 
77  static ioBinVec bin_pt {{ // track pT bins
78  0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2,
79  1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4,
80  2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.5, 3.6,
81  3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 4.7, 4.9, 5.1,
82  5.3, 5.5, 5.7, 5.9, 6.1, 6.3, 6.5, 6.7, 6.9, 7.2, 7.5, 7.8,
83  8.1, 8.4, 8.7, 9.1, 9.5, 9.9, 10.3, 10.8, 11.3, 11.9, 12.5, 13.2,
84  14.0, 15.0 }}; // 74 bins
85 
86  static ioBinVec bin_trigEt {{ 0., 4., 8., 12., 30. }};
87  static ioBinVec bin_dca {{ 0., 0., 15, 3. }};
88  static ioBinVec bin_hitsfit {{ 0.5, 0.5, 47., 47.5 }};
89  static ioBinVec bin_fitrat {{ 0.02, 0.02, 10, 1.02 }};
90  static ioBinVec bin_eta {{ -1., -0.3, 0.3, 1. }};
91  static ioBinVec bin_eta_p9 {{ -0.9, -0.3, 0.3, 0.9 }};
92  static ioBinVec bin_eta_p20 {{ -1., -1., 20, 1. }};
93  static ioBinVec bin_dummy {{ 0., 1. }}; // a dummy bin for non-used axes
94 #endif
95 
96 #ifdef noiBinVec__h
97  static tuBinVec bin_id {{ -0.5, -0.5, 6, 5.5 }};
98  static tuBinVec bin_bbc {{ 0., 2802.38, 5460.32, 8420.66, 11722, 15410.1, 19584.9, 24371, 30105.3, 37666.3, 64000}};
99  static tuBinVec bin_bbc3070 {{ 0., 8420.66, 24371, 64000}};
100  static tuBinVec bin_zdcx {{ 0., 0., 30, 30000. }};
101 
102  static tuBinVec bin_pt {{ // track pT bins
103  0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2,
104  1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4,
105  2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.5, 3.6,
106  3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 4.7, 4.9, 5.1,
107  5.3, 5.5, 5.7, 5.9, 6.1, 6.3, 6.5, 6.7, 6.9, 7.2, 7.5, 7.8,
108  8.1, 8.4, 8.7, 9.1, 9.5, 9.9, 10.3, 10.8, 11.3, 11.9, 12.5, 13.2,
109  14.0, 15.0 }};
110 
111  static tuBinVec bin_pt_v0 {{ // track pT bins
112  0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4,
113  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,
114  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,
115  4.5, 4.7, 4.9, 5.1, 5.3, 5.5, 5.7, 5.9, 6.1, 6.3, 6.5, 6.7, 6.9, 7.2, 7.5,
116  7.8, 8.1, 8.4, 8.7, 9.1, 9.5, 9.9, 10.4, 10.9, 11.4, 12.0, 12.6, 13.3, 15.0
117  }};
118 
119  static tuBinVec bin_trigEt {{ 0., 4., 8., 12., 30. }};
120  static tuBinVec bin_dca {{ 0., 0., 15, 3. }};
121  static tuBinVec bin_hitsfit {{ 0.5, 0.5, 47., 47.5 }};
122  static tuBinVec bin_fitrat {{ 0.02, 0.02, 10, 1.02 }};
123  static tuBinVec bin_eta {{ -1., -0.3, 0.3, 1. }};
124  static tuBinVec bin_eta_p9 {{ -0.9, -0.3, 0.3, 0.9 }};
125  static tuBinVec bin_eta_p20 {{ -1., -1., 20, 1. }};
126  static tuBinVec bin_dummy {{ 0., 1. }}; // a dummy bin for non-used axes
127 #endif
128 
129  static map<string, array<double,4>> TsallisParams {
130  {"Tsallis_pp_pi", {{ 0.135000, 5.5037062019, 0.1277746601, 9.7595242212 }}},
131  {"Tsallis_pp_antipi", {{ 0.135000, 5.5865021613, 0.1270067077, 9.7325864609 }}},
132  {"Tsallis_pp_K", {{ 0.493677, 0.1386983308, 0.1929685766, 11.8247896746 }}},
133  {"Tsallis_pp_antiK", {{ 0.493677, 0.1386983308, 0.1929685766, 11.8247896746 }}},
134  {"Tsallis_pp_p", {{ 0.938272, 0.0749981986, 0.1758814595, 10.5452494593 }}},
135  {"Tsallis_pp_pbar", {{ 0.938272, 0.0642725426, 0.1691919677, 10.0635969124 }}},
136 
137  {"Tsallis_dAu_pi", {{ 0.135000, 13.1377700117, 0.1525501142, 10.1869046869 }}},
138  {"Tsallis_dAu_antipi", {{ 0.135000, 13.6251153259, 0.1505198390, 10.0949606259 }}},
139  {"Tsallis_dAu_K", {{ 0.493677, 0.4615505225, 0.2215729988, 11.6137583852 }}},
140  {"Tsallis_dAu_antiK", {{ 0.493677, 0.4675391506, 0.2192846303, 11.4475846094 }}},
141  {"Tsallis_dAu_p", {{ 0.938272, 0.2873036847, 0.2114254603, 11.4410990835 }}},
142  {"Tsallis_dAu_pbar", {{ 0.938272, 0.2257531137, 0.2195209347, 13.0260305553 }}}
143  };
144  //
145  static constexpr double _n10xPyth6[9] {
146  3920155 , 2101168 , 1187058 , 1695141 ,
147  4967075 , 1797387 , 260676 , 261926 , 262366
148  };
149  static constexpr double _XsecPyth6[9] {
150  0.107509, 0.0190967, 0.00475202,
151  0.00198812, 0.000361282, 9.65463E-06,
152  4.71077E-07, 2.68464E-08, 1.38211E-09
153  }; //
154  static constexpr double _pthat9[9] {
155  0.107509 / 3920155,
156  0.0190967 / 2101168,
157  0.00475202 / 1187058,
158  0.00198812 / 1695141,
159  0.000361282 / 4967075,
160  9.65463E-06 / 1797387,
161  4.71077E-07 / 260676,
162  2.68464E-08 / 261926,
163  1.38211E-09 / 262366
164  };
165  //
166 #endif // pAu_bins__h
167