check_mc_tree_sums()
{
  gROOT->LoadMacro("$VMCWORKDIR/macro/mvd/Tools.C");
  gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
  gROOT->Macro("~/rootlogon.C");

  TTree *abc=cbmsim;

    TClonesArray* mctrack_array=new TClonesArray("PndMCTrack");
  abc->SetBranchAddress("MCTrack",&mctrack_array);

  Int_t ncounts = abc->GetEntries();
  for(int k = 0; k < 2/*ncounts*/; k++){
    abc->GetEntry(k);

    cout << "============="<<endl;
    for(Int_t i=0; i<mctrack_array->GetEntries(); i++){
      mo=(PndMCTrack *)mctrack_array->At(i);
      if(mo->GetMotherID()==-1 && mo->GetSecondMotherID()==-1) propagate_mclist(i, mctrack_array);
    }
    cout << "============="<<endl;
   }
}

void add_mccode(TString &mctxt,Int_t id)
{
	switch(id)
	{
	case 22: mctxt+="#gamma" ; break;
	case 11: mctxt+="e^{-}" ; break;
	case -11: mctxt+="e^{+}" ; break;
	case 13: mctxt+="#mu^{-}" ; break;
	case -13: mctxt+="#mu^{+}" ; break;
	case 111: mctxt+="#pi^{0}" ; break;
	case 211: mctxt+="#pi^{+}" ; break;
	case -211: mctxt+="#pi^{-}" ; break;
	case 310: mctxt+="K^{0}_{S}" ; break;
	case 130: mctxt+="K^{0}_{L}" ; break;
	case 321: mctxt+="K^{+}" ; break;
	case -321: mctxt+="K^{-}" ; break;

	case 12: mctxt+="#nu_{e}"; break; 
	case 14: mctxt+="#nu_{#mu}"; break; 
	case -12: mctxt+="#bar{#nu}_{e}"; break; 
	case -14: mctxt+="#bar{#nu}_{#mu}"; break; 

	case 2112: mctxt+="n" ; break;
	case 2212: mctxt+="p" ; break;

	case 3122: mctxt+="#Lambda"; break; 
	case 3124: mctxt+="#Lambda(1520)"; break; 
	case 13122: mctxt+="#Lambda(1405)"; break; 

	case 3222: mctxt+="#Sigma^{+}"; break; 
	case 3224: mctxt+="#Sigma^{*+}"; break; 
	case 3212: mctxt+="#Sigma^{0}"; break; 
	case 3214: mctxt+="#Sigma^{*0}"; break; 
	case 3112: mctxt+="#Sigma^{-}"; break; 
	case 3114: mctxt+="#Sigma^{*-}"; break; 
	
	case 3312: mctxt+="#Xi^{-}"; break; 
	case 3322: mctxt+="#Xi^{0}"; break; 
	case 3314: mctxt+="#Xi^{*-}"; break; 
	case 3324: mctxt+="#Xi^{*0}"; break; 
	
	case 3334: mctxt+="#Omega^{-}"; break; 
	
	case -2112: mctxt+="#bar{n}" ; break;
	case -2212: mctxt+="#bar{p}" ; break;

	case -3122: mctxt+="#bar{#Lambda}"; break; 
	case -3124: mctxt+="#bar{#Lambda}(1520)"; break; 
	case -13122: mctxt+="#bar{#Lambda}(1405)"; break; 
	
	case -3222: mctxt+="#bar{#Sigma}^{-}"; break; 
	case -3224: mctxt+="#bar{#Sigma}^{*-}"; break; 
	case -3212: mctxt+="#bar{#Sigma}^{0}"; break; 
	case -3214: mctxt+="#bar{#Sigma}^{*0}"; break; 
	case -3112: mctxt+="#bar{#Sigma}^{+}"; break; 
	case -3114: mctxt+="#bar{#Sigma}^{*+}"; break; 

	case -3312: mctxt+="#bar{#Xi}^{+}"; break; 
	case -3322: mctxt+="#bar{#Xi}^{0}"; break; 
	case -3314: mctxt+="#bar{#Xi}^{*+}"; break; 
	case -3324: mctxt+="#bar{#Xi}^{*0}"; break; 
	
	case -3334: mctxt+="#bar{#Omega}^{+}"; break; 
	
	case 443: mctxt+="J/#Psi"; break; 
	case 100443: mctxt+="#Psi(2S)"; break; 
	case 441: mctxt+="#eta_c(1S)"; break; 
        case 10441: mctxt+="#xi_{c0}(1P)"; break; 
        case 20443: mctxt+="#xi_{c1}(1P)"; break; 
        case 445: mctxt+="#xi_{c2}(1P)"; break; 
        case 100441: mctxt+="#eta_c(2S)"; break; 
	case 10443: mctxt+="h_c(1P)"; break; 

	default: mctxt+="?" ; break;
	}
}


TString propagate_mclist(Int_t id,TClonesArray * mctrack_array)
{
    PndMCTrack *mct;
    PndMCTrack *mo;
    
    TLorentzVector sum, tmp;
    TString chtxt="";
    
    double fak=0.3;

    if(id<0) return "none";
    mct=(PndMCTrack *)mctrack_array->At(id);
    sum=mct->Get4Momentum();
    tmp=sum;

    if(mct->GetMotherID()==-1){// show only first geant generation
      for(Int_t i=0; i<mctrack_array->GetEntries(); i++){
	mo=(PndMCTrack *)mctrack_array->At(i);
	
	if(mo->GetMotherID()==id){
	   sum-=mo->Get4Momentum();
	  chtxt+=propagate_mclist( i, mctrack_array);
	}else if(mo->GetMotherID()==-1 && mo->GetSecondMotherID()==id){
	   sum-=mo->Get4Momentum();
	  chtxt+=propagate_mclist( i, mctrack_array);
	}
      }
    }

    TString mctxt;
    mctxt="";
    add_mccode(mctxt,mct->GetPdgCode());

      if( sum.P()>1e-4 || sum.E()>1e-4)
    cout << "id "<<id<<": "<<mctxt<<" : "<<
      tmp.Px()<<" , "<<
      tmp.Py()<<" , "<<
      tmp.Pz()<<" , "<<
      tmp.E()<<" , "<<
      sum.Px()<<" , "<<
      sum.Py()<<" , "<<
      sum.Pz()<<" , "<<
      sum.E()<<" , "<<
     chtxt<< endl;
   return (mctxt+" ");
}