PndEmcErrorMatrix.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// Description:
//      Class PndEmcErrorMatrix
//      Calculate Error Matrix for the given EmcCluster
//      with parametrization defined by the given parameter PndEmcErrorMatrixPar
//
//------------------------------------------------------------------------

//#pragma once
#ifndef PNDEMCERRORMATRIX_H
#define PNDEMCERRORMATRIX_H

#include "TMatrixD.h"
#include "PndEmcCluster.h"
#include "TString.h"
#include <cstdlib>

class PndEmcErrorMatrixPar;
class PndEmcErrorMatrixParObject;

class PndEmcErrorMatrix 
{
        public:
 
        // Constructors
        PndEmcErrorMatrix();
        // Destructor
        ~PndEmcErrorMatrix();
        
        void Init(PndEmcErrorMatrixParObject *par);
        void InitFromFile(Int_t geomVersion);
        PndEmcErrorMatrixParObject* GetParObject();
        
        // Error Matrix for Cluster
        // 4x4 matrix in E, theta, phi, R
        TMatrixD GetErrorMatrix(const PndEmcCluster &cluster) const;
                
        // 4x4 matrix in (px, py, pz, E)
        TMatrixD Get4MomentumErrorMatrix(const PndEmcCluster &cluster) const;

        // 7x7 matrix in (x, y, z, px, py, pz, E)
        TMatrixD GetErrorP7(const PndEmcCluster &cluster) const;
         
        private:
                PndEmcErrorMatrixParObject *fErrorMatrixParObject;

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

// helper function to convert error matrix
// copied from BbrGeom/BbrError.cc
TMatrixD similarityWith(const TMatrixD& mat, const TMatrixD& m1);

#endif