ecalLightMap.h

Go to the documentation of this file.

#ifndef ECALLIGHTMAP_H
#define ECALLIGHTMAP_H
 
#include "TNamed.h"
 
class ecalLightMap : public TNamed
{
public:
 ecalLightMap() : TNamed(), fSSide(0.), fS(0), fSize(0), fData(NULL) {};
  ecalLightMap(const char* fname, const char* title="Light collection efficiency map");
  Double_t Data(Double_t x, Double_t y)
    {Int_t n=GetNum(x,y); if (n<0) return n; return fData[n];}
  Int_t GetNum(Double_t x, Double_t y)
  {
    Double_t lx=x+0.5; Double_t ly=y+0.5;
    if (lx<0) lx=0; if (ly<0) ly=0;
    if (lx>1) lx=1; if (ly>1) ly=1;
    Int_t ix=(Int_t)(lx/fSSide);
    Int_t iy=(Int_t)(ly/fSSide);
    return iy*fS+ix;
  }
  virtual ~ecalLightMap() {delete fData;}
private:
  void Init(const char* filename);
  void FillGaps();
  void Normalize();
  Double_t fSSide;      
  Int_t fS;         
  Int_t fSize;          
  Double_t* fData;      
 
  ecalLightMap(const ecalLightMap&);
  ecalLightMap& operator=(const ecalLightMap&);
 
  ClassDef(ecalLightMap, 1)
};
 
#endif