#ifndef ___CONVERT_ROOT___
#define ___CONVERT_ROOT___


//############################ CONVERTERS ################################
Int_t convert_ROOT(objects& obj)
{
    // converts SMASH root file to HGeant ascii file.
    // Its assumed that the reaction is in lab frame already (see the smash settings)
    // K0_S and K0_L are equaly distributed 50:50

    TFile*  f = new TFile(obj.input.Data(),"READ");

    if(!f){
       cerr<<"ERROR: Could not open infile : "<<obj.input.Data()<<"!"<<endl;
       return 1;
    }
    f->cd();

    TTree* T = (TTree*) gDirectory->Get("particles");
    if(!T){
	cerr<<"ERROR: Could not reed TTee \"particles\" from infile : "<<obj.input.Data()<<"!"<<endl;
	return 1;
    }

    obj.evt.flag = 1;

#define NPAR 1000

    //Declaration of leaves types
    Int_t           npart;
    Double_t        impact_b;
    Int_t           ev;
    Int_t           tcounter;
    Int_t           pdgcode[NPAR];
    Double_t        p0[NPAR];
    Double_t        px[NPAR];
    Double_t        py[NPAR];
    Double_t        pz[NPAR];
    Double_t        t [NPAR];
    Double_t        x [NPAR];
    Double_t        y [NPAR];
    Double_t        z [NPAR];

    // Set branch addresses.
    T->SetBranchAddress("npart",&npart);
    T->SetBranchAddress("impact_b",&impact_b);
    T->SetBranchAddress("ev",&ev);
    T->SetBranchAddress("tcounter",&tcounter);
    T->SetBranchAddress("pdgcode",pdgcode);
    T->SetBranchAddress("p0",p0);
    T->SetBranchAddress("px",px);
    T->SetBranchAddress("py",py);
    T->SetBranchAddress("pz",pz);
    T->SetBranchAddress("t",t);
    T->SetBranchAddress("x",x);
    T->SetBranchAddress("y",y);
    T->SetBranchAddress("z",z);


    Long64_t nentries = T->GetEntries();

    if(nentries > obj.nEvts && obj.nEvts != -1) nentries = obj.nEvts;


    Long64_t nbytes = 0;
    smash_p p;
    Int_t evtCt     = 0;
    for (Long64_t i=0; i<nentries;i++) {
	nbytes = T->GetEntry(i);
        if(nbytes<=0) break;

	obj.evt.evtNr     = ev+1;
	obj.evt.nParticle = npart;
	obj.evt.impact    = impact_b;

	//if(obj.debug){ cout<<"Event : "<<setw(5)<<obj.evt.evtNr<<"  npart "<<setw(5)<<obj.evt.nParticle <<" impact "<<setw(10)<<obj.evt.impact<<endl;}

	for(Int_t j = 0; j < npart; j++){
	    p.t      = (Int_t) t[j];
	    p.x      = x[j];
	    p.y      = y[j];
	    p.z      = z[j];
	    p.p0     = p0[j];
	    p.px     = px[j];
	    p.py     = py[j];
	    p.pz     = pz[j];
	    p.pdg    = (Int_t) pdgcode[j];

	    mapPdgPIDtoGeant3(obj.PDG,p); // mass + charge + weight + geantID of p are set

	    //if(obj.debug){cout<<setw(4)<<j<<"    : "<<p.t<< " "<<setw(12)<<p.x<<" "<<setw(12)<<p.y<<" "<<setw(12)<<p.z<<" "<<setw(12)<<p.p0<<" "<<setw(12)<<p.px<<" "<<setw(12)<<p.py<<" "<<setw(12)<<p.pz<<" "<<setw(6)<<p.pdg<<endl;}

	    p.ID     = -1;

	    p.v.SetXYZM(p.px,p.py,p.pz,p.mass);

	    p.E = p.p0;
            p.state = 1;
	    obj.particles.push_back(p);
	}

	if(obj.debug) obj.printEvent("after event processing (no boost)");

	if(params.doBoost){
	    obj.boost();
	}

#ifdef WITHPLUTO
	//-----------------------------------------------------------
	// loop over event and fill PLUTO tree

	if(params.doTree == 1)
	{
	    fillPutoTree(obj,0.);
	}
	//-----------------------------------------------------------
#endif
	write_Header(obj);
	write_Event (obj);
	evtCt++;
	if(evtCt == obj.nEvts) {
	    break;
	}
    }

    obj.out.close();
    cout<<"Number of events read from file : "<<setw(6)<<evtCt<<" (requested "<<obj.nEvts<<")"<<endl;

    return 0;
}

#endif