@@ -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 , 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