Skip to content

[PWGLF] Add jet pt histograms #18284

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
alcaliva wants to merge 2 commits into
base: master
Choose a base branch
from
Open

[PWGLF] Add jet pt histograms #18284

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
6e7c820
[PWGLF] Add jet pt histograms
alcaliva Oct 14, 2026
e9f9e54
restarting the CIs
alcaliva Oct 15, 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
23 changes: 22 additions & 1 deletion PWGLF/Tasks/Nuspex/antinucleiInJets.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -282,7 +282,7 @@ struct AntinucleiInJets {
getReweightingHistograms(ccdb, TString(pathToFile));
}

// Binning
// pt binning
double min = 0.0;
double max = 6.0;
int nbins = 120;
Expand Down Expand Up @@ -326,6 +326,9 @@ struct AntinucleiInJets {
registryData.add("jetEffectiveAreaOverPiR2", "jet effective area / piR^2", HistType::kTH1F, {{2000, 0, 2, "Area/#piR^{2}"}});
registryData.add("jetArea", "jetArea", HistType::kTH1F, {{5000, 0, 5, "area"}});

// Jet pt distribution in data
registryData.add("jetPtData", "jet pt distribution in data", HistType::kTH1F, {{2000, 0, 200, "#it{p}^{jet}_{T} (GeV/#it{c})"}});

// angle between track and jet axis
registryData.add("theta_track_jet", "theta_track_jet", HistType::kTH2F, {{100, 0, 100, "#it{p}^{jet}_{T} (GeV/#it{c})"}, {400, 0, 20.0, "#theta_{track-jet} (deg)"}});

Expand Down Expand Up @@ -379,6 +382,9 @@ struct AntinucleiInJets {
registryMC.add("genEvents", "number of generated events in mc", HistType::kTH1F, {{10, 0, 10, "counter"}});
registryMC.add("genJets", "number of generated jets", HistType::kTH1F, {{10, 0, 10, "counter"}});

// Jet pt distribution in gen MC
registryMC.add("jetPtMCgenerated", "jet pt distribution in gen MC", HistType::kTH1F, {{2000, 0, 200, "#it{p}^{jet}_{T} (GeV/#it{c})"}});

// Size of particle array
registryMC.add("sizeParticleArray", "number of particles", HistType::kTH1F, {{1000, 0, 1000, "#it{N}_{part}"}});

Expand Down Expand Up @@ -415,6 +421,9 @@ struct AntinucleiInJets {
registryMC.add("recEvents", "number of reconstructed events in mc", HistType::kTH1F, {{20, 0, 20, "counter"}});
registryMC.add("recJets", "number of reconstructed jets", HistType::kTH1F, {{10, 0, 10, "counter"}});

// Jet pt distribution in rec MC
registryMC.add("jetPtMCreconstructed", "jet pt distribution in rec MC", HistType::kTH1F, {{2000, 0, 200, "#it{p}^{jet}_{T} (GeV/#it{c})"}});

// Configuration
registryMC.add("settingMC", "settingMC", HistType::kTH2F, {{100, 0.0, 50.0, "min #it{p}^{jet}_{T} (GeV/#it{c})"}, {20, 0.0, 1.0, "#it{R}_{jet}"}});

Expand Down Expand Up @@ -1303,6 +1312,10 @@ struct AntinucleiInJets {
// Jet pt must be larger than threshold
auto jetForSub = jet;
fastjet::PseudoJet jetMinusBkg = backgroundSub.doRhoAreaSub(jetForSub, rhoPerp, rhoMPerp);

// Fill jet pt distribution
registryData.fill(HIST("jetPtData"), jetMinusBkg.pt());

if (isppRefAnalysis && (jet.pt() < minJetPt || jet.pt() > maxJetPt))
continue;
if (!isppRefAnalysis && (jetMinusBkg.pt() < minJetPt || jetMinusBkg.pt() > maxJetPt))
Expand Down Expand Up @@ -2141,6 +2154,10 @@ struct AntinucleiInJets {
// Jet pt must be larger than threshold
auto jetForSub = jet;
fastjet::PseudoJet jetMinusBkg = backgroundSub.doRhoAreaSub(jetForSub, rhoPerp, rhoMPerp);

// Fill jet pt distribution
registryMC.fill(HIST("jetPtMCgenerated"), jetMinusBkg.pt());

if (isppRefAnalysis && (jet.pt() < minJetPt || jet.pt() > maxJetPt))
continue;
if (!isppRefAnalysis && (jetMinusBkg.pt() < minJetPt || jetMinusBkg.pt() > maxJetPt))
Expand Down Expand Up @@ -2440,6 +2457,10 @@ struct AntinucleiInJets {
// Jet pt must be larger than threshold
auto jetForSub = jet;
fastjet::PseudoJet jetMinusBkg = backgroundSub.doRhoAreaSub(jetForSub, rhoPerp, rhoMPerp);

// Fill jet pt distribution
registryMC.fill(HIST("jetPtMCreconstructed"), jetMinusBkg.pt());

if (isppRefAnalysis && (jet.pt() < minJetPt || jet.pt() > maxJetPt))
continue;
if (!isppRefAnalysis && (jetMinusBkg.pt() < minJetPt || jetMinusBkg.pt() > maxJetPt))
Expand Down
Loading