class HParticleT0Reco: public HReconstructor

 This function fills the particle buffers for TOF and RPC.
 It suppose that the meta matching was already recalculated!!!!!
 It uses tracks with some minimum quality criteria.
 Also it does a mild cut on vertex : dist>(1./mom)*1200.+10.

 Function which loops over candidates onces corrected by T0 and sets the best hypo just to the best traces.

 hist are store as arrays:

    void readHistWalk(){

    TFile* fileWC = new TFile("walk_par_histos3.root","READ");

    TH1D* h [2][16];
    TH1D* he[2][16];

    TArrayD  astartdt [2][16];
    TArrayD  astartsdt[2][16];

    TArrayD  astartdtback [2][16];
    TArrayD  astartsdtback[2][16];

    for(Int_t i=0;i<2;i++) {
	for(Int_t j=0;j<16;j++) {
	    h [i][j] = (TH1D*)fileWC->Get(Form("h_q_maxdt_mod%i_strip%i",i,j));
	    he[i][j] = (TH1D*)fileWC->Get(Form("h_q_sdt_mod%i_strip%i"  ,i,j));
	    HHistConverter::fillArray(h [i][j],astartdt [i][j],Form("h_q_maxdt_mod%i_strip%i",i,j),10,10, kFALSE);
	    HHistConverter::fillArray(he[i][j],astartsdt[i][j],Form("h_q_sdt_mod%i_strip%i"  ,i,j),10,10, kFALSE);
	}
    }




    ofstream out;
    out.open("myhists.txt");

    for(Int_t i=0;i<2;i++) {
	for(Int_t j=0;j<16;j++) {
	    HHistConverter::writeArray(out,h [i][j]->GetName(),astartdt [i][j],10);
	}
    }
    for(Int_t i=0;i<2;i++) {
	for(Int_t j=0;j<16;j++) {
	    HHistConverter::writeArray(out,he[i][j]->GetName(),astartsdt[i][j],10);
	}
    }
    out.close();


    out.open("myhists_array.txt");

    Int_t imax = 2;
    Int_t jmax = 16;
    Int_t kmax = astartdt[0][0].GetSize();

    out<<Form("Double_t astartdt[%i][%i][%i]={ ",imax,jmax,kmax);
    for(Int_t i=0;i<imax;i++) {
	out<<"{ ";
	for(Int_t j=0;j<jmax;j++) {
	    if(j>0)out<<endl;
	    out<<"{ ";
	    for(Int_t k=0;k<kmax;k++) {
		if(k<kmax-1) out<<astartdt[i][j][k]<<" , ";
       	else         out<<astartdt[i][j][k];
		if(k>0&&k<kmax&&k%10==0) out<<endl;
	    }
	    if(j<jmax-1) out<<"} , ";
            else         out<<"} ";
	}
	if(i<imax-1)out<<"} , "<<endl;
        else        out<<"} ";
    }
    out<<"};"<<endl;

    out<<endl;

    out<<Form("Double_t astartsdt[%i][%i][%i]={ ",imax,jmax,kmax);
    for(Int_t i=0;i<imax;i++) {
	out<<"{ ";
	for(Int_t j=0;j<jmax;j++) {
	    if(j>0)out<<endl;
	    out<<"{ ";
	    for(Int_t k=0;k<kmax;k++) {
		if(k<kmax-1) out<<astartsdt[i][j][k]<<" , ";
		else         out<<astartsdt[i][j][k];
       	if(k>0&&k<kmax&&k%10==0) out<<endl;
	    }
	    if(j<jmax-1) out<<"} , ";
            else         out<<"} ";
	}
	if(i<imax-1)out<<"} , "<<endl;
        else        out<<"} ";
    }
    out<<"};"<<endl;

    out.close();


}