/**
 * \file TrkParams.cpp
 * \author Elena Vannuccini
 */

#include <TrkParams.h>
#include <iostream>
using namespace std;


// -------------------------
// CALIBRATION (PED,SIG,BAD)
// -------------------------
GL_TRK_CALIB TrkParams::calib      = GL_TRK_CALIB();
TString      TrkParams::calibpatht = "";
TString      TrkParams::calibpathf = "";
Bool_t       TrkParams::calib104   = 0;
Bool_t       TrkParams::calibload  = false;
UInt_t       TrkParams::trkparamtype[] = {1,2,3,4,5,6,7};
GL_PARAM     TrkParams::gl[]           = {GL_PARAM(),GL_PARAM(),GL_PARAM(),GL_PARAM(),GL_PARAM(),GL_PARAM(),GL_PARAM()};
TString      TrkParams::glpath[]       = {"","","","","","",""};
Bool_t       TrkParams::glload[]       = {false,false,false,false,false,false,false}; 

TrkParams::TrkParams(){

    cout<<"TrkParams::TrkParams()"<<endl;

}
float aint(float f){
    int   ret  = (int)f;
    float ret1 = 0;
    ret1=(float)ret;
    return ret1;	
}

float anint(float f){
    int   ret  = (int)(f+0.5);
    float ret1 = 0;
    ret1=(float)ret;
    return ret1;	
}
/**
 * Static method to set the calibration to be loaded. 

 * It retrieves the calibration associated to a given run, from 
 * a given DB. It fills the F77 common.
 * @param glrun Pointer to a GL_RUN row 
 * @param dbc   Pointer to a DB connection 
 */
Bool_t TrkParams::SetCalib( GL_RUN *glrun, TSQLServer* dbc ){

    GL_TRK_CALIB q2;
    GL_ROOT q3;
    GL_PARAM q4;

    if( q2.Query_GL_TRK_CALIB(glrun->RUNHEADER_TIME,dbc) )return false;
    if(q2.EV_ROOT_CALIBTRK1 != q2.EV_ROOT_CALIBTRK2) 
	printf("WARNING!! ---> EV_ROOT_CALIBTRK1=%d it's different from EV_ROOT_CALIBTRK2=%d \n\n",q2.EV_ROOT_CALIBTRK1,q2.EV_ROOT_CALIBTRK2);	
    if( q3.Query_GL_ROOT(q2.ID_ROOT_L0,dbc) )return false;

    calib      = q2;
    calibpathf = q3.PATH + q3.NAME;
    calib104   = (glrun->TRK_CALIB_USED==104);
    calibload  = kFALSE;
    if(calib104){
	if (q4.Query_GL_PARAM(glrun->RUNHEADER_TIME,7,dbc) )return false; 
	calibpatht = q4.PATH + q4.NAME;
    }


//    return TrkParams::LoadCalib();
    return true;
};

/**
 * Method to fill the tracker calibration-parameter struct from on-line calibrations
 */
void TrkParams::FillACalibFrom(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){

    cTrkCalib *ca = &pedsigbad_;

    TTree *tr1 = (TTree*)f0->Get("CalibTrk1");  
    TTree *tr2 = (TTree*)f0->Get("CalibTrk2"); 

    CalibTrk1Event *calibdata1 = 0;
    CalibTrk2Event *calibdata2 = 0;
    tr1->SetBranchAddress("CalibTrk1",&calibdata1);
    tr2->SetBranchAddress("CalibTrk2",&calibdata2);
    tr1->GetEntry(ev_reg1);
    tr2->GetEntry(ev_reg2); 
    Int_t dspn1(0),dspn2(0);
	
    for(Int_t i=0;i<6;i++){
	dspn1=calibdata1->DSPnumber[i]-1;
	dspn2=calibdata2->DSPnumber[i]-1;
	/*     cout<<"dspn1= "<<dspn1<<"   dspn2= "<<dspn2<<endl; */
	for(Int_t ii=0;ii<3072;ii++){
	    Int_t j=(ii/128);
	    Int_t jj=ii-j*128; 
	    ca->pedestal[jj][j][dspn1]=calibdata1->DSPped_par[i][ii];
	    ca->pedestal[jj][j][dspn2]=calibdata2->DSPped_par[i][ii];
	    ca->pedestal_t[jj][j][dspn1]=aint(calibdata1->DSPped_par[i][ii]);
	    ca->pedestal_t[jj][j][dspn2]=aint(calibdata2->DSPped_par[i][ii]);
	    ca->sigma[jj][j][dspn1]=calibdata1->DSPsig_par[i][ii];
	    ca->sigma[jj][j][dspn2]=calibdata2->DSPsig_par[i][ii];
	    ca->sigma_t[jj][j][dspn1]=anint(calibdata1->DSPsig_par[i][ii]);
	    ca->sigma_t[jj][j][dspn2]=anint(calibdata2->DSPsig_par[i][ii]);
	    //	      cout<<"sigma trk1= "<<sigma[jj][j][dspn1]<<endl;
	    //        cout<<"sigma_t trk1= "<<sigma_t[jj][j][dspn1]<<endl;
	    //        cout<<"sigma trk2= "<<sigma[jj][j][dspn2]<<endl;
	    //        cout<<"sigma_t trk2= "<<sigma_t[jj][j][dspn2]<<endl;
	
	    //
	    // invert the bad
	    //
	    if(calibdata1->DSPbad_par[i][ii]==0) 
		ca->bad[jj][j][dspn1]=1;
	    else if(calibdata1->DSPbad_par[i][ii]==1)
		ca->bad[jj][j][dspn1]=0;
	    if(calibdata2->DSPbad_par[i][ii]==0)
		ca->bad[jj][j][dspn2]=1;
	    else if(calibdata2->DSPbad_par[i][ii]==1)
		ca->bad[jj][j][dspn2]=0;
	};
    };
	
    tr1->Delete();
    tr2->Delete();
};
/**
 * Method to fill the tracker calibration-parameter struct from on-line calibrations
 */
void TrkParams::FillFCalibFrom(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){

    cTrkCalib *ca = &pedsigbad_;

    TTree *tr1 = (TTree*)f0->Get("CalibTrk1");  
    TTree *tr2 = (TTree*)f0->Get("CalibTrk2"); 

    CalibTrk1Event *calibdata1 = 0;
    CalibTrk2Event *calibdata2 = 0;
    tr1->SetBranchAddress("CalibTrk1",&calibdata1);
    tr2->SetBranchAddress("CalibTrk2",&calibdata2);
    tr1->GetEntry(ev_reg1);
    tr2->GetEntry(ev_reg2); 
    Int_t dspn1(0),dspn2(0);
	
    for(Int_t i=0;i<6;i++){
	dspn1=calibdata1->DSPnumber[i]-1;
	dspn2=calibdata2->DSPnumber[i]-1;
	/*     cout<<"dspn1= "<<dspn1<<"   dspn2= "<<dspn2<<endl; */
	for(Int_t ii=0;ii<3072;ii++){
	    Int_t j=(ii/128);
	    Int_t jj=ii-j*128; 
	    ca->pedestal[jj][j][dspn1]=calibdata1->DSPped_par[i][ii];
	    ca->pedestal[jj][j][dspn2]=calibdata2->DSPped_par[i][ii];
	    ca->sigma[jj][j][dspn1]=calibdata1->DSPsig_par[i][ii];
	    ca->sigma[jj][j][dspn2]=calibdata2->DSPsig_par[i][ii];
	    //	      cout<<"sigma trk1= "<<sigma[jj][j][dspn1]<<endl;
	    //        cout<<"sigma_t trk1= "<<sigma_t[jj][j][dspn1]<<endl;
	    //        cout<<"sigma trk2= "<<sigma[jj][j][dspn2]<<endl;
	    //        cout<<"sigma_t trk2= "<<sigma_t[jj][j][dspn2]<<endl;
	
	    //
	    // invert the bad
	    //
	    if(calibdata1->DSPbad_par[i][ii]==0)
		ca->bad[jj][j][dspn1]=1;
	    else if(calibdata1->DSPbad_par[i][ii]==1)
		ca->bad[jj][j][dspn1]=0;
	    if(calibdata2->DSPbad_par[i][ii]==0)
		ca->bad[jj][j][dspn2]=1;
	    else if(calibdata2->DSPbad_par[i][ii]==1)
		ca->bad[jj][j][dspn2]=0;
	};
    };
	
    tr1->Delete();
    tr2->Delete();
};
/**
 * Method to fill the tracker calibration-parameter struct from on-line calibrations
 */
void TrkParams::FillTCalibFrom(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){

    cTrkCalib *ca = &pedsigbad_;

    TTree *tr1 = (TTree*)f0->Get("CalibTrk1");  
    TTree *tr2 = (TTree*)f0->Get("CalibTrk2"); 

    CalibTrk1Event *calibdata1 = 0;
    CalibTrk2Event *calibdata2 = 0;
    tr1->SetBranchAddress("CalibTrk1",&calibdata1);
    tr2->SetBranchAddress("CalibTrk2",&calibdata2);
    tr1->GetEntry(ev_reg1);
    tr2->GetEntry(ev_reg2); 
    Int_t dspn1(0),dspn2(0);
	
    for(Int_t i=0;i<6;i++){
	dspn1=calibdata1->DSPnumber[i]-1;
	dspn2=calibdata2->DSPnumber[i]-1;
	/*     cout<<"dspn1= "<<dspn1<<"   dspn2= "<<dspn2<<endl; */
	for(Int_t ii=0;ii<3072;ii++){
	    Int_t j=(ii/128);
	    Int_t jj=ii-j*128; 
	    ca->pedestal_t[jj][j][dspn1]=aint(calibdata1->DSPped_par[i][ii]);
	    ca->pedestal_t[jj][j][dspn2]=aint(calibdata2->DSPped_par[i][ii]);
	    ca->sigma_t[jj][j][dspn1]=anint(calibdata1->DSPsig_par[i][ii]);
	    ca->sigma_t[jj][j][dspn2]=anint(calibdata2->DSPsig_par[i][ii]);	
	};
    };
	
    tr1->Delete();
    tr2->Delete();
};
/**
 * Method to fill the tracker calibration-parameter struct from default calibration
 */
void TrkParams::FillTCalibFrom(TString path){
    
    path_.FillWith(path);
    fillpedsigfromdefault_();
    if(path_.error) throw -216;

};
/**
 * Method to fill the tracker VA1-mask struct from calibration
 */
void TrkParams::FillMask(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){
    
    cTrkMask *ma = &mask_;

    TrkCalib *calib = new TrkCalib();

    for(Int_t i=0; i<2; i++){
	TTree *tr=0; 
	if(i==0){
	    //-------------------------
	    //read calibration packet 1
	    //-------------------------
	    tr = (TTree*)f0->Get("CalibTrk1");  
	    tr->SetBranchAddress("CalibTrk1",calib->GetPointerTo());
	    tr->GetEntry(ev_reg1);
	}
	if(i==1){
	    //-------------------------
	    //read calibration packet 2
	    //-------------------------
	    tr = (TTree*)f0->Get("CalibTrk2"); 
	    tr->SetBranchAddress("CalibTrk2",calib->GetPointerTo());
	    tr->GetEntry(ev_reg2); 
	}
	for(Int_t ipkt=0; ipkt<NPLANE; ipkt++){
	    for(Int_t ivk=0; ivk<NVK; ivk++){
		Int_t view   = calib->GetView(ipkt);
		Float_t sigm = calib->GetMean("SIG-BAD",ipkt,ivk+1);
		Float_t sigmin = 0.;
		if(view%2)sigmin = 6.5;
		else      sigmin = 2.5;
		Int_t mask = 1;
		if( sigm<sigmin )mask = 0;
	    
		ma->mask_vk_run[ivk][view-1] = mask;
	    
//		if(mask==0)cout<< " VIEW "<<view<<" VK "<<ivk+1<<" -- <SIG> "<<sigm<<" ---> MASKED!"<<endl;
	    }	
	}
	if(tr)tr->Delete();
    }

    if(calib)calib->Delete();

};

/**
 * Static method to load calibration. 
 * Calibration is loaded ONLY IF TrkParams::CalibLoaded()==kTRUE, which appens the 
 * first time the method is called (but it does not work becouse I didn't set any 
 * default calibration for the moment) or  after the method TrkParams::SetCalib(***) is called.
 */
Bool_t TrkParams::LoadCalib( ){

    if( TrkParams::CalibIsLoaded() )return false;

//    cTrkCalib *ca = new cTrkCalib; //&pedsigbad_;

//     extern cTrkCalib   pedsigbad_;
//     extern cTrkMask    mask_; 

    cout << "--------------------------------------"<<endl;
    cout << "Retrieve calibration (PED-SIG-BAD) ..."<<endl;
    // =============================================================
    // retrieve calibration file needed to reduce data 
    // =============================================================
    // if run OBT is > last calibration "expiration date"
    //  - search for new calibration packet
    //  - load calibration parameters (full + truncated)    
    cout << "Full pedestals for cluster finding:";
    cout << " >> Loading from LEVEL0 file: "<< calibpathf << endl;
    FileStat_t t;
    if( gSystem->GetPathInfo(calibpathf.Data(),t) )return false;
   
    TFile *f0_c = new TFile(calibpathf);
    if ( !f0_c ) return false;
    cout << " calibration entries "<< calib.EV_ROOT_CALIBTRK1 << " " << calib.EV_ROOT_CALIBTRK2; 
    cout << " (from time "<< calib.FROM_TIME <<" to time "<< calib.TO_TIME <<")"<<endl;
//    ca->FillACalibFrom(f0_c,calib.EV_ROOT_CALIBTRK1,calib.EV_ROOT_CALIBTRK2);
    TrkParams::FillACalibFrom(f0_c,calib.EV_ROOT_CALIBTRK1,calib.EV_ROOT_CALIBTRK2);
    TrkParams::FillMask(f0_c,calib.EV_ROOT_CALIBTRK1,calib.EV_ROOT_CALIBTRK2);
//    mask_.Set(f0_c,calib.EV_ROOT_CALIBTRK1,calib.EV_ROOT_CALIBTRK2);

//    };
    // =============================================================
    // retrieve calibration file needed to uncompress data
    // =============================================================
    // if the run was compressed using default calib
    // load truncated pedestals from default
    // otherwise reload them from on-line calibration
    cout << "Truncated pedestals for uncompression:";
    if( calib104 ){

 	cout << " >> Loading default calibration: "<< calibpatht << endl; 
// 	ca->FillTCalibFrom(calibpatht);
	TrkParams::FillTCalibFrom(calibpatht);
	    
    }else{

	if ( !f0_c ) return false;
	cout << ">> Loading on-line calibration " << endl;
//	ca->FillTCalibFrom(f0_c,calib.EV_ROOT_CALIBTRK1,calib.EV_ROOT_CALIBTRK2);
	TrkParams::FillTCalibFrom(f0_c,calib.EV_ROOT_CALIBTRK1,calib.EV_ROOT_CALIBTRK2);
    };
    f0_c->Close();
//    delete f0_c;
    cout << "--------------------------------------"<<endl;
    
    calibload = kTRUE;
    return true;
};

Int_t TrkParams::GetIndex( UInt_t type ){
    Int_t index = -1;
    for (Int_t i=0; i< NTYPES; i++)if( type == trkparamtype[i] )index=i;
    return index;
};

/**
 * Set tracker parameters proper for a run querying the DB.
 * @param glrun Input run info 
 * @param dbc   Database 
 * @param type  Parameter type
 * @return FALSE if parameter set fails.
 * If glrun==NULL the last parameter entry from GL_PARAM table is set. 
 * If dbc==NULL, try to open a connection taking DB-connection
 * parameters from PAMELA environment variables.
 */
Bool_t TrkParams::Set( GL_RUN* glrun , TSQLServer* dbc , UInt_t type){

//    if( !dbc )return false;


    if( !dbc || (dbc && !dbc->IsConnected()) ){

	cout <<" Missing DB connection -- check PAMELA environment variables "<<endl; 
	TString host = "mysql://localhost/pamelaprod";
	TString user = "anonymous";
	TString psw = "";
	const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
	const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
	const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
	if ( !pamdbhost ) pamdbhost = "";
	if ( !pamdbuser ) pamdbuser = "";
	if ( !pamdbpsw ) pamdbpsw = "";
	if ( strcmp(pamdbhost,"") ) host = pamdbhost;
	if ( strcmp(pamdbuser,"") ) user = pamdbuser;
	if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
	cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl; 
	cout<<"Connecting to DB"<<endl;
	cout<<"HOST "<<host<<endl;
	cout<<"USER "<<user<<endl;
	cout<<"PSW  "<<psw<<endl;
	dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
	if( !dbc )return false;
	if( !dbc->IsConnected() )return false;    
	cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl; 

    }

    UInt_t index = TrkParams::GetIndex( type );
    if(index < 0)return false;

    UInt_t runheadtime = numeric_limits<UInt_t>::max();
    if(glrun)runheadtime = glrun->RUNHEADER_TIME;
    else cout <<"Missing run info -- taking last entry "<<endl;

    GL_PARAM q4;
    if( q4.Query_GL_PARAM(runheadtime,type,dbc) )return false;

    gl[index]     = q4;
    glpath[index] = q4.PATH+q4.NAME;
    glload[index] = false;
    cout << "<< set TrkParams type "<<type<<" >>"<<endl;
    return true;
//    return TrkParams::Load(type);

};
/**
 * Set tracker parameters proper for the run, from the from input.
 * @param path  Path to parameters directory
 * @param type  Parameter type
 * @return FALSE if parameter set fails.
 */
Bool_t TrkParams::Set( TString path , UInt_t type){

    UInt_t index = TrkParams::GetIndex( type );
    if(index < 0)return false;

    gl[index].TYPE= type;
    glpath[index] = path;
    cout << "<< set TrkParams type "<<type<<" >>"<<endl;
    glload[index] = false;

//    return TrkParams::Load(type);
    return true;

};

/**
 * Set and load tracker parameters.
 * @param type  Parameter type
 * If no path as been yet specified, it set a default path.
 * @return FALSE if parameter set fails.
 */
Bool_t TrkParams::Set(UInt_t type){

    UInt_t index = TrkParams::GetIndex( type );
    if(index < 0)return false;

    gl[index].TYPE= type;
    glload[index] = false;

    // -----------------------------------------
    // check if enviroment variables are defined
    // -----------------------------------------
    glpath[index] = gSystem->Getenv("PAM_CALIB");

    if( !glpath[index].IsNull() ){
    
	switch(type){
	case 1:
	    glpath[index].Append("/trk-param/field_param-0/");
	    break;
	case 2:
	    glpath[index].Append("/trk-param/mip_param-0/");
	    break;
	case 3:
	    glpath[index].Append("/trk-param/charge_param-1/");
	    break;
	case 4:
	    glpath[index].Append("/trk-param/eta_param-0/");
	    break;
	case 5:
	    glpath[index].Append("/trk-param/align_param-0/");
	    break;
	case 6:
	    glpath[index].Append("/trk-param/mask_param-1/");
	    break;
	case 7:
	    glpath[index].Append("/trk-param/default_calib/");
	    break;
	}
	
    }else{
	cout << " TrkParams::Set( UInt_t type) ==> No PAMELA environment variables defined "<<endl;
	return false;
	
    }

    cout << "<< set TrkParams type "<<type<<" >>"<<endl;
//    return TrkParams::Load(type);
    return true;

};
/**
 * Set tracker parameters.
 * If no path as been yet specified, try to read from DB 
 * otherwise it set a default path.
 * @return FALSE if parameter set fails.
 */
Bool_t  TrkParams::Set( ){

    Bool_t connected = false;

    cout <<" Missing DB connection -- check PAMELA environment variables "<<endl; 
    TString host = "mysql://localhost/pamelaprod";
    TString user = "anonymous";
    TString psw = "";
    const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
    const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
    const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
    if ( !pamdbhost ) pamdbhost = "";
    if ( !pamdbuser ) pamdbuser = "";
    if ( !pamdbpsw ) pamdbpsw = "";
    if ( strcmp(pamdbhost,"") ) host = pamdbhost;
    if ( strcmp(pamdbuser,"") ) user = pamdbuser;
    if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
    cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl; 
    cout<<"Connecting to DB"<<endl;
    cout<<"HOST "<<host<<endl;
    cout<<"USER "<<user<<endl;
    cout<<"PSW  "<<psw<<endl;
    TSQLServer *dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
    if( dbc || dbc->IsConnected() )connected = true;
    else cout << " >> FAILED!!!"<<endl;
    cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl; 
    
    if(connected){
	if ( !TrkParams::Set(0,dbc) )return false;
	dbc->Close();
    }else for(Int_t i=0; i<NTYPES; i++)if ( !TrkParams::Set(trkparamtype[i])) return false;
}

/**
 * Load tracker parameter, if not already loaded.
 * @param type Parameter type
 * ( Any call to TrkParams::Set(...) forces reloading )
 */
Bool_t TrkParams::Load(UInt_t type){
    
    Int_t index = TrkParams::GetIndex( type );
    if( index<0 )return false;
    if( glload[index] )return true;

    if( glpath[index].IsNull() )TrkParams::Set(type);

    switch( type ){
    //------------------------------------------
    case 1:	
	cout <<type<<" -- Loading magnetic field "<<glpath[index]<<endl;
	strcpy(path_.path,glpath[index].Data());
	path_.pathlen = glpath[index].Length();
	path_.error   = 0;
	readb_();
	if(path_.error) return false;
	glload[index]=true;
	break;
    //------------------------------------------
    case 2:	
	cout <<type<<" -- Loading ADC-to-MIP conversion parameters "<<glpath[index]<<endl;
	strcpy(path_.path,glpath[index].Data());
	path_.pathlen = glpath[index].Length();
	path_.error   = 0;
	readmipparam_();
	if(path_.error) return false;
	glload[index]=true;
	break;
    //------------------------------------------
    case 3:	
	cout <<type<<" -- Loading charge-correlation parameters "<<glpath[index]<<endl;
	strcpy(path_.path,glpath[index].Data());
	path_.pathlen = glpath[index].Length();
	path_.error   = 0;
	readchargeparam_();
	if(path_.error) return false;
	glload[index]=true;
	break;
    //------------------------------------------
    case 4:	
	cout <<type<<" -- Loading p.f.a. parameters "<<glpath[index]<<endl;
	strcpy(path_.path,glpath[index].Data());
	path_.pathlen = glpath[index].Length();
	path_.error   = 0;
	readetaparam_();
	if(path_.error) return false;
	glload[index]=true;
	break;
    //------------------------------------------
    case 5:	
	cout <<type<<" -- Loading alignment parameters "<<glpath[index]<<endl;
	strcpy(path_.path,glpath[index].Data());
	path_.pathlen = glpath[index].Length();
	path_.error   = 0;
	readalignparam_();
	if(path_.error) return false;
	glload[index]=true;
	break;
    //------------------------------------------
    case 6:	
	cout <<type<<" -- Loading VK-mask parameters "<<glpath[index]<<endl;
	strcpy(path_.path,glpath[index].Data());
	path_.pathlen = glpath[index].Length();
	path_.error   = 0;
	readvkmask_();
	if(path_.error) return false;
	glload[index]=true;
	break;
    //------------------------------------------
    case 7:	
	cout <<type<<" -- Loading default calibration "<<glpath[index]<<endl;
	strcpy(path_.path,glpath[index].Data());
	path_.pathlen = glpath[index].Length();
	path_.error   = 0;
//	fillpedsigfromdefault_();
	cout << ">> not implemented <<"<<endl;
	if(path_.error) return false;
	glload[index]=true;
	break;
    default:
	cout << " Bool_t TrkParams::Load(UInt_t type) -- type "<< type << "not yet implemented"<<endl;
	return false;
    };

    return true;
}

Bool_t TrkParams::Load( ){
    for(Int_t i=0 ; i<NTYPES; i++)if(!glload[i] && trkparamtype[i]>0 )if( !TrkParams::Load(trkparamtype[i]) )return false;
    return true;
};
Bool_t TrkParams::IsLoaded( ){
    for(Int_t i=0 ; i<NTYPES; i++)if(!glload[i] && trkparamtype[i]>0 )return false;
    return true;
};

ClassImp(TrkParams);