Skip to content

Commit 8863b2d

Browse files
committed
Merge remote-tracking branch 'upstream/master' into fixes-PWGUD
2 parents 6fac105 + 33bb1c0 commit 8863b2d

7 files changed

Lines changed: 75 additions & 47 deletions

File tree

PWGCF/EbyEFluctuations/Tasks/radialFlowDecorr.cxx

Lines changed: 58 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -76,7 +76,7 @@ struct RadialFlowDecorr {
7676
static constexpr int KnFt0cCell = 96;
7777
static constexpr int KIntM = 3;
7878
static constexpr int KIntK = 3;
79-
static constexpr int KNEta = 17;
79+
static constexpr int KNEta = 9;
8080
static constexpr float KFloatEpsilon = 1e-6f;
8181
static constexpr int KPiPlus = 211;
8282
static constexpr int KKPlus = 321;
@@ -151,10 +151,10 @@ struct RadialFlowDecorr {
151151
static constexpr float KinvalidCentrality = -1.0f;
152152
inline static const std::vector<float> etaLw = {
153153
-0.8,
154-
-0.8, -0.7, -0.6, -0.5, -0.4, -0.3, -0.2, -0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7};
154+
-0.8, -0.6, -0.4, -0.2, 0.0, 0.2, 0.4, 0.6};
155155
inline static const std::vector<float> etaUp = {
156156
0.8,
157-
-0.7, -0.6, -0.5, -0.4, -0.3, -0.2, -0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8};
157+
-0.6, -0.4, -0.2, 0.0, 0.2, 0.4, 0.6, 0.8};
158158

159159
Configurable<float> cfgVtxZCut{"cfgVtxZCut", 10.f, "z-vertex range"};
160160
Configurable<float> cfgPtMin{"cfgPtMin", 0.2f, "min pT"};
@@ -187,6 +187,7 @@ struct RadialFlowDecorr {
187187
Configurable<float> cfgCutPtUpper{"cfgCutPtUpper", 10.0f, "Higher pT cut for inclusive hadron analysis"};
188188
Configurable<float> cfgCutPtUpperPID{"cfgCutPtUpperPID", 6.0f, "Higher pT cut for identified particle analysis"};
189189
Configurable<float> cfgCutEta{"cfgCutEta", 0.8f, "absolute Eta cut"};
190+
Configurable<int> cfgMinTracksPerEtaBin{"cfgMinTracksPerEtaBin", 0, "Min weighted-track sum required in every narrow eta bin for inclusive species (0 = disabled)"};
190191
Configurable<int> cfgNsubsample{"cfgNsubsample", 10, "Number of subsamples"};
191192
Configurable<int> cfgCentralityChoice{"cfgCentralityChoice", 1, "Which centrality estimator? 1-->FT0C, 2-->FT0M, 3-->FDDM, 4-->FV0A"};
192193
Configurable<bool> cfgEvSelNoSameBunchPileup{"cfgEvSelNoSameBunchPileup", true, "Pileup removal"};
@@ -201,7 +202,7 @@ struct RadialFlowDecorr {
201202
Configurable<float> cfgLinPupParam2{"cfgLinPupParam2", 3.0f, "(Lower) Linear Pileup Cut Const"};
202203
Configurable<float> cfgLinPupParam3{"cfgLinPupParam3", 3.0f, "(Lower) Linear Pileup Slope"};
203204

204-
Configurable<int> cfgNchPbMax{"cfgNchPbMax", 4000, "Max Nch range for PbPb collisions"};
205+
Configurable<int> cfgNchPbMax{"cfgNchPbMax", 5000, "Max Nch range for PbPb collisions"};
205206
Configurable<int> cfgNchOMax{"cfgNchOMax", 800, "Max Nch range for OO collisions"};
206207

207208
Configurable<int> cfgSys{"cfgSys", 1, "Efficiency to be used for which system? 1-->PbPb, 2-->NeNe, 3-->OO, 4-->pp"};
@@ -222,21 +223,17 @@ struct RadialFlowDecorr {
222223
const AxisSpec vzAxis{5, -12.5, 12.5, "Vz"};
223224
const AxisSpec chgAxis{3, -1.5, 1.5};
224225
const AxisSpec pTAxis{{0.0, 0.2, 0.4, 0.6, 0.8, 1, 3, 5, 7, 10}, "pT Axis"};
225-
const AxisSpec etaAxis{{-0.9, -0.8, -0.7, -0.6, -0.5, -0.4, -0.3, -0.2, -0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9}, "Eta"};
226+
const AxisSpec etaAxis{{-0.8, -0.6, -0.4, -0.2, 0.0, 0.2, 0.4, 0.6, 0.8}, "Eta"};
226227
const AxisSpec phiAxis{KNbinsPhi, KPhiMin, TwoPI, "#phi"};
227228
const AxisSpec etaBinAxis{KNEta + 1, -0.5, KNEta + 0.5, "#eta bin Number"};
228229
const AxisSpec spBinAxis{KNsp + 1, -KBinOffset, static_cast<float>(KNsp) + KBinOffset, "species index Number"};
229230

230-
const AxisSpec gapAxis{{-1.55, -1.45, -1.35, -1.25, -1.15, -1.05, -0.95, -0.85,
231-
-0.75, -0.65, -0.55, -0.45, -0.35, -0.25, -0.15, -0.05,
232-
0.05, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75,
233-
0.85, 0.95, 1.05, 1.15, 1.25, 1.35, 1.45, 1.55},
231+
const AxisSpec gapAxis{{-1.5, -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1,
232+
0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5},
234233
"Gap"};
235234

236-
const AxisSpec sumAxis{{-1.55, -1.45, -1.35, -1.25, -1.15, -1.05, -0.95, -0.85,
237-
-0.75, -0.65, -0.55, -0.45, -0.35, -0.25, -0.15, -0.05,
238-
0.05, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75,
239-
0.85, 0.95, 1.05, 1.15, 1.25, 1.35, 1.45, 1.55},
235+
const AxisSpec sumAxis{{-1.5, -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1,
236+
0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5},
240237
"Sum"};
241238

242239
Configurable<bool> cfgRunMCGetNSig{"cfgRunMCGetNSig", false, "Run MC pass to get mean of Nsig Plots"};
@@ -552,6 +549,32 @@ struct RadialFlowDecorr {
552549
return true;
553550
}
554551

552+
template <std::size_t NspT, std::size_t NetaT, std::size_t NkT>
553+
bool hasMinTracksInAllEtaBins(const double (&sw)[NspT][NetaT][NkT])
554+
{
555+
const int minTracks = cfgMinTracksPerEtaBin;
556+
if (minTracks <= 0)
557+
return true;
558+
for (std::size_t ieta = 1; ieta < NetaT; ++ieta) {
559+
if (sw[kInclusiveIdx][ieta][1] < static_cast<double>(minTracks))
560+
return false;
561+
}
562+
return true;
563+
}
564+
565+
template <std::size_t NspT, std::size_t NetaT>
566+
bool hasMinTracksInAllEtaBins(const double (&sw)[NspT][NetaT])
567+
{
568+
const int minTracks = cfgMinTracksPerEtaBin;
569+
if (minTracks <= 0)
570+
return true;
571+
for (std::size_t ieta = 1; ieta < NetaT; ++ieta) {
572+
if (sw[kInclusiveIdx][ieta] < static_cast<double>(minTracks))
573+
return false;
574+
}
575+
return true;
576+
}
577+
555578
template <typename T>
556579
bool isTrackSelected(const T& trk)
557580
{
@@ -1996,6 +2019,9 @@ struct RadialFlowDecorr {
19962019
}
19972020
}
19982021

2022+
if (!hasMinTracksInAllEtaBins(sumWiTruth) || !hasMinTracksInAllEtaBins(sumWiReco))
2023+
return;
2024+
19992025
for (int isp = 0; isp < KNsp; ++isp) {
20002026
histos.fill(HIST("MCReco/Prof_Cent_Nsp_Nchrec"), cent, isp, sumWiReco[isp][0]);
20012027
histos.fill(HIST("MCReco/Prof_Mult_Nsp_Nchrec"), multPV, isp, sumWiReco[isp][0]);
@@ -2352,6 +2378,9 @@ struct RadialFlowDecorr {
23522378
}
23532379
} // trkslice
23542380

2381+
if (!hasMinTracksInAllEtaBins(sumWkTru) || !hasMinTracksInAllEtaBins(sumWkReco))
2382+
return;
2383+
23552384
for (int ieta = 0; ieta < KNEta; ++ieta) {
23562385
const int ibx = state.pmeanTruNchEtabinSpbinStep2->GetXaxis()->FindBin(mcCollision.multNTracksPV());
23572386
const int iby = ieta + 1;
@@ -2546,9 +2575,12 @@ struct RadialFlowDecorr {
25462575
float sum = (etaValA + etaValB);
25472576
for (int isp = 0; isp < KNsp; ++isp) {
25482577

2549-
float c2SubTru = p1kBarTru[isp][ietaA] * p1kBarTru[isp][ietaC];
2550-
float c2SubReco = p1kBarReco[isp][ietaA] * p1kBarReco[isp][ietaC];
2551-
float c2SubRecoEffCor = p1kBarRecoEffCor[isp][ietaA] * p1kBarRecoEffCor[isp][ietaC];
2578+
float c2SubTru = (ietaA == ietaC) ? static_cast<float>(c2Tru[isp][ietaA])
2579+
: p1kBarTru[isp][ietaA] * p1kBarTru[isp][ietaC];
2580+
float c2SubReco = (ietaA == ietaC) ? static_cast<float>(c2Reco[isp][ietaA])
2581+
: p1kBarReco[isp][ietaA] * p1kBarReco[isp][ietaC];
2582+
float c2SubRecoEffCor = (ietaA == ietaC) ? static_cast<float>(c2RecoEffCor[isp][ietaA])
2583+
: p1kBarRecoEffCor[isp][ietaA] * p1kBarRecoEffCor[isp][ietaC];
25522584

25532585
float covTru = p1kBarTruMult[isp][ietaA] * p1kBarTru[isp][ietaC];
25542586
float covReco = p1kBarRecoMult[isp][ietaA] * p1kBarReco[isp][ietaC];
@@ -3200,6 +3232,9 @@ struct RadialFlowDecorr {
32003232
}
32013233
}
32023234

3235+
if (!hasMinTracksInAllEtaBins(sumWi))
3236+
return;
3237+
32033238
for (int isp = 0; isp < KNsp; ++isp) {
32043239
if (sumWi[isp][0] < 1.0f)
32053240
continue;
@@ -3373,6 +3408,9 @@ struct RadialFlowDecorr {
33733408
}
33743409
}
33753410

3411+
if (!hasMinTracksInAllEtaBins(sumwk))
3412+
return;
3413+
33763414
double amplFT0A = 0, amplFT0C = 0;
33773415
if (coll.has_foundFT0()) {
33783416
const auto& ft0 = coll.foundFT0();
@@ -3440,8 +3478,8 @@ struct RadialFlowDecorr {
34403478
histos.fill(HIST("Prof_Cov_Mult_etabin_spbin"), coll.multNTracksPV(), ietaA, isp, covAC);
34413479
}
34423480
if (std::isfinite(covCA)) {
3443-
histos.fill(HIST("Prof_Cov_Cent_etabin_spbin"), cent, ietaA, isp, covCA);
3444-
histos.fill(HIST("Prof_Cov_Mult_etabin_spbin"), coll.multNTracksPV(), ietaA, isp, covCA);
3481+
histos.fill(HIST("Prof_Cov_Cent_etabin_spbin"), cent, ietaC, isp, covCA);
3482+
histos.fill(HIST("Prof_Cov_Mult_etabin_spbin"), coll.multNTracksPV(), ietaC, isp, covCA);
34453483
}
34463484
if (std::isfinite(covFT0A)) {
34473485
histos.fill(HIST("Prof_CovFT0A_Cent_etabin_spbin"), cent, ietaA, isp, covFT0A);
@@ -3464,7 +3502,8 @@ struct RadialFlowDecorr {
34643502

34653503
for (int isp = 0; isp < KNsp; ++isp) {
34663504

3467-
float c2Sub = p1kBar[isp][ietaA] * p1kBar[isp][ietaC];
3505+
float c2Sub = (ietaA == ietaC) ? static_cast<float>(c2[isp][ietaA])
3506+
: p1kBar[isp][ietaA] * p1kBar[isp][ietaC];
34683507
float cov = p1kBarMult[isp][ietaA] * p1kBar[isp][ietaC];
34693508
float covFT0A = p1kBarFt0A * p1kBar[isp][ietaC];
34703509
float covFT0C = p1kBarFt0C * p1kBar[isp][ietaA];

PWGHF/D2H/Macros/runMassFitter.C

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -38,6 +38,7 @@
3838
#include <RtypesCore.h>
3939

4040
#include <algorithm>
41+
#include <array>
4142
#include <cmath>
4243
#include <cstdio>
4344
#include <functional>

PWGHF/D2H/Tasks/taskFlowCharmHadrons.cxx

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -52,6 +52,7 @@
5252
#include <cmath>
5353
#include <cstdint>
5454
#include <cstdlib>
55+
#include <numeric>
5556
#include <string>
5657
#include <vector>
5758

PWGHF/D2H/Tasks/taskLc.cxx

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -32,7 +32,6 @@
3232
#include "Common/Core/RecoDecay.h"
3333
#include "Common/DataModel/Centrality.h"
3434
#include "Common/DataModel/EventSelection.h"
35-
#include "Common/DataModel/Multiplicity.h"
3635

3736
#include <CCDB/BasicCCDBManager.h>
3837
#include <CommonConstants/PhysicsConstants.h>
@@ -45,7 +44,6 @@
4544
#include <Framework/HistogramSpec.h>
4645
#include <Framework/InitContext.h>
4746
#include <Framework/Logger.h>
48-
#include <Framework/OutputObjHeader.h>
4947
#include <Framework/runDataProcessing.h>
5048

5149
#include <THnSparse.h>

PWGHF/D2H/Tasks/taskUpcLc.cxx

Lines changed: 11 additions & 26 deletions
Original file line numberDiff line numberDiff line change
@@ -17,7 +17,6 @@
1717
/// \author Ran Tu <ran.tu@cern.ch>, Fudan University
1818

1919
#include "PWGHF/Core/CentralityEstimation.h"
20-
#include "PWGHF/Core/DecayChannels.h"
2120
#include "PWGHF/Core/HfHelper.h"
2221
#include "PWGHF/Core/SelectorCuts.h"
2322
#include "PWGHF/DataModel/AliasTables.h"
@@ -29,12 +28,11 @@
2928
#include "PWGUD/Core/SGSelector.h"
3029
#include "PWGUD/Core/UPCHelpers.h"
3130

32-
#include "Common/Core/RecoDecay.h"
3331
#include "Common/DataModel/EventSelection.h"
3432
#include "Common/DataModel/Multiplicity.h"
3533

3634
#include <CCDB/BasicCCDBManager.h>
37-
#include <CommonConstants/PhysicsConstants.h>
35+
#include <CommonDataFormat/TimeStamp.h>
3836
#include <Framework/ASoA.h>
3937
#include <Framework/AnalysisDataModel.h>
4038
#include <Framework/AnalysisHelpers.h>
@@ -44,12 +42,9 @@
4442
#include <Framework/HistogramSpec.h>
4543
#include <Framework/InitContext.h>
4644
#include <Framework/Logger.h>
47-
#include <Framework/OutputObjHeader.h>
4845
#include <Framework/runDataProcessing.h>
4946
#include <ReconstructionDataFormats/Vertex.h>
5047

51-
#include <TPDGCode.h>
52-
5348
#include <array>
5449
#include <cmath>
5550
#include <numeric>
@@ -71,18 +66,15 @@ namespace full
7166
DECLARE_SOA_COLUMN(M, m, float);
7267
DECLARE_SOA_COLUMN(Pt, pt, float);
7368
DECLARE_SOA_COLUMN(BkgScore, bkgScore, float);
74-
DECLARE_SOA_COLUMN(PromptScore, promptScore, float);
75-
DECLARE_SOA_COLUMN(FdScore, fdScore, float);
7669
DECLARE_SOA_COLUMN(PtProng0, ptProng0, float);
7770
DECLARE_SOA_COLUMN(PtProng1, ptProng1, float);
7871
DECLARE_SOA_COLUMN(PtProng2, ptProng2, float);
7972
DECLARE_SOA_COLUMN(Chi2PCA, chi2PCA, float);
8073
DECLARE_SOA_COLUMN(DecayLength, decayLength, float);
8174
DECLARE_SOA_COLUMN(Cpa, cpa, float);
82-
DECLARE_SOA_COLUMN(PvContributors, pvContributors, float);
83-
DECLARE_SOA_COLUMN(Multiplicity, multiplicity, float);
75+
DECLARE_SOA_COLUMN(PvContributors, pvContributors, int);
76+
DECLARE_SOA_COLUMN(Multiplicity, multiplicity, int);
8477
DECLARE_SOA_COLUMN(Vtz, vtz, float);
85-
DECLARE_SOA_COLUMN(AmpFV0A, ampFV0A, float);
8678
DECLARE_SOA_COLUMN(AmpFT0A, ampFT0A, float);
8779
DECLARE_SOA_COLUMN(AmpFT0C, ampFT0C, float);
8880
DECLARE_SOA_COLUMN(ZdcTimeZNA, zdcTimeZNA, float);
@@ -92,7 +84,6 @@ DECLARE_SOA_TABLE(HfUpcQa, "AOD", "HFUPCQA",
9284
full::PvContributors,
9385
full::Multiplicity,
9486
full::Vtz,
95-
full::AmpFV0A,
9687
full::AmpFT0A,
9788
full::AmpFT0C,
9889
full::ZdcTimeZNA,
@@ -102,9 +93,6 @@ DECLARE_SOA_TABLE(HfUpcLcBdtInfos, "AOD", "HFUPCLCBDTINFOS",
10293
full::M,
10394
full::Pt,
10495
full::BkgScore,
105-
full::PromptScore,
106-
full::FdScore,
107-
full::AmpFV0A,
10896
full::AmpFT0A,
10997
full::AmpFT0C,
11098
full::ZdcTimeZNA,
@@ -119,7 +107,6 @@ DECLARE_SOA_TABLE(HfUpcLcInfos, "AOD", "HFUPCLCINFOS",
119107
full::Chi2PCA,
120108
full::DecayLength,
121109
full::Cpa,
122-
full::AmpFV0A,
123110
full::AmpFT0A,
124111
full::AmpFT0C,
125112
full::ZdcTimeZNA,
@@ -256,9 +243,9 @@ struct HfTaskUpcLc {
256243
if (gap == o2::aod::sgselector::TrueGap::SingleGapA || gap == o2::aod::sgselector::TrueGap::SingleGapC) {
257244
registry.fill(HIST("Data/hUpcMulti"), collision.multNTracksPV());
258245
registry.fill(HIST("Data/hUpcVtz"), collision.posZ());
259-
}
260-
if (fillTreeUpcQa) {
261-
rowUpcQa(numPvContributors, collision.multNTracksPV(), collision.posZ(), fitInfo.ampFV0A, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
246+
if (fillTreeUpcQa) {
247+
rowUpcQa(numPvContributors, collision.multNTracksPV(), collision.posZ(), fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
248+
}
262249
}
263250

264251
for (const auto& candidate : groupedLcCandidates) {
@@ -276,7 +263,7 @@ struct HfTaskUpcLc {
276263
const auto chi2PCA = candidate.chi2PCA();
277264
const auto cpa = candidate.cpa();
278265

279-
double outputBkg(-1), outputPrompt(-1), outputFD(-1);
266+
double outputBkg(-1);
280267

281268
auto fillTHnData = [&](bool isPKPi) {
282269
const auto massLc = isPKPi ? HfHelper::invMassLcToPKPi(candidate) : HfHelper::invMassLcToPiKP(candidate);
@@ -285,25 +272,23 @@ struct HfTaskUpcLc {
285272
const auto& mlProb = isPKPi ? candidate.mlProbLcToPKPi() : candidate.mlProbLcToPiKP();
286273
if (mlProb.size() == NumberOfMlClasses) {
287274
outputBkg = mlProb[MlClassBackground]; /// bkg score
288-
outputPrompt = mlProb[MlClassPrompt]; /// prompt score
289-
outputFD = mlProb[MlClassNonPrompt]; /// non-prompt score
290275
}
291276
/// Fill the ML outputScores and variables of candidate
292277
if (fillTreeOnlySingleGap) {
293278
if (gap == o2::aod::sgselector::TrueGap::SingleGapA || gap == o2::aod::sgselector::TrueGap::SingleGapC) {
294-
rowCandUpcBdt(massLc, pt, outputBkg, outputPrompt, outputFD, fitInfo.ampFV0A, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
279+
rowCandUpcBdt(massLc, pt, outputBkg, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
295280
}
296281
} else {
297-
rowCandUpcBdt(massLc, pt, outputBkg, outputPrompt, outputFD, fitInfo.ampFV0A, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
282+
rowCandUpcBdt(massLc, pt, outputBkg, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
298283
}
299284

300285
} else {
301286
if (fillTreeOnlySingleGap) {
302287
if (gap == o2::aod::sgselector::TrueGap::SingleGapA || gap == o2::aod::sgselector::TrueGap::SingleGapC) {
303-
rowCandUpc(massLc, pt, ptProng0, ptProng1, ptProng2, chi2PCA, decayLength, cpa, fitInfo.ampFV0A, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
288+
rowCandUpc(massLc, pt, ptProng0, ptProng1, ptProng2, chi2PCA, decayLength, cpa, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
304289
}
305290
} else {
306-
rowCandUpc(massLc, pt, ptProng0, ptProng1, ptProng2, chi2PCA, decayLength, cpa, fitInfo.ampFV0A, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
291+
rowCandUpc(massLc, pt, ptProng0, ptProng1, ptProng2, chi2PCA, decayLength, cpa, fitInfo.ampFT0A, fitInfo.ampFT0C, zdcTimeZNA, zdcTimeZNC);
307292
}
308293
}
309294
};

PWGHF/HFC/TableProducer/correlatorLcScHadrons.cxx

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -54,6 +54,8 @@
5454
#include <Framework/runDataProcessing.h>
5555
#include <ReconstructionDataFormats/PID.h>
5656

57+
#include <Math/Vector4D.h> // IWYU pragma: keep (do not replace with Math/Vector4Dfwd.h)
58+
#include <Math/Vector4Dfwd.h>
5759
#include <TPDGCode.h>
5860
#include <TRandom3.h>
5961

PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
#include "PWGHF/DataModel/TrackIndexSkimmingTables.h"
2929
#include "PWGHF/Utils/utilsBfieldCCDB.h"
3030
#include "PWGHF/Utils/utilsEvSelHf.h"
31+
//
3132
#include "PWGLF/DataModel/LFStrangenessTables.h"
3233
#include "PWGLF/DataModel/mcCentrality.h"
3334

@@ -57,6 +58,7 @@
5758
#include <Framework/RunningWorkflowInfo.h>
5859
#include <Framework/runDataProcessing.h>
5960
#include <ReconstructionDataFormats/DCA.h>
61+
#include <ReconstructionDataFormats/PID.h>
6062
#include <ReconstructionDataFormats/Track.h>
6163

6264
#include <TH1.h>

0 commit comments

Comments
 (0)