@@ -265,7 +265,8 @@ struct HfDataCreatorCharmHadPiReduced {
265265 void fillMcRecoInfo (const PParticles& particlesMc,
266266 const std::vector<TTrack>& vecDaughtersB,
267267 int & indexHfCandCharm,
268- std::map<int64_t , int64_t > selectedTracksPion) {
268+ std::map<int64_t , int64_t > selectedTracksPion)
269+ {
269270
270271 // we check the MC matching to be stored
271272 int8_t sign{0 };
@@ -369,8 +370,7 @@ struct HfDataCreatorCharmHadPiReduced {
369370 }
370371 rowHfDPiMcCheckReduced (pdgCodeBeautyMother, pdgCodeCharmMother, pdgCodeProng0, pdgCodeProng1, pdgCodeProng2, pdgCodeProng3);
371372 }
372- }
373- else if constexpr (decChannel == DecayChannel::BplusToD0barPi) {
373+ } else if constexpr (decChannel == DecayChannel::BplusToD0barPi) {
374374 // B+ → D0(bar) π+ → (K+ π-) π+
375375 auto indexRec = RecoDecay::getMatchedMCRec (particlesMc, std::array{vecDaughtersB[0 ], vecDaughtersB[1 ], vecDaughtersB[2 ]}, Pdg::kBPlus , std::array{+kPiPlus , +kKPlus , -kPiPlus }, true , &sign, 2 );
376376 if (indexRec > -1 ) {
@@ -436,7 +436,7 @@ struct HfDataCreatorCharmHadPiReduced {
436436 indexHfCandCharm = hfCand3Prong.lastIndex () + 1 ;
437437 invMassC0 = hfHelper.invMassDplusToPiKPi (candC);
438438 registry.fill (HIST (" hMassCharmHad" ), invMassC0);
439- } else if constexpr (decChannel == DecayChannel::BplusToD0barPi){
439+ } else if constexpr (decChannel == DecayChannel::BplusToD0barPi) {
440440 indexHfCandCharm = hfCand2Prong.lastIndex () + 1 ;
441441 if (candC.isSelD0 () >= selectionFlagD0) {
442442 invMassC0 = hfHelper.invMassD0ToPiK (candC);
@@ -501,7 +501,7 @@ struct HfDataCreatorCharmHadPiReduced {
501501 df3.getTrack (1 ).getPxPyPzGlo (pVec1);
502502 df3.getTrack (2 ).getPxPyPzGlo (pVec2);
503503 pVecCharm = RecoDecay::pVec (pVec0, pVec1, pVec2);
504- trackParCovCharmHad.setAbsCharge (charmHadDauTracks[1 ].sign ()); // to be sure
504+ trackParCovCharmHad.setAbsCharge (charmHadDauTracks[1 ].sign ()); // to be sure
505505 } else if constexpr (decChannel == DecayChannel::BplusToD0barPi) { // D0(bar) → K± π∓
506506 if (df2.process (trackParCov0, trackParCov1) == 0 ) {
507507 continue ;
@@ -543,7 +543,7 @@ struct HfDataCreatorCharmHadPiReduced {
543543 if (!isPionSelected (trackPion, trackParCovPion, dcaPion, charmHadDauTracks)) {
544544 continue ;
545545 }
546-
546+
547547 registry.fill (HIST (" hPtPion" ), trackParCovPion.getPt ());
548548 // compute invariant mass square and apply selection
549549 auto invMass2DPi = RecoDecay::m2 (std::array{pVecCharm, pVecPion}, std::array{massC, massPi});
@@ -579,15 +579,15 @@ struct HfDataCreatorCharmHadPiReduced {
579579 fillMcRecoInfo<decChannel>(particlesMc, beautyHadDauTracks, indexHfCandCharm, selectedTracksPion);
580580 }
581581 fillHfCandCharm = true ;
582- } // pion loop
583- if (fillHfCandCharm) { // fill candCplus table only once per D candidate
582+ } // pion loop
583+ if (fillHfCandCharm) { // fill candCplus table only once per D candidate
584584 if constexpr (decChannel == DecayChannel::B0ToDminusPi) { // D∓ → π∓ K± π∓
585585 hfCand3Prong (charmHadDauTracks[0 ].globalIndex (), charmHadDauTracks[1 ].globalIndex (), charmHadDauTracks[2 ].globalIndex (),
586- indexHfReducedCollision,
587- trackParCovCharmHad.getX (), trackParCovCharmHad.getAlpha (),
588- trackParCovCharmHad.getY (), trackParCovCharmHad.getZ (), trackParCovCharmHad.getSnp (),
589- trackParCovCharmHad.getTgl (), trackParCovCharmHad.getQ2Pt (),
590- candC.xSecondaryVertex (), candC.ySecondaryVertex (), candC.zSecondaryVertex (), invMassC0);
586+ indexHfReducedCollision,
587+ trackParCovCharmHad.getX (), trackParCovCharmHad.getAlpha (),
588+ trackParCovCharmHad.getY (), trackParCovCharmHad.getZ (), trackParCovCharmHad.getSnp (),
589+ trackParCovCharmHad.getTgl (), trackParCovCharmHad.getQ2Pt (),
590+ candC.xSecondaryVertex (), candC.ySecondaryVertex (), candC.zSecondaryVertex (), invMassC0);
591591 hfCand3ProngCov (trackParCovCharmHad.getSigmaY2 (), trackParCovCharmHad.getSigmaZY (), trackParCovCharmHad.getSigmaZ2 (),
592592 trackParCovCharmHad.getSigmaSnpY (), trackParCovCharmHad.getSigmaSnpZ (),
593593 trackParCovCharmHad.getSigmaSnp2 (), trackParCovCharmHad.getSigmaTglY (), trackParCovCharmHad.getSigmaTglZ (),
@@ -599,11 +599,11 @@ struct HfDataCreatorCharmHadPiReduced {
599599 }
600600 } else if constexpr (decChannel == DecayChannel::BplusToD0barPi) { // D0(bar) → K± π∓
601601 hfCand2Prong (charmHadDauTracks[0 ].globalIndex (), charmHadDauTracks[1 ].globalIndex (),
602- indexHfReducedCollision,
603- trackParCovCharmHad.getX (), trackParCovCharmHad.getAlpha (),
604- trackParCovCharmHad.getY (), trackParCovCharmHad.getZ (), trackParCovCharmHad.getSnp (),
605- trackParCovCharmHad.getTgl (), trackParCovCharmHad.getQ2Pt (),
606- candC.xSecondaryVertex (), candC.ySecondaryVertex (), candC.zSecondaryVertex (), invMassC0, invMassC1);
602+ indexHfReducedCollision,
603+ trackParCovCharmHad.getX (), trackParCovCharmHad.getAlpha (),
604+ trackParCovCharmHad.getY (), trackParCovCharmHad.getZ (), trackParCovCharmHad.getSnp (),
605+ trackParCovCharmHad.getTgl (), trackParCovCharmHad.getQ2Pt (),
606+ candC.xSecondaryVertex (), candC.ySecondaryVertex (), candC.zSecondaryVertex (), invMassC0, invMassC1);
607607 hfCand2ProngCov (trackParCovCharmHad.getSigmaY2 (), trackParCovCharmHad.getSigmaZY (), trackParCovCharmHad.getSigmaZ2 (),
608608 trackParCovCharmHad.getSigmaSnpY (), trackParCovCharmHad.getSigmaSnpZ (),
609609 trackParCovCharmHad.getSigmaSnp2 (), trackParCovCharmHad.getSigmaTglY (), trackParCovCharmHad.getSigmaTglZ (),
@@ -613,10 +613,10 @@ struct HfDataCreatorCharmHadPiReduced {
613613 if constexpr (withMl) {
614614 std::array<float , 6 > bdtScores = {-1 .f , -1 .f , -1 .f , -1 .f , -1 .f , -1 .f };
615615 if (candC.mlProbD0 ().size () == 3 ) {
616- std::copy (candC.mlProbD0 ().begin (), candC.mlProbD0 ().end (), bdtScores.begin ());
616+ std::copy (candC.mlProbD0 ().begin (), candC.mlProbD0 ().end (), bdtScores.begin ());
617617 }
618618 if (candC.mlProbD0bar ().size () == 3 ) {
619- std::copy (candC.mlProbD0bar ().begin (), candC.mlProbD0bar ().end (), bdtScores.begin () + 3 );
619+ std::copy (candC.mlProbD0bar ().begin (), candC.mlProbD0bar ().end (), bdtScores.begin () + 3 );
620620 }
621621
622622 hfCand2ProngMl (bdtScores[0 ], bdtScores[1 ], bdtScores[2 ], bdtScores[3 ], bdtScores[4 ], bdtScores[5 ]);
@@ -639,7 +639,7 @@ struct HfDataCreatorCharmHadPiReduced {
639639 bz);
640640 }
641641
642- template <uint8_t decayChannel>
642+ template <uint8_t decayChannel>
643643 void runMcGen (aod::McParticles const & particlesMc)
644644 {
645645 // Match generated particles.
@@ -680,7 +680,7 @@ struct HfDataCreatorCharmHadPiReduced {
680680 ptProngs[0 ], yProngs[0 ], etaProngs[0 ],
681681 ptProngs[1 ], yProngs[1 ], etaProngs[1 ]);
682682 } else if constexpr (decayChannel == DecayChannel::BplusToD0barPi) {
683- // B+ → D0bar π+
683+ // B+ → D0bar π+
684684 if (RecoDecay::isMatchedMCGen (particlesMc, particle, Pdg::kBPlus , std::array{static_cast <int >(Pdg::kD0 ), +kPiPlus }, true )) {
685685 // Match D0bar -> π- K+
686686 auto candD0MC = particlesMc.rawIteratorAt (particle.daughtersIds ().front ());
0 commit comments