TrigFTF_GNNR3_Layer Class Reference

#include <GNNR3_Geometry.h>

Collaboration diagram for TrigFTF_GNNR3_Layer:

Public Member Functions
	TrigFTF_GNNR3_Layer (const TrigInDetSiLayer &, float, int)
	~TrigFTF_GNNR3_Layer ()
int	getEtaBin (float, float) const
float	getMinBinRadius (int) const
float	getMaxBinRadius (int) const
int	num_bins () const
bool	verifyBin (const TrigFTF_GNNR3_Layer *, int, int, float, float) const

Public Attributes
const TrigInDetSiLayer &	m_layer
std::vector< int >	m_bins
std::vector< float >	m_minRadius
std::vector< float >	m_maxRadius
std::vector< float >	m_minBinCoord
std::vector< float >	m_maxBinCoord
float	m_minEta
float	m_maxEta
float	m_etaBin

Protected Attributes
float	m_etaBinWidth
float	m_r1
float	m_z1
float	m_r2
float	m_z2
int	m_nBins

Detailed Description

Definition at line 16 of file GNNR3_Geometry.h.

Constructor & Destructor Documentation

TrigFTF_GNNR3_Layer()

TrigFTF_GNNR3_Layer::TrigFTF_GNNR3_Layer	(	const TrigInDetSiLayer &	ls,
		float	ew,
		int	bin0 )

Definition at line 13 of file GNNR3_Geometry.cxx.

                                                                                       : m_layer(ls), m_etaBinWidth(ew) {
 
  if(m_layer.m_type == 0) {//barrel
    m_r1 = m_layer.m_refCoord;
    m_r2 = m_layer.m_refCoord;
    m_z1 = m_layer.m_minBound;
    m_z2 = m_layer.m_maxBound;
  }
  else {//endcap
    m_r1 = m_layer.m_minBound;
    m_r2 = m_layer.m_maxBound;
    m_z1 = m_layer.m_refCoord;
    m_z2 = m_layer.m_refCoord;
  }
 
  float t1   = m_z1/m_r1;
  float eta1 = -std::log(std::sqrt(1+t1*t1)-t1);
 
 
  float t2   = m_z2/m_r2;
  float eta2 = -std::log(std::sqrt(1+t2*t2)-t2);
 
  m_minEta = eta1;
  m_maxEta = eta2;
  
  if(m_maxEta < m_minEta) {
    m_minEta = eta2;
    m_maxEta = eta1;
  }
 
  m_maxEta += 1e-6;//increasing them slightly to avoid range_check exceptions
  m_minEta -= 1e-6;
 
  float deltaEta = m_maxEta - m_minEta;
 
  int binCounter = bin0;
 
  if(deltaEta < m_etaBinWidth) {
    m_nBins = 1;
    m_bins.push_back(binCounter++);
    m_etaBin = deltaEta;
    if(m_layer.m_type == 0) {//barrel
      m_minRadius.push_back(m_layer.m_refCoord - 2.0);
      m_maxRadius.push_back(m_layer.m_refCoord + 2.0);
      m_minBinCoord.push_back(m_layer.m_minBound);
      m_maxBinCoord.push_back(m_layer.m_maxBound);
    }
    else {//endcap
      m_minRadius.push_back(m_layer.m_minBound - 2.0);
      m_maxRadius.push_back(m_layer.m_maxBound + 2.0);
      m_minBinCoord.push_back(m_layer.m_minBound);
      m_maxBinCoord.push_back(m_layer.m_maxBound);
    }
  }
  else {
    float nB = static_cast<int>(deltaEta/m_etaBinWidth);
    m_nBins = nB;
    if(deltaEta - m_etaBinWidth*nB > 0.5*m_etaBinWidth) m_nBins++;
 
    m_etaBin = deltaEta/m_nBins;
 
    if(m_nBins == 1) {
      m_bins.push_back(binCounter++);
      if(m_layer.m_type == 0) {//barrel
    m_minRadius.push_back(m_layer.m_refCoord - 2.0);
        m_maxRadius.push_back(m_layer.m_refCoord + 2.0);
    m_minBinCoord.push_back(m_layer.m_minBound);
    m_maxBinCoord.push_back(m_layer.m_maxBound);
      }
      else {//endcap
    m_minRadius.push_back(m_layer.m_minBound - 2.0);
        m_maxRadius.push_back(m_layer.m_maxBound + 2.0);
    m_minBinCoord.push_back(m_layer.m_minBound);
    m_maxBinCoord.push_back(m_layer.m_maxBound);
      }
    }
    else {
 
      float eta = m_minEta+0.5*m_etaBin;
 
      for(int i=1;i<=m_nBins;i++) {
 
    m_bins.push_back(binCounter++);
    
    float e1 = eta - 0.5*m_etaBin;
    float e2 = eta + 0.5*m_etaBin;
    
    if(m_layer.m_type == 0) {//barrel
      m_minRadius.push_back(m_layer.m_refCoord - 2.0);
          m_maxRadius.push_back(m_layer.m_refCoord + 2.0);
      float z1 = m_layer.m_refCoord*std::sinh(e1);
      m_minBinCoord.push_back(z1);
      float z2 = m_layer.m_refCoord*std::sinh(e2);
      m_maxBinCoord.push_back(z2);
    }
    else {//endcap
      float r = m_layer.m_refCoord/std::sinh(e1);
      m_minBinCoord.push_back(r);
      m_minRadius.push_back(r - 2.0);
      r = m_layer.m_refCoord/std::sinh(e2);
      m_maxBinCoord.push_back(r);
      m_maxRadius.push_back(r + 2.0);
      
    }
    
    eta += m_etaBin;
      }
    }
  }
}

~TrigFTF_GNNR3_Layer()

TrigFTF_GNNR3_Layer::~TrigFTF_GNNR3_Layer

(

)

Definition at line 215 of file GNNR3_Geometry.cxx.

                                          {
  m_bins.clear();
}

Member Function Documentation

getEtaBin()

int TrigFTF_GNNR3_Layer::getEtaBin	(	float	zh,
		float	rh ) const

Definition at line 184 of file GNNR3_Geometry.cxx.

                                                           {
  
  if(m_bins.size() == 1) return m_bins.at(0);
 
  float t1   = zh/rh;
  float eta = -std::log(std::sqrt(1+t1*t1)-t1);
  
 
  int idx = static_cast<int>((eta-m_minEta)/m_etaBin);
 
  
  if(idx < 0) idx = 0;
  if(idx >= static_cast<int>(m_bins.size())) idx = static_cast<int>(m_bins.size())-1;
 
  return m_bins.at(idx);//index in the global storage
}

getMaxBinRadius()

float TrigFTF_GNNR3_Layer::getMaxBinRadius

(

int

idx

)

const

Definition at line 208 of file GNNR3_Geometry.cxx.

                                                        {
  if(idx >= static_cast<int>(m_maxRadius.size())) idx = idx-1;
  if(idx < 0) idx = 0;
  
  return m_maxRadius.at(idx);
}

getMinBinRadius()

float TrigFTF_GNNR3_Layer::getMinBinRadius

(

int

idx

)

const

Definition at line 201 of file GNNR3_Geometry.cxx.

                                                        {
  if(idx >= static_cast<int>(m_minRadius.size())) idx = idx-1;
  if(idx < 0) idx = 0;
  
  return m_minRadius.at(idx);
}

num_bins()

int TrigFTF_GNNR3_Layer::num_bins ( ) const

inline

Definition at line 26 of file GNNR3_Geometry.h.

{return m_bins.size();} 

verifyBin()

bool TrigFTF_GNNR3_Layer::verifyBin	(	const TrigFTF_GNNR3_Layer *	pL,
		int	b1,
		int	b2,
		float	min_z0,
		float	max_z0 ) const

Definition at line 124 of file GNNR3_Geometry.cxx.

                                                                                                                   {
 
  float z1min = m_minBinCoord.at(b1);
  float z1max = m_maxBinCoord.at(b1);
  float r1 = m_layer.m_refCoord;
 
  if(m_layer.m_type == 0 && pL->m_layer.m_type == 0) {//barrel -> barrel
 
    const float tol = 5.0;
 
    float min_b2 = pL->m_minBinCoord.at(b2);
    float max_b2 = pL->m_maxBinCoord.at(b2);
 
    float r2 = pL->m_layer.m_refCoord;
 
    float A = r2/(r2-r1);
    float B = r1/(r2-r1);
 
    float z0_min = z1min*A - max_b2*B;
    float z0_max = z1max*A - min_b2*B;
    
    if(z0_max < min_z0-tol || z0_min > max_z0+tol) return false;
 
    return true;
  }
 
  if(m_layer.m_type == 0 && pL->m_layer.m_type != 0) {//barrel -> endcap
 
    const float tol = 10.0;
 
    float z2 = pL->m_layer.m_refCoord;
    float r2max = pL->m_maxBinCoord.at(b2);
    float r2min = pL->m_minBinCoord.at(b2);
 
    if(r2max <= r1) return false;
 
    if(r2min <= r1) {
      r2min = r1 + 1e-3;
    }
 
    float z0_max = 0.0;
    float z0_min = 0.0;
    
    if(z2 > 0) { 
      z0_max = (z1max*r2max - z2*r1)/(r2max-r1);
      z0_min = (z1min*r2min - z2*r1)/(r2min-r1);
    }
    else {
      z0_max = (z1max*r2min - z2*r1)/(r2min-r1);
      z0_min = (z1min*r2max - z2*r1)/(r2max-r1);
    }
 
    if(z0_max < min_z0-tol || z0_min > max_z0+tol) return false;
    return true;
  }
 
   return true;
}

Member Data Documentation

m_bins

std::vector<int> TrigFTF_GNNR3_Layer::m_bins

Definition at line 31 of file GNNR3_Geometry.h.

m_etaBin

float TrigFTF_GNNR3_Layer::m_etaBin

Definition at line 37 of file GNNR3_Geometry.h.

m_etaBinWidth

float TrigFTF_GNNR3_Layer::m_etaBinWidth

protected

Definition at line 41 of file GNNR3_Geometry.h.

m_layer

const TrigInDetSiLayer& TrigFTF_GNNR3_Layer::m_layer

Definition at line 30 of file GNNR3_Geometry.h.

m_maxBinCoord

std::vector<float> TrigFTF_GNNR3_Layer::m_maxBinCoord

Definition at line 35 of file GNNR3_Geometry.h.

m_maxEta

float TrigFTF_GNNR3_Layer::m_maxEta

Definition at line 37 of file GNNR3_Geometry.h.

m_maxRadius

std::vector<float> TrigFTF_GNNR3_Layer::m_maxRadius

Definition at line 33 of file GNNR3_Geometry.h.

m_minBinCoord

std::vector<float> TrigFTF_GNNR3_Layer::m_minBinCoord

Definition at line 34 of file GNNR3_Geometry.h.

m_minEta

float TrigFTF_GNNR3_Layer::m_minEta

Definition at line 37 of file GNNR3_Geometry.h.

m_minRadius

std::vector<float> TrigFTF_GNNR3_Layer::m_minRadius

Definition at line 32 of file GNNR3_Geometry.h.

m_nBins

int TrigFTF_GNNR3_Layer::m_nBins

protected

Definition at line 44 of file GNNR3_Geometry.h.

m_r1

float TrigFTF_GNNR3_Layer::m_r1

protected

Definition at line 43 of file GNNR3_Geometry.h.

m_r2

float TrigFTF_GNNR3_Layer::m_r2

protected

Definition at line 43 of file GNNR3_Geometry.h.

m_z1

float TrigFTF_GNNR3_Layer::m_z1

protected

Definition at line 43 of file GNNR3_Geometry.h.

m_z2

float TrigFTF_GNNR3_Layer::m_z2

protected

Definition at line 43 of file GNNR3_Geometry.h.

---

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

• GNNR3_Geometry.h
• GNNR3_Geometry.cxx