fix array

Author: cnkoster

PWGCF/Flow/Tasks/flowSP.cxx

@@ -700,10 +700,10 @@ struct FlowSP {
       fMultCutHigh = new TF1("fMultCutHigh", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x", 0, 100);
 
       std::vector<double> paramsMultPVCut;
-      std::vector<double> paramsMultCut;
-      int y2023 = 2023;
-      int y2024 = 2024;
-      std::array<int> nSigma = {1, 2, 3};
+      std::vector<double> paramsMultCut; 
+      int y2023 = 2023; 
+      int y2024 = 2024; 
+      std::array<int,3> nSigma = {1,2,3}; 
 
       if (cfg.cUsePredeFinedSigma) {
         if (cfg.cUsePredeFinedSigmaYear == y2023) {

Tutorials/PWGCF/EventPlane/src/spectatorPlaneTutorial.cxx

@@ -102,8 +102,8 @@ struct SpectatorPlaneTutorial {
     AxisSpec axisEta = {64, -1.6, 1.6, "#eta"};
     AxisSpec axisEtaVn = {8, -.8, .8, "#eta"};
     AxisSpec axisCent = {90, 0, 90, "Centrality(%)"};
-    AxisSpec axisPhiPlane = {100, -constants::math::PI, constants::math::PI, "#Psi"};
-    AxisSpec axisQQ = {100, -0.2, 0.2, "#LT Q_{X}^{A}Q_{Y}^{C} #GT"};
+    AxisSpec axisPhiPlane = {40, -constants::math::PI, constants::math::PI, "#Psi"};
+    AxisSpec axisQQ = {10, -0.2, 0.2, "#LT Q_{X}^{A}Q_{Y}^{C} #GT"};
 
     std::vector<double> ptbinning = {0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2, 2.2, 2.4, 2.6, 2.8, 3, 3.5, 4, 5, 6, 8, 10};
     AxisSpec axisPt = {ptbinning, "#it{p}_{T} GeV/#it{c}"};
@@ -127,13 +127,13 @@ struct SpectatorPlaneTutorial {
     registry.add("CalibHistos/hEvPlaneRes", "Event Plane Resolution; #Events; Event Plane Resolution", {HistType::kTProfile, {axisCent}});
 
     // Flow Histograms
-    registry.add("flow/v1A", "", {HistType::kTProfile, {axisPt}});
-    registry.add("flow/v1C", "", {HistType::kTProfile, {axisPt}});
+    registry.add("flow/v1A", "", {HistType::kTProfile, {axisEtaVn}});
+    registry.add("flow/v1C", "", {HistType::kTProfile, {axisEtaVn}});
 
-    registry.add("flow/vnAxCxUxMH", "", {HistType::kTProfile, {axisCent}});
-    registry.add("flow/vnAyCyUxMH", "", {HistType::kTProfile, {axisCent}});
-    registry.add("flow/vnAxCyUyMH", "", {HistType::kTProfile, {axisCent}});
-    registry.add("flow/vnAyCxUyMH", "", {HistType::kTProfile, {axisCent}});
+    registry.add("flow/v2AxCxUxMH", "", {HistType::kTProfile, {axisCent}});
+    registry.add("flow/v2AyCyUxMH", "", {HistType::kTProfile, {axisCent}});
+    registry.add("flow/v2AxCyUyMH", "", {HistType::kTProfile, {axisCent}});
+    registry.add("flow/v2AyCxUyMH", "", {HistType::kTProfile, {axisCent}});
   }
 
   void process(ZDCCollisions::iterator const& collision, aod::BCsWithTimestamps const&, UsedTracks const& tracks)
@@ -147,7 +147,7 @@ struct SpectatorPlaneTutorial {
     float centMin = 0;
     float centMax = 80;
 
-    if (centrality > centMax || centrality < centMin)
+    if (centrality >= centMax || centrality < centMin)
       return;
 
     if (collision.isSelected() == false)
@@ -174,10 +174,10 @@ struct SpectatorPlaneTutorial {
     // Fill the q-vector correlations
     registry.fill(HIST("qAqCXY"), centrality, qxA * qxC + qyA * qyC);
 
-    double corrQQx = 1;
-    double corrQQy = 1;
-    double corrQQ = 1;
-    double evPlaneRes = 1;
+    double corrQQx = 1.0;
+    double corrQQy = 1.0;
+    double corrQQ = 1.0 ;
+    double evPlaneRes = 1.0;
 
     // Get QQ-correlations from CCDB
     if (cfgCCDBdir_QQ.value.empty() == false) {
@@ -200,7 +200,6 @@ struct SpectatorPlaneTutorial {
     }
 
     for (const auto& track : tracks) {
-
       // constrain angle to 0 -> [0,0+2pi]
       auto phi = RecoDecay::constrainAngle(track.phi(), 0);
 
@@ -210,6 +209,9 @@ struct SpectatorPlaneTutorial {
       double uxMH = std::cos(2 * phi);
       double uyMH = std::sin(2 * phi);
 
+      if(corrQQ==0 || corrQQx==0 || corrQQy==0)
+        continue;
+
       double v1A = (uy * qyA + ux * qxA) / std::sqrt(std::fabs(corrQQ));
       double v1C = (uy * qyC + ux * qxC) / std::sqrt(std::fabs(corrQQ));