PndEmcClusterSimpleCalibrator Class Reference

#include <PndEmcClusterCalibrator.h>

Inheritance diagram for PndEmcClusterSimpleCalibrator:

Public Member Functions
	PndEmcClusterSimpleCalibrator (Int_t version=1)
virtual	~PndEmcClusterSimpleCalibrator ()
void	Init ()
virtual Double_t	Energy (PndEmcCluster *clust, Int_t pid=22)
virtual TVector3	Where (PndEmcCluster *clust, Int_t pid=22)

Protected Attributes
Int_t	fVersion
TString	fPath

Private Member Functions
	PndEmcClusterSimpleCalibrator (const PndEmcClusterSimpleCalibrator &L)
PndEmcClusterSimpleCalibrator &	operator= (const PndEmcClusterSimpleCalibrator &)

Private Attributes
PndEmcClusterCalibrationParObject *	fParObject

Detailed Description

Definition at line 101 of file PndEmcClusterCalibrator.h.

Constructor & Destructor Documentation

PndEmcClusterSimpleCalibrator::PndEmcClusterSimpleCalibrator

(

Int_t

version = 1

)

Definition at line 465 of file PndEmcClusterCalibrator.cxx.

References fParObject, PndEmcAbsClusterCalibrator::fPath, and TString.

                                                                         :
fParObject(new PndEmcClusterCalibrationParObject())
{
        fPath = getenv("VMCWORKDIR");
        fPath += "/macro/params/";
 
        TString fileNamePhoton;
        fileNamePhoton.Form("emc_correction_par_gamma_%i.root",version);
        fileNamePhoton=fPath+fileNamePhoton;
        TFile *fPhoton;
        
//      TString fileNameElectron;
//      TString fileNamePion;
//      fileNameElectron.Form("emc_correction_par_electron_%i.root",version);
//      fileNamePion.Form("emc_correction_par_pion_%i.root",version);
//      fileNameElectron=fPath+fileNameElectron;
//      fileNamePion=fPath+fileNamePion;
//      TFile *fElectron, *fPion;

        fPhoton = new TFile(fileNamePhoton,"READ");
        if (fPhoton->IsZombie())
        { 
                std::cout<<"EMC cluster correction file "<<fileNamePhoton<<" for photons does not exist for given version"<<std::endl;
        } else 
        {
                // Photon
                fPhoton->GetObject("PndEmcClusterCalibrationParObject",fParObject);
                if(fParObject == NULL){
                        std::cout << "-E- PndEmcClusterSimpleCalibrator: Could not get Emc cluster calibration information from file " << fPhoton << std::endl;
                } else 
                {
                        std::cout<<"EMC photon calibration parameters are read succesfully"<<std::endl;
                }
        }
//      fElectron = new TFile(fileNameElectron,"READ");
//      if (fElectron->IsZombie())
//      {
//              std::cout<<"EMC cluster correction file "<<fileNameElectron<<" for electrons does not exist for given version"<<std::endl;
//      } else
//      {
//              // Photon
//              fElectron->GetObject("PndEmcClusterCalibrationParObject",fParObject);
//              if(fParObject == NULL){
//                      std::cout << "-E- PndEmcClusterSimpleCalibrator: Could not get Emc cluster calibration information from file " << fElectron<< std::endl;
//              }
//      }
//      
//      fPion = new TFile(fileNamePion,"READ");
//      if (fPion->IsZombie())
//      {
//              std::cout<<"EMC cluster correction file "<<fileNamePion<<" for pions does not exist for given version"<<std::endl;
//      } else
//      {
//              // Photon
//              fPion->GetObject("PndEmcClusterCalibrationParObject",fParObject);
//              if(fParObject == NULL){
//                      std::cout << "-E- PndEmcClusterSimpleCalibrator: Could not get Emc cluster calibration information from file " << fPion << std::endl;
//              }
//      }
}

PndEmcClusterSimpleCalibrator::~PndEmcClusterSimpleCalibrator ( )

virtual

Definition at line 628 of file PndEmcClusterCalibrator.cxx.

{
}

PndEmcClusterSimpleCalibrator::PndEmcClusterSimpleCalibrator ( const PndEmcClusterSimpleCalibrator &  L )

private

Member Function Documentation

Double_t PndEmcClusterSimpleCalibrator::Energy ( PndEmcCluster *  clust, Int_t  pid = 22  )

virtual

Implements PndEmcAbsClusterCalibrator.

Definition at line 531 of file PndEmcClusterCalibrator.cxx.

References cos(), Double_t, PndEmcCluster::energy(), exp(), fParObject, PndEmcClusterCalibrationParObject::GetCalibrationPar(), log(), p1, p2, Pi, sqrt(), and PndEmcCluster::where().

{
        if (pid!=22)
        {
                std::cout<<"PndEmcClusterSimpleCalibrator:: EMC Cluster energy correction is not defined for pid="<<pid<<". Photon correction is used"<<std::endl;
        }
        
        Double_t e=clust->energy();
        TVector3 clusterPosition= clust->where();
        Double_t theta_cluster=clusterPosition.Theta();
        
        Double_t e1=e;
        Double_t theta1=theta_cluster;
        
        if ( (clusterPosition.Z() < 180.0)&&(theta_cluster<141.*TMath::Pi()/180.))
        {
                if (e<0.03) e1 = 0.03;
                if (e>8.0)  e1 = 8.0 ;
        }
        
        if ( (clusterPosition.Z() < 180.0)&&(theta_cluster>141.*TMath::Pi()/180.))
        {
                if (e<0.03) e1 = 0.03;
                if (e>2.0)  e1 = 2.0 ;
        }
        
        if (clusterPosition.Z() > 180.0)
        {
                if (e<0.01) e1 = 0.01;
                if (e>16.0)  e1 = 16.0 ;
        }

        Int_t iParSet; // EMC component and energy range (1 - barrel (below 1 GeV, 2 -barrel (above 1 geV, 3 -forward endcap, 4 -backward endcap, 5 -shashlyk) )
        
        if ( clusterPosition.Z() > 500.0)
        {
                iParSet=5;
        }
        else if ( (clusterPosition.Z() < 180.0)&&(theta_cluster>141.*TMath::Pi()/180.))
        {
                iParSet=4;
        }
        else  if ( (clusterPosition.Z() < 180.0)&&(theta_cluster<141.*TMath::Pi()/180.))
        {
                if (e1<1.0)
                        iParSet=1;
                else
                        iParSet=2;
        }
        else
        { 
                iParSet=3;
        }
        
        Double_t pars[10];
        fParObject->GetCalibrationPar(iParSet, pars);
        
        double p0, p1, p2, p3, p4, p5, p6, p7, p8, p9;
        double eout;
        
        if (iParSet==5)
        {
                p0=pars[0];
                p1=pars[1];
                p2=pars[2];
                p3=pars[3];
                p4=pars[4];
                eout=(p0-p1/sqrt(e1)+p2/e1+p3/(e1*e1)-p4/(e1*sqrt(e1)))*e;
        } else
        { 
                p0=pars[0];
                p1=pars[1];
                p2=pars[2];
                p3=pars[3];
                p4=pars[4];
                p5=pars[5];
                p6=pars[6];
                p7=pars[7];
                p8=pars[8];
                p9=pars[9];
                
                double 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);
                eout=e* exp(factor1);
        }
        
        return eout;
}

void PndEmcClusterSimpleCalibrator::Init ( )

inline

Definition at line 107 of file PndEmcClusterCalibrator.h.

{};

PndEmcClusterSimpleCalibrator& PndEmcClusterSimpleCalibrator::operator= ( const PndEmcClusterSimpleCalibrator &  )

inlineprivate

Definition at line 117 of file PndEmcClusterCalibrator.h.

{return *this;};

TVector3 PndEmcClusterSimpleCalibrator::Where ( PndEmcCluster *  clust, Int_t  pid = 22  )

virtual

Implements PndEmcAbsClusterCalibrator.

Definition at line 526 of file PndEmcClusterCalibrator.cxx.

References PndEmcCluster::where().

{
        return clust->where();
}

Member Data Documentation

PndEmcClusterCalibrationParObject* PndEmcClusterSimpleCalibrator::fParObject

private

Definition at line 114 of file PndEmcClusterCalibrator.h.

Referenced by Energy(), and PndEmcClusterSimpleCalibrator().

TString PndEmcAbsClusterCalibrator::fPath

protectedinherited

Definition at line 50 of file PndEmcClusterCalibrator.h.

Referenced by PndEmcClusterHistCalibrator::PndEmcClusterHistCalibrator(), and PndEmcClusterSimpleCalibrator().

Int_t PndEmcAbsClusterCalibrator::fVersion

protectedinherited

Definition at line 49 of file PndEmcClusterCalibrator.h.

---

The documentation for this class was generated from the following files:

• PndEmcClusterCalibrator.h
• PndEmcClusterCalibrator.cxx