1111// / \author Maxim Virta (maxim.virta@cern.ch)
1212// / \since Jul 2024
1313
14- #include < string>
15- #include < vector>
16-
1714#include " Framework/AnalysisTask.h"
1815#include " Framework/RunningWorkflowInfo.h"
1916#include " Framework/HistogramRegistry.h"
2623#include " Common/DataModel/Qvectors.h"
2724#include " Common/Core/EventPlaneHelper.h"
2825
29- #include " FlowJHistManager.h"
30- #include " JEPFlowAnalysis.h"
31-
3226#include " CCDB/CcdbApi.h"
3327#include " CCDB/BasicCCDBManager.h"
3428
29+ #include " FlowJHistManager.h"
30+ #include " JEPFlowAnalysis.h"
31+
3532using namespace o2 ;
3633using namespace o2 ::framework;
3734using namespace o2 ::framework::expressions;
@@ -47,11 +44,9 @@ struct jEPFlowAnalysis {
4744 EventPlaneHelper helperEP;
4845 FlowJHistManager histManager;
4946 Bool_t debug = kFALSE ;
50-
5147 Service<o2::ccdb::BasicCCDBManager> ccdb;
5248 o2::ccdb::CcdbApi ccdbApi;
5349
54-
5550 // Set Configurables here
5651 struct : ConfigurableGroup {
5752 Configurable<float > cfgPtMin{" cfgPtMin" , 0 .2f , " Minimum pT used for track selection." };
@@ -125,7 +120,7 @@ struct jEPFlowAnalysis {
125120 Float_t EPs[3 ] = {0 .};
126121
127122 if (cfgShiftCorr) {
128- auto bc = collision .bc_as <aod::BCsWithTimestamps>();
123+ auto bc = coll .bc_as <aod::BCsWithTimestamps>();
129124 currentRunNumber = bc.runNumber ();
130125 if (currentRunNumber != lastRunNumber) {
131126 shiftprofile.clear ();
@@ -154,12 +149,12 @@ struct jEPFlowAnalysis {
154149
155150 if (cfgShiftCorr) {
156151 for (int ishift = 1 ; ishift <= 10 ; ishift++) {
157- auto coeffshiftxDet = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (centrality , 0.5 , ishift - 0.5 ));
158- auto coeffshiftyDet = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (centrality , 1.5 , ishift - 0.5 ));
159- auto coeffshiftxRefA = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (centrality , 2.5 , ishift - 0.5 ));
160- auto coeffshiftyRefA = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (centrality , 3.5 , ishift - 0.5 ));
161- auto coeffshiftxRefB = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (centrality , 4.5 , ishift - 0.5 ));
162- auto coeffshiftyRefB = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (centrality , 5.5 , ishift - 0.5 )); // currently only FT0C/TPCpos/TPCneg
152+ auto coeffshiftxDet = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (cent , 0.5 , ishift - 0.5 ));
153+ auto coeffshiftyDet = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (cent , 1.5 , ishift - 0.5 ));
154+ auto coeffshiftxRefA = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (cent , 2.5 , ishift - 0.5 ));
155+ auto coeffshiftyRefA = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (cent , 3.5 , ishift - 0.5 ));
156+ auto coeffshiftxRefB = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (cent , 4.5 , ishift - 0.5 ));
157+ auto coeffshiftyRefB = shiftprofile.at (i - 2 )->GetBinContent (shiftprofile.at (i - 2 )->FindBin (cent , 5.5 , ishift - 0.5 )); // currently only FT0C/TPCpos/TPCneg
163158
164159 deltapsiDet += ((1 / (1.0 * ishift)) * (-coeffshiftxDet * TMath::Cos (ishift * static_cast <float >(i) * EPs[0 ]) + coeffshiftyDet * TMath::Sin (ishift * static_cast <float >(i) * EPs[0 ])));
165160 deltapsiRefA += ((1 / (1.0 * ishift)) * (-coeffshiftxRefA * TMath::Cos (ishift * static_cast <float >(i) * EPs[1 ]) + coeffshiftyRefA * TMath::Sin (ishift * static_cast <float >(i) * EPs[1 ])));
0 commit comments