Skip to content

Commit 8972c5f

Browse files
author
Shirajum Monira
committed
fixed bugs
1 parent cd7b04c commit 8972c5f

File tree

1 file changed

+21
-20
lines changed

1 file changed

+21
-20
lines changed

‎PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Extended.cxx‎

Lines changed: 21 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -506,60 +506,61 @@ struct FemtoUniversePairTaskTrackV0Extended {
506506
for (const auto& part : groupPartsTwo) {
507507
if (!invMLambda(part.mLambda(), part.mAntiLambda()))
508508
continue;
509-
const auto& posChild = parts.iteratorAt(part.globalIndex() - 2 - parts.begin().globalIndex());
510-
const auto& negChild = parts.iteratorAt(part.globalIndex() - 1 - parts.begin().globalIndex());
511-
509+
const auto& posChild1 = parts.iteratorAt(part.globalIndex() - 2 - parts.begin().globalIndex());
510+
const auto& negChild1 = parts.iteratorAt(part.globalIndex() - 1 - parts.begin().globalIndex());
512511
/// Check daughters of first V0 particle
513512
if constexpr (std::experimental::is_detected<hasSigma, typename PartType::iterator>::value) {
514-
if (!isParticleTPC(posChild, V0ChildTable[ConfV0Selection.confV0Type1][0]) || !isParticleTPC(negChild, V0ChildTable[ConfV0Selection.confV0Type1][1]))
513+
if (!isParticleTPC(posChild1, V0ChildTable[ConfV0Selection.confV0Type1][0]) || !isParticleTPC(negChild1, V0ChildTable[ConfV0Selection.confV0Type1][1]))
515514
continue;
516-
if (!isParticleTOF(posChild, V0ChildTable[ConfV0Selection.confV0Type1][0]) || !isParticleTOF(negChild, V0ChildTable[ConfV0Selection.confV0Type1][1]))
515+
if (!isParticleTOF(posChild1, V0ChildTable[ConfV0Selection.confV0Type1][0]) || !isParticleTOF(negChild1, V0ChildTable[ConfV0Selection.confV0Type1][1]))
517516
continue;
518517
trackHistoV0Type1.fillQABase<false, true>(part, HIST("V0Type1"));
519-
posChildV0Type1.fillQABase<false, true>(posChild, HIST("posChildV0Type1"));
520-
negChildV0Type1.fillQABase<false, true>(negChild, HIST("negChildV0Type1"));
518+
posChildV0Type1.fillQABase<false, true>(posChild1, HIST("posChildV0Type1"));
519+
negChildV0Type1.fillQABase<false, true>(negChild1, HIST("negChildV0Type1"));
521520
} else {
522-
if ((posChild.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type1][0])) == 0 || (negChild.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type1][1])) == 0)
521+
if ((posChild1.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type1][0])) == 0 || (negChild1.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type1][1])) == 0)
523522
continue;
524523
if (ConfV0Selection.confUseStrangenessTOF) {
525524
if (((ConfV0Selection.confV0Type1 == 0) && (part.pidCut() & 3) != 3) || ((ConfV0Selection.confV0Type1 == 1) && (part.pidCut() & 12) != 12) || ((ConfV0Selection.confV0Type1 == 2) && (part.pidCut() & 48) != 48))
526525
continue;
527526
} else {
528-
if ((posChild.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type1][0])) == 0 || (negChild.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type1][1])) == 0)
527+
if ((posChild1.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type1][0])) == 0 || (negChild1.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type1][1])) == 0)
529528
continue;
530529
}
531530
trackHistoV0Type1.fillQABase<false, false>(part, HIST("V0Type1"));
532-
posChildV0Type1.fillQABase<false, false>(posChild, HIST("posChildV0Type1"));
533-
negChildV0Type1.fillQABase<false, false>(negChild, HIST("negChildV0Type1"));
531+
posChildV0Type1.fillQABase<false, false>(posChild1, HIST("posChildV0Type1"));
532+
negChildV0Type1.fillQABase<false, false>(negChild1, HIST("negChildV0Type1"));
534533
}
535534
qaRegistry.fill(HIST("V0Type1/hInvMassLambdaVsCent"), multCol, part.mLambda());
536535
qaRegistry.fill(HIST("V0Type1/hInvMassAntiLambdaVsCent"), multCol, part.mAntiLambda());
537536
if constexpr (isMC) {
538537
effCorrection.fillRecoHist<ParticleNo::ONE, FilteredFDCollisions>(part, kLambda0);
539538
}
540539

540+
const auto& posChild2 = parts.iteratorAt(part.globalIndex() - 2 - parts.begin().globalIndex());
541+
const auto& negChild2 = parts.iteratorAt(part.globalIndex() - 1 - parts.begin().globalIndex());
541542
/// Check daughters of second V0 particle
542543
if constexpr (std::experimental::is_detected<hasSigma, typename PartType::iterator>::value) {
543-
if (!isParticleTPC(posChild, V0ChildTable[ConfV0Selection.confV0Type2][0]) || !isParticleTPC(negChild, V0ChildTable[ConfV0Selection.confV0Type2][1]))
544+
if (!isParticleTPC(posChild2, V0ChildTable[ConfV0Selection.confV0Type2][0]) || !isParticleTPC(negChild2, V0ChildTable[ConfV0Selection.confV0Type2][1]))
544545
continue;
545-
if (!isParticleTOF(posChild, V0ChildTable[ConfV0Selection.confV0Type2][0]) || !isParticleTOF(negChild, V0ChildTable[ConfV0Selection.confV0Type2][1]))
546+
if (!isParticleTOF(posChild2, V0ChildTable[ConfV0Selection.confV0Type2][0]) || !isParticleTOF(negChild2, V0ChildTable[ConfV0Selection.confV0Type2][1]))
546547
continue;
547548
trackHistoV0Type2.fillQABase<false, true>(part, HIST("V0Type2"));
548-
posChildV0Type2.fillQABase<false, true>(posChild, HIST("posChildV0Type2"));
549-
negChildV0Type2.fillQABase<false, true>(negChild, HIST("negChildV0Type2"));
549+
posChildV0Type2.fillQABase<false, true>(posChild2, HIST("posChildV0Type2"));
550+
negChildV0Type2.fillQABase<false, true>(negChild2, HIST("negChildV0Type2"));
550551
} else {
551-
if ((posChild.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type1][0])) == 0 || (negChild.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type1][1])) == 0)
552+
if ((posChild2.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type2][0])) == 0 || (negChild2.pidCut() & (1u << V0ChildTable[ConfV0Selection.confV0Type2][1])) == 0)
552553
continue;
553554
if (ConfV0Selection.confUseStrangenessTOF) {
554-
if (((ConfV0Selection.confV0Type1 == 0) && (part.pidCut() & 3) != 3) || ((ConfV0Selection.confV0Type1 == 1) && (part.pidCut() & 12) != 12) || ((ConfV0Selection.confV0Type1 == 2) && (part.pidCut() & 48) != 48))
555+
if (((ConfV0Selection.confV0Type2 == 0) && (part.pidCut() & 3) != 3) || ((ConfV0Selection.confV0Type2 == 1) && (part.pidCut() & 12) != 12) || ((ConfV0Selection.confV0Type2 == 2) && (part.pidCut() & 48) != 48))
555556
continue;
556557
} else {
557-
if ((posChild.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type1][0])) == 0 || (negChild.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type1][1])) == 0)
558+
if ((posChild2.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type2][0])) == 0 || (negChild2.pidCut() & (8u << V0ChildTable[ConfV0Selection.confV0Type2][1])) == 0)
558559
continue;
559560
}
560561
trackHistoV0Type2.fillQABase<false, false>(part, HIST("V0Type2"));
561-
posChildV0Type2.fillQABase<false, false>(posChild, HIST("posChildV0Type2"));
562-
negChildV0Type2.fillQABase<false, false>(negChild, HIST("negChildV0Type2"));
562+
posChildV0Type2.fillQABase<false, false>(posChild2, HIST("posChildV0Type2"));
563+
negChildV0Type2.fillQABase<false, false>(negChild2, HIST("negChildV0Type2"));
563564
}
564565
qaRegistry.fill(HIST("V0Type2/hInvMassLambdaVsCent"), multCol, part.mLambda());
565566
qaRegistry.fill(HIST("V0Type2/hInvMassAntiLambdaVsCent"), multCol, part.mAntiLambda());

0 commit comments

Comments
 (0)