|
17 | 17 | #include "Common/DataModel/Centrality.h" |
18 | 18 | #include "Common/DataModel/EventSelection.h" |
19 | 19 | #include "Common/DataModel/Multiplicity.h" |
| 20 | +#include "Common/DataModel/PIDResponseITS.h" |
20 | 21 | #include "Common/DataModel/PIDResponseTOF.h" |
21 | 22 | #include "Common/DataModel/PIDResponseTPC.h" |
22 | | -#include "Common/DataModel/PIDResponseITS.h" |
23 | 23 | #include "Common/DataModel/Qvectors.h" |
24 | 24 | #include "Common/DataModel/TrackSelectionTables.h" |
25 | 25 |
|
@@ -202,27 +202,27 @@ struct F0980pbpbanalysis { |
202 | 202 | int nmode = FlowConfig.cfgQvecNMods; |
203 | 203 | static constexpr double QvecAmpMin = 1e-4; |
204 | 204 |
|
205 | | - int detId; |
206 | | - int refAId; |
207 | | - int refBId; |
| 205 | + int detId = 0; |
| 206 | + int refAId = 4; |
| 207 | + int refBId = 5; |
208 | 208 |
|
209 | | - int qVecDetInd; |
210 | | - int qVecRefAInd; |
211 | | - int qVecRefBInd; |
| 209 | + int qVecDetInd = 0; |
| 210 | + int qVecRefAInd = 0; |
| 211 | + int qVecRefBInd = 0; |
212 | 212 |
|
213 | | - double eventPlaneDet; |
214 | | - double eventPlaneRefA; |
215 | | - double eventPlaneRefB; |
| 213 | + double eventPlaneDet = 0.; |
| 214 | + double eventPlaneRefA = 0.; |
| 215 | + double eventPlaneRefB = 0.; |
216 | 216 |
|
217 | | - float centrality; |
| 217 | + float centrality =. 1.; |
218 | 218 |
|
219 | | - double angle; |
220 | | - double relPhi; |
221 | | - double relPhiRot; |
222 | | - double relPhiMix; |
| 219 | + double angle = -1.; |
| 220 | + double relPhi = -1.; |
| 221 | + double relPhiRot = -1.; |
| 222 | + double relPhiMix = -1.; |
223 | 223 |
|
224 | 224 | // double massPi = o2::constants::physics::MassPionCharged; |
225 | | - double massPtl; |
| 225 | + double massPtl = 0.; |
226 | 226 |
|
227 | 227 | int nTotalEvents = 0; |
228 | 228 |
|
@@ -597,15 +597,15 @@ struct F0980pbpbanalysis { |
597 | 597 | template <typename TrackType1, typename TrackType2> |
598 | 598 | bool pairAngleSelection(const TrackType1& track1, const TrackType2& track2) |
599 | 599 | { |
600 | | - double pt1, pt2, pz1, pz2, p1, p2, angle; |
| 600 | + double pt1, pt2, pz1, pz2, p1, p2, angle12; |
601 | 601 | pt1 = track1.pt(); |
602 | 602 | pt2 = track2.pt(); |
603 | 603 | pz1 = track1.pz(); |
604 | 604 | pz2 = track2.pz(); |
605 | 605 | p1 = track1.p(); |
606 | 606 | p2 = track2.p(); |
607 | | - angle = std::acos((pt1 * pt2 + pz1 * pz2) / (p1 * p2)); |
608 | | - if (PhiTestConfig.cfgPhiDeepAngleSel && angle < PhiTestConfig.cfgPhiDeepAngle) { |
| 607 | + angle12 = std::acos((pt1 * pt2 + pz1 * pz2) / (p1 * p2)); |
| 608 | + if (PhiTestConfig.cfgPhiDeepAngleSel && angle12 < PhiTestConfig.cfgPhiDeepAngle) { |
609 | 609 | return 0; |
610 | 610 | } |
611 | 611 | return 1; |
@@ -736,7 +736,7 @@ struct F0980pbpbanalysis { |
736 | 736 | if (c1.bcId() == c2.bcId()) { |
737 | 737 | continue; |
738 | 738 | } |
739 | | - double eventPlaneDet = std::atan2(c1.qvecIm()[qVecDetInd], c1.qvecRe()[qVecDetInd]) / static_cast<float>(nmode); |
| 739 | + eventPlaneDet = std::atan2(c1.qvecIm()[qVecDetInd], c1.qvecRe()[qVecDetInd]) / static_cast<float>(nmode); |
740 | 740 |
|
741 | 741 | for (const auto& trk1 : t1) { |
742 | 742 | if (!trackSelected(trk1, false)) { |
|
0 commit comments