From e7f3d4f82aa07e877d5fb3d334271db0489e33ac Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Tue, 27 May 2025 10:23:30 +0200 Subject: [PATCH] PWGHF: Exclude parton event in Ds reso triggered config --- ..._D2H_ccbar_and_bbbar_gap5_DResoTrigger.ini | 2 +- ...HF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.C | 21 +------------------ 2 files changed, 2 insertions(+), 21 deletions(-) diff --git a/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.ini b/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.ini index 5391b6f4d..8c5f26dd1 100644 --- a/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.ini +++ b/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 diff --git a/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.C b/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.C index 0fda5cf6b..008ef62f8 100644 --- a/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.C +++ b/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";