From fb8a2351d7a8ddcea98bcb6e565c43b88319b93a Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 9 May 2024 05:55:39 +0000 Subject: [PATCH] Please consider the following formatting changes --- PWGJE/Tasks/phiInJets.cxx | 57 ++++++++++++++++++--------------------- 1 file changed, 26 insertions(+), 31 deletions(-) diff --git a/PWGJE/Tasks/phiInJets.cxx b/PWGJE/Tasks/phiInJets.cxx index 52611da091d..29f816014ec 100644 --- a/PWGJE/Tasks/phiInJets.cxx +++ b/PWGJE/Tasks/phiInJets.cxx @@ -132,10 +132,10 @@ struct phiInJets { JEhistos.add("nEvents_MCGen", "nEvents_MCGen", kTH1F, {{4, 0.0, 4.0}}); JEhistos.add("nEvents_MCRec_MATCHED", "nEvents_MCRec_MATCHED", kTH1F, {{4, 0.0, 4.0}}); JEhistos.add("nEvents_MCGen_MATCHED", "nEvents_MCGen_MATCHED", kTH1F, {{4, 0.0, 4.0}}); - - JEhistos.add("JetVsPhi_GEN","JetVsPhi_GEN",kTH2F,{{4000, 0., 200.},{200, 0, 20.0}}); - JEhistos.add("JetVsPhi_REC","JetVsPhi_REC",kTH2F,{{4000, 0., 200.},{200, 0, 20.0}}); - + + JEhistos.add("JetVsPhi_GEN", "JetVsPhi_GEN", kTH2F, {{4000, 0., 200.}, {200, 0, 20.0}}); + JEhistos.add("JetVsPhi_REC", "JetVsPhi_REC", kTH2F, {{4000, 0., 200.}, {200, 0, 20.0}}); + JEhistos.add("nJetsPerEvent", "nJetsPerEvent", kTH1F, {{10, 0.0, 10.0}}); JEhistos.add("hDCArToPv", "DCArToPv", kTH1F, {{300, 0.0, 3.0}}); @@ -462,7 +462,7 @@ struct phiInJets { if (!trackSelection(originalTrack)) continue; - + JEhistos.fill(HIST("etaHistogram"), trackC.eta()); JEhistos.fill(HIST("phiHistogram"), trackC.phi()); } // JTrack Loop @@ -625,13 +625,12 @@ struct phiInJets { JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_OUTSIDE"), 1.0, lResonance.Pt(), lResonance.M()); } //! jetflag - if (hasJets){ + if (hasJets) { JEhistos.fill(HIST("ptJEHistogramPhi_JetTrigger"), lResonance.Pt()); - auto triggerjet = std::min_element(mcd_pt.begin(), mcd_pt.end()); - double triggerjet_pt = *triggerjet; - JEhistos.fill(HIST("JetVsPhi_REC"), triggerjet_pt, lResonance.Pt()); - - } + auto triggerjet = std::min_element(mcd_pt.begin(), mcd_pt.end()); + double triggerjet_pt = *triggerjet; + JEhistos.fill(HIST("JetVsPhi_REC"), triggerjet_pt, lResonance.Pt()); + } JEhistos.fill(HIST("minvJEHistogramPhi"), lResonance.M()); } // mcpart check } // tracks2 @@ -644,8 +643,8 @@ struct phiInJets { //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// int nprocessSimEvents = 0; // Preslice slice = o2::aod::JCollision::collisionId; - void processSim(o2::aod::JMcCollision const& collision, soa::SmallGroups> const& recocolls, aod::JMcParticles const& mcParticles, soa::Filtered const& mcpjets) - { + void processSim(o2::aod::JMcCollision const& collision, soa::SmallGroups> const& recocolls, aod::JMcParticles const& mcParticles, soa::Filtered const& mcpjets) + { if (cDebugLevel > 0) { nprocessSimEvents++; if ((nprocessSimEvents + 1) % 10000 == 0) @@ -654,17 +653,15 @@ struct phiInJets { JEhistos.fill(HIST("nEvents_MCGen"), 0.5); - if(recocolls.size()<=0)//not reconstructed + if (recocolls.size() <= 0) // not reconstructed return; - - for(auto& recocoll : recocolls){ //poorly reconstructed + + for (auto& recocoll : recocolls) { // poorly reconstructed if (!jetderiveddatautilities::selectCollision(recocoll, jetderiveddatautilities::JCollisionSel::sel8)) - return; + return; } - - - if (fabs(collision.posZ()) > 10)//bad vertex + if (fabs(collision.posZ()) > 10) // bad vertex return; bool INELgt0 = false; for (const auto& mcParticle : mcParticles) { @@ -673,10 +670,9 @@ struct phiInJets { break; } } - if (!INELgt0) //not INEL + if (!INELgt0) // not INEL return; - std::vector mcp_pt{}; std::vector mcp_phi{}; std::vector mcp_eta{}; @@ -753,9 +749,9 @@ struct phiInJets { ////////////////////////////Phi found if (hasJets) { JEhistos.fill(HIST("ptGeneratedPhi_JetTrigger"), mcParticle.pt()); - auto triggerjet = std::min_element(mcp_pt.begin(), mcp_pt.end()); - double triggerjet_pt = *triggerjet; - JEhistos.fill(HIST("JetVsPhi_GEN"), triggerjet_pt, mcParticle.pt()); + auto triggerjet = std::min_element(mcp_pt.begin(), mcp_pt.end()); + double triggerjet_pt = *triggerjet; + JEhistos.fill(HIST("JetVsPhi_GEN"), triggerjet_pt, mcParticle.pt()); } // check for jets } // check for phi @@ -773,7 +769,7 @@ struct phiInJets { // void processMatchedGen(o2::aod::JMcCollision const& collision, aod::JMcParticles const& mcParticles, soa::Filtered const& mcpjets) int nprocessSimJEEvents = 0; void processMatchedGen(aod::JMcCollision const& collision, - soa::SmallGroups> const& recocolls, + soa::SmallGroups> const& recocolls, JetMCDTable const& mcdjets, JetMCPTable const& mcpjets, aod::JMcParticles const& mcParticles) @@ -790,14 +786,13 @@ struct phiInJets { if (fabs(collision.posZ()) > 10) return; - if(recocolls.size()<=0)//not reconstructed + if (recocolls.size() <= 0) // not reconstructed return; - for(auto& recocoll : recocolls){ //poorly reconstructed + for (auto& recocoll : recocolls) { // poorly reconstructed if (!jetderiveddatautilities::selectCollision(recocoll, jetderiveddatautilities::JCollisionSel::sel8)) - return; + return; } - bool INELgt0 = false; for (const auto& mcParticle : mcParticles) { if (TMath::Abs(mcParticle.eta()) < cfgtrkMaxEta) { @@ -847,7 +842,7 @@ struct phiInJets { for (const auto& mcParticle : mcParticles) { if (fabs(mcParticle.eta()) > cfgtrkMaxEta) continue; - + if (fabs(mcParticle.pdgCode()) == 333) { TLorentzVector lResonance;