@@ -55,10 +55,10 @@ struct HfTaskBs {
5555 void init (o2::framework::InitContext&)
5656 {
5757 auto vbins = (std::vector<double >)bins;
58- registry.add (" hMass" , " B_{s} candidates;inv. mass #Lambda_{c }^{#plus}#pi^{#minus} (GeV/#it{c}^{2});#it{p}_{T} (GeV/#it{c}); centrality" , {HistType::kTH3F , {{500 , 0 ., 10 .}, {vbins, " #it{p}_{T} (GeV/#it{c})" }, {100 , 0 ., 100 .}}});
58+ registry.add (" hMass" , " B_{s} candidates;inv. mass D_{s }^{#plus}#pi^{#minus} (GeV/#it{c}^{2});#it{p}_{T} (GeV/#it{c}); centrality" , {HistType::kTH3F , {{500 , 0 ., 10 .}, {vbins, " #it{p}_{T} (GeV/#it{c})" }, {100 , 0 ., 100 .}}});
5959 registry.add (" hDecLength" , " B_{s} candidates;decay length (cm);entries" , {HistType::kTH2F , {{200 , 0 ., 0.4 }, {vbins, " #it{p}_{T} (GeV/#it{c})" }}});
6060 registry.add (" hDecLengthXY" , " B_{s} candidates;decay length xy (cm);entries" , {HistType::kTH2F , {{200 , 0 ., 0.4 }, {vbins, " #it{p}_{T} (GeV/#it{c})" }}});
61- registry.add (" hd0Prong0" , " B_{s} candidates;prong 0 (#Lambda_{c }^{#plus}) DCAxy to prim. vertex (cm);entries" , {HistType::kTH2F , {{100 , -0.05 , 0.05 }, {vbins, " #it{p}_{T} (GeV/#it{c})" }}});
61+ registry.add (" hd0Prong0" , " B_{s} candidates;prong 0 (D_{s }^{#plus}) DCAxy to prim. vertex (cm);entries" , {HistType::kTH2F , {{100 , -0.05 , 0.05 }, {vbins, " #it{p}_{T} (GeV/#it{c})" }}});
6262 registry.add (" hd0Prong1" , " B_{s} candidates;prong 1 (#pi^{#minus}) DCAxy to prim. vertex (cm);entries" , {HistType::kTH2F , {{100 , -0.05 , 0.05 }, {vbins, " #it{p}_{T} (GeV/#it{c})" }}});
6363 registry.add (" hCPA" , " B_{s} candidates;B_{s} candidate cosine of pointing angle;entries" , {HistType::kTH2F , {{110 , -1.1 , 1.1 }, {vbins, " #it{p}_{T} (GeV/#it{c})" }}});
6464 registry.add (" hEta" , " B_{s} candidates;B_{s} candidate #it{#eta};entries" , {HistType::kTH2F , {{100 , -2 ., 2 .}, {vbins, " #it{p}_{T} (GeV/#it{c})" }}});
@@ -190,7 +190,7 @@ struct HfTaskBsMc {
190190 auto candDs = candidate.index0_as <aod::HfCandProng3>();
191191 if (std::abs (candidate.flagMCMatchRec ()) == 1 << hf_cand_bs::DecayType::BsToDsPi) {
192192
193- auto indexMother = RecoDecay::getMother (particlesMC, candidate.index1_as <aod::BigTracksMC>().mcParticle_as <soa::Join<aod::McParticles, aod::HfCandBsMCGen>>(), pdg::Code::kLambdaB0 , true );
193+ auto indexMother = RecoDecay::getMother (particlesMC, candidate.index1_as <aod::BigTracksMC>().mcParticle_as <soa::Join<aod::McParticles, aod::HfCandBsMCGen>>(), pdg::Code::kBs , true );
194194 auto particleMother = particlesMC.iteratorAt (indexMother);
195195 registry.fill (HIST (" hPtGenSig" ), particleMother.pt ());
196196 registry.fill (HIST (" hPtRecSig" ), candidate.pt ());
@@ -238,7 +238,7 @@ struct HfTaskBsMc {
238238 for (auto & particle : particlesMC) {
239239 if (std::abs (particle.flagMCMatchGen ()) == 1 << hf_cand_bs::DecayType::BsToDsPi) {
240240
241- auto yParticle = RecoDecay::Y (array{particle.px (), particle.py (), particle.pz ()}, RecoDecay::getMassPDG (pdg::Code::kLambdaB0 ));
241+ auto yParticle = RecoDecay::Y (array{particle.px (), particle.py (), particle.pz ()}, RecoDecay::getMassPDG (pdg::Code::kBs ));
242242 if (cutYCandMax >= 0 . && std::abs (yParticle) > cutYCandMax) {
243243 continue ;
244244 }
0 commit comments