break;
}
}
+
+
+ if(!aod) {delete [] pdgdaughters;return;}
+
Bool_t isSimpleMode=fSimpleMode;
if(!arrayProng) {
AliInfo("Branch not found! The output will contain only trak related histograms\n");
return;
}
}
- if(!aod) {delete [] pdgdaughters;return;}
// fix for temporary bug in ESDfilter
// the AODs with null vertex pointer didn't pass the PhysSel
- if(!aod->GetPrimaryVertex() || TMath::Abs(aod->GetMagneticField())<0.001) return;
+ if(!aod->GetPrimaryVertex() || TMath::Abs(aod->GetMagneticField())<0.001) {
+ delete [] pdgdaughters;
+ return;
+ }
// count event
fNEntries->Fill(0);
AliAODMCParticle* mcPartDaughterD0 = dynamic_cast<AliAODMCParticle*>(fmcArray->At(daughter0ds));
AliAODMCParticle* mcPartDaughterPis = dynamic_cast<AliAODMCParticle*>(fmcArray->At(daughter1ds));
+ if(!mcPartDaughterD0 || !mcPartDaughterPis) return kFALSE;
+
Double_t vtx1[3] = {0,0,0}; // primary vertex
Double_t vtx2daughter0[3] = {0,0,0}; // secondary vertex from daughter 0
Double_t vtx2daughter1[3] = {0,0,0}; // secondary vertex from daughter 1
AliAODMCParticle* mcPartDaughter0 = dynamic_cast<AliAODMCParticle*>(fmcArray->At(daughter0)); //D0
AliAODMCParticle* mcPartDaughter1 = dynamic_cast<AliAODMCParticle*>(fmcArray->At(daughter1)); //pis
+ if(!mcPartDaughter0 || !mcPartDaughter1) return kFALSE;
+
// getting vertex from daughters
mcPartDaughter0->XvYvZv(vtx2daughter0); // cm
mcPartDaughter1->XvYvZv(vtx2daughter1); //cm