Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-06 08:14:12

0001 #ifndef pAu_bins__h
0002 #define pAu_bins__h
0003 
0004 #include "noiBinVec.h"
0005 
0006 const int ITER_BASE = 6;
0007 const int ITER_LO   = 4;
0008 const int ITER_HI   = 8;
0009 
0010 static const int _trigEt   { 0 };
0011 static const int _geantid  { 0 };
0012 static const int _bbc      { 1 };
0013 static const int _zdcx     { 2 };
0014 static const int _pt       { 3 };
0015 static const int _dphi     { 4 }; // used for jets relative to trigger
0016 static const int _eta      { 4 };
0017 static const int _dca      { 5 };
0018 static const int _ptmatch  { 5 }; // used for Aj jet pairs
0019 static const int _nhitfit  { 6 };
0020 static const int _Aj       { 6 }; // used for Aj jet pairs
0021 static const int _nhitrat  { 7 };
0022 static const int _ptreco   { 8 };
0023 
0024 static const int _rec_absdphi_event { 3 };
0025 static const int _rec_eta_event     { 4 };
0026 static const int _lead_eta_event    { 5 };
0027 static const int _lead_deta_event   { 6 };
0028 
0029 static const int _rec_absdphi_track { 8 };
0030 static const int _rec_eta_track     { 9 };
0031 static const int _lead_eta_track    { 10 };
0032 static const int _lead_deta_track   { 11 };
0033 
0034 // used for Aj and acoplanarity 
0035 /* static const int _dphiAj   { 6 }; */
0036 
0037 static const int _etaAu  { 0 }; // use like (*data)[_etaAu][_PUfit_m]
0038 static const int _etaMid { 1 };
0039 static const int _etaPP  { 2 };
0040 static const int _etaAll { 3 };
0041 
0042 static const int _phi_trig   { 1 };
0043 static const int _phi_trans  { 2 };
0044 static const int _phi_recoil { 3 };
0045 
0046 static const int _PUfit_b { 0 };
0047 static const int _PUfit_m { 1 };
0048 
0049 // bin boundaries for 6 bins of ZDCx (the 7th is for all ZDCx)
0050 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
0051 static const array<int, 9> _1_zdcx_3kHz { 3, 6, 9, 12, 14, 17, 20, 23, 23 };
0052 static const array<int,9> _i0_zdcx { 6, 9,  12, 15, 18, 21, 24, 27,  6 }; // binning by 3 for zdcx
0053 static const array<int,9> _i1_zdcx { 8, 11, 14, 17, 20, 23, 26, 29,  29};
0054 
0055 static array<int,5> _50_zdcx { 6, 9,  12, 15, 18 }; // binning by 3 for zdcx
0056 static array<int,5> _51_zdcx { 8, 11, 14, 17, 20 };
0057 
0058 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]" };
0059 // bin kHz for the various bins (bin 1 at 4 kHz, etc...)
0060 static const double kHz_ZDCx (int bin) { return (float)bin+3; };
0061 
0062 
0063 #ifdef ioClass__h
0064     static ioBinVec bin_id      {{ -0.5,  -0.5,   6,      5.5     }};
0065     static ioBinVec bin_bbc     {{ 0.,  2802.38, 5460.32, 8420.66, 11722, 15410.1, 19584.9, 24371, 30105.3, 37666.3, 64000}};
0066     static ioBinVec bin_bbc3070 {{ 0.,  8420.66,  24371, 64000}};
0067     static ioBinVec bin_zdcx    {{ 0., 0.,  30,     30000.  }};
0068 
0069     static ioBinVec bin_pt_v0   {{ // track pT bins                                               
0070          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,
0071      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,
0072      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,
0073      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,
0074      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
0075     }};
0076 
0077     static ioBinVec bin_pt {{ // track pT bins
0078      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,
0079      1.3,  1.4,  1.5, 1.6, 1.7, 1.8, 1.9,  2.0,  2.1,  2.2,  2.3,  2.4,
0080      2.5,  2.6,  2.7, 2.8, 2.9, 3.0, 3.1,  3.2,  3.3,  3.4,  3.5,  3.6,
0081      3.7,  3.8,  3.9, 4.0, 4.1, 4.2, 4.3,  4.4,  4.5,  4.7,  4.9,  5.1,
0082      5.3,  5.5,  5.7, 5.9, 6.1, 6.3, 6.5,  6.7,  6.9,  7.2,  7.5,  7.8,
0083      8.1,  8.4,  8.7, 9.1, 9.5, 9.9, 10.3, 10.8, 11.3, 11.9, 12.5, 13.2,
0084      14.0, 15.0 }}; // 74 bins
0085 
0086     static ioBinVec bin_trigEt  {{ 0.,    4.,     8.,    12.,  30. }};
0087     static ioBinVec bin_dca     {{ 0.,    0.,     15,     3.       }};
0088     static ioBinVec bin_hitsfit {{ 0.5,   0.5,    47.,    47.5     }};
0089     static ioBinVec bin_fitrat  {{ 0.02,  0.02,   10,     1.02     }};
0090     static ioBinVec bin_eta     {{ -1.,  -0.3,   0.3,     1.       }};
0091     static ioBinVec bin_eta_p9  {{ -0.9,  -0.3,   0.3,     0.9     }};
0092     static ioBinVec bin_eta_p20 {{ -1., -1., 20, 1. }};
0093     static ioBinVec bin_dummy   {{ 0., 1. }}; // a dummy bin for non-used axes
0094 #endif
0095 
0096 #ifdef noiBinVec__h
0097     static tuBinVec bin_id      {{ -0.5,  -0.5,   6,      5.5     }};
0098     static tuBinVec bin_bbc     {{ 0.,  2802.38, 5460.32, 8420.66, 11722, 15410.1, 19584.9, 24371, 30105.3, 37666.3, 64000}};
0099     static tuBinVec bin_bbc3070 {{ 0.,  8420.66,  24371, 64000}};
0100     static tuBinVec bin_zdcx    {{ 0., 0.,  30,     30000.  }};
0101 
0102     static tuBinVec bin_pt {{ // track pT bins
0103      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,
0104      1.3,  1.4,  1.5, 1.6, 1.7, 1.8, 1.9,  2.0,  2.1,  2.2,  2.3,  2.4,
0105      2.5,  2.6,  2.7, 2.8, 2.9, 3.0, 3.1,  3.2,  3.3,  3.4,  3.5,  3.6,
0106      3.7,  3.8,  3.9, 4.0, 4.1, 4.2, 4.3,  4.4,  4.5,  4.7,  4.9,  5.1,
0107      5.3,  5.5,  5.7, 5.9, 6.1, 6.3, 6.5,  6.7,  6.9,  7.2,  7.5,  7.8,
0108      8.1,  8.4,  8.7, 9.1, 9.5, 9.9, 10.3, 10.8, 11.3, 11.9, 12.5, 13.2,
0109      14.0, 15.0 }};
0110 
0111     static tuBinVec bin_pt_v0   {{ // track pT bins                                               
0112          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,
0113      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,
0114      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,
0115      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,
0116      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
0117     }};
0118 
0119     static tuBinVec bin_trigEt  {{ 0.,    4.,     8.,    12.,  30. }};
0120     static tuBinVec bin_dca     {{ 0.,    0.,     15,     3.       }};
0121     static tuBinVec bin_hitsfit {{ 0.5,   0.5,    47.,    47.5     }};
0122     static tuBinVec bin_fitrat  {{ 0.02,  0.02,   10,     1.02     }};
0123     static tuBinVec bin_eta     {{ -1.,  -0.3,   0.3,     1.       }};
0124     static tuBinVec bin_eta_p9  {{ -0.9,  -0.3,   0.3,     0.9     }};
0125     static tuBinVec bin_eta_p20 {{ -1., -1., 20, 1. }};
0126     static tuBinVec bin_dummy   {{ 0., 1. }}; // a dummy bin for non-used axes
0127 #endif
0128 
0129     static map<string, array<double,4>> TsallisParams {
0130         {"Tsallis_pp_pi",       {{ 0.135000,  5.5037062019,  0.1277746601,  9.7595242212 }}},
0131         {"Tsallis_pp_antipi",   {{ 0.135000,  5.5865021613,  0.1270067077,  9.7325864609 }}},
0132         {"Tsallis_pp_K",        {{ 0.493677,  0.1386983308,  0.1929685766, 11.8247896746 }}},
0133         {"Tsallis_pp_antiK",    {{ 0.493677,  0.1386983308,  0.1929685766, 11.8247896746 }}},
0134         {"Tsallis_pp_p",        {{ 0.938272,  0.0749981986,  0.1758814595, 10.5452494593 }}},
0135         {"Tsallis_pp_pbar",     {{ 0.938272,  0.0642725426,  0.1691919677, 10.0635969124 }}},
0136 
0137         {"Tsallis_dAu_pi",      {{ 0.135000, 13.1377700117,  0.1525501142, 10.1869046869 }}},
0138         {"Tsallis_dAu_antipi",  {{ 0.135000, 13.6251153259,  0.1505198390, 10.0949606259 }}},
0139         {"Tsallis_dAu_K",       {{ 0.493677,  0.4615505225,  0.2215729988, 11.6137583852 }}},
0140         {"Tsallis_dAu_antiK",   {{ 0.493677,  0.4675391506,  0.2192846303, 11.4475846094 }}},
0141         {"Tsallis_dAu_p",       {{ 0.938272,  0.2873036847,  0.2114254603, 11.4410990835 }}},
0142         {"Tsallis_dAu_pbar",    {{ 0.938272,  0.2257531137,  0.2195209347, 13.0260305553 }}}
0143     };
0144                                               //
0145     static constexpr double _n10xPyth6[9] {
0146         3920155 , 2101168 , 1187058 , 1695141 ,
0147         4967075 , 1797387 , 260676 , 261926 , 262366
0148     };
0149     static constexpr double _XsecPyth6[9] {
0150         0.107509,   0.0190967,  0.00475202,
0151         0.00198812, 0.000361282, 9.65463E-06,
0152         4.71077E-07, 2.68464E-08, 1.38211E-09
0153     };                                             //
0154     static constexpr double _pthat9[9] {
0155         0.107509    / 3920155,
0156         0.0190967   / 2101168,
0157         0.00475202  / 1187058,
0158         0.00198812  / 1695141,
0159         0.000361282 / 4967075,
0160         9.65463E-06 / 1797387,
0161         4.71077E-07 / 260676,
0162         2.68464E-08 / 261926,
0163         1.38211E-09 / 262366
0164     };
0165                                               //
0166 #endif // pAu_bins__h
0167