Skip to content

Commit a6f5cda

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent ddf3b45 commit a6f5cda

File tree

1 file changed

+105
-105
lines changed

1 file changed

+105
-105
lines changed

‎PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx‎

Lines changed: 105 additions & 105 deletions
Original file line numberDiff line numberDiff line change
@@ -601,59 +601,59 @@ struct ThreeParticleCorrelations {
601601

602602
// Start of the MC Track QA
603603
for (const auto& track : groupMCAssociates) {
604-
if (track.isPhysicalPrimary()) {
605-
606-
if (track.pdgCode() > 0) {
607-
assocSign = 1;
608-
} else if (track.pdgCode() < 0) {
609-
assocSign = -1;
610-
}
611-
612-
if (std::abs(track.pdgCode()) == kPiPlus) { // Pions
613-
rQARegistry.fill(HIST("hPtPion_MC"), track.pt(), collision.bestCollisionCentFT0C(), assocSign);
614-
} else if (std::abs(track.pdgCode()) == kKPlus) { // Kaons
615-
rQARegistry.fill(HIST("hPtKaon_MC"), track.pt(), collision.bestCollisionCentFT0C(), assocSign);
616-
} else if (std::abs(track.pdgCode()) == kProton) { // Protons
617-
rQARegistry.fill(HIST("hPtProton_MC"), track.pt(), collision.bestCollisionCentFT0C(), assocSign);
618-
}
619-
}
604+
if (track.isPhysicalPrimary()) {
605+
606+
if (track.pdgCode() > 0) {
607+
assocSign = 1;
608+
} else if (track.pdgCode() < 0) {
609+
assocSign = -1;
610+
}
611+
612+
if (std::abs(track.pdgCode()) == kPiPlus) { // Pions
613+
rQARegistry.fill(HIST("hPtPion_MC"), track.pt(), collision.bestCollisionCentFT0C(), assocSign);
614+
} else if (std::abs(track.pdgCode()) == kKPlus) { // Kaons
615+
rQARegistry.fill(HIST("hPtKaon_MC"), track.pt(), collision.bestCollisionCentFT0C(), assocSign);
616+
} else if (std::abs(track.pdgCode()) == kProton) { // Protons
617+
rQARegistry.fill(HIST("hPtProton_MC"), track.pt(), collision.bestCollisionCentFT0C(), assocSign);
618+
}
619+
}
620620
}
621621
// End of the MC Track QA
622622

623623
// Start of the MC Same-Event correlations
624624
for (const auto& trigger : groupMCTriggers) {
625-
if (trigger.isPhysicalPrimary()) {
626-
627-
if (trigger.pdgCode() > 0) {
628-
triggSign = 1;
629-
} else if (trigger.pdgCode() < 0) {
630-
triggSign = -1;
631-
}
632-
rQARegistry.fill(HIST("hPtV0_MC"), trigger.pt(), collision.bestCollisionCentFT0C(), triggSign);
633-
rQARegistry.fill(HIST("hNLambdas"), triggSign, trigger.pt(), collision.bestCollisionCentFT0C());
634-
635-
for (const auto& associate : groupMCAssociates) {
636-
if (associate.isPhysicalPrimary()) {
637-
638-
if (associate.pdgCode() > 0) {
639-
assocSign = 1;
640-
} else if (associate.pdgCode() < 0) {
641-
assocSign = -1;
642-
}
643-
644-
deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf);
645-
deltaEta = trigger.eta() - associate.eta();
646-
647-
if (std::abs(associate.pdgCode()) == kPiPlus) {
648-
rSECorrRegistry.fill(HIST("hSameLambdaPion_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
649-
} else if (std::abs(associate.pdgCode()) == kKPlus) {
650-
rSECorrRegistry.fill(HIST("hSameLambdaKaon_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
651-
} else if (std::abs(associate.pdgCode()) == kProton) {
652-
rSECorrRegistry.fill(HIST("hSameLambdaProton_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
653-
}
654-
}
655-
}
656-
}
625+
if (trigger.isPhysicalPrimary()) {
626+
627+
if (trigger.pdgCode() > 0) {
628+
triggSign = 1;
629+
} else if (trigger.pdgCode() < 0) {
630+
triggSign = -1;
631+
}
632+
rQARegistry.fill(HIST("hPtV0_MC"), trigger.pt(), collision.bestCollisionCentFT0C(), triggSign);
633+
rQARegistry.fill(HIST("hNLambdas"), triggSign, trigger.pt(), collision.bestCollisionCentFT0C());
634+
635+
for (const auto& associate : groupMCAssociates) {
636+
if (associate.isPhysicalPrimary()) {
637+
638+
if (associate.pdgCode() > 0) {
639+
assocSign = 1;
640+
} else if (associate.pdgCode() < 0) {
641+
assocSign = -1;
642+
}
643+
644+
deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf);
645+
deltaEta = trigger.eta() - associate.eta();
646+
647+
if (std::abs(associate.pdgCode()) == kPiPlus) {
648+
rSECorrRegistry.fill(HIST("hSameLambdaPion_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
649+
} else if (std::abs(associate.pdgCode()) == kKPlus) {
650+
rSECorrRegistry.fill(HIST("hSameLambdaKaon_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
651+
} else if (std::abs(associate.pdgCode()) == kProton) {
652+
rSECorrRegistry.fill(HIST("hSameLambdaProton_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
653+
}
654+
}
655+
}
656+
}
657657
}
658658
// End of the MC Same-Event Correlations
659659
}
@@ -668,35 +668,35 @@ struct ThreeParticleCorrelations {
668668
auto recCollsA2 = recCollisions.sliceBy(perMCCol, coll_2.globalIndex());
669669
if (recCollsA1.size() == 1 && recCollsA2.size() == 1 && acceptEvent(recCollsA1.begin(), false) && acceptEvent(recCollsA2.begin(), false)) {
670670

671-
LOGF(info, "Size_1 = %i, Size_2 = %i", recCollsA1.size(), recCollsA1.size());
672-
auto groupMCTriggers = mcTriggers->sliceByCached(aod::mcparticle::mcCollisionId, coll_1.globalIndex(), cache);
673-
auto groupMCAssociates = mcAssociates->sliceByCached(aod::mcparticle::mcCollisionId, coll_2.globalIndex(), cache);
674-
for (const auto& [trigger, associate] : soa::combinations(soa::CombinationsFullIndexPolicy(groupMCTriggers, groupMCAssociates))) {
675-
if (trigger.isPhysicalPrimary() && associate.isPhysicalPrimary()) {
676-
677-
if (trigger.pdgCode() > 0) {
678-
triggSign = 1;
679-
} else if (trigger.pdgCode() < 0) {
680-
triggSign = -1;
681-
}
682-
if (associate.pdgCode() > 0) {
683-
assocSign = 1;
684-
} else if (associate.pdgCode() < 0) {
685-
assocSign = -1;
686-
}
687-
688-
deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf);
689-
deltaEta = trigger.eta() - associate.eta();
690-
691-
if (std::abs(associate.pdgCode()) == kPiPlus) {
692-
rMECorrRegistry.fill(HIST("hMixLambdaPion_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
693-
} else if (std::abs(associate.pdgCode()) == kKPlus) {
694-
rMECorrRegistry.fill(HIST("hMixLambdaKaon_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
695-
} else if (std::abs(associate.pdgCode()) == kProton) {
696-
rMECorrRegistry.fill(HIST("hMixLambdaProton_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
697-
}
698-
}
699-
}
671+
LOGF(info, "Size_1 = %i, Size_2 = %i", recCollsA1.size(), recCollsA1.size());
672+
auto groupMCTriggers = mcTriggers->sliceByCached(aod::mcparticle::mcCollisionId, coll_1.globalIndex(), cache);
673+
auto groupMCAssociates = mcAssociates->sliceByCached(aod::mcparticle::mcCollisionId, coll_2.globalIndex(), cache);
674+
for (const auto& [trigger, associate] : soa::combinations(soa::CombinationsFullIndexPolicy(groupMCTriggers, groupMCAssociates))) {
675+
if (trigger.isPhysicalPrimary() && associate.isPhysicalPrimary()) {
676+
677+
if (trigger.pdgCode() > 0) {
678+
triggSign = 1;
679+
} else if (trigger.pdgCode() < 0) {
680+
triggSign = -1;
681+
}
682+
if (associate.pdgCode() > 0) {
683+
assocSign = 1;
684+
} else if (associate.pdgCode() < 0) {
685+
assocSign = -1;
686+
}
687+
688+
deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf);
689+
deltaEta = trigger.eta() - associate.eta();
690+
691+
if (std::abs(associate.pdgCode()) == kPiPlus) {
692+
rMECorrRegistry.fill(HIST("hMixLambdaPion_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
693+
} else if (std::abs(associate.pdgCode()) == kKPlus) {
694+
rMECorrRegistry.fill(HIST("hMixLambdaKaon_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
695+
} else if (std::abs(associate.pdgCode()) == kProton) {
696+
rMECorrRegistry.fill(HIST("hMixLambdaProton_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
697+
}
698+
}
699+
}
700700
}
701701
}
702702
// End of the MC Mixed-events Correlations
@@ -712,41 +712,41 @@ struct ThreeParticleCorrelations {
712712

713713
// Start of the Monte-Carlo generated QA
714714
for (const auto& particle : groupMCTracks) {
715-
if (particle.isPhysicalPrimary()) {
716-
717-
// Track efficiency - Generated
718-
rMCRegistry.fill(HIST("hGenerated"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
719-
if (particle.pdgCode() == kPiPlus) { // Pos pions
720-
rMCRegistry.fill(HIST("hGenPionP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
721-
} else if (particle.pdgCode() == kPiMinus) { // Neg pions
722-
rMCRegistry.fill(HIST("hGenPionN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
723-
} else if (particle.pdgCode() == kKPlus) { // Pos kaons
724-
rMCRegistry.fill(HIST("hGenKaonP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
725-
} else if (particle.pdgCode() == kKMinus) { // Neg kaons
726-
rMCRegistry.fill(HIST("hGenKaonN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
727-
} else if (particle.pdgCode() == kProton) { // Pos protons
728-
rMCRegistry.fill(HIST("hGenProtonP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
729-
} else if (particle.pdgCode() == kProtonBar) { // Neg protons
730-
rMCRegistry.fill(HIST("hGenProtonN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
731-
}
732-
}
715+
if (particle.isPhysicalPrimary()) {
716+
717+
// Track efficiency - Generated
718+
rMCRegistry.fill(HIST("hGenerated"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
719+
if (particle.pdgCode() == kPiPlus) { // Pos pions
720+
rMCRegistry.fill(HIST("hGenPionP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
721+
} else if (particle.pdgCode() == kPiMinus) { // Neg pions
722+
rMCRegistry.fill(HIST("hGenPionN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
723+
} else if (particle.pdgCode() == kKPlus) { // Pos kaons
724+
rMCRegistry.fill(HIST("hGenKaonP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
725+
} else if (particle.pdgCode() == kKMinus) { // Neg kaons
726+
rMCRegistry.fill(HIST("hGenKaonN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
727+
} else if (particle.pdgCode() == kProton) { // Pos protons
728+
rMCRegistry.fill(HIST("hGenProtonP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
729+
} else if (particle.pdgCode() == kProtonBar) { // Neg protons
730+
rMCRegistry.fill(HIST("hGenProtonN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
731+
}
732+
}
733733
}
734734

735735
for (const auto& particle : groupMCV0s) {
736-
if (particle.isPhysicalPrimary()) {
737-
738-
// V0 efficiency - Generated
739-
if (particle.pdgCode() == kLambda0) { // Lambdas
740-
rMCRegistry.fill(HIST("hGenLambdaP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
741-
} else if (particle.pdgCode() == kLambda0Bar) { // AntiLambdas
742-
rMCRegistry.fill(HIST("hGenLambdaN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
743-
}
744-
}
736+
if (particle.isPhysicalPrimary()) {
737+
738+
// V0 efficiency - Generated
739+
if (particle.pdgCode() == kLambda0) { // Lambdas
740+
rMCRegistry.fill(HIST("hGenLambdaP"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
741+
} else if (particle.pdgCode() == kLambda0Bar) { // AntiLambdas
742+
rMCRegistry.fill(HIST("hGenLambdaN"), particle.pt(), particle.eta(), collision.bestCollisionCentFT0C());
743+
}
744+
}
745745
}
746746
// End of the Monte-Carlo generated QA
747747
}
748748
}
749-
749+
750750
void processMCRec(MyFilteredMCRecCollisions::iterator const& collision, MyMCV0s const& v0s, MyFilteredMCTracks const& tracks, aod::McCollisions const&, aod::McParticles const&)
751751
{
752752

0 commit comments

Comments
 (0)