hcalStructure Class Reference

#include <hcalStructure.h>

Inheritance diagram for hcalStructure:

Collaboration diagram for hcalStructure:

Public Member Functions
	hcalStructure (hcalInf *hcalinf)
void	SetUseMC (Int_t mc=0)
Int_t	GetUseMC () const
void	Construct ()
Int_t	GetNumber (Int_t x, Int_t y) const
Bool_t	AddEnergy (Float_t x, Float_t y, Float_t energy, Float_t energy2)
Float_t	GetEnergy (Float_t x, Float_t y, Int_t section=0) const
hcalModule *	GetModule (Float_t x, Float_t y) const
Int_t	GetModuleNumber (Float_t x, Float_t y) const
Float_t	GetX1 () const
Float_t	GetY1 () const
Float_t	GetX2 () const
Float_t	GetY2 () const
hcalInf *	GetHcalInf () const
void	GetStructure (std::vector< hcalModule * > &stru) const
void	GetModules (std::list< hcalModule * > &mdls) const
void	ResetModules ()
hcalModule *	CreateModule (char type, Int_t number, Float_t x1, Float_t y1, Float_t x2, Float_t y2)
Int_t	GetHitId (Float_t x, Float_t y) const
void	GetHitXY (const Int_t hitId, Float_t &x, Float_t &y) const
hcalModule *	GetModule (Int_t fVolId, Int_t &section)
hcalModule *	GetHitModule (const Int_t hitId) const

Private Member Functions
Int_t	GetNum (Int_t x, Int_t y) const
void	Serialize ()
	hcalStructure (const hcalStructure &)
hcalStructure &	operator= (const hcalStructure &)
	ClassDef (hcalStructure, 1)

Private Attributes
Int_t	fUseMC
Float_t	fX1
Float_t	fY1
hcalInf *	fHcalInf
std::vector< hcalModule * >	fStructure
std::list< hcalModule * >	fModules
std::vector< hcalModule * >	fHash

Detailed Description

Definition at line 20 of file hcalStructure.h.

Constructor & Destructor Documentation

hcalStructure() [1/2]

hcalStructure::hcalStructure

(

hcalInf *

hcalinf

)

Definition at line 69 of file hcalStructure.cxx.

  : TNamed("hcalStructure", "Hadron calorimeter structure"), 
    fUseMC(0),
    fX1(0.),
    fY1(0.),
    fHcalInf(hcalinf),
    fStructure(),
    fModules(),
    fHash()
{
  fX1=fHcalInf->GetXPos()-fHcalInf->GetModuleSize()*fHcalInf->GetXSize()/2.0;
  fY1=fHcalInf->GetYPos()-fHcalInf->GetModuleSize()*fHcalInf->GetYSize()/2.0;
}

hcalStructure() [2/2]

hcalStructure::hcalStructure ( const hcalStructure &  )

private

Member Function Documentation

AddEnergy()

Bool_t hcalStructure::AddEnergy ( Float_t  x, Float_t  y, Float_t  energy, Float_t  energy2  )

inline

Add preshower or calorimeter energy to a module with coordinate (x,y)

Definition at line 125 of file hcalStructure.h.

{
  hcalModule* mdl=GetModule(x,y);
  if (mdl)
  {
    mdl->AddEnergy(energy);
    mdl->AddEnergy2(energy2);
  }
  else
    return kFALSE;
  return kTRUE;
}

ClassDef()

hcalStructure::ClassDef ( hcalStructure  , 1    )

private

Construct()

void hcalStructure::Construct

(

)

Definition at line 84 of file hcalStructure.cxx.

{
  if (!fHcalInf) return;
 
  Float_t x1=GetX1();
  Float_t y1=GetY1();
  Float_t x;
  Float_t y;
  Float_t dx;
  Float_t dy;
  Int_t i;
  Int_t j;
  Int_t k;
  Int_t number;
  char type;
 
  fStructure.resize(fHcalInf->GetXSize()*fHcalInf->GetYSize(), NULL);
    
  dx=fHcalInf->GetModuleSize();
  dy=fHcalInf->GetModuleSize();
  //Creating ECAL Matrix
  for(i=0;i<fHcalInf->GetXSize();i++)
    for(j=0;j<fHcalInf->GetYSize();j++) {
      type=fHcalInf->GetType(i,j);
      if (type) {
    x=x1+i*dx;
    y=y1+j*dy;
    number=(i*100+j)*100;
    fStructure[GetNum(i,j)]=CreateModule(type,number,x,y,x+dx,y+dy);
      }
      else
    fStructure[GetNum(i,j)]=NULL;
    }
#ifdef _DECALSTRUCT
  Info("Construct()", "Calorimeter matrix created.");
#endif
  //Now HCAL matrix created
  list<hcalModule*> neib;
  vector<hcalModule*> cl;
  vector<hcalModule*>::const_iterator pcl;
  
  Int_t num;
  //We want neighbors for ecalModules be ecalModules
  for(i=0;i<fHcalInf->GetXSize();i++)
    for(j=0;j<fHcalInf->GetYSize();j++)
      if (fStructure[GetNum(i,j)]) {
    neib.clear();
    
    num=GetNumber(i-1,j);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
            
    num=GetNumber(i-1,j+1);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
    
    num=GetNumber(i,j+1);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
    
    num=GetNumber(i+1,j+1);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
    
    num=GetNumber(i+1,j);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
            
    num=GetNumber(i+1,j-1);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
    
    num=GetNumber(i,j-1);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
    
    num=GetNumber(i-1,j-1);
    if (-1!=num) {
      neib.push_back(fStructure[num]);
    }
    
    num=GetNumber(i,j);
    fStructure[num]->SetNeighborsList(neib);
      }
  Serialize();
}

CreateModule()

hcalModule * hcalStructure::CreateModule	(	char	type,
		Int_t	number,
		Float_t	x1,
		Float_t	y1,
		Float_t	x2,
		Float_t	y2
	)

Definition at line 55 of file hcalStructure.cxx.

{
  if (type!=1)
  {
    Fatal("CreateModule", "All modules in hcal should have type 1");
    return NULL;
  }
  if (fUseMC)
    return new hcalModuleMC(number, x1, y1, x2, y2);
  else
    return new hcalModule(number, x1, y1, x2, y2);
}

GetEnergy()

Float_t hcalStructure::GetEnergy ( Float_t  x, Float_t  y, Int_t  section = 0  ) const

inline

Definition at line 139 of file hcalStructure.h.

{
  hcalModule* mdl=GetModule(x,y);
  if (mdl)
  {
    if (section==1)
      return mdl->GetEnergy2();
    else
      return mdl->GetEnergy();
  }
  return -1111;
}

GetHcalInf()

hcalInf * hcalStructure::GetHcalInf ( ) const

inline

Definition at line 38 of file hcalStructure.h.

{return fHcalInf;}

GetHitId()

Int_t hcalStructure::GetHitId ( Float_t  x, Float_t  y  ) const

inline

Definition at line 153 of file hcalStructure.h.

{
  hcalModule* mdl=GetModule(x,y);
  if (mdl)
    return mdl->GetNumber();
  else
    return -1111;
}

GetHitModule()

hcalModule * hcalStructure::GetHitModule

(

const Int_t

hitId

)

const

Hit Id -> Cell

Definition at line 215 of file hcalStructure.cxx.

{
  // Some translation from x*100+y  to y*sizex+x coding...
 
  Int_t mnum=hitId/10;
  Int_t cellx = mnum/100;
  Int_t celly = mnum%100;
  mnum = GetNum(cellx, celly);
  
  // end translation
 
  return fStructure[mnum];
}

GetHitXY()

void hcalStructure::GetHitXY	(	const Int_t	hitId,
		Float_t &	x,
		Float_t &	y
	)		const

Hit Id -> (x,y)

Definition at line 195 of file hcalStructure.cxx.

{
  // Some translation from x*100+y  to y*sizex+x coding...
 
  Int_t mnum=hitId/10;
  Int_t cellx = mnum/100;
  Int_t celly = mnum%100;
  mnum = GetNum(cellx, celly);
  
  // end translation
 
  hcalModule* module=fStructure[mnum];
  if (module==NULL) return;
  x=module->GetCenterX();
  y=module->GetCenterY();
}

GetModule() [1/2]

hcalModule * hcalStructure::GetModule ( Float_t  x, Float_t  y  ) const

inline

get ECAL module by known module center coordinate (x,y)

Definition at line 82 of file hcalStructure.h.

{
  Int_t num=GetModuleNumber(x,y);
  if (-1==num) return NULL; else return fStructure[num]; 
}

GetModule() [2/2]

hcalModule * hcalStructure::GetModule	(	Int_t	fVolId,
		Int_t &	section
	)

Definition at line 20 of file hcalStructure.cxx.

{
  UInt_t i;
  static Int_t volidmax = 0;
  volidmax=10000000;
 
  if ((Int_t)fHash.size()<volidmax)
  {
    fHash.resize(volidmax);
    for(i=0;i<fHash.size();i++)
      fHash[i]=NULL;
  }
  if (volId>volidmax)
    return NULL;
  if (fHash[volId]==NULL)
  {
    Float_t x;
    Float_t y;
    hcal::GetCellCoordInf(volId, x, y, section);
    fHash[volId]=GetModule(x+0.025,y+0.025);
  }
  section=volId%10;
  return fHash[volId];
}

GetModuleNumber()

Int_t hcalStructure::GetModuleNumber ( Float_t  x, Float_t  y  ) const

inline

get ECAL module by known module center coordinate (x,y)

Definition at line 89 of file hcalStructure.h.

{
  Int_t ix=(Int_t)TMath::Floor((x-GetX1())/fHcalInf->GetModuleSize());
  Int_t iy=(Int_t)TMath::Floor((y-GetY1())/fHcalInf->GetModuleSize());
  return GetNumber(ix,iy);
}

GetModules()

void hcalStructure::GetModules ( std::list< hcalModule * > &  mdls ) const

inline

Definition at line 40 of file hcalStructure.h.

{mdls=fModules;}

GetNum()

Int_t hcalStructure::GetNum ( Int_t  x, Int_t  y  ) const

inlineprivate

get ECAL absolute module number by known module relative number (x,y)

Definition at line 107 of file hcalStructure.h.

{
  return y*fHcalInf->GetXSize()+x;
}

GetNumber()

Int_t hcalStructure::GetNumber ( Int_t  x, Int_t  y  ) const

inline

get ECAL absolute module number by known module relative number (x,y) with check for the ECAL boundaries

Definition at line 97 of file hcalStructure.h.

{
  if (x>-1&&y>-1)
    if (x<fHcalInf->GetXSize()&&y<fHcalInf->GetYSize())
      return GetNum(x,y);
  return -1;
}

GetStructure()

void hcalStructure::GetStructure ( std::vector< hcalModule * > &  stru ) const

inline

Definition at line 39 of file hcalStructure.h.

{stru=fStructure;}

GetUseMC()

Int_t hcalStructure::GetUseMC ( ) const

inline

Definition at line 25 of file hcalStructure.h.

{return fUseMC;}

GetX1()

Float_t hcalStructure::GetX1 ( ) const

inline

Definition at line 34 of file hcalStructure.h.

{return fX1;};

GetX2()

Float_t hcalStructure::GetX2 ( ) const

inline

get ECAL right edge coordinate in cm

Definition at line 113 of file hcalStructure.h.

{
  return fHcalInf->GetXPos()+fHcalInf->GetModuleSize()*fHcalInf->GetXSize()/2.0;
}

GetY1()

Float_t hcalStructure::GetY1 ( ) const

inline

Definition at line 35 of file hcalStructure.h.

{return fY1;};

GetY2()

Float_t hcalStructure::GetY2 ( ) const

inline

get ECAL upper edge coordinate in cm

Definition at line 119 of file hcalStructure.h.

{
  return fHcalInf->GetYPos()+fHcalInf->GetModuleSize()*fHcalInf->GetYSize()/2.0;
}

operator=()

hcalStructure & hcalStructure::operator= ( const hcalStructure &  )

private

ResetModules()

void hcalStructure::ResetModules

(

)

Definition at line 179 of file hcalStructure.cxx.

{
  list<hcalModule*>::const_iterator p=fModules.begin();
  if (fUseMC==0)
  {
    for(;p!=fModules.end();++p)
      (*p)->ResetEnergyFast();
  }
  else
  {
    for(;p!=fModules.end();++p)
    ((hcalModuleMC*)(*p))->ResetEnergy();
  }
}

Serialize()

void hcalStructure::Serialize ( )

private

Creates modules lists

Definition at line 46 of file hcalStructure.cxx.

{
  fModules.clear();
  for(UInt_t i=0;i<fStructure.size();i++)
  if (fStructure[i])
    fModules.push_back(fStructure[i]);
}

SetUseMC()

void hcalStructure::SetUseMC ( Int_t  mc = 0 )

inline

Definition at line 24 of file hcalStructure.h.

{fUseMC=mc;}

Member Data Documentation

fHash

std::vector<hcalModule*> hcalStructure::fHash

private

MCPoint id -> HCAL mdoule

Definition at line 74 of file hcalStructure.h.

fHcalInf

hcalInf* hcalStructure::fHcalInf

private

HCAL geometry container

Definition at line 68 of file hcalStructure.h.

fModules

std::list<hcalModule*> hcalStructure::fModules

private

All ECAL modules

Definition at line 72 of file hcalStructure.h.

fStructure

std::vector<hcalModule*> hcalStructure::fStructure

private

vector of HCAL modules

Definition at line 70 of file hcalStructure.h.

fUseMC

Int_t hcalStructure::fUseMC

private

Use store MC information in modules

Definition at line 62 of file hcalStructure.h.

fX1

Float_t hcalStructure::fX1

private

X coordibate of left bottom angle of ECAL

Definition at line 64 of file hcalStructure.h.

fY1

Float_t hcalStructure::fY1

private

Y coordibate of left bottom angle of ECAL

Definition at line 66 of file hcalStructure.h.

---

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

• hcal/hcalStructure.h
• hcal/hcalStructure.cxx