PWGHF: Exclude parton event in Ds reso triggered config

Author: fgrosa

MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.ini

@@ -5,4 +5,4 @@ funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(5, -1.5, 1.5, -1.5, 1.5, {10
 
 [GeneratorPythia8]
 config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_DResoTrigger.cfg
-includePartonEvent=true
+includePartonEvent=false # not needed for jet studies, hence no need to keep parton event

MC/config/PWGHF/ini/tests/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.C

@@ -37,7 +37,7 @@ int External() {
     tree->SetBranchAddress("MCEventHeader.", &eventHeader);
 
     int nEventsMB{}, nEventsInjOne{}, nEventsInjTwo{};
-    int nQuarksOne{}, nQuarksTwo{}, nSignals{}, nSignalGoodDecay{};
+    int nSignals{}, nSignalGoodDecay{};
     auto nEvents = tree->GetEntries();
 
     for (int i = 0; i < nEvents; i++) {
@@ -60,14 +60,6 @@ int External() {
         for (auto &track : *tracks) {
             auto pdg = track.GetPdgCode();
             auto absPdg = std::abs(pdg);
-            if (absPdg == checkPdgQuarkOne) {
-                nQuarksOne++;
-                continue;
-            }
-            if (absPdg == checkPdgQuarkTwo) {
-                nQuarksTwo++;
-                continue;
-            }
             if (std::find(checkPdgHadron.begin(), checkPdgHadron.end(), absPdg) != checkPdgHadron.end()) { // found signal
                 nSignals++; // count signal PDG
 
@@ -115,8 +107,6 @@ int External() {
     std::cout << "# MB events: " << nEventsMB << "\n";
     std::cout << Form("# events injected with %d quark pair: ", checkPdgQuarkOne) << nEventsInjOne << "\n";
     std::cout << Form("# events injected with %d quark pair: ", checkPdgQuarkTwo) << nEventsInjTwo << "\n";
-    std::cout << Form("# %d (anti)quarks: ", checkPdgQuarkOne) << nQuarksOne << "\n";
-    std::cout << Form("# %d (anti)quarks: ", checkPdgQuarkTwo) << nQuarksTwo << "\n";
     std::cout <<"# signal hadrons: " << nSignals << "\n";
     std::cout <<"# signal hadrons decaying in the correct channel: " << nSignalGoodDecay << "\n";
 
@@ -133,15 +123,6 @@ int External() {
         return 1;
     }
 
-    if (nQuarksOne < nEvents * ratioTrigger) { // we expect anyway more because the same quark is repeated several time, after each gluon radiation
-        std::cerr << "Number of generated (anti)quarks " << checkPdgQuarkOne << " lower than expected\n";
-        return 1;
-    }
-    if (nQuarksTwo < nEvents * ratioTrigger) { // we expect anyway more because the same quark is repeated several time, after each gluon radiation
-        std::cerr << "Number of generated (anti)quarks " << checkPdgQuarkTwo << " lower than expected\n";
-        return 1;
-    }
-
     float fracForcedDecays = float(nSignalGoodDecay) / nSignals;
     if (fracForcedDecays < 0.9) { // we put some tolerance (e.g. due to oscillations which might change the final state)
         std::cerr << "Fraction of signals decaying into the correct channel " << fracForcedDecays << " lower than expected\n";