- j++;
-
- // debug info to check fro charm daugthers
- if((TMath::Abs(part->GetPdgCode()))==411){
- iCharm++;
- Int_t d0 = mcpart->GetFirstDaughter();
- Int_t d1 = mcpart->GetLastDaughter();
- if(d0>0&&d1>0){
- for(int id = d0;id <= d1;id++){
- iAll++;
- if(mcH->IsParticleSelected(id))iTaken++;
+ j++;
+ }
+ }
+
+ /*
+
+ // check varibales for charm need all daughters
+ static int iTaken = 0;
+ static int iAll = 0;
+ static int iCharm = 0;
+
+ for (Int_t ip = 0; ip < np; ip++){
+ AliMCParticle* mcpart = (AliMCParticle*) mcE->GetTrack(ip);
+ TParticle* part = mcpart->Particle();
+ // debug info to check fro charm daugthers
+ if((TMath::Abs(part->GetPdgCode()))==411){
+ iCharm++;
+ Int_t d0 = part->GetFirstDaughter();
+ Int_t d1 = part->GetLastDaughter();
+ if(d0>0&&d1>0){
+ for(int id = d0;id <= d1;id++){
+ iAll++;
+ if(mcH->IsParticleSelected(id)){
+ iTaken++;
+ Printf("> %d/%d Taken",id,nprim);
+ PrintMCParticle( (AliMCParticle*)mcE->GetTrack(id),id);
+ }
+ else{
+ Printf("> %d/%d NOT Taken",id,nprim);
+ PrintMCParticle( (AliMCParticle*)mcE->GetTrack(id),id);