Skip to content

Commit 7b877c4

Browse files
author
miranov25
committed
O2-6235 - change pt range and bug fix
kMaxInvPt to 2 (pt>0.5) make the inv_pt independent of the pt scale
1 parent 0f6d137 commit 7b877c4

File tree

1 file changed

+3
-2
lines changed

1 file changed

+3
-2
lines changed

MC/config/common/external/generator/performanceGenerator.C

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -208,7 +208,7 @@ namespace o2
208208
};
209209

210210
// pT bounds: Max pT ~5 TeV (ALICE Pb-Pb energy)
211-
const float kMaxInvPt = 1.0f; // Min pT = 1 GeV
211+
const float kMaxInvPt = 2.f; // Min pT = 0.5 GeV
212212
const float kBaseMinInvPt = 2e-4f; // Max pT = 5000 GeV (unscaled)
213213

214214
// Check if particle is a parton (quark/gluon, status=11)
@@ -262,7 +262,8 @@ namespace o2
262262
int status = isParton(pdgCode) ? 11 : 1;
263263
// 4. Kinematics (flat 1/pT, max ~5000 GeV / pTScale)
264264
float min_inv_pt = kBaseMinInvPt / pTScale; // E.g., max pT=40,000 GeV for b quarks
265-
float inv_pt = (gRandom->Rndm() / pTScale) * (kMaxInvPt - min_inv_pt) + min_inv_pt;
265+
float max_inv_pt = kBaseMaxInvPt / pTScale; // E.g., max pT=40,000 GeV for b quarks
266+
float inv_pt = (gRandom->Rndm() / pTScale) * (max_inv_pt - min_inv_pt) + min_inv_pt;
266267
float pt = 1.0f / inv_pt;
267268
float phi = gRandom->Rndm() * 2.0f * TMath::Pi();
268269
float eta = gRandom->Rndm() * 3.0f - 1.5f; // ALICE TPC: -1.5 to 1.5

0 commit comments

Comments
 (0)