ecalLightMap Class Reference

#include <ecalLightMap.h>

Inheritance diagram for ecalLightMap:

Collaboration diagram for ecalLightMap:

Public Member Functions
	ecalLightMap ()
	ecalLightMap (const char *fname, const char *title="Light collection efficiency map")
Double_t	Data (Double_t x, Double_t y)
Int_t	GetNum (Double_t x, Double_t y)
virtual	~ecalLightMap ()

Private Member Functions
void	Init (const char *filename)
void	FillGaps ()
void	Normalize ()
	ecalLightMap (const ecalLightMap &)
ecalLightMap &	operator= (const ecalLightMap &)

Private Attributes
Double_t	fSSide
Int_t	fS
Int_t	fSize
Double_t *	fData

Detailed Description

Definition at line 6 of file ecalLightMap.h.

Constructor & Destructor Documentation

ecalLightMap() [1/3]

ecalLightMap::ecalLightMap ( )

inline

Definition at line 9 of file ecalLightMap.h.

: TNamed(), fSSide(0.), fS(0), fSize(0), fData(NULL) {};

ecalLightMap() [2/3]

ecalLightMap::ecalLightMap	(	const char *	fname,
		const char *	title = "Light collection efficiency map"
	)

Definition at line 13 of file ecalLightMap.cxx.

  : TNamed(fname, title),
    fSSide(-1111.), fS(-1111), fSize(0), fData(NULL)
{
  Init(fname);
}

~ecalLightMap()

virtual ecalLightMap::~ecalLightMap ( )

inlinevirtual

Definition at line 22 of file ecalLightMap.h.

{delete fData;}

ecalLightMap() [3/3]

ecalLightMap::ecalLightMap ( const ecalLightMap &  )

private

Member Function Documentation

Data()

Double_t ecalLightMap::Data ( Double_t  x, Double_t  y  )

inline

Definition at line 11 of file ecalLightMap.h.

    {Int_t n=GetNum(x,y); if (n<0) return n; return fData[n];}

FillGaps()

void ecalLightMap::FillGaps ( )

private

Fix a light collection map

Definition at line 99 of file ecalLightMap.cxx.

{
  Int_t i;
  Int_t j;
  Double_t x;
  Double_t y;
  Int_t n[3];
  Double_t v;
  Int_t d;
 
  for(i=0;i<fSize;i++)
  if (fData[i]<0)
  {
    x=(i%fS)*fSSide-0.5+fSSide/2.0;
    y=(i/fS)*fSSide-0.5+fSSide/2.0;
    n[0]=GetNum( x,-y);
    n[1]=GetNum(-x, y);
    n[2]=GetNum(-x,-y);
    v=0; d=0;
    for(j=0;j<3;j++)
    if (n[j]>=0&&fData[n[j]]>=0)
      { d++; v+=fData[n[j]]; }
    if (d>0)
    {
      v/=d; fData[i]=v;
      if (fData[n[0]]<=0) fData[n[0]]=v;
      if (fData[n[1]]<=0) fData[n[1]]=v;
      if (fData[n[2]]<=0) fData[n[2]]=v;
    }
    else
      Info("FillGaps","No data for (%f,%f)", x, y);
  }
}

GetNum()

Int_t ecalLightMap::GetNum ( Double_t  x, Double_t  y  )

inline

Definition at line 13 of file ecalLightMap.h.

  {
    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;
  }

Init()

void ecalLightMap::Init ( const char *  filename )

private

Read information from file

Definition at line 21 of file ecalLightMap.cxx.

{
  TString fn=filename;
  gSystem->ExpandPathName(fn);
  ifstream f(fn);
  list<Double_t> lst;
  string buf;
  string token;
  string message;
  if (!f)
  {
    Fatal("Init","Can't open file %s.", filename);
    return;
  }
 
  while(getline(f, buf))
  {
    if (buf.empty()) continue;
    if (buf.find_first_not_of("     ")==string::npos) continue;
    //Skipping initial spaces
    message=buf.substr(buf.find_first_not_of("  "));
    if (message.empty()) continue;
    //Removing comments
    message=message.substr(0,message.find("#"));
    if (message.empty()) continue;
    while(!message.empty())
    {
      token=message.substr(0,message.find_first_of("    "));
      if (token.empty()) break;
      lst.push_back(atof(token.c_str()));
      if (token==message) break;
      token=message.substr(message.find_first_of("  "));
      if (token.empty()) break;
      if (token.find_first_not_of("     ")==string::npos) break;
      message=token.substr(token.find_first_not_of("    "));
    }
  }
  f.close();
 
  list<Double_t>::const_iterator p=lst.begin();
  Double_t xsize;
  Double_t ysize;
  Double_t sqside;
  Double_t l;
  Double_t x;
  Double_t y;
  Double_t z;
  Double_t v;
  Int_t n;
  Int_t i;
 
  xsize=(*p); ++p; ysize=(*p); ++p; sqside=(*p); ++p;
  fS=(Int_t)((xsize+0.00001)/sqside);
  fSize=fS*fS; fSSide=sqside/xsize;
  fData=new Double_t[fSize];
  for(i=0;i<fSize;i++)
    fData[i]=-1111;
  for(;;)
  {
    if (p==lst.end()) break; l=(*p); ++p;
    if (p==lst.end()) break; x=(*p); ++p;
    if (p==lst.end()) break; y=(*p); ++p;
    if (p==lst.end()) break; z=(*p); ++p;
    if (p==lst.end()) break; v=(*p); ++p;
    x/=xsize; y/=ysize;
    n=GetNum(x ,y);
    if (n>=fSize)
    {
      Info("Init","Data is not selfconsistent (%f, %f), %d", x, y, fSize);
    }
    fData[n]=v;
  }
  lst.clear();
  FillGaps();
  Normalize();
}

Normalize()

void ecalLightMap::Normalize ( )

private

Set average efficiency of light collection to 1.0

Definition at line 134 of file ecalLightMap.cxx.

{
  Int_t i;
  Int_t n=0;
  Double_t v=0;
  for(i=0;i<fSize;i++)
    if (fData[i]>0)
      { v+=fData[i]; n++; }
  v/=n;
  for(i=0;i<fSize;i++)
    if (fData[i]>0)
      fData[i]/=v;
}

operator=()

ecalLightMap & ecalLightMap::operator= ( const ecalLightMap &  )

private

Member Data Documentation

fData

Double_t* ecalLightMap::fData

private

Light collection efficiency map

Definition at line 37 of file ecalLightMap.h.

fS

Int_t ecalLightMap::fS

private

Size of side of the light map in steps

Definition at line 33 of file ecalLightMap.h.

fSize

Int_t ecalLightMap::fSize

private

Size of the light map

Definition at line 35 of file ecalLightMap.h.

fSSide

Double_t ecalLightMap::fSSide

private

Step of the light map

Definition at line 31 of file ecalLightMap.h.

---

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

• ecal/ecalLightMap.h
• ecal/ecalLightMap.cxx