RegSelLayer Class Reference

#include <RegSelLayer.h>

Inheritance diagram for RegSelLayer:

Collaboration diagram for RegSelLayer:

Classes
class	PhiCluster

Public Member Functions
	RegSelLayer ()
	RegSelLayer (double rmin, double rmax, double zmin, double zmax)
void	reset ()
void	addModule (const RegSelModule &module)
void	setupDisabledModuleList ()
void	getModules (const RegSelRoI &roi, std::vector< const RegSelModule * > &modules) const
	Fixme: This is the completely correct solution, where each module is allowed to be at a completely different radius.
void	getModules (std::vector< const RegSelModule * > &modules) const
void	getDisabledModules (const RegSelRoI &roi, std::vector< const RegSelModule * > &modules) const
int	ID () const
int	ID (int id)
int	size () const
void	createMaps ()
double	phiMin () const
double	phiMax () const
double	phiMin (double phimin)
double	phiMax (double phimax)
virtual bool	inRoI (const RegSelRoI &roi) const
double	rMin () const
double	rMin (double rmin)
double	rMax () const
double	rMax (double rmax)
double	zMin () const
double	zMin (double zmin)
double	zMax () const
double	zMax (double zmax)
double	z2Min () const
double	z2Min (double zmin)
double	z2Max () const
double	z2Max (double zmax)
bool	enabled () const
bool	enable ()
bool	disable ()
double	_etaMin (double z=0) const
double	_etaMax (double z=0) const

Protected Attributes
double	m_rMin
double	m_rMax
double	m_zMin
double	m_zMax
double	m_z2Min
double	m_z2Max
bool	m_enabled

Private Member Functions
int	clusterModules () const

Private Attributes
bool	m_set
int	m_ID
std::vector< const RegSelModule * >	m_modules
std::vector< const RegSelModule * >	m_disabled
int	m_Nphi
double	m_ideltaphi
double	m_phiMin
double	m_phiMax
std::vector< std::vector< const RegSelModule * > >	m_phimaps

Detailed Description

Definition at line 39 of file RegSelLayer.h.

Constructor & Destructor Documentation

RegSelLayer() [1/2]

RegSelLayer::RegSelLayer ( )

inline

Definition at line 43 of file RegSelLayer.h.

                : 
    ZRObject(0,0,0,0), 
    m_set(false), m_ID(0),
    m_modules(0),   m_disabled(0), m_Nphi(0),
    m_ideltaphi(0), m_phiMin(0), m_phiMax(0), m_phimaps(0) { }   

RegSelLayer() [2/2]

RegSelLayer::RegSelLayer ( double rmin, double rmax, double zmin, double zmax )

inline

Definition at line 49 of file RegSelLayer.h.

                                                                  :
    ZRObject(rmin, rmax, zmin, zmax), 
    m_set(true), m_ID(0), m_modules(0), m_disabled(0), m_Nphi(0),
    m_ideltaphi(0), m_phiMin(0), m_phiMax(0), m_phimaps(0) { }    

Member Function Documentation

_etaMax()

double ZRObject::_etaMax ( double z = 0 ) const

inlineinherited

Definition at line 109 of file ZRObject.h.

                                   { 
    double eta1 = calceta( rMin(), zMax()-z);
    double eta2 = calceta( rMax(), z2Max()-z);
    if ( eta1>eta2 ) return eta1;
    else             return eta2;
  }

_etaMin()

double ZRObject::_etaMin ( double z = 0 ) const

inlineinherited

Definition at line 102 of file ZRObject.h.

                                   { 
    double eta1 = calceta( rMin(), zMin()-z);
    double eta2 = calceta( rMax(), z2Min()-z);
    if ( eta1<eta2 ) return eta1;
    else             return eta2;
  }

addModule()

void RegSelLayer::addModule

(

const RegSelModule &

module

)

Definition at line 130 of file RegSelLayer.cxx.

                                                      { 
  m_modules.push_back(&module);
 
  //  if ( module.enabled() ) m_enabledmodules.push_back(&module);
  //  m_enabledmodules.push_back(&module);
 
  //  std::cout << module 
  //       << "\t rmin=" << rMin() << "\t module.rMin=" << module.rMin() 
  //       << "\t zmin=" << zMin() << "\t module.zMin=" << module.zMin() << std::endl;
 
  if ( m_set ) {
    if ( module.rMin()<rMin()) rMin(module.rMin());
    if ( module.rMax()>rMax()) rMax(module.rMax());
    
    if ( module.zMin()<zMin()) zMin(module.zMin());
    if ( module.zMax()>zMax()) zMax(module.zMax());
 
    if ( module.phiMin()<phiMin()) phiMin(module.phiMin());
    if ( module.phiMax()>phiMax()) phiMax(module.phiMax());
  }
  else { 
    m_set = true;
    rMin(module.rMin());
    rMax(module.rMax());
    zMin(module.zMin());
    zMax(module.zMax());
    phiMin(module.phiMin());
    phiMax(module.phiMax());
  }
 
}

clusterModules()

int RegSelLayer::clusterModules ( ) const

private

Definition at line 50 of file RegSelLayer.cxx.

                                      { 
  
  // if multiple radial "layers" in this actual layer (as in the TRT 
  // barrel, or the SCT endcaps) just find the inner most layer, as 
  // this will have the largest segmentation in phi
 
  if ( m_modules.size()==0 ) return 0;
  
  double r_in  = m_modules.back()->rMax();
 
  for ( int i=m_modules.size()-1 ; i-- ; ) {
    double mod_rmax = m_modules[i]->rMax();
    if ( mod_rmax<r_in )  r_in=mod_rmax;
  }
 
  // now cluster the ends of these innermost modules to calcute how
  // many phi segments to use
 
  // cluster vectors
  PhiCluster cluster;
 
  // use the minimum rMax+4mm just in case the detector has "moved"
  double r2_in = 2*(r_in+4);
 
  for ( int i=m_modules.size() ; i-- ; ) { 
    if ( (m_modules[i]->rMin()+m_modules[i]->rMax())<r2_in )  cluster.addModule(*m_modules[i]);    
  }
 
  return cluster.size();
}

createMaps()

void RegSelLayer::createMaps

(

)

can this really ever happen?

Definition at line 334 of file RegSelLayer.cxx.

                             { 
 
  //  std::cout << "RegSelLayer::createMaps() " << *this << std::endl; 
  
  m_Nphi = clusterModules();
 
  if ( m_Nphi == 0 ) return; 
 
  //  std::cout << "\t\t\tRegSelLayer::createMaps() layer with " << m_Nphi << " phi segments" << std::endl; 
 
  m_phimaps.resize(m_Nphi);
  
  m_phiMin = -M_PI+0.001;
  m_phiMax =  M_PI+0.001;
  
  m_ideltaphi = 0.5*m_Nphi/M_PI;
  
  //  std::map<const RegSelModule*, const RegSelModule*> tmap;
  std::set<const RegSelModule*> tmap;
  
  double trmax = 0;
 
  std::vector<const RegSelModule*>::const_iterator mptr(m_modules.begin());
  for ( ; mptr!=m_modules.end() ; ++mptr ) {
    
    double phicent = 0.5*((*mptr)->phiMax()+(*mptr)->phiMin());
    
    if ( (*mptr)->rMax()>trmax ) trmax = (*mptr)->rMax();
 
    // if cell spans pi boundary
    if ( (*mptr)->phiMax()<(*mptr)->phiMin() ) phicent += M_PI;
    
    if ( phicent>M_PI  ) phicent -= 2*M_PI;
    if ( phicent<-M_PI ) phicent += 2*M_PI;
    
    int phibin = ((int)((phicent-phiMin())*m_ideltaphi))%m_Nphi;
    
    //   std::cout << "phibin=" << (int)phibin 
    //             << "\tphicent=" << phicent 
    //             << "\t  [ " << (*mptr)->phiMin() << " - " << (*mptr)->phiMax() << " ]" << std::endl;
    
    // check module isn't already in another phi segment 
    if ( tmap.find(*mptr)!=tmap.end() ) continue;
    
    //    tmap[*mptr] = *mptr;
    tmap.insert(*mptr);
    m_phimaps[phibin].push_back(*mptr);
    
  }
  
  //  std::cout << "RegSelLayer::createMaps() exiting" << std::endl;   
}

disable()

bool ZRObject::disable ( )

inlineinherited

Definition at line 100 of file ZRObject.h.

{ return m_enabled=false; } 

enable()

bool ZRObject::enable ( )

inlineinherited

Definition at line 99 of file ZRObject.h.

{ return m_enabled=true;  } 

enabled()

bool ZRObject::enabled ( ) const

inlineinherited

Definition at line 98 of file ZRObject.h.

{ return m_enabled; } 

getDisabledModules()

void RegSelLayer::getDisabledModules	(	const RegSelRoI &	roi,
		std::vector< const RegSelModule * > &	modules ) const

Definition at line 314 of file RegSelLayer.cxx.

                                                                                                        { 
 
  //  std::cout << "RegSelLayer::getDisabledModules() " << *this << std::endl;
  //  std::cout << "RegSelLayer::getDisabledModules() m_disabled.size()=" << m_disabled.size() < std::endl;
 
  std::vector<const RegSelModule*>::const_iterator mptr(m_disabled.begin());
  //  std::vector<const RegSelModule*>::const_iterator mptr(m_modules.begin());
 
  for ( ; mptr!=m_disabled.end() ; ++mptr ) { 
    // is the module inside the RoI ???
    if ( (*mptr)->enabled()==false && (*mptr)->inRoI(roi) ) {
      // then add them to the list 
      modules.push_back(*mptr);
    }
  }
 
}

getModules() [1/2]

void RegSelLayer::getModules	(	const RegSelRoI &	roi,
		std::vector< const RegSelModule * > &	modules ) const

Fixme: This is the completely correct solution, where each module is allowed to be at a completely different radius.

If we want to approximate that all the modules are at the same radius within a layer, then we only have to calculate the RoI z end points in this layer once and then just do the comparisons. it will save 2 multiplications and two additions per module, but only for those in the RoI in the phi region, since those outside the RoI in phi aren't even considered at the moment. Tried this, hardly quicker.

NB. Using segmentation in phi and (or) z should speed things up. This implements segmentation in phi only as the z didn't appear much faster.

Definition at line 197 of file RegSelLayer.cxx.

                                                                                                { 
 
  //  std::cout << "RegSelLayer::getModules() " << *this << std::endl;
 
#if 0
 
  // look through all the modules in this layer and 
  // don't use the phi segmentation speed up
 
  std::vector<const RegSelModule*>::const_iterator mptr(m_modules.begin());
 
  for ( ; mptr!=m_modules.end() ; mptr++ ) { 
    // is the module inside the RoI ???
    if ( (*mptr)->enabled() && (*mptr)->inRoI(roi) ) {
      // then add them to the list 
      modules.push_back(*mptr);
    }
  }
 
#else
    
  //  std::cout << "\tRegSelLayer::getModules() " << roi << std::endl;
 
  // get the list of all modules in the phi segments 
  // corresponding to this RoI
    
  int roi_phiMin = (int) ((roi.getphiMin()-phiMin())*m_ideltaphi-1); // look at adjacent segments
  int roi_phiMax = (int) ((roi.getphiMax()-phiMin())*m_ideltaphi+1); // also, hence the +-1
 
  int roi_phiWidth = roi_phiMax-roi_phiMin;
 
  // have to see if the RoI spans the phi=pi boundary, 
  // and also whether we need to check segments that span the booundary
  bool noswap = true;
 
 
  // there MUST be a simpler way to do it than all these add hoc conditions
  // there just has to be, this is *awful* code
  if ( roi_phiWidth>=m_Nphi ) {  
    if ( roi_phiMin>roi_phiMax ) noswap = false;
    roi_phiMin = 0;
    roi_phiMax = m_Nphi-1;
  }
  else { 
    
    if ( roi.getphiWidth()*m_ideltaphi>=m_Nphi-2 ) roi_phiMax = roi_phiMin+m_Nphi-1;
 
    if ( roi_phiMax>=m_Nphi )   { roi_phiMax -= m_Nphi; noswap=false; }
    if ( roi_phiMin<0 )         { roi_phiMin += m_Nphi; noswap=false; }
    if ( roi_phiMin==0 && roi_phiMax<m_Nphi-1 && roi_phiMax!=0 )   { roi_phiMin = m_Nphi-1; noswap=false; }
    //    if ( roi_phiMin==0 && roi_phiMax==0 )   { roi_phiMin = m_Nphi-1; noswap=false; }
 
  }
 
  //    std::cout << "roi  phimin=" << roi.getphiMin() << "\tphimax=" << roi.getphiMax() << std::endl;
  //    std::cout << "roi iphimin=" << roi_phiMin << "\tiphimax=" << roi_phiMax << std::endl;
 
  //  roi_phiWidth = roi_phiMax-roi_phiMin;
 
  // if roi not in phi boundary
  if ( roi.getphiMin()<=roi.getphiMax() && noswap ) { 
    for ( int i=roi_phiMin ; i<=roi_phiMax ; i++ ) {
      std::vector<const RegSelModule*>::const_iterator mptr(m_phimaps[i].begin());  
      std::vector<const RegSelModule*>::const_iterator eptr(m_phimaps[i].end());  
      for ( ; mptr!=eptr ; ++mptr ) {
    if ( (*mptr)->enabled() && (*mptr)->inRoI(roi) ) modules.push_back(*mptr); 
      }
    }
  }
  else { 
    // roi spans phi=pi boundary
    // do phi<0 part of roi
    for ( int i=0 ; i<=roi_phiMax ; i++ ) { 
      std::vector<const RegSelModule*>::const_iterator mptr(m_phimaps[i].begin());  
      std::vector<const RegSelModule*>::const_iterator eptr(m_phimaps[i].end());  
      for ( ; mptr!=eptr ; ++mptr ) {
    if ( (*mptr)->enabled() && (*mptr)->inRoI(roi) )  modules.push_back(*mptr); 
      }
    }
    // do phi>0 part of roi
    // this is needed in case the max and min phi slice are the same
    // so the modules aren't added twice
    int firstbin = (roi_phiMin<roi_phiMax?roi_phiMax+1:roi_phiMin);
    for ( int i=firstbin ; i<m_Nphi ; i++ ) { 
      std::vector<const RegSelModule*>::const_iterator mptr(m_phimaps[i].begin());  
      std::vector<const RegSelModule*>::const_iterator eptr(m_phimaps[i].end());  
      for ( ; mptr!=eptr ; ++mptr ) {
    if ( (*mptr)->enabled() && (*mptr)->inRoI(roi) )  modules.push_back(*mptr); 
      }
    }
  }
  
#endif
 
}

getModules() [2/2]

void RegSelLayer::getModules

(

std::vector< const RegSelModule * > &

modules

)

const

Definition at line 295 of file RegSelLayer.cxx.

                                                                           { 
  std::vector<const RegSelModule*>::const_iterator mptr(m_modules.begin());  
  std::vector<const RegSelModule*>::const_iterator eptr(m_modules.end());  
  for ( ; mptr!=eptr ; ++mptr ) if ( (*mptr)->enabled() ) modules.push_back(*mptr); 
}

ID() [1/2]

int RegSelLayer::ID ( ) const

inline

Definition at line 64 of file RegSelLayer.h.

{ return m_ID; } 

ID() [2/2]

int RegSelLayer::ID ( int id )

inline

Definition at line 65 of file RegSelLayer.h.

{ m_ID = id; return m_ID; } 

inRoI()

virtual bool ZRObject::inRoI ( const RegSelRoI & roi ) const

inlinevirtualinherited

Reimplemented in RegSelModule.

Definition at line 71 of file ZRObject.h.

                                                  {   
    if ( !enabled() ) return false; 
    if ( (roi.getzMax(rMin())<zMin() && roi.getzMax(rMax())<z2Min())  || 
     (roi.getzMin(rMin())>zMax() && roi.getzMin(rMax())>z2Max()) ) {
      return false; 
    }
    return true;
  }

phiMax() [1/2]

double RegSelLayer::phiMax ( ) const

inline

Definition at line 72 of file RegSelLayer.h.

{ return m_phiMax; } 

phiMax() [2/2]

double RegSelLayer::phiMax ( double phimax )

inline

Definition at line 75 of file RegSelLayer.h.

{ return m_phiMax=phimax; } 

phiMin() [1/2]

double RegSelLayer::phiMin ( ) const

inline

Definition at line 71 of file RegSelLayer.h.

{ return m_phiMin; } 

phiMin() [2/2]

double RegSelLayer::phiMin ( double phimin )

inline

Definition at line 74 of file RegSelLayer.h.

{ return m_phiMin=phimin; } 

reset()

void RegSelLayer::reset

(

)

Definition at line 38 of file RegSelLayer.cxx.

                        { 
  m_modules.clear();
  m_disabled.clear();
  std::vector<std::vector<const RegSelModule*> >::iterator itr(m_phimaps.begin()) ; 
  for ( ; itr!=m_phimaps.end() ; ++itr ) (*itr).clear();
  m_set = false;
}

rMax() [1/2]

double ZRObject::rMax ( ) const

inlineinherited

Definition at line 81 of file ZRObject.h.

{ return m_rMax; } 

rMax() [2/2]

double ZRObject::rMax ( double rmax )

inlineinherited

Definition at line 90 of file ZRObject.h.

{ m_rMax = rmax; return m_rMax; } 

rMin() [1/2]

double ZRObject::rMin ( ) const

inlineinherited

Definition at line 80 of file ZRObject.h.

{ return m_rMin; } 

rMin() [2/2]

double ZRObject::rMin ( double rmin )

inlineinherited

Definition at line 89 of file ZRObject.h.

{ m_rMin = rmin; return m_rMin; } 

setupDisabledModuleList()

void RegSelLayer::setupDisabledModuleList

(

)

Definition at line 166 of file RegSelLayer.cxx.

                                          { 
  //  std::cout << "RegSelLayer::setupDisabledModuleList() clear() m_disabled.size()=" << m_disabled.size() << std::endl;
  m_disabled.clear();
  //  std::cout << "RegSelLayer::setupDisabledModuleList() going into loop m_modules.size()=" << m_modules.size() << std::endl;
  std::vector<const RegSelModule*>::const_iterator mptr(m_modules.begin());
  for ( ; mptr!=m_modules.end() ; ++mptr ) { 
    if ( !(*mptr)->enabled() ) m_disabled.push_back(*mptr);
  }
  //  std::cout << "RegSelLayer::setupDisabledModuleList() disabled " << m_disabled.size() << " modules" << std::endl;
}

size()

int RegSelLayer::size ( ) const

inline

Definition at line 67 of file RegSelLayer.h.

{ return m_modules.size(); } 

z2Max() [1/2]

double ZRObject::z2Max ( ) const

inlineinherited

Definition at line 87 of file ZRObject.h.

{ return m_z2Max; } 

z2Max() [2/2]

double ZRObject::z2Max ( double zmax )

inlineinherited

Definition at line 96 of file ZRObject.h.

{ m_z2Max = zmax; return m_z2Max; } 

z2Min() [1/2]

double ZRObject::z2Min ( ) const

inlineinherited

Definition at line 86 of file ZRObject.h.

{ return m_z2Min; } 

z2Min() [2/2]

double ZRObject::z2Min ( double zmin )

inlineinherited

Definition at line 95 of file ZRObject.h.

{ m_z2Min = zmin; return m_z2Min; } 

zMax() [1/2]

double ZRObject::zMax ( ) const

inlineinherited

Definition at line 84 of file ZRObject.h.

{ return m_zMax; } 

zMax() [2/2]

double ZRObject::zMax ( double zmax )

inlineinherited

Definition at line 93 of file ZRObject.h.

{ m_zMax = zmax; return m_zMax; } 

zMin() [1/2]

double ZRObject::zMin ( ) const

inlineinherited

Definition at line 83 of file ZRObject.h.

{ return m_zMin; } 

zMin() [2/2]

double ZRObject::zMin ( double zmin )

inlineinherited

Definition at line 92 of file ZRObject.h.

{ m_zMin = zmin; return m_zMin; } 

Member Data Documentation

m_disabled

std::vector<const RegSelModule*> RegSelLayer::m_disabled

private

Definition at line 103 of file RegSelLayer.h.

m_enabled

bool ZRObject::m_enabled

protectedinherited

Definition at line 122 of file ZRObject.h.

m_ID

int RegSelLayer::m_ID

private

Definition at line 100 of file RegSelLayer.h.

m_ideltaphi

double RegSelLayer::m_ideltaphi

private

Definition at line 110 of file RegSelLayer.h.

m_modules

std::vector<const RegSelModule*> RegSelLayer::m_modules

private

Definition at line 102 of file RegSelLayer.h.

m_Nphi

int RegSelLayer::m_Nphi

private

Definition at line 108 of file RegSelLayer.h.

m_phimaps

std::vector<std::vector<const RegSelModule*> > RegSelLayer::m_phimaps

private

Definition at line 115 of file RegSelLayer.h.

m_phiMax

double RegSelLayer::m_phiMax

private

Definition at line 113 of file RegSelLayer.h.

m_phiMin

double RegSelLayer::m_phiMin

private

Definition at line 112 of file RegSelLayer.h.

m_rMax

double ZRObject::m_rMax

protectedinherited

Definition at line 118 of file ZRObject.h.

m_rMin

double ZRObject::m_rMin

protectedinherited

Definition at line 118 of file ZRObject.h.

m_set

bool RegSelLayer::m_set

private

Definition at line 99 of file RegSelLayer.h.

m_z2Max

double ZRObject::m_z2Max

protectedinherited

Definition at line 120 of file ZRObject.h.

m_z2Min

double ZRObject::m_z2Min

protectedinherited

Definition at line 120 of file ZRObject.h.

m_zMax

double ZRObject::m_zMax

protectedinherited

Definition at line 119 of file ZRObject.h.

m_zMin

double ZRObject::m_zMin

protectedinherited

Definition at line 119 of file ZRObject.h.

---

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

• RegSelLayer.h
• RegSelLayer.cxx