From d5e1c729edbe44e16376038a7858f987ba94436f Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Mon, 24 Feb 2025 10:49:24 +0000 Subject: [PATCH] Please consider the following formatting changes --- .../Tasks/threeParticleCorrelations.cxx | 130 +++++++++--------- 1 file changed, 65 insertions(+), 65 deletions(-) diff --git a/PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx b/PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx index 1c0d50ecf87..d0746966a3d 100644 --- a/PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx +++ b/PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx @@ -68,7 +68,7 @@ struct ThreeParticleCorrelations { // Table aliases - MC using MyFilteredMCGenCollisions = soa::Filtered>; - using MyFilteredMCGenCollision = MyFilteredMCGenCollisions::iterator; + using MyFilteredMCGenCollision = MyFilteredMCGenCollisions::iterator; using MyFilteredMCParticles = soa::Filtered; using MyFilteredMCRecCollision = soa::Filtered>::iterator; using MyFilteredMCTracks = soa::Filtered mcTriggers = (aod::mcparticle::pdgCode == static_cast(kLambda0) || aod::mcparticle::pdgCode == static_cast(kLambda0Bar)) && aod::mcparticle::pt > 0.6f && aod::mcparticle::pt < 12.0f && nabs(aod::mcparticle::eta) < 0.72f; + Partition mcTriggers = (aod::mcparticle::pdgCode == static_cast(kLambda0) || aod::mcparticle::pdgCode == static_cast(kLambda0Bar)) && aod::mcparticle::pt > 0.6f && aod::mcparticle::pt < 12.0f && nabs(aod::mcparticle::eta) < 0.72f; Partition mcAssociates = (((aod::mcparticle::pdgCode == static_cast(kPiPlus) || aod::mcparticle::pdgCode == static_cast(kPiMinus)) && aod::mcparticle::pt > 0.3f && aod::mcparticle::pt < 2.3f) || - ((aod::mcparticle::pdgCode == static_cast(kKPlus) || aod::mcparticle::pdgCode == static_cast(kKMinus)) && aod::mcparticle::pt > 0.5f && aod::mcparticle::pt < 2.5f) || - ((aod::mcparticle::pdgCode == static_cast(kProton) || aod::mcparticle::pdgCode == static_cast(kProtonBar)) && aod::mcparticle::pt > 0.5f)); + ((aod::mcparticle::pdgCode == static_cast(kKPlus) || aod::mcparticle::pdgCode == static_cast(kKMinus)) && aod::mcparticle::pt > 0.5f && aod::mcparticle::pt < 2.5f) || + ((aod::mcparticle::pdgCode == static_cast(kProton) || aod::mcparticle::pdgCode == static_cast(kProtonBar)) && aod::mcparticle::pt > 0.5f)); mcTriggers.bindTable(particles); mcAssociates.bindTable(particles); - + // Start of the MC Same-Event Correlations for (const auto& trigger : mcTriggers) { if (trigger.isPhysicalPrimary()) { - - if (trigger.pdgCode() > 0) { - triggSign = 1; - rQARegistry.fill(HIST("hNLambdas_MC"), 1, trigger.pt(), collision.centFT0M()); - } else if (trigger.pdgCode() < 0) { - triggSign = -1; - rQARegistry.fill(HIST("hNLambdas_MC"), -1, trigger.pt(), collision.centFT0M()); - } - - for (const auto& associate : mcAssociates) { - if (associate.isPhysicalPrimary()) { - - if (associate.pdgCode() > 0) { - assocSign = 1; - } else if (associate.pdgCode() < 0) { - assocSign = -1; - } - - deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf); - deltaEta = trigger.eta() - associate.eta(); - - if (std::abs(associate.pdgCode()) == kPiPlus) { - rSECorrRegistry.fill(HIST("hSameLambdaPion_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign); - } else if (std::abs(associate.pdgCode()) == kKPlus) { - rSECorrRegistry.fill(HIST("hSameLambdaKaon_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign); - } else if (std::abs(associate.pdgCode()) == kProton) { - rSECorrRegistry.fill(HIST("hSameLambdaProton_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign); - } - } - } + + if (trigger.pdgCode() > 0) { + triggSign = 1; + rQARegistry.fill(HIST("hNLambdas_MC"), 1, trigger.pt(), collision.centFT0M()); + } else if (trigger.pdgCode() < 0) { + triggSign = -1; + rQARegistry.fill(HIST("hNLambdas_MC"), -1, trigger.pt(), collision.centFT0M()); + } + + for (const auto& associate : mcAssociates) { + if (associate.isPhysicalPrimary()) { + + if (associate.pdgCode() > 0) { + assocSign = 1; + } else if (associate.pdgCode() < 0) { + assocSign = -1; + } + + deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf); + deltaEta = trigger.eta() - associate.eta(); + + if (std::abs(associate.pdgCode()) == kPiPlus) { + rSECorrRegistry.fill(HIST("hSameLambdaPion_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign); + } else if (std::abs(associate.pdgCode()) == kKPlus) { + rSECorrRegistry.fill(HIST("hSameLambdaKaon_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign); + } else if (std::abs(associate.pdgCode()) == kProton) { + rSECorrRegistry.fill(HIST("hSameLambdaProton_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign); + } + } + } } } // End of the MC Same-Event Correlations @@ -407,49 +407,49 @@ struct ThreeParticleCorrelations { // Start of the MC Mixed-events Correlations for (const auto& [coll_1, v0_1, coll_2, track_2] : pairMC) { - Partition mcTriggers = (aod::mcparticle::pdgCode == static_cast(kLambda0) || aod::mcparticle::pdgCode == static_cast(kLambda0Bar)) && aod::mcparticle::pt > 0.6f && aod::mcparticle::pt < 12.0f && nabs(aod::mcparticle::eta) < 0.72f; + Partition mcTriggers = (aod::mcparticle::pdgCode == static_cast(kLambda0) || aod::mcparticle::pdgCode == static_cast(kLambda0Bar)) && aod::mcparticle::pt > 0.6f && aod::mcparticle::pt < 12.0f && nabs(aod::mcparticle::eta) < 0.72f; Partition mcAssociates = (((aod::mcparticle::pdgCode == static_cast(kPiPlus) || aod::mcparticle::pdgCode == static_cast(kPiMinus)) && aod::mcparticle::pt > 0.3f && aod::mcparticle::pt < 2.3f) || - ((aod::mcparticle::pdgCode == static_cast(kKPlus) || aod::mcparticle::pdgCode == static_cast(kKMinus)) && aod::mcparticle::pt > 0.5f && aod::mcparticle::pt < 2.5f) || - ((aod::mcparticle::pdgCode == static_cast(kProton) || aod::mcparticle::pdgCode == static_cast(kProtonBar)) && aod::mcparticle::pt > 0.5f)); + ((aod::mcparticle::pdgCode == static_cast(kKPlus) || aod::mcparticle::pdgCode == static_cast(kKMinus)) && aod::mcparticle::pt > 0.5f && aod::mcparticle::pt < 2.5f) || + ((aod::mcparticle::pdgCode == static_cast(kProton) || aod::mcparticle::pdgCode == static_cast(kProtonBar)) && aod::mcparticle::pt > 0.5f)); mcTriggers.bindTable(v0_1); mcAssociates.bindTable(track_2); - + for (const auto& [trigger, associate] : soa::combinations(soa::CombinationsFullIndexPolicy(mcTriggers, mcAssociates))) { - if (trigger.isPhysicalPrimary() && associate.isPhysicalPrimary()) { - - if (trigger.pdgCode() > 0) { - triggSign = 1; - } else if (trigger.pdgCode() < 0) { - triggSign = -1; - } - if (associate.pdgCode() > 0) { - assocSign = 1; - } else if (associate.pdgCode() < 0) { - assocSign = -1; - } - - deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf); - deltaEta = trigger.eta() - associate.eta(); - - if (std::abs(associate.pdgCode()) == kPiPlus) { - rMECorrRegistry.fill(HIST("hMixLambdaPion_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign); - } else if (std::abs(associate.pdgCode()) == kKPlus) { - rMECorrRegistry.fill(HIST("hMixLambdaKaon_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign); - } else if (std::abs(associate.pdgCode()) == kProton) { - rMECorrRegistry.fill(HIST("hMixLambdaProton_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign); - } - } + if (trigger.isPhysicalPrimary() && associate.isPhysicalPrimary()) { + + if (trigger.pdgCode() > 0) { + triggSign = 1; + } else if (trigger.pdgCode() < 0) { + triggSign = -1; + } + if (associate.pdgCode() > 0) { + assocSign = 1; + } else if (associate.pdgCode() < 0) { + assocSign = -1; + } + + deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf); + deltaEta = trigger.eta() - associate.eta(); + + if (std::abs(associate.pdgCode()) == kPiPlus) { + rMECorrRegistry.fill(HIST("hMixLambdaPion_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign); + } else if (std::abs(associate.pdgCode()) == kKPlus) { + rMECorrRegistry.fill(HIST("hMixLambdaKaon_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign); + } else if (std::abs(associate.pdgCode()) == kProton) { + rMECorrRegistry.fill(HIST("hMixLambdaProton_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign); + } + } } } // End of the MC Mixed-events Correlations - } + } void processMCGen(MyFilteredMCGenCollision const&, MyFilteredMCParticles const& particles) { Partition mcParticles = aod::mcparticle::pt > 0.2f && aod::mcparticle::pt < 3.0f; mcParticles.bindTable(particles); - + // Start of the Monte-Carlo generated QA for (const auto& particle : mcParticles) { if (particle.isPhysicalPrimary()) {