Fix alternate injection for more than two hadrons

atriolo · gconesab · commit 3a712a8274ae · 2023-11-17T09:48:23.000+01:00
diff --git a/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C b/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C
@@ -96,7 +96,7 @@ protected:
       // Alternate hadrons if enabled (with the same ratio)
       if (mHadronPdgList.size() >= 1)
       {
-        int iHadron = (nInjectedEvents / std::max(mHadronPdgList.size(), 1ul)) % mHadronPdgList.size();
+        int iHadron = (nInjectedEvents / std::max(mQuarkPdgList.size(), 1ul)) % mHadronPdgList.size();
         mHadronPdg = mHadronPdgList[iHadron];
         LOG(debug)<<"SELECTED hadron: "<<mHadronPdgList[iHadron];
       }

Original file line number	Diff line number	Diff line change
`@@ -96,7 +96,7 @@ protected:`
`96`	`96`	`// Alternate hadrons if enabled (with the same ratio)`
`97`	`97`	`if (mHadronPdgList.size() >= 1)`
`98`	`98`	`{`
`99`		`- int iHadron = (nInjectedEvents / std::max(mHadronPdgList.size(), 1ul)) % mHadronPdgList.size();`
	`99`	`+ int iHadron = (nInjectedEvents / std::max(mQuarkPdgList.size(), 1ul)) % mHadronPdgList.size();`
`100`	`100`	`mHadronPdg = mHadronPdgList[iHadron];`
`101`	`101`	`LOG(debug)<<"SELECTED hadron: "<<mHadronPdgList[iHadron];`
`102`	`102`	`}`