|
42 | 42 |
|
43 | 43 | #include <fmt/format.h> |
44 | 44 |
|
45 | | -#include <algorithm> |
46 | 45 | #include <array> |
47 | 46 | #include <cmath> |
48 | 47 | #include <cstdlib> |
@@ -192,7 +191,6 @@ struct FoxWolframCorrelation { |
192 | 191 | Configurable<float> vtxRange{"vtxRange", 10.0f, "Vertex Z range to consider"}; |
193 | 192 | Configurable<float> etaRange{"etaRange", 0.8f, "eta range to consider"}; |
194 | 193 | Configurable<float> dcaZ{"dcaZ", 0.2f, "custom DCA Z cut (ignored if negative)"}; |
195 | | - Configurable<bool> doMultipSampling{"doMultipSampling", false, "doMultipSampling"}; |
196 | 194 | Configurable<bool> doSameEvent{"doSameEvent", true, "Fill same-event TRI correlations"}; |
197 | 195 | Configurable<bool> doFWMCorrelations{"doFWMCorrelations", false, |
198 | 196 | "Fill H_i-H_j correlation histograms for RANGE1-RANGE3"}; |
@@ -248,45 +246,6 @@ struct FoxWolframCorrelation { |
248 | 246 |
|
249 | 247 | Configurable<int> noOfMultipBins{"noOfMultipBins", 10, "noOfMultipBins"}; |
250 | 248 |
|
251 | | - // Provide reference multiplicity weights for optional sampling. |
252 | | - std::array<int, 500> limMultip{ |
253 | | - 7058, 43307, 85345, 110727, 118857, 118652, 115912, 112011, 109254, 106044, |
254 | | - 102421, 100599, 97952, 95908, 93883, 91890, 90264, 88097, 86388, 84582, |
255 | | - 82542, 81212, 79240, 77903, 75562, 74099, 71973, 70777, 69100, 67880, |
256 | | - 66162, 64827, 63513, 62101, 61110, 59307, 58419, 57719, 56344, 55228, |
257 | | - 54615, 53847, 52387, 51925, 51560, 50217, 49962, 48766, 48005, 47651, |
258 | | - 46758, 46046, 45531, 45063, 44271, 43530, 43151, 42489, 42347, 41576, |
259 | | - 40991, 40166, 39959, 39494, 39367, 38830, 38225, 37781, 37435, 36934, |
260 | | - 36501, 35980, 35846, 35312, 34926, 34640, 33693, 33721, 33296, 32570, |
261 | | - 32272, 32122, 31690, 31390, 31125, 30532, 29855, 29623, 29486, 29032, |
262 | | - 28534, 28140, 27491, 27430, 26760, 26601, 26217, 25625, 25464, 25067, |
263 | | - 24395, 24271, 23722, 23240, 22953, 22474, 22166, 21621, 21312, 20935, |
264 | | - 20444, 20157, 19714, 19077, 18774, 18284, 18005, 17391, 16963, 16864, |
265 | | - 16184, 15722, 15456, 14956, 14386, 14248, 13853, 13419, 13132, 12602, |
266 | | - 12192, 11667, 11397, 10967, 10701, 10316, 9831, 9758, 9371, 8905, |
267 | | - 8759, 8255, 7992, 7603, 7363, 7150, 6816, 6477, 6168, 5940, |
268 | | - 5638, 5502, 5154, 5014, 4686, 4434, 4263, 3998, 3959, 3623, |
269 | | - 3479, 3323, 3196, 3005, 2923, 2770, 2607, 2354, 2226, 2187, |
270 | | - 2058, 1980, 1794, 1739, 1601, 1515, 1419, 1321, 1292, 1175, |
271 | | - 1144, 1048, 995, 909, 860, 783, 783, 745, 680, 571, |
272 | | - 566, 524, 533, 462, 444, 427, 367, 320, 303, 317, |
273 | | - 270, 250, 227, 221, 207, 193, 188, 163, 131, 151, |
274 | | - 139, 104, 105, 70, 94, 83, 68, 62, 85, 70, |
275 | | - 46, 38, 43, 47, 39, 33, 38, 23, 28, 27, |
276 | | - 24, 17, 24, 17, 14, 27, 22, 14, 11, 11, |
277 | | - 10, 11, 7, 6, 5, 6, 6, 7, 2, 5, |
278 | | - 6, 4, 2, 5, 3, 2, 3, 1, 0, 0, |
279 | | - 2, 0, 1, 3, 2, 0, 1, 3, 1, 1, |
280 | | - 2, 0, 1, 1, 0, 0, 0, 1, 2, 0, |
281 | | - 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, |
282 | | - 0, 0, 1, 0, 0, 0, 3, 0, 0, 0, |
283 | | - 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, |
284 | | - 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, |
285 | | - 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, |
286 | | - 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, |
287 | | - 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}; |
288 | | - std::array<int, 500> countMultip{}; |
289 | | - |
290 | 249 | static constexpr std::array<std::string_view, 10> binNames = { |
291 | 250 | "BIN1", "BIN2", "BIN3", "BIN4", "BIN5", |
292 | 251 | "BIN6", "BIN7", "BIN8", "BIN9", "BIN10"}; |
@@ -777,16 +736,6 @@ struct FoxWolframCorrelation { |
777 | 736 | (topologyMode.value == kTriEvents && isTRI) || |
778 | 737 | (topologyMode.value == kIsoEvents && isISO); |
779 | 738 |
|
780 | | - // Optionally cap the sampled events at each multiplicity. |
781 | | - int bin = std::clamp<int>(nGlobalTracks, 0, 499); |
782 | | - if (doMultipSampling) { |
783 | | - if (countMultip[bin] >= limMultip[bin]) { |
784 | | - return; |
785 | | - } |
786 | | - |
787 | | - ++countMultip[bin]; |
788 | | - } |
789 | | - |
790 | 739 | // Fill inclusive event multiplicity and recoil. |
791 | 740 | histos.fill(HIST("Multip/ITSTPC_Multiplicity"), nGlobalTracks); |
792 | 741 | histos.fill(HIST("pTSpectra/RecoilALL"), nGlobalTracks, nRecoil); |
|
0 commit comments