//*-- Author  : Rafał Lalik
//*-- Created : 01.06.2016

//_HADES_CLASS_DESCRIPTION
/////////////////////////////////////////////////////////////
//  HFRpcCalSim
//
//  This class contains Forward Rpc detector Cal data
//
//  Containing calibrated Time and energy loss
//
/////////////////////////////////////////////////////////////

#include "hfrpccal.h"

#include <hgeomvector.h>

HFRpcCal::HFRpcCal()
{
    fSector = -1;
    fTimeN = fTimeF = -10.F;
    fChargeN = fChargeF = -10.F;
    fTof = -1.0;
    fV = -1000.0;
    fX = fY = -1000.0;
    fZ = -1000.0;
}

void HFRpcCal::reconstructHits(Float_t tl, Float_t tr, Float_t ql, Float_t qr, Float_t z, Float_t v, Float_t L,
                               Float_t po, Float_t to, HGeomVector *pos, Float_t cosa, Float_t sina, Float_t twc0,
                               Float_t twc1, Float_t twc2)
{
    fTimeN = tl;
    fTimeF = tr;
    fChargeN = ql;
    fChargeF = qr;

    Float_t twc = twc0 * ((Float_t) TMath::Exp(-twc1 * ((fChargeN + fChargeF) - twc2))); // time walk correction
    fTof = (fTimeN + fTimeF - L / v) / 2.0 - to - twc;                    // tof
    fV = (fTimeF - fTimeN) * v / 2. - po;                                 // pos along time axis
    fX = pos->X() - fV * sina;
    fY = pos->Y() + fV * cosa;
    fZ = z;
}

void HFRpcCal::print() const
{
    printf(" RPC cal hit: sec=%d strip=%2d:  t_n=%7.3f t_f=%7.3f q_n=%6.2f q_f=%6.2f "
           " tof=%7.3f  v=%8.2f  x,y=%7.2f,%7.2f  z=%6.2f\n",
           fSector, fStrip, fTimeN, fTimeF, fChargeN, fChargeF, fTof, fV, fX, fY, fZ);
}