CaloGeometryLookup Class Reference

#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.

{return m_cell_grid_eta;};

cell_grid_phi()

int CaloGeometryLookup::cell_grid_phi ( ) const

inline

Definition at line 83 of file CaloGeometryLookup.h.

{return m_cell_grid_phi;};

deta()

const MeanAndRMS & CaloGeometryLookup::deta ( )

inline

Definition at line 68 of file CaloGeometryLookup.h.

{return m_deta;};

dphi()

const MeanAndRMS & CaloGeometryLookup::dphi ( )

inline

Definition at line 69 of file CaloGeometryLookup.h.

{return m_dphi;};

dx()

const MeanAndRMS & CaloGeometryLookup::dx ( )

inline

Definition at line 72 of file CaloGeometryLookup.h.

{return m_deta;};

dy()

const MeanAndRMS & CaloGeometryLookup::dy ( )

inline

Definition at line 73 of file CaloGeometryLookup.h.

{return m_dphi;};

eta_correction()

const MeanAndRMS & CaloGeometryLookup::eta_correction ( )

inline

Definition at line 77 of file CaloGeometryLookup.h.

{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.

{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.

{return m_maxeta;};

maxeta_raw()

float CaloGeometryLookup::maxeta_raw ( ) const

inline

Definition at line 46 of file CaloGeometryLookup.h.

{return m_maxeta_raw;};

maxphi()

float CaloGeometryLookup::maxphi ( ) const

inline

Definition at line 43 of file CaloGeometryLookup.h.

{return m_maxphi;};

maxphi_raw()

float CaloGeometryLookup::maxphi_raw ( ) const

inline

Definition at line 48 of file CaloGeometryLookup.h.

{return m_maxphi_raw;};

maxx()

float CaloGeometryLookup::maxx ( ) const

inline

Definition at line 59 of file CaloGeometryLookup.h.

{return m_maxeta;};

maxx_raw()

float CaloGeometryLookup::maxx_raw ( ) const

inline

Definition at line 64 of file CaloGeometryLookup.h.

{return m_maxeta_raw;};

maxy()

float CaloGeometryLookup::maxy ( ) const

inline

Definition at line 61 of file CaloGeometryLookup.h.

{return m_maxphi;};

maxy_raw()

float CaloGeometryLookup::maxy_raw ( ) const

inline

Definition at line 66 of file CaloGeometryLookup.h.

{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.

{return m_mindeta;};

mindphi()

float CaloGeometryLookup::mindphi ( )

inline

Definition at line 71 of file CaloGeometryLookup.h.

{return m_mindphi;};

mindx()

float CaloGeometryLookup::mindx ( )

inline

Definition at line 74 of file CaloGeometryLookup.h.

{return m_mindeta;};

mindy()

float CaloGeometryLookup::mindy ( )

inline

Definition at line 75 of file CaloGeometryLookup.h.

{return m_mindphi;};

mineta()

float CaloGeometryLookup::mineta ( ) const

inline

Definition at line 40 of file CaloGeometryLookup.h.

{return m_mineta;};

mineta_raw()

float CaloGeometryLookup::mineta_raw ( ) const

inline

Definition at line 45 of file CaloGeometryLookup.h.

{return m_mineta_raw;};

minphi()

float CaloGeometryLookup::minphi ( ) const

inline

Definition at line 42 of file CaloGeometryLookup.h.

{return m_minphi;};

minphi_raw()

float CaloGeometryLookup::minphi_raw ( ) const

inline

Definition at line 47 of file CaloGeometryLookup.h.

{return m_minphi_raw;};

minx()

float CaloGeometryLookup::minx ( ) const

inline

Definition at line 58 of file CaloGeometryLookup.h.

{return m_mineta;};

minx_raw()

float CaloGeometryLookup::minx_raw ( ) const

inline

Definition at line 63 of file CaloGeometryLookup.h.

{return m_mineta_raw;};

miny()

float CaloGeometryLookup::miny ( ) const

inline

Definition at line 60 of file CaloGeometryLookup.h.

{return m_minphi;};

miny_raw()

float CaloGeometryLookup::miny_raw ( ) const

inline

Definition at line 65 of file CaloGeometryLookup.h.

{return m_minphi_raw;};

neta()

Int_t CaloGeometryLookup::neta ( )

inlineprotected

Definition at line 99 of file CaloGeometryLookup.h.

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

neta_double()

float CaloGeometryLookup::neta_double ( )

inlineprotected

Definition at line 97 of file CaloGeometryLookup.h.

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

nphi()

Int_t CaloGeometryLookup::nphi ( )

inlineprotected

Definition at line 100 of file CaloGeometryLookup.h.

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

nphi_double()

float CaloGeometryLookup::nphi_double ( )

inlineprotected

Definition at line 98 of file CaloGeometryLookup.h.

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

nx()

Int_t CaloGeometryLookup::nx ( )

inlineprotected

Definition at line 105 of file CaloGeometryLookup.h.

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

nx_double()

float CaloGeometryLookup::nx_double ( )

inlineprotected

Definition at line 103 of file CaloGeometryLookup.h.

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

ny()

Int_t CaloGeometryLookup::ny ( )

inlineprotected

Definition at line 106 of file CaloGeometryLookup.h.

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

ny_double()

float CaloGeometryLookup::ny_double ( )

inlineprotected

Definition at line 104 of file CaloGeometryLookup.h.

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

phi_correction()

const MeanAndRMS & CaloGeometryLookup::phi_correction ( )

inline

Definition at line 78 of file CaloGeometryLookup.h.

{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.

{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.

{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.

{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.

{m_xy_grid_adjustment_factor=factor;};

size()

t_cellmap::size_type CaloGeometryLookup::size ( ) const

inline

Definition at line 32 of file CaloGeometryLookup.h.

{return m_cells.size();};

x_correction()

const MeanAndRMS & CaloGeometryLookup::x_correction ( )

inline

Definition at line 79 of file CaloGeometryLookup.h.

{return m_eta_correction;};

y_correction()

const MeanAndRMS & CaloGeometryLookup::y_correction ( )

inline

Definition at line 80 of file CaloGeometryLookup.h.

{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:

• CaloGeometryLookup.h
• CaloGeometryLookup.cxx