Skip to content

Commit 617939f

Browse files
MRazza879Marta Razza
MRazza879
and
Marta Razza
authored
[PWGHF] Fix DCAxy selection and code-check errors/warnings in taskDeuteronFromLb (#17098)
Co-authored-by: Marta Razza <marta.razza@cern.ch>
1 parent cd60106 commit 617939f

1 file changed

Lines changed: 87 additions & 71 deletions

File tree

‎PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx‎

Lines changed: 87 additions & 71 deletions
Original file line numberDiff line numberDiff line change
@@ -78,10 +78,10 @@ struct HfTaskDeuteronFromLb {
7878
Configurable<int> pdgCodeDaughter{"pdgCodeDaughter", -1000010020, "PDG code of the daughter particle (default: anti-deuteron)"};
7979

8080
int mRunNumber = 0;
81-
float d_bz = 0.f;
81+
float dBz = 0.f;
8282
int mCurrentRun = -1;
8383

84-
framework::Service<ccdb::BasicCCDBManager> ccdb;
84+
framework::Service<ccdb::BasicCCDBManager> ccdb{};
8585

8686
using CollisionCandidates = o2::soa::Join<o2::aod::Collisions, o2::aod::EvSels>;
8787
using MCTrackCandidates = o2::soa::Join<o2::aod::TracksIU, o2::aod::TracksExtra, o2::aod::TracksDCA, o2::aod::McTrackLabels>;
@@ -95,7 +95,7 @@ struct HfTaskDeuteronFromLb {
9595
ConfigurableAxis dcaXyAxis{"dcaXyAxis", {1000, -0.2f, 0.2f}, "DCA xy (cm)"};
9696
ConfigurableAxis dcaZAxis{"dcaZAxis", {1000, -0.2f, 0.2f}, "DCA z (cm)"};
9797

98-
HistogramRegistry QAHistos{"QAHistos", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};
98+
HistogramRegistry qaHistos{"QAHistos", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};
9999
HistogramRegistry registry{"registry", {}, OutputObjHandlingPolicy::AnalysisObject};
100100

101101
OutputObj<ZorroSummary> zorroSummary{"zorroSummary"};
@@ -129,38 +129,38 @@ struct HfTaskDeuteronFromLb {
129129
zorroSummary.setObject(zorro.getZorroSummary());
130130
}
131131

132-
QAHistos.add("Data/hDCAxyVsPt", "DCAxy #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaXyAxis}});
133-
QAHistos.add("Data/hDCAzVsPt", "DCAz #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaZAxis}});
134-
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}});
135-
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}});
136-
QAHistos.add("Data/ptAntiDeuteron", "ptAntiDeuteron", {HistType::kTH1F, {ptAxis}});
137-
QAHistos.add("Data/etaAntideuteron", "etaAntideuteron", {HistType::kTH1F, {{100, -1.0f, 1.0f, "eta #bar{d}"}}});
138-
QAHistos.add("Data/hVtxZ", "Z-Vertex distribution after selection;Z (cm)", HistType::kTH1F, {{100, -50, 50}});
132+
qaHistos.add("Data/hDCAxyVsPt", "DCAxy #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaXyAxis}});
133+
qaHistos.add("Data/hDCAzVsPt", "DCAz #bar{d} vs p_{T}", {HistType::kTH2D, {ptAxis, dcaZAxis}});
134+
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}});
135+
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}});
136+
qaHistos.add("Data/ptAntiDeuteron", "ptAntiDeuteron", {HistType::kTH1F, {ptAxis}});
137+
qaHistos.add("Data/etaAntideuteron", "etaAntideuteron", {HistType::kTH1F, {{100, -1.0f, 1.0f, "eta #bar{d}"}}});
138+
qaHistos.add("Data/hVtxZ", "Z-Vertex distribution after selection;Z (cm)", HistType::kTH1F, {{100, -50, 50}});
139139
// MC generated-level histograms
140-
QAHistos.add("MCGen/ptGeneratedBminus", "p_{T} generated B^{-};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
141-
QAHistos.add("MCGen/ptGeneratedAntiLambdaB", "p_{T} generated #bar{#Lambda}_{b};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
140+
qaHistos.add("MCGen/ptGeneratedBminus", "p_{T} generated B^{-};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
141+
qaHistos.add("MCGen/ptGeneratedAntiLambdaB", "p_{T} generated #bar{#Lambda}_{b};p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
142142

143-
QAHistos.add("MCGen/ptAntiDeuteronPrimary", "p_{T} #bar{d} primary gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
144-
QAHistos.add("MCGen/ptAntiDeuteronFromBeautyHadron", "p_{T} #bar{d} from beauty gen hadrons;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
143+
qaHistos.add("MCGen/ptAntiDeuteronPrimary", "p_{T} #bar{d} primary gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
144+
qaHistos.add("MCGen/ptAntiDeuteronFromBeautyHadron", "p_{T} #bar{d} from beauty gen hadrons;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
145145

146-
QAHistos.add("MCGen/hMotherPdgCode", "PDG code of mother, gen level;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}});
147-
QAHistos.add("MCGen/ctauBminus", "ctau of B^{-}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
148-
QAHistos.add("MCGen/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
146+
qaHistos.add("MCGen/hMotherPdgCode", "PDG code of mother, gen level;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}});
147+
qaHistos.add("MCGen/ctauBminus", "ctau of B^{-}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
148+
qaHistos.add("MCGen/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, gen level;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
149149
// MC reco/MC-anchored histograms
150-
QAHistos.add("MCReco/ptAntiDeuteronFromBminus",
150+
qaHistos.add("MCReco/ptAntiDeuteronFromBminus",
151151
"p_{T} #bar{d} from B^{-} reco/MC anchored;p_{T} (GeV/c);Counts",
152152
HistType::kTH1F, {ptAxis});
153153

154-
QAHistos.add("MCReco/ptAntiDeuteronFromAntiLambdaB",
154+
qaHistos.add("MCReco/ptAntiDeuteronFromAntiLambdaB",
155155
"p_{T} #bar{d} from #bar{#Lambda}_{b} reco/MC anchored;p_{T} (GeV/c);Counts",
156156
HistType::kTH1F, {ptAxis});
157-
QAHistos.add("MCGen/ptAntiDeuteronFromBminus", "p_{T} #bar{d} from B^{-} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
158-
QAHistos.add("MCGen/ptAntiDeuteronFromAntiLambdaB", "p_{T} #bar{d} from #bar{#Lambda}_{b} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
159-
QAHistos.add("MCReco/hDCAxy-Primary", "DCAxy primary reco/MC anchored;DCA xy (cm);Counts", {HistType::kTH1D, {{400, -0.2f, 0.2f, "DCA xy (cm)"}}});
160-
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)"}}});
161-
QAHistos.add("MCReco/hMotherPdgCode", "PDG code of mother, reco/MC anchored;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}});
162-
QAHistos.add("MCReco/ctauBminus", "ctau of B^{-}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
163-
QAHistos.add("MCReco/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
157+
qaHistos.add("MCGen/ptAntiDeuteronFromBminus", "p_{T} #bar{d} from B^{-} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
158+
qaHistos.add("MCGen/ptAntiDeuteronFromAntiLambdaB", "p_{T} #bar{d} from #bar{#Lambda}_{b} gen;p_{T} (GeV/c);Counts", HistType::kTH1F, {ptAxis});
159+
qaHistos.add("MCReco/hDCAxy-Primary", "DCAxy primary reco/MC anchored;DCA xy (cm);Counts", {HistType::kTH1D, {{400, -0.2f, 0.2f, "DCA xy (cm)"}}});
160+
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)"}}});
161+
qaHistos.add("MCReco/hMotherPdgCode", "PDG code of mother, reco/MC anchored;PDG code;Counts", HistType::kTH1I, {{12000, -6000, 6000}});
162+
qaHistos.add("MCReco/ctauBminus", "ctau of B^{-}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
163+
qaHistos.add("MCReco/ctauAntiLambdaB", "ctau of #bar{#Lambda}_{b}, reco/MC anchored;ctau (#mu m);Counts", HistType::kTH1F, {{25, 0., 2000.f}});
164164

165165
hProcessedEvents->GetXaxis()->SetBinLabel(1, "Events processed");
166166
hProcessedEvents->GetXaxis()->SetBinLabel(2, "ZORRO");
@@ -173,7 +173,7 @@ struct HfTaskDeuteronFromLb {
173173
if (mRunNumber == bc.runNumber()) {
174174
return;
175175
}
176-
d_bz = 0.f;
176+
dBz = 0.f;
177177

178178
if (applySkimming) {
179179
zorro.initCCDB(ccdb.service, bc.runNumber(), bc.timestamp(), cfgSkimming.value);
@@ -184,30 +184,39 @@ struct HfTaskDeuteronFromLb {
184184
template <typename T1>
185185
bool passedSingleTrackSelection(const T1& track)
186186
{
187-
if (std::abs(track.eta()) > cfgEta)
187+
if (std::abs(track.eta()) > cfgEta) {
188188
return false;
189-
if (std::abs(track.dcaXY()) < cfgDCAmin || std::abs(track.dcaXY()) > cfgDCAmax)
190-
return false;
191-
if (!track.hasITS())
189+
}
190+
if (!track.hasITS()) {
192191
return false;
193-
if (!track.hasTPC())
192+
}
193+
if (!track.hasTPC()) {
194194
return false;
195-
if (!track.hasTOF())
195+
}
196+
if (!track.hasTOF()) {
196197
return false;
197-
if (track.tpcNClsFound() < cfgTPCNclsFound)
198+
}
199+
if (track.tpcNClsFound() < cfgTPCNclsFound) {
198200
return false;
199-
if (track.tpcChi2NCl() > cfgTPCChi2Ncl)
201+
}
202+
if (track.tpcChi2NCl() > cfgTPCChi2Ncl) {
200203
return false;
201-
if (track.itsChi2NCl() > cfgITSChi2Ncl)
204+
}
205+
if (track.itsChi2NCl() > cfgITSChi2Ncl) {
202206
return false;
203-
if (track.itsNCls() < cfgITScls)
207+
}
208+
if (track.itsNCls() < cfgITScls) {
204209
return false;
205-
if (track.pt() > cfgMaxPt)
210+
}
211+
if (track.pt() > cfgMaxPt) {
206212
return false;
207-
if (track.pt() < cfgMinPt)
213+
}
214+
if (track.pt() < cfgMinPt) {
208215
return false;
209-
if (track.sign() > 0)
216+
}
217+
if (track.sign() > 0) {
210218
return false;
219+
}
211220

212221
return true;
213222
}
@@ -219,7 +228,7 @@ struct HfTaskDeuteronFromLb {
219228
{
220229
for (const auto& collision : collisions) {
221230
if (mCurrentRun != collision.bc_as<o2::aod::BCsWithTimestamps>().runNumber()) {
222-
o2::parameters::GRPMagField* grpo = ccdb->getForTimeStamp<o2::parameters::GRPMagField>("GLO/Config/GRPMagField", collision.bc_as<o2::aod::BCsWithTimestamps>().timestamp());
231+
auto* grpo = ccdb->getForTimeStamp<o2::parameters::GRPMagField>("GLO/Config/GRPMagField", collision.bc_as<o2::aod::BCsWithTimestamps>().timestamp());
223232
o2::base::Propagator::initFieldFromGRP(grpo);
224233
mCurrentRun = collision.bc_as<o2::aod::BCsWithTimestamps>().runNumber();
225234
}
@@ -241,7 +250,7 @@ struct HfTaskDeuteronFromLb {
241250
continue;
242251
}
243252
hProcessedEvents->Fill(3.5);
244-
QAHistos.fill(HIST("Data/hVtxZ"), collision.posZ());
253+
qaHistos.fill(HIST("Data/hVtxZ"), collision.posZ());
245254

246255
const auto& trackIdsThisCollision = trackIndices.sliceBy(trackIndicesPerCollision, collision.globalIndex());
247256

@@ -254,6 +263,10 @@ struct HfTaskDeuteronFromLb {
254263
o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackPar, 2.f, noMatCorr, &dca);
255264
}
256265

266+
if (std::abs(dca[0]) < cfgDCAmin || std::abs(dca[0]) > cfgDCAmax) {
267+
continue;
268+
}
269+
257270
if (!passedSingleTrackSelection(track)) {
258271
continue;
259272
}
@@ -263,21 +276,21 @@ struct HfTaskDeuteronFromLb {
263276

264277
if (track.pt() < ptThresholdPid) {
265278
if (isTPCDe) {
266-
QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt());
267-
QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta());
268-
QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]);
269-
QAHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]);
270-
QAHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe());
271-
QAHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe());
279+
qaHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt());
280+
qaHistos.fill(HIST("Data/etaAntideuteron"), track.eta());
281+
qaHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]);
282+
qaHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]);
283+
qaHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe());
284+
qaHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe());
272285
}
273286
} else {
274287
if (isTPCDe && isTOFDe) {
275-
QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt());
276-
QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta());
277-
QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]);
278-
QAHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]);
279-
QAHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe());
280-
QAHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe());
288+
qaHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt());
289+
qaHistos.fill(HIST("Data/etaAntideuteron"), track.eta());
290+
qaHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]);
291+
qaHistos.fill(HIST("Data/hDCAzVsPt"), track.pt(), dca[1]);
292+
qaHistos.fill(HIST("Data/hnSigmaTPCVsPt"), track.pt(), track.tpcNSigmaDe());
293+
qaHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe());
281294
}
282295
}
283296
}
@@ -288,6 +301,9 @@ struct HfTaskDeuteronFromLb {
288301
void processMC(MCCollisionCandidates::iterator const& collision, MCTrackCandidates const& tracks, o2::aod::McParticles const&)
289302
{
290303
for (const auto& track : tracks) {
304+
if (std::abs(track.dcaXY()) < cfgDCAmin || std::abs(track.dcaXY()) > cfgDCAmax) {
305+
continue;
306+
}
291307
if (!passedSingleTrackSelection(track)) {
292308
continue;
293309
}
@@ -305,7 +321,7 @@ struct HfTaskDeuteronFromLb {
305321
int motherPdg = 0;
306322
if (separateAntideuterons) {
307323
for (const auto& mom : mcParticle.mothers_as<o2::aod::McParticles>()) {
308-
QAHistos.fill(HIST("MCReco/hMotherPdgCode"), mom.pdgCode());
324+
qaHistos.fill(HIST("MCReco/hMotherPdgCode"), mom.pdgCode());
309325
if (mom.pdgCode() == pdgCodeBeautyMeson || mom.pdgCode() == pdgCodeBeautyBaryon) {
310326
motherPdg = mom.pdgCode();
311327
double dx = mcParticle.vx() - collision.posX();
@@ -320,24 +336,24 @@ struct HfTaskDeuteronFromLb {
320336
ctauBeauty *= 1.e4;
321337

322338
if (motherPdg == pdgCodeBeautyMeson) {
323-
QAHistos.fill(HIST("MCReco/ctauBminus"), ctauBeauty);
339+
qaHistos.fill(HIST("MCReco/ctauBminus"), ctauBeauty);
324340
} else if (motherPdg == pdgCodeBeautyBaryon) {
325-
QAHistos.fill(HIST("MCReco/ctauAntiLambdaB"), ctauBeauty);
341+
qaHistos.fill(HIST("MCReco/ctauAntiLambdaB"), ctauBeauty);
326342
}
327343
}
328344
break;
329345
}
330346
}
331347
}
332348
if (motherPdg == pdgCodeBeautyMeson) {
333-
QAHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY());
334-
QAHistos.fill(HIST("MCReco/ptAntiDeuteronFromBminus"), track.pt());
349+
qaHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY());
350+
qaHistos.fill(HIST("MCReco/ptAntiDeuteronFromBminus"), track.pt());
335351
} else if (motherPdg == pdgCodeBeautyBaryon) {
336-
QAHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY());
337-
QAHistos.fill(HIST("MCReco/ptAntiDeuteronFromAntiLambdaB"), track.pt());
352+
qaHistos.fill(HIST("MCReco/hDCAxy-FromBeautyHadron"), track.dcaXY());
353+
qaHistos.fill(HIST("MCReco/ptAntiDeuteronFromAntiLambdaB"), track.pt());
338354
} else {
339-
QAHistos.fill(HIST("MCReco/hDCAxy-Primary"), track.dcaXY());
340-
QAHistos.fill(HIST("MCReco/ptAntiDeuteronPrimary"), track.pt());
355+
qaHistos.fill(HIST("MCReco/hDCAxy-Primary"), track.dcaXY());
356+
qaHistos.fill(HIST("MCReco/ptAntiDeuteronPrimary"), track.pt());
341357
}
342358
}
343359
}
@@ -355,9 +371,9 @@ struct HfTaskDeuteronFromLb {
355371
if (isSelectedBeautyHadron(mcParticle.pdgCode())) {
356372
if (std::abs(mcParticle.y()) <= rapidityCut) {
357373
if (mcParticle.pdgCode() == pdgCodeBeautyMeson) {
358-
QAHistos.fill(HIST("MCGen/ptGeneratedBminus"), mcParticle.pt());
374+
qaHistos.fill(HIST("MCGen/ptGeneratedBminus"), mcParticle.pt());
359375
} else if (mcParticle.pdgCode() == pdgCodeBeautyBaryon) {
360-
QAHistos.fill(HIST("MCGen/ptGeneratedAntiLambdaB"), mcParticle.pt());
376+
qaHistos.fill(HIST("MCGen/ptGeneratedAntiLambdaB"), mcParticle.pt());
361377
}
362378
}
363379

@@ -376,9 +392,9 @@ struct HfTaskDeuteronFromLb {
376392
ctauBeauty *= 1.e4; // cm -> micrometers
377393

378394
if (mcParticle.pdgCode() == pdgCodeBeautyMeson) {
379-
QAHistos.fill(HIST("MCGen/ctauBminus"), ctauBeauty);
395+
qaHistos.fill(HIST("MCGen/ctauBminus"), ctauBeauty);
380396
} else if (mcParticle.pdgCode() == pdgCodeBeautyBaryon) {
381-
QAHistos.fill(HIST("MCGen/ctauAntiLambdaB"), ctauBeauty);
397+
qaHistos.fill(HIST("MCGen/ctauAntiLambdaB"), ctauBeauty);
382398
}
383399
}
384400
break;
@@ -396,7 +412,7 @@ struct HfTaskDeuteronFromLb {
396412
if (mcParticle.has_mothers()) {
397413
for (const auto& mom : mcParticle.mothers_as<o2::aod::McParticles>()) {
398414

399-
QAHistos.fill(HIST("MCGen/hMotherPdgCode"), mom.pdgCode());
415+
qaHistos.fill(HIST("MCGen/hMotherPdgCode"), mom.pdgCode());
400416

401417
if (mom.pdgCode() == pdgCodeBeautyMeson) {
402418
motherPdg = mom.pdgCode();
@@ -411,11 +427,11 @@ struct HfTaskDeuteronFromLb {
411427
}
412428

413429
if (motherPdg == pdgCodeBeautyMeson) {
414-
QAHistos.fill(HIST("MCGen/ptAntiDeuteronFromBminus"), mcParticle.pt());
430+
qaHistos.fill(HIST("MCGen/ptAntiDeuteronFromBminus"), mcParticle.pt());
415431
} else if (motherPdg == pdgCodeBeautyBaryon) {
416-
QAHistos.fill(HIST("MCGen/ptAntiDeuteronFromAntiLambdaB"), mcParticle.pt());
432+
qaHistos.fill(HIST("MCGen/ptAntiDeuteronFromAntiLambdaB"), mcParticle.pt());
417433
} else if (mcParticle.isPhysicalPrimary()) {
418-
QAHistos.fill(HIST("MCGen/ptAntiDeuteronPrimary"), mcParticle.pt());
434+
qaHistos.fill(HIST("MCGen/ptAntiDeuteronPrimary"), mcParticle.pt());
419435
}
420436
}
421437
}

0 commit comments

Comments
 (0)