@@ -159,7 +159,7 @@ struct HfDataCreatorJpsiHadReduced {
159159 Configurable<LabeledArray<double >> cuts{" cuts" , {hf_cuts_jpsi_to_mu_mu::Cuts[0 ], hf_cuts_jpsi_to_mu_mu::NBinsPt, hf_cuts_jpsi_to_mu_mu::NCutVars, hf_cuts_jpsi_to_mu_mu::labelsPt, hf_cuts_jpsi_to_mu_mu::labelsCutVar}, " J/Psi candidate selection per pT bin" };
160160 Configurable<double > invMassWindowJpsiHad{" invMassWindowJpsiHad" , 0.3 , " invariant-mass window for Jpsi-Had pair preselections (GeV/c2)" };
161161 Configurable<double > deltaMPhiMax{" deltaMPhiMax" , 0.02 , " invariant-mass window for phi preselections (GeV/c2) (only for Bs->J/PsiPhi)" };
162- Configurable<double > deltaMK0StarMax{" deltaMK0StarMax" , 0.02 , " invariant-mass window for K*0 preselections (GeV/c2) (only for B0->J/PsiK0*)" };
162+ Configurable<double > deltaMK0StarMax{" deltaMK0StarMax" , 0.15 , " invariant-mass window for K*0 preselections (GeV/c2) (only for B0->J/PsiK0*)" };
163163 Configurable<double > cpaMin{" cpaMin" , 0 ., " Minimum cosine of pointing angle for B candidates" };
164164 Configurable<double > decLenMin{" decLenMin" , 0 ., " Minimum decay length for B candidates" };
165165
@@ -647,7 +647,7 @@ struct HfDataCreatorJpsiHadReduced {
647647 }
648648 // Check K*0 -> K+pi-
649649 if (RecoDecay::isMatchedMCGen (particlesMc, candK0StarMC, static_cast <int >(Pdg::kK0Star892 ), std::array{+kKPlus , -kPiPlus }, true )) {
650- channel = o2::hf_decay::hf_cand_beauty::B0ToJpsiPiK ;
650+ channel = o2::hf_decay::hf_cand_beauty::B0ToJpsiKstar0 ;
651651 }
652652 }
653653
@@ -927,6 +927,9 @@ struct HfDataCreatorJpsiHadReduced {
927927 } else if constexpr (DecChannel == DecayChannel::B0ToJpsiK0Star) {
928928 for (auto trackBachId2 = trackId + 1 ; trackBachId2 != trackIndices.end (); ++trackBachId2) {
929929 auto trackBach2 = trackBachId2.template track_as <TTracks>();
930+ if (trackBach.sign () == trackBach2.sign ()) {
931+ continue ;
932+ }
930933 auto trackBach2ParCov = getTrackParCov (trackBach2);
931934
932935 std::array<float , 2 > dcaBach2{trackBach2.dcaXY (), trackBach2.dcaZ ()};
@@ -1054,6 +1057,9 @@ struct HfDataCreatorJpsiHadReduced {
10541057 } else if constexpr (DecChannel == DecayChannel::BsToJpsiPhi) {
10551058 for (auto trackBachId2 = trackId + 1 ; trackBachId2 != trackIndices.end (); ++trackBachId2) {
10561059 auto trackBach2 = trackBachId2.template track_as <TTracks>();
1060+ if (trackBach.sign () == trackBach2.sign ()) {
1061+ continue ;
1062+ }
10571063 auto trackBach2ParCov = getTrackParCov (trackBach2);
10581064
10591065 std::array<float , 2 > dcaBach2{trackBach2.dcaXY (), trackBach2.dcaZ ()};
0 commit comments