#include <CaloGeometryLookup.h>

Collaboration diagram for CaloGeometryLookup:

Public Member Functions
	CaloGeometryLookup (int ind=0)

virtual	~CaloGeometryLookup ()

bool	IsCompatible (const CaloDetDescrElement *cell)

void	add (const CaloDetDescrElement *cell)

t_cellmap::size_type	size () const

int	index () const

void	set_index (int ind)

void	post_process ()

bool	has_overlap (CaloGeometryLookup *ref)

void	merge_into_ref (CaloGeometryLookup *ref)

float	mineta () const

float	maxeta () const

float	minphi () const

float	maxphi () const

float	mineta_raw () const

float	maxeta_raw () const

float	minphi_raw () const

float	maxphi_raw () const

float	minx () const

float	maxx () const

float	miny () const

float	maxy () const

float	minx_raw () const

float	maxx_raw () const

float	miny_raw () const

float	maxy_raw () const

const MeanAndRMS &	deta ()

const MeanAndRMS &	dphi ()

float	mindeta ()

float	mindphi ()

const MeanAndRMS &	dx ()

const MeanAndRMS &	dy ()

float	mindx ()

float	mindy ()

const MeanAndRMS &	eta_correction ()

const MeanAndRMS &	phi_correction ()

const MeanAndRMS &	x_correction ()

const MeanAndRMS &	y_correction ()

int	cell_grid_eta () const

int	cell_grid_phi () const

void	set_xy_grid_adjustment_factor (float factor)

virtual const CaloDetDescrElement *	getDDE (float eta, float phi, float distance=0, int steps=0)

Protected Member Functions
float	neta_double ()

float	nphi_double ()

Int_t	neta ()

Int_t	nphi ()

float	nx_double ()

float	ny_double ()

Int_t	nx ()

Int_t	ny ()

int	raw_eta_position_to_index (float eta_raw) const

int	raw_phi_position_to_index (float phi_raw) const

bool	index_range_adjust (int &ieta, int &iphi) const

float	calculate_distance_eta_phi (const CaloDetDescrElement *DDE, float eta, float phi, float &dist_eta0, float &dist_phi0) const

Protected Attributes
float	m_xy_grid_adjustment_factor

int	m_index

t_cellmap	m_cells

std::vector< std::vector< const CaloDetDescrElement * > >	m_cell_grid

int	m_cell_grid_eta

int	m_cell_grid_phi

float	m_mineta

float	m_maxeta

float	m_minphi

float	m_maxphi

float	m_mineta_raw

float	m_maxeta_raw

float	m_minphi_raw

float	m_maxphi_raw

float	m_mineta_correction

float	m_maxeta_correction

float	m_minphi_correction

float	m_maxphi_correction

MeanAndRMS	m_deta

MeanAndRMS	m_dphi

MeanAndRMS	m_eta_correction

MeanAndRMS	m_phi_correction

float	m_mindeta

float	m_mindphi

float	m_deta_double

float	m_dphi_double

Detailed Description

Definition at line 25 of file CaloGeometryLookup.h.

Constructor & Destructor Documentation

◆ CaloGeometryLookup()

CaloGeometryLookup::CaloGeometryLookup ( int ind = 0 )

Definition at line 24 of file CaloGeometryLookup.cxx.

                                              :m_xy_grid_adjustment_factor(0.75),m_index(ind)
 {
   m_mineta=+10000;
   m_maxeta=-10000;
   m_minphi=+10000;
   m_maxphi=-10000;
  
   m_mineta_raw=+10000;
   m_maxeta_raw=-10000;
   m_minphi_raw=+10000;
   m_maxphi_raw=-10000;
  
   m_mindeta=10000;
   m_mindphi=10000;
   
   m_mineta_correction=+10000;
   m_maxeta_correction=-10000;
   m_minphi_correction=+10000;
   m_maxphi_correction=-10000;
  
   m_cell_grid_eta=0.;
   m_cell_grid_phi=0.;
   m_deta_double  =0.;
   m_dphi_double  =0.;
 }

◆ ~CaloGeometryLookup()

CaloGeometryLookup::~CaloGeometryLookup ( )

virtualdefault

Member Function Documentation

◆ add()

void CaloGeometryLookup::add ( const CaloDetDescrElement * cell )

Definition at line 122 of file CaloGeometryLookup.cxx.

 {
   if(cell->getSampling()<21) {
     m_deta.add(cell->deta());
     m_dphi.add(cell->dphi());
     m_mindeta=TMath::Min(cell->deta(),m_mindeta);
     m_mindphi=TMath::Min(cell->dphi(),m_mindphi);
  
     m_eta_correction.add(cell->eta_raw()-cell->eta());
     m_phi_correction.add(cell->phi_raw()-cell->phi());
     m_mineta_correction=TMath::Min(cell->eta_raw()-cell->eta() , m_mineta_correction);
     m_maxeta_correction=TMath::Max(cell->eta_raw()-cell->eta() , m_maxeta_correction);
     m_minphi_correction=TMath::Min(cell->phi_raw()-cell->phi() , m_minphi_correction);
     m_maxphi_correction=TMath::Max(cell->phi_raw()-cell->phi() , m_maxphi_correction);
  
     m_mineta=TMath::Min(cell->eta()-cell->deta()/2,m_mineta);
     m_maxeta=TMath::Max(cell->eta()+cell->deta()/2,m_maxeta);
     m_minphi=TMath::Min(cell->phi()-cell->dphi()/2,m_minphi);
     m_maxphi=TMath::Max(cell->phi()+cell->dphi()/2,m_maxphi);
  
     m_mineta_raw=TMath::Min(cell->eta_raw()-cell->deta()/2,m_mineta_raw);
     m_maxeta_raw=TMath::Max(cell->eta_raw()+cell->deta()/2,m_maxeta_raw);
     m_minphi_raw=TMath::Min(cell->phi_raw()-cell->dphi()/2,m_minphi_raw);
     m_maxphi_raw=TMath::Max(cell->phi_raw()+cell->dphi()/2,m_maxphi_raw);
   } else {
     m_deta.add(cell->dx());
     m_dphi.add(cell->dy());
     m_mindeta=TMath::Min(cell->dx(),m_mindeta);
     m_mindphi=TMath::Min(cell->dy(),m_mindphi);
  
     m_eta_correction.add(cell->x_raw()-cell->x());
     m_phi_correction.add(cell->y_raw()-cell->y());
     m_mineta_correction=TMath::Min(cell->x_raw()-cell->x() , m_mineta_correction);
     m_maxeta_correction=TMath::Max(cell->x_raw()-cell->x() , m_maxeta_correction);
     m_minphi_correction=TMath::Min(cell->y_raw()-cell->y() , m_minphi_correction);
     m_maxphi_correction=TMath::Max(cell->y_raw()-cell->y() , m_maxphi_correction);
  
     m_mineta=TMath::Min(cell->x()-cell->dx()/2,m_mineta);
     m_maxeta=TMath::Max(cell->x()+cell->dx()/2,m_maxeta);
     m_minphi=TMath::Min(cell->y()-cell->dy()/2,m_minphi);
     m_maxphi=TMath::Max(cell->y()+cell->dy()/2,m_maxphi);
  
     m_mineta_raw=TMath::Min(cell->x_raw()-cell->dx()/2,m_mineta_raw);
     m_maxeta_raw=TMath::Max(cell->x_raw()+cell->dx()/2,m_maxeta_raw);
     m_minphi_raw=TMath::Min(cell->y_raw()-cell->dy()/2,m_minphi_raw);
     m_maxphi_raw=TMath::Max(cell->y_raw()+cell->dy()/2,m_maxphi_raw);
   }  
   m_cells[cell->identify()]=cell;
 }

◆ calculate_distance_eta_phi()

float CaloGeometryLookup::calculate_distance_eta_phi	(	const CaloDetDescrElement *	DDE,
		float	eta,
		float	phi,
		float &	dist_eta0,
		float &	dist_phi0
	)		const

protected

Definition at line 249 of file CaloGeometryLookup.cxx.

 {
   dist_eta0=(eta - DDE->eta())/m_deta_double;
   dist_phi0=(TVector2::Phi_mpi_pi(phi - DDE->phi()))/m_dphi_double;
   float abs_dist_eta0=TMath::Abs(dist_eta0);
   float abs_dist_phi0=TMath::Abs(dist_phi0);
   return TMath::Max(abs_dist_eta0,abs_dist_phi0)-0.5;
 }

◆ cell_grid_eta()

int CaloGeometryLookup::cell_grid_eta ( ) const

inline

Definition at line 82 of file CaloGeometryLookup.h.

82 {return m_cell_grid_eta;};

◆ cell_grid_phi()

int CaloGeometryLookup::cell_grid_phi ( ) const

inline

Definition at line 83 of file CaloGeometryLookup.h.

83 {return m_cell_grid_phi;};

◆ deta()

const MeanAndRMS& CaloGeometryLookup::deta ( )

inline

Definition at line 68 of file CaloGeometryLookup.h.

68 {return m_deta;};

◆ dphi()

const MeanAndRMS& CaloGeometryLookup::dphi ( )

inline

Definition at line 69 of file CaloGeometryLookup.h.

69 {return m_dphi;};

◆ dx()

const MeanAndRMS& CaloGeometryLookup::dx ( )

inline

Definition at line 72 of file CaloGeometryLookup.h.

72 {return m_deta;};

◆ dy()

const MeanAndRMS& CaloGeometryLookup::dy ( )

inline

Definition at line 73 of file CaloGeometryLookup.h.

73 {return m_dphi;};

◆ eta_correction()

const MeanAndRMS& CaloGeometryLookup::eta_correction ( )

inline

Definition at line 77 of file CaloGeometryLookup.h.

77 {return m_eta_correction;};

◆ getDDE()

const CaloDetDescrElement * CaloGeometryLookup::getDDE	(	float	eta,
		float	phi,
		float *	distance = `0`,
		int *	steps = `0`
	)

virtual

Definition at line 258 of file CaloGeometryLookup.cxx.

 {
   float dist = 0;
   const CaloDetDescrElement* bestDDE=nullptr;
   if(!distance) distance=&dist;
   *distance=+10000000;
   int intsteps=0;
   if(!steps) steps=&intsteps;
   
   float best_eta_corr=m_eta_correction;
   float best_phi_corr=m_phi_correction;
   
   float raw_eta=eta+best_eta_corr;
   float raw_phi=phi+best_phi_corr;
  
   int ieta=raw_eta_position_to_index(raw_eta);
   int iphi=raw_phi_position_to_index(raw_phi);
   index_range_adjust(ieta,iphi);
   
   const CaloDetDescrElement* newDDE=m_cell_grid[ieta][iphi];
   float bestdist=+10000000;
   ++(*steps);
   int nsearch=0;
   while(newDDE && nsearch<3) {
     float dist_eta0,dist_phi0;
     *distance=calculate_distance_eta_phi(newDDE,eta,phi,dist_eta0,dist_phi0);
     bestDDE=newDDE;
     bestdist=*distance;
  
     if(CaloGeometry::m_debug || CaloGeometry::m_debug_identify==newDDE->identify()) {
       cout<<"CaloGeometryLookup::getDDE, search iteration "<<nsearch<<endl;
       cout<<"cell id="<<newDDE->identify()<<" steps "<<*steps<<" steps, eta="<<eta<<" phi="<<phi<<" dist="<<*distance<<" cell_grid_eta="<<cell_grid_eta()<<" cell_grid_phi="<<cell_grid_phi()<<" m_deta_double="<<m_deta_double<<" m_dphi_double="<<m_dphi_double<<" 1st_eta_corr="<<best_eta_corr<<" 1st_phi_corr="<<best_phi_corr<<endl;
       cout<<"  input sampling="<<newDDE->getSampling()<<" eta="<<newDDE->eta()<<" eta_raw="<<newDDE->eta_raw()<<" deta="<<newDDE->deta()<<" ("<<(newDDE->eta_raw()-newDDE->eta())/newDDE->deta()<<") phi="<<newDDE->phi()<<" phi_raw="<<newDDE->phi_raw()<<" dphi="<<newDDE->dphi()<<" ("<<(newDDE->phi_raw()-newDDE->phi())/newDDE->dphi()<<")"<<endl;
       cout<<"  ieta="<<ieta<<" iphi="<<iphi<<endl;
       cout<<"  dist_eta0="<<dist_eta0<<" ,dist_phi0="<<dist_phi0<<endl;
     }
  
     if(*distance<0) return newDDE;
     //correct ieta and iphi by the observed difference to the hit cell
     ieta+=TMath::Nint(dist_eta0);
     iphi+=TMath::Nint(dist_phi0);
     index_range_adjust(ieta,iphi);
     const CaloDetDescrElement* oldDDE=newDDE;
     newDDE=m_cell_grid[ieta][iphi];
     ++(*steps);
     ++nsearch;
     if(oldDDE==newDDE) break;
   } 
   if(CaloGeometry::m_debug && !newDDE) {
     cout<<"CaloGeometryLookup::getDDE, direct search ieta="<<ieta<<" iphi="<<iphi<<" is empty"<<endl;
   }
   float minieta=ieta+TMath::Nint(TMath::Floor(m_mineta_correction/cell_grid_eta()));
   float maxieta=ieta+TMath::Nint(TMath::Ceil(m_maxeta_correction/cell_grid_eta()));
   float miniphi=iphi+TMath::Nint(TMath::Floor(m_minphi_correction/cell_grid_phi()));
   float maxiphi=iphi+TMath::Nint(TMath::Ceil(m_maxphi_correction/cell_grid_phi()));
   if(minieta<0) minieta=0;
   if(maxieta>=m_cell_grid_eta) maxieta=m_cell_grid_eta-1;
   for(int iieta=minieta;iieta<=maxieta;++iieta) {
     for(int iiphi=miniphi;iiphi<=maxiphi;++iiphi) {
       ieta=iieta;
       iphi=iiphi;
       index_range_adjust(ieta,iphi);
       newDDE=m_cell_grid[ieta][iphi];
       ++(*steps);
       if(newDDE) {
         float dist_eta0,dist_phi0;
         *distance=calculate_distance_eta_phi(newDDE,eta,phi,dist_eta0,dist_phi0);
  
         if(CaloGeometry::m_debug || CaloGeometry::m_debug_identify==newDDE->identify()) {
           cout<<"CaloGeometryLookup::getDDE, windows search! minieta="<<m_mineta_correction/cell_grid_eta()<<" maxieta="<<m_maxeta_correction/cell_grid_eta()<<" miniphi="<<m_minphi_correction/cell_grid_phi()<<" maxiphi="<<m_maxphi_correction/cell_grid_phi()<<endl;
           cout<<"cell id="<<newDDE->identify()<<" steps "<<*steps<<" steps, eta="<<eta<<" phi="<<phi<<" dist="<<*distance<<endl;
           cout<<"  input sampling="<<newDDE->getSampling()<<" eta="<<newDDE->eta()<<" eta_raw="<<newDDE->eta_raw()<<" deta="<<newDDE->deta()<<" ("<<(newDDE->eta_raw()-newDDE->eta())/newDDE->deta()<<") phi="<<newDDE->phi()<<" phi_raw="<<newDDE->phi_raw()<<" dphi="<<newDDE->dphi()<<" ("<<(newDDE->phi_raw()-newDDE->phi())/newDDE->dphi()<<")"<<endl;
           cout<<"  ieta="<<ieta<<" iphi="<<iphi<<endl;
           cout<<"  dist_eta0="<<dist_eta0<<" ,dist_phi0="<<dist_phi0<<endl;
         }
  
         if(*distance<0) return newDDE;
         if(*distance<bestdist) {
           bestDDE=newDDE;
           bestdist=*distance;
         }
       } else if(CaloGeometry::m_debug) {
         cout<<"CaloGeometryLookup::getDDE, windows search ieta="<<ieta<<" iphi="<<iphi<<" is empty"<<endl;
       }
     }
   }
   *distance=bestdist;
   return bestDDE;
 }

◆ has_overlap()

bool CaloGeometryLookup::has_overlap ( CaloGeometryLookup * ref )

Definition at line 53 of file CaloGeometryLookup.cxx.

 {
   if(m_cells.empty()) return false;
   for(t_cellmap::iterator ic=m_cells.begin();ic!=m_cells.end();++ic) {
     const CaloDetDescrElement* cell=ic->second;
     if(ref->IsCompatible(cell)) return true;
   }
   return false;
 }

◆ index()

int CaloGeometryLookup::index ( ) const

inline

Definition at line 33 of file CaloGeometryLookup.h.

33 {return m_index;};

◆ index_range_adjust()

bool CaloGeometryLookup::index_range_adjust	(	int &	ieta,
		int &	iphi
	)		const

protected

Definition at line 172 of file CaloGeometryLookup.cxx.

 {
   while(iphi<0) {iphi+=m_cell_grid_phi;};
   while(iphi>=m_cell_grid_phi) {iphi-=m_cell_grid_phi;};
   if(ieta<0) {
     ieta=0;
     return false;
   }
   if(ieta>=m_cell_grid_eta) {
     ieta=m_cell_grid_eta-1;
     return false;
   }
   return true;
 }

◆ IsCompatible()

bool CaloGeometryLookup::IsCompatible ( const CaloDetDescrElement * cell )

Definition at line 72 of file CaloGeometryLookup.cxx.

 {
   if(m_cells.empty()) return true;
   t_cellmap::iterator ic=m_cells.begin();
   const CaloDetDescrElement* refcell=ic->second;
   int sampling=refcell->getSampling();
   if(cell->getSampling()!=sampling) return false;
   if(cell->eta_raw()*refcell->eta_raw()<0) return false;
   if(sampling<21) { // keep only cells reasonable close in eta to each other
     if(TMath::Abs(cell->deta()-refcell->deta())>0.001) return false;
     if(TMath::Abs(cell->dphi()-refcell->dphi())>0.001) return false;
     
     if( cell->eta()<mineta()-2.1*cell->deta() ) return false;  
     if( cell->eta()>maxeta()+2.1*cell->deta() ) return false;  
  
     double delta_eta=(cell->eta_raw()-refcell->eta_raw())/cell->deta();
     double delta_phi=(cell->phi_raw()-refcell->phi_raw())/cell->dphi();
     if(TMath::Abs(delta_eta-TMath::Nint(delta_eta)) > 0.01) return false;
     if(TMath::Abs(delta_phi-TMath::Nint(delta_phi)) > 0.01) return false;
  
     /*
     cout<<"Compatible: s="<<cell->getSampling()<<"~"<<refcell->getSampling()<<"; ";
     cout<<"eta="<<cell->eta_raw()<<","<<refcell->eta_raw()<<"; ";
     cout<<"phi="<<cell->phi_raw()<<","<<refcell->phi_raw()<<"; ";
     cout<<"deta="<<cell->deta()<<"~"<<refcell->deta()<<"; ";
     cout<<"dphi="<<cell->dphi()<<"~"<<refcell->dphi()<<";";
     cout<<"mineta="<<mineta()<<", maxeta="<<maxeta()<<"; ";
     cout<<endl;
     */
   } else {
     //FCal is not sufficiently regular to use submaps with regular mapping
     return true;
     /*
     if(TMath::Abs(cell->dx()-refcell->dx())>0.001) return false;
     if(TMath::Abs(cell->dy()-refcell->dy())>0.001) return false;
     if( cell->x()<minx()-20*cell->dx() ) return false;
     if( cell->x()>maxx()+20*cell->dx() ) return false;
     if( cell->y()<miny()-20*cell->dy() ) return false;
     if( cell->y()>maxy()+20*cell->dy() ) return false;
  
     double delta_x=(cell->x_raw()-refcell->x_raw())/cell->dx();
     double delta_y=(cell->y_raw()-refcell->y_raw())/cell->dy();
     if(TMath::Abs(delta_x-TMath::Nint(delta_x)) > 0.01) return false;
     if(TMath::Abs(delta_y-TMath::Nint(delta_y)) > 0.01) return false;
     */
   }  
  
   return true;  
 }

◆ maxeta()

float CaloGeometryLookup::maxeta ( ) const

inline

Definition at line 41 of file CaloGeometryLookup.h.

41 {return m_maxeta;};

◆ maxeta_raw()

float CaloGeometryLookup::maxeta_raw ( ) const

inline

Definition at line 46 of file CaloGeometryLookup.h.

46 {return m_maxeta_raw;};

◆ maxphi()

float CaloGeometryLookup::maxphi ( ) const

inline

Definition at line 43 of file CaloGeometryLookup.h.

43 {return m_maxphi;};

◆ maxphi_raw()

float CaloGeometryLookup::maxphi_raw ( ) const

inline

Definition at line 48 of file CaloGeometryLookup.h.

48 {return m_maxphi_raw;};

◆ maxx()

float CaloGeometryLookup::maxx ( ) const

inline

Definition at line 59 of file CaloGeometryLookup.h.

59 {return m_maxeta;};

◆ maxx_raw()

float CaloGeometryLookup::maxx_raw ( ) const

inline

Definition at line 64 of file CaloGeometryLookup.h.

64 {return m_maxeta_raw;};

◆ maxy()

float CaloGeometryLookup::maxy ( ) const

inline

Definition at line 61 of file CaloGeometryLookup.h.

61 {return m_maxphi;};

◆ maxy_raw()

float CaloGeometryLookup::maxy_raw ( ) const

inline

Definition at line 66 of file CaloGeometryLookup.h.

66 {return m_maxphi_raw;};

◆ merge_into_ref()

void CaloGeometryLookup::merge_into_ref ( CaloGeometryLookup * ref )

Definition at line 63 of file CaloGeometryLookup.cxx.

 {
   for(t_cellmap::iterator ic=m_cells.begin();ic!=m_cells.end();++ic) {
     const CaloDetDescrElement* cell=ic->second;
     ref->add(cell);
   }
 }

◆ mindeta()

float CaloGeometryLookup::mindeta ( )

inline

Definition at line 70 of file CaloGeometryLookup.h.

70 {return m_mindeta;};

◆ mindphi()

float CaloGeometryLookup::mindphi ( )

inline

Definition at line 71 of file CaloGeometryLookup.h.

71 {return m_mindphi;};

◆ mindx()

float CaloGeometryLookup::mindx ( )

inline

Definition at line 74 of file CaloGeometryLookup.h.

74 {return m_mindeta;};

◆ mindy()

float CaloGeometryLookup::mindy ( )

inline

Definition at line 75 of file CaloGeometryLookup.h.

75 {return m_mindphi;};

◆ mineta()

float CaloGeometryLookup::mineta ( ) const

inline

Definition at line 40 of file CaloGeometryLookup.h.

40 {return m_mineta;};

◆ mineta_raw()

float CaloGeometryLookup::mineta_raw ( ) const

inline

Definition at line 45 of file CaloGeometryLookup.h.

45 {return m_mineta_raw;};

◆ minphi()

float CaloGeometryLookup::minphi ( ) const

inline

Definition at line 42 of file CaloGeometryLookup.h.

42 {return m_minphi;};

◆ minphi_raw()

float CaloGeometryLookup::minphi_raw ( ) const

inline

Definition at line 47 of file CaloGeometryLookup.h.

47 {return m_minphi_raw;};

◆ minx()

float CaloGeometryLookup::minx ( ) const

inline

Definition at line 58 of file CaloGeometryLookup.h.

58 {return m_mineta;};

◆ minx_raw()

float CaloGeometryLookup::minx_raw ( ) const

inline

Definition at line 63 of file CaloGeometryLookup.h.

63 {return m_mineta_raw;};

◆ miny()

float CaloGeometryLookup::miny ( ) const

inline

Definition at line 60 of file CaloGeometryLookup.h.

60 {return m_minphi;};

◆ miny_raw()

float CaloGeometryLookup::miny_raw ( ) const

inline

Definition at line 65 of file CaloGeometryLookup.h.

65 {return m_minphi_raw;};

◆ neta()

Int_t CaloGeometryLookup::neta ( )

inlineprotected

Definition at line 99 of file CaloGeometryLookup.h.

99 {return TMath::Nint( neta_double() );};

◆ neta_double()

float CaloGeometryLookup::neta_double ( )

inlineprotected

Definition at line 97 of file CaloGeometryLookup.h.

97 {return (maxeta_raw()-mineta_raw())/deta().mean();};

◆ nphi()

Int_t CaloGeometryLookup::nphi ( )

inlineprotected

Definition at line 100 of file CaloGeometryLookup.h.

100 {return TMath::Nint( nphi_double() );};

◆ nphi_double()

float CaloGeometryLookup::nphi_double ( )

inlineprotected

Definition at line 98 of file CaloGeometryLookup.h.

98 {return (maxphi_raw()-minphi_raw())/dphi().mean();};

◆ nx()

Int_t CaloGeometryLookup::nx ( )

inlineprotected

Definition at line 105 of file CaloGeometryLookup.h.

105 {return TMath::Nint(TMath::Ceil( nx_double() ));};

◆ nx_double()

float CaloGeometryLookup::nx_double ( )

inlineprotected

Definition at line 103 of file CaloGeometryLookup.h.

103 {return (maxx_raw()-minx_raw())/mindx();};

◆ ny()

Int_t CaloGeometryLookup::ny ( )

inlineprotected

Definition at line 106 of file CaloGeometryLookup.h.

106 {return TMath::Nint(TMath::Ceil( ny_double() ));};

◆ ny_double()

float CaloGeometryLookup::ny_double ( )

inlineprotected

Definition at line 104 of file CaloGeometryLookup.h.

104 {return (maxy_raw()-miny_raw())/mindy();};

◆ phi_correction()

const MeanAndRMS& CaloGeometryLookup::phi_correction ( )

inline

Definition at line 78 of file CaloGeometryLookup.h.

78 {return m_phi_correction;};

◆ post_process()

void CaloGeometryLookup::post_process ( )

Definition at line 187 of file CaloGeometryLookup.cxx.

 {
   if(size()==0) return;
   t_cellmap::iterator ic=m_cells.begin();
   const CaloDetDescrElement* refcell=ic->second;
   int sampling=refcell->getSampling();
   if(sampling<21) {
     double rneta=neta_double()-neta();
     double rnphi=nphi_double()-nphi();
     if(TMath::Abs(rneta)>0.05 || TMath::Abs(rnphi)>0.05) {
       cout<<"ERROR: can't map cells into a regular grid for sampling "<<sampling<<", map index "<<index()<<endl;
       return;
     }
     m_cell_grid_eta=neta();
     m_cell_grid_phi=nphi();
     m_deta_double=deta();
     m_dphi_double=dphi();
   } else {
     //double rnx=nx_double()-nx();
     //double rny=ny_double()-ny();
     //if(TMath::Abs(rnx)>0.05 || TMath::Abs(rny)>0.05) {
     //  cout<<"Grid: Sampling "<<sampling<<"_"<<index()<<": mapping cells into a regular grid, although cells don't fit well"<<endl;
     //}
  
     m_cell_grid_eta=TMath::Nint(TMath::Ceil(nx_double()/m_xy_grid_adjustment_factor));
     m_cell_grid_phi=TMath::Nint(TMath::Ceil(ny_double()/m_xy_grid_adjustment_factor));
     m_deta_double=mindx()*m_xy_grid_adjustment_factor;
     m_dphi_double=mindy()*m_xy_grid_adjustment_factor;
   }
   
   m_cell_grid.resize(m_cell_grid_eta); 
   for(int ieta=0;ieta<m_cell_grid_eta;++ieta) {
     m_cell_grid[ieta].resize(m_cell_grid_phi,nullptr); 
   }  
   
   for(ic=m_cells.begin();ic!=m_cells.end();++ic) {
     refcell=ic->second;
     int ieta,iphi;
     if(sampling<21) {
       ieta=raw_eta_position_to_index(refcell->eta_raw());
       iphi=raw_phi_position_to_index(refcell->phi_raw());
     } else {
       ieta=raw_eta_position_to_index(refcell->x_raw());
       iphi=raw_phi_position_to_index(refcell->y_raw());
     }
     if(index_range_adjust(ieta,iphi)) {
       if(m_cell_grid[ieta][iphi]) {
         cout<<"Grid: Sampling "<<sampling<<"_"<<index()<<": Already cell found at pos ("<<ieta<<","<<iphi<<"): id="
             <<m_cell_grid[ieta][iphi]->identify()<<"->"<<refcell->identify()<<endl;
         cout<<"    x="<<m_cells[m_cell_grid[ieta][iphi]->identify()]->x_raw()<<" -> "<<refcell->x_raw();
         cout<<" , y="<<m_cells[m_cell_grid[ieta][iphi]->identify()]->y_raw()<<" -> "<<refcell->y_raw();
         cout<<" mindx="<<m_deta_double<<" mindy="<<m_dphi_double<<endl;
         cout<<endl;
       }      
       m_cell_grid[ieta][iphi]=refcell;
     } else {
       cout<<"Grid: Sampling "<<sampling<<"_"<<index()<<": Cell at pos ("<<ieta<<","<<iphi<<"): id="
           <<refcell->identify()<<" outside eta range"<<endl;
     }
   }
 }

◆ raw_eta_position_to_index()

int CaloGeometryLookup::raw_eta_position_to_index ( float eta_raw ) const

inlineprotected

Definition at line 110 of file CaloGeometryLookup.h.

110 {return TMath::Floor((eta_raw-mineta_raw())/m_deta_double);};

◆ raw_phi_position_to_index()

int CaloGeometryLookup::raw_phi_position_to_index ( float phi_raw ) const

inlineprotected

Definition at line 111 of file CaloGeometryLookup.h.

111 {return TMath::Floor((phi_raw-minphi_raw())/m_dphi_double);};

◆ set_index()

void CaloGeometryLookup::set_index ( int ind )

inline

Definition at line 34 of file CaloGeometryLookup.h.

34 {m_index=ind;};

◆ set_xy_grid_adjustment_factor()

void CaloGeometryLookup::set_xy_grid_adjustment_factor ( float factor )

inline

Definition at line 84 of file CaloGeometryLookup.h.

84 {m_xy_grid_adjustment_factor=factor;};

◆ size()

t_cellmap::size_type CaloGeometryLookup::size ( ) const

inline

Definition at line 32 of file CaloGeometryLookup.h.

32 {return m_cells.size();};

◆ x_correction()

const MeanAndRMS& CaloGeometryLookup::x_correction ( )

inline

Definition at line 79 of file CaloGeometryLookup.h.

79 {return m_eta_correction;};

◆ y_correction()

const MeanAndRMS& CaloGeometryLookup::y_correction ( )

inline

Definition at line 80 of file CaloGeometryLookup.h.

80 {return m_phi_correction;};

Member Data Documentation

◆ m_cell_grid

std::vector< std::vector< const CaloDetDescrElement* > > CaloGeometryLookup::m_cell_grid

protected

Definition at line 117 of file CaloGeometryLookup.h.

◆ m_cell_grid_eta

int CaloGeometryLookup::m_cell_grid_eta

protected

Definition at line 118 of file CaloGeometryLookup.h.

◆ m_cell_grid_phi

int CaloGeometryLookup::m_cell_grid_phi

protected

Definition at line 118 of file CaloGeometryLookup.h.

◆ m_cells

t_cellmap CaloGeometryLookup::m_cells

protected

Definition at line 116 of file CaloGeometryLookup.h.

◆ m_deta

MeanAndRMS CaloGeometryLookup::m_deta

protected

Definition at line 123 of file CaloGeometryLookup.h.

◆ m_deta_double

float CaloGeometryLookup::m_deta_double

protected

Definition at line 125 of file CaloGeometryLookup.h.

◆ m_dphi

MeanAndRMS CaloGeometryLookup::m_dphi

protected

Definition at line 123 of file CaloGeometryLookup.h.

◆ m_dphi_double

float CaloGeometryLookup::m_dphi_double

protected

Definition at line 125 of file CaloGeometryLookup.h.

◆ m_eta_correction

MeanAndRMS CaloGeometryLookup::m_eta_correction

protected

Definition at line 123 of file CaloGeometryLookup.h.

◆ m_index

int CaloGeometryLookup::m_index

protected

Definition at line 115 of file CaloGeometryLookup.h.

◆ m_maxeta

float CaloGeometryLookup::m_maxeta

protected

Definition at line 119 of file CaloGeometryLookup.h.

◆ m_maxeta_correction

float CaloGeometryLookup::m_maxeta_correction

protected

Definition at line 121 of file CaloGeometryLookup.h.

◆ m_maxeta_raw

float CaloGeometryLookup::m_maxeta_raw

protected

Definition at line 120 of file CaloGeometryLookup.h.

◆ m_maxphi

float CaloGeometryLookup::m_maxphi

protected

Definition at line 119 of file CaloGeometryLookup.h.

◆ m_maxphi_correction

float CaloGeometryLookup::m_maxphi_correction

protected

Definition at line 121 of file CaloGeometryLookup.h.

◆ m_maxphi_raw

float CaloGeometryLookup::m_maxphi_raw

protected

Definition at line 120 of file CaloGeometryLookup.h.

◆ m_mindeta

float CaloGeometryLookup::m_mindeta

protected

Definition at line 124 of file CaloGeometryLookup.h.

◆ m_mindphi

float CaloGeometryLookup::m_mindphi

protected

Definition at line 124 of file CaloGeometryLookup.h.

◆ m_mineta

float CaloGeometryLookup::m_mineta

protected

Definition at line 119 of file CaloGeometryLookup.h.

◆ m_mineta_correction

float CaloGeometryLookup::m_mineta_correction

protected

Definition at line 121 of file CaloGeometryLookup.h.

◆ m_mineta_raw

float CaloGeometryLookup::m_mineta_raw

protected

Definition at line 120 of file CaloGeometryLookup.h.

◆ m_minphi

float CaloGeometryLookup::m_minphi

protected

Definition at line 119 of file CaloGeometryLookup.h.

◆ m_minphi_correction

float CaloGeometryLookup::m_minphi_correction

protected

Definition at line 121 of file CaloGeometryLookup.h.

◆ m_minphi_raw

float CaloGeometryLookup::m_minphi_raw

protected

Definition at line 120 of file CaloGeometryLookup.h.

◆ m_phi_correction

MeanAndRMS CaloGeometryLookup::m_phi_correction

protected

Definition at line 123 of file CaloGeometryLookup.h.

◆ m_xy_grid_adjustment_factor

float CaloGeometryLookup::m_xy_grid_adjustment_factor

protected

Definition at line 108 of file CaloGeometryLookup.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CaloGeometryLookup()

◆ ~CaloGeometryLookup()

Member Function Documentation

◆ add()

◆ calculate_distance_eta_phi()

◆ cell_grid_eta()

◆ cell_grid_phi()

◆ deta()

◆ dphi()

◆ dx()

◆ dy()

◆ eta_correction()

◆ getDDE()

◆ has_overlap()

◆ index()

◆ index_range_adjust()

◆ IsCompatible()

◆ maxeta()

◆ maxeta_raw()

◆ maxphi()

◆ maxphi_raw()

◆ maxx()

◆ maxx_raw()

◆ maxy()

◆ maxy_raw()

◆ merge_into_ref()

◆ mindeta()

◆ mindphi()

◆ mindx()

◆ mindy()

◆ mineta()

◆ mineta_raw()

◆ minphi()

◆ minphi_raw()

◆ minx()

◆ minx_raw()

◆ miny()

◆ miny_raw()

◆ neta()

◆ neta_double()

◆ nphi()

◆ nphi_double()

◆ nx()

◆ nx_double()

◆ ny()

◆ ny_double()

◆ phi_correction()

◆ post_process()

◆ raw_eta_position_to_index()

◆ raw_phi_position_to_index()

◆ set_index()

◆ set_xy_grid_adjustment_factor()

◆ size()

◆ x_correction()

◆ y_correction()

Member Data Documentation

◆ m_cell_grid

◆ m_cell_grid_eta

◆ m_cell_grid_phi

◆ m_cells

◆ m_deta

◆ m_deta_double

◆ m_dphi

◆ m_dphi_double

◆ m_eta_correction

◆ m_index

◆ m_maxeta

◆ m_maxeta_correction

◆ m_maxeta_raw

◆ m_maxphi

◆ m_maxphi_correction

◆ m_maxphi_raw

◆ m_mindeta

◆ m_mindphi

◆ m_mineta