Skip to content

[PWGJE] add occupancy cut and PbPb checks #15564

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

Merged
nzardosh merged 1 commit into from
Mar 30, 2026
+19 −2
Merged

[PWGJE] add occupancy cut and PbPb checks #15564

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
21 changes: 19 additions & 2 deletions PWGJE/Tasks/jetHadronRecoil.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -70,6 +70,7 @@ struct JetHadronRecoil {
Configurable<float> centralityMin{"centralityMin", -999.0, "minimum centrality"};
Configurable<float> centralityMax{"centralityMax", 999.0, "maximum centrality"};
Configurable<float> vertexZCut{"vertexZCut", 10.0f, "Accepted z-vertex range"};
Configurable<int> trackOccupancyInTimeRangeMax{"trackOccupancyInTimeRangeMax", 999999, "maximum track occupancy in time range"};
Configurable<float> ptTTrefMin{"ptTTrefMin", 5, "reference minimum trigger track pt"};
Configurable<float> ptTTrefMax{"ptTTrefMax", 7, "reference maximum trigger track pt"};
Configurable<float> ptTTsigMin{"ptTTsigMin", 20, "signal minimum trigger track pt"};
Expand All @@ -93,7 +94,7 @@ struct JetHadronRecoil {
Filter jetCuts = aod::jet::r == nround(jetR.node() * 100.0f);
Filter trackCuts = (aod::jtrack::pt >= trackPtMin && aod::jtrack::pt < trackPtMax && aod::jtrack::eta > trackEtaMin && aod::jtrack::eta < trackEtaMax);
Filter particleCuts = (aod::jmcparticle::pt >= trackPtMin && aod::jmcparticle::pt < trackPtMax && aod::jmcparticle::eta > trackEtaMin && aod::jmcparticle::eta < trackEtaMax);
Filter eventCuts = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centFT0M >= centralityMin && aod::jcollision::centFT0M < centralityMax);
Filter eventCuts = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centFT0M >= centralityMin && aod::jcollision::centFT0M < centralityMax && aod::jcollision::trackOccupancyInTimeRange < trackOccupancyInTimeRangeMax);

std::vector<double> ptBinningPart = {0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0,
15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0, 55.0, 60.0,
Expand Down Expand Up @@ -147,8 +148,11 @@ struct JetHadronRecoil {
jetReclusterer.recombScheme = fastjet::WTA_pt_scheme;

registry.add("hZvtxSelected", "Z vertex position;Z_{vtx};entries", {HistType::kTH1F, {{80, -20, 20}}}, doSumw);
registry.add("hCentOccupancy", "centrality vs track occupancy;centrality (%);track occupancy", {HistType::kTH2F, {{100, 0, 100}, {200, 0, 10000}}}, doSumw);
registry.add("hCentRho", "centrality vs rho;centrality (%);#rho", {HistType::kTH2F, {{100, 0, 100}, {125, 0, 250}}}, doSumw);

if (doprocessData || doprocessDataWithRhoSubtraction || doprocessMCD || doprocessMCDWithRhoSubtraction || doprocessMCDWeighted || doprocessMCDWeightedWithRhoSubtraction || doprocessMCP || doprocessMCPWeighted || doprocessMCPOnTheFly || doprocessMCPOnTheFlyWeighted) {
registry.add("hNtracksRho", "track multiplicity vs #rho;N_{tracks};#rho", {HistType::kTH2F, {{125, 0, 5000}, {125, 0, 250}}}, doSumw);
registry.add("hNtrig", "number of triggers;trigger type;entries", {HistType::kTH1F, {{2, 0, 2}}}, doSumw);
registry.add("hSignalTriggersPtHard", "Signal triggers vs PtHard", {HistType::kTH1F, {pThatAxis}}, doSumw);
registry.add("hReferenceTriggersPtHard", "Reference triggers vs PtHard", {HistType::kTH1F, {pThatAxis}}, doSumw);
Expand All @@ -168,6 +172,7 @@ struct JetHadronRecoil {
registry.add("hJetEta", "jet #eta;#eta_{jet};entries", {HistType::kTH1F, {{100, -1.0, 1.0}}}, doSumw);
registry.add("hJetPhi", "jet #phi;#phi_{jet};entries", {HistType::kTH1F, {{100, 0.0, o2::constants::math::TwoPI}}}, doSumw);
registry.add("hJet3D", "3D jet distribution;p_{T};#eta;#phi", {HistType::kTH3F, {{500, -100, 400}, {100, -1.0, 1.0}, {100, 0.0, o2::constants::math::TwoPI}}}, doSumw);
registry.add("hJetPtArea", "jet p_{T}vs area;jet p_{T};A_{jet}", {HistType::kTH2F, {{500, -100, 400}, {100, 0, 2}}}, doSumw);
}

if (doprocessData || doprocessDataWithRhoSubtraction || doprocessMCD || doprocessMCDWithRhoSubtraction || doprocessMCDWeighted || doprocessMCDWeightedWithRhoSubtraction) {
Expand Down Expand Up @@ -253,6 +258,7 @@ struct JetHadronRecoil {
double phiTT = 0;
double ptTT = 0;
int nTT = 0;
int nTracks = 0;
double leadingPT = 0;
double leadingTrackPt = 0;
double leadingJetPt = 0;
Expand Down Expand Up @@ -295,7 +301,9 @@ struct JetHadronRecoil {
registry.fill(HIST("hPhiTrack"), track.phi(), weight);
registry.fill(HIST("hTrack3D"), track.pt(), track.eta(), track.phi(), weight);
registry.fill(HIST("hPtTrackPtHard"), track.pt() / pTHat, track.pt(), weight);
nTracks++;
}
registry.fill(HIST("hNtracksRho"), nTracks, rho);
if (nTT > 0) {
int trigNumber = rand->Integer(nTT);
phiTT = phiTTAr[trigNumber];
Expand Down Expand Up @@ -340,6 +348,7 @@ struct JetHadronRecoil {
registry.fill(HIST("hJetEta"), jet.eta(), weight);
registry.fill(HIST("hJetPhi"), jet.phi(), weight);
registry.fill(HIST("hJet3D"), jet.pt() - (rho * jet.area()), jet.eta(), jet.phi(), weight);
registry.fill(HIST("hJetPtArea"), jet.pt(), jet.area(), weight);

if (nTT > 0) {
float dphi = RecoDecay::constrainAngle(jet.phi() - phiTT);
Expand Down Expand Up @@ -388,6 +397,7 @@ struct JetHadronRecoil {
double phiTT = 0;
double ptTT = 0;
int nTT = 0;
int nTracks = 0;
double leadingPT = 0;
double leadingTrackPt = 0;
double leadingJetPt = 0;
Expand Down Expand Up @@ -430,7 +440,9 @@ struct JetHadronRecoil {
registry.fill(HIST("hPhiTrack"), track.phi(), weight);
registry.fill(HIST("hTrack3D"), track.pt(), track.eta(), track.phi(), weight);
registry.fill(HIST("hPtTrackPtHard"), track.pt() / pTHat, track.pt(), weight);
nTracks++;
}
registry.fill(HIST("hNtracksRho"), nTracks, rho);
if (nTT > 0) {
int trigNumber = rand->Integer(nTT);
phiTT = phiTTAr[trigNumber];
Expand Down Expand Up @@ -641,7 +653,6 @@ struct JetHadronRecoil {
const auto jetsTagColl = jetsTag.sliceBy(perColJetMC, collisionID);

for (const auto& jetTag : jetsTagColl) {

if (jetTag.pt() > pTHatMaxMCP * pTHat) {
if (outlierRejectEvent) {
return;
Expand Down Expand Up @@ -859,6 +870,8 @@ struct JetHadronRecoil {
return;
}
registry.fill(HIST("hZvtxSelected"), collision.posZ());
registry.fill(HIST("hCentOccupancy"), collision.centFT0M(), collision.trackOccupancyInTimeRange());
registry.fill(HIST("hCentRho"), collision.centFT0M(), collision.rho());
fillHistograms(jets, tracks, 1.0, collision.rho());
}
PROCESS_SWITCH(JetHadronRecoil, processDataWithRhoSubtraction, "process data with rho subtraction", false);
Expand Down Expand Up @@ -903,6 +916,8 @@ struct JetHadronRecoil {
return;
}
registry.fill(HIST("hZvtxSelected"), collision.posZ());
registry.fill(HIST("hCentOccupancy"), collision.centFT0M(), collision.trackOccupancyInTimeRange());
registry.fill(HIST("hCentRho"), collision.centFT0M(), collision.rho());
fillHistogramsMCD(jets, tracks, 1.0, collision.rho(), collision.mcCollision().ptHard());
}
PROCESS_SWITCH(JetHadronRecoil, processMCDWithRhoSubtraction, "process MC detector level with rho subtraction", false);
Expand Down Expand Up @@ -947,6 +962,8 @@ struct JetHadronRecoil {
return;
}
registry.fill(HIST("hZvtxSelected"), collision.posZ(), collision.mcCollision().weight());
registry.fill(HIST("hCentOccupancy"), collision.centFT0M(), collision.trackOccupancyInTimeRange(), collision.mcCollision().weight());
registry.fill(HIST("hCentRho"), collision.centFT0M(), collision.rho());
fillHistogramsMCD(jets, tracks, collision.mcCollision().weight(), collision.rho(), collision.mcCollision().ptHard());
}
PROCESS_SWITCH(JetHadronRecoil, processMCDWeightedWithRhoSubtraction, "process MC detector level with event weights and rho subtraction", false);
Expand Down
Loading