Skip to content

Commit 6e05d66

Browse files
zchochulZuzannaZuzanna Chochulska
authored
[PWGCF] FemtoUniverse -- corrected mc truth filling for the secondary contamination estimation (#11987)
Co-authored-by: Zuzanna <01150674@pw.edu.pl> Co-authored-by: Zuzanna Chochulska <zchochul@cern.ch>
1 parent 6f75b98 commit 6e05d66

File tree

1 file changed

+5
-15
lines changed

1 file changed

+5
-15
lines changed

PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackPhi.cxx

Lines changed: 5 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -513,13 +513,6 @@ struct FemtoUniversePairTaskTrackPhi {
513513
if constexpr (isMC) {
514514
// reco
515515
effCorrection.fillRecoHist<ParticleNo::ONE, FilteredFDCollisions>(phicandidate, 333);
516-
// truth
517-
auto mcPartId1 = phicandidate.fdMCParticleId();
518-
auto const& mcpart1 = mcParts.iteratorAt(mcPartId1);
519-
if (mcpart1.pdgMCTruth() != 333) {
520-
continue;
521-
}
522-
effCorrection.fillTruthHist<ParticleNo::ONE, FilteredFDCollisions>(phicandidate);
523516
}
524517
}
525518

@@ -560,14 +553,6 @@ struct FemtoUniversePairTaskTrackPhi {
560553

561554
if constexpr (isMC) {
562555
effCorrection.fillRecoHist<ParticleNo::TWO, FilteredFDCollisions>(track, ConfTrackPDGCode);
563-
564-
// truth
565-
auto mcPartId2 = track.fdMCParticleId();
566-
auto const& mcpart2 = mcParts.iteratorAt(mcPartId2);
567-
if (mcpart2.pdgMCTruth() != ConfTrackPDGCode) {
568-
continue;
569-
}
570-
effCorrection.fillTruthHist<ParticleNo::TWO, FilteredFDCollisions>(track);
571556
}
572557
}
573558

@@ -728,6 +713,10 @@ struct FemtoUniversePairTaskTrackPhi {
728713
continue;
729714
}
730715

716+
if (pdgCode == ConfTrackPDGCode) {
717+
effCorrection.fillTruthHist<ParticleNo::TWO, FilteredFDCollisions>(part);
718+
}
719+
731720
// charge +
732721
if (pdgParticle->Charge() > 0.0) {
733722
registryMCtruth.fill(HIST("MCtruthAllPositivePt"), part.pt());
@@ -745,6 +734,7 @@ struct FemtoUniversePairTaskTrackPhi {
745734
if (pdgCode == 333) {
746735
registryMCtruth.fill(HIST("MCtruthPhi"), part.pt(), part.eta());
747736
registryMCtruth.fill(HIST("MCtruthPhiPt"), part.pt());
737+
effCorrection.fillTruthHist<ParticleNo::ONE, FilteredFDCollisions>(part);
748738
continue;
749739
}
750740

0 commit comments

Comments
 (0)