@@ -90,44 +90,40 @@ const std::vector<int> particleCharge{1, 1, 1, 1, 2, 2};
9090constexpr int NBetheParams = 6 ;
9191const std::vector<int > hfMothCodes = {511 , 521 , 531 , 541 , 5122 }; // b-mesons + Lambda_b
9292const std::vector<std::string> betheBlochParNames{" p0" , " p1" , " p2" , " p3" , " p4" , " resolution" };
93- const std::array<std::array<double , NBetheParams>, NParticles> kBetheBlochDefault {{
94- {13.611469 , 3.598765 , -0.021138 , 2.039562 , 0.651040 , 0.09 }, // pion
95- {5.393020 , 7.859534 , 0.004048 , 2.323197 , 1.609307 , 0.09 }, // proton
96- {5.393020 , 7.859534 , 0.004048 , 2.323197 , 1.609307 , 0.09 }, // deuteron
97- {5.393020 , 7.859534 , 0.004048 , 2.323197 , 1.609307 , 0.09 }, // triton
98- {-126.557359 , -0.858569 , 1.111643 , 1.210323 , 2.656374 , 0.09 }, // helion
99- {-126.557359 , -0.858569 , 1.111643 , 1.210323 , 2.656374 , 0.09 }}}; // alpha
93+ const std::array<std::array<double , NBetheParams>, NParticles> kBetheBlochDefault {{{13.611469 , 3.598765 , -0.021138 , 2.039562 , 0.651040 , 0.09 }, // pion
94+ {5.393020 , 7.859534 , 0.004048 , 2.323197 , 1.609307 , 0.09 }, // proton
95+ {5.393020 , 7.859534 , 0.004048 , 2.323197 , 1.609307 , 0.09 }, // deuteron
96+ {5.393020 , 7.859534 , 0.004048 , 2.323197 , 1.609307 , 0.09 }, // triton
97+ {-126.557359 , -0.858569 , 1.111643 , 1.210323 , 2.656374 , 0.09 }, // helion
98+ {-126.557359 , -0.858569 , 1.111643 , 1.210323 , 2.656374 , 0.09 }}}; // alpha
10099
101100// ---- Track PID/quality selection defaults (per species) --------------------
102101constexpr int NTrkSettings = 13 ;
103102const std::vector<std::string> trackPIDsettingsNames{" useBBparams" , " minITSnCls" , " minITSnClscos" , " minTPCnCls" , " maxTPCchi2" , " minTPCchi2" , " maxITSchi2" , " maxTPCnSigma" , " maxDcaXY" , " maxDcaZ" , " minITSclsSize" , " minTPCnClsCrossedRows" , " minReqClusterITSib" };
104- const std::array<std::array<double , NTrkSettings>, NParticles> kTrackPIDSettings {{
105- {0 , 0 , 4 , 60 , 4.0 , 0.5 , 100 , 2.5 , 2 ., 2 ., 0 ., 70 , 1 },
106- {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 3.0 , 2 ., 2 ., 0 ., 70 , 1 },
107- {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 3.0 , 2 ., 2 ., 0 ., 70 , 1 },
108- {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 3.0 , 2 ., 2 ., 0 ., 70 , 1 },
109- {1 , 0 , 4 , 75 , 4.0 , 0.5 , 100 , 5.0 , 2 ., 2 ., 0 ., 70 , 1 },
110- {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 5.0 , 2 ., 2 ., 0 ., 70 , 1 }}};
103+ const std::array<std::array<double , NTrkSettings>, NParticles> kTrackPIDSettings {{{0 , 0 , 4 , 60 , 4.0 , 0.5 , 100 , 2.5 , 2 ., 2 ., 0 ., 70 , 1 },
104+ {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 3.0 , 2 ., 2 ., 0 ., 70 , 1 },
105+ {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 3.0 , 2 ., 2 ., 0 ., 70 , 1 },
106+ {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 3.0 , 2 ., 2 ., 0 ., 70 , 1 },
107+ {1 , 0 , 4 , 75 , 4.0 , 0.5 , 100 , 5.0 , 2 ., 2 ., 0 ., 70 , 1 },
108+ {1 , 0 , 4 , 70 , 4.0 , 0.5 , 100 , 5.0 , 2 ., 2 ., 0 ., 70 , 1 }}};
111109
112110constexpr int NTrkSettings2 = 7 ;
113111const std::vector<std::string> trackPIDsettingsNames2{" useITSnsigma" , " minITSnsigma" , " maxITSnsigma" , " fillsparsh" , " useTPCnsigmaTOF" , " maxTPCnsigmaTOF" , " maxTOFnsigma" };
114- const std::array<std::array<double , NTrkSettings2>, NParticles> kTrackPIDSettings2 {{
115- {1 , -5 , 4 , 0 , 1 , 2 , 2 },
116- {1 , -5 , 4 , 0 , 1 , 2 , 2 },
117- {1 , -5 , 4 , 0 , 1 , 2 , 2 },
118- {1 , -5 , 4 , 1 , 1 , 2 , 2 },
119- {1 , -5 , 4 , 1 , 1 , 2 , 2 },
120- {1 , -5 , 4 , 1 , 1 , 2 , 2 }}};
112+ const std::array<std::array<double , NTrkSettings2>, NParticles> kTrackPIDSettings2 {{{1 , -5 , 4 , 0 , 1 , 2 , 2 },
113+ {1 , -5 , 4 , 0 , 1 , 2 , 2 },
114+ {1 , -5 , 4 , 0 , 1 , 2 , 2 },
115+ {1 , -5 , 4 , 1 , 1 , 2 , 2 },
116+ {1 , -5 , 4 , 1 , 1 , 2 , 2 },
117+ {1 , -5 , 4 , 1 , 1 , 2 , 2 }}};
121118
122119// ---- pT-dependent reconstruction correction parameters for He3/He4 ---------
123120constexpr int NFittingParticle = 4 ;
124121constexpr int NFittingParameters = 4 ;
125122const std::vector<std::string> trackcorrectionNames{" a" , " b" , " c" , " d" };
126- const std::array<std::array<double , NFittingParameters>, NFittingParticle> ktrackcorrection{{
127- {0.464215 , 0.195771 , 0.0183111 , 0.0 }, // He3
128- {0.464215 , 0.195771 , 0.0183111 , 0.0 }, // anti-He3
129- {0.00765 , 0.503791 , -1.10517 , 0.0 }, // He4
130- {0.00765 , 0.503791 , -1.10517 , 0.0 }}}; // anti-He4
123+ const std::array<std::array<double , NFittingParameters>, NFittingParticle> ktrackcorrection{{{0.464215 , 0.195771 , 0.0183111 , 0.0 }, // He3
124+ {0.464215 , 0.195771 , 0.0183111 , 0.0 }, // anti-He3
125+ {0.00765 , 0.503791 , -1.10517 , 0.0 }, // He4
126+ {0.00765 , 0.503791 , -1.10517 , 0.0 }}}; // anti-He4
131127
132128// ---- pT-dependent DCA cut parameters ----------------------------------------
133129// DCAxy(pt) = p0 * exp( p1 * pt) + p2
@@ -136,9 +132,8 @@ constexpr int NDCATypes = 2;
136132constexpr int NDCAParameters = 3 ;
137133const std::vector<std::string> dcaNames{" p0" , " p1" , " p2" };
138134const std::vector<std::string> dcaTypeNames{" DCAxy" , " DCAz" };
139- const std::array<std::array<double , NDCAParameters>, NDCATypes> kDCAcorrection {{
140- {0.0118 , -0.6889 , 0.0017 },
141- {0.1014 , 1.7512 , 0.0024 }}};
135+ const std::array<std::array<double , NDCAParameters>, NDCATypes> kDCAcorrection {{{0.0118 , -0.6889 , 0.0017 },
136+ {0.1014 , 1.7512 , 0.0024 }}};
142137
143138// ---- MC decay-origin classification -----------------------------------------
144139constexpr int DecayTypePrimary = 0 ;
0 commit comments