From e572aee1a5f94af75e75f3743f79cd984f695dc3 Mon Sep 17 00:00:00 2001 From: aimeric-landou Date: Sun, 26 Jan 2025 21:43:47 +0100 Subject: [PATCH 1/2] PWWGGAJE - add configuration for pp sim anchored to Pb-Pb @5.36TeV --- .../GeneratorJE_gapgen5_hook_pp5360GeV.ini | 5 ++++ .../generator/pythia8_jet_pp5360GeV.cfg | 25 +++++++++++++++++++ .../generator/pythia8_minbias_pp5360GeV.cfg | 18 +++++++++++++ 3 files changed, 48 insertions(+) create mode 100644 MC/config/PWGGAJE/ini/GeneratorJE_gapgen5_hook_pp5360GeV.ini create mode 100644 MC/config/PWGGAJE/pythia8/generator/pythia8_jet_pp5360GeV.cfg create mode 100644 MC/config/PWGGAJE/pythia8/generator/pythia8_minbias_pp5360GeV.cfg diff --git a/MC/config/PWGGAJE/ini/GeneratorJE_gapgen5_hook_pp5360GeV.ini b/MC/config/PWGGAJE/ini/GeneratorJE_gapgen5_hook_pp5360GeV.ini new file mode 100644 index 000000000..dc65d729d --- /dev/null +++ b/MC/config/PWGGAJE/ini/GeneratorJE_gapgen5_hook_pp5360GeV.ini @@ -0,0 +1,5 @@ +### jet-jet production with MB Gap 5 +### The external generator derives from GeneratorPythia8. +[GeneratorExternal] +fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGGAJE/external/generator/generator_pythia8_gaptrigger_jets_hook.C +funcName = getGeneratorPythia8GapGenJE(5,"${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGGAJE/pythia8/generator/pythia8_minbias_pp5360GeV.cfg","${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGGAJE/pythia8/generator/pythia8_jet_pp5360GeV.cfg") diff --git a/MC/config/PWGGAJE/pythia8/generator/pythia8_jet_pp5360GeV.cfg b/MC/config/PWGGAJE/pythia8/generator/pythia8_jet_pp5360GeV.cfg new file mode 100644 index 000000000..59f936556 --- /dev/null +++ b/MC/config/PWGGAJE/pythia8/generator/pythia8_jet_pp5360GeV.cfg @@ -0,0 +1,25 @@ + +# ### random +# # if we want to input a spectific seed: + +# Random:setSeed = on +# Random:seed = seed + +### beams +Beams:idA = 2212 # proton +Beams:idB = 2212 # proton +Beams:eCM = 5360. # GeV + +### processes +SoftQCD:inelastic = off +HardQCD:all = on + +### decays +ParticleDecays:limitTau0 = on +ParticleDecays:tau0Max = 10. + +### phase space cuts +PhaseSpace:pTHatMin = 5 +PhaseSpace:pTHatMax = 600 +PhaseSpace:bias2Selection = on +PhaseSpace:bias2SelectionPow = 4 diff --git a/MC/config/PWGGAJE/pythia8/generator/pythia8_minbias_pp5360GeV.cfg b/MC/config/PWGGAJE/pythia8/generator/pythia8_minbias_pp5360GeV.cfg new file mode 100644 index 000000000..c0397d1d6 --- /dev/null +++ b/MC/config/PWGGAJE/pythia8/generator/pythia8_minbias_pp5360GeV.cfg @@ -0,0 +1,18 @@ + +# ### random +# # if we want to input a spectific seed: + +# Random:setSeed = on +# Random:seed = seed + +### beams +Beams:idA = 2212 # proton +Beams:idB = 2212 # proton +Beams:eCM = 5360. # GeV + +### processes +SoftQCD:inelastic = on + +### decays +ParticleDecays:limitTau0 = on +ParticleDecays:tau0Max = 10. From 76d9dbd8c0cb95a51633c1b86730f153245327fe Mon Sep 17 00:00:00 2001 From: aimeric-landou Date: Sun, 26 Jan 2025 23:41:48 +0100 Subject: [PATCH 2/2] add test --- .../GeneratorJE_gapgen5_hook_pp5360GeV.C | 71 +++++++++++++++++++ 1 file changed, 71 insertions(+) create mode 100644 MC/config/PWGGAJE/ini/tests/GeneratorJE_gapgen5_hook_pp5360GeV.C diff --git a/MC/config/PWGGAJE/ini/tests/GeneratorJE_gapgen5_hook_pp5360GeV.C b/MC/config/PWGGAJE/ini/tests/GeneratorJE_gapgen5_hook_pp5360GeV.C new file mode 100644 index 000000000..5ac7f5173 --- /dev/null +++ b/MC/config/PWGGAJE/ini/tests/GeneratorJE_gapgen5_hook_pp5360GeV.C @@ -0,0 +1,71 @@ +int External() { + std::string path{"o2sim_Kine.root"}; + + float ratioTrigger = 1./5; // one event triggered out of 5 + + + TFile file(path.c_str(), "READ"); + if (file.IsZombie()) { + std::cerr << "Cannot open ROOT file " << path << "\n"; + return 1; + } + + auto tree = (TTree *)file.Get("o2sim"); + std::vector *tracks{}; + tree->SetBranchAddress("MCTrack", &tracks); + o2::dataformats::MCEventHeader *eventHeader = nullptr; + tree->SetBranchAddress("MCEventHeader.", &eventHeader); + + int nEventsMB{}, nEventsJetJet{}; + float sumWeightsMB{}, sumWeightsJetJet{}; + int sumTracks{}; + auto nEvents = tree->GetEntries(); + + for (int i = 0; i < nEvents; i++) { + tree->GetEntry(i); + + // check subgenerator information and event weights + if (eventHeader->hasInfo(o2::mcgenid::GeneratorProperty::SUBGENERATORID)) { + bool isValid = false; + int subGeneratorId = eventHeader->getInfo(o2::mcgenid::GeneratorProperty::SUBGENERATORID, isValid); + if (eventHeader->hasInfo(o2::dataformats::MCInfoKeys::weight)) { + float weight = eventHeader->getInfo(o2::dataformats::MCInfoKeys::weight,isValid); + if (subGeneratorId == 0) { + nEventsMB++; + sumWeightsMB += weight; + } + else if (subGeneratorId == 1) { + nEventsJetJet++; + sumWeightsJetJet += weight; + } + } + } + sumTracks += tracks->size(); + } + + std::cout << "--------------------------------\n"; + std::cout << "# Events: " << nEvents << "\n"; + std::cout << "# MB events: " << nEventsMB << "\n"; + std::cout << " sum of weights for MB events: " << sumWeightsMB << "\n"; + std::cout << "# Jet-jet events " << nEventsJetJet << "\n"; + std::cout << " sum of weights jet-jet events: " << sumWeightsJetJet << "\n"; + std::cout << "# tracks summed over all events (jet-jet + MB): " << sumTracks << "\n"; + + if (nEventsMB < nEvents * (1 - ratioTrigger) * 0.95 || nEventsMB > nEvents * (1 - ratioTrigger) * 1.05) { // we put some tolerance since the number of generated events is small + std::cerr << "Number of generated MB events different than expected\n"; + return 1; + } + if (nEventsJetJet < nEvents * ratioTrigger * 0.95 || nEventsJetJet > nEvents * ratioTrigger * 1.05) { + std::cerr << "Number of jet-jet generated events different than expected\n"; + return 1; + } + if(nEventsMB < sumWeightsMB * 0.95 || nEventsMB > sumWeightsMB * 1.05) { + std::cerr << "Weights of MB events do not = 1 as expected\n"; + return 1; + } + if(sumTracks < 1) { + std::cerr << "No tracks in simulated events\n"; + return 1; + } + return 0; +}