Skip to content

Commit a23acba

Browse files
authored
[PWGLF] Fix missing Bz initialisation in the DCAFitter (#9432)
1 parent c6b3382 commit a23acba

File tree

1 file changed

+7
-6
lines changed

1 file changed

+7
-6
lines changed

PWGLF/Tasks/Strangeness/nonPromptCascade.cxx

Lines changed: 7 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -194,7 +194,7 @@ struct NonPromptCascadeTask {
194194

195195
Service<o2::ccdb::BasicCCDBManager> ccdb;
196196
int mRunNumber = 0;
197-
float bz = 0.f;
197+
float mBz = 0.f;
198198

199199
HistogramRegistry registry{
200200
"registry",
@@ -252,6 +252,7 @@ struct NonPromptCascadeTask {
252252

253253
if (o2::parameters::GRPMagField* grpmag = ccdb->getForRun<o2::parameters::GRPMagField>(cfgGRPmagPath, mRunNumber)) {
254254
o2::base::Propagator::initFieldFromGRP(grpmag);
255+
mBz = static_cast<float>(grpmag->getNominalL3Field());
255256
}
256257

257258
if (static_cast<o2::base::Propagator::MatCorrType>(cfgMaterialCorrection.value) == o2::base::Propagator::MatCorrType::USEMatCorrLUT) {
@@ -297,7 +298,7 @@ struct NonPromptCascadeTask {
297298
auto trackCovTrk = getTrackParCov(track);
298299
o2::dataformats::DCA impactParameterTrk;
299300

300-
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovTrk, bz, 2.f, matCorr, &impactParameterTrk)) {
301+
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovTrk, mBz, 2.f, matCorr, &impactParameterTrk)) {
301302
if (protonTrack.hasTPC() && pionTrack.hasTPC()) {
302303
if (isOmega) {
303304
registry.fill(HIST("h_dca_Omega"), std::sqrt(impactParameterTrk.getR2()));
@@ -327,7 +328,7 @@ struct NonPromptCascadeTask {
327328

328329
auto trackCovBach = getTrackParCov(bachelor);
329330
o2::dataformats::DCA impactParameterBach;
330-
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovBach, bz, 2.f, matCorr, &impactParameterBach)) {
331+
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovBach, mBz, 2.f, matCorr, &impactParameterBach)) {
331332
if (isOmega) {
332333
if (bachelor.sign() < 0) {
333334
registry.fill(HIST("h_bachdcaxyM_Omega"), impactParameterBach.getY());
@@ -354,7 +355,7 @@ struct NonPromptCascadeTask {
354355

355356
auto trackCovNtrack = getTrackParCov(pionTrack);
356357
o2::dataformats::DCA impactParameterPiontrack;
357-
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovNtrack, bz, 2.f, matCorr, &impactParameterPiontrack)) {
358+
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovNtrack, mBz, 2.f, matCorr, &impactParameterPiontrack)) {
358359
if (isOmega) {
359360
registry.fill(HIST("h_ntrackdcavspt_Omega"), impactParameterPiontrack.getY(), pionTrack.pt());
360361
}
@@ -363,7 +364,7 @@ struct NonPromptCascadeTask {
363364

364365
auto trackCovPtrack = getTrackParCov(protonTrack);
365366
o2::dataformats::DCA impactParameterProtontrack;
366-
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovPtrack, bz, 2.f, matCorr, &impactParameterProtontrack)) {
367+
if (o2::base::Propagator::Instance()->propagateToDCA(primaryVertex, trackCovPtrack, mBz, 2.f, matCorr, &impactParameterProtontrack)) {
367368
if (isOmega) {
368369
registry.fill(HIST("h_ptrackdcavspt_Omega"), impactParameterProtontrack.getY(), protonTrack.pt());
369370
}
@@ -391,7 +392,7 @@ struct NonPromptCascadeTask {
391392
const auto primaryVertex = getPrimaryVertex(collision);
392393

393394
o2::vertexing::DCAFitterN<2> df2;
394-
df2.setBz(bz);
395+
df2.setBz(mBz);
395396
df2.setPropagateToPCA(propToDCA);
396397
df2.setMaxR(maxR);
397398
df2.setMaxDZIni(maxDZIni);

0 commit comments

Comments
 (0)