@@ -267,7 +267,7 @@ struct HfDataCreatorCharmHadPiReduced {
267267
268268 std::array<int , 2 > arrPDGResonantDsPhiPi = {kPhi , kPiPlus }; // Ds± → Phi π±
269269 std::array<int , 2 > arrPDGResonantDKstarK = {kK0Star892 , kKPlus }; // Ds± → K*(892)0bar K± and D± → K*(892)0bar K±
270- std::array<int , 2 > arrPDGResonantDstarD0Pi = {kD0 , kPiPlus }; // D*± → D0 π±
270+ std::array<int , 2 > arrPDGResonantDstarD0Pi = {kD0 , kPiPlus }; // D*± → D0 π±
271271
272272 void init (InitContext& initContext)
273273 {
@@ -388,10 +388,10 @@ struct HfDataCreatorCharmHadPiReduced {
388388
389389 // init HF event selection helper
390390 hfEvSel.init (registry);
391- if (doprocessDplusPiMc || doprocessDplusPiMcWithMl ||
392- doprocessDsPiMc || doprocessDsPiMcWithMl ||
393- doprocessD0PiMc || doprocessD0PiMcWithMl ||
394- doprocessLcPiMc || doprocessLcPiMcWithMl ||
391+ if (doprocessDplusPiMc || doprocessDplusPiMcWithMl ||
392+ doprocessDsPiMc || doprocessDsPiMcWithMl ||
393+ doprocessD0PiMc || doprocessD0PiMcWithMl ||
394+ doprocessLcPiMc || doprocessLcPiMcWithMl ||
395395 doprocessDstarPiMc || doprocessDstarPiMcWithMl) {
396396 const auto & workflows = initContext.services ().get <RunningWorkflowInfo const >();
397397 for (const DeviceSpec& device : workflows.devices ) {
@@ -1213,8 +1213,8 @@ struct HfDataCreatorCharmHadPiReduced {
12131213 trackParCovCharmHad = df2.createParentTrackParCov ();
12141214 trackParCovCharmHad.setAbsCharge (0 ); // to be sure
12151215 } else if constexpr (decChannel == DecayChannel::B0ToDstarPi) {
1216- hCandidatesDstar->Fill (SVFitting::BeforeFit);
1217-
1216+ hCandidatesDstar->Fill (SVFitting::BeforeFit);
1217+
12181218 df2.process (trackParCov0, trackParCov1); // D0 vertex
12191219 df2.propagateTracksToVertex ();
12201220 std::array<float , 3 > pVecPosVtx{}, pVecNegVtx{};
@@ -1241,7 +1241,7 @@ struct HfDataCreatorCharmHadPiReduced {
12411241 df2.getTrack (0 ).getPxPyPzGlo (pVecD0); // D0
12421242 pVecCharm = RecoDecay::pVec (pVecD0, pVec2);
12431243 trackParCovCharmHad = df2.createParentTrackParCov ();
1244- trackParCovCharmHad.setAbsCharge (charmHadDauTracks[2 ].sign ()); // sign of soft pion
1244+ trackParCovCharmHad.setAbsCharge (charmHadDauTracks[2 ].sign ()); // sign of soft pion
12451245 }
12461246
12471247 float ptDauMin = 1 .e6 , etaDauMin = 999 .f , chi2TpcDauMax = -1 .f ;
@@ -1390,18 +1390,18 @@ struct HfDataCreatorCharmHadPiReduced {
13901390 }
13911391 } else if constexpr (decChannel == DecayChannel::BplusToD0barPi) { // D0(bar) → K± π∓
13921392 tables.hfCand2Prong (charmHadDauTracks[0 ].globalIndex (), charmHadDauTracks[1 ].globalIndex (),
1393- indexHfReducedCollision,
1394- trackParCovCharmHad.getX (), trackParCovCharmHad.getAlpha (),
1395- trackParCovCharmHad.getY (), trackParCovCharmHad.getZ (), trackParCovCharmHad.getSnp (),
1396- trackParCovCharmHad.getTgl (), trackParCovCharmHad.getQ2Pt (),
1397- candC.xSecondaryVertex (), candC.ySecondaryVertex (), candC.zSecondaryVertex (), invMassC0, invMassC1,
1398- ptDauMin, etaDauMin, nItsClsDauMin, nTpcCrossRowsDauMin, chi2TpcDauMax);
1393+ indexHfReducedCollision,
1394+ trackParCovCharmHad.getX (), trackParCovCharmHad.getAlpha (),
1395+ trackParCovCharmHad.getY (), trackParCovCharmHad.getZ (), trackParCovCharmHad.getSnp (),
1396+ trackParCovCharmHad.getTgl (), trackParCovCharmHad.getQ2Pt (),
1397+ candC.xSecondaryVertex (), candC.ySecondaryVertex (), candC.zSecondaryVertex (), invMassC0, invMassC1,
1398+ ptDauMin, etaDauMin, nItsClsDauMin, nTpcCrossRowsDauMin, chi2TpcDauMax);
13991399 tables.hfCand2ProngCov (trackParCovCharmHad.getSigmaY2 (), trackParCovCharmHad.getSigmaZY (), trackParCovCharmHad.getSigmaZ2 (),
1400- trackParCovCharmHad.getSigmaSnpY (), trackParCovCharmHad.getSigmaSnpZ (),
1401- trackParCovCharmHad.getSigmaSnp2 (), trackParCovCharmHad.getSigmaTglY (), trackParCovCharmHad.getSigmaTglZ (),
1402- trackParCovCharmHad.getSigmaTglSnp (), trackParCovCharmHad.getSigmaTgl2 (),
1403- trackParCovCharmHad.getSigma1PtY (), trackParCovCharmHad.getSigma1PtZ (), trackParCovCharmHad.getSigma1PtSnp (),
1404- trackParCovCharmHad.getSigma1PtTgl (), trackParCovCharmHad.getSigma1Pt2 ());
1400+ trackParCovCharmHad.getSigmaSnpY (), trackParCovCharmHad.getSigmaSnpZ (),
1401+ trackParCovCharmHad.getSigmaSnp2 (), trackParCovCharmHad.getSigmaTglY (), trackParCovCharmHad.getSigmaTglZ (),
1402+ trackParCovCharmHad.getSigmaTglSnp (), trackParCovCharmHad.getSigmaTgl2 (),
1403+ trackParCovCharmHad.getSigma1PtY (), trackParCovCharmHad.getSigma1PtZ (), trackParCovCharmHad.getSigma1PtSnp (),
1404+ trackParCovCharmHad.getSigma1PtTgl (), trackParCovCharmHad.getSigma1Pt2 ());
14051405 tables.hfCandPidProng0 (candC.nSigTpcPi0 (), candC.nSigTofPi0 (), candC.nSigTpcKa0 (), candC.nSigTofKa0 (), 0 ., 0 ., charmHadDauTracks[0 ].hasTOF (), charmHadDauTracks[0 ].hasTPC ());
14061406 tables.hfCandPidProng1 (candC.nSigTpcPi1 (), candC.nSigTofPi1 (), candC.nSigTpcKa1 (), candC.nSigTofKa1 (), 0 ., 0 ., charmHadDauTracks[1 ].hasTOF (), charmHadDauTracks[1 ].hasTPC ());
14071407 if constexpr (withMl) {
@@ -1597,8 +1597,8 @@ struct HfDataCreatorCharmHadPiReduced {
15971597 counter++;
15981598 }
15991599 tables.rowHfBsMcGenReduced (flag, -1 /* channel*/ , ptParticle, yParticle, etaParticle,
1600- ptProngs[0 ], yProngs[0 ], etaProngs[0 ],
1601- ptProngs[1 ], yProngs[1 ], etaProngs[1 ], hfRejMap, centFT0C, centFT0M);
1600+ ptProngs[0 ], yProngs[0 ], etaProngs[0 ],
1601+ ptProngs[1 ], yProngs[1 ], etaProngs[1 ], hfRejMap, centFT0C, centFT0M);
16021602 } else if constexpr (decayChannel == DecayChannel::BplusToD0barPi) {
16031603 // B+ → D0bar π+
16041604 if (RecoDecay::isMatchedMCGen (particlesMc, particle, Pdg::kBPlus , std::array{-static_cast <int >(Pdg::kD0 ), +kPiPlus }, true )) {
@@ -1630,8 +1630,8 @@ struct HfDataCreatorCharmHadPiReduced {
16301630 counter++;
16311631 }
16321632 tables.rowHfBpMcGenReduced (flag, -1 /* channel*/ , ptParticle, yParticle, etaParticle,
1633- ptProngs[0 ], yProngs[0 ], etaProngs[0 ],
1634- ptProngs[1 ], yProngs[1 ], etaProngs[1 ], hfRejMap, centFT0C, centFT0M);
1633+ ptProngs[0 ], yProngs[0 ], etaProngs[0 ],
1634+ ptProngs[1 ], yProngs[1 ], etaProngs[1 ], hfRejMap, centFT0C, centFT0M);
16351635 } else if constexpr (decayChannel == DecayChannel::LbToLcplusPi) {
16361636 // Lb → Lc+ π-
16371637 if (RecoDecay::isMatchedMCGen<true >(particlesMc, particle, Pdg::kLambdaB0 , std::array{static_cast <int >(Pdg::kLambdaCPlus ), -kPiPlus }, true )) {
0 commit comments