From 9e1514bbc0b4d71d523e0d78289a350528b787e9 Mon Sep 17 00:00:00 2001 From: Marta Razza Date: Fri, 10 Jul 2026 18:49:00 +0200 Subject: [PATCH 1/3] PWGHF/D2H: PID selections changed in taskDeuteronFromLb --- PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx | 18 ++++++++---------- 1 file changed, 8 insertions(+), 10 deletions(-) diff --git a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx index db0af2da598..2c050c1e9c2 100644 --- a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx +++ b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx @@ -66,10 +66,9 @@ struct HfTaskDeuteronFromLb { Configurable cfgITScls{"cfgITScls", 2, "Minimum ITS clusters"}; Configurable cfgMaxPt{"cfgMaxPt", 5.0f, "Maximum pT cut"}; Configurable cfgMinPt{"cfgMinPt", 0.5f, "Minimum pT cut"}; - Configurable cfgTPCNsigma{"cfgTPCNsigma", 4.0f, "TPC n sigma for deuteron PID"}; - Configurable cfgTofNsigmaMin{"cfgTofNsigmaMin", 3.0f, "TOF n sigma min for deuteron PID"}; - Configurable cfgTofNsigmaMax{"cfgTofNsigmaMax", 4.0f, "TOF n sigma max for deuteron PID"}; - Configurable ptThresholdPid{"ptThresholdPid", 0.5f, "pT threshold to switch between 4 and 3 sigmas for TOF PID"}; + Configurable cfgTPCNsigma{"cfgTPCNsigma", 3.0f, "TPC n sigma for deuteron PID"}; + Configurable cfgTofNsigma{"cfgTofNsigma", 3.0f, "TOF n sigma for deuteron PID"}; + Configurable ptThresholdPid{"ptThresholdPid", 0.5f, "pT threshold to switch between TPC and TPC+TOF PID"}; Configurable cfgDCAmin{"cfgDCAmin", 0.05f, "Minimum DCA for deuteron PID"}; Configurable cfgDCAmax{"cfgDCAmax", 1000.0f, "Maximum DCA for deuteron PID"}; Configurable rapidityCut{"rapidityCut", 0.5f, "Rapidity cut"}; @@ -260,11 +259,10 @@ struct HfTaskDeuteronFromLb { } const bool isTPCDe = std::abs(track.tpcNSigmaDe()) < cfgTPCNsigma; - const bool isTOFDe_min = std::abs(track.tofNSigmaDe()) > cfgTofNsigmaMin; - const bool isTOFDe_max = std::abs(track.tofNSigmaDe()) < cfgTofNsigmaMax; + const bool isTOFDe = std::abs(track.tofNSigmaDe()) < cfgTofNsigma; if (track.pt() < ptThresholdPid) { - if (isTPCDe && isTOFDe_max) { + if (isTPCDe) { QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); @@ -273,7 +271,7 @@ struct HfTaskDeuteronFromLb { QAHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe()); } } else { - if (isTPCDe && isTOFDe_min && isTOFDe_max) { + if (isTPCDe && isTOFDe) { QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); @@ -285,7 +283,7 @@ struct HfTaskDeuteronFromLb { } } } - PROCESS_SWITCH(HfTaskDeuteronFromLb, processData, "processData", false); + PROCESS_SWITCH(HfTaskDeuteronFromLb, processData, "processData", true); void processMC(MCCollisionCandidates::iterator const& collision, MCTrackCandidates const& tracks, o2::aod::McParticles const&) { @@ -423,7 +421,7 @@ struct HfTaskDeuteronFromLb { } } - PROCESS_SWITCH(HfTaskDeuteronFromLb, processGen, "processGen", true); + PROCESS_SWITCH(HfTaskDeuteronFromLb, processGen, "processGen", false); }; WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) From 85fe02c283561ce2926722214bcb6fd5a4570756 Mon Sep 17 00:00:00 2001 From: Marta Razza Date: Thu, 16 Jul 2026 20:05:38 +0200 Subject: [PATCH 2/3] Fix DCAxy selection and Errors/Warnings in taskDeuteronFromLb --- PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx | 46 +++++++++++++++++--------- 1 file changed, 31 insertions(+), 15 deletions(-) diff --git a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx index 2c050c1e9c2..02c0f9638ce 100644 --- a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx +++ b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx @@ -81,7 +81,7 @@ struct HfTaskDeuteronFromLb { float d_bz = 0.f; int mCurrentRun = -1; - framework::Service ccdb; + framework::Service ccdb{}; using CollisionCandidates = o2::soa::Join; using MCTrackCandidates = o2::soa::Join; @@ -184,30 +184,39 @@ struct HfTaskDeuteronFromLb { template bool passedSingleTrackSelection(const T1& track) { - if (std::abs(track.eta()) > cfgEta) + if (std::abs(track.eta()) > cfgEta) { return false; - if (std::abs(track.dcaXY()) < cfgDCAmin || std::abs(track.dcaXY()) > cfgDCAmax) - return false; - if (!track.hasITS()) + } + if (!track.hasITS()) { return false; - if (!track.hasTPC()) + } + if (!track.hasTPC()) { return false; - if (!track.hasTOF()) + } + if (!track.hasTOF()) { return false; - if (track.tpcNClsFound() < cfgTPCNclsFound) + } + if (track.tpcNClsFound() < cfgTPCNclsFound) { return false; - if (track.tpcChi2NCl() > cfgTPCChi2Ncl) + } + if (track.tpcChi2NCl() > cfgTPCChi2Ncl) { return false; - if (track.itsChi2NCl() > cfgITSChi2Ncl) + } + if (track.itsChi2NCl() > cfgITSChi2Ncl) { return false; - if (track.itsNCls() < cfgITScls) + } + if (track.itsNCls() < cfgITScls) { return false; - if (track.pt() > cfgMaxPt) + } + if (track.pt() > cfgMaxPt) { return false; - if (track.pt() < cfgMinPt) + } + if (track.pt() < cfgMinPt) { return false; - if (track.sign() > 0) + } + if (track.sign() > 0) { return false; + } return true; } @@ -219,7 +228,7 @@ struct HfTaskDeuteronFromLb { { for (const auto& collision : collisions) { if (mCurrentRun != collision.bc_as().runNumber()) { - o2::parameters::GRPMagField* grpo = ccdb->getForTimeStamp("GLO/Config/GRPMagField", collision.bc_as().timestamp()); + auto* grpo = ccdb->getForTimeStamp("GLO/Config/GRPMagField", collision.bc_as().timestamp()); o2::base::Propagator::initFieldFromGRP(grpo); mCurrentRun = collision.bc_as().runNumber(); } @@ -254,6 +263,10 @@ struct HfTaskDeuteronFromLb { o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackPar, 2.f, noMatCorr, &dca); } + if (std::abs(dca[0]) < cfgDCAmin || std::abs(dca[0]) > cfgDCAmax) { + continue; + } + if (!passedSingleTrackSelection(track)) { continue; } @@ -288,6 +301,9 @@ struct HfTaskDeuteronFromLb { void processMC(MCCollisionCandidates::iterator const& collision, MCTrackCandidates const& tracks, o2::aod::McParticles const&) { for (const auto& track : tracks) { + if (std::abs(track.dcaXY()) < cfgDCAmin || std::abs(track.dcaXY()) > cfgDCAmax) { + continue; + } if (!passedSingleTrackSelection(track)) { continue; } From 482a2ba3b114f5f7bcb63c298b5179d5cd5e3b36 Mon Sep 17 00:00:00 2001 From: Marta Razza Date: Thu, 16 Jul 2026 20:25:09 +0200 Subject: [PATCH 3/3] Fix for O2 linter warinings --- PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx | 112 ++++++++++++------------- 1 file changed, 56 insertions(+), 56 deletions(-) diff --git a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx index 02c0f9638ce..09ddb9cba1f 100644 --- a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx +++ b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx @@ -78,7 +78,7 @@ struct HfTaskDeuteronFromLb { Configurable pdgCodeDaughter{"pdgCodeDaughter", -1000010020, "PDG code of the daughter particle (default: anti-deuteron)"}; int mRunNumber = 0; - float d_bz = 0.f; + float dBz = 0.f; int mCurrentRun = -1; framework::Service ccdb{}; @@ -95,7 +95,7 @@ struct HfTaskDeuteronFromLb { ConfigurableAxis dcaXyAxis{"dcaXyAxis", {1000, -0.2f, 0.2f}, "DCA xy (cm)"}; ConfigurableAxis dcaZAxis{"dcaZAxis", {1000, -0.2f, 0.2f}, "DCA z (cm)"}; - HistogramRegistry QAHistos{"QAHistos", {}, OutputObjHandlingPolicy::AnalysisObject, false, true}; + HistogramRegistry qaHistos{"QAHistos", {}, OutputObjHandlingPolicy::AnalysisObject, false, true}; HistogramRegistry registry{"registry", {}, OutputObjHandlingPolicy::AnalysisObject}; OutputObj zorroSummary{"zorroSummary"}; @@ -129,38 +129,38 @@ struct HfTaskDeuteronFromLb { zorroSummary.setObject(zorro.getZorroSummary()); } - QAHistos.add("Data/hDCAxyVsPt", "DCAxy #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaXyAxis}}); - QAHistos.add("Data/hDCAzVsPt", "DCAz #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaZAxis}}); - QAHistos.add("Data/hnSigmaTPCVsPt", "n#sigma TPC vs p_{T} for #bar{d} hypothesis; p_{T} (GeV/c); n#sigma TPC", {HistType::kTH2D, {ptAxis, nSigmaAxis}}); - QAHistos.add("Data/hnSigmaTOFVsPt", "n#sigma TOF vs p_{T} for #bar{d} hypothesis; p_{T} (GeV/c); n#sigma TOF", {HistType::kTH2D, {ptAxis, nSigmaAxis}}); - QAHistos.add("Data/ptAntiDeuteron", "ptAntiDeuteron", {HistType::kTH1F, {ptAxis}}); - QAHistos.add("Data/etaAntideuteron", "etaAntideuteron", {HistType::kTH1F, {{100, -1.0f, 1.0f, "eta #bar{d}"}}}); - QAHistos.add("Data/hVtxZ", "Z-Vertex distribution after selection;Z (cm)", HistType::kTH1F, {{100, -50, 50}}); + qaHistos.add("Data/hDCAxyVsPt", "DCAxy #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaXyAxis}}); + qaHistos.add("Data/hDCAzVsPt", "DCAz #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaZAxis}}); + qaHistos.add("Data/hnSigmaTPCVsPt", "n#sigma TPC vs p_{T} for #bar{d} hypothesis; p_{T} (GeV/c); n#sigma TPC", {HistType::kTH2D, {ptAxis, nSigmaAxis}}); + qaHistos.add("Data/hnSigmaTOFVsPt", "n#sigma TOF vs p_{T} for #bar{d} hypothesis; p_{T} (GeV/c); n#sigma TOF", {HistType::kTH2D, {ptAxis, nSigmaAxis}}); + qaHistos.add("Data/ptAntiDeuteron", "ptAntiDeuteron", {HistType::kTH1F, {ptAxis}}); + qaHistos.add("Data/etaAntideuteron", "etaAntideuteron", {HistType::kTH1F, {{100, -1.0f, 1.0f, "eta #bar{d}"}}}); + qaHistos.add("Data/hVtxZ", "Z-Vertex distribution after selection;Z (cm)", HistType::kTH1F, {{100, -50, 50}}); // MC generated-level histograms - QAHistos.add("MCGen/ptGeneratedBminus", "p_{T} generated B^{-};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCGen/ptGeneratedAntiLambdaB", "p_{T} generated #bar{#Lambda}_{b};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); + qaHistos.add("MCGen/ptGeneratedBminus", "p_{T} generated B^{-};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); + qaHistos.add("MCGen/ptGeneratedAntiLambdaB", "p_{T} generated #bar{#Lambda}_{b};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCGen/ptAntiDeuteronPrimary", "p_{T} #bar{d} primary gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCGen/ptAntiDeuteronFromBeautyHadron", "p_{T} #bar{d} from beauty gen hadrons;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); + qaHistos.add("MCGen/ptAntiDeuteronPrimary", "p_{T} #bar{d} primary gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); + qaHistos.add("MCGen/ptAntiDeuteronFromBeautyHadron", "p_{T} #bar{d} from beauty gen hadrons;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCGen/hMotherPdgCode", "PDG code of mother, gen level;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}}); - QAHistos.add("MCGen/ctauBminus", "ctau of B^{-}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); - QAHistos.add("MCGen/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); + qaHistos.add("MCGen/hMotherPdgCode", "PDG code of mother, gen level;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}}); + qaHistos.add("MCGen/ctauBminus", "ctau of B^{-}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); + qaHistos.add("MCGen/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); // MC reco/MC-anchored histograms - QAHistos.add("MCReco/ptAntiDeuteronFromBminus", + qaHistos.add("MCReco/ptAntiDeuteronFromBminus", "p_{T} #bar{d} from B^{-} reco/MC anchored;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCReco/ptAntiDeuteronFromAntiLambdaB", + qaHistos.add("MCReco/ptAntiDeuteronFromAntiLambdaB", "p_{T} #bar{d} from #bar{#Lambda}_{b} reco/MC anchored;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCGen/ptAntiDeuteronFromBminus", "p_{T} #bar{d} from B^{-} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCGen/ptAntiDeuteronFromAntiLambdaB", "p_{T} #bar{d} from #bar{#Lambda}_{b} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); - QAHistos.add("MCReco/hDCAxy-Primary", "DCAxy primary reco/MC anchored;DCA xy (cm);Counts", {HistType::kTH1D, {{400, -0.2f, 0.2f, "DCA xy (cm)"}}}); - QAHistos.add("MCReco/hDCAxy-FromBeautyHadron", "DCAxy from beauty reco/MC anchored;DCA xy (cm);Counts", {HistType::kTH1D, {{400, -0.2f, 0.2f, "DCA xy (cm)"}}}); - QAHistos.add("MCReco/hMotherPdgCode", "PDG code of mother, reco/MC anchored;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}}); - QAHistos.add("MCReco/ctauBminus", "ctau of B^{-}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); - QAHistos.add("MCReco/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); + qaHistos.add("MCGen/ptAntiDeuteronFromBminus", "p_{T} #bar{d} from B^{-} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); + qaHistos.add("MCGen/ptAntiDeuteronFromAntiLambdaB", "p_{T} #bar{d} from #bar{#Lambda}_{b} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis}); + qaHistos.add("MCReco/hDCAxy-Primary", "DCAxy primary reco/MC anchored;DCA xy (cm);Counts", {HistType::kTH1D, {{400, -0.2f, 0.2f, "DCA xy (cm)"}}}); + qaHistos.add("MCReco/hDCAxy-FromBeautyHadron", "DCAxy from beauty reco/MC anchored;DCA xy (cm);Counts", {HistType::kTH1D, {{400, -0.2f, 0.2f, "DCA xy (cm)"}}}); + qaHistos.add("MCReco/hMotherPdgCode", "PDG code of mother, reco/MC anchored;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}}); + qaHistos.add("MCReco/ctauBminus", "ctau of B^{-}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); + qaHistos.add("MCReco/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}}); hProcessedEvents->GetXaxis()->SetBinLabel(1, "Events processed"); hProcessedEvents->GetXaxis()->SetBinLabel(2, "ZORRO"); @@ -173,7 +173,7 @@ struct HfTaskDeuteronFromLb { if (mRunNumber == bc.runNumber()) { return; } - d_bz = 0.f; + dBz = 0.f; if (applySkimming) { zorro.initCCDB(ccdb.service, bc.runNumber(), bc.timestamp(), cfgSkimming.value); @@ -250,7 +250,7 @@ struct HfTaskDeuteronFromLb { continue; } hProcessedEvents->Fill(3.5); - QAHistos.fill(HIST("Data/hVtxZ"), collision.posZ()); + qaHistos.fill(HIST("Data/hVtxZ"), collision.posZ()); const auto& trackIdsThisCollision = trackIndices.sliceBy(trackIndicesPerCollision, collision.globalIndex()); @@ -276,21 +276,21 @@ struct HfTaskDeuteronFromLb { if (track.pt() < ptThresholdPid) { if (isTPCDe) { - QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); - QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); - QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); - QAHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]); - QAHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe()); - QAHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe()); + qaHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); + qaHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); + qaHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); + qaHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]); + qaHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe()); + qaHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe()); } } else { if (isTPCDe && isTOFDe) { - QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); - QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); - QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); - QAHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]); - QAHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe()); - QAHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe()); + qaHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); + qaHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); + qaHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); + qaHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]); + qaHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe()); + qaHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe()); } } } @@ -321,7 +321,7 @@ struct HfTaskDeuteronFromLb { int motherPdg = 0; if (separateAntideuterons) { for (const auto& mom : mcParticle.mothers_as()) { - QAHistos.fill(HIST("MCReco/hMotherPdgCode"), mom.pdgCode()); + qaHistos.fill(HIST("MCReco/hMotherPdgCode"), mom.pdgCode()); if (mom.pdgCode() == pdgCodeBeautyMeson || mom.pdgCode() == pdgCodeBeautyBaryon) { motherPdg = mom.pdgCode(); double dx = mcParticle.vx() - collision.posX(); @@ -336,9 +336,9 @@ struct HfTaskDeuteronFromLb { ctauBeauty *= 1.e4; if (motherPdg == pdgCodeBeautyMeson) { - QAHistos.fill(HIST("MCReco/ctauBminus"), ctauBeauty); + qaHistos.fill(HIST("MCReco/ctauBminus"), ctauBeauty); } else if (motherPdg == pdgCodeBeautyBaryon) { - QAHistos.fill(HIST("MCReco/ctauAntiLambdaB"), ctauBeauty); + qaHistos.fill(HIST("MCReco/ctauAntiLambdaB"), ctauBeauty); } } break; @@ -346,14 +346,14 @@ struct HfTaskDeuteronFromLb { } } if (motherPdg == pdgCodeBeautyMeson) { - QAHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY()); - QAHistos.fill(HIST("MCReco/ptAntiDeuteronFromBminus"), track.pt()); + qaHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY()); + qaHistos.fill(HIST("MCReco/ptAntiDeuteronFromBminus"), track.pt()); } else if (motherPdg == pdgCodeBeautyBaryon) { - QAHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY()); - QAHistos.fill(HIST("MCReco/ptAntiDeuteronFromAntiLambdaB"), track.pt()); + qaHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY()); + qaHistos.fill(HIST("MCReco/ptAntiDeuteronFromAntiLambdaB"), track.pt()); } else { - QAHistos.fill(HIST("MCReco/hDCAxy-Primary"), track.dcaXY()); - QAHistos.fill(HIST("MCReco/ptAntiDeuteronPrimary"), track.pt()); + qaHistos.fill(HIST("MCReco/hDCAxy-Primary"), track.dcaXY()); + qaHistos.fill(HIST("MCReco/ptAntiDeuteronPrimary"), track.pt()); } } } @@ -371,9 +371,9 @@ struct HfTaskDeuteronFromLb { if (isSelectedBeautyHadron(mcParticle.pdgCode())) { if (std::abs(mcParticle.y()) <= rapidityCut) { if (mcParticle.pdgCode() == pdgCodeBeautyMeson) { - QAHistos.fill(HIST("MCGen/ptGeneratedBminus"), mcParticle.pt()); + qaHistos.fill(HIST("MCGen/ptGeneratedBminus"), mcParticle.pt()); } else if (mcParticle.pdgCode() == pdgCodeBeautyBaryon) { - QAHistos.fill(HIST("MCGen/ptGeneratedAntiLambdaB"), mcParticle.pt()); + qaHistos.fill(HIST("MCGen/ptGeneratedAntiLambdaB"), mcParticle.pt()); } } @@ -392,9 +392,9 @@ struct HfTaskDeuteronFromLb { ctauBeauty *= 1.e4; // cm -> micrometers if (mcParticle.pdgCode() == pdgCodeBeautyMeson) { - QAHistos.fill(HIST("MCGen/ctauBminus"), ctauBeauty); + qaHistos.fill(HIST("MCGen/ctauBminus"), ctauBeauty); } else if (mcParticle.pdgCode() == pdgCodeBeautyBaryon) { - QAHistos.fill(HIST("MCGen/ctauAntiLambdaB"), ctauBeauty); + qaHistos.fill(HIST("MCGen/ctauAntiLambdaB"), ctauBeauty); } } break; @@ -412,7 +412,7 @@ struct HfTaskDeuteronFromLb { if (mcParticle.has_mothers()) { for (const auto& mom : mcParticle.mothers_as()) { - QAHistos.fill(HIST("MCGen/hMotherPdgCode"), mom.pdgCode()); + qaHistos.fill(HIST("MCGen/hMotherPdgCode"), mom.pdgCode()); if (mom.pdgCode() == pdgCodeBeautyMeson) { motherPdg = mom.pdgCode(); @@ -427,11 +427,11 @@ struct HfTaskDeuteronFromLb { } if (motherPdg == pdgCodeBeautyMeson) { - QAHistos.fill(HIST("MCGen/ptAntiDeuteronFromBminus"), mcParticle.pt()); + qaHistos.fill(HIST("MCGen/ptAntiDeuteronFromBminus"), mcParticle.pt()); } else if (motherPdg == pdgCodeBeautyBaryon) { - QAHistos.fill(HIST("MCGen/ptAntiDeuteronFromAntiLambdaB"), mcParticle.pt()); + qaHistos.fill(HIST("MCGen/ptAntiDeuteronFromAntiLambdaB"), mcParticle.pt()); } else if (mcParticle.isPhysicalPrimary()) { - QAHistos.fill(HIST("MCGen/ptAntiDeuteronPrimary"), mcParticle.pt()); + qaHistos.fill(HIST("MCGen/ptAntiDeuteronPrimary"), mcParticle.pt()); } } }