Skip to content

[PWGJE] Minimum pT hat cut feature in jetSpectraEseTask #18372

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
joachimckh wants to merge 3 commits into
base: master
Choose a base branch
from
+39 −1
Open

[PWGJE] Minimum pT hat cut feature in jetSpectraEseTask #18372

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
054984c
minimum pT hat feature in jetSpectraEseTask
joachimckh Oct 20, 2026
f0ef073
remove comment
joachimckh Oct 20, 2026
9d566d1
safe against 0/inf weights
joachimckh Oct 21, 2026
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
40 changes: 39 additions & 1 deletion PWGJE/Tasks/jetSpectraEseTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -92,6 +92,7 @@ struct JetSpectraEseTask {
Configurable<float> pTHatMaxMCD{"pTHatMaxMCD", 999.0, "maximum fraction of hard scattering for jet acceptance in detector MC"};
Configurable<float> pTHatMaxMCP{"pTHatMaxMCP", 999.0, "maximum fraction of hard scattering for jet acceptance in particle MC"};
Configurable<float> pTHatExponent{"pTHatExponent", 6.0, "exponent of the event weight for the calculation of pTHat"};
Configurable<float> pTHatAbsoluteMin{"pTHatAbsoluteMin", -99.0f, "minimum value of pTHat"};

Configurable<bool> fLeadJetPtCut{"fLeadJetPtCut", false, "Flag for leading jet pT cut"};
Configurable<float> leadJetPtMin{"leadJetPtMin", 20.0, "minimum leading jet pT cut"};
Expand Down Expand Up @@ -540,9 +541,13 @@ struct JetSpectraEseTask {
registry.add("h3CentdeltapTRndmConePhi_localrhovsphi", "centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho}; #Delta#varphi_{jet}", {HistType::kTH3F, {{100, 0.0, 100.0}, {400, -200.0, 200.0}, {dPhiAxis}}});
registry.add("h3CentdeltapTRndmConePhi_rhovsphi", "centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho}; #Delta#varphi_{jet}", {HistType::kTH3F, {{100, 0.0, 100.0}, {400, -200.0, 200.0}, {dPhiAxis}}});
}
// Normalize MC spectra with hUsedEventWeights bin 1 content, not GetEntries().
if (doprocessMCParticleLevel) {
LOGF(info, "JetSpectraEseTask::init() - MC Particle level");
registry.add("mcp/hEventCounter", "event status;event status;entries", {HistType::kTH1F, {{10, 0.0, 10.0}}});
registry.add("mcp/hUsedEventWeights", "Used MC collisions;Used events;Sum of event weights", {HistType::kTH1D, {{1, 0.0, 1.0}}});
registry.get<TH1>(HIST("mcp/hUsedEventWeights"))->Sumw2();
registry.get<TH1>(HIST("mcp/hUsedEventWeights"))->GetXaxis()->SetBinLabel(1, "Sum of weights of used events");
registry.add("mcp/hCentralitySel", ";centr;entries", {HistType::kTH1F, {{centAxis}}});
// registry.add("mcp/hJetSparse", ";Centrality;#it{p}_{T,jet part}; #eta; #phi", {HistType::kTHnSparseF, {{centAxis}, {jetPtAxis}, {etaAxis}, {phiAxis}}});
registry.add("mcp/hJetSparse", ";Centrality;#it{p}_{T,jet part}; #eta; #phi", {HistType::kTHnSparseF, {{centAxis}, {jetPtAxis}, {etaAxis}, {phiAxis}}});
Expand All @@ -560,6 +565,9 @@ struct JetSpectraEseTask {
if (doprocessMCDetectorLevel) {
LOGF(info, "JetSpectraEseTask::init() - MC Detector level");
registry.add("mcd/hEventCounter", "event status;event status;entries", {HistType::kTH1F, {{10, 0.0, 10.0}}});
registry.add("mcd/hUsedEventWeights", "Used reconstructed collisions;Used events;Sum of event weights", {HistType::kTH1D, {{1, 0.0, 1.0}}});
registry.get<TH1>(HIST("mcd/hUsedEventWeights"))->Sumw2();
registry.get<TH1>(HIST("mcd/hUsedEventWeights"))->GetXaxis()->SetBinLabel(1, "Sum of weights of used events");
registry.add("mcd/hCentralitySel", ";centr;entries", {HistType::kTH1F, {{centAxis}}});
// registry.add("mcd/hJetPt", "particle level jet pT;#it{p}_{T,jet part} (GeV/#it{c});entries", {HistType::kTH1F, {{jetPtAxis}}});
// registry.add("mcd/hJetSparse", ";Centrality;#it{p}_{T,jet det}; #eta; #phi", {HistType::kTHnSparseF, {{centAxis}, {jetPtAxis}, {etaAxis}, {phiAxis}}});
Expand All @@ -582,6 +590,9 @@ struct JetSpectraEseTask {

registry.add("mcm/hMCEventCounter", "event status;event status;entries", {HistType::kTH1F, {{10, 0.0, 10.0}}});
registry.add("mcm/hMCDMatchedEventCounter", "event status;event status;entries", {HistType::kTH1F, {{10, 0.0, 10.0}}});
registry.add("mcm/hUsedEventWeights", "Used reconstructed collisions with exactly one reconstructed collision per MC collision;Used events;Sum of event weights", {HistType::kTH1D, {{1, 0.0, 1.0}}});
registry.get<TH1>(HIST("mcm/hUsedEventWeights"))->Sumw2();
registry.get<TH1>(HIST("mcm/hUsedEventWeights"))->GetXaxis()->SetBinLabel(1, "Sum of weights of used events");
registry.add("mcm/hCentralityAnalyzed", ";Centrality;entries", {HistType::kTH1F, {{centAxis}}});

registry.add("mcm/hMatchedJetsPtDelta", "#it{p}_{T,jet part}; det - part", {HistType::kTH2F, {{jetPtAxis}, {100, -20., 20.0}}});
Expand Down Expand Up @@ -1184,14 +1195,18 @@ struct JetSpectraEseTask {
registry.fill(HIST("mcp/hCentralitySel"), centrality);
const float eventWeight = cfgUseMCEventWeights ? mcCollision.weight() : 1.0f;
const float pTHat = getPTHat(eventWeight);
if (!(pTHat < pTHatAbsoluteMin)) {
registry.fill(HIST("mcp/hUsedEventWeights"), 0.5, eventWeight);
}
jetLoopMCP<aod::JetParticles>(jets, centrality, mcCollision.rho(), eventWeight, pTHat);
}
PROCESS_SWITCH(JetSpectraEseTask, processMCParticleLevel, "jets on particle level MC", false);

void processMCDetectorLevel(soa::Filtered<soa::Join<aod::JetCollisionsMCD, aod::BkgChargedRhos>>::iterator const& collision,
ChargedMCDJets const& mcdjets,
aod::JetTracks const&,
aod::JetParticles const&)
aod::JetParticles const&,
aod::JetMcCollisions const&)
{
float counter{0.5f};
registry.fill(HIST("mcd/hEventCounter"), counter++);
Expand All @@ -1212,6 +1227,17 @@ struct JetSpectraEseTask {
auto centrality = cfgisPbPb ? getCentrality(collision) : -1;

registry.fill(HIST("mcd/hCentralitySel"), centrality);
float eventWeight = 1.0f;
if (cfgUseMCEventWeights) {
eventWeight = collision.has_mcCollision() ? collision.mcCollision_as<aod::JetMcCollisions>().weight() : -1.0f;
}
if (eventWeight <= 0.0f) {
return;
}
const float pTHat = getPTHat(eventWeight);
if (!(pTHat < pTHatAbsoluteMin)) {
registry.fill(HIST("mcd/hUsedEventWeights"), 0.5, eventWeight);
}
jetLoopMCD<aod::JetTracks>(mcdjets, centrality, collision.rho());
}
PROCESS_SWITCH(JetSpectraEseTask, processMCDetectorLevel, "jets on detector level", false);
Expand Down Expand Up @@ -1264,6 +1290,9 @@ struct JetSpectraEseTask {
registry.fill(HIST("mcm/hCentralityAnalyzed"), centrality);
float eventWeight = cfgUseMCEventWeights ? mcCol.weight() : 1.0;
float pTHat = getPTHat(eventWeight);
if (!(pTHat < pTHatAbsoluteMin)) {
registry.fill(HIST("mcm/hUsedEventWeights"), 0.5, eventWeight);
}
matchedJetLoop<JetMCPTable, aod::JetTracks>(mcdjets.sliceBy(mcdjetsPerJCollision, collision.globalIndex()), centrality, collision.rho(), mcCol.rho(), eventWeight, pTHat);

registry.fill(HIST("mcm/hMCDMatchedEventCounter"), secCount++);
Expand Down Expand Up @@ -1914,6 +1943,9 @@ struct JetSpectraEseTask {
weight = jet.eventWeight();
}
const float pTHat = getPTHat(weight);
if (pTHat < pTHatAbsoluteMin) {
continue;
}
if (jet.pt() > pTHatMaxMCD * pTHat) {
continue;
}
Expand All @@ -1924,6 +1956,9 @@ struct JetSpectraEseTask {
template <typename JTracks, typename Jets>
void jetLoopMCP(const Jets& jets, const float& centrality, const float& rho, float weight = 1.0f, float pTHat = 999.0f)
{
if (pTHat < pTHatAbsoluteMin) {
return;
}
bool mcLevelIsParticleLevel = true;
for (const auto& jet : jets) {
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackCuts.trackEtaMin, trackCuts.trackEtaMax)) {
Expand All @@ -1946,6 +1981,9 @@ struct JetSpectraEseTask {
template <typename MCPTab, typename JTracks, typename Jets>
void matchedJetLoop(const Jets& jets, const float& centrality, const float& rho, const float& rho2, float weight = 1.0, float pTHat = 999.0)
{
if (pTHat < pTHatAbsoluteMin) {
return;
}
for (const auto& jet : jets) {
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackCuts.trackEtaMin, trackCuts.trackEtaMax)) {
continue;
Expand Down
Loading