@@ -78,7 +78,7 @@ 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
8484 framework::Service<ccdb::BasicCCDBManager> ccdb{};
@@ -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 );
@@ -250,7 +250,7 @@ struct HfTaskDeuteronFromLb {
250250 continue ;
251251 }
252252 hProcessedEvents->Fill (3.5 );
253- QAHistos .fill (HIST (" Data/hVtxZ" ), collision.posZ ());
253+ qaHistos .fill (HIST (" Data/hVtxZ" ), collision.posZ ());
254254
255255 const auto & trackIdsThisCollision = trackIndices.sliceBy (trackIndicesPerCollision, collision.globalIndex ());
256256
@@ -276,21 +276,21 @@ struct HfTaskDeuteronFromLb {
276276
277277 if (track.pt () < ptThresholdPid) {
278278 if (isTPCDe) {
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 ());
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 ());
285285 }
286286 } else {
287287 if (isTPCDe && isTOFDe) {
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 ());
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 ());
294294 }
295295 }
296296 }
@@ -321,7 +321,7 @@ struct HfTaskDeuteronFromLb {
321321 int motherPdg = 0 ;
322322 if (separateAntideuterons) {
323323 for (const auto & mom : mcParticle.mothers_as <o2::aod::McParticles>()) {
324- QAHistos .fill (HIST (" MCReco/hMotherPdgCode" ), mom.pdgCode ());
324+ qaHistos .fill (HIST (" MCReco/hMotherPdgCode" ), mom.pdgCode ());
325325 if (mom.pdgCode () == pdgCodeBeautyMeson || mom.pdgCode () == pdgCodeBeautyBaryon) {
326326 motherPdg = mom.pdgCode ();
327327 double dx = mcParticle.vx () - collision.posX ();
@@ -336,24 +336,24 @@ struct HfTaskDeuteronFromLb {
336336 ctauBeauty *= 1 .e4 ;
337337
338338 if (motherPdg == pdgCodeBeautyMeson) {
339- QAHistos .fill (HIST (" MCReco/ctauBminus" ), ctauBeauty);
339+ qaHistos .fill (HIST (" MCReco/ctauBminus" ), ctauBeauty);
340340 } else if (motherPdg == pdgCodeBeautyBaryon) {
341- QAHistos .fill (HIST (" MCReco/ctauAntiLambdaB" ), ctauBeauty);
341+ qaHistos .fill (HIST (" MCReco/ctauAntiLambdaB" ), ctauBeauty);
342342 }
343343 }
344344 break ;
345345 }
346346 }
347347 }
348348 if (motherPdg == pdgCodeBeautyMeson) {
349- QAHistos .fill (HIST (" MCReco/hDCAxy-FromBeautyHadron" ), track.dcaXY ());
350- QAHistos .fill (HIST (" MCReco/ptAntiDeuteronFromBminus" ), track.pt ());
349+ qaHistos .fill (HIST (" MCReco/hDCAxy-FromBeautyHadron" ), track.dcaXY ());
350+ qaHistos .fill (HIST (" MCReco/ptAntiDeuteronFromBminus" ), track.pt ());
351351 } else if (motherPdg == pdgCodeBeautyBaryon) {
352- QAHistos .fill (HIST (" MCReco/hDCAxy-FromBeautyHadron" ), track.dcaXY ());
353- QAHistos .fill (HIST (" MCReco/ptAntiDeuteronFromAntiLambdaB" ), track.pt ());
352+ qaHistos .fill (HIST (" MCReco/hDCAxy-FromBeautyHadron" ), track.dcaXY ());
353+ qaHistos .fill (HIST (" MCReco/ptAntiDeuteronFromAntiLambdaB" ), track.pt ());
354354 } else {
355- QAHistos .fill (HIST (" MCReco/hDCAxy-Primary" ), track.dcaXY ());
356- QAHistos .fill (HIST (" MCReco/ptAntiDeuteronPrimary" ), track.pt ());
355+ qaHistos .fill (HIST (" MCReco/hDCAxy-Primary" ), track.dcaXY ());
356+ qaHistos .fill (HIST (" MCReco/ptAntiDeuteronPrimary" ), track.pt ());
357357 }
358358 }
359359 }
@@ -371,9 +371,9 @@ struct HfTaskDeuteronFromLb {
371371 if (isSelectedBeautyHadron (mcParticle.pdgCode ())) {
372372 if (std::abs (mcParticle.y ()) <= rapidityCut) {
373373 if (mcParticle.pdgCode () == pdgCodeBeautyMeson) {
374- QAHistos .fill (HIST (" MCGen/ptGeneratedBminus" ), mcParticle.pt ());
374+ qaHistos .fill (HIST (" MCGen/ptGeneratedBminus" ), mcParticle.pt ());
375375 } else if (mcParticle.pdgCode () == pdgCodeBeautyBaryon) {
376- QAHistos .fill (HIST (" MCGen/ptGeneratedAntiLambdaB" ), mcParticle.pt ());
376+ qaHistos .fill (HIST (" MCGen/ptGeneratedAntiLambdaB" ), mcParticle.pt ());
377377 }
378378 }
379379
@@ -392,9 +392,9 @@ struct HfTaskDeuteronFromLb {
392392 ctauBeauty *= 1 .e4 ; // cm -> micrometers
393393
394394 if (mcParticle.pdgCode () == pdgCodeBeautyMeson) {
395- QAHistos .fill (HIST (" MCGen/ctauBminus" ), ctauBeauty);
395+ qaHistos .fill (HIST (" MCGen/ctauBminus" ), ctauBeauty);
396396 } else if (mcParticle.pdgCode () == pdgCodeBeautyBaryon) {
397- QAHistos .fill (HIST (" MCGen/ctauAntiLambdaB" ), ctauBeauty);
397+ qaHistos .fill (HIST (" MCGen/ctauAntiLambdaB" ), ctauBeauty);
398398 }
399399 }
400400 break ;
@@ -412,7 +412,7 @@ struct HfTaskDeuteronFromLb {
412412 if (mcParticle.has_mothers ()) {
413413 for (const auto & mom : mcParticle.mothers_as <o2::aod::McParticles>()) {
414414
415- QAHistos .fill (HIST (" MCGen/hMotherPdgCode" ), mom.pdgCode ());
415+ qaHistos .fill (HIST (" MCGen/hMotherPdgCode" ), mom.pdgCode ());
416416
417417 if (mom.pdgCode () == pdgCodeBeautyMeson) {
418418 motherPdg = mom.pdgCode ();
@@ -427,11 +427,11 @@ struct HfTaskDeuteronFromLb {
427427 }
428428
429429 if (motherPdg == pdgCodeBeautyMeson) {
430- QAHistos .fill (HIST (" MCGen/ptAntiDeuteronFromBminus" ), mcParticle.pt ());
430+ qaHistos .fill (HIST (" MCGen/ptAntiDeuteronFromBminus" ), mcParticle.pt ());
431431 } else if (motherPdg == pdgCodeBeautyBaryon) {
432- QAHistos .fill (HIST (" MCGen/ptAntiDeuteronFromAntiLambdaB" ), mcParticle.pt ());
432+ qaHistos .fill (HIST (" MCGen/ptAntiDeuteronFromAntiLambdaB" ), mcParticle.pt ());
433433 } else if (mcParticle.isPhysicalPrimary ()) {
434- QAHistos .fill (HIST (" MCGen/ptAntiDeuteronPrimary" ), mcParticle.pt ());
434+ qaHistos .fill (HIST (" MCGen/ptAntiDeuteronPrimary" ), mcParticle.pt ());
435435 }
436436 }
437437 }
0 commit comments