vis_mc_ids_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);

  TClonesArray* hit_array=new TClonesArray("PndEmcHit");		
  abc->SetBranchAddress("EmcHit",&hit_array);

//   TClonesArray* digi_array=new TClonesArray("PndEmcDigi");		
//   abc->SetBranchAddress("EmcDigi",&digi_array);

  TClonesArray* cluster_array=new TClonesArray("PndEmcCluster");		
  abc->SetBranchAddress("EmcCluster",&cluster_array);

  TCanvas* c1 = new TCanvas("c1", "c1", 100, 100, 800, 800);
  TLatex *lt;
  lt=new TLatex();
  lt->SetTextFont(32);
  lt->SetTextAlign(22);

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

    cout << "Hits: "<<endl;
    lt->SetTextColor(1);
    lt->SetTextSize(0.015);
    Int_t nhit = hit_array->GetEntries();
    for (Int_t i=0; i<nhit; i++)
    {
	    PndEmcHit *hit;
	    hit=(PndEmcHit*)hit_array->At(i);
//         cout <<"Hit "<<i<<" energy: "<<hit->GetEnergy()<<endl;
    	Double_t x,y;
    	x=hit->GetX();
    	y=hit->GetY();
    	TLine *l;
    	Double_t fak;
    	fak=4.1e-3+1.5e-3*hit->GetEnergy();
	    l=new TLine(x*4e-3+.5,y*4e-3+.5,x*fak+.5,y*fak+.5);
	    l->SetLineWidth(5);
	    l->SetLineColor(6);
	    l->Draw();

	    TString mctxt;
	    mctxt="";
	    if( hit->GetMcList().size()==0) cout << "Zero List??"<<endl;
		for(Int_t j=0; j<hit->GetMcList().size(); j++){
			Int_t id;
			id=hit->GetMcList()[j];
			propagate_mclist(id, mctrack_array,x,y);
			if(j!=0) mctxt+=",";
			add_mccode(mctxt,((PndMCTrack *)mctrack_array->At(id))->GetPdgCode());
//                        cout <<"MCIndex: "<<id<<": "<<((PndMCTrack *)mctrack_array->At(id))->Get4Momentum().E()<<endl;
 		}

	    lt->DrawLatex(x*(fak+.8e-3)+.5,y*(fak+.8e-3)+.5,mctxt.Data());
    }

    cout << "Clusters: "<<endl;
    lt->SetTextColor(2);
    lt->SetTextSize(0.03);
    Int_t nclu = cluster_array->GetEntries();
    for (Int_t i=0; i<nclu; i++)
    {
    	PndEmcCluster *clu;
    	clu=(PndEmcCluster*)cluster_array->At(i);
	cout <<"Cluster "<<i<<" energy: "<<clu->energy()<< " , "<< clu->GetEnergyCorrected()<<endl;
    	Double_t x,y;
    	x=clu->x();
    	y=clu->y();
    	TLine *l;
    	Double_t fak;
    	fak=4.1e-3+1.5e-3*clu->energy();
	    l=new TLine(x*4e-3+.5,y*4e-3+.5,x*fak+.5,y*fak+.5);
	    l->SetLineWidth(9);
	    l->Draw();

	    TString mctxt;
	    mctxt="";
	    if( clu->GetMcList().size()==0) cout << "Zero List??"<<endl;
		for(Int_t j=0; j<clu->GetMcList().size(); j++){
			Int_t id;
			id=clu->GetMcList()[j];
			propagate_mclist(id, mctrack_array,x,y);
			if(j!=0) mctxt+=",";
			add_mccode(mctxt,((PndMCTrack *)mctrack_array->At(id))->GetPdgCode());
			cout <<"MCIndex: "<<id<<": "<<((PndMCTrack *)mctrack_array->At(id))->Get4Momentum().E()<<endl;
		}
	    lt->DrawLatex(x*(fak+.8e-3)+.5,y*(fak+.8e-3)+.5,mctxt.Data());
    }


    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)
{
	PndMCTrack *mct;
	PndMCTrack *mo;
	TLine *l;
	TVector3 v;
	Double_t x1,y1,x2,y2;

	if(id<0) return;
	mct=(PndMCTrack *)mctrack_array->At(id);
	v=mct->GetStartVertex();
	x1=v.x();
	y1=v.y();

    l=new TLine(x1*4e-3+.5,y1*4e-3+.5,x*4e-3+.5,y*4e-3+.5);
	l->SetLineWidth(1);
	l->SetLineColor(3);
    l->Draw();

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

	    cout <<"Loop Start: "<<id;
	Bool_t loop=false;
	while(mct){
		int mid;
		mid=mct->GetMotherID();
		cout <<" ("<< mct->GetPdgCode() <<")";
		if(mid<0){
		  cout <<"*";
		  mid=mct->GetSecondMotherID();
		}
		cout <<"=> "<<mid;
		if(mid<0) break;

		mo=((PndMCTrack *))mctrack_array->At(mid);

		v=mct->GetStartVertex();
		x1=v.x();
		y1=v.y();
		v=mo->GetStartVertex();
		x2=v.x();
		y2=v.y();

	    l=new TLine(x1*4e-3+.5,y1*4e-3+.5,x2*4e-3+.5,y2*4e-3+.5);
	    if(loop){
	    	l->SetLineWidth(1);
	    	l->SetLineColor(7);
	    }else{
	    	l->SetLineWidth(3);
	    	l->SetLineColor(2);
	    }
	    l->Draw();
	    
	    mct=mo;
	    loop=true;
	}
	cout <<endl;
}