// ---------------------------------------------------------------------
// -----               KrattaCalibrator source file                -----
// -----            Created 10/11/13  by P. Pawlowski              -----
// ---------------------------------------------------------------------



#include <iostream>
#include <fstream>

#include "KrattaCalibrator.hxx"

using namespace Kratta;


void KrattaCalibrator::LoadRun(int irun){
  partype dummy;
  int runmin, runmax;
  for(int mod=0; mod<NMODULE; mod++){/// Loop over modules
    //partype * p = GetModuleData(mod);
    for(int i=0; i<NPAR; i++) par[mod][i] = 0;
    
    TString fileName = Form(fFILE_MASK,mod);
    ifstream f( fileName.Data() );
    while(f){/// Loop over run ranges
      f >> runmin >> runmax;
      if (!f) break;
      if(irun < runmin || irun > runmax){
	for(int i=0; i<NPAR; i++)
	  f >> dummy >> dummy >> dummy;
      }
      else{
	for(int i=0; i<NPAR; i++){
	  f >> dummy >> par[mod][i] >> dummy;
	}
	break;
      }
  }/// End loop over run ranges
  f.close();

  }/// End loop over modules
}



void KrattaCalibrator::Init(const char * mask){
  
      fFILE_MASK = TString(mask);
      Run        = 0;
      LoadRun(0);
      
}


void KrattaCalibrator::ls()const{
  for(int mod=0; mod<Kratta::NMODULE; mod++){
    std::cout << mod;
     const partype * p = GetModuleData(mod);
    for(int i=0; i<NPAR; i++) std::cout << "\t"<< p[i];
    std::cout << std::endl;
  }
}

float GetEnergy(const float * par, int Z, int A, ASYFadcPeak * peak){
  return 0;
}

void KrattaCalibrator::Calibrate(TRootKRATParticle * part, ASYFadcPeak * peak)const{
  const partype * p = GetModuleData(part->Module);
  part->Energy = GetEnergy(p, part->Z, part->A, peak);
}