Skip to content

Commit 0f553c0

Browse files
prchakraprchakra
authored
[PWGCF] FemtoUniverse: Fixing issues in upper limits of Cent. (#8333)
1 parent a8bbad7 commit 0f553c0

File tree

1 file changed

+32
-32
lines changed

1 file changed

+32
-32
lines changed

‎PWGCF/FemtoUniverse/Core/FemtoUniversePairWithCentMultKt.h‎

Lines changed: 32 additions & 32 deletions
Original file line numberDiff line numberDiff line change
@@ -105,57 +105,57 @@ class PairWithCentMultKt
105105
void fill(t1 kstar_value, t1 cent_mult_value, t1 kt_value)
106106
{
107107

108-
if (cent_mult_value > CentMultBins[CentMultBins.size() - 1] || cent_mult_value < CentMultBins[0]) {
108+
if (cent_mult_value >= CentMultBins[CentMultBins.size() - 1] || cent_mult_value < CentMultBins[0]) {
109109
PairWithCentMultKtRegistry->fill(HIST("Beyond_Max"), kstar_value);
110-
} else if (cent_mult_value <= CentMultBins[1]) {
110+
} else if (cent_mult_value < CentMultBins[1]) {
111111
PairWithCentMultKtRegistry->fill(HIST("mult_0_1/kstar"), kstar_value);
112112
if (UseKt) {
113113
auto histMultFolder = HIST("mult_0_1/");
114114
fill_kT(kstar_value, kt_value, histMultFolder);
115115
}
116-
} else if (cent_mult_value <= CentMultBins[2]) {
116+
} else if (cent_mult_value < CentMultBins[2]) {
117117
PairWithCentMultKtRegistry->fill(HIST("mult_1_2/kstar"), kstar_value);
118118
if (UseKt) {
119119
auto histMultFolder = HIST("mult_1_2/");
120120
fill_kT(kstar_value, kt_value, histMultFolder);
121121
}
122-
} else if (cent_mult_value <= CentMultBins[3]) {
122+
} else if (cent_mult_value < CentMultBins[3]) {
123123
PairWithCentMultKtRegistry->fill(HIST("mult_2_3/kstar"), kstar_value);
124124
if (UseKt) {
125125
auto histMultFolder = HIST("mult_2_3/");
126126
fill_kT(kstar_value, kt_value, histMultFolder);
127127
}
128-
} else if (cent_mult_value <= CentMultBins[4]) {
128+
} else if (cent_mult_value < CentMultBins[4]) {
129129
PairWithCentMultKtRegistry->fill(HIST("mult_3_4/kstar"), kstar_value);
130130
if (UseKt) {
131131
auto histMultFolder = HIST("mult_3_4/");
132132
fill_kT(kstar_value, kt_value, histMultFolder);
133133
}
134-
} else if (cent_mult_value <= CentMultBins[5]) {
134+
} else if (cent_mult_value < CentMultBins[5]) {
135135
PairWithCentMultKtRegistry->fill(HIST("mult_4_5/kstar"), kstar_value);
136136
if (UseKt) {
137137
auto histMultFolder = HIST("mult_4_5/");
138138
fill_kT(kstar_value, kt_value, histMultFolder);
139139
}
140-
} else if (cent_mult_value <= CentMultBins[6]) {
140+
} else if (cent_mult_value < CentMultBins[6]) {
141141
PairWithCentMultKtRegistry->fill(HIST("mult_5_6/kstar"), kstar_value);
142142
if (UseKt) {
143143
auto histMultFolder = HIST("mult_5_6/");
144144
fill_kT(kstar_value, kt_value, histMultFolder);
145145
}
146-
} else if (cent_mult_value <= CentMultBins[7]) {
146+
} else if (cent_mult_value < CentMultBins[7]) {
147147
PairWithCentMultKtRegistry->fill(HIST("mult_6_7/kstar"), kstar_value);
148148
if (UseKt) {
149149
auto histMultFolder = HIST("mult_6_7/");
150150
fill_kT(kstar_value, kt_value, histMultFolder);
151151
}
152-
} else if (cent_mult_value <= CentMultBins[8]) {
152+
} else if (cent_mult_value < CentMultBins[8]) {
153153
PairWithCentMultKtRegistry->fill(HIST("mult_7_8/kstar"), kstar_value);
154154
if (UseKt) {
155155
auto histMultFolder = HIST("mult_7_8/");
156156
fill_kT(kstar_value, kt_value, histMultFolder);
157157
}
158-
} else if (cent_mult_value <= CentMultBins[9]) {
158+
} else if (cent_mult_value < CentMultBins[9]) {
159159
PairWithCentMultKtRegistry->fill(HIST("mult_8_9/kstar"), kstar_value);
160160
if (UseKt) {
161161
auto histMultFolder = HIST("mult_8_9/");
@@ -173,23 +173,23 @@ class PairWithCentMultKt
173173
template <typename t1, typename t2>
174174
void fill_kT(t1 kstar_value, t1 kt_value, t2 folder)
175175
{
176-
if (kt_value <= KtBins[1]) {
176+
if (kt_value < KtBins[1]) {
177177
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_0_1"), kstar_value);
178-
} else if (kt_value <= KtBins[2]) {
178+
} else if (kt_value < KtBins[2]) {
179179
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_1_2"), kstar_value);
180-
} else if (kt_value <= KtBins[3]) {
180+
} else if (kt_value < KtBins[3]) {
181181
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_2_3"), kstar_value);
182-
} else if (kt_value <= KtBins[4]) {
182+
} else if (kt_value < KtBins[4]) {
183183
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_3_4"), kstar_value);
184-
} else if (kt_value <= KtBins[5]) {
184+
} else if (kt_value < KtBins[5]) {
185185
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_4_5"), kstar_value);
186-
} else if (kt_value <= KtBins[6]) {
186+
} else if (kt_value < KtBins[6]) {
187187
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_5_6"), kstar_value);
188-
} else if (kt_value <= KtBins[7]) {
188+
} else if (kt_value < KtBins[7]) {
189189
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_6_7"), kstar_value);
190-
} else if (kt_value <= KtBins[8]) {
190+
} else if (kt_value < KtBins[8]) {
191191
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_7_8"), kstar_value);
192-
} else if (kt_value <= KtBins[9]) {
192+
} else if (kt_value < KtBins[9]) {
193193
PairWithCentMultKtRegistry->fill(folder + HIST("kstar_kt_8_9"), kstar_value);
194194
}
195195
}
@@ -204,57 +204,57 @@ class PairWithCentMultKt
204204
void fill_3D(t1 qout_value, t1 qside_value, t1 qlong_value, t1 cent_mult_value, t1 kt_value)
205205
{
206206

207-
if (cent_mult_value > CentMultBins[CentMultBins.size() - 1] || cent_mult_value < CentMultBins[0]) {
207+
if (cent_mult_value >= CentMultBins[CentMultBins.size() - 1] || cent_mult_value < CentMultBins[0]) {
208208
PairWithCentMultKtRegistry->fill(HIST("Beyond_Max_3D"), qout_value, qside_value, qlong_value);
209-
} else if (cent_mult_value <= CentMultBins[1]) {
209+
} else if (cent_mult_value < CentMultBins[1]) {
210210
PairWithCentMultKtRegistry->fill(HIST("mult_0_1/q3D"), qout_value, qside_value, qlong_value);
211211
if (Use3D) {
212212
auto histMultFolder = HIST("mult_0_1/");
213213
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
214214
}
215-
} else if (cent_mult_value <= CentMultBins[2]) {
215+
} else if (cent_mult_value < CentMultBins[2]) {
216216
// PairWithCentMultKtRegistry->fill(HIST("mult_1_2/q3D"), qout_value, qside_value, qlong_value);
217217
if (Use3D) {
218218
auto histMultFolder = HIST("mult_1_2/");
219219
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
220220
}
221-
} else if (cent_mult_value <= CentMultBins[3]) {
221+
} else if (cent_mult_value < CentMultBins[3]) {
222222
// PairWithCentMultKtRegistry->fill(HIST("mult_2_3/q3D"), qout_value, qside_value, qlong_value);
223223
if (Use3D) {
224224
auto histMultFolder = HIST("mult_2_3/");
225225
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
226226
}
227-
} else if (cent_mult_value <= CentMultBins[4]) {
227+
} else if (cent_mult_value < CentMultBins[4]) {
228228
// PairWithCentMultKtRegistry->fill(HIST("mult_3_4/q3D"), qout_value, qside_value, qlong_value);
229229
if (Use3D) {
230230
auto histMultFolder = HIST("mult_3_4/");
231231
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
232232
}
233-
} else if (cent_mult_value <= CentMultBins[5]) {
233+
} else if (cent_mult_value < CentMultBins[5]) {
234234
// PairWithCentMultKtRegistry->fill(HIST("mult_4_5/q3D"), qout_value, qside_value, qlong_value);
235235
if (Use3D) {
236236
auto histMultFolder = HIST("mult_4_5/");
237237
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
238238
}
239-
} else if (cent_mult_value <= CentMultBins[6]) {
239+
} else if (cent_mult_value < CentMultBins[6]) {
240240
// PairWithCentMultKtRegistry->fill(HIST("mult_5_6/q3D"), qout_value, qside_value, qlong_value);
241241
if (Use3D) {
242242
auto histMultFolder = HIST("mult_5_6/");
243243
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
244244
}
245-
} else if (cent_mult_value <= CentMultBins[7]) {
245+
} else if (cent_mult_value < CentMultBins[7]) {
246246
// PairWithCentMultKtRegistry->fill(HIST("mult_6_7/q3D"), qout_value, qside_value, qlong_value);
247247
if (Use3D) {
248248
auto histMultFolder = HIST("mult_6_7/");
249249
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
250250
}
251-
} else if (cent_mult_value <= CentMultBins[8]) {
251+
} else if (cent_mult_value < CentMultBins[8]) {
252252
// PairWithCentMultKtRegistry->fill(HIST("mult_7_8/q3D"), qout_value, qside_value, qlong_value);
253253
if (Use3D) {
254254
auto histMultFolder = HIST("mult_7_8/");
255255
fill_kT_3d(qout_value, qside_value, qlong_value, kt_value, histMultFolder);
256256
}
257-
} else if (cent_mult_value <= CentMultBins[9]) {
257+
} else if (cent_mult_value < CentMultBins[9]) {
258258
// PairWithCentMultKtRegistry->fill(HIST("mult_8_9/q3D"), qout_value, qside_value, qlong_value);
259259
if (Use3D) {
260260
auto histMultFolder = HIST("mult_8_9/");
@@ -273,11 +273,11 @@ class PairWithCentMultKt
273273
template <typename t1, typename t2>
274274
void fill_kT_3d(t1 qout_value, t1 qside_value, t1 qlong_value, t1 kt_value, t2 folder)
275275
{
276-
if (kt_value <= KtBins[1]) {
276+
if (kt_value < KtBins[1]) {
277277
PairWithCentMultKtRegistry->fill(folder + HIST("q3D_kt_0_1"), qout_value, qside_value, qlong_value);
278-
} else if (kt_value <= KtBins[2]) {
278+
} else if (kt_value < KtBins[2]) {
279279
PairWithCentMultKtRegistry->fill(folder + HIST("q3D_kt_1_2"), qout_value, qside_value, qlong_value);
280-
} else if (kt_value <= KtBins[3]) {
280+
} else if (kt_value < KtBins[3]) {
281281
PairWithCentMultKtRegistry->fill(folder + HIST("q3D_kt_2_3"), qout_value, qside_value, qlong_value);
282282
}
283283
}

0 commit comments

Comments
 (0)