//----------------------------------------------------------------------
// File and Version Information:
//      $Id: //
// Description:
//      Class EmcWaveformToDigi. Module to take the ADC waveforms and produces digi.
// 
//	 Software developed for the BaBar Detector at the SLAC B-Factory.
// Adapted for the PANDA experiment at GSI		
//		
// Author List:
//      Phil Strother                  Original Author
// Dima Melnichuk - adaption for PANDA		
// Copyright Information:
//      Copyright (C) 1996             Imperial College
//
//----------------------------------------------------------------------

#include "TClonesArray.h"
//#include "TObjArray.h"		
		
#include "EmcWaveformToDigi.h"
#include "EmcWaveform.h"
#include "EmcDigi.h"
#include "EmcDigiPar.h"		

#include "CbmRootManager.h"
#include "CbmRunAna.h"
#include "CbmRuntimeDb.h"

//#include "TStopwatch.h"

#include <iostream>
#include <map>

using std::endl;
using std::fstream;

EmcWaveformToDigi::EmcWaveformToDigi()
{
	fDigiPosMethod="surface";// "surface" or "depth"
	fEmcDigiRescaleFactor=1.08;
	fEmcDigiPositionDepth=12.5;
}

//--------------
// Destructor --
//--------------

EmcWaveformToDigi::~EmcWaveformToDigi()
{
}


InitStatus EmcWaveformToDigi::Init()
{
	// Get RootManager
	CbmRootManager* ioman = CbmRootManager::Instance();
	if ( ! ioman )
	{
		cout << "-E- EmcWaveformToDigi::Init: "
		<< "RootManager not instantiated!" << endl;
		return kFATAL;
	}
	
	// Get input array
	fWaveformArray = (TClonesArray*) ioman->GetObject("EmcWaveform");
	if ( ! fWaveformArray ) {
		cout << "-W- EmcWaveformToDigi::Init: "
		<< "No EmcWaveform array!" << endl;
		return kERROR;
	}
	
	// Create and register output array
	fDigiArray = new TClonesArray("EmcDigi");

	ioman->Register("EmcDigi","Emc",fDigiArray,kTRUE);
	sampleRate=fDigiPar->GetSampleRate();
	threshold=fDigiPar->GetThreshold();
	
	if (!fDigiPosMethod.compare("surface"))
	{
		EmcDigi::selectDigiPositionMethod( EmcDigi::surface, 1., 0. );
	}
	else if (!fDigiPosMethod.compare("depth"))
	{
	   EmcDigi::selectDigiPositionMethod( EmcDigi::depth, 
			 fEmcDigiRescaleFactor, fEmcDigiPositionDepth);
	}
	else 
	{
		cout << "-W- EmcWaveformToDigi::Init: "
		<< "Unknown digi position method!" << endl;
		return kERROR;
	}

	cout << "-I- EmcWaveformToDigi: Intialization successfull" << endl;
	
	return kSUCCESS;
}

void EmcWaveformToDigi::Exec(Option_t* opt)
{
	// Reset output array
	if ( ! fDigiArray ) Fatal("Exec", "No Digi Array");
  	fDigiArray->Delete();
	double peakPosition;
	double max_energy;
	double digi_time;
	int i_digi=0; //index of digi in TClonesArray
	Int_t nWaveforms = fWaveformArray->GetEntriesFast();
	//cout<<"EmcWaveformToDigi: "<<nWaveforms<<" waveforms to convert"<<endl;
	for (Int_t iWaveform=0; iWaveform<nWaveforms; iWaveform++) {
		EmcWaveform* theWaveform = (EmcWaveform*) fWaveformArray->At(iWaveform);
		hitList=theWaveform->itsHitList();
		int detId=theWaveform->GetDetectorId();
		int trackId=theWaveform->GetTrackId();
		theWaveform->fitPeak(max_energy,peakPosition);
		max_energy/=theWaveform->get_scale();
		digi_time=peakPosition/sampleRate;
		if (max_energy>threshold)
		{
			new((*fDigiArray)[i_digi]) EmcDigi(trackId,detId, max_energy, peakPosition,hitList);
			i_digi++;
			
		}
		hitList.clear();
	}

}

void EmcWaveformToDigi::SetParContainers() {

  // Get run and runtime database
  CbmRunAna* run = CbmRunAna::Instance();
  if ( ! run ) Fatal("SetParContainers", "No analysis run");

  CbmRuntimeDb* db = run->GetRuntimeDb();
  if ( ! db ) Fatal("SetParContainers", "No runtime database");

  // Get Emc digitisation parameter container
  fDigiPar = (EmcDigiPar*) db->getContainer("EmcDigiPar");
 
}

ClassImp(EmcWaveformToDigi)