@@ -262,7 +262,7 @@ struct HfDataCreatorCharmResoReduced {
262262 registry.add (" hMCGenOrigin" , " Origin of Generated particles" , {HistType::kTH1F , {{3 , -0.5 , 2.5 }}});
263263 registry.add (" hMCOriginCounterWrongDecay" , " Origin of Generated particles in Wrong decay" , {HistType::kTH1F , {{3 , -0.5 , 2.5 }}});
264264
265- if (doMcRecQa){
265+ if (doMcRecQa) {
266266 registry.add (" hMassVsPtK0Matched" , " K0s candidates Matched ;#it{p}_{T} (GeV/#it{c});inv. mass (GeV/#it{c}^{2});entries" , {HistType::kTH2F , {axisPt, axisMassKzero}});
267267 registry.add (" hMassVsPtD0Matched" , " D0 candidates Matched ;#it{p}_{T} (GeV/#it{c});inv. mass (GeV/#it{c}^{2});entries" , {HistType::kTH2F , {axisPt, axisMassDplus}});
268268 registry.add (" hMassVsPtDstarMatched" , " Dstar candidates Matched ;#it{p}_{T} (GeV/#it{c});inv. mass (GeV/#it{c}^{2});entries" , {HistType::kTH2F , {axisPt, axisMassDstar}});
@@ -562,7 +562,7 @@ struct HfDataCreatorCharmResoReduced {
562562 indexRecReso = RecoDecay::getMatchedMCRec<false , true , false , true >(particlesMc, std::array{vecDaughtersReso[0 ], vecDaughtersReso[1 ], vecDaughtersReso[2 ], vecDaughtersReso[3 ], vecDaughtersReso[4 ]}, Pdg::kDS1 , std::array{+kPiPlus , -kKPlus , +kPiPlus , +kPiPlus , -kPiPlus }, true , &sign, 3 , &nPiToMuReso, &nKaToPiReso);
563563 if (indexRecReso > -1 && nPiToMuReso == 0 && nKaToPiReso == 0 ) {
564564 flag = sign * DecayTypeMc::Ds1ToDStarK0ToD0PiK0s;
565- } else if (indexRecReso > -1 && nPiToMuReso >= 1 && nKaToPiReso == 0 ){
565+ } else if (indexRecReso > -1 && nPiToMuReso >= 1 && nKaToPiReso == 0 ) {
566566 flag = sign * DecayTypeMc::Ds1ToDStarK0ToD0PiK0sOneMu;
567567 }
568568 }
@@ -600,57 +600,57 @@ struct HfDataCreatorCharmResoReduced {
600600 indexParticle = indexRecReso;
601601 } else if (indexRecResoPartReco > -1 ) {
602602 indexParticle = indexRecResoPartReco;
603- }
603+ }
604604 auto particleReso = particlesMc.iteratorAt (indexParticle);
605605 origin = RecoDecay::getCharmHadronOrigin (particlesMc, particleReso, false , &idxBhadMothers);
606606 motherPt = particleReso.pt ();
607607 }
608- if (doMcRecQa){
609- if (indexRecReso > -1 ){
610- if (nPiToMuReso == 0 && nKaToPiReso == 0 ){
608+ if (doMcRecQa) {
609+ if (indexRecReso > -1 ) {
610+ if (nPiToMuReso == 0 && nKaToPiReso == 0 ) {
611611 registry.fill (HIST (" hMassVsPtDs1Matched" ), varUtils.ptD , varUtils.invMassKPiPiV0 - varUtils.invMassD );
612612 }
613- if (nPiToMuReso >= 1 ){
613+ if (nPiToMuReso >= 1 ) {
614614 registry.fill (HIST (" hMassVsPtDs1MatchedPiToMu" ), varUtils.ptD , varUtils.invMassKPiPiV0 - varUtils.invMassD );
615615 }
616- if (nKaToPiReso >= 1 ){
616+ if (nKaToPiReso >= 1 ) {
617617 registry.fill (HIST (" hMassVsPtDs1MatchedKaToPi" ), varUtils.ptD , varUtils.invMassKPiPiV0 - varUtils.invMassD );
618618 }
619619 }
620- if (indexRecD0 > -1 ){
621- if (nPiToMuD0 == 0 && nKaToPiD0 == 0 ){
620+ if (indexRecD0 > -1 ) {
621+ if (nPiToMuD0 == 0 && nKaToPiD0 == 0 ) {
622622 registry.fill (HIST (" hMassVsPtD0Matched" ), varUtils.ptD , varUtils.invMassDdau );
623623 }
624- if (nPiToMuD0 >= 1 ){
624+ if (nPiToMuD0 >= 1 ) {
625625 registry.fill (HIST (" hMassVsPtD0MatchedPiToMu" ), varUtils.ptD , varUtils.invMassDdau );
626626 }
627- if (nKaToPiD0 >= 1 ){
627+ if (nKaToPiD0 >= 1 ) {
628628 registry.fill (HIST (" hMassVsPtD0MatchedKaToPi" ), varUtils.ptD , varUtils.invMassDdau );
629629 }
630630 }
631- if (indexRecDstar > -1 ){
632- if (nPiToMuDstar == 0 && nKaToPiDstar == 0 ){
631+ if (indexRecDstar > -1 ) {
632+ if (nPiToMuDstar == 0 && nKaToPiDstar == 0 ) {
633633 registry.fill (HIST (" hMassVsPtDstarMatched" ), varUtils.ptD , varUtils.invMassD - varUtils.invMassDdau );
634634 }
635- if (nPiToMuDstar >= 1 ){
635+ if (nPiToMuDstar >= 1 ) {
636636 registry.fill (HIST (" hMassVsPtDstarMatchedPiToMu" ), varUtils.ptD , varUtils.invMassD - varUtils.invMassDdau );
637637 }
638- if (nKaToPiDstar >= 1 ){
638+ if (nKaToPiDstar >= 1 ) {
639639 registry.fill (HIST (" hMassVsPtDstarMatchedKaToPi" ), varUtils.ptD , varUtils.invMassD - varUtils.invMassDdau );
640640 }
641641 }
642- if (indexRecK0 > -1 ){
643- if (nPiToMuV0 == 0 && nKaToPiV0 == 0 ){
644- registry.fill (HIST (" hMassVsPtK0Matched" ), candidateV0.pT , candidateV0.mK0Short );
642+ if (indexRecK0 > -1 ) {
643+ if (nPiToMuV0 == 0 && nKaToPiV0 == 0 ) {
644+ registry.fill (HIST (" hMassVsPtK0Matched" ), candidateV0.pT , candidateV0.mK0Short );
645645 }
646- if (nPiToMuV0 >= 1 ){
646+ if (nPiToMuV0 >= 1 ) {
647647 registry.fill (HIST (" hMassVsPtK0MatchedPiToMu" ), candidateV0.pT , candidateV0.mK0Short );
648648 }
649- if (nKaToPiV0 >= 1 ){
649+ if (nKaToPiV0 >= 1 ) {
650650 registry.fill (HIST (" hMassVsPtK0MatchedKaToPi" ), candidateV0.pT , candidateV0.mK0Short );
651651 }
652652 }
653- }
653+ }
654654 } else if constexpr (decChannel == DecayChannel::DplusV0) {
655655 // Ds2Star → D+ K0 → (π+K-π+) K0s → (π+K-π+)(π+π-)
656656 indexRecK0 = RecoDecay::getMatchedMCRec<false , true , false , true >(particlesMc, std::array{vecDaughtersReso[3 ], vecDaughtersReso[4 ]}, kK0 , std::array{+kPiPlus , -kPiPlus }, true , &signV0, 2 , &nPiToMuV0, &nKaToPiV0);
@@ -659,60 +659,60 @@ struct HfDataCreatorCharmResoReduced {
659659 indexRecReso = RecoDecay::getMatchedMCRec<false , false , false , true >(particlesMc, std::array{vecDaughtersReso[0 ], vecDaughtersReso[1 ], vecDaughtersReso[2 ], vecDaughtersReso[3 ], vecDaughtersReso[4 ]}, Pdg::kDS2Star , std::array{+kPiPlus , -kKPlus , +kPiPlus , +kPiPlus , -kPiPlus }, true , &sign, 3 , &nPiToMuReso, &nKaToPiReso);
660660 if (indexRecReso > -1 && nPiToMuReso == 0 && nKaToPiReso == 0 ) {
661661 flag = sign * DecayTypeMc::Ds2StarToDplusK0sToPiKaPiPiPi;
662- } else if (indexRecReso > -1 && nPiToMuReso >= 1 && nKaToPiReso == 0 ){
662+ } else if (indexRecReso > -1 && nPiToMuReso >= 1 && nKaToPiReso == 0 ) {
663663 flag = sign * DecayTypeMc::Ds2StarToDplusK0sOneMu;
664- } else if (indexRecReso < 0 ) {
664+ } else if (indexRecReso < 0 ) {
665665 // Verify partly reconstructed decay Ds1 -> D* K0s -> D+ π0 K0s
666666 indexRecDstar = RecoDecay::getMatchedMCRec<false , false , true , true >(particlesMc, std::array{vecDaughtersReso[0 ], vecDaughtersReso[1 ], vecDaughtersReso[2 ]}, Pdg::kDStar , std::array{-kKPlus , +kPiPlus , +kPiPlus }, true , &signDStar, 2 );
667- if (indexRecDstar > -1 ){
667+ if (indexRecDstar > -1 ) {
668668 indexRecReso = RecoDecay::getMatchedMCRec<false , true , true , true >(particlesMc, std::array{vecDaughtersReso[0 ], vecDaughtersReso[1 ], vecDaughtersReso[2 ], vecDaughtersReso[3 ], vecDaughtersReso[4 ]}, Pdg::kDS1 , std::array{+kPiPlus , -kKPlus , +kPiPlus , +kPiPlus , -kPiPlus }, true , &sign, 3 );
669- if (indexRecReso > -1 ){
669+ if (indexRecReso > -1 ) {
670670 flag = sign * DecayTypeMc::Ds1ToDStarK0ToDPlusPi0K0s;
671671 }
672672 }
673673 }
674- }
674+ }
675675 if (flag != 0 ) {
676676 auto particleReso = particlesMc.iteratorAt (indexRecReso);
677677 origin = RecoDecay::getCharmHadronOrigin (particlesMc, particleReso, false , &idxBhadMothers);
678678 motherPt = particleReso.pt ();
679679 }
680- if (doMcRecQa){
681- if (indexRecReso > -1 ){
682- if (nPiToMuReso == 0 && nKaToPiReso == 0 ){
680+ if (doMcRecQa) {
681+ if (indexRecReso > -1 ) {
682+ if (nPiToMuReso == 0 && nKaToPiReso == 0 ) {
683683 registry.fill (HIST (" hMassVsPtDs2StarMatched" ), varUtils.ptD , varUtils.invMassKPiPiV0 - varUtils.invMassD );
684684 }
685- if (nPiToMuReso >= 1 ){
685+ if (nPiToMuReso >= 1 ) {
686686 registry.fill (HIST (" hMassVsPtDs2StarMatchedPiToMu" ), varUtils.ptD , varUtils.invMassKPiPiV0 - varUtils.invMassD );
687687 }
688- if (nKaToPiReso >= 1 ){
688+ if (nKaToPiReso >= 1 ) {
689689 registry.fill (HIST (" hMassVsPtDs2StarMatchedKaToPi" ), varUtils.ptD , varUtils.invMassKPiPiV0 - varUtils.invMassD );
690690 }
691691 }
692- if (indexRecDplus > -1 ){
693- if (nPiToMuDplus == 0 && nKaToPiDplus == 0 ){
692+ if (indexRecDplus > -1 ) {
693+ if (nPiToMuDplus == 0 && nKaToPiDplus == 0 ) {
694694 registry.fill (HIST (" hMassVsPtDplusMatched" ), varUtils.ptD , varUtils.invMassD );
695695 }
696- if (nPiToMuDplus >= 1 ){
696+ if (nPiToMuDplus >= 1 ) {
697697 registry.fill (HIST (" hMassVsPtDplusMatchedPiToMu" ), varUtils.ptD , varUtils.invMassD );
698698 }
699- if (nKaToPiDplus >= 1 ){
699+ if (nKaToPiDplus >= 1 ) {
700700 registry.fill (HIST (" hMassVsPtDplusMatchedKaToPi" ), varUtils.ptD , varUtils.invMassD );
701701 }
702702 }
703- if (indexRecK0 > -1 ){
704- if (nPiToMuV0 == 0 && nKaToPiV0 == 0 ){
705- registry.fill (HIST (" hMassVsPtK0Matched" ), candidateV0.pT , candidateV0.mK0Short );
703+ if (indexRecK0 > -1 ) {
704+ if (nPiToMuV0 == 0 && nKaToPiV0 == 0 ) {
705+ registry.fill (HIST (" hMassVsPtK0Matched" ), candidateV0.pT , candidateV0.mK0Short );
706706 }
707- if (nPiToMuV0 >= 1 ){
707+ if (nPiToMuV0 >= 1 ) {
708708 registry.fill (HIST (" hMassVsPtK0MatchedPiToMu" ), candidateV0.pT , candidateV0.mK0Short );
709709 }
710- if (nKaToPiV0 >= 1 ){
710+ if (nKaToPiV0 >= 1 ) {
711711 registry.fill (HIST (" hMassVsPtK0MatchedKaToPi" ), candidateV0.pT , candidateV0.mK0Short );
712712 }
713713 }
714- }
715- } // DecayChannel::DplusV0
714+ }
715+ } // DecayChannel::DplusV0
716716 if (flag != 0 ) {
717717 registry.fill (HIST (" hMCRecCounter" ), flag);
718718 registry.fill (HIST (" hMCRecOrigin" ), origin);
0 commit comments