Skip to content
Merged
Changes from 2 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
5 changes: 2 additions & 3 deletions MC/config/common/external/generator/performanceGenerator.C
Original file line number Diff line number Diff line change
Expand Up @@ -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)
Expand Down Expand Up @@ -261,8 +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 inv_pt = (gRandom->Rndm() / pTScale) * (kMaxInvPt - min_inv_pt) + min_inv_pt;
ffloat inv_pt = gRandom->Rndm() * (kBaseMaxInvPt - kBaseMinInvPt) + kBaseMinInvPt;
Comment thread
miranov25 marked this conversation as resolved.
Outdated
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
Expand Down