From 35efc45e7ebb46e37d7664c8cf47c5738915e7fc Mon Sep 17 00:00:00 2001 From: sawan Date: Thu, 8 Oct 2026 18:44:10 +0530 Subject: [PATCH] made the mass and pt cut configurable --- .../GeneratorLF_doublephi_triggerMasspTcut.ini | 2 +- .../generator_pythia8_twophi_triggerMassCut.C | 17 +++++++++++------ 2 files changed, 12 insertions(+), 7 deletions(-) diff --git a/MC/config/PWGLF/ini/GeneratorLF_doublephi_triggerMasspTcut.ini b/MC/config/PWGLF/ini/GeneratorLF_doublephi_triggerMasspTcut.ini index a51496b44..753c01d31 100644 --- a/MC/config/PWGLF/ini/GeneratorLF_doublephi_triggerMasspTcut.ini +++ b/MC/config/PWGLF/ini/GeneratorLF_doublephi_triggerMasspTcut.ini @@ -1,6 +1,6 @@ [GeneratorExternal] fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGLF/pythia8/generator_pythia8_twophi_triggerMassCut.C -funcName=generateDoublePhi(0, 0.0, 100.0, 0.8) +funcName=generateDoublePhi(0, 0.0, 100.0, 0.8, 2.4, 6.0) [GeneratorPythia8] config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGLF/pythia8/generator/pythia8_inel_136tev.cfg diff --git a/MC/config/PWGLF/pythia8/generator_pythia8_twophi_triggerMassCut.C b/MC/config/PWGLF/pythia8/generator_pythia8_twophi_triggerMassCut.C index 4cb177bc3..4851840c8 100644 --- a/MC/config/PWGLF/pythia8/generator_pythia8_twophi_triggerMassCut.C +++ b/MC/config/PWGLF/pythia8/generator_pythia8_twophi_triggerMassCut.C @@ -25,14 +25,17 @@ class GeneratorPythia8DoublePhi : public o2::eventgen::GeneratorPythia8 { public: /// Constructor - GeneratorPythia8DoublePhi(int gapSize = 0, double minPt = 0.0, double maxPt = 100.0, double maxEta = 0.8) + GeneratorPythia8DoublePhi(int gapSize = 0, double minPt = 0.0, double maxPt = 100.0, double maxEta = 0.8, double minPairMass = 2.4, double minPairPt = 6.0) : o2::eventgen::GeneratorPythia8(), mGapSize(gapSize), mMinPt(minPt), mMaxPt(maxPt), - mMaxEta(maxEta) + mMaxEta(maxEta), + mMinPairMass(minPairMass), + mMinPairPt(minPairPt) { - fmt::printf(">> Pythia8 generator: two phi(1020) mesons, gap = %d, minPtPhi = %f, maxPtPhi = %f, |etaPhi| < %f\n", gapSize, minPt, maxPt, maxEta); + fmt::printf(">> Pythia8 generator: two phi(1020) mesons, gap = %d, minPtPhi = %f, maxPtPhi = %f, |etaPhi| < %f, minPairMass = %f, minPairPt = %f\n", + gapSize, minPt, maxPt, maxEta, minPairMass, minPairPt); } /// Destructor ~GeneratorPythia8DoublePhi() = default; @@ -156,7 +159,7 @@ protected: double mass = phiPhi.M(); double pt = phiPhi.Pt(); - if (mass > 2.4 && pt > 4.0) + if (mass > mMinPairMass && pt > mMinPairPt) return true; } } @@ -169,13 +172,15 @@ private: double mMinPt{0.0}; double mMaxPt{100.0}; double mMaxEta{0.8}; + double mMinPairMass{2.4}; + double mMinPairPt{6.0}; uint64_t mGeneratedEvents{0}; }; ///___________________________________________________________ -FairGenerator *generateDoublePhi(int gap = 0, double minPt = 0.0, double maxPt = 100.0, double maxEta = 0.8) +FairGenerator *generateDoublePhi(int gap = 0, double minPt = 0.0, double maxPt = 100.0, double maxEta = 0.8, double minPairMass = 2.4, double minPairPt = 6.0) { - auto myGenerator = new GeneratorPythia8DoublePhi(gap, minPt, maxPt, maxEta); + auto myGenerator = new GeneratorPythia8DoublePhi(gap, minPt, maxPt, maxEta, minPairMass, minPairPt); myGenerator->readString("333:onMode = off"); myGenerator->readString("333:onIfMatch = 321 -321");