diff --git a/PWGCF/Flow/Tasks/flowSP.cxx b/PWGCF/Flow/Tasks/flowSP.cxx index 1156a8f4f43..a50fa73859e 100644 --- a/PWGCF/Flow/Tasks/flowSP.cxx +++ b/PWGCF/Flow/Tasks/flowSP.cxx @@ -14,36 +14,38 @@ /// \since 01/12/2024 /// \brief task to evaluate flow with respect to spectator plane. -#include -#include -#include -#include -#include +#include "GFWWeights.h" -#include "Framework/runDataProcessing.h" -#include "Framework/AnalysisTask.h" -#include "Framework/ASoAHelpers.h" -#include "Framework/RunningWorkflowInfo.h" -#include "Framework/HistogramRegistry.h" -#include "Framework/O2DatabasePDGPlugin.h" -#include "CCDB/BasicCCDBManager.h" -#include "DataFormatsParameters/GRPObject.h" -#include "DataFormatsParameters/GRPMagField.h" -#include "DataFormatsParameters/GRPLHCIFData.h" +#include "PWGCF/DataModel/SPTableZDC.h" -#include "Common/DataModel/EventSelection.h" +#include "Common/Core/RecoDecay.h" #include "Common/Core/TrackSelection.h" -#include "Common/DataModel/TrackSelectionTables.h" -#include "Common/DataModel/Multiplicity.h" #include "Common/DataModel/Centrality.h" -#include "Common/Core/RecoDecay.h" +#include "Common/DataModel/EventSelection.h" +#include "Common/DataModel/Multiplicity.h" #include "Common/DataModel/PIDResponse.h" +#include "Common/DataModel/TrackSelectionTables.h" + +#include "CCDB/BasicCCDBManager.h" +#include "DataFormatsParameters/GRPLHCIFData.h" +#include "DataFormatsParameters/GRPMagField.h" +#include "DataFormatsParameters/GRPObject.h" +#include "Framework/ASoAHelpers.h" +#include "Framework/AnalysisTask.h" +#include "Framework/HistogramRegistry.h" +#include "Framework/O2DatabasePDGPlugin.h" +#include "Framework/RunningWorkflowInfo.h" +#include "Framework/runDataProcessing.h" -#include "PWGCF/DataModel/SPTableZDC.h" -#include "GFWWeights.h" #include "TF1.h" #include "TPDGCode.h" +#include +#include +#include +#include +#include + using namespace o2; using namespace o2::framework; using namespace o2::framework::expressions; @@ -376,7 +378,7 @@ struct FlowSP { registry.add("incl/QA/after/hCrossedRows_pt", "", {HistType::kTH2D, {axisPt, axisCl}}); registry.add("incl/QA/after/hCrossedRows_vs_SharedClusters", "", {HistType::kTH2D, {axisCl, axisShCl}}); - if(cfgTrackSelDoTrackQAvsCent){ + if (cfgTrackSelDoTrackQAvsCent) { registry.add("incl/QA/after/hPt", "", kTH2D, {axisPt, axisCent}); registry.add("incl/QA/after/hPt_forward", "", kTH2D, {axisPt, axisCent}); registry.add("incl/QA/after/hPt_forward_uncorrected", "", kTH2D, {axisPt, axisCent}); @@ -386,8 +388,8 @@ struct FlowSP { registry.add("incl/QA/after/hPhi_uncorrected", "", kTH2D, {axisPhi, axisCent}); registry.add("incl/QA/after/hEta", "", kTH2D, {axisEta, axisCent}); registry.add("incl/QA/after/hEta_uncorrected", "", kTH2D, {axisEta, axisCent}); - } else { - registry.add("incl/QA/after/hPt", "", kTH1D, {axisPt}); + } else { + registry.add("incl/QA/after/hPt", "", kTH1D, {axisPt}); registry.add("incl/QA/after/hPt_forward", "", kTH1D, {axisPt}); registry.add("incl/QA/after/hPt_forward_uncorrected", "", kTH1D, {axisPt}); registry.add("incl/QA/after/hPt_backward", "", kTH1D, {axisPt}); @@ -656,7 +658,7 @@ struct FlowSP { return -1; } - int getMagneticField(uint64_t timestamp) + int getMagneticField(uint64_t timestamp) { // TODO done only once (and not per run). Will be replaced by CCDBConfigurable static o2::parameters::GRPMagField* grpo = nullptr; @@ -674,14 +676,14 @@ struct FlowSP { std::pair getCrossingAngleCCDB(uint64_t timestamp) { // TODO done only once (and not per run). Will be replaced by CCDBConfigurable - auto grpo = ccdb->getForTimeStamp("GLO/Config/GRPLHCIF", timestamp); - if (grpo == nullptr) { - LOGF(fatal, "GRP object for Crossing Angle not found for timestamp %llu", timestamp); - return {0, 0}; - } - float crossingAngle = grpo->getCrossingAngle(); - long crossingAngleTime = grpo->getCrossingAngleTime(); - return {crossingAngle, crossingAngleTime}; + auto grpo = ccdb->getForTimeStamp("GLO/Config/GRPLHCIF", timestamp); + if (grpo == nullptr) { + LOGF(fatal, "GRP object for Crossing Angle not found for timestamp %llu", timestamp); + return {0, 0}; + } + float crossingAngle = grpo->getCrossingAngle(); + long crossingAngleTime = grpo->getCrossingAngleTime(); + return {crossingAngle, crossingAngleTime}; } // From Generic Framework @@ -1131,7 +1133,6 @@ struct FlowSP { registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("QA/") + HIST(Time[ft]) + HIST("hCrossedRows_vs_SharedClusters"), track.tpcNClsFound(), track.tpcFractionSharedCls(), wacc * weff); } - template inline void fillPIDQA(TrackObject track) { @@ -1201,21 +1202,21 @@ struct FlowSP { template void fillAllQA(TrackObject track, double vtxz, double centrality, bool pos, float wacc = 1, float weff = 1, float waccP = 1, float weffP = 1, float waccN = 1, float weffN = 1) { - if(!cfgTrackSelDoTrackQAvsCent) { + if (!cfgTrackSelDoTrackQAvsCent) { fillTrackQA(track, vtxz, wacc, weff); } else { fillTrackQA(track, vtxz, centrality, wacc, weff); - } + } fillPIDQA(track); if (pos) { - if(!cfgTrackSelDoTrackQAvsCent) { - fillTrackQA(track, vtxz, waccP, weffP); - } else { - fillTrackQA(track, vtxz, centrality, waccP, weffP); - } + if (!cfgTrackSelDoTrackQAvsCent) { + fillTrackQA(track, vtxz, waccP, weffP); + } else { + fillTrackQA(track, vtxz, centrality, waccP, weffP); + } fillPIDQA(track); } else { - if(!cfgTrackSelDoTrackQAvsCent) { + if (!cfgTrackSelDoTrackQAvsCent) { fillTrackQA(track, vtxz, waccN, weffN); } else { fillTrackQA(track, vtxz, centrality, waccN, weffN);