message: Fix D0 MC selection
diff: @@ -467,6 +467,10 @@ void AliAnalysisTaskSEDstarPolarization::UserExec(Option_t * /*option*/)
  else
  {
  if(labD > 0) {
  + //check if reflected signal
  + int labDauFirst = dauPi->GetLabel();
  + AliAODMCParticle* dauFirst = dynamic_cast<AliAODMCParticle *>(arrayMC->At(labDauFirst));
  + if(TMath::Abs(dauFirst->GetPdgCode()) == 211) {
  if(orig == 4) {
  fnSparseReco[1]->Fill(var4nSparse.data());
  fnSparseRecoThetaPhiStar[1]->Fill(var4nSparseThetaPhiStar.data());
  @@ -476,6 +480,7 @@ void AliAnalysisTaskSEDstarPolarization::UserExec(Option_t * /*option*/)
  fnSparseRecoThetaPhiStar[2]->Fill(var4nSparseThetaPhiStar.data());
  }
  }
  + }
  else {
  fnSparseReco[3]->Fill(var4nSparse.data());
  fnSparseRecoThetaPhiStar[3]->Fill(var4nSparseThetaPhiStar.data());
  @@ -517,6 +522,10 @@ void AliAnalysisTaskSEDstarPolarization::UserExec(Option_t * /*option*/)
  else
  {
  if(labD > 0) {
  + //check if reflected signal
  + int labDauFirst = dauPi->GetLabel();
  + AliAODMCParticle* dauFirst = dynamic_cast<AliAODMCParticle *>(arrayMC->At(labDauFirst));
  + if(TMath::Abs(dauFirst->GetPdgCode()) == 211) {
  if(orig == 4) {
  fnSparseReco[1]->Fill(var4nSparse.data());
  fnSparseRecoThetaPhiStar[1]->Fill(var4nSparseThetaPhiStar.data());
  @@ -526,6 +535,7 @@ void AliAnalysisTaskSEDstarPolarization::UserExec(Option_t * /*option*/)
  fnSparseRecoThetaPhiStar[2]->Fill(var4nSparseThetaPhiStar.data());
  }
  }
  + }
  else {
  fnSparseReco[3]->Fill(var4nSparse.data());
  fnSparseRecoThetaPhiStar[3]->Fill(var4nSparseThetaPhiStar.data());
  @@ -731,11 +741,14 @@ void AliAnalysisTaskSEDstarPolarization::FillMCGenAccHistos(TClonesArray *arrayM
  bool isFidAcc = false;
  bool isDaugInAcc = false;
  int nDau = 3;
  + if(fDecChannel == kDstartoD0pi)
  deca = AliVertexingHFUtils::CheckDstarDecay(arrayMC, mcPart, labDau);
  + else
  + deca = AliVertexingHFUtils::CheckD0Decay(arrayMC, mcPart, labDau);
  if (labDau[0] == -1)
  continue; //protection against unfilled array of labels
  - if (deca > 0)
  + if (deca == 1)
  isGoodDecay = true;
  if (isGoodDecay)
  