message: select only particles with status 1 in primary isolation cone loop
diff: @@ -7914,7 +7914,7 @@ void AliAnaParticleIsolation::FillAcceptanceHistograms()
  // Consider only final state particles, but this depends on generator,
  // status 1 is the usual one, in case of not being ok, leave the possibility
  // to not consider this.
  - if( partInConeStatus > 1 &&
  + if( partInConeStatus != 1 &&
  GetMCAnalysisUtils()->GetMCGenerator()!= AliMCAnalysisUtils::kBoxLike ) continue ;
  // Protection against floating point exception
  @@ -7966,8 +7966,9 @@ void AliAnaParticleIsolation::FillAcceptanceHistograms()
  }
  }
  - //printf("status %d, E %f, pt %f, pz %f, pdg %d, gamma eta %f, phi %f, par in cone eta %f phi %f\n",
  - // partInConeStatus, partInConeE, partInConePt, mcisopAOD->Pz(), partInConePDG, photonEta, photonPhi, partInConeEta, partInConePhi);
  +// printf("Input status %d, E %f, pt %f, pz %f, pdg %d, gamma eta %f, phi %f, par in cone eta %f phi %f\n",
  +// partInConeStatus, partInConeE, partInConePt, mcisop->Pz(), partInConePDG,
  +// photonEta, photonPhi, partInConeEta, partInConePhi);
  dR = GetIsolationCut()->Radius(photonEta, photonPhi, partInConeEta, partInConePhi);
  @@ -7992,12 +7993,15 @@ void AliAnaParticleIsolation::FillAcceptanceHistograms()
  //printf("Selected particles pdg %d, status %d\n", partInConePDG, partInConeStatus);
  }
  -
  if ( !physPrimary ) continue ;
  sumPtInCone += partInConePt;
  if(partInConePt > GetIsolationCut()->GetPtThreshold() &&
  partInConePt < GetIsolationCut()->GetPtThresholdMax()) npart++;
  +
  +// printf("Selected: status %d, E %f, pt %f, pz %f, pdg %d, eta %f, phi %f, par in cone eta %f phi %f\n",
  +// partInConeStatus, partInConeE, partInConePt, mcisop->Pz(), partInConePDG,
  +// photonEta, photonPhi, partInConeEta, partInConePhi);
  }
  ///////END ISO MC/////////////////////////
  @@ -8161,6 +8165,8 @@ void AliAnaParticleIsolation::FillAcceptanceHistograms()
  sumPtInCone > GetIsolationCut()->GetSumPtThresholdMax()))
  isolated = kTRUE;
  + //printf("Primary in cone n Part %d sum pT in cone %2.2f\n",npart, sumPtInCone);
  +
  if(GetIsolationCut()->GetICMethod() == AliIsolationCut::kPtThresIC &&
  npart == 0)
  isolated = kTRUE;
  