// ---------------------------------------------------------------------------------
// -----                           PndSttUnpack                                -----
// -----                        Author: S.Costanza                             -----
// ---------------------------------------------------------------------------------

//ROOT headers
#include "TClonesArray.h"
#include "TH1F.h"
#include "TMath.h"
//Fair heaeders
#include "FairRootManager.h"
#include "FairRunOnline.h"
#include "FairLogger.h"
#include "FairSttSource.h"

//Stt headers
#include "PndSttRawHit.h"
#include "PndSttUnpack.h"

#include <iostream>

using std::cout;
using std::endl;

//PndSttUnpack: Constructor
PndSttUnpack::PndSttUnpack(char *strCalDir,
			   Short_t type, Short_t subType,
			   Short_t procId,
			   Short_t subCrate, Short_t control)
  : FairUnpack(type, subType, procId, subCrate, control),
    fRawData(new TClonesArray("PndSttRawHit"))
{
}
	
PndSttUnpack::PndSttUnpack()
  : fRawData(new TClonesArray("PndSttRawHit"))
{
}
	
//Virtual PndSttUnpack: Public method
PndSttUnpack::~PndSttUnpack()
{
  delete fRawData;
}



//Init: Public method
Bool_t PndSttUnpack::Init()
{
  Register();
  return kTRUE;
}
	
	
	
//Register: Protected method
void PndSttUnpack::Register()
{
  FairRootManager *fMan = FairRootManager::Instance();
  if(! fMan) {
    return;
  }
  fMan->Register("SttRawHit", "Raw data from Stt", fRawData, kTRUE);
}
	
	
	
//DoUnpack: Public method
Bool_t PndSttUnpack::DoUnpack(Int_t *data, Int_t size)
{
  return kTRUE;
}


Bool_t PndSttUnpack::DoUnpack(Int_t nhits, TMatrixT<Float_t> HitInfo)
{
  
  Int_t evNum;        // global event number
  Int_t layer;        // layer number
  Int_t tube;         // tube number
  Float_t raw_led;    // raw leading edge discriminator [ns]
  Float_t raw_zct;    // raw zero crossing time [ns]
  Float_t raw_cfd;    // raw constant fraction discriminator [ns]
  Float_t raw_ampl;   // raw pulse amplitude (first maximum) [bins - a.u.]
  Float_t raw_de;     // raw energy loss (pulse area) [bins - a.u.]
  Float_t fpga_led;    // fpga leading edge discriminator [ns]
  Float_t fpga_zct;    // fpga zero crossing time [ns]
  Float_t fpga_cfd;    // fpga constant fraction discriminator [ns]
  Float_t fpga_ampl;   // fpga pulse amplitude (first maximum) [bins - a.u.]
  Float_t fpga_de;     // fpga energy loss (pulse area) [bins - a.u.]
  
  for (Int_t i=0; i<nhits; i++) {
    evNum     = (Int_t)HitInfo[i][0];
    layer     = (Int_t)HitInfo[i][1];
    tube      = (Int_t)HitInfo[i][2];
    tube     += (layer-1)*16;
    raw_led   = HitInfo[i][3];
    raw_zct   = HitInfo[i][4];
    raw_cfd   = HitInfo[i][5];
    raw_ampl  = HitInfo[i][6];
    raw_de    = HitInfo[i][7];
    fpga_led  = HitInfo[i][8];
    fpga_zct  = HitInfo[i][9];
    fpga_cfd  = HitInfo[i][10];
    fpga_ampl = HitInfo[i][11];
    fpga_de   = HitInfo[i][12];

    new ((*fRawData)[i]) PndSttRawHit(evNum, layer, tube, raw_led, raw_zct, raw_cfd, raw_ampl, raw_de, fpga_led, fpga_zct, fpga_cfd, fpga_ampl, fpga_de);
  }
  
  return kTRUE;
}
	 
//Reset: Public method
void PndSttUnpack::Reset()
{
  fRawData->Clear();
}
	
	
ClassImp(PndSttUnpack)