1414// / \since Apr/2/2026
1515// / \brief flow efficiency analysis on UPC MC
1616
17- #include < CCDB/BasicCCDBManager.h>
1817#include " PWGUD/Core/SGSelector.h"
1918#include " PWGUD/DataModel/UDTables.h"
20- #include < vector>
21- #include < string>
22- #include " Framework/AnalysisDataModel.h"
23- #include " Framework/AnalysisTask.h"
24- #include " Framework/HistogramRegistry.h"
25- #include " Framework/runDataProcessing.h"
26- #include " Framework/ASoAHelpers.h"
27- #include " Framework/RunningWorkflowInfo.h"
28- #include " Common/DataModel/TrackSelectionTables.h"
19+
2920#include " Common/Core/RecoDecay.h"
3021#include " Common/Core/TrackSelection.h"
3122#include " Common/Core/TrackSelectionDefaults.h"
3223#include " Common/Core/trackUtilities.h"
24+ #include " Common/DataModel/TrackSelectionTables.h"
25+
26+ #include " Framework/ASoAHelpers.h"
27+ #include " Framework/AnalysisDataModel.h"
28+ #include " Framework/AnalysisTask.h"
29+ #include " Framework/HistogramRegistry.h"
30+ #include " Framework/RunningWorkflowInfo.h"
31+ #include " Framework/runDataProcessing.h"
3332#include " ReconstructionDataFormats/Track.h"
33+ #include < CCDB/BasicCCDBManager.h>
34+
35+ #include < TF1.h>
36+ #include < TPDGCode.h>
3437#include < TProfile.h>
3538#include < TRandom3.h>
36- #include < TPDGCode.h>
37- #include < TF1.h>
39+
40+ #include < string>
41+ #include < vector>
3842
3943using namespace o2 ;
4044using namespace o2 ::framework;
@@ -54,7 +58,6 @@ struct FlowMcUpc {
5458 O2_DEFINE_CONFIGURABLE (cfgCutDCAxy, float , 0 .2f , " DCAxy cut for tracks" )
5559 O2_DEFINE_CONFIGURABLE (cfgDcaxy, bool , true , " choose dcaxy" )
5660
57-
5861 ConfigurableAxis axisB{" axisB" , {100 , 0 .0f , 20 .0f }, " " };
5962 ConfigurableAxis axisPt{" axisPt" , {VARIABLE_WIDTH, 0 .0f , 0 .1f , 0 .2f , 0 .3f , 0 .4f , 0 .5f , 0 .6f , 0 .7f , 0 .8f , 0 .9f , 1 .0f , 1 .1f , 1 .2f , 1 .3f , 1 .4f , 1 .5f , 1 .6f , 1 .7f , 1 .8f , 1 .9f , 2 .0f , 2 .2f , 2 .4f , 2 .6f , 2 .8f , 3 .0f , 3 .2f , 3 .4f , 3 .6f , 3 .8f , 4 .0f , 4 .4f , 4 .8f , 5 .2f , 5 .6f , 6 .0f , 6 .5f , 7 .0f , 7 .5f , 8 .0f , 9 .0f , 10 .0f , 11 .0f , 12 .0f }, " pt axis" };
6063 // Connect to ccdb
@@ -65,7 +68,6 @@ struct FlowMcUpc {
6568
6669 using McParticles = soa::Join<aod::UDMcParticles, aod::UDMcTrackLabels>;
6770
68-
6971 void init (InitContext&)
7072 {
7173 ccdb->setURL (ccdbUrl.value );
@@ -85,7 +87,6 @@ struct FlowMcUpc {
8587 histos.add <TH3>(" hEtaPtVtxzMCGen" , " Monte Carlo Truth; #eta; p_{T} (GeV/c); V_{z} (cm);" , {HistType::kTH3D , {axisEta, axisPt, axisVertex}});
8688 histos.add <TH1>(" hPtReco" , " Monte Carlo Reco Global; pT (GeV/c);" , {HistType::kTH1D , {axisPt}});
8789 histos.add <TH3>(" hEtaPtVtxzMCReco" , " Monte Carlo Global; #eta; p_{T} (GeV/c); V_{z} (cm);" , {HistType::kTH3D , {axisEta, axisPt, axisVertex}});
88-
8990 }
9091
9192 // template <typename TCollision>
@@ -136,12 +137,11 @@ struct FlowMcUpc {
136137 continue ;
137138 // if (std::fabs(mcParticle.eta()) > cfgCutEta) // main acceptance
138139 // continue;
139-
140+
140141 histos.fill (HIST (" hPtMCGen" ), pt);
141142 histos.fill (HIST (" hEtaPtVtxzMCGen" ), eta, pt, vtxz);
142143 }
143144 }
144-
145145 }
146146 PROCESS_SWITCH (FlowMcUpc, processMCTrue, " process pure simulation information" , true );
147147
@@ -165,9 +165,9 @@ struct FlowMcUpc {
165165
166166 for (const auto & track : tracks) {
167167 auto momentum = std::array<double , 3 >{track.px (), track.py (), track.pz ()};
168- double pt = RecoDecay::pt (momentum);
169- double phi = RecoDecay::phi (momentum);
170- double eta = RecoDecay::eta (momentum);
168+ double pt = RecoDecay::pt (momentum);
169+ double phi = RecoDecay::phi (momentum);
170+ double eta = RecoDecay::eta (momentum);
171171 if (!trackSelected (track) || (!track.has_udMcParticle ()))
172172 continue ;
173173 auto mcParticle = track.udMcParticle ();
@@ -178,7 +178,7 @@ struct FlowMcUpc {
178178 // continue;
179179 if (!mcParticle.isPhysicalPrimary ())
180180 continue ;
181-
181+
182182 histos.fill (HIST (" hPtReco" ), pt);
183183 histos.fill (HIST (" hEtaPtVtxzMCReco" ), eta, pt, vtxz);
184184 }
0 commit comments