From 6197da4c7f9b41ebe50148caee685cea0e671f5c Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 28 Mar 2024 16:09:53 +0000 Subject: [PATCH] Please consider the following formatting changes --- PWGJE/Tasks/jetChCorr.cxx | 328 +++++++++++++++++++------------------- 1 file changed, 162 insertions(+), 166 deletions(-) diff --git a/PWGJE/Tasks/jetChCorr.cxx b/PWGJE/Tasks/jetChCorr.cxx index 53064dc16ad..a652c2f489b 100644 --- a/PWGJE/Tasks/jetChCorr.cxx +++ b/PWGJE/Tasks/jetChCorr.cxx @@ -54,11 +54,11 @@ int trackSelection = -1; int trackL = -1; int n_trackL = -1; -TVector3 v1, v2, vR; // leading and next to leading vector -int ch_l =0; // charge leading -int ch_nl =0; // charge next to leading +TVector3 v1, v2, vR; // leading and next to leading vector +int ch_l = 0; // charge leading +int ch_nl = 0; // charge next to leading -int ch_mult = 0; // charge next to leading +int ch_mult = 0; // charge next to leading #include "Framework/runDataProcessing.h" @@ -92,7 +92,6 @@ struct JetChCorr { registry.add("h2_jet_pt_jet_rg_eventwiseconstituentsubtracted", ";#it{p}_{T,jet} (GeV/#it{c});#it{R}_{g}", {HistType::kTH2F, {{200, 0., 200.}, {22, 0.0, 1.1}}}); registry.add("h2_jet_pt_jet_nsd_eventwiseconstituentsubtracted", ";#it{p}_{T,jet} (GeV/#it{c});#it{n}_{SD}", {HistType::kTH2F, {{200, 0., 200.}, {15, -0.5, 14.5}}}); - //// For L and NL particles //// For L and NL particles @@ -101,7 +100,7 @@ struct JetChCorr { registry.add("h_ch_s_kt", ";#it{k}_{T,pair} (GeV/#it{c});N_{ch} (hh)", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("h_ch_d_kt", ";#it{k}_{T,pair} (GeV/#it{c});N_{ch} (h#bar{h})", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("h_ch_s_z", ";#it{z}_{pair} ;N_{ch} (hh)", {HistType::kTH1F, {{20, 0., .5}}}); - registry.add("h_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}}", {HistType::kTH1F, {{20, 0.,.5}}}); + registry.add("h_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}}", {HistType::kTH1F, {{20, 0., .5}}}); registry.add("h_ch_s_ft", ";#it{fm}_{time} (GeV/#it{c});N_{ch} (hh)", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("h_ch_d_ft", ";#it{fm}_{time} (GeV/#it{c});N_{ch} (h#bar{h})", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("h_ch_s_mult", ";#it{N}_{ch,jet} (GeV/#it{c});N_{ch} (hh)", {HistType::kTH1F, {{15, 0., 30.}}}); @@ -112,7 +111,7 @@ struct JetChCorr { registry.add("hr1_ch_s_kt", ";#it{k}_{T,pair} (GeV/#it{c}) {nsd1};N_{ch} (hh)", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("hr1_ch_d_kt", ";#it{k}_{T,pair} (GeV/#it{c}) {nsd1};N_{ch} (h#bar{h})", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("hr1_ch_s_z", ";#it{z}_{pair} ;N_{ch} (hh) {nsd1}", {HistType::kTH1F, {{20, 0., .5}}}); - registry.add("hr1_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}} {nsd1}", {HistType::kTH1F, {{20, 0.,.5}}}); + registry.add("hr1_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}} {nsd1}", {HistType::kTH1F, {{20, 0., .5}}}); registry.add("hr1_ch_s_ft", ";#it{fm}_{time} (GeV/#it{c}) {nsd1};N_{ch} (hh)", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("hr1_ch_d_ft", ";#it{fm}_{time} (GeV/#it{c}) {nsd1};N_{ch} (h#bar{h})", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("hr1_ch_s_mult", ";#it{N}_{ch,jet} (GeV/#it{c}) {nsd1};N_{ch} (hh)", {HistType::kTH1F, {{15, 0., 30.}}}); @@ -129,7 +128,7 @@ struct JetChCorr { registry.add("hr2_ch_s_kt", ";#it{k}_{T,pair} (GeV/#it{c}) {nsd2};N_{ch} (hh)", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("hr2_ch_d_kt", ";#it{k}_{T,pair} (GeV/#it{c}) {nsd2};N_{ch} (h#bar{h})", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("hr2_ch_s_z", ";#it{z}_{pair} ;N_{ch} (hh) {nsd2}", {HistType::kTH1F, {{20, 0., .5}}}); - registry.add("hr2_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}} {nsd2}", {HistType::kTH1F, {{20, 0.,.5}}}); + registry.add("hr2_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}} {nsd2}", {HistType::kTH1F, {{20, 0., .5}}}); registry.add("hr2_ch_s_ft", ";#it{fm}_{time} (GeV/#it{c}) {nsd2};N_{ch} (hh)", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("hr2_ch_d_ft", ";#it{fm}_{time} (GeV/#it{c}) {nsd2};N_{ch} (h#bar{h})", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("hr2_ch_s_mult", ";#it{N}_{ch,jet} (GeV/#it{c}) {nsd2};N_{ch} (hh)", {HistType::kTH1F, {{15, 0., 30.}}}); @@ -146,7 +145,7 @@ struct JetChCorr { registry.add("hr3_ch_s_kt", ";#it{k}_{T,pair} (GeV/#it{c}) {nsd3};N_{ch} (hh)", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("hr3_ch_d_kt", ";#it{k}_{T,pair} (GeV/#it{c}) {nsd3};N_{ch} (h#bar{h})", {HistType::kTH1F, {{50, 0., 5.}}}); registry.add("hr3_ch_s_z", ";#it{z}_{pair} ;N_{ch} (hh) {nsd3}", {HistType::kTH1F, {{20, 0., .5}}}); - registry.add("hr3_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}} {nsd3}", {HistType::kTH1F, {{20, 0.,.5}}}); + registry.add("hr3_ch_d_z", ";#it{z}_{pair} ;ch_{h#bar{h}} {nsd3}", {HistType::kTH1F, {{20, 0., .5}}}); registry.add("hr3_ch_s_ft", ";#it{fm}_{time} (GeV/#it{c}) {nsd3};N_{ch} (hh)", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("hr3_ch_d_ft", ";#it{fm}_{time} (GeV/#it{c}) {nsd3};N_{ch} (h#bar{h})", {HistType::kTH1F, {{20, -1., 4.}}}); registry.add("hr3_ch_s_mult", ";#it{N}_{ch,jet} (GeV/#it{c}) {nsd3};N_{ch} (hh)", {HistType::kTH1F, {{15, 0., 30.}}}); @@ -162,7 +161,7 @@ struct JetChCorr { jetReclusterer.algorithm = fastjet::JetAlgorithm::cambridge_algorithm; } - //template + // template template void jetReclustering(T const& jet, U& outputTable) { @@ -183,7 +182,6 @@ struct JetChCorr { float ptJet = daughterSubJet.pt(); - while (daughterSubJet.has_parents(parentSubJet1, parentSubJet2)) { if (parentSubJet1.perp() < parentSubJet2.perp()) { std::swap(parentSubJet1, parentSubJet2); @@ -191,22 +189,23 @@ struct JetChCorr { vector constituents1 = parentSubJet1.constituents(); vector constituents2 = parentSubJet2.constituents(); - //cout<<"in subJET1 ********************************************* "<= zCut * TMath::Power(theta / (jet.r() / 100.f), beta)) { - if(found1>0&&found2>0){ // found leading and next-to-leading in seperate prongs - if (!softDropped) { - zg = z; - rg = theta; - if constexpr (!isSubtracted && !isMCP) { - registry.fill(HIST("h2_jet_pt_jet_zg"), jet.pt(), zg); - registry.fill(HIST("h2_jet_pt_jet_rg"), jet.pt(), rg); - } - if constexpr (!isSubtracted && isMCP) { - registry.fill(HIST("h2_jet_pt_part_jet_zg_part"), jet.pt(), zg); - registry.fill(HIST("h2_jet_pt_part_jet_rg_part"), jet.pt(), rg); - } - if constexpr (isSubtracted && !isMCP) { - registry.fill(HIST("h2_jet_pt_jet_zg_eventwiseconstituentsubtracted"), jet.pt(), zg); - registry.fill(HIST("h2_jet_pt_jet_rg_eventwiseconstituentsubtracted"), jet.pt(), rg); - } - softDropped = true; - - - softDropped = true; - - v1.SetXYZ(parentSubJet1.px(), parentSubJet1.py(), parentSubJet1.pz()); - v2.SetXYZ(parentSubJet2.px(), parentSubJet2.py(), parentSubJet2.pz()); - - vR = v1 + v2; - float z = v2.Perp(vR.Orthogonal())/(v1.Perp(vR.Orthogonal())+v2.Perp(vR.Orthogonal())); - float fT = ((2.*z*(1-z)*vR.Mag())/v1.Perp2(vR))/6.; - float kt_p = v1.Perp(vR); - - //////// for first matched split - if(nsd==0){ - if(ch_l == ch_nl) { - registry.fill(HIST("hr1_ch_s_pt"), ptJet); - registry.fill(HIST("hr1_ch_s_kt"), kt_p); - registry.fill(HIST("hr1_ch_s_z"), z); - registry.fill(HIST("hr1_ch_s_ft"), log10(fT)); - registry.fill(HIST("hr1_ch_s_mult"), ch_mult); - registry.fill(HIST("hr1_ch_s_zg"), zg); - registry.fill(HIST("hr1_ch_s_rg"), rg); - registry.fill(HIST("hr1_ch_s_lu"), log(1/rg),log(kt_p)); - } - if(ch_l != ch_nl) { - registry.fill(HIST("hr1_ch_d_pt"), ptJet); - registry.fill(HIST("hr1_ch_d_kt"), kt_p); - registry.fill(HIST("hr1_ch_d_z"), z); - registry.fill(HIST("hr1_ch_d_ft"), log10(fT)); - registry.fill(HIST("hr1_ch_d_mult"), ch_mult); - registry.fill(HIST("hr1_ch_d_zg"), zg); - registry.fill(HIST("hr1_ch_d_rg"), rg); - registry.fill(HIST("hr1_ch_d_lu"), log(1/rg),log(kt_p)); - } - } - //////// for 2nd matched split - if(nsd==1){ - if(ch_l == ch_nl) { - registry.fill(HIST("hr2_ch_s_pt"), ptJet); - registry.fill(HIST("hr2_ch_s_kt"), kt_p); - registry.fill(HIST("hr2_ch_s_z"), z); - registry.fill(HIST("hr2_ch_s_ft"), log10(fT)); - registry.fill(HIST("hr2_ch_s_mult"), ch_mult); - registry.fill(HIST("hr2_ch_s_zg"), zg); - registry.fill(HIST("hr2_ch_s_rg"), rg); - registry.fill(HIST("hr2_ch_s_lu"), log(1/rg),log(kt_p)); - } - if(ch_l != ch_nl) { - registry.fill(HIST("hr2_ch_d_pt"), ptJet); - registry.fill(HIST("hr2_ch_d_kt"), kt_p); - registry.fill(HIST("hr2_ch_d_z"), z); - registry.fill(HIST("hr2_ch_d_ft"), log10(fT)); - registry.fill(HIST("hr2_ch_d_mult"), ch_mult); - registry.fill(HIST("hr2_ch_d_zg"), zg); - registry.fill(HIST("hr2_ch_d_rg"), rg); - registry.fill(HIST("hr2_ch_d_lu"), log(1/rg),log(kt_p)); - } - } - //////// for 3rd marched split - if(nsd>=2){ - if(ch_l == ch_nl) { - registry.fill(HIST("hr3_ch_s_pt"), ptJet); - registry.fill(HIST("hr3_ch_s_kt"), kt_p); - registry.fill(HIST("hr3_ch_s_z"), z); - registry.fill(HIST("hr3_ch_s_ft"), log10(fT)); - registry.fill(HIST("hr3_ch_s_mult"), ch_mult); - registry.fill(HIST("hr3_ch_s_zg"), zg); - registry.fill(HIST("hr3_ch_s_rg"), rg); - registry.fill(HIST("hr3_ch_s_lu"), log(1/rg),log(kt_p)); - } - if(ch_l != ch_nl) { - registry.fill(HIST("hr3_ch_d_pt"), ptJet); - registry.fill(HIST("hr3_ch_d_kt"), kt_p); - registry.fill(HIST("hr3_ch_d_z"), z); - registry.fill(HIST("hr3_ch_d_ft"), log10(fT)); - registry.fill(HIST("hr3_ch_d_mult"), ch_mult); - registry.fill(HIST("hr3_ch_d_zg"), zg); - registry.fill(HIST("hr3_ch_d_rg"), rg); - registry.fill(HIST("hr3_ch_d_lu"), log(1/rg),log(kt_p)); - } - } - } - } + if (found1 > 0 && found2 > 0) { // found leading and next-to-leading in seperate prongs + if (!softDropped) { + zg = z; + rg = theta; + if constexpr (!isSubtracted && !isMCP) { + registry.fill(HIST("h2_jet_pt_jet_zg"), jet.pt(), zg); + registry.fill(HIST("h2_jet_pt_jet_rg"), jet.pt(), rg); + } + if constexpr (!isSubtracted && isMCP) { + registry.fill(HIST("h2_jet_pt_part_jet_zg_part"), jet.pt(), zg); + registry.fill(HIST("h2_jet_pt_part_jet_rg_part"), jet.pt(), rg); + } + if constexpr (isSubtracted && !isMCP) { + registry.fill(HIST("h2_jet_pt_jet_zg_eventwiseconstituentsubtracted"), jet.pt(), zg); + registry.fill(HIST("h2_jet_pt_jet_rg_eventwiseconstituentsubtracted"), jet.pt(), rg); + } + softDropped = true; + + softDropped = true; + + v1.SetXYZ(parentSubJet1.px(), parentSubJet1.py(), parentSubJet1.pz()); + v2.SetXYZ(parentSubJet2.px(), parentSubJet2.py(), parentSubJet2.pz()); + + vR = v1 + v2; + float z = v2.Perp(vR.Orthogonal()) / (v1.Perp(vR.Orthogonal()) + v2.Perp(vR.Orthogonal())); + float fT = ((2. * z * (1 - z) * vR.Mag()) / v1.Perp2(vR)) / 6.; + float kt_p = v1.Perp(vR); + + //////// for first matched split + if (nsd == 0) { + if (ch_l == ch_nl) { + registry.fill(HIST("hr1_ch_s_pt"), ptJet); + registry.fill(HIST("hr1_ch_s_kt"), kt_p); + registry.fill(HIST("hr1_ch_s_z"), z); + registry.fill(HIST("hr1_ch_s_ft"), log10(fT)); + registry.fill(HIST("hr1_ch_s_mult"), ch_mult); + registry.fill(HIST("hr1_ch_s_zg"), zg); + registry.fill(HIST("hr1_ch_s_rg"), rg); + registry.fill(HIST("hr1_ch_s_lu"), log(1 / rg), log(kt_p)); + } + if (ch_l != ch_nl) { + registry.fill(HIST("hr1_ch_d_pt"), ptJet); + registry.fill(HIST("hr1_ch_d_kt"), kt_p); + registry.fill(HIST("hr1_ch_d_z"), z); + registry.fill(HIST("hr1_ch_d_ft"), log10(fT)); + registry.fill(HIST("hr1_ch_d_mult"), ch_mult); + registry.fill(HIST("hr1_ch_d_zg"), zg); + registry.fill(HIST("hr1_ch_d_rg"), rg); + registry.fill(HIST("hr1_ch_d_lu"), log(1 / rg), log(kt_p)); + } + } + //////// for 2nd matched split + if (nsd == 1) { + if (ch_l == ch_nl) { + registry.fill(HIST("hr2_ch_s_pt"), ptJet); + registry.fill(HIST("hr2_ch_s_kt"), kt_p); + registry.fill(HIST("hr2_ch_s_z"), z); + registry.fill(HIST("hr2_ch_s_ft"), log10(fT)); + registry.fill(HIST("hr2_ch_s_mult"), ch_mult); + registry.fill(HIST("hr2_ch_s_zg"), zg); + registry.fill(HIST("hr2_ch_s_rg"), rg); + registry.fill(HIST("hr2_ch_s_lu"), log(1 / rg), log(kt_p)); + } + if (ch_l != ch_nl) { + registry.fill(HIST("hr2_ch_d_pt"), ptJet); + registry.fill(HIST("hr2_ch_d_kt"), kt_p); + registry.fill(HIST("hr2_ch_d_z"), z); + registry.fill(HIST("hr2_ch_d_ft"), log10(fT)); + registry.fill(HIST("hr2_ch_d_mult"), ch_mult); + registry.fill(HIST("hr2_ch_d_zg"), zg); + registry.fill(HIST("hr2_ch_d_rg"), rg); + registry.fill(HIST("hr2_ch_d_lu"), log(1 / rg), log(kt_p)); + } + } + //////// for 3rd marched split + if (nsd >= 2) { + if (ch_l == ch_nl) { + registry.fill(HIST("hr3_ch_s_pt"), ptJet); + registry.fill(HIST("hr3_ch_s_kt"), kt_p); + registry.fill(HIST("hr3_ch_s_z"), z); + registry.fill(HIST("hr3_ch_s_ft"), log10(fT)); + registry.fill(HIST("hr3_ch_s_mult"), ch_mult); + registry.fill(HIST("hr3_ch_s_zg"), zg); + registry.fill(HIST("hr3_ch_s_rg"), rg); + registry.fill(HIST("hr3_ch_s_lu"), log(1 / rg), log(kt_p)); + } + if (ch_l != ch_nl) { + registry.fill(HIST("hr3_ch_d_pt"), ptJet); + registry.fill(HIST("hr3_ch_d_kt"), kt_p); + registry.fill(HIST("hr3_ch_d_z"), z); + registry.fill(HIST("hr3_ch_d_ft"), log10(fT)); + registry.fill(HIST("hr3_ch_d_mult"), ch_mult); + registry.fill(HIST("hr3_ch_d_zg"), zg); + registry.fill(HIST("hr3_ch_d_rg"), rg); + registry.fill(HIST("hr3_ch_d_lu"), log(1 / rg), log(kt_p)); + } + } + } + } nsd++; } daughterSubJet = parentSubJet1; @@ -336,62 +334,61 @@ struct JetChCorr { { jetConstituents.clear(); - //int leadpt =0; - //int n_leadpt =0; - int nn=0; - ch_mult = 0;// jet.size(); - ch_mult = jet.tracksIds().size(); - //for (auto& jetConstituent : jet.template tracks_as()) { - for (auto& jetConstituent : jet.template tracks_as()){ + // int leadpt =0; + // int n_leadpt =0; + int nn = 0; + ch_mult = 0; // jet.size(); + ch_mult = jet.tracksIds().size(); + // for (auto& jetConstituent : jet.template tracks_as()) { + for (auto& jetConstituent : jet.template tracks_as()) { fastjetutilities::fillTracks(jetConstituent, jetConstituents, jetConstituent.globalIndex()); - - //int tr = jetConstituent.globalIndex(); - //auto const& trackSel = tracks.iteratorAt(tr); + // int tr = jetConstituent.globalIndex(); + // auto const& trackSel = tracks.iteratorAt(tr); TVector3 p_leading(jetConstituent.px(), jetConstituent.py(), jetConstituent.pz()); - if(nn==0){ - //leadpt = p_leading.Perp(); - trackL = jetConstituent.globalIndex(); - //ch_l = trackSel.sign(); - ch_l = jetConstituent.sign();//signed1Pt();//pt(); - v1.SetXYZ(jetConstituent.pt(), jetConstituent.py(), jetConstituent.pz()); + if (nn == 0) { + // leadpt = p_leading.Perp(); + trackL = jetConstituent.globalIndex(); + // ch_l = trackSel.sign(); + ch_l = jetConstituent.sign(); // signed1Pt();//pt(); + v1.SetXYZ(jetConstituent.pt(), jetConstituent.py(), jetConstituent.pz()); } - if(nn==1){ - //n_leadpt = p_leading.Perp(); - n_trackL = jetConstituent.globalIndex(); - ch_nl = jetConstituent.sign();//signed1Pt();//pt(); - v2.SetXYZ(jetConstituent.px(), jetConstituent.py(), jetConstituent.pz()); - - vR = v1 + v2; - float z = v2.Perp(vR.Orthogonal())/(v1.Perp(vR.Orthogonal())+v2.Perp(vR.Orthogonal())); - float fT = ((2.*z*(1-z)*vR.Mag())/v1.Perp2(vR))/6.; - float kt_p = v1.Perp(vR); - // float th_p = v1.Angle(v2); - // cout<1) jetReclustering(jet, outputTable); + if (nn > 1) + jetReclustering(jet, outputTable); } void processDummy(JetTracks const& tracks) @@ -399,7 +396,6 @@ struct JetChCorr { } PROCESS_SWITCH(JetChCorr, processDummy, "Dummy process function turned on by default", true); - void processChargedJetsData(soa::Join::iterator const& jet, JetTracks const& tracks) // void processChargedJetsData(soa::Join::iterator const& jet, soa::Join const& tracks) { @@ -438,5 +434,5 @@ WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) { return WorkflowSpec{adaptAnalysisTask( - cfgc, TaskName{"jet-ch-corr"})}; + cfgc, TaskName{"jet-ch-corr"})}; }