Skip to content

Commit 0ae1a11

Browse files
prchakraprchakra
authored
[PWGCF] FemtoUniverse: Quick fix for centrality bining in storing CFs (#9178)
1 parent 8aad3e6 commit 0ae1a11

File tree

2 files changed

+21
-21
lines changed

2 files changed

+21
-21
lines changed

‎PWGCF/FemtoUniverse/Core/FemtoUniversePairSHCentMultKt.h‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -167,11 +167,11 @@ class PairSHCentMultKt
167167

168168
if ((absmultval >= centMultBins[0]) && (absmultval < centMultBins[1])) {
169169
multbinval = 0;
170-
} else if (absmultval < centMultBins[2]) {
170+
} else if ((absmultval >= centMultBins[1]) && (absmultval < centMultBins[2])) {
171171
multbinval = 1;
172-
} else if (absmultval < centMultBins[3]) {
172+
} else if ((absmultval >= centMultBins[2]) && (absmultval < centMultBins[3])) {
173173
multbinval = 2;
174-
} else if (absmultval < centMultBins[4]) {
174+
} else if ((absmultval >= centMultBins[3]) && (absmultval < centMultBins[4])) {
175175
multbinval = 3;
176176
} else {
177177
return;

‎PWGCF/FemtoUniverse/Core/FemtoUniversePairWithCentMultKt.h‎

Lines changed: 18 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -106,55 +106,55 @@ class FemtoUniversePairWithCentMultKt
106106

107107
if (cent_mult_value >= centMultBins[centMultBins.size() - 1] || cent_mult_value < centMultBins[0]) {
108108
pairWithCentMultKtRegistry->fill(HIST("Beyond_Max"), kstar_value);
109-
} else if (cent_mult_value < centMultBins[1]) {
109+
} else if (cent_mult_value >= centMultBins[0] && cent_mult_value < centMultBins[1]) {
110110
pairWithCentMultKtRegistry->fill(HIST("mult_0_1/kstar"), kstar_value);
111111
if (useKt) {
112112
auto histMultFolder = HIST("mult_0_1/");
113113
fillkT(kstar_value, kt_value, histMultFolder);
114114
}
115-
} else if (cent_mult_value < centMultBins[2]) {
115+
} else if (cent_mult_value >= centMultBins[1] && cent_mult_value < centMultBins[2]) {
116116
pairWithCentMultKtRegistry->fill(HIST("mult_1_2/kstar"), kstar_value);
117117
if (useKt) {
118118
auto histMultFolder = HIST("mult_1_2/");
119119
fillkT(kstar_value, kt_value, histMultFolder);
120120
}
121-
} else if (cent_mult_value < centMultBins[3]) {
121+
} else if (cent_mult_value >= centMultBins[2] && cent_mult_value < centMultBins[3]) {
122122
pairWithCentMultKtRegistry->fill(HIST("mult_2_3/kstar"), kstar_value);
123123
if (useKt) {
124124
auto histMultFolder = HIST("mult_2_3/");
125125
fillkT(kstar_value, kt_value, histMultFolder);
126126
}
127-
} else if (cent_mult_value < centMultBins[4]) {
127+
} else if (cent_mult_value >= centMultBins[3] && cent_mult_value < centMultBins[4]) {
128128
pairWithCentMultKtRegistry->fill(HIST("mult_3_4/kstar"), kstar_value);
129129
if (useKt) {
130130
auto histMultFolder = HIST("mult_3_4/");
131131
fillkT(kstar_value, kt_value, histMultFolder);
132132
}
133-
} else if (cent_mult_value < centMultBins[5]) {
133+
} else if (cent_mult_value >= centMultBins[4] && cent_mult_value < centMultBins[5]) {
134134
pairWithCentMultKtRegistry->fill(HIST("mult_4_5/kstar"), kstar_value);
135135
if (useKt) {
136136
auto histMultFolder = HIST("mult_4_5/");
137137
fillkT(kstar_value, kt_value, histMultFolder);
138138
}
139-
} else if (cent_mult_value < centMultBins[6]) {
139+
} else if (cent_mult_value >= centMultBins[5] && cent_mult_value < centMultBins[6]) {
140140
pairWithCentMultKtRegistry->fill(HIST("mult_5_6/kstar"), kstar_value);
141141
if (useKt) {
142142
auto histMultFolder = HIST("mult_5_6/");
143143
fillkT(kstar_value, kt_value, histMultFolder);
144144
}
145-
} else if (cent_mult_value < centMultBins[7]) {
145+
} else if (cent_mult_value >= centMultBins[6] && cent_mult_value < centMultBins[7]) {
146146
pairWithCentMultKtRegistry->fill(HIST("mult_6_7/kstar"), kstar_value);
147147
if (useKt) {
148148
auto histMultFolder = HIST("mult_6_7/");
149149
fillkT(kstar_value, kt_value, histMultFolder);
150150
}
151-
} else if (cent_mult_value < centMultBins[8]) {
151+
} else if (cent_mult_value >= centMultBins[7] && cent_mult_value < centMultBins[8]) {
152152
pairWithCentMultKtRegistry->fill(HIST("mult_7_8/kstar"), kstar_value);
153153
if (useKt) {
154154
auto histMultFolder = HIST("mult_7_8/");
155155
fillkT(kstar_value, kt_value, histMultFolder);
156156
}
157-
} else if (cent_mult_value < centMultBins[9]) {
157+
} else if (cent_mult_value >= centMultBins[8] && cent_mult_value < centMultBins[9]) {
158158
pairWithCentMultKtRegistry->fill(HIST("mult_8_9/kstar"), kstar_value);
159159
if (useKt) {
160160
auto histMultFolder = HIST("mult_8_9/");
@@ -205,55 +205,55 @@ class FemtoUniversePairWithCentMultKt
205205

206206
if (cent_mult_value >= centMultBins[centMultBins.size() - 1] || cent_mult_value < centMultBins[0]) {
207207
pairWithCentMultKtRegistry->fill(HIST("Beyond_Max_3D"), qout_value, qside_value, qlong_value);
208-
} else if (cent_mult_value < centMultBins[1]) {
208+
} else if (cent_mult_value >= centMultBins[0] && cent_mult_value < centMultBins[1]) {
209209
pairWithCentMultKtRegistry->fill(HIST("mult_0_1/q3D"), qout_value, qside_value, qlong_value);
210210
if (use3D) {
211211
auto histMultFolder = HIST("mult_0_1/");
212212
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
213213
}
214-
} else if (cent_mult_value < centMultBins[2]) {
214+
} else if (cent_mult_value >= centMultBins[1] && cent_mult_value < centMultBins[2]) {
215215
// pairWithCentMultKtRegistry->fill(HIST("mult_1_2/q3D"), qout_value, qside_value, qlong_value);
216216
if (use3D) {
217217
auto histMultFolder = HIST("mult_1_2/");
218218
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
219219
}
220-
} else if (cent_mult_value < centMultBins[3]) {
220+
} else if (cent_mult_value >= centMultBins[2] && cent_mult_value < centMultBins[3]) {
221221
// pairWithCentMultKtRegistry->fill(HIST("mult_2_3/q3D"), qout_value, qside_value, qlong_value);
222222
if (use3D) {
223223
auto histMultFolder = HIST("mult_2_3/");
224224
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
225225
}
226-
} else if (cent_mult_value < centMultBins[4]) {
226+
} else if (cent_mult_value >= centMultBins[3] && cent_mult_value < centMultBins[4]) {
227227
// pairWithCentMultKtRegistry->fill(HIST("mult_3_4/q3D"), qout_value, qside_value, qlong_value);
228228
if (use3D) {
229229
auto histMultFolder = HIST("mult_3_4/");
230230
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
231231
}
232-
} else if (cent_mult_value < centMultBins[5]) {
232+
} else if (cent_mult_value >= centMultBins[4] && cent_mult_value < centMultBins[5]) {
233233
// pairWithCentMultKtRegistry->fill(HIST("mult_4_5/q3D"), qout_value, qside_value, qlong_value);
234234
if (use3D) {
235235
auto histMultFolder = HIST("mult_4_5/");
236236
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
237237
}
238-
} else if (cent_mult_value < centMultBins[6]) {
238+
} else if (cent_mult_value >= centMultBins[5] && cent_mult_value < centMultBins[6]) {
239239
// pairWithCentMultKtRegistry->fill(HIST("mult_5_6/q3D"), qout_value, qside_value, qlong_value);
240240
if (use3D) {
241241
auto histMultFolder = HIST("mult_5_6/");
242242
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
243243
}
244-
} else if (cent_mult_value < centMultBins[7]) {
244+
} else if (cent_mult_value >= centMultBins[6] && cent_mult_value < centMultBins[7]) {
245245
// pairWithCentMultKtRegistry->fill(HIST("mult_6_7/q3D"), qout_value, qside_value, qlong_value);
246246
if (use3D) {
247247
auto histMultFolder = HIST("mult_6_7/");
248248
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
249249
}
250-
} else if (cent_mult_value < centMultBins[8]) {
250+
} else if (cent_mult_value >= centMultBins[7] && cent_mult_value < centMultBins[8]) {
251251
// pairWithCentMultKtRegistry->fill(HIST("mult_7_8/q3D"), qout_value, qside_value, qlong_value);
252252
if (use3D) {
253253
auto histMultFolder = HIST("mult_7_8/");
254254
fillkT3D(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
255255
}
256-
} else if (cent_mult_value < centMultBins[9]) {
256+
} else if (cent_mult_value >= centMultBins[8] && cent_mult_value < centMultBins[9]) {
257257
// pairWithCentMultKtRegistry->fill(HIST("mult_8_9/q3D"), qout_value, qside_value, qlong_value);
258258
if (use3D) {
259259
auto histMultFolder = HIST("mult_8_9/");

0 commit comments

Comments
 (0)