//interpolation usage for to find the right residual energy and right A, using
//the final interval returned by bisection
void CsICalib::Interpolate(){
//L->Log<<"esi1 : "<<esi1<<" ecsi1 : "<<ecsi1<<endl;
//L->Log<<"esi2 : "<<esi2<<" ecsi2 : "<<ecsi2<<endl;
sRefECsI= GetInterpolationD(esi1,ecsi1,esi2,ecsi2,eEnergySi);
//L->Log<<"sRefECsI : "<<sRefECsI<<endl;
iA=GetInterpolationD(esi1,(Double_t)left,esi2,(Double_t)right,eEnergySi);
//L->Log<<"iA : "<<iA<<endl;
}
//interpolation usage for to find the right residual energy and right A, using
//the final interval returned by bisection
void CsICalib::Interpolate()
{
sRefECsI = GetInterpolationD(esi1, ecsi1, esi2, ecsi2, eEnergySi);
if ((sRefECsI < 0.) || (sRefECsI >= 3000.)) {
std::cerr << "CsICalib::Interpolate Error Report:" << std::endl;
std::cerr << "sRefECsI: " << sRefECsI << std::endl;
std::cerr << "esi1: " << esi1 << std::endl;
std::cerr << "ecsi1: " << ecsi1 << std::endl;
std::cerr << "esi2: " << esi2 << std::endl;
std::cerr << "ecsi2: " << ecsi2 << std::endl;
std::cerr << "eEnergySi: " << eEnergySi << std::endl;
}
assert((sRefECsI >= 0.) && (sRefECsI < 3000.));
iA = GetInterpolationD(esi1, (Double_t)left, esi2, (Double_t)right, eEnergySi);
}
