PndEmcClusterCalibrationPar.h

Go to the documentation of this file.

// 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():fClusterCalibrationPars() { fClusterCalibrationPars.clear();};
                ~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;
                PndEmcClusterCalibrationPar(const  PndEmcClusterCalibrationPar& L);
                PndEmcClusterCalibrationPar& operator= (const  PndEmcClusterCalibrationPar&) {return *this;};
                ClassDef(PndEmcClusterCalibrationPar,1);
};

#endif