message: Fix theta* computation
diff: @@ -296,7 +296,7 @@ void AliAnalysisTaskSEDstarPolarization::UserExec(Option_t * /*option*/)
  double cosThetaStarProd = TMath::Abs(normalVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  double cosThetaStarHelicity = TMath::Abs(helicityVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  double cosThetaStarBeam = TMath::Abs(beamVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  - double thetaStarBeam = TMath::ACos(cosThetaStarBeam);
  + double thetaStarBeam = TMath::ACos(beamVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  double phiStarBeam = TMath::ATan2(threeVecPiCM.Y(), threeVecPiCM.X());
  std::vector<double> var4nSparse = {mass, ptCand, cosThetaStarBeam, cosThetaStarProd, cosThetaStarHelicity, centrality};
  @@ -491,7 +491,7 @@ void AliAnalysisTaskSEDstarPolarization::FillMCGenAccHistos(TClonesArray *arrayM
  double cosThetaStarProd = TMath::Abs(normalVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  double cosThetaStarHelicity = TMath::Abs(helicityVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  double cosThetaStarBeam = TMath::Abs(beamVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  - double thetaStarBeam = TMath::ACos(cosThetaStarBeam);
  + double thetaStarBeam = TMath::ACos(beamVec.Dot(threeVecPiCM) / TMath::Sqrt(threeVecPiCM.Mag2()));
  double phiStarBeam = TMath::ATan2(threeVecPiCM.Y(), threeVecPiCM.X());
  if (orig == 4 && !isParticleFromOutOfBunchPileUpEvent)
  