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 static const int _trigEt   { 0 };
0007 static const int _geantid  { 0 };
0008 static const int _bbc      { 1 };
0009 static const int _zdcx     { 2 };
0010 static const int _pt       { 3 };
0011 static const int _dphi     { 4 }; // used for jets relative to trigger
0012 static const int _eta      { 4 };
0013 static const int _dphi     { 4 }; // used for jets relative to trigger
0014 static const int _dca      { 5 };
0015 static const int _ptmatch  { 5 }; // used for Aj jet pairs
0016 static const int _nhitfit  { 6 };
0017 static const int _Aj       { 6 }; // used for Aj jet pairs
0018 static const int _nhitrat  { 7 };
0019 static const int _ptreco   { 8 };
0020 
0021 // used for Aj and acoplanarity 
0022 /* static const int _dphiAj   { 6 }; */
0023 
0024 static const int _etaAu  { 0 }; // use like (*data)[_etaAu][_PUfit_m]
0025 static const int _etaMid { 1 };
0026 static const int _etaPP  { 2 };
0027 static const int _etaAll { 3 };
0028 
0029 static const int _phi_trig   { 1 };
0030 static const int _phi_trans  { 2 };
0031 static const int _phi_recoil { 3 };
0032 
0033 static const int _PUfit_b { 0 };
0034 static const int _PUfit_m { 1 };
0035 
0036 // bin boundaries for 6 bins of ZDCx (the 7th is for all ZDCx)
0037 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
0038 static const array<int, 9> _1_zdcx_3kHz { 3, 6, 9, 12, 14, 17, 20, 23, 23 };
0039 array<int,9> _i0_zdcx { 5, 8,  11, 14, 17, 20, 23, 26,  5 }; // binning by 3 for zdcx
0040 array<int,9> _i1_zdcx { 7, 10, 13, 16, 19, 22, 25, 28,  28};
0041 // bin kHz for the various bins (bin 1 at 4 kHz, etc...)
0042 static const double kHz_ZDCx (int bin) { return (float)bin+3; };
0043 
0044 
0045 #ifdef ioClass__h
0046     static ioBinVec bin_id      {{ -0.5,  -0.5,   6,      5.5     }};
0047     static ioBinVec bin_bbc     {{ 0.,  2802.38, 5460.32, 8420.66, 11722, 15410.1, 19584.9, 24371, 30105.3, 37666.3, 64000}};
0048     static ioBinVec bin_bbc3070 {{ 0.,  8420.66,  24371, 64000}};
0049     static ioBinVec bin_zdcx    {{ 0., 0.,  30,     30000.  }};
0050     static ioBinVec bin_pt   {{ // track pT bins                                               
0051          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,
0052      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,
0053      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,
0054      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,
0055      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
0056     }};
0057     static ioBinVec bin_trigEt  {{ 0.,    4.,     8.,    12.,  30. }};
0058     static ioBinVec bin_dca     {{ 0.,    0.,     15,     3.       }};
0059     static ioBinVec bin_hitsfit {{ 0.5,   0.5,    47.,    47.5     }};
0060     static ioBinVec bin_fitrat  {{ 0.02,  0.02,   10,     1.02     }};
0061     static ioBinVec bin_eta     {{ -1.,  -0.3,   0.3,     1.       }};
0062     static ioBinVec bin_dummy   {{ 0., 1. }}; // a dummy bin for non-used axes
0063 #endif
0064 
0065 #ifdef noiBinVec__h
0066     static tuBinVec bin_id      {{ -0.5,  -0.5,   6,      5.5     }};
0067     static tuBinVec bin_bbc     {{ 0.,  2802.38, 5460.32, 8420.66, 11722, 15410.1, 19584.9, 24371, 30105.3, 37666.3, 64000}};
0068     static tuBinVec bin_bbc3070 {{ 0.,  8420.66,  24371, 64000}};
0069     static tuBinVec bin_zdcx    {{ 0., 0.,  30,     30000.  }};
0070     static tuBinVec bin_pt   {{ // track pT bins                                               
0071          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,
0072      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,
0073      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,
0074      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,
0075      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
0076     }};
0077     static tuBinVec bin_trigEt  {{ 0.,    4.,     8.,    12.,  30. }};
0078     static tuBinVec bin_dca     {{ 0.,    0.,     15,     3.       }};
0079     static tuBinVec bin_hitsfit {{ 0.5,   0.5,    47.,    47.5     }};
0080     static tuBinVec bin_fitrat  {{ 0.02,  0.02,   10,     1.02     }};
0081     static tuBinVec bin_eta     {{ -1.,  -0.3,   0.3,     1.       }};
0082     static tuBinVec bin_dummy   {{ 0., 1. }}; // a dummy bin for non-used axes
0083 #endif
0084 
0085     static map<string, array<double,4>> TsallisParams {
0086         {"Tsallis_pp_pi",       {{ 0.135000,  5.5037062019,  0.1277746601,  9.7595242212 }}},
0087         {"Tsallis_pp_antipi",   {{ 0.135000,  5.5865021613,  0.1270067077,  9.7325864609 }}},
0088         {"Tsallis_pp_K",        {{ 0.493677,  0.1386983308,  0.1929685766, 11.8247896746 }}},
0089         {"Tsallis_pp_antiK",    {{ 0.493677,  0.1386983308,  0.1929685766, 11.8247896746 }}},
0090         {"Tsallis_pp_p",        {{ 0.938272,  0.0749981986,  0.1758814595, 10.5452494593 }}},
0091         {"Tsallis_pp_pbar",     {{ 0.938272,  0.0642725426,  0.1691919677, 10.0635969124 }}},
0092 
0093         {"Tsallis_dAu_pi",      {{ 0.135000, 13.1377700117,  0.1525501142, 10.1869046869 }}},
0094         {"Tsallis_dAu_antipi",  {{ 0.135000, 13.6251153259,  0.1505198390, 10.0949606259 }}},
0095         {"Tsallis_dAu_K",       {{ 0.493677,  0.4615505225,  0.2215729988, 11.6137583852 }}},
0096         {"Tsallis_dAu_antiK",   {{ 0.493677,  0.4675391506,  0.2192846303, 11.4475846094 }}},
0097         {"Tsallis_dAu_p",       {{ 0.938272,  0.2873036847,  0.2114254603, 11.4410990835 }}},
0098         {"Tsallis_dAu_pbar",    {{ 0.938272,  0.2257531137,  0.2195209347, 13.0260305553 }}}
0099     };
0100                                               //
0101     static constexpr double _n10xPyth6[9] {
0102         3920155 , 2101168 , 1187058 , 1695141 ,
0103         4967075 , 1797387 , 260676 , 261926 , 262366
0104     };
0105     static constexpr double _XsecPyth6[9] {
0106         0.107509,   0.0190967,  0.00475202,
0107         0.00198812, 0.000361282, 9.65463E-06,
0108         4.71077E-07, 2.68464E-08, 1.38211E-09
0109     };                                             //
0110     static constexpr double _pthat9[9] {
0111         0.107509    / 3920155,
0112         0.0190967   / 2101168,
0113         0.00475202  / 1187058,
0114         0.00198812  / 1695141,
0115         0.000361282 / 4967075,
0116         9.65463E-06 / 1797387,
0117         4.71077E-07 / 260676,
0118         2.68464E-08 / 261926,
0119         1.38211E-09 / 262366
0120     };
0121                                               //
0122 #endif // pAu_bins__h
0123        //
0124