Skip to content

Commit db5c993

Browse files
vkuceravkucera
committed
Fix const
1 parent 4be695a commit db5c993

File tree

141 files changed

+1887
-1887
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

141 files changed

+1887
-1887
lines changed

‎PWGHF/Core/HfHelper.h‎

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -304,12 +304,12 @@ class HfHelper
304304
momK2 = candidate.pVectorProng2();
305305
}
306306

307-
ROOT::Math::PxPyPzMVector vecPi(momPi[0], momPi[1], momPi[2], o2::constants::physics::MassPiPlus);
308-
ROOT::Math::PxPyPzMVector vecK1(momK1[0], momK1[1], momK1[2], o2::constants::physics::MassKPlus);
309-
ROOT::Math::PxPyPzMVector vecK2(momK2[0], momK2[1], momK2[2], o2::constants::physics::MassKPlus);
310-
ROOT::Math::PxPyPzMVector vecPhi = vecK1 + vecK2;
307+
ROOT::Math::PxPyPzMVector const vecPi(momPi[0], momPi[1], momPi[2], o2::constants::physics::MassPiPlus);
308+
ROOT::Math::PxPyPzMVector const vecK1(momK1[0], momK1[1], momK1[2], o2::constants::physics::MassKPlus);
309+
ROOT::Math::PxPyPzMVector const vecK2(momK2[0], momK2[1], momK2[2], o2::constants::physics::MassKPlus);
310+
ROOT::Math::PxPyPzMVector const vecPhi = vecK1 + vecK2;
311311

312-
ROOT::Math::Boost boostToPhiRestFrame(vecPhi.BoostToCM());
312+
ROOT::Math::Boost const boostToPhiRestFrame(vecPhi.BoostToCM());
313313
auto momPiPhiRestFrame = boostToPhiRestFrame(vecPi).Vect();
314314
auto momK1PhiRestFrame = boostToPhiRestFrame(vecK1).Vect();
315315

@@ -522,8 +522,8 @@ class HfHelper
522522
phiPi = RecoDecay::phi(candidate.pxProng2(), candidate.pyProng2());
523523
}
524524

525-
double deltaEta = etaJpsi - etaPi;
526-
double deltaPhi = RecoDecay::constrainAngle(phiJpsi - phiPi, -o2::constants::math::PI);
525+
double const deltaEta = etaJpsi - etaPi;
526+
double const deltaPhi = RecoDecay::constrainAngle(phiJpsi - phiPi, -o2::constants::math::PI);
527527

528528
return RecoDecay::sqrtSumOfSquares(deltaEta, deltaPhi);
529529
}

‎PWGHF/D2H/Macros/HFInvMassFitter.cxx‎

Lines changed: 35 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -334,7 +334,7 @@ void HFInvMassFitter::doFit()
334334
}
335335
mTotalPdf->plotOn(mInvMassFrame, Name("Tot_c"));
336336
mReflPdf = new RooAddPdf("mReflPdf", "reflection fit function", RooArgList(*reflPdf), RooArgList(*mRooNRefl));
337-
RooAddPdf reflBkgPdf("reflBkgPdf", "reflBkgPdf", RooArgList(*bkgPdf, *reflPdf), RooArgList(*mRooNBkg, *mRooNRefl));
337+
RooAddPdf const reflBkgPdf("reflBkgPdf", "reflBkgPdf", RooArgList(*bkgPdf, *reflPdf), RooArgList(*mRooNBkg, *mRooNRefl));
338338
reflBkgPdf.plotOn(mInvMassFrame, Normalization(1.0, RooAbsReal::RelativeExpected), LineStyle(7), LineColor(kRed + 1), Name("ReflBkg_c"));
339339
plotBkg(mTotalPdf); // plot bkg pdf in total pdf
340340
plotRefl(mTotalPdf); // plot reflection in total pdf
@@ -385,30 +385,30 @@ void HFInvMassFitter::fillWorkspace(RooWorkspace& workspace) const
385385
workspace.import(*bkgFuncExpo);
386386
delete bkgFuncExpo;
387387
// bkg poly1
388-
RooRealVar polyParam0("polyParam0", "Parameter of Poly function", 0.5, -5., 5.);
389-
RooRealVar polyParam1("polyParam1", "Parameter of Poly function", 0.2, -5., 5.);
388+
RooRealVar const polyParam0("polyParam0", "Parameter of Poly function", 0.5, -5., 5.);
389+
RooRealVar const polyParam1("polyParam1", "Parameter of Poly function", 0.2, -5., 5.);
390390
RooAbsPdf* bkgFuncPoly1 = new RooPolynomial("bkgFuncPoly1", "background fit function", mass, RooArgSet(polyParam0, polyParam1));
391391
workspace.import(*bkgFuncPoly1);
392392
delete bkgFuncPoly1;
393393
// bkg poly2
394-
RooRealVar polyParam2("polyParam2", "Parameter of Poly function", 0.2, -5., 5.);
394+
RooRealVar const polyParam2("polyParam2", "Parameter of Poly function", 0.2, -5., 5.);
395395
RooAbsPdf* bkgFuncPoly2 = new RooPolynomial("bkgFuncPoly2", "background fit function", mass, RooArgSet(polyParam0, polyParam1, polyParam2));
396396
workspace.import(*bkgFuncPoly2);
397397
delete bkgFuncPoly2;
398398
// bkg poly3
399-
RooRealVar polyParam3("polyParam3", "Parameter of Poly function", 0.2, -1., 1.);
399+
RooRealVar const polyParam3("polyParam3", "Parameter of Poly function", 0.2, -1., 1.);
400400
RooAbsPdf* bkgFuncPoly3 = new RooPolynomial("bkgFuncPoly3", "background pdf", mass, RooArgSet(polyParam0, polyParam1, polyParam2, polyParam3));
401401
workspace.import(*bkgFuncPoly3);
402402
delete bkgFuncPoly3;
403403
// bkg power law
404-
RooRealVar powParam1("powParam1", "Parameter of Pow function", TDatabasePDG::Instance()->GetParticle("pi+")->Mass());
405-
RooRealVar powParam2("powParam2", "Parameter of Pow function", 1., -10, 10);
404+
RooRealVar const powParam1("powParam1", "Parameter of Pow function", TDatabasePDG::Instance()->GetParticle("pi+")->Mass());
405+
RooRealVar const powParam2("powParam2", "Parameter of Pow function", 1., -10, 10);
406406
RooAbsPdf* bkgFuncPow = new RooGenericPdf("bkgFuncPow", "bkgFuncPow", "(mass-powParam1)^powParam2", RooArgSet(mass, powParam1, powParam2));
407407
workspace.import(*bkgFuncPow);
408408
delete bkgFuncPow;
409409
// pow * exp
410-
RooRealVar powExpoParam1("powExpoParam1", "Parameter of PowExpo function", 1. / 2.);
411-
RooRealVar powExpoParam2("powExpoParam2", "Parameter of PowExpo function", 1, -10, 10);
410+
RooRealVar const powExpoParam1("powExpoParam1", "Parameter of PowExpo function", 1. / 2.);
411+
RooRealVar const powExpoParam2("powExpoParam2", "Parameter of PowExpo function", 1, -10, 10);
412412
RooRealVar massPi("massPi", "mass of pion", TDatabasePDG::Instance()->GetParticle("pi+")->Mass());
413413
RooFormulaVar powExpoParam3("powExpoParam3", "powExpoParam1 + 1", RooArgList(powExpoParam1));
414414
RooFormulaVar powExpoParam4("powExpoParam4", "1./powExpoParam2", RooArgList(powExpoParam2));
@@ -453,8 +453,8 @@ void HFInvMassFitter::fillWorkspace(RooWorkspace& workspace) const
453453
sigmaDoubleGaus.setVal(mSigmaSgnDoubleGaus);
454454
sigmaDoubleGaus.setConstant(kTRUE);
455455
}
456-
RooGaussian gaus1("gaus1", "gaus1", mass, mean, sigma);
457-
RooGaussian gaus2("gaus2", "gaus2", mass, mean, sigmaDoubleGaus);
456+
RooGaussian const gaus1("gaus1", "gaus1", mass, mean, sigma);
457+
RooGaussian const gaus2("gaus2", "gaus2", mass, mean, sigmaDoubleGaus);
458458
RooRealVar fracDoubleGaus("fracDoubleGaus", "frac of two gauss", mFracDoubleGaus, 0, 1.);
459459
if (mFixedFracDoubleGaus) {
460460
fracDoubleGaus.setVal(mFracDoubleGaus);
@@ -478,8 +478,8 @@ void HFInvMassFitter::fillWorkspace(RooWorkspace& workspace) const
478478
sigma.setMin(mSigmaSgn * (1 - mParamSgn));
479479
}
480480
RooRealVar sigmaDoubleGausRatio("sigmaDoubleGausRatio", "sigmaDoubleGausRatio", sigma.getVal() * ratio.getVal());
481-
RooGaussian gausRatio1("gausRatio1", "gausratio1", mass, mean, sigma);
482-
RooGaussian gausRatio2("gausRatio2", "gausratio2", mass, mean, sigmaDoubleGausRatio);
481+
RooGaussian const gausRatio1("gausRatio1", "gausratio1", mass, mean, sigma);
482+
RooGaussian const gausRatio2("gausRatio2", "gausratio2", mass, mean, sigmaDoubleGausRatio);
483483
RooRealVar fracDoubleGausRatio("fracDoubleGausRatio", "fraction of two gauss ratio", 0.5, 0, 1.);
484484
if (mFixedFracDoubleGaus) {
485485
fracDoubleGausRatio.setVal(mFracDoubleGaus);
@@ -507,9 +507,9 @@ void HFInvMassFitter::fillWorkspace(RooWorkspace& workspace) const
507507
sigmaSec.setMax(mSecSigma * (1 + mParamSgn));
508508
sigmaSec.setMin(mSecSigma * (1 - mParamSgn));
509509
}
510-
RooGaussian gausSec1("gausSec1", "gausSec1", mass, mean, sigmaSec);
511-
RooGaussian gausSec2("gausSec2", "gausSec2", mass, meanSec, sigmaSec);
512-
RooRealVar fracSec("fracSec", "frac of two peak", 0.5, 0, 1.);
510+
RooGaussian const gausSec1("gausSec1", "gausSec1", mass, mean, sigmaSec);
511+
RooGaussian const gausSec2("gausSec2", "gausSec2", mass, meanSec, sigmaSec);
512+
RooRealVar const fracSec("fracSec", "frac of two peak", 0.5, 0, 1.);
513513
RooAbsPdf* sgnFuncDoublePeak = new RooAddPdf("sgnFuncDoublePeak", "signal pdf", RooArgList(gausSec1, gausSec2), fracSec);
514514
workspace.import(*sgnFuncDoublePeak);
515515
delete sgnFuncDoublePeak;
@@ -530,24 +530,24 @@ void HFInvMassFitter::fillWorkspace(RooWorkspace& workspace) const
530530
meanReflDoubleGaus.setMin(mMassReflLowLimit);
531531
}
532532
RooRealVar sigmaReflDoubleGaus("sigmaReflDoubleGaus", "sigmaReflDoubleGaus", 0.012, 0.0, 0.25);
533-
RooGaussian gausRefl1("gausRefl1", "gausRefl1", mass, meanRefl, sigmaRefl);
534-
RooGaussian gausRefl2("gausRefl2", "gausRefl2", mass, meanReflDoubleGaus, sigmaReflDoubleGaus);
535-
RooRealVar fracRefl("fracRefl", "frac of two gauss", 0.5, 0, 1.);
533+
RooGaussian const gausRefl1("gausRefl1", "gausRefl1", mass, meanRefl, sigmaRefl);
534+
RooGaussian const gausRefl2("gausRefl2", "gausRefl2", mass, meanReflDoubleGaus, sigmaReflDoubleGaus);
535+
RooRealVar const fracRefl("fracRefl", "frac of two gauss", 0.5, 0, 1.);
536536
RooAbsPdf* reflFuncDoubleGaus = new RooAddPdf("reflFuncDoubleGaus", "reflection pdf", RooArgList(gausRefl1, gausRefl2), fracRefl);
537537
workspace.import(*reflFuncDoubleGaus);
538538
delete reflFuncDoubleGaus;
539539
// reflection poly3
540-
RooRealVar polyReflParam0("polyReflParam0", "polyReflParam0", 0.5, -1., 1.);
541-
RooRealVar polyReflParam1("polyReflParam1", "polyReflParam1", 0.2, -1., 1.);
542-
RooRealVar polyReflParam2("polyReflParam2", "polyReflParam2", 0.2, -1., 1.);
543-
RooRealVar polyReflParam3("polyReflParam3", "polyReflParam3", 0.2, -1., 1.);
540+
RooRealVar const polyReflParam0("polyReflParam0", "polyReflParam0", 0.5, -1., 1.);
541+
RooRealVar const polyReflParam1("polyReflParam1", "polyReflParam1", 0.2, -1., 1.);
542+
RooRealVar const polyReflParam2("polyReflParam2", "polyReflParam2", 0.2, -1., 1.);
543+
RooRealVar const polyReflParam3("polyReflParam3", "polyReflParam3", 0.2, -1., 1.);
544544
RooAbsPdf* reflFuncPoly3 = new RooPolynomial("reflFuncPoly3", "reflection PDF", mass, RooArgSet(polyReflParam0, polyReflParam1, polyReflParam2, polyReflParam3));
545545
workspace.import(*reflFuncPoly3);
546546
delete reflFuncPoly3;
547547
// reflection poly6
548-
RooRealVar polyReflParam4("polyReflParam4", "polyReflParam4", 0.2, -1., 1.);
549-
RooRealVar polyReflParam5("polyReflParam5", "polyReflParam5", 0.2, -1., 1.);
550-
RooRealVar polyReflParam6("polyReflParam6", "polyReflParam6", 0.2, -1., 1.);
548+
RooRealVar const polyReflParam4("polyReflParam4", "polyReflParam4", 0.2, -1., 1.);
549+
RooRealVar const polyReflParam5("polyReflParam5", "polyReflParam5", 0.2, -1., 1.);
550+
RooRealVar const polyReflParam6("polyReflParam6", "polyReflParam6", 0.2, -1., 1.);
551551
RooAbsPdf* reflFuncPoly6 = new RooPolynomial("reflFuncPoly6", "reflection pdf", mass, RooArgSet(polyReflParam0, polyReflParam1, polyReflParam2, polyReflParam3, polyReflParam4, polyReflParam5, polyReflParam6));
552552
workspace.import(*reflFuncPoly6);
553553
delete reflFuncPoly6;
@@ -643,8 +643,8 @@ void HFInvMassFitter::highlightPeakRegion(const RooPlot* plot, Color_t color, Wi
643643
if (!mHighlightPeakRegion) {
644644
return;
645645
}
646-
double yMin = plot->GetMinimum();
647-
double yMax = plot->GetMaximum();
646+
double const yMin = plot->GetMinimum();
647+
double const yMax = plot->GetMaximum();
648648
const Double_t mean = mRooMeanSgn->getVal();
649649
const Double_t sigma = mRooSigmaSgn->getVal();
650650
const Double_t minForSgn = mean - mNSigmaForSidebands * sigma;
@@ -716,20 +716,20 @@ void HFInvMassFitter::calculateSignificance(Double_t& significance, Double_t& er
716716
calculateSignal(signal, errSignal);
717717
Double_t bkg, errBkg;
718718
calculateBackground(bkg, errBkg);
719-
Double_t sgnErrSquare = errSignal * errSignal;
720-
Double_t bkgErrSquare = errBkg * errBkg;
721-
Double_t totalSgnBkg = signal + bkg;
719+
Double_t const sgnErrSquare = errSignal * errSignal;
720+
Double_t const bkgErrSquare = errBkg * errBkg;
721+
Double_t const totalSgnBkg = signal + bkg;
722722
significance = signal / std::sqrt(signal + bkg);
723723
errSignificance = significance * std::sqrt((sgnErrSquare + bkgErrSquare) / (mNSigmaForSidebands * totalSgnBkg * totalSgnBkg) + (bkg / totalSgnBkg) * (sgnErrSquare / signal / signal));
724724
}
725725

726726
// estimate Signal
727727
void HFInvMassFitter::checkForSignal(Double_t& estimatedSignal)
728728
{
729-
Double_t minForSgn = mMass - 4 * mSigmaSgn;
730-
Double_t maxForSgn = mMass + 4 * mSigmaSgn;
731-
Int_t binForMinSgn = mHistoInvMass->FindBin(minForSgn);
732-
Int_t binForMaxSgn = mHistoInvMass->FindBin(maxForSgn);
729+
Double_t const minForSgn = mMass - 4 * mSigmaSgn;
730+
Double_t const maxForSgn = mMass + 4 * mSigmaSgn;
731+
Int_t const binForMinSgn = mHistoInvMass->FindBin(minForSgn);
732+
Int_t const binForMaxSgn = mHistoInvMass->FindBin(maxForSgn);
733733

734734
Double_t sum = 0;
735735
for (Int_t i = binForMinSgn; i <= binForMaxSgn; i++) {

‎PWGHF/D2H/Macros/runMassFitter.C‎

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -62,7 +62,7 @@ void readArray(const Value& jsonArray, std::vector<ValueType>& output)
6262

6363
void parseStringArray(const Value& jsonArray, std::vector<std::string>& output)
6464
{
65-
size_t arrayLength = jsonArray.Size();
65+
size_t const arrayLength = jsonArray.Size();
6666
for (size_t i = 0; i < arrayLength; i++) {
6767
if (jsonArray[i].IsString()) {
6868
output.emplace_back(jsonArray[i].GetString());
@@ -87,11 +87,11 @@ int runMassFitter(const TString& configFileName)
8787
config.ParseStream(is);
8888
fclose(configFile);
8989

90-
Bool_t isMc = config["IsMC"].GetBool();
91-
TString inputFileName = config["InFileName"].GetString();
92-
TString reflFileName = config["ReflFileName"].GetString();
90+
Bool_t const isMc = config["IsMC"].GetBool();
91+
TString const inputFileName = config["InFileName"].GetString();
92+
TString const reflFileName = config["ReflFileName"].GetString();
9393
TString outputFileName = config["OutFileName"].GetString();
94-
TString particleName = config["Particle"].GetString();
94+
TString const particleName = config["Particle"].GetString();
9595

9696
std::vector<std::string> inputHistoName;
9797
std::vector<std::string> promptHistoName;
@@ -166,8 +166,8 @@ int runMassFitter(const TString& configFileName)
166166
const Value& rebinValue = config["Rebin"];
167167
readArray(rebinValue, nRebin);
168168

169-
bool includeSecPeak = config["InclSecPeak"].GetBool();
170-
bool useLikelihood = config["UseLikelihood"].GetBool();
169+
bool const includeSecPeak = config["InclSecPeak"].GetBool();
170+
bool const useLikelihood = config["UseLikelihood"].GetBool();
171171

172172
const Value& bkgFuncValue = config["BkgFunc"];
173173
readArray(bkgFuncValue, bkgFuncConfig);
@@ -349,7 +349,7 @@ int runMassFitter(const TString& configFileName)
349349
canvasSize[1] = 500;
350350
}
351351

352-
Int_t nCanvasesMax = 20; // do not put more than 20 bins per canvas to make them visible
352+
Int_t const nCanvasesMax = 20; // do not put more than 20 bins per canvas to make them visible
353353
const Int_t nCanvases = std::ceil(static_cast<float>(nSliceVarBins) / nCanvasesMax);
354354
std::vector<TCanvas*> canvasMass(nCanvases);
355355
std::vector<TCanvas*> canvasResiduals(nCanvases);
@@ -372,7 +372,7 @@ int runMassFitter(const TString& configFileName)
372372
const Int_t iCanvas = std::floor(static_cast<float>(iSliceVar) / nCanvasesMax);
373373

374374
hMassForFit[iSliceVar] = hMass[iSliceVar]->Rebin(nRebin[iSliceVar]);
375-
TString ptTitle =
375+
TString const ptTitle =
376376
Form("%0.2f < " + sliceVarName + " < %0.2f " + sliceVarUnit, sliceVarMin[iSliceVar], sliceVarMax[iSliceVar]);
377377
hMassForFit[iSliceVar]->SetTitle(Form("%s;%s;Counts per %0.1f MeV/#it{c}^{2}",
378378
ptTitle.Data(), massAxisTitle.Data(),

‎PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx‎

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -352,8 +352,8 @@ struct HfCandidateCreatorB0Reduced {
352352
}
353353
// invMassWindowDPiTolerance is used to apply a slightly tighter cut than in DPi pair preselection
354354
// to avoid accepting DPi pairs that were not formed in DPi pair creator
355-
float invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
356-
float invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
355+
float const invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
356+
float const invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
357357

358358
for (const auto& collisionCounter : collisionsCounter) {
359359
registry.fill(HIST("hEvents"), 1, collisionCounter.originalCollisionCount());
@@ -386,8 +386,8 @@ struct HfCandidateCreatorB0Reduced {
386386
}
387387
// invMassWindowDPiTolerance is used to apply a slightly tighter cut than in DPi pair preselection
388388
// to avoid accepting DPi pairs that were not formed in DPi pair creator
389-
float invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
390-
float invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
389+
float const invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
390+
float const invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
391391

392392
for (const auto& collisionCounter : collisionsCounter) {
393393
registry.fill(HIST("hEvents"), 1, collisionCounter.originalCollisionCount());
@@ -421,8 +421,8 @@ struct HfCandidateCreatorB0Reduced {
421421
}
422422
// invMassWindowDPiTolerance is used to apply a slightly tighter cut than in DPi pair preselection
423423
// to avoid accepting DPi pairs that were not formed in DPi pair creator
424-
float invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
425-
float invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
424+
float const invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
425+
float const invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
426426

427427
for (const auto& collisionCounter : collisionsCounter) {
428428
registry.fill(HIST("hEvents"), 1, collisionCounter.originalCollisionCount());
@@ -456,8 +456,8 @@ struct HfCandidateCreatorB0Reduced {
456456
}
457457
// invMassWindowDPiTolerance is used to apply a slightly tighter cut than in DPi pair preselection
458458
// to avoid accepting DPi pairs that were not formed in DPi pair creator
459-
float invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
460-
float invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
459+
float const invMass2DPiMin = (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 - myInvMassWindowDPi + invMassWindowDPiTolerance);
460+
float const invMass2DPiMax = (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance) * (o2::constants::physics::MassB0 + myInvMassWindowDPi - invMassWindowDPiTolerance);
461461

462462
for (const auto& collisionCounter : collisionsCounter) {
463463
registry.fill(HIST("hEvents"), 1, collisionCounter.originalCollisionCount());

0 commit comments

Comments
 (0)