PndEmcClusterCalibrator.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// Description:
//      Class PndEmcClusterCalibrator
//      Do an energy and position corrections for EMC cluster
//      PndEmcClusterCalibrator is a factory which call specific method depending on the input parameter "method"
//      method=1 - PndEmcClusterHistCalibrator is called, which perform correction from the interpolation on 2-dimensional histogram on (Energy, theta)
//      method=2 - PndEmcClusterSimpleCalibrator is called, where energy correction is parametrised as a function of (Energy, theta)
//      PndEmcAbsClusterCalibrator - abstract interface class
//
// Author List:
//      D.Melnychuk
//      PndEmcClusterHistCalibrator class is based on code PndEmcMakeCorr.cxx
//      (A. Biegun, M. Babai)
//      PndEmcClusterSimpleCalibrator class is based on EmcPhotonSimpleCalib class in Babar framework (Jan Zhong)
//------------------------------------------------------------------------
//#pragma once
#ifndef PNDEMCCLUSTERCALIBRATOR_H
#define PNDEMCCLUSTERCALIBRATOR_H

#include "TObject.h"
#include "TVector3.h"
#include <TH2.h>
#include <TFile.h>
#include <TString.h>
#include <cstdlib>

class PndEmcCluster;
class PndEmcDigi;
class PndEmcClusterCalibrationParObject;

// Abstract interface class for EMC Cluster calibrator
class PndEmcAbsClusterCalibrator
{
  
        public:
  
  // Constructors
                PndEmcAbsClusterCalibrator(Int_t version=1);

  // Destructor
                virtual ~PndEmcAbsClusterCalibrator();

  // Methods 
                virtual Double_t Energy(PndEmcCluster *clust, Int_t pid=22) = 0;
                virtual TVector3 Where(PndEmcCluster *clust, Int_t pid=22) = 0;
  
        protected: 
        
                Int_t fVersion;
                TString fPath;
};

// Factory for creation EMC Cluster calibrator
class PndEmcClusterCalibrator
{
  
        public:
  
  // Constructors
                PndEmcClusterCalibrator();

  // Destructor
                virtual ~PndEmcClusterCalibrator();
  
                static PndEmcAbsClusterCalibrator *MakeEmcClusterCalibrator(Int_t method, Int_t version=1);

        private:

                PndEmcClusterCalibrator(const  PndEmcClusterCalibrator& L);
                PndEmcClusterCalibrator& operator= (const  PndEmcClusterCalibrator&) {return *this;};

};

class PndEmcClusterHistCalibrator: public PndEmcAbsClusterCalibrator
{
        public:
                PndEmcClusterHistCalibrator(Int_t version=1);
                virtual ~PndEmcClusterHistCalibrator();
                
  // Methods 
        virtual Double_t Energy(PndEmcCluster *clust, Int_t pid=22);
        virtual TVector3 Where(PndEmcCluster *clust, Int_t pid=22);
        
        private:
        Int_t FindTheBin(TH2* lookup_table, Float_t value_x, Float_t value_y, Int_t &bin_x, Int_t &bin_y);

        Double_t GetValueInZ(TH2 *lookup_table, Float_t value_x, Float_t value_y, Bool_t use_interpolation = kFALSE);
        
        TFile *fPhoton, *fElectron, *fPion;
        TH2F *fHEnergyRatioBarrelPhoton, *fHEnergyRatioFwdPhoton, *fHEnergyRatioBwdPhoton, *fHEnergyRatioShashlykPhoton;
        TH2F *fHThetaDiffBarrelPhoton, *fHThetaDiffFwdPhoton, *fHThetaDiffBwdPhoton, *fHThetaDiffShashlykPhoton;
//      TH2F *fHEnergyRatioBarrelElectron, *fHEnergyRatioFwdElectron, *fHEnergyRatioBwdElectron, *fHEnergyRatioShashlykElectron;
//      TH2F *fHThetaDiffBarrelElectron, *fHThetaDiffFwdElectron, *fHThetaDiffBwdElectron, *fHThetaDiffShashlykElectron;
//      TH2F *fHEnergyRatioBarrelPion, *fHEnergyRatioFwdPion, *fHEnergyRatioBwdPion, *fHEnergyRatioShashlykPion;
//      TH2F *fHThetaDiffBarrelPion, *fHThetaDiffFwdPion, *fHThetaDiffBwdPion, *fHThetaDiffShashlykPion;

        PndEmcClusterHistCalibrator(const  PndEmcClusterHistCalibrator& L);
        PndEmcClusterHistCalibrator& operator= (const  PndEmcClusterHistCalibrator&) {return *this;};
};

class PndEmcClusterSimpleCalibrator: public PndEmcAbsClusterCalibrator
{
        public:
                PndEmcClusterSimpleCalibrator(Int_t version=1);
                virtual ~PndEmcClusterSimpleCalibrator();
                
                void Init() {};
                
  // Methods 
                virtual Double_t Energy(PndEmcCluster *clust, Int_t pid=22);
                virtual TVector3 Where(PndEmcCluster *clust, Int_t pid=22);
        
        private: 
                PndEmcClusterCalibrationParObject *fParObject;

                PndEmcClusterSimpleCalibrator(const  PndEmcClusterSimpleCalibrator& L);
                PndEmcClusterSimpleCalibrator& operator= (const  PndEmcClusterSimpleCalibrator&) {return *this;};
};
#endif // PndEmcClusterCalibrator_HH