//---------------------------------------------
// Go4 Tasca analysis
// Author: Hans G. Essel
//         H.Essel@gsi.de
// GSI, Experiment Electronics, Data Processing
// modified by r.maendl
// July 2021 - cleanup by JAM (j.adamczewski@gsi.de)
//---------------------------------------------

#include "TascaCalibration.h"

#include "TMath.h"
#include <iostream>
#include <fstream>
using namespace std;

#include "TGo4Log.h"

//***********************************************************
TascaCalibration::TascaCalibration() : TGo4Parameter()
{}
//***********************************************************
TascaCalibration::TascaCalibration(const char* name) : TGo4Parameter(name)
{
    cout << "Tasca> TascaCalibration: " << name << " created" << endl;
    UseCalibration=kFALSE;
}
//***********************************************************
TascaCalibration::~TascaCalibration(){
    cout << "Tasca> TascaCalibration: " << GetName() << " deleted" << endl;
}
//-----------------------------------------------------------
Int_t TascaCalibration::PrintParameter(){
  return 0;
  cout << "Calibration " << GetName()<<":" <<endl;
  return 0;
}
//-----------------------------------------------------------
Bool_t TascaCalibration::Preset(){
	Int_t i;


	for(i=0;i<TASCA_CALI_NUMSTOP;i++){
	  fdStopXL_a0[i]=0.;fdStopYL_a0[i]=0.;
	  fdStopXL_a1[i]=1.;fdStopYL_a1[i]=1.;
	  fdStopXL_a2[i]=0.;fdStopYL_a2[i]=0.;
	}
	
	for(i=0;i<TASCA_CALI_NUMGAMMA;i++){
	    fdGammaE_a0[i]=0.;  fdGammaT_a0[i]=0.;
		fdGammaE_a1[i]=1.;  fdGammaT_a1[i]=1.;
		fdGammaE_a2[i]=0.;  fdGammaT_a2[i]=0.;


	}

	return kTRUE;
  }

//-----------------------------------------------------------
void TascaCalibration::ReadSingleCoefficients(const char * file
		, UInt_t size, Double_t *a0, Double_t *a1, Double_t *a2){
	Text_t line[128];
	Double_t v0,v1,v2;
	UInt_t ii;
	//File=file;
	char name[32];
	std::ifstream database(file);
    if(!database.is_open()){
	    cout << "Tasca> TascaCalibration: Error open " << file << endl;
	    return;
	}
	cout << "Tasca> TascaCalibration: Coefficients from " << file<< endl;
	while(1){
		database.getline(line,127,'\n' ); // read whole line
		if(database.eof() || !database.good()) break;
		if(strstr(line,"#") || strstr(line,"!")){
			cout << line << endl;
			continue;// skip any comments
		}
		sscanf(line,"%s %d %lf %le %le",name,&ii,&v0,&v1,&v2);
//		cout << name<<" "<<ii<<" "
//			<< v0<<" "<<v1<<" "<<v2<<endl;
		if(ii>=size) cout << "Index error: " << ii << " should less equal "<<size<< endl;
		a0[ii]=v0;
		a1[ii]=v1;
		a2[ii]=v2;
	}
	database.close();
	return;

}
void TascaCalibration::ReadCoefficients(){
	TString full;
	TString pref=prefix;

	UInt_t size=TASCA_CALI_NUMSTOP;
	full.Form("%s_StopXL.txt",prefix.Data());
	ReadSingleCoefficients(full.Data(),size,fdStopXL_a0,fdStopXL_a1,fdStopXL_a2);
	full.Form("%s_StopYL.txt",prefix.Data());
	ReadSingleCoefficients(full.Data(),size,fdStopYL_a0,fdStopYL_a1,fdStopYL_a2);
	size=TASCA_CALI_NUMGAMMA;
	full.Form("%s_GammaE.txt",prefix.Data());
	ReadSingleCoefficients(full.Data(),size,fdGammaE_a0,fdGammaE_a1,fdGammaE_a2);
	full.Form("%s_GammaT.txt",prefix.Data());
	ReadSingleCoefficients(full.Data(),size,fdGammaT_a0,fdGammaT_a1,fdGammaT_a2);
	//	for(UInt_t i=0;i<size;i++)
//	cout <<fdStopXL_a0[i]<<" "<< fdStopXL_a1[i]<<" "<<  fdStopXL_a2[i]<<endl;
}