InDet::TRT_DetElementsLayer_xk Class Reference

#include <TRT_DetElementsLayer_xk.h>

Collaboration diagram for InDet::TRT_DetElementsLayer_xk:

Public Member Functions
	TRT_DetElementsLayer_xk ()
	TRT_DetElementsLayer_xk (double, double, double, double, double)
	TRT_DetElementsLayer_xk (const TRT_DetElementsLayer_xk &)=default
	~TRT_DetElementsLayer_xk ()
TRT_DetElementsLayer_xk &	operator= (const TRT_DetElementsLayer_xk &)=default
TRT_DetElementsLayer_xk &	operator= (TRT_DetElementsLayer_xk &&)=default
const float &	r () const
const float &	dr () const
const float &	z () const
const float &	dz () const
const float &	dfe () const
void	set (double, double, double, double, double, double, double)
void	add (TRT_DetElementLink_xk &&)
int	nElements () const
	Get number of links to detector elements.
void	getBarrelDetElementsATL (const float *, const float *, std::vector< std::pair< const InDet::TRT_DetElementLink_xk *, float > > &, std::vector< InDet::TRT_DetElementLink_xk::Used_t > &used) const
void	getBarrelDetElementsCTB (const float *, const float *, std::vector< std::pair< const InDet::TRT_DetElementLink_xk *, float > > &, std::vector< InDet::TRT_DetElementLink_xk::Used_t > &used) const
void	getEndcapDetElements (const float *, const float *, std::vector< std::pair< const InDet::TRT_DetElementLink_xk *, float > > &, std::vector< InDet::TRT_DetElementLink_xk::Used_t > &used) const
void	reserve (std::size_t n_elements)

Protected Attributes
float	m_z
float	m_dz
float	m_r
float	m_dr
float	m_dfe
float	m_f0
float	m_sfi
float	m_wf
float	m_wz
std::vector< TRT_DetElementLink_xk >	m_elements

Detailed Description

Definition at line 27 of file TRT_DetElementsLayer_xk.h.

Constructor & Destructor Documentation

TRT_DetElementsLayer_xk() [1/3]

InDet::TRT_DetElementsLayer_xk::TRT_DetElementsLayer_xk ( )

inline

Definition at line 125 of file TRT_DetElementsLayer_xk.h.

  {
    m_z   = 0.;
    m_dz  = 0.;
    m_r   = 0.;
    m_dr  = 0.;
    m_dfe = 0.;
    m_f0  = 0.;
    m_sfi = 0.;
    m_wf  = 0.;
    m_wz  = 0.;
  } 

TRT_DetElementsLayer_xk() [2/3]

InDet::TRT_DetElementsLayer_xk::TRT_DetElementsLayer_xk ( double r, double dr, double z, double dz, double df )

inline

Definition at line 138 of file TRT_DetElementsLayer_xk.h.

    {
      m_r   = float(r );
      m_dr  = float(dr);
      m_z   = float(z );
      m_dz  = float(dz);
      m_dfe = float(df);
      m_f0  = 0.       ;
      m_sfi = 0.       ;
      m_wf  = 0.       ;
      m_wz  = 0.       ;
    } 

TRT_DetElementsLayer_xk() [3/3]

InDet::TRT_DetElementsLayer_xk::TRT_DetElementsLayer_xk ( const TRT_DetElementsLayer_xk & )

default

~TRT_DetElementsLayer_xk()

InDet::TRT_DetElementsLayer_xk::~TRT_DetElementsLayer_xk ( )

inline

Definition at line 152 of file TRT_DetElementsLayer_xk.h.

{}

Member Function Documentation

add()

void InDet::TRT_DetElementsLayer_xk::add ( TRT_DetElementLink_xk && link )

inline

Definition at line 154 of file TRT_DetElementsLayer_xk.h.

    {
       m_elements.push_back(std::move(link));
    }

dfe()

const float & InDet::TRT_DetElementsLayer_xk::dfe ( ) const

inline

Definition at line 54 of file TRT_DetElementsLayer_xk.h.

{return this->m_dfe;}

dr()

const float & InDet::TRT_DetElementsLayer_xk::dr ( ) const

inline

Definition at line 51 of file TRT_DetElementsLayer_xk.h.

{return this->m_dr ;}

dz()

const float & InDet::TRT_DetElementsLayer_xk::dz ( ) const

inline

Definition at line 53 of file TRT_DetElementsLayer_xk.h.

{return this->m_dz ;}

getBarrelDetElementsATL()

void InDet::TRT_DetElementsLayer_xk::getBarrelDetElementsATL	(	const float *	P,
		const float *	A,
		std::vector< std::pair< const InDet::TRT_DetElementLink_xk *, float > > &	lDE,
		std::vector< InDet::TRT_DetElementLink_xk::Used_t > &	used ) const

Parameters

?
?
will	be filled with pairs of detector elements and way(?)
array	with size matching number of elements, which tells which elements are already used.

Definition at line 55 of file TRT_DetElementsLayer_xk.cxx.

{
  // Tell clang to optimize assuming that FP exceptions can trap.
  // Otherwise, it can vectorize the division, which can lead to
  // spurious division-by-zero traps from unused vector lanes.
  CXXUTILS_TRAPPING_FP;
  float a    = (A[0]*P[0]+A[1]*P[1])*2.; 
  float s    = 0. ;
  if(a!=0.) {
    float d  = (m_r-P[0]-P[1])*(m_r+P[0]+P[1])+2.*P[0]*P[1];
    float b  = 2.*(A[0]*A[0]+A[1]*A[1]);
    float sq = a*a+2.*d*b;  sq>0. ? sq=std::sqrt(sq) : sq=0.;
    float s1 =-(a+sq)/b;
    float s2 =-(a-sq)/b; 
 
    if((s1*s2) > 0.) {std::abs(s1) < std::abs(s2) ? s = s1 : s = s2;}
    else             {     s1  > 0.       ? s = s1 : s = s2;}  
 
  }
  float zc  = P[2]+A[2]*s ;  if(std::abs(zc-m_z)>(m_dz+P[4])) return;
  float zc0 = zc-P[4]-m_z ;
  float zc1 = zc+P[4]-m_z ;
  float yc  = P[1]+A[1]*s ;
  float xc  = P[0]+A[0]*s ; s+=P[5];
  float fc  = std::atan2(yc,xc);
  int m     = m_elements.size()-2;
  int n     = int((fc-m_f0)*m_sfi)*2; if(n<0) n=0; else if(n>m) n=m;
 
  float dF = fc-m_elements[n].phi();
  float dx = std::abs(yc*m_elements[n].cos()-xc*m_elements[n].sin()-m_elements[n].centerf());
  if((dx-m_wf-P[4]) <= 0.) {
 
    assert( used.size() > static_cast<unsigned int>(n));
    if(zc0 <= 0 && !used[n].used()) {
       lDE.emplace_back(&m_elements[n],s); used[n].setUsed();
    }
 
    int k = n+1;
    assert( used.size() > static_cast<unsigned int>(k));
    if(zc1 >= 0. && !used[k].used() ) {
       lDE.emplace_back(&m_elements[k],s); used[k].setUsed();
    }
  }
  if(dF>0.) { n!=62 ? n+=2 : n=0 ;}
  else      { n!= 0 ? n-=2 : n=62;}
 
  dx = std::abs(yc*m_elements[n].cos()-xc*m_elements[n].sin()-m_elements[n].centerf());
  if((dx-m_wf-P[4]) <= 0.) {
 
    assert( used.size() > static_cast<unsigned int>(n));
    if(zc0 <= 0. && !used[n].used()) {
       lDE.emplace_back(&m_elements[n],s); used[n].setUsed();
    }
 
    int k = n+1;
    assert( used.size() > static_cast<unsigned int>(k));
    if(zc1 >= 0. && !used[k].used()) {
       lDE.emplace_back(&m_elements[k],s); used[k].setUsed();
    }
  }
}

getBarrelDetElementsCTB()

void InDet::TRT_DetElementsLayer_xk::getBarrelDetElementsCTB	(	const float *	P,
		const float *	A,
		std::vector< std::pair< const InDet::TRT_DetElementLink_xk *, float > > &	lDE,
		std::vector< InDet::TRT_DetElementLink_xk::Used_t > &	used ) const

Parameters

?
?
will	be filled with pairs of detector elements and way(?)
array	with size matching number of elements, which tells which elements are already used.

Definition at line 131 of file TRT_DetElementsLayer_xk.cxx.

{
  float a    = (A[0]*P[0]+A[1]*P[1])*2.; 
  float s    = 0. ;
  if(a!=0.) {
    float d  = (m_r-P[0]-P[1])*(m_r+P[0]+P[1])+2.*P[0]*P[1];
    float b  = 2.*(A[0]*A[0]+A[1]*A[1]);
    float sq = a*a+2.*d*b;  sq>0. ? sq=std::sqrt(sq) : sq=0.;
    float s1 =-(a+sq)/b;
    float s2 =-(a-sq)/b; 
 
    if((s1*s2) > 0.) {std::abs(s1) < std::abs(s2) ? s = s1 : s = s2;}
    else             {     s1  > 0.       ? s = s1 : s = s2;}  
 
  }
  float zc  = P[2]+A[2]*s ;  if(std::abs(zc-m_z)>(m_dz+P[4])) return;
  float zc0 = zc-P[4]-m_z ;
  float zc1 = zc+P[4]-m_z ;
  float yc  = P[1]+A[1]*s ;
  float xc  = P[0]+A[0]*s ; s+=P[5];
  
  for(size_t n = 0; n < m_elements.size(); n+=2) {
 
    float dx = std::abs(yc*m_elements[n].cos()-xc*m_elements[n].sin()-m_elements[n].centerf());
    if((dx-m_wf-P[4]) <= 0.) {
 
      assert( used.size() > static_cast<unsigned int>(n));
      if(zc0 <= 0. &&  !used[n].used()) {
         lDE.emplace_back(&m_elements[n],s); used[n].setUsed();
      }
 
      assert( used.size() > static_cast<unsigned int>(n+1));
      if(zc1 >= 0. && !used[n+1].used()) {
 
         lDE.emplace_back(&m_elements[n+1],s); used[n+1].setUsed();
      }
    }
  }
}

getEndcapDetElements()

void InDet::TRT_DetElementsLayer_xk::getEndcapDetElements	(	const float *	P,
		const float *	A,
		std::vector< std::pair< const InDet::TRT_DetElementLink_xk *, float > > &	lDE,
		std::vector< InDet::TRT_DetElementLink_xk::Used_t > &	used ) const

Parameters

?
?
will	be filled with pairs of detector elements and way(?)
array	with size matching number of elements, which tells which elements are already used.

Definition at line 185 of file TRT_DetElementsLayer_xk.cxx.

{
  const float pi2=2.*M_PI;
 
  float s   = (m_z-P[2])/A[2];
  float xc  = P[0]+A[0]*s;
  float yc  = P[1]+A[1]*s; s+=P[5]; 
  float rc  = std::sqrt(xc*xc+yc*yc); if(std::abs(rc-m_r)>m_dr+P[4]) return;
  float fc  = std::atan2(yc,xc);
  float sf  = 0.09817477+P[4]/rc;
 
  int m     = m_elements.size()-1;
  int n     = int((fc-m_f0)*m_sfi); if(n<0) n=0; else if(n>m) n=m;
  float dF  = fc-m_elements[n].phi();
 
  assert( used.size() > static_cast<unsigned int>(n));
  if(std::abs(dF)                     < sf && !used[n].used()) {
     lDE.emplace_back(&m_elements[n],s); used[n].setUsed();
  }
 
  //  if(dF>0.) {if(n!=63) ++n; else {n=0 ; fc-=pi2;}}
  //  else      {if(n!=0 ) --n; else {n=63; fc+=pi2;}}
  if(dF>0.) {if(n!=m) ++n; else {n=0 ; fc-=pi2;}}
  else      {if(n!=0) --n; else {n=m ; fc+=pi2;}}
 
  assert( used.size() > static_cast<size_t>(n));
  if(std::abs(fc-m_elements[n].phi()) < sf && !used[n].used()) {
     lDE.emplace_back(&m_elements[n],s); used[n].setUsed();
  }
}

nElements()

int InDet::TRT_DetElementsLayer_xk::nElements ( ) const

inline

Get number of links to detector elements.

Definition at line 159 of file TRT_DetElementsLayer_xk.h.

    {
      return m_elements.size();
    }  

operator=() [1/2]

TRT_DetElementsLayer_xk & InDet::TRT_DetElementsLayer_xk::operator= ( const TRT_DetElementsLayer_xk & )

default

operator=() [2/2]

TRT_DetElementsLayer_xk & InDet::TRT_DetElementsLayer_xk::operator= ( TRT_DetElementsLayer_xk && )

default

r()

const float & InDet::TRT_DetElementsLayer_xk::r ( ) const

inline

Definition at line 50 of file TRT_DetElementsLayer_xk.h.

{return this->m_r  ;}

reserve()

void InDet::TRT_DetElementsLayer_xk::reserve ( std::size_t n_elements )

inline

Definition at line 96 of file TRT_DetElementsLayer_xk.h.

                                         {
         m_elements.reserve(n_elements);
      }

set()

void InDet::TRT_DetElementsLayer_xk::set	(	double	r,
		double	dr,
		double	z,
		double	dz,
		double	df,
		double	wf,
		double	wz )

Definition at line 27 of file TRT_DetElementsLayer_xk.cxx.

{
  m_r   = float(r ) ;
  m_dr  = float(dr) ;
  m_z   = float(z ) ;
  m_dz  = float(dz) ;
  m_dfe = float(df) ;
  m_wf  = float(wf) ;
  m_wz  = float(wz) ;
  m_f0  = m_elements[0].phi();
  m_sfi = 0.        ;
  if(m_elements.size()<=1) return;
  m_sfi = m_elements[1].phi()-m_elements[0].phi();
  if(m_sfi < .001) m_sfi = m_elements[2].phi()-m_elements[0].phi();
  m_sfi = 1./m_sfi;
} 

z()

const float & InDet::TRT_DetElementsLayer_xk::z ( ) const

inline

Definition at line 52 of file TRT_DetElementsLayer_xk.h.

{return this->m_z  ;}

Member Data Documentation

m_dfe

float InDet::TRT_DetElementsLayer_xk::m_dfe

protected

Definition at line 109 of file TRT_DetElementsLayer_xk.h.

m_dr

float InDet::TRT_DetElementsLayer_xk::m_dr

protected

Definition at line 108 of file TRT_DetElementsLayer_xk.h.

m_dz

float InDet::TRT_DetElementsLayer_xk::m_dz

protected

Definition at line 106 of file TRT_DetElementsLayer_xk.h.

m_elements

std::vector<TRT_DetElementLink_xk> InDet::TRT_DetElementsLayer_xk::m_elements

protected

Definition at line 114 of file TRT_DetElementsLayer_xk.h.

m_f0

float InDet::TRT_DetElementsLayer_xk::m_f0

protected

Definition at line 110 of file TRT_DetElementsLayer_xk.h.

m_r

float InDet::TRT_DetElementsLayer_xk::m_r

protected

Definition at line 107 of file TRT_DetElementsLayer_xk.h.

m_sfi

float InDet::TRT_DetElementsLayer_xk::m_sfi

protected

Definition at line 111 of file TRT_DetElementsLayer_xk.h.

m_wf

float InDet::TRT_DetElementsLayer_xk::m_wf

protected

Definition at line 112 of file TRT_DetElementsLayer_xk.h.

m_wz

float InDet::TRT_DetElementsLayer_xk::m_wz

protected

Definition at line 113 of file TRT_DetElementsLayer_xk.h.

m_z

float InDet::TRT_DetElementsLayer_xk::m_z

protected

Definition at line 105 of file TRT_DetElementsLayer_xk.h.

---

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

• TRT_DetElementsLayer_xk.h
• TRT_DetElementsLayer_xk.cxx