@@ -121,6 +121,7 @@ struct AntinucleiInJets {
121121 Configurable<double> maxNormalizedJetArea{"maxNormalizedJetArea", 1.0, "area cut"};
122122 Configurable<double> deltaEtaEdge{"deltaEtaEdge", 0.05, "eta gap from the edge"};
123123 Configurable<int> nSyst{"nSyst", 50, "number of systematic variations"};
124+ Configurable<int> nSubsamples{"nSubsamples", 50, "number of subsamples"};
124125
125126 // Track quality, kinematic, and PID selection parameters
126127 Configurable<bool> requirePvContributor{"requirePvContributor", false, "require that the track is a PV contributor"};
@@ -472,7 +473,7 @@ struct AntinucleiInJets {
472473 const AxisSpec nBarD2Axis{100, 0.0, 100.0, "N_{#bar{d}}^{i} #times N_{#bar{d}}^{j}"};
473474 const AxisSpec nBarP2Axis{100, 0.0, 100.0, "N_{#bar{p}}^{i} #times N_{#bar{p}}^{j}"};
474475 const AxisSpec nBarDnBarPAxis{100, 0.0, 100.0, "N_{#bar{d}}^{i} #times N_{#bar{p}}^{j}"};
475- const AxisSpec subsampleAxis{20 , 0, 20 , "Subsample Index"};
476+ const AxisSpec subsampleAxis{nSubsamples , 0, static_cast<double>(nSubsamples) , "Subsample Index"};
476477
477478 // Event counter
478479 registryCorr.add("eventCounter", "number of events", HistType::kTH1F, {{20, 0, 20, "counter"}});
@@ -2684,7 +2685,7 @@ struct AntinucleiInJets {
26842685 registryCorr.fill(HIST("eventCounter"), 7.5);
26852686
26862687 // Assign event to a random subsample (0-19)
2687- double sampleId = mRand.Integer(20 ) + 0.5;
2688+ double sampleId = mRand.Integer(nSubsamples ) + 0.5;
26882689
26892690 // Multiplicity percentile
26902691 const float multiplicity = collision.centFT0M();
0 commit comments