|
11 | 11 | /// \author Maxim Virta (maxim.virta@cern.ch) |
12 | 12 | /// \since Jul 2024 |
13 | 13 |
|
14 | | -#include <string> |
15 | | -#include <vector> |
| 14 | +#include "JEPFlowAnalysis.h" |
16 | 15 |
|
17 | | -#include "Framework/AnalysisTask.h" |
18 | | -#include "Framework/RunningWorkflowInfo.h" |
19 | | -#include "Framework/HistogramRegistry.h" |
| 16 | +#include "FlowJHistManager.h" |
20 | 17 |
|
21 | | -#include "Common/DataModel/EventSelection.h" |
| 18 | +#include "Common/Core/EventPlaneHelper.h" |
22 | 19 | #include "Common/Core/TrackSelection.h" |
23 | | -#include "Framework/runDataProcessing.h" |
24 | | -#include "Common/DataModel/TrackSelectionTables.h" |
25 | | - |
| 20 | +#include "Common/DataModel/EventSelection.h" |
26 | 21 | #include "Common/DataModel/Qvectors.h" |
27 | | -#include "Common/Core/EventPlaneHelper.h" |
| 22 | +#include "Common/DataModel/TrackSelectionTables.h" |
28 | 23 |
|
29 | | -#include "CCDB/CcdbApi.h" |
30 | 24 | #include "CCDB/BasicCCDBManager.h" |
| 25 | +#include "CCDB/CcdbApi.h" |
| 26 | +#include "Framework/AnalysisTask.h" |
| 27 | +#include "Framework/HistogramRegistry.h" |
| 28 | +#include "Framework/RunningWorkflowInfo.h" |
| 29 | +#include "Framework/runDataProcessing.h" |
31 | 30 |
|
32 | | -#include "FlowJHistManager.h" |
33 | | -#include "JEPFlowAnalysis.h" |
| 31 | +#include <string> |
| 32 | +#include <vector> |
34 | 33 |
|
35 | 34 | using namespace o2; |
36 | 35 | using namespace o2::framework; |
@@ -118,9 +117,9 @@ struct jEPFlowAnalysis { |
118 | 117 | if (cfgAddEvtSel && (!coll.sel8() || !coll.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) || !coll.selection_bit(aod::evsel::kNoSameBunchPileup))) |
119 | 118 | return; |
120 | 119 |
|
121 | | - Float_t cent = coll.cent(); |
| 120 | + float cent = coll.cent(); |
122 | 121 | EPFlowHistograms.fill(HIST("FullCentrality"), cent); |
123 | | - Float_t EPs[3] = {0.}; |
| 122 | + float EPs[3] = {0.}; |
124 | 123 |
|
125 | 124 | if (cfgShiftCorr) { |
126 | 125 | auto bc = coll.bc_as<aod::BCsWithTimestamps>(); |
@@ -157,28 +156,29 @@ struct jEPFlowAnalysis { |
157 | 156 | auto coeffshiftxRefA = shiftprofile.at(i - 2)->GetBinContent(shiftprofile.at(i - 2)->FindBin(cent, 2.5, ishift - 0.5)); |
158 | 157 | auto coeffshiftyRefA = shiftprofile.at(i - 2)->GetBinContent(shiftprofile.at(i - 2)->FindBin(cent, 3.5, ishift - 0.5)); |
159 | 158 | auto coeffshiftxRefB = shiftprofile.at(i - 2)->GetBinContent(shiftprofile.at(i - 2)->FindBin(cent, 4.5, ishift - 0.5)); |
160 | | - auto coeffshiftyRefB = shiftprofile.at(i - 2)->GetBinContent(shiftprofile.at(i - 2)->FindBin(cent, 5.5, ishift - 0.5)); //currently only FT0C/TPCpos/TPCneg |
| 159 | + auto coeffshiftyRefB = shiftprofile.at(i - 2)->GetBinContent(shiftprofile.at(i - 2)->FindBin(cent, 5.5, ishift - 0.5)); // currently only FT0C/TPCpos/TPCneg |
161 | 160 |
|
162 | | - 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]))); |
163 | | - 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]))); |
164 | | - deltapsiRefB += ((1 / (1.0 * ishift)) * (-coeffshiftxRefB * TMath::Cos(ishift * static_cast<float>(i) * EPs[2]) + coeffshiftyRefB * TMath::Sin(ishift * static_cast<float>(i) * EPs[2]))); |
| 161 | + deltapsiDet += ((1 / (1.0 * ishift)) * (-coeffshiftxDet * std::cos(ishift * static_cast<float>(i) * EPs[0]) + coeffshiftyDet * std::sin(ishift * static_cast<float>(i) * EPs[0]))); |
| 162 | + deltapsiRefA += ((1 / (1.0 * ishift)) * (-coeffshiftxRefA * std::cos(ishift * static_cast<float>(i) * EPs[1]) + coeffshiftyRefA * std::sin(ishift * static_cast<float>(i) * EPs[1]))); |
| 163 | + deltapsiRefB += ((1 / (1.0 * ishift)) * (-coeffshiftxRefB * std::cos(ishift * static_cast<float>(i) * EPs[2]) + coeffshiftyRefB * std::sin(ishift * static_cast<float>(i) * EPs[2]))); |
165 | 164 | } |
166 | 165 |
|
167 | 166 | EPs[0] += deltapsiDet; |
168 | 167 | EPs[1] += deltapsiRefA; |
169 | 168 | EPs[2] += deltapsiRefB; |
170 | 169 | } |
171 | 170 |
|
172 | | - if (cfgSPmethod) weight *= sqrt(pow(coll.qvecRe()[DetId + harmInd], 2) + pow(coll.qvecIm()[DetId + harmInd], 2)); |
| 171 | + if (cfgSPmethod) |
| 172 | + weight *= sqrt(pow(coll.qvecRe()[DetId + harmInd], 2) + pow(coll.qvecIm()[DetId + harmInd], 2)); |
173 | 173 |
|
174 | | - Float_t resNumA = helperEP.GetResolution(EPs[0], EPs[1], i); |
175 | | - Float_t resNumB = helperEP.GetResolution(EPs[0], EPs[2], i); |
176 | | - Float_t resDenom = helperEP.GetResolution(EPs[1], EPs[2], i); |
| 174 | + float resNumA = helperEP.GetResolution(EPs[0], EPs[1], i); |
| 175 | + float resNumB = helperEP.GetResolution(EPs[0], EPs[2], i); |
| 176 | + float resDenom = helperEP.GetResolution(EPs[1], EPs[2], i); |
177 | 177 | epAnalysis.FillResolutionHistograms(cent, static_cast<float>(i), resNumA, resNumB, resDenom); |
178 | 178 | for (uint j = 0; j < 3; j++) { // loop over detectors used |
179 | 179 | for (auto& track : tracks) { |
180 | | - Float_t vn = TMath::Cos((i) * (track.phi() - EPs[j])); |
181 | | - Float_t vn_sin = TMath::Sin((i) * (track.phi() - EPs[j])); |
| 180 | + float vn = std::cos((i) * (track.phi() - EPs[j])); |
| 181 | + float vn_sin = std::sin((i) * (track.phi() - EPs[j])); |
182 | 182 | epAnalysis.FillVnHistograms(i, cent, static_cast<float>(j + 1), track.pt(), vn * weight, vn_sin * weight); |
183 | 183 | } |
184 | 184 | } |
|
0 commit comments