1414// / \since 01/12/2024
1515// / \brief task to evaluate flow with respect to spectator plane.
1616
17- #include < algorithm>
18- #include < numeric>
19- #include < vector>
20- #include < string>
21- #include < unordered_map>
17+ #include " GFWWeights.h"
2218
23- #include " Framework/runDataProcessing.h"
24- #include " Framework/AnalysisTask.h"
25- #include " Framework/ASoAHelpers.h"
26- #include " Framework/RunningWorkflowInfo.h"
27- #include " Framework/HistogramRegistry.h"
28- #include " Framework/O2DatabasePDGPlugin.h"
29- #include " CCDB/BasicCCDBManager.h"
30- #include " DataFormatsParameters/GRPObject.h"
31- #include " DataFormatsParameters/GRPMagField.h"
32- #include " DataFormatsParameters/GRPLHCIFData.h"
19+ #include " PWGCF/DataModel/SPTableZDC.h"
3320
34- #include " Common/DataModel/EventSelection .h"
21+ #include " Common/Core/RecoDecay .h"
3522#include " Common/Core/TrackSelection.h"
36- #include " Common/DataModel/TrackSelectionTables.h"
37- #include " Common/DataModel/Multiplicity.h"
3823#include " Common/DataModel/Centrality.h"
39- #include " Common/Core/RecoDecay.h"
24+ #include " Common/DataModel/EventSelection.h"
25+ #include " Common/DataModel/Multiplicity.h"
4026#include " Common/DataModel/PIDResponse.h"
27+ #include " Common/DataModel/TrackSelectionTables.h"
28+
29+ #include " CCDB/BasicCCDBManager.h"
30+ #include " DataFormatsParameters/GRPLHCIFData.h"
31+ #include " DataFormatsParameters/GRPMagField.h"
32+ #include " DataFormatsParameters/GRPObject.h"
33+ #include " Framework/ASoAHelpers.h"
34+ #include " Framework/AnalysisTask.h"
35+ #include " Framework/HistogramRegistry.h"
36+ #include " Framework/O2DatabasePDGPlugin.h"
37+ #include " Framework/RunningWorkflowInfo.h"
38+ #include " Framework/runDataProcessing.h"
4139
42- #include " PWGCF/DataModel/SPTableZDC.h"
43- #include " GFWWeights.h"
4440#include " TF1.h"
4541#include " TPDGCode.h"
4642
43+ #include < algorithm>
44+ #include < numeric>
45+ #include < string>
46+ #include < unordered_map>
47+ #include < vector>
48+
4749using namespace o2 ;
4850using namespace o2 ::framework;
4951using namespace o2 ::framework::expressions;
@@ -376,7 +378,7 @@ struct FlowSP {
376378 registry.add (" incl/QA/after/hCrossedRows_pt" , " " , {HistType::kTH2D , {axisPt, axisCl}});
377379 registry.add (" incl/QA/after/hCrossedRows_vs_SharedClusters" , " " , {HistType::kTH2D , {axisCl, axisShCl}});
378380
379- if (cfgTrackSelDoTrackQAvsCent){
381+ if (cfgTrackSelDoTrackQAvsCent) {
380382 registry.add <TH2>(" incl/QA/after/hPt" , " " , kTH2D , {axisPt, axisCent});
381383 registry.add <TH2>(" incl/QA/after/hPt_forward" , " " , kTH2D , {axisPt, axisCent});
382384 registry.add <TH2>(" incl/QA/after/hPt_forward_uncorrected" , " " , kTH2D , {axisPt, axisCent});
@@ -386,8 +388,8 @@ struct FlowSP {
386388 registry.add <TH2>(" incl/QA/after/hPhi_uncorrected" , " " , kTH2D , {axisPhi, axisCent});
387389 registry.add <TH2>(" incl/QA/after/hEta" , " " , kTH2D , {axisEta, axisCent});
388390 registry.add <TH2>(" incl/QA/after/hEta_uncorrected" , " " , kTH2D , {axisEta, axisCent});
389- } else {
390- registry.add <TH1>(" incl/QA/after/hPt" , " " , kTH1D , {axisPt});
391+ } else {
392+ registry.add <TH1>(" incl/QA/after/hPt" , " " , kTH1D , {axisPt});
391393 registry.add <TH1>(" incl/QA/after/hPt_forward" , " " , kTH1D , {axisPt});
392394 registry.add <TH1>(" incl/QA/after/hPt_forward_uncorrected" , " " , kTH1D , {axisPt});
393395 registry.add <TH1>(" incl/QA/after/hPt_backward" , " " , kTH1D , {axisPt});
@@ -656,7 +658,7 @@ struct FlowSP {
656658 return -1 ;
657659 }
658660
659- int getMagneticField (uint64_t timestamp)
661+ int getMagneticField (uint64_t timestamp)
660662 {
661663 // TODO done only once (and not per run). Will be replaced by CCDBConfigurable
662664 static o2::parameters::GRPMagField* grpo = nullptr ;
@@ -674,14 +676,14 @@ struct FlowSP {
674676 std::pair<float , long > getCrossingAngleCCDB (uint64_t timestamp)
675677 {
676678 // TODO done only once (and not per run). Will be replaced by CCDBConfigurable
677- auto grpo = ccdb->getForTimeStamp <o2::parameters::GRPLHCIFData>(" GLO/Config/GRPLHCIF" , timestamp);
678- if (grpo == nullptr ) {
679- LOGF (fatal, " GRP object for Crossing Angle not found for timestamp %llu" , timestamp);
680- return {0 , 0 };
681- }
682- float crossingAngle = grpo->getCrossingAngle ();
683- long crossingAngleTime = grpo->getCrossingAngleTime ();
684- return {crossingAngle, crossingAngleTime};
679+ auto grpo = ccdb->getForTimeStamp <o2::parameters::GRPLHCIFData>(" GLO/Config/GRPLHCIF" , timestamp);
680+ if (grpo == nullptr ) {
681+ LOGF (fatal, " GRP object for Crossing Angle not found for timestamp %llu" , timestamp);
682+ return {0 , 0 };
683+ }
684+ float crossingAngle = grpo->getCrossingAngle ();
685+ long crossingAngleTime = grpo->getCrossingAngleTime ();
686+ return {crossingAngle, crossingAngleTime};
685687 }
686688
687689 // From Generic Framework
@@ -1131,7 +1133,6 @@ struct FlowSP {
11311133 registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_vs_SharedClusters" ), track.tpcNClsFound (), track.tpcFractionSharedCls (), wacc * weff);
11321134 }
11331135
1134-
11351136 template <FillType ft, ChargeType ct, typename TrackObject>
11361137 inline void fillPIDQA (TrackObject track)
11371138 {
@@ -1201,21 +1202,21 @@ struct FlowSP {
12011202 template <FillType ft, ParticleType ct, typename TrackObject>
12021203 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 )
12031204 {
1204- if (!cfgTrackSelDoTrackQAvsCent) {
1205+ if (!cfgTrackSelDoTrackQAvsCent) {
12051206 fillTrackQA<ft, kInclusive , ct>(track, vtxz, wacc, weff);
12061207 } else {
12071208 fillTrackQA<ft, kInclusive , ct>(track, vtxz, centrality, wacc, weff);
1208- }
1209+ }
12091210 fillPIDQA<ft, kInclusive >(track);
12101211 if (pos) {
1211- if (!cfgTrackSelDoTrackQAvsCent) {
1212- fillTrackQA<ft, kPositive , ct>(track, vtxz, waccP, weffP);
1213- } else {
1214- fillTrackQA<ft, kPositive , ct>(track, vtxz, centrality, waccP, weffP);
1215- }
1212+ if (!cfgTrackSelDoTrackQAvsCent) {
1213+ fillTrackQA<ft, kPositive , ct>(track, vtxz, waccP, weffP);
1214+ } else {
1215+ fillTrackQA<ft, kPositive , ct>(track, vtxz, centrality, waccP, weffP);
1216+ }
12161217 fillPIDQA<ft, kPositive >(track);
12171218 } else {
1218- if (!cfgTrackSelDoTrackQAvsCent) {
1219+ if (!cfgTrackSelDoTrackQAvsCent) {
12191220 fillTrackQA<ft, kNegative , ct>(track, vtxz, waccN, weffN);
12201221 } else {
12211222 fillTrackQA<ft, kNegative , ct>(track, vtxz, centrality, waccN, weffN);
0 commit comments