MapEta Class Reference

#include <MapEta.h>

Collaboration diagram for MapEta:

Public Member Functions
	~MapEta ()
	MapEta (int isampling)
void	SetDirectory (const std::string &dir)
void	Initialize (int isampling)
int	Nx () const
int	Ny () const
float	Xmin () const
float	Xmax () const
float	Ymin () const
float	Ymax () const
void	GetData (double x, double y, double *resp, double *xt0, double *xt1, double *xt2) const
void	GetData0 (double x, double y, double *resp) const

Private Attributes
std::string	m_directory
int	m_init
int	m_nx
int	m_ny
float	m_xmin = 0.0F
float	m_xmax = 0.0F
float	m_ymin = 0.0F
float	m_ymax = 0.0F
float	m_deltax = 0.0F
float	m_deltay = 0.0F
float *	m_resp
float *	m_xt0
float *	m_xt1
float *	m_xt2

Detailed Description

Definition at line 12 of file MapEta.h.

Constructor & Destructor Documentation

~MapEta()

MapEta::~MapEta

(

)

Definition at line 29 of file MapEta.cxx.

{
  if (m_resp) delete [] m_resp;
  if (m_xt0) delete [] m_xt0;
  if (m_xt1) delete [] m_xt1;
  if (m_xt2) delete [] m_xt2;
}

MapEta()

MapEta::MapEta

(

int

isampling

)

Definition at line 16 of file MapEta.cxx.

  : m_directory ("/afs/cern.ch/atlas/offline/data/lar/calo_data"),
    m_init (0),
    m_nx (0),
    m_ny (0),
    m_resp (nullptr),
    m_xt0 (nullptr),
    m_xt1 (nullptr),
    m_xt2 (nullptr)
{
  Initialize(isampling);
}

Member Function Documentation

GetData()

void MapEta::GetData	(	double	x,
		double	y,
		double *	resp,
		double *	xt0,
		double *	xt1,
		double *	xt2 ) const

Definition at line 111 of file MapEta.cxx.

{
  *resp=1;
  *xt0=1;
  *xt1=0;
  *xt2=0;
  if (m_nx==0 || m_ny ==0) return;
  if (x<m_xmin ) x=m_xmin+0.01*m_deltax;
  if (x>=m_xmax) x=m_xmax-0.01*m_deltax;
  if (y<m_ymin ) y=m_ymin+0.01*m_deltay;
  if (y>=m_ymax) y=m_ymax-0.01*m_deltay;
 
  int ix,iy;
  ix = (int) ((x-m_xmin)/m_deltax);
  iy = (int) ((y-m_ymin)/m_deltay);
  float x0 = ((float) ix)*m_deltax+m_xmin;
  float x1 = x0+m_deltax;
  float y0 = ((float) iy)*m_deltay+m_ymin;
  float y1 = y0+m_deltay;
 
 
  if (ix<0 || ix+1 >= m_nx || iy<0 || iy+1 >= m_ny) {
    std::cout << "MapEta:  Out of range " << ix << " " << iy << std::endl;
    return;
  }
  double w[4];
  w[0]=(x1-x)*(y1-y);
  w[1]=(x-x0)*(y1-y);
  w[2]=(x1-x)*(y-y0);
  w[3]=(x-x0)*(y-y0);
 
  double sumw=0.;
  *resp=0;
  *xt0=0;
  for (int i=0;i<2;i++) {
    for (int j=0;j<2;j++) {
      int n=ix+i+(iy+j)*m_nx;
      int m=i+2*j;
      sumw +=w[m];
      *resp += m_resp[n]*w[m];
      *xt0 += m_xt0[n]*w[m];
      *xt1 += m_xt1[n]*w[m];
      *xt2 += m_xt2[n]*w[m];
    }
  }
  const double inv_sumw = 1. / sumw;
  if (sumw>0.) {
    *resp = *resp*inv_sumw;
    *xt0 = *xt0*inv_sumw;
    *xt1 = *xt1*inv_sumw;
    *xt2 = *xt2*inv_sumw;
  }
 
}

GetData0()

void MapEta::GetData0	(	double	x,
		double	y,
		double *	resp ) const

Definition at line 165 of file MapEta.cxx.

{
  *resp=1;
  if (m_nx==0 || m_ny ==0) return;
  if (x<m_xmin ) x=m_xmin+0.01*m_deltax;
  if (x>=m_xmax) x=m_xmax-0.01*m_deltax;
  if (y<m_ymin ) y=m_ymin+0.01*m_deltay;
  if (y>=m_ymax) y=m_ymax-0.01*m_deltay;
 
  int ix,iy;
  ix = (int) ((x-m_xmin)/m_deltax);
  iy = (int) ((y-m_ymin)/m_deltay);
  float x0 = ((float) ix)*m_deltax+m_xmin;
  float x1 = x0+m_deltax;
  float y0 = ((float) iy)*m_deltay+m_ymin;
  float y1 = y0+m_deltay;
 
  if (ix<0 || ix+1 >= m_nx || iy<0 || iy+1 >= m_ny) {
    std::cout << "MapEta:  Out of range " << ix << " " << iy << std::endl;
    return;
  }
  double w[4];
  w[0]=(x1-x)*(y1-y);
  w[1]=(x-x0)*(y1-y);
  w[2]=(x1-x)*(y-y0);
  w[3]=(x-x0)*(y-y0);
 
  double sumw=0.;
  *resp = 0;
  for (int i=0;i<2;i++) {
    for (int j=0;j<2;j++) {
      int n=ix+i+(iy+j)*m_nx;
      int m=i+2*j;
      sumw +=w[m];
      *resp += m_resp[n]*w[m];
    }
  }
  if (sumw>0.) *resp = *resp/sumw;
  else *resp = 1;
 
}

Initialize()

void MapEta::Initialize

(

int

isampling

)

Definition at line 43 of file MapEta.cxx.

{
  if (m_init==1) return;
  if (isampling < 1 || isampling >3) return;
  std::string filename{};
  if (isampling==1) filename = "deta_strip.map";
  if (isampling==2) filename = "deta_middle.map";
  if (isampling==3) filename = "eta_trans.map";
  std::string fileLocation;
#ifdef LARG4_STAND_ALONE
  // The stand-alone program expects to find the file via AFS.
  fileLocation = m_directory + "/" + filename;
#else
  // In Athena, the PathResolver tool will find the file for us.
  std::cout << "filename " << filename << std::endl;
  //std::string larLocation = PathResolver::find_directory("lar","DATAPATH");
  //fileLocation=larLocation+"/calo_data/"+filename;
  std::string larLocation = PathResolver::find_directory("LArG4Barrel","ATLASCALDATA");
  fileLocation=larLocation+"/"+filename;
  std::cout << "fileLocation " << fileLocation << std::endl;
#endif
 
  std::ifstream in(fileLocation);
  if (in)
    {
#ifndef LARG4_STAND_ALONE
      // Tell clang to optimize assuming that FP operations may trap.
      CXXUTILS_TRAPPING_FP;
#endif
 
      in>>m_nx>>m_xmin>>m_xmax>>m_ny>>m_ymin>>m_ymax;
      if(m_nx>0 && m_ny>0 && m_nx<10000 && m_ny<10000){//coverity issue. This is a tainted variable protection, 10000 can be changed if required.
        m_deltax=(m_xmax-m_xmin)/((float) m_nx);
        m_deltay=(m_ymax-m_ymin)/((float) m_ny);
        // what is written as xmax in the map is x of last point + delta x
        // (such that npoint = (xmax-xmin)/deltax
        // to get the real last point in the map we should subtract deltax
        m_xmax = m_xmax - m_deltax;
        m_ymax = m_ymax - m_deltay;
        m_resp = new float[m_ny*m_nx];
        m_xt0  = new float[m_ny*m_nx];
        m_xt1  = new float[m_ny*m_nx];
        m_xt2  = new float[m_ny*m_nx];
        int ii,jj;
        for (int ix=0;ix<m_nx;ix++) {
          for (int iy=0;iy<m_ny;iy++) {
            in>>ii>>jj>>m_resp[ix+iy*m_nx]>>m_xt0[ix+iy*m_nx]>>m_xt1[ix+iy*m_nx]>>m_xt2[ix+iy*m_nx];
            if (ii != ix || jj != iy) std::cout << "MapEta: inconsistency when reading map..." << std::endl;
          }
        }
 
        std::cout << "Eta Map read " << " Nbins " << m_nx << " " << m_ny
                  << " bin size*1000 " << 1000*m_deltax << " " << 1000*m_deltay
                  << " X range  " << m_xmin   << " " << m_xmax
                  << " Y range  " << m_ymin   << " " << m_ymax << std::endl;
        m_init=1;
      }
      else{
        std::cout << "Error in MapEta::Initialize: nx or ny out of limits." << std::endl;
        return;
      }
    }
  else{
    std::cout << "Error in MapEta::Initialize: The file could not be open." << std::endl;
    return;
  }
}

Nx()

int MapEta::Nx ( ) const

inline

Definition at line 28 of file MapEta.h.

{return m_nx;}

Ny()

int MapEta::Ny ( ) const

inline

Definition at line 29 of file MapEta.h.

{return m_ny;}

SetDirectory()

void MapEta::SetDirectory

(

const std::string &

dir

)

Definition at line 38 of file MapEta.cxx.

{
  m_directory=dir;
}

Xmax()

float MapEta::Xmax ( ) const

inline

Definition at line 31 of file MapEta.h.

{return m_xmax;}

Xmin()

float MapEta::Xmin ( ) const

inline

Definition at line 30 of file MapEta.h.

{return m_xmin;}

Ymax()

float MapEta::Ymax ( ) const

inline

Definition at line 33 of file MapEta.h.

{return m_ymax;}

Ymin()

float MapEta::Ymin ( ) const

inline

Definition at line 32 of file MapEta.h.

{return m_ymin;}

Member Data Documentation

m_deltax

float MapEta::m_deltax = 0.0F

private

Definition at line 18 of file MapEta.h.

m_deltay

float MapEta::m_deltay = 0.0F

private

Definition at line 18 of file MapEta.h.

m_directory

std::string MapEta::m_directory

private

Definition at line 14 of file MapEta.h.

m_init

int MapEta::m_init

private

Definition at line 15 of file MapEta.h.

m_nx

int MapEta::m_nx

private

Definition at line 16 of file MapEta.h.

m_ny

int MapEta::m_ny

private

Definition at line 16 of file MapEta.h.

m_resp

float* MapEta::m_resp

private

Definition at line 19 of file MapEta.h.

m_xmax

float MapEta::m_xmax = 0.0F

private

Definition at line 17 of file MapEta.h.

m_xmin

float MapEta::m_xmin = 0.0F

private

Definition at line 17 of file MapEta.h.

m_xt0

float* MapEta::m_xt0

private

Definition at line 20 of file MapEta.h.

m_xt1

float* MapEta::m_xt1

private

Definition at line 21 of file MapEta.h.

m_xt2

float* MapEta::m_xt2

private

Definition at line 22 of file MapEta.h.

m_ymax

float MapEta::m_ymax = 0.0F

private

Definition at line 17 of file MapEta.h.

m_ymin

float MapEta::m_ymin = 0.0F

private

Definition at line 17 of file MapEta.h.

---

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

• MapEta.h
• MapEta.cxx