@@ -172,7 +172,7 @@ class FemtoUniversePairWithCentMultKt
172172 template <typename t1, typename t2>
173173 void fillkT (t1 kstar_value, t1 kt_value, t2 folder)
174174 {
175- if (kt_value < ktBins[1 ]) {
175+ if (kt_value >= ktBins[ 0 ] && kt_value < ktBins[1 ]) {
176176 pairWithCentMultKtRegistry->fill (folder + HIST (" kstar_kt_0_1" ), kstar_value);
177177 } else if (kt_value < ktBins[2 ]) {
178178 pairWithCentMultKtRegistry->fill (folder + HIST (" kstar_kt_1_2" ), kstar_value);
@@ -272,12 +272,18 @@ class FemtoUniversePairWithCentMultKt
272272 template <typename t1, typename t2>
273273 void fillkT3D (t1 qout_value, t1 qside_value, t1 qlong_value, t1 kt_value, t2 folder)
274274 {
275- if (kt_value < ktBins[1 ]) {
275+ if (kt_value >= ktBins[ 0 ] && kt_value < ktBins[1 ]) {
276276 pairWithCentMultKtRegistry->fill (folder + HIST (" q3D_kt_0_1" ), qout_value, qside_value, qlong_value);
277277 } else if (kt_value < ktBins[2 ]) {
278278 pairWithCentMultKtRegistry->fill (folder + HIST (" q3D_kt_1_2" ), qout_value, qside_value, qlong_value);
279279 } else if (kt_value < ktBins[3 ]) {
280280 pairWithCentMultKtRegistry->fill (folder + HIST (" q3D_kt_2_3" ), qout_value, qside_value, qlong_value);
281+ } else if (kt_value < ktBins[4 ]) {
282+ pairWithCentMultKtRegistry->fill (folder + HIST (" q3D_kt_3_4" ), qout_value, qside_value, qlong_value);
283+ } else if (kt_value < ktBins[5 ]) {
284+ pairWithCentMultKtRegistry->fill (folder + HIST (" q3D_kt_4_5" ), qout_value, qside_value, qlong_value);
285+ } else if (kt_value < ktBins[6 ]) {
286+ pairWithCentMultKtRegistry->fill (folder + HIST (" q3D_kt_5_6" ), qout_value, qside_value, qlong_value);
281287 }
282288 }
283289
0 commit comments