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

//  TTree *sim1, *sim2;
//  TFile *_file0 = TFile::Open("full_emc_test_gamma.root");
//  sim1=cbmsim;
//  sim1->SetName("sim1");
  TChain *abc=new TChain("cbmsim");
  abc->Add("emc_complete.root");
//  abc->Add("sim_emc_test_gamma.root");
//  abc->AddFriend("sim1");

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

  TCanvas* c1 = new TCanvas("c1", "c1", 100, 100, 800, 800);

  Int_t ncounts = abc->GetEntries();
  for(int k = 0; k < ncounts; k++){
    abc->GetEntry(k);
    c1->Clear();

    TLatex *lt;
    lt=new TLatex();
    lt->SetTextFont(32);
    lt->SetTextAlign(22);
    lt->SetTextSize(0.03);

    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,0,0,lt);
    }

    TString bild="bild.ps[";
    c1->Print(bild);
   }
   c1->Print("bild.ps]");
}

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;
	}
}


void propagate_mclist(Int_t id,TClonesArray * mctrack_array,Double_t x,Double_t y,TLatex *lt)
{
    PndMCTrack *mct;
    PndMCTrack *mo;
    TLine *l;
    TVector3 v;
    Double_t x1,y1;
    
    double fak=0.3;

    if(id<0) return;
    mct=(PndMCTrack *)mctrack_array->At(id);
    x1=mct->Get4Momentum().Px();
    y1=mct->Get4Momentum().Py();
    x1/=mct->Get4Momentum().E();
    y1/=mct->Get4Momentum().E();
    x1+=x;
    y1+=y;

    
    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){
	  propagate_mclist( i, mctrack_array, x1, y1,lt);
	}else if(mo->GetMotherID()==-1 && mo->GetSecondMotherID()==id){
	  propagate_mclist( i, mctrack_array, x1, y1,lt);
	}
      }
    }

    l=new TLine(x1*fak+.5,y1*fak+.5,x*fak+.5,y*fak+.5);
    l->SetLineWidth(1);

    if( mct->IsGeneratorCreated()){
      if(mct->IsGeneratorDecayed()){
	l->SetLineColor(4);
      }else{
	l->SetLineColor(2);
//    This is identical to 	
//    if( mct->IsGeneratorLast()){
//	l->SetLineColor(6);
//    }
      }
    }else{
      l->SetLineColor(3);
    }
    
    l->Draw();

    TString mctxt;
    mctxt="";
    add_mccode(mctxt,mct->GetPdgCode());
    lt->DrawLatex((x1+x)/2*fak+.5,(y1+y)/2*fak+.5,mctxt.Data());    
}