From 9ce92df06c7358b192226b319fda68475fbf124b Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Mon, 13 Jul 2026 17:04:35 +0200 Subject: [PATCH 1/9] MC HMPID table changes --- .../AOTTrack/PID/HMPID/hmpidTableProducer.cxx | 66 +++++++++++++------ DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h | 4 +- 2 files changed, 49 insertions(+), 21 deletions(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx index 0e7c1e3b464..f7243e1a5bd 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx +++ b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx @@ -24,13 +24,8 @@ #include #include #include -#include #include -#include #include -#include -#include -#include #include #include @@ -72,6 +67,10 @@ struct HmpidTableProducer { Configurable requireTPC{"requireTPC", true, "Require TPC track"}; Configurable requireTOF{"requireTOF", true, "Require TOF track"}; + // Absorbers position + Configurable absorberEdgeRich2{"absorberEdgeRich2", 435.5f, "Raggio esterno (bordo) del target Al davanti a rich2 [cm]"}; + Configurable absorberEdgeRich4{"absorberEdgeRich4", 435.0f, "Raggio esterno (bordo) del target Al davanti a rich4 [cm]"}; + using CollisionCandidates = o2::soa::Join; using TrackCandidates = soa::Join mCollisionsWithHmpid; + const int rich2 = 2, rich4 = 4; + const float kRich2AbsorberThickness = 4.0f; + const float kRich4AbsorberThickness = 8.0f; + void init(o2::framework::InitContext&) { ccdb->setURL(ccdbConfig.ccdbUrl); @@ -305,14 +308,15 @@ struct HmpidTableProducer { } PROCESS_SWITCH(HmpidTableProducer, processEvent, "Process event level", true); - template + template void runHmpidAnalysis( aod::HMPIDs const& hmpids, TTrackTable const&, CollisionCandidates const&, - aod::BCsWithTimestamps const&) + aod::BCsWithTimestamps const&, + TMcParticles const& mcParticles = 0) { - for (auto const& t : hmpids) { + for (auto const& t : hmpids) { // begin loop over hmpids const auto& globalTrack = t.template track_as(); @@ -366,17 +370,17 @@ struct HmpidTableProducer { histos.fill(HIST("hChamberM1"), chamberM1); - if (!isCorrupt) { + if (!isCorrupt) { // begin if(!isCorrupt) - fill M1vsM2 histos.fill(HIST("hChamberM1vsM2"), chamberM2, chamberM1); - } + } // end if(!isCorrupt) - fill M1vsM2 // --- M3: hybrid --- int chamberM3 = -1; - if (!isCorrupt) { + if (!isCorrupt) { // begin if/else - M3 assignment chamberM3 = chamberM2; } else { chamberM3 = chamberM1; - } + } // end if/else - M3 assignment // Legend: // bin 0 = clusSize > 0, chamber found (M2 ok) @@ -400,8 +404,9 @@ struct HmpidTableProducer { continue; float hmpidPhotsCharge2[o2::aod::kDimPhotonsCharge]; - for (int i = 0; i < o2::aod::kDimPhotonsCharge; i++) + for (int i = 0; i < o2::aod::kDimPhotonsCharge; i++) // begin for - copy photon charges hmpidPhotsCharge2[i] = t.hmpidPhotsCharge()[i]; + // end for - copy photon charges // fill hmpid table hmpidAnalysis( @@ -427,14 +432,35 @@ struct HmpidTableProducer { if constexpr (isMC) { if (globalTrack.has_mcParticle()) { const auto& mc = globalTrack.mcParticle(); - hmpidAnalysisMC(mc.pdgCode(), mc.vx(), mc.vy(), mc.vz(), mc.isPhysicalPrimary(), mc.getProcess()); + + bool interactionInAbsorber = false; + + if (mc.has_daughters()) { + auto dIds = mc.daughtersIds(); + for (int32_t idx = dIds.front(); idx <= dIds.back(); ++idx) { + auto daughter = mcParticles.rawIteratorAt(idx); + double r = std::hypot(daughter.vx(), daughter.vy()); + + if (chamberM3 == rich2 && r >= absorberEdgeRich2 - kRich2AbsorberThickness && r <= absorberEdgeRich2) { + interactionInAbsorber = true; + break; + } + if (chamberM3 == rich4 && r >= absorberEdgeRich4 - kRich4AbsorberThickness && r <= absorberEdgeRich4) { + interactionInAbsorber = true; + break; + } + } // end loop daughters + } // end if has_daughters + + hmpidAnalysisMC(mc.pdgCode(), mc.vx(), mc.vy(), mc.vz(), + mc.isPhysicalPrimary(), mc.getProcess(), interactionInAbsorber); } else { - hmpidAnalysisMC(-1, 0.f, 0.f, 0.f, false, -100); + hmpidAnalysisMC(-1, 0.f, 0.f, 0.f, false, -100, false); } - } + } // end if constexpr (isMC) - } // end HMPID loop - } + } // end for - loop over hmpids + } // end runHmpidAnalysis // process real data void processHmpid(aod::HMPIDs const& hmpids, @@ -452,9 +478,9 @@ struct HmpidTableProducer { TrackCandidatesMC const& tracks, CollisionCandidates const& cols, aod::BCsWithTimestamps const& bcs, - aod::McParticles const&) + aod::McParticles const& mcParticles) { - runHmpidAnalysis(hmpids, tracks, cols, bcs); + runHmpidAnalysis(hmpids, tracks, cols, bcs, mcParticles); } PROCESS_SWITCH(HmpidTableProducer, processHmpidMC, "Process HMPID MC entries", true); }; diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h index 1fcf9980763..f3f6a93f843 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h +++ b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h @@ -103,6 +103,7 @@ DECLARE_SOA_COLUMN(McVy, mcVy, float); DECLARE_SOA_COLUMN(McVz, mcVz, float); DECLARE_SOA_COLUMN(IsPhysPrimary, isPhysPrimary, bool); DECLARE_SOA_COLUMN(ProcessCode, processCode, int); +DECLARE_SOA_COLUMN(HasInteractedInAbsorber, hasInteractedInAbsorber, bool); } // namespace hmpid_mc @@ -112,7 +113,8 @@ DECLARE_SOA_TABLE(HmpidAnalysisMC, "AOD", "HMPIDANALYSISMC", hmpid_mc::McVy, hmpid_mc::McVz, hmpid_mc::IsPhysPrimary, - hmpid_mc::ProcessCode); + hmpid_mc::ProcessCode, + hmpid_mc::HasInteractedInAbsorber); } // namespace o2::aod From e74e8bae60374d7907d998534c53dfbb4b4ff305 Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Mon, 13 Jul 2026 17:18:11 +0200 Subject: [PATCH 2/9] Fixing MegaLinter --- DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx index f7243e1a5bd..c4bf56ab7e5 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx +++ b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx @@ -252,11 +252,14 @@ struct HmpidTableProducer { } int getHmpidChamber( - std::array x, - std::array p, + std::array xIn, + std::array pIn, double bz, int charge) { + auto x = xIn; + auto p = pIn; + auto* param = o2::hmpid::Param::instance(); for (int ch = o2::hmpid::Param::kMinCh; ch <= o2::hmpid::Param::kMaxCh; ++ch) { From e1ea7b25f4f884b1aa934bc0d07adc12b517b900 Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Mon, 13 Jul 2026 18:27:07 +0200 Subject: [PATCH 3/9] Fixing build --- .../AOTTrack/PID/HMPID/hmpidTableProducer.cxx | 40 ++++++++++--------- 1 file changed, 21 insertions(+), 19 deletions(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx index c4bf56ab7e5..8b42193b7ae 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx +++ b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx @@ -50,7 +50,7 @@ struct HmpidTableProducer { HistogramRegistry histos{"histos", {}, OutputObjHandlingPolicy::AnalysisObject}; - const AxisSpec axisEvtCounter{1, 0, +1, ""}; + AxisSpec axisEvtCounter{1, 0, +1, ""}; Service ccdb; struct : ConfigurableGroup { @@ -68,8 +68,8 @@ struct HmpidTableProducer { Configurable requireTOF{"requireTOF", true, "Require TOF track"}; // Absorbers position - Configurable absorberEdgeRich2{"absorberEdgeRich2", 435.5f, "Raggio esterno (bordo) del target Al davanti a rich2 [cm]"}; - Configurable absorberEdgeRich4{"absorberEdgeRich4", 435.0f, "Raggio esterno (bordo) del target Al davanti a rich4 [cm]"}; + Configurable absorberEdgeRich2{"absorberEdgeRich2", 435.5f, "Raggio esterno (bordo) del target Al davanti a Rich2 [cm]"}; + Configurable absorberEdgeRich4{"absorberEdgeRich4", 435.0f, "Raggio esterno (bordo) del target Al davanti a Rich4 [cm]"}; using CollisionCandidates = o2::soa::Join; @@ -90,9 +90,9 @@ struct HmpidTableProducer { std::unordered_set mCollisionsWithHmpid; - const int rich2 = 2, rich4 = 4; - const float kRich2AbsorberThickness = 4.0f; - const float kRich4AbsorberThickness = 8.0f; + static constexpr int Rich2 = 2, Rich4 = 4; + static constexpr float Rich2AbsorberThickness = 4.0f; + static constexpr float Rich4AbsorberThickness = 8.0f; void init(o2::framework::InitContext&) { @@ -252,13 +252,16 @@ struct HmpidTableProducer { } int getHmpidChamber( - std::array xIn, - std::array pIn, + const std::array& xIn, + const std::array& pIn, double bz, int charge) { - auto x = xIn; - auto p = pIn; + std::array x{}; + std::array p{}; + + x = xIn; + p = pIn; auto* param = o2::hmpid::Param::instance(); @@ -274,21 +277,21 @@ struct HmpidTableProducer { param->norm(ch, nPc.data()); // Intersection track - radiator plane - std::array xRad, pAtRad; + std::array xRad{}, pAtRad{}; if (!intersectHelixPlane(bz, charge, x, p, pRad, nRad, xRad, pAtRad)) continue; // Intersection track - PC plane - std::array xPc, pAtPc; + std::array xPc{}, pAtPc{}; if (!intersectHelixPlane(bz, charge, xRad, pAtRad, pPc, nPc, xPc, pAtPc)) continue; - double theta, phi; + double theta = 0., phi = 0.; param->mars2LorsVec(ch, pAtRad.data(), theta, phi); - double xL, yL; + double xL = 0., yL = 0.; param->mars2Lors(ch, xPc.data(), xL, yL); // Use isInside to check Chamber intersected @@ -367,7 +370,7 @@ struct HmpidTableProducer { int16_t charge = globalTrack.sign(); auto prop = o2::base::Propagator::Instance(); - double bz = static_cast(prop->getNominalBz()); + auto bz = static_cast(prop->getNominalBz()); int chamberM1 = getHmpidChamber(x, p, bz, charge); @@ -407,9 +410,8 @@ struct HmpidTableProducer { continue; float hmpidPhotsCharge2[o2::aod::kDimPhotonsCharge]; - for (int i = 0; i < o2::aod::kDimPhotonsCharge; i++) // begin for - copy photon charges + for (int i = 0; i < o2::aod::kDimPhotonsCharge; i++) hmpidPhotsCharge2[i] = t.hmpidPhotsCharge()[i]; - // end for - copy photon charges // fill hmpid table hmpidAnalysis( @@ -444,11 +446,11 @@ struct HmpidTableProducer { auto daughter = mcParticles.rawIteratorAt(idx); double r = std::hypot(daughter.vx(), daughter.vy()); - if (chamberM3 == rich2 && r >= absorberEdgeRich2 - kRich2AbsorberThickness && r <= absorberEdgeRich2) { + if (chamberM3 == Rich2 && r >= absorberEdgeRich2 - Rich2AbsorberThickness && r <= absorberEdgeRich2) { interactionInAbsorber = true; break; } - if (chamberM3 == rich4 && r >= absorberEdgeRich4 - kRich4AbsorberThickness && r <= absorberEdgeRich4) { + if (chamberM3 == Rich4 && r >= absorberEdgeRich4 - Rich4AbsorberThickness && r <= absorberEdgeRich4) { interactionInAbsorber = true; break; } From 7ce6db031bdb6ddaf11c20df2180a2e07c5991e8 Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Mon, 13 Jul 2026 19:33:11 +0200 Subject: [PATCH 4/9] Fixing includes --- DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx index 8b42193b7ae..eb089da60c4 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx +++ b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx @@ -24,8 +24,13 @@ #include #include #include +#include #include +#include #include +#include +#include +#include #include #include @@ -222,7 +227,8 @@ struct HmpidTableProducer { const double kConvThr = 0.00001; const int kMaxIter = 100; - std::array x, p; + std::array x{}; + std::array p{}; for (int iter = 0; iter < kMaxIter; ++iter) { From 4d708db7a97aa6888ff946c1f726c44074db1192 Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Mon, 13 Jul 2026 20:30:49 +0200 Subject: [PATCH 5/9] Change of type, photon charge --- DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx | 6 ++++-- DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h | 4 +++- 2 files changed, 7 insertions(+), 3 deletions(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx index eb089da60c4..f267b25b443 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx +++ b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx @@ -415,9 +415,11 @@ struct HmpidTableProducer { if (chamberM3 < 0) continue; - float hmpidPhotsCharge2[o2::aod::kDimPhotonsCharge]; - for (int i = 0; i < o2::aod::kDimPhotonsCharge; i++) + std::array hmpidPhotsCharge2{}; + + for (int i = 0; i < o2::aod::kDimPhotonsCharge; i++) { hmpidPhotsCharge2[i] = t.hmpidPhotsCharge()[i]; + } // fill hmpid table hmpidAnalysis( diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h index f3f6a93f843..30b1707189e 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h +++ b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h @@ -15,6 +15,8 @@ #include #include +#include + namespace o2::aod { @@ -33,7 +35,7 @@ DECLARE_SOA_COLUMN(NPhotons, nPhotons, float); DECLARE_SOA_COLUMN(ChargeMip, chargeMip, float); DECLARE_SOA_COLUMN(ClusterSize, clusterSize, float); DECLARE_SOA_COLUMN(Chamber, chamber, float); -DECLARE_SOA_COLUMN(PhotonsCharge, photonsCharge, float[kDimPhotonsCharge]); +DECLARE_SOA_COLUMN(PhotonsCharge, photonsCharge, (std::array)); DECLARE_SOA_COLUMN(EtaTrack, etaTrack, float); DECLARE_SOA_COLUMN(PhiTrack, phiTrack, float); DECLARE_SOA_COLUMN(Px, px, float); From 0f80ae50983706e638cc67414c23344f0cdfbb35 Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Tue, 14 Jul 2026 13:13:39 +0200 Subject: [PATCH 6/9] Fixing table --- DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h index 30b1707189e..282638206b6 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h +++ b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h @@ -24,6 +24,8 @@ inline constexpr int kDimPhotonsCharge = 10; namespace variables_table { +using PhotonsChargeArray = std::array; + DECLARE_SOA_COLUMN(ChAngle, chAngle, float); DECLARE_SOA_COLUMN(MomentumHmpid, momentumHmpid, float); DECLARE_SOA_COLUMN(MomentumTrack, momentumTrack, float); @@ -35,7 +37,7 @@ DECLARE_SOA_COLUMN(NPhotons, nPhotons, float); DECLARE_SOA_COLUMN(ChargeMip, chargeMip, float); DECLARE_SOA_COLUMN(ClusterSize, clusterSize, float); DECLARE_SOA_COLUMN(Chamber, chamber, float); -DECLARE_SOA_COLUMN(PhotonsCharge, photonsCharge, (std::array)); +DECLARE_SOA_COLUMN(PhotonsCharge, photonsCharge, PhotonsChargeArray); DECLARE_SOA_COLUMN(EtaTrack, etaTrack, float); DECLARE_SOA_COLUMN(PhiTrack, phiTrack, float); DECLARE_SOA_COLUMN(Px, px, float); From 377c077cc4b1d1f413416bd0c914c91c235f8b0a Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Tue, 14 Jul 2026 14:01:58 +0200 Subject: [PATCH 7/9] Fixing column type --- DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h index 282638206b6..90df4e273a6 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h +++ b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h @@ -24,7 +24,7 @@ inline constexpr int kDimPhotonsCharge = 10; namespace variables_table { -using PhotonsChargeArray = std::array; +using PhotonsChargeArray = o2::soa::Array; DECLARE_SOA_COLUMN(ChAngle, chAngle, float); DECLARE_SOA_COLUMN(MomentumHmpid, momentumHmpid, float); From e3e3ce67ee32a75dd7002839a3ee52ecc2ef914d Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Tue, 14 Jul 2026 17:04:10 +0200 Subject: [PATCH 8/9] Fixing column type (vector) --- DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx | 2 +- DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h | 5 ++--- 2 files changed, 3 insertions(+), 4 deletions(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx index f267b25b443..a643b313331 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx +++ b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx @@ -415,7 +415,7 @@ struct HmpidTableProducer { if (chamberM3 < 0) continue; - std::array hmpidPhotsCharge2{}; + std::vector hmpidPhotsCharge2(o2::aod::kDimPhotonsCharge, 0.f); for (int i = 0; i < o2::aod::kDimPhotonsCharge; i++) { hmpidPhotsCharge2[i] = t.hmpidPhotsCharge()[i]; diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h index 90df4e273a6..c1fc0365b8f 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h +++ b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h @@ -24,8 +24,7 @@ inline constexpr int kDimPhotonsCharge = 10; namespace variables_table { -using PhotonsChargeArray = o2::soa::Array; - + DECLARE_SOA_COLUMN(ChAngle, chAngle, float); DECLARE_SOA_COLUMN(MomentumHmpid, momentumHmpid, float); DECLARE_SOA_COLUMN(MomentumTrack, momentumTrack, float); @@ -37,7 +36,7 @@ DECLARE_SOA_COLUMN(NPhotons, nPhotons, float); DECLARE_SOA_COLUMN(ChargeMip, chargeMip, float); DECLARE_SOA_COLUMN(ClusterSize, clusterSize, float); DECLARE_SOA_COLUMN(Chamber, chamber, float); -DECLARE_SOA_COLUMN(PhotonsCharge, photonsCharge, PhotonsChargeArray); +DECLARE_SOA_COLUMN(PhotonsCharge, photonsCharge, std::vector); DECLARE_SOA_COLUMN(EtaTrack, etaTrack, float); DECLARE_SOA_COLUMN(PhiTrack, phiTrack, float); DECLARE_SOA_COLUMN(Px, px, float); From 899cb2dee575073f6323e682ca060ab4244ec425 Mon Sep 17 00:00:00 2001 From: Rocco Liotino Date: Tue, 14 Jul 2026 17:38:06 +0200 Subject: [PATCH 9/9] Fixing format --- DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx | 1 + DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h | 4 ++-- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx index a643b313331..d79cfdc7b98 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx +++ b/DPG/Tasks/AOTTrack/PID/HMPID/hmpidTableProducer.cxx @@ -42,6 +42,7 @@ #include #include #include +#include using namespace o2; using namespace o2::framework; diff --git a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h index c1fc0365b8f..1e2445d5b05 100644 --- a/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h +++ b/DPG/Tasks/AOTTrack/PID/HMPID/tableHMPID.h @@ -15,7 +15,7 @@ #include #include -#include +#include namespace o2::aod { @@ -24,7 +24,7 @@ inline constexpr int kDimPhotonsCharge = 10; namespace variables_table { - + DECLARE_SOA_COLUMN(ChAngle, chAngle, float); DECLARE_SOA_COLUMN(MomentumHmpid, momentumHmpid, float); DECLARE_SOA_COLUMN(MomentumTrack, momentumTrack, float);