From 7b877c48454d6f727b9e3d680a0001e4c591a108 Mon Sep 17 00:00:00 2001 From: miranov25 Date: Fri, 13 Feb 2026 13:07:59 +0100 Subject: [PATCH 1/3] O2-6235 - change pt range and bug fix kMaxInvPt to 2 (pt>0.5) make the inv_pt independent of the pt scale --- MC/config/common/external/generator/performanceGenerator.C | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/MC/config/common/external/generator/performanceGenerator.C b/MC/config/common/external/generator/performanceGenerator.C index ce59cdfaa..b7deef243 100644 --- a/MC/config/common/external/generator/performanceGenerator.C +++ b/MC/config/common/external/generator/performanceGenerator.C @@ -208,7 +208,7 @@ namespace o2 }; // pT bounds: Max pT ~5 TeV (ALICE Pb-Pb energy) - const float kMaxInvPt = 1.0f; // Min pT = 1 GeV + const float kMaxInvPt = 2.f; // Min pT = 0.5 GeV const float kBaseMinInvPt = 2e-4f; // Max pT = 5000 GeV (unscaled) // Check if particle is a parton (quark/gluon, status=11) @@ -262,7 +262,8 @@ namespace o2 int status = isParton(pdgCode) ? 11 : 1; // 4. Kinematics (flat 1/pT, max ~5000 GeV / pTScale) float min_inv_pt = kBaseMinInvPt / pTScale; // E.g., max pT=40,000 GeV for b quarks - float inv_pt = (gRandom->Rndm() / pTScale) * (kMaxInvPt - min_inv_pt) + min_inv_pt; + float max_inv_pt = kBaseMaxInvPt / pTScale; // E.g., max pT=40,000 GeV for b quarks + float inv_pt = (gRandom->Rndm() / pTScale) * (max_inv_pt - min_inv_pt) + min_inv_pt; float pt = 1.0f / inv_pt; float phi = gRandom->Rndm() * 2.0f * TMath::Pi(); float eta = gRandom->Rndm() * 3.0f - 1.5f; // ALICE TPC: -1.5 to 1.5 From eb27358c6cb55ce9ccccc323d3278185eb9e3b96 Mon Sep 17 00:00:00 2001 From: miranov25 Date: Fri, 13 Feb 2026 13:53:41 +0100 Subject: [PATCH 2/3] Make range simple --- MC/config/common/external/generator/performanceGenerator.C | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/MC/config/common/external/generator/performanceGenerator.C b/MC/config/common/external/generator/performanceGenerator.C index b7deef243..a8fc1f583 100644 --- a/MC/config/common/external/generator/performanceGenerator.C +++ b/MC/config/common/external/generator/performanceGenerator.C @@ -261,9 +261,7 @@ namespace o2 // 3. Status: 11 for partons (jets), 1 for final-state int status = isParton(pdgCode) ? 11 : 1; // 4. Kinematics (flat 1/pT, max ~5000 GeV / pTScale) - float min_inv_pt = kBaseMinInvPt / pTScale; // E.g., max pT=40,000 GeV for b quarks - float max_inv_pt = kBaseMaxInvPt / pTScale; // E.g., max pT=40,000 GeV for b quarks - float inv_pt = (gRandom->Rndm() / pTScale) * (max_inv_pt - min_inv_pt) + min_inv_pt; + ffloat inv_pt = gRandom->Rndm() * (kBaseMaxInvPt - kBaseMinInvPt) + kBaseMinInvPt; float pt = 1.0f / inv_pt; float phi = gRandom->Rndm() * 2.0f * TMath::Pi(); float eta = gRandom->Rndm() * 3.0f - 1.5f; // ALICE TPC: -1.5 to 1.5 From 5075cb5a3a9c1b8ecff587aeb13572c0a3aef207 Mon Sep 17 00:00:00 2001 From: miranov25 Date: Mon, 16 Feb 2026 14:34:39 +0100 Subject: [PATCH 3/3] Correct names for ranges --- MC/config/common/external/generator/performanceGenerator.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/MC/config/common/external/generator/performanceGenerator.C b/MC/config/common/external/generator/performanceGenerator.C index a8fc1f583..d592d6f25 100644 --- a/MC/config/common/external/generator/performanceGenerator.C +++ b/MC/config/common/external/generator/performanceGenerator.C @@ -209,7 +209,7 @@ namespace o2 // pT bounds: Max pT ~5 TeV (ALICE Pb-Pb energy) const float kMaxInvPt = 2.f; // Min pT = 0.5 GeV - const float kBaseMinInvPt = 2e-4f; // Max pT = 5000 GeV (unscaled) + const float kMinInvPt = 2e-4f; // Max pT = 5000 GeV (unscaled) // Check if particle is a parton (quark/gluon, status=11) bool isParton(int& pdgCode) @@ -261,7 +261,7 @@ namespace o2 // 3. Status: 11 for partons (jets), 1 for final-state int status = isParton(pdgCode) ? 11 : 1; // 4. Kinematics (flat 1/pT, max ~5000 GeV / pTScale) - ffloat inv_pt = gRandom->Rndm() * (kBaseMaxInvPt - kBaseMinInvPt) + kBaseMinInvPt; + float inv_pt = gRandom->Rndm() * (kMaxInvPt - kMinInvPt) + kMinInvPt; float pt = 1.0f / inv_pt; float phi = gRandom->Rndm() * 2.0f * TMath::Pi(); float eta = gRandom->Rndm() * 3.0f - 1.5f; // ALICE TPC: -1.5 to 1.5