From 24d6b5b9461da3d9c1132e588c5bb1b223452c56 Mon Sep 17 00:00:00 2001 From: Joonsuk Date: Sun, 24 Aug 2025 02:16:26 +0900 Subject: [PATCH 1/2] PWGJE: QA process for track resolution added --- PWGJE/Tasks/jetBackgroundAnalysis.cxx | 30 ++++++++++++++++++++++++++- 1 file changed, 29 insertions(+), 1 deletion(-) diff --git a/PWGJE/Tasks/jetBackgroundAnalysis.cxx b/PWGJE/Tasks/jetBackgroundAnalysis.cxx index 8e93fc305d0..26097abcbde 100644 --- a/PWGJE/Tasks/jetBackgroundAnalysis.cxx +++ b/PWGJE/Tasks/jetBackgroundAnalysis.cxx @@ -82,6 +82,7 @@ struct JetBackgroundAnalysisTask { // histogram definitions if (doprocessRho) { + registry.add("h2_leadingjet_pt_rho", "; #it{p}_{T, leading jet} (GeV/#it{c}); #it{rho} (GeV/area);", {HistType::kTH2F, {{200, 0., 200.0}, {400, 0.0, 400.0}}}); registry.add("h2_centrality_ntracks", "; centrality; N_{tracks};", {HistType::kTH2F, {{1100, 0., 110.0}, {10000, 0.0, 10000.0}}}); registry.add("h2_ntracks_rho", "; N_{tracks}; #it{rho} (GeV/area);", {HistType::kTH2F, {{10000, 0.0, 10000.0}, {400, 0.0, 400.0}}}); registry.add("h2_ntracks_rhom", "; N_{tracks}; #it{rho}_{m} (GeV/area);", {HistType::kTH2F, {{10000, 0.0, 10000.0}, {100, 0.0, 100.0}}}); @@ -90,9 +91,14 @@ struct JetBackgroundAnalysisTask { } if (doprocessBkgFluctuationsData || doprocessBkgFluctuationsMCD) { + registry.add("h2_leadingjet_pt_rhorandomcone", "; #it{p}_{T, leading jet} (GeV/#it{c}); #it{rho} (GeV/area);", {HistType::kTH2F, {{200, 0., 200.0}, {400, 0.0, 400.0}}}); registry.add("h2_centrality_rhorandomcone", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, bkgFluctuationsAxis}}); + registry.add("h2_leadingjet_pt_rhorandomconerandomtrackdirection", "; #it{p}_{T, leading jet} (GeV/#it{c}); #it{rho} (GeV/area);", {HistType::kTH2F, {{200, 0., 200.0}, {400, 0.0, 400.0}}}); registry.add("h2_centrality_rhorandomconerandomtrackdirection", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, bkgFluctuationsAxis}}); + registry.add("h2_leadingjet_pt_rhorandomconewithoutleadingjet", "; #it{p}_{T, leading jet} (GeV/#it{c}); #it{rho} (GeV/area);", {HistType::kTH2F, {{200, 0., 200.0}, {400, 0.0, 400.0}}}); registry.add("h2_centrality_rhorandomconewithoutleadingjet", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, bkgFluctuationsAxis}}); + registry.add("h2_leadingjet_pt_rhorandomconerandomtrackdirectionwithoutoneleadingjets", "; #it{p}_{T, leading jet} (GeV/#it{c}); #it{rho} (GeV/area);", {HistType::kTH2F, {{200, 0., 200.0}, {400, 0.0, 400.0}}}); + registry.add("h2_leadingjet_pt_rhorandomconerandomtrackdirectionwithouttwoleadingjets", "; #it{p}_{T, leading jet} (GeV/#it{c}); #it{rho} (GeV/area);", {HistType::kTH2F, {{200, 0., 200.0}, {400, 0.0, 400.0}}}); registry.add("h2_centrality_rhorandomconerandomtrackdirectionwithoutoneleadingjets", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, bkgFluctuationsAxis}}); registry.add("h2_centrality_rhorandomconerandomtrackdirectionwithouttwoleadingjets", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, bkgFluctuationsAxis}}); } @@ -128,6 +134,13 @@ struct JetBackgroundAnalysisTask { } } } + + // Check if jets exist before accessing them + std::cout << "bkg fluctuations jets size: " << jets.size() << std::endl; + if (jets.size() > 0) { + std::cout << "leading jet pt random cone: " << jets.iteratorAt(0).pt() << std::endl; + registry.fill(HIST("h2_leadingjet_pt_rhorandomcone"), jets.iteratorAt(0).pt(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho()); + } registry.fill(HIST("h2_centrality_rhorandomcone"), collision.centFT0M(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho()); // randomised eta,phi for tracks, to assess part of fluctuations coming from statistically independently emitted particles @@ -141,6 +154,9 @@ struct JetBackgroundAnalysisTask { } } } + if (jets.size() > 0) { + registry.fill(HIST("h2_leadingjet_pt_rhorandomconerandomtrackdirection"), jets.iteratorAt(0).pt(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho()); + } registry.fill(HIST("h2_centrality_rhorandomconerandomtrackdirection"), collision.centFT0M(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho()); // removing the leading jet from the random cone @@ -169,6 +185,9 @@ struct JetBackgroundAnalysisTask { } } } + if (jets.size() > 0) { + registry.fill(HIST("h2_leadingjet_pt_rhorandomconewithoutleadingjet"), jets.iteratorAt(0).pt(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho()); + } registry.fill(HIST("h2_centrality_rhorandomconewithoutleadingjet"), collision.centFT0M(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho()); // randomised eta,phi for tracks, to assess part of fluctuations coming from statistically independently emitted particles, removing tracks from 2 leading jets @@ -190,11 +209,15 @@ struct JetBackgroundAnalysisTask { } } } + if (jets.size() > 0) { + registry.fill(HIST("h2_leadingjet_pt_rhorandomconerandomtrackdirectionwithoutoneleadingjets"), jets.iteratorAt(0).pt(), randomConePtWithoutOneLeadJet - M_PI * randomConeR * randomConeR * collision.rho()); + registry.fill(HIST("h2_leadingjet_pt_rhorandomconerandomtrackdirectionwithouttwoleadingjets"), jets.iteratorAt(0).pt(), randomConePtWithoutTwoLeadJet - M_PI * randomConeR * randomConeR * collision.rho()); + } registry.fill(HIST("h2_centrality_rhorandomconerandomtrackdirectionwithoutoneleadingjets"), collision.centFT0M(), randomConePtWithoutOneLeadJet - M_PI * randomConeR * randomConeR * collision.rho()); registry.fill(HIST("h2_centrality_rhorandomconerandomtrackdirectionwithouttwoleadingjets"), collision.centFT0M(), randomConePtWithoutTwoLeadJet - M_PI * randomConeR * randomConeR * collision.rho()); } - void processRho(soa::Filtered>::iterator const& collision, soa::Filtered const& tracks) + void processRho(soa::Filtered>::iterator const& collision, soa::Filtered const& tracks, aod::ChargedJets const& jets) { if (!jetderiveddatautilities::selectCollision(collision, eventSelectionBits, skipMBGapEvents)) { return; @@ -208,6 +231,11 @@ struct JetBackgroundAnalysisTask { nTracks++; } } + std::cout << "process rho jets size: " << jets.size() << std::endl; + if (jets.size() > 0) { + std::cout << "leading jet pt rho: " << jets.iteratorAt(0).pt() << std::endl; + registry.fill(HIST("h2_leadingjet_pt_rho"), jets.iteratorAt(0).pt(), collision.rho()); + } registry.fill(HIST("h2_centrality_ntracks"), collision.centFT0M(), nTracks); registry.fill(HIST("h2_ntracks_rho"), nTracks, collision.rho()); registry.fill(HIST("h2_ntracks_rhom"), nTracks, collision.rhoM()); From 28336607ad05970a446827780ae8d75bd2097363 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Sat, 23 Aug 2025 17:19:08 +0000 Subject: [PATCH 2/2] Please consider the following formatting changes --- PWGJE/Tasks/jetBackgroundAnalysis.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGJE/Tasks/jetBackgroundAnalysis.cxx b/PWGJE/Tasks/jetBackgroundAnalysis.cxx index 26097abcbde..ceb73649361 100644 --- a/PWGJE/Tasks/jetBackgroundAnalysis.cxx +++ b/PWGJE/Tasks/jetBackgroundAnalysis.cxx @@ -134,7 +134,7 @@ struct JetBackgroundAnalysisTask { } } } - + // Check if jets exist before accessing them std::cout << "bkg fluctuations jets size: " << jets.size() << std::endl; if (jets.size() > 0) {