for(int ii=0;ii<mcarray->GetEntries();ii++){
AliAODMCParticle *mctrack = (AliAODMCParticle*) mcarray->At(ii);
if(mctrack->GetPdgCode()!=13) continue;
for(int ii=0;ii<mcarray->GetEntries();ii++){
AliAODMCParticle *mctrack = (AliAODMCParticle*) mcarray->At(ii);
if(mctrack->GetPdgCode()!=13) continue;