Skip to content

Commit aaeb1d4

Browse files
a-m-andrushkoa-m-andrushko
committed
Improve K0Short histogramming
1 parent e326226 commit aaeb1d4

1 file changed

Lines changed: 40 additions & 39 deletions

File tree

‎PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx‎

Lines changed: 40 additions & 39 deletions
Original file line numberDiff line numberDiff line change
@@ -519,47 +519,48 @@ struct FemtoUniversePairTaskTrackV0Helicity {
519519

520520
/// Histogramming for MC Reco to calculate fraction of K0Short in the Lambda sample for each helicity bin
521521
if constexpr (confIsMC) {
522-
if (part.has_fdMCParticle())
522+
if (part.has_fdMCParticle()) {
523523
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
524524
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMC"), part.mLambda());
525-
else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short)
526-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMC"), part.mKaon());
527-
528-
if (cosineTheta <= helicityconfigs.confLimitHel0 && cosineTheta >= helicityconfigs.confLimitHel1) {
529-
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
530-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel1"), part.mLambda());
531-
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
532-
if (helicityconfigs.isRealK0Short)
533-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel1"), part.mKaon());
534-
else
535-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel1"), part.mLambda());
536-
}
537-
} else if (cosineTheta < helicityconfigs.confLimitHel1 && cosineTheta >= helicityconfigs.confLimitHel2) {
538-
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
539-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel2"), part.mLambda());
540-
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
541-
if (helicityconfigs.isRealK0Short)
542-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel2"), part.mKaon());
543-
else
544-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel2"), part.mLambda());
545-
}
546-
} else if (cosineTheta < helicityconfigs.confLimitHel2 && cosineTheta >= helicityconfigs.confLimitHel3) {
547-
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
548-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel3"), part.mLambda());
549-
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
550-
if (helicityconfigs.isRealK0Short)
551-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel3"), part.mKaon());
552-
else
553-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel3"), part.mLambda());
554-
}
555-
} else if (cosineTheta < helicityconfigs.confLimitHel3 && cosineTheta >= helicityconfigs.confLimitHel4) {
556-
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
557-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel4"), part.mLambda());
558-
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
559-
if (helicityconfigs.isRealK0Short)
560-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel4"), part.mKaon());
561-
else
562-
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel4"), part.mLambda());
525+
else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short)
526+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMC"), part.mKaon());
527+
528+
if (cosineTheta <= helicityconfigs.confLimitHel0 && cosineTheta >= helicityconfigs.confLimitHel1) {
529+
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
530+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel1"), part.mLambda());
531+
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
532+
if (helicityconfigs.isRealK0Short)
533+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel1"), part.mKaon());
534+
else
535+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel1"), part.mLambda());
536+
}
537+
} else if (cosineTheta < helicityconfigs.confLimitHel1 && cosineTheta >= helicityconfigs.confLimitHel2) {
538+
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
539+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel2"), part.mLambda());
540+
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
541+
if (helicityconfigs.isRealK0Short)
542+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel2"), part.mKaon());
543+
else
544+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel2"), part.mLambda());
545+
}
546+
} else if (cosineTheta < helicityconfigs.confLimitHel2 && cosineTheta >= helicityconfigs.confLimitHel3) {
547+
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
548+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel3"), part.mLambda());
549+
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
550+
if (helicityconfigs.isRealK0Short)
551+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel3"), part.mKaon());
552+
else
553+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel3"), part.mLambda());
554+
}
555+
} else if (cosineTheta < helicityconfigs.confLimitHel3 && cosineTheta >= helicityconfigs.confLimitHel4) {
556+
if ((part.fdMCParticle()).pdgMCTruth() == kLambda0) {
557+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassLambdaMCHel4"), part.mLambda());
558+
} else if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
559+
if (helicityconfigs.isRealK0Short)
560+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel4"), part.mKaon());
561+
else
562+
thetaRegistry.fill(HIST("Theta/Mother/MCReco/hInvMassK0ShortMCHel4"), part.mLambda());
563+
}
563564
}
564565
}
565566
}

0 commit comments

Comments
 (0)