PndEmcErrorMatrixPar.cxx

Go to the documentation of this file.

// PndEmcErrorMatrixPar
//
// Container class for EMC error matrix parameter
// class is inherited from FairParGenericSet
//
// Parametrization of EMC error matrix
// functions used for parameterization
// Energy: Delta(E)/E = (a^2/E^power) + const^2 + (quadr/E)^2
// position: Delta(x)=(a*a/E^power) + const^2
// There are 10 parameters (engParA, engPower, engConst, engQuadr, pos1ParA, pos1Power,
// pos1Const, pos2ParA, pos2Power, pos2Const) for each component of EMC (barrel,
// fwd endcap, bwd endcap, shashlyk)

#include "PndEmcErrorMatrixPar.h"
#include <iostream>
#include "assert.h"

ClassImp(PndEmcErrorMatrixParObject);

void PndEmcErrorMatrixParObject::SetErrorMatrix(Int_t detectorComponent, Double_t *pars){
        std::vector<Double_t> parvec;
        parvec.assign(pars,pars+10);
        fErrorMatrixPars[detectorComponent]=parvec;
}

void PndEmcErrorMatrixParObject::GetErrorMatrix(Int_t detectorComponent, Double_t *pars){
        pars[0]=fErrorMatrixPars[detectorComponent][0];
        pars[1]=fErrorMatrixPars[detectorComponent][1];
        pars[2]=fErrorMatrixPars[detectorComponent][2];
        pars[3]=fErrorMatrixPars[detectorComponent][3];
        pars[4]=fErrorMatrixPars[detectorComponent][4];
        pars[5]=fErrorMatrixPars[detectorComponent][5];
        pars[6]=fErrorMatrixPars[detectorComponent][6];
        pars[7]=fErrorMatrixPars[detectorComponent][7];
        pars[8]=fErrorMatrixPars[detectorComponent][8];
        pars[9]=fErrorMatrixPars[detectorComponent][9];
}


ClassImp(PndEmcErrorMatrixPar)

PndEmcErrorMatrixPar::PndEmcErrorMatrixPar(const char* name,const char* title,
                             const char* context)
    : FairParGenericSet(name,title,context), fParObject(new PndEmcErrorMatrixParObject()) 
{
        fParObject = NULL;
}
Bool_t PndEmcErrorMatrixPar::IsValid()
{
        if (fParObject == NULL) return false;
        
        return true;
}

void PndEmcErrorMatrixPar::GetErrorMatrixParameters(Int_t detectorComponent, Double_t *pars){

        if(fParObject !=0){
                fParObject->GetErrorMatrix(detectorComponent,pars);
        } else {
                std::cout<<"Wrong detector component in PndEmcErrorMatrixPar"<<std::endl;
                abort();
        }

}

void PndEmcErrorMatrixPar::putParams(FairParamList* list)
{
        if(!list) return;
        list->addObject("PndEmcErrorMatrix",fParObject);
}

Bool_t PndEmcErrorMatrixPar::getParams(FairParamList* list)
{
        if (!list) return kFALSE;
        if(!list->fillObject("PndEmcErrorMatrix",fParObject)) return kFALSE;
        return kTRUE;
}