/////////////////////////////////////////////////////////////
// PndEmcClusterCalibrationPar
//
// Container class for parametrization of EMC cluster energy correction 
// class is inherited from FairParGenericSet
//
// Parametrization of EMC energy correction is a function of energy and polar angle and has different parameters in each of 4 regions:
// barrel, forward and backward endcap, shashlyk. Two different parameters set are used in barrel for energy below and above 1 GeV
// For the barrel , forward and backward endcap the following function is used with 10 parameters
// eout1=e* exp(factor1);
// factor1=
// p0
// +p1*log(e1)
// +p2*log(e1)*log(e1)
// +p3*log(e1)*log(e1)*log(e1)
// +p4*cos(theta1)
// +p5*cos(theta1)*cos(theta1)
// +p6*cos(theta1)*cos(theta1)*cos(theta1)
// +p7*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
// +p8*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
// +p9*log(e1)*cos(theta1);
// For shashlyk the following parametrization is used with 5 parameters:
// eout=(3.31694-0.0183379/sqrt(e1)+0.0327113/e1+0.00040156/(e1*e1)
// -0.00641305/(e1*sqrt(e1)))*e1;
// 
// The index for parameter set within
// map<Int_t,std::vector<Double_t> > fClusterCalibrationPars
// are: 1 - barrel (below 1 GeV), 2 - barrel (above 1 GeV)
// 3 - forward endcap, 4 -backward endcap, 5 - shashlyk 
/////////////////////////////////////////////////////////////

#pragma once
#ifndef PNDEMCCLUSTERCALIBRATIONPAR_H
#define PNDEMCCLUSTERCALIBRATIONPAR_H

#include "FairParGenericSet.h"
#include "FairParamList.h"
#include <cstdlib>
#include <map>
#include <vector>

class PndEmcClusterCalibrationParObject: public TObject
{
	public:
		void SetCalibrationPar(Int_t iParSet, Double_t *pars);
		void GetCalibrationPar(Int_t iParSet, Double_t *pars);
		
		PndEmcClusterCalibrationParObject(){};
		~PndEmcClusterCalibrationParObject(){};
	private:
		std::map<Int_t,std::vector<Double_t> > fClusterCalibrationPars;
		ClassDef(PndEmcClusterCalibrationParObject,1);
};


class PndEmcClusterCalibrationPar : public FairParGenericSet
{
	public:
		void GetClusterCalibrationParameters(Int_t iParSet, Double_t *pars);
		void SetClusterCalibrationObject(PndEmcClusterCalibrationParObject *ParObject){fParObject = ParObject;};
		PndEmcClusterCalibrationParObject *GetParObject() {return fParObject;};
		PndEmcClusterCalibrationPar (const char* name="PndEmcClusterCalibrationPar",
			    const char* title="Emc cluster calibration parameter",
			    const char* context="TestDefaultContext");
		~PndEmcClusterCalibrationPar(void){};
		
		Bool_t IsValid();
		
		void clear(void){};
		
		void putParams(FairParamList* list);
		Bool_t getParams(FairParamList* list);
	
	private :
		PndEmcClusterCalibrationParObject *fParObject;
		ClassDef(PndEmcClusterCalibrationPar,1);
};

#endif