T2GeometryTool Class Reference

#include <T2GeometryTool.h>

Inheritance diagram for T2GeometryTool:

Collaboration diagram for T2GeometryTool:

Public Member Functions
	T2GeometryTool (const std::string &type, const std::string &name, const IInterface *parent)
	Constructor.
virtual StatusCode	initialize () override
	Initialize method.
virtual int	EtaPhiRange (const int nCaloType, const int nCaloSamp, const double &energyEta, const double &energyPhi) override
	This prepares internal tables based on the cluster position.
virtual bool	CellInNorCluster (const int nCaloType, const int nCaloSamp, const double &etaCell, const double &phiCell) override
	Finds whether a cell is in a Normal Sized Cluster.
virtual bool	CellInNarCluster (const int nCaloType, const int nCaloSamp, const double &etaCell, const double &phiCell) override
	Finds whether a cell is in a Narrow Sized Cluster.
virtual bool	CellInWidCluster (const int nCaloType, const int nCaloSamp, const double &etaCell, const double &phiCell) override
	Finds whether a cell is in a Wide Sized Cluster.
virtual int	nStripetaEMnar (const int layer) const override
	Get number of cells in eta for Narrow Sized Cluster in EM calorimeter.
virtual int	nStripphiEMnar (const int layer) const override
	Get number of cells in phi for Narrow Sized Cluster in EM calorimeter.
virtual int	nStripetaEMwid (const int layer) const override
	Get number of cells in eta for Wide Sized Cluster in EM calorimeter.
virtual int	nStripphiEMwid (const int layer) const override
	Get number of cells in phi for Wide Sized Cluster in EM calorimeter.
virtual int	nStripetaEMnor (const int layer) const override
	Get number of cells in eta for Normal Sized Cluster in EM calorimeter.
virtual int	nStripphiEMnor (const int layer) const override
	Get number of cells in phi for Normal Sized Cluster in EM calorimeter.
virtual int	nStripetaHADnar (const int layer) const override
	Get number of cells in eta for Narrow Sized Cluster in HAD calorimeter.
virtual int	nStripphiHADnar (const int layer) const override
	Get number of cells in phi for Narrow Sized Cluster in HAD calorimeter.
virtual int	nStripetaHADwid (const int layer) const override
	Get number of cells in eta for Wide Sized Cluster in HADcalorimeter.
virtual int	nStripphiHADwid (const int layer) const override
	Get number of cells in phi for Wide Sized Cluster in HAD calorimeter.
virtual int	nStripetaHADnor (const int layer) const override
	Get number of cells in eta for Normal Sized Cluster in HAD calorimeter.
virtual int	nStripphiHADnor (const int layer) const override
	Get number of cells in phi for Normal Sized Cluster in HAD calorimeter.
	DeclareInterfaceID (IT2GeometryTool, 1, 0)
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
virtual int	EtaPhiRangeInt (const int layer, const int &nStripEta, const int &nStripPhi, const double &energyEta, const double &energyPhi, const int index_size) override
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
double	m_CellEtaNorMin [7] {}
	Normal cluster variables.
double	m_CellEtaNorMax [7] {}
double	m_CellPhiNorMin [7] {}
double	m_CellPhiNorMax [7] {}
double	m_CellEtaNarMin [7] {}
	Narrow cluster variables.
double	m_CellEtaNarMax [7] {}
double	m_CellPhiNarMin [7] {}
double	m_CellPhiNarMax [7] {}
double	m_CellEtaWidMin [7] {}
	Wider cluster variables.
double	m_CellEtaWidMax [7] {}
double	m_CellPhiWidMin [7] {}
double	m_CellPhiWidMax [7] {}
int	m_nStripetaEMnar [4] {}
	Reference changed to vectors.
int	m_nStripphiEMnar [4] {}
int	m_nStripetaEMwid [4] {}
int	m_nStripphiEMwid [4] {}
int	m_nStripetaEMnor [4] {}
int	m_nStripphiEMnor [4] {}
int	m_nStripetaHADnar [3] {}
int	m_nStripphiHADnar [3] {}
int	m_nStripetaHADwid [3] {}
int	m_nStripphiHADwid [3] {}
int	m_nStripetaHADnor [3] {}
int	m_nStripphiHADnor [3] {}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Static Private Attributes
static const int	m_netareg [7] = {3,7,4,3,2,2,3}
	Granularity and eta ranges relevant for above energy sums.
static const double	m_etareg [7][8]
static const double	m_etaGran [7][7]
static const int	m_nphireg [7] = {1,3,2,1,2,2,2}
static const double	m_phireg [7][4]
static const double	m_phiGran [7][3]

Detailed Description

Definition at line 25 of file T2GeometryTool.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

T2GeometryTool()

T2GeometryTool::T2GeometryTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Constructor.

Loads cluster size definition tables.

Definition at line 9 of file T2GeometryTool.cxx.

                                              : AthAlgTool(type,name,parent)
{
 
  declareProperty( "nStripEtaEM0nar" , m_nStripetaEMnar[0] = 4); // 0.1
  declareProperty( "nStripPhiEM0nar" , m_nStripphiEMnar[0] = 1); // 0.098
  declareProperty( "nStripEtaEM1nar" , m_nStripetaEMnar[1] = 32); // 0.1
  declareProperty( "nStripPhiEM1nar" , m_nStripphiEMnar[1] = 1); // 0.098
  declareProperty( "nStripEtaEM2nar" , m_nStripetaEMnar[2] = 4); // 0.1
  declareProperty( "nStripPhiEM2nar" , m_nStripphiEMnar[2] = 4); // 0.098
  declareProperty( "nStripEtaEM3nar" , m_nStripetaEMnar[3] = 2); // 0.1
  declareProperty( "nStripPhiEM3nar" , m_nStripphiEMnar[3] = 4); // 0.098
 
  declareProperty( "nStripEtaEM0wid" , m_nStripetaEMwid[0] = 8); // 0.2
  declareProperty( "nStripPhiEM0wid" , m_nStripphiEMwid[0] = 2); // 0.196
  declareProperty( "nStripEtaEM1wid" , m_nStripetaEMwid[1] = 64); // 0.2
  declareProperty( "nStripPhiEM1wid" , m_nStripphiEMwid[1] = 2); // 0.196
  declareProperty( "nStripEtaEM2wid" , m_nStripetaEMwid[2] = 8); // 0.2
  declareProperty( "nStripPhiEM2wid" , m_nStripphiEMwid[2] = 8); // 0.196
  declareProperty( "nStripEtaEM3wid" , m_nStripetaEMwid[3] = 4); // 0.2
  declareProperty( "nStripPhiEM3wid" , m_nStripphiEMwid[3] = 8); // 0.196
 
  declareProperty( "nStripEtaEM0nor" , m_nStripetaEMnor[0] = 16); // 0.4
  declareProperty( "nStripPhiEM0nor" , m_nStripphiEMnor[0] = 4); // 0.39
  declareProperty( "nStripEtaEM1nor" , m_nStripetaEMnor[1] = 128); // 0.4
  declareProperty( "nStripPhiEM1nor" , m_nStripphiEMnor[1] = 4); // 0.39
  declareProperty( "nStripEtaEM2nor" , m_nStripetaEMnor[2] = 16); // 0.4
  declareProperty( "nStripPhiEM2nor" , m_nStripphiEMnor[2] = 16); // 0.39
  declareProperty( "nStripEtaEM3nor" , m_nStripetaEMnor[3] = 8); // 0.4
  declareProperty( "nStripPhiEM3nor" , m_nStripphiEMnor[3] = 16); // 0.39
 
 
  declareProperty( "nStripEtaHAD0nar" , m_nStripetaHADnar[0] = 1); // 0.1
  declareProperty( "nStripPhiHAD0nar" , m_nStripphiHADnar[0] = 1); // 0.098
  declareProperty( "nStripEtaHAD1nar" , m_nStripetaHADnar[1] = 1); // 0.1
  declareProperty( "nStripPhiHAD1nar" , m_nStripphiHADnar[1] = 1); // 0.098
  declareProperty( "nStripEtaHAD2nar" , m_nStripetaHADnar[2] = 1); // 0.2
  declareProperty( "nStripPhiHAD2nar" , m_nStripphiHADnar[2] = 1); // 0.098
 
  declareProperty( "nStripEtaHAD0wid" , m_nStripetaHADwid[0] = 2); // 0.2
  declareProperty( "nStripPhiHAD0wid" , m_nStripphiHADwid[0] = 2); // 0.196
  declareProperty( "nStripEtaHAD1wid" , m_nStripetaHADwid[1] = 2); // 0.2
  declareProperty( "nStripPhiHAD1wid" , m_nStripphiHADwid[1] = 2); // 0.196
  declareProperty( "nStripEtaHAD2wid" , m_nStripetaHADwid[2] = 1); // 0.2
  declareProperty( "nStripPhiHAD2wid" , m_nStripphiHADwid[2] = 2); // 0.196
 
  declareProperty( "nStripEtaHAD0nor" , m_nStripetaHADnor[0] = 4); // 0.4
  declareProperty( "nStripPhiHAD0nor" , m_nStripphiHADnor[0] = 4); // 0.39
  declareProperty( "nStripEtaHAD1nor" , m_nStripetaHADnor[1] = 4); // 0.4
  declareProperty( "nStripPhiHAD1nor" , m_nStripphiHADnor[1] = 4); // 0.39
  declareProperty( "nStripEtaHAD2nor" , m_nStripetaHADnor[2] = 2); // 0.4
  declareProperty( "nStripPhiHAD2nor" , m_nStripphiHADnor[2] = 4); // 0.39
 
 
 
}

Member Function Documentation

CellInNarCluster()

bool T2GeometryTool::CellInNarCluster ( const int nCaloType, const int nCaloSamp, const double & etaCell, const double & phiCell )

inlineoverridevirtual

Finds whether a cell is in a Narrow Sized Cluster.

Parameters

[in]	nCaloType,nCaloSamp	: See T2GeometryTool::EtaPhiRange.
[in]	etaCell	: Cell eta position.
[in]	phiCell	: Cell phi position.

Returns

cell is or not in this cluster.

Implements IT2GeometryTool.

Definition at line 426 of file T2GeometryTool.h.

                                                                         {
  // test for invalid calorimeter type
  if (nCaloType < 0 || nCaloType > 1) {
    ATH_MSG_ERROR("Invalid CaloType");
    return false;
  }
  if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
    ATH_MSG_ERROR("Invalid CaloSamp");
    return false;
  }
  int layer = nCaloType*4+nCaloSamp;
 
        int IetaPass=0;
        int IphiPass=0;
        if ( (etaCell > m_CellEtaNarMin[layer]) &&
           (etaCell < m_CellEtaNarMax[layer]) ) IetaPass=1;
        if ( ( m_CellPhiNarMax[layer] > m_CellPhiNarMin[layer] &&
           phiCell > m_CellPhiNarMin[layer]  &&
           phiCell < m_CellPhiNarMax[layer]  ) ||
         ( ( m_CellPhiNarMax[layer] < m_CellPhiNarMin[layer]  &&
            phiCell > 0 && phiCell > m_CellPhiNarMin[layer]  &&
            phiCell -2*M_PI < m_CellPhiNarMax[layer] ) || 
           ( phiCell < 0 && phiCell < m_CellPhiNarMax[layer] && 
             phiCell +2*M_PI > m_CellPhiNarMin[layer] ) )         
          )  IphiPass=1;
 
      if ( IetaPass == 1 && IphiPass == 1 ) return true;
        else return false;
} // end of CellInCluster

CellInNorCluster()

bool T2GeometryTool::CellInNorCluster ( const int nCaloType, const int nCaloSamp, const double & etaCell, const double & phiCell )

inlineoverridevirtual

Finds whether a cell is in a Normal Sized Cluster.

Parameters

[in]	nCaloType,nCaloSamp	: See T2GeometryTool::EtaPhiRange.
[in]	etaCell	: Cell eta position.
[in]	phiCell	: Cell phi position.

Returns

cell is or not in this cluster.

Implements IT2GeometryTool.

Definition at line 393 of file T2GeometryTool.h.

                                                                         {
 
  // test for invalid calorimeter type
  if (nCaloType < 0 || nCaloType > 1) {
    ATH_MSG_ERROR("Invalid CaloType");
    return false;
  }
  if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
    ATH_MSG_ERROR("Invalid CaloSamp");
    return false;
  }
  
  int layer = nCaloType*4+nCaloSamp;
 
        int IetaPass=0;
        int IphiPass=0;
        if ( (etaCell > m_CellEtaNorMin[layer]) &&
           (etaCell < m_CellEtaNorMax[layer] ) ) IetaPass=1;
        if ( ((m_CellPhiNorMax[layer] > m_CellPhiNorMin[layer]) &&
           (phiCell > m_CellPhiNorMin[layer]) &&
           (phiCell < m_CellPhiNorMax[layer]) ) ||
        ( (m_CellPhiNorMax[layer] < m_CellPhiNorMin[layer]  &&
            phiCell > 0 && phiCell > m_CellPhiNorMin[layer]  &&
            phiCell -2*M_PI < m_CellPhiNorMax[layer] ) ||
           ( phiCell < 0 && phiCell < m_CellPhiNorMax[layer] &&
             phiCell +2*M_PI > m_CellPhiNorMin[layer] ) )
          )  IphiPass=1;
 
      if ( IetaPass == 1 && IphiPass == 1 ) return true;
        else return false;
} // end of CellInCluster

CellInWidCluster()

bool T2GeometryTool::CellInWidCluster ( const int nCaloType, const int nCaloSamp, const double & etaCell, const double & phiCell )

inlineoverridevirtual

Finds whether a cell is in a Wide Sized Cluster.

Parameters

[in]	nCaloType,nCaloSamp	: See T2GeometryTool::EtaPhiRange.
[in]	etaCell	: Cell eta position.
[in]	phiCell	: Cell phi position.

Returns

cell is or not in this cluster.

Implements IT2GeometryTool.

Definition at line 457 of file T2GeometryTool.h.

                                                                         {
  // test for invalid calorimeter type
  if (nCaloType < 0 || nCaloType > 1) {
    ATH_MSG_ERROR("Invalid CaloType");
    return false;
  }
  if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
    ATH_MSG_ERROR("Invalid CaloSamp");
    return false;
  }
  int layer = nCaloType*4+nCaloSamp;
 
 
  int IetaPass=0;
  int IphiPass=0;
  if ( (etaCell > m_CellEtaWidMin[layer] ) &&
       (etaCell < m_CellEtaWidMax[layer] ) ) IetaPass=1;
  if ( ( m_CellPhiWidMax[layer] > m_CellPhiWidMin[layer]  &&
       phiCell > m_CellPhiWidMin[layer]  &&
       phiCell < m_CellPhiWidMax[layer]  ) ||
       ( ( m_CellPhiWidMax[layer] < m_CellPhiWidMin[layer]  &&
          phiCell > 0 && phiCell > m_CellPhiWidMin[layer]  &&
          phiCell -2*M_PI < m_CellPhiWidMax[layer] ) ||
         ( phiCell < 0 && phiCell < m_CellPhiWidMax[layer] &&
           phiCell +2*M_PI > m_CellPhiWidMin[layer] ) )
        )  IphiPass=1;
 
  if ( IetaPass == 1 && IphiPass == 1 ) return true;
  else return false;
} // end of CellInWidCluster

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

DeclareInterfaceID()

IT2GeometryTool::DeclareInterfaceID ( IT2GeometryTool , 1 , 0  )

inherited

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

EtaPhiRange()

int T2GeometryTool::EtaPhiRange ( const int nCaloType, const int nCaloSamp, const double & energyEta, const double & energyPhi )

inlineoverridevirtual

This prepares internal tables based on the cluster position.

Such tables should be used to determine whether a cells is inside clusters with different sizes

Parameters

[in]	nCaloType	should be 0 for LAr and 1 for TileCal.
[in]	nCaloSamp	should be calorimeter sample.
[in]	energyEta	: eta of hottest cell.
[in]	energyPhi	: phi of hottest cell.

Implements IT2GeometryTool.

Definition at line 330 of file T2GeometryTool.h.

                                                    {
 
  // test for invalid calorimeter type
  if (nCaloType < 0 || nCaloType > 1) {
    ATH_MSG_ERROR("Invalid CaloType");
    return 0;
  }
  if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
    ATH_MSG_ERROR("Invalid CaloSamp");
    return 0;
  }
  int layer = nCaloType*4+nCaloSamp;
 
  if( nCaloType == 0)
  {
 
    EtaPhiRangeInt(layer,m_nStripetaEMnar[nCaloSamp], 
        m_nStripphiEMnar[nCaloSamp],energyEta,energyPhi,0);
    EtaPhiRangeInt(layer,m_nStripetaEMwid[nCaloSamp],
        m_nStripphiEMwid[nCaloSamp],energyEta,energyPhi,1);
    EtaPhiRangeInt(layer,m_nStripetaEMnor[nCaloSamp],
        m_nStripphiEMnor[nCaloSamp],energyEta,energyPhi,2);
  }
  else 
  {
    EtaPhiRangeInt(layer,m_nStripetaHADnar[nCaloSamp], 
        m_nStripphiHADnar[nCaloSamp],energyEta,energyPhi,0);
    EtaPhiRangeInt(layer,m_nStripetaHADwid[nCaloSamp],
        m_nStripphiHADwid[nCaloSamp],energyEta,energyPhi,1);
    EtaPhiRangeInt(layer,m_nStripetaHADnor[nCaloSamp],
        m_nStripphiHADnor[nCaloSamp],energyEta,energyPhi,2);
 
  }
     
  // log removed, this is not an AlgTool
  ATH_MSG_DEBUG( "At end of EtaPhiRange. layer " << layer 
                 << " CellEtaNorMin : " << m_CellEtaNorMin[layer] << " CellEtaNorMax : " 
                 << m_CellEtaNorMax[layer] << " CellPhiNorMin : " << m_CellPhiNorMin[layer] 
                 << " CellPhiNorMax : " << m_CellPhiNorMax[layer] << " deta : " 
                 << m_CellEtaNorMin[layer] - energyEta << " dphi : "  
                 << m_CellPhiNorMin[layer] - energyPhi );
 
  ATH_MSG_DEBUG( " Nar variables ----> CellEtaNarMin : " 
                 << m_CellEtaNarMin[layer]
                 << " CellEtaNarMax : " << m_CellEtaNarMax[layer]  << " CellPhiNarMin : " 
                 << m_CellPhiNarMin[layer] << " CellPhiNarMax : " << m_CellPhiNarMax[layer]
                 << " detaNar : " << m_CellEtaNarMin[layer] - energyEta << " dphiNar : "
                 << m_CellPhiNarMin[layer] - energyPhi );
 
  ATH_MSG_DEBUG( " Wid variables ----> CellEtaWidMin : " 
                 << m_CellEtaWidMin[layer] << " CellEtaWidMax : " << m_CellEtaWidMax[layer]
                 << " CellPhiWidMin : "  << m_CellPhiWidMin[layer] << " CellPhiWidMax : " 
                 << m_CellPhiWidMax[layer] << " detaWid : " 
                 << m_CellEtaWidMin[layer] - energyEta << " dphiWid : "
                 << m_CellPhiWidMin[layer] - energyPhi );
 
 
 
 
    return 0;
} // End of EtaPhiRange

EtaPhiRangeInt()

int T2GeometryTool::EtaPhiRangeInt ( const int layer, const int & nStripEta, const int & nStripPhi, const double & energyEta, const double & energyPhi, const int index_size )

inlineoverrideprivatevirtual

Implements IT2GeometryTool.

Definition at line 202 of file T2GeometryTool.h.

                                                                               {
 
  double aeta = fabs(energyEta);
  double dgra = m_etaGran[layer][0]; // eta granularity at eta=0
  int netaregions=m_netareg[layer]; // number of different eta regions
  double CellEtaMin=0.;
  double CellEtaMax=0.;
  double CellPhiMin=0.;
  double CellPhiMax=0.;
 
  CellEtaMin = aeta - dgra*nStripEta*0.5;
  int EtaMinSign = 1;
  if ( CellEtaMin < 0. ) {
    EtaMinSign = -1;
    CellEtaMin = fabs(CellEtaMin );  // need a non-negative value
  }
  ATH_MSG_DEBUG( "At start of EtaPhiRange. layer " << layer <<
                 " aeta " << aeta << " dgra " << dgra << " netaregions " << netaregions
                 << " First CellEtaMin " << CellEtaMin );
 
  // Go by cell edges not centres since boundaries are at edges.
  // Determine correct cell edges below, if seed is exactly at a cell
  // edge rounding uncertainties could theoretically result in wrong
  // number of cells being used, this is not checked.
 
  // When dealing with a non-integer number of cells (due to including cells
  // with a different granularity then the central region) round up to the next
  // integer number of cells if at least 30% of the cell would be in the
  // desired range (i.e. error on the side of including too much rather
  // then too little).  This is an arbitrary choice and not proven to be ideal.
 
  int irmin;
  int irmax;
  double dgra1;
 
if( CellEtaMin < m_etareg[layer][netaregions]){
  for (irmin=0;irmin<netaregions;irmin++) {
 
    if (CellEtaMin >= m_etareg[layer][irmin] &&
        CellEtaMin < m_etareg[layer][irmin+1]) {
 
      dgra1 = m_etaGran[layer][irmin];
      CellEtaMin = (CellEtaMin-m_etareg[layer][irmin])/dgra1; //number of cells
      if(CellEtaMin - trunc(CellEtaMin) > 0.7 ) CellEtaMin++; //round down
      CellEtaMin = dgra1*trunc(CellEtaMin) + m_etareg[layer][irmin];
 
      CellEtaMin = CellEtaMin*EtaMinSign; // correct for edge that spans 0
      CellEtaMax = aeta + dgra*nStripEta * 0.5;
 
 
        for (irmax=irmin;irmax<netaregions;irmax++) {
          if (CellEtaMax > m_etareg[layer][irmax] &&
              CellEtaMax <= m_etareg[layer][irmax+1]) {
            dgra1 =  m_etaGran[layer][irmax];
            CellEtaMax = (CellEtaMax-m_etareg[layer][irmax])/dgra1; //number of cells
            if (CellEtaMax-trunc(CellEtaMax) >= 0.3) CellEtaMax++; // round up
            CellEtaMax = dgra1*trunc(CellEtaMax)+m_etareg[layer][irmax];
          }
        }  // for (irmax=irmin+1;irmax<netaregions;irmax++)
    }  // if (*CellEtaMin >= m_etareg[layer][irmin]...
  }  // for (irmin=0;irmin<=netaregions;irmin++)
} //  if (CellEtaMin < m_etareg[layer][netaregions])
 
    // set to correct sign
  double tempeta;
  if (energyEta < 0) {
    tempeta=-CellEtaMin;
    CellEtaMin=-CellEtaMax;
    CellEtaMax=tempeta;
 
  }
 
  // set phi range based only on phi granularity in the eta region
  // of CellEtaMin, this may not be the optimum method when the eta
  // range spans a change in phi granularity but the only place that
  // might be a concern is the crack region of EM samp 1 since the
  // other phi gran. changes occur at eta > 2.5
 
  double dgrap= m_phiGran[layer][0]; // phi granularity at eta=0
  int nphiregions=m_nphireg[layer]; // number of different phi regions
  CellPhiMin = energyPhi - dgrap*nStripPhi*0.5;
 
if (fabs(CellEtaMin) < m_phireg[layer][nphiregions]){
  for (irmin=0;irmin<nphiregions;irmin++) {
 
    if (fabs(CellEtaMin) >= m_phireg[layer][irmin] &&
        fabs(CellEtaMin) < m_phireg[layer][irmin+1]) {
 
      dgra1 = m_phiGran[layer][irmin];
      CellPhiMin = dgra1*rint(CellPhiMin/dgra1);
      CellPhiMax = CellPhiMin + dgrap*nStripPhi;
    }  // if *(CellEtaMin >= m_phireg[layer][irmin]...
  }  // for (irmin=0;irmin<=nphiregions;irmin++)
} // if (fabs(CellEtaMin) < m_phireg[layer][nphiregions]){
 
  if (CellPhiMax >   M_PI) CellPhiMax -= 2*M_PI;
  if (CellPhiMin <= -M_PI) CellPhiMin += 2*M_PI;
  if (CellPhiMin >   M_PI) CellPhiMin -= 2*M_PI;
  if (CellPhiMax <= -M_PI) CellPhiMax += 2*M_PI;
 
  if( index_size == 0)
  {
    m_CellEtaNarMax[layer] = CellEtaMax;
    m_CellEtaNarMin[layer] = CellEtaMin;
    m_CellPhiNarMax[layer] = CellPhiMax;
    m_CellPhiNarMin[layer] = CellPhiMin;
  }
  else if (index_size == 1)
  {
    m_CellEtaWidMax[layer] = CellEtaMax;
    m_CellEtaWidMin[layer] = CellEtaMin;
    m_CellPhiWidMax[layer] = CellPhiMax;
    m_CellPhiWidMin[layer] = CellPhiMin;
  }
  else if (index_size == 2)
  {
    m_CellEtaNorMax[layer] = CellEtaMax;
    m_CellEtaNorMin[layer] = CellEtaMin;
    m_CellPhiNorMax[layer] = CellPhiMax;
    m_CellPhiNorMin[layer] = CellPhiMin;
  }
 
 
  return 0;
} // End of EtaPhiRangeInt method

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

initialize()

StatusCode T2GeometryTool::initialize ( )

overridevirtual

Initialize method.

Definition at line 67 of file T2GeometryTool.cxx.

{
  ATH_MSG_DEBUG("in initialize() by T2GeometryTool");
 return StatusCode::SUCCESS; 
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

nStripetaEMnar()

int T2GeometryTool::nStripetaEMnar ( const int layer ) const

inlineoverridevirtual

Get number of cells in eta for Narrow Sized Cluster in EM calorimeter.

Parameters

[in]

layer

integer to identify the EM sampling

Implements IT2GeometryTool.

Definition at line 490 of file T2GeometryTool.h.

        { return m_nStripetaEMnar[layer]; }

nStripetaEMnor()

int T2GeometryTool::nStripetaEMnor ( const int layer ) const

inlineoverridevirtual

Get number of cells in eta for Normal Sized Cluster in EM calorimeter.

Parameters

[in]

layer

integer to identify the EM sampling

Implements IT2GeometryTool.

Definition at line 500 of file T2GeometryTool.h.

        { return m_nStripetaEMnor[layer]; }

nStripetaEMwid()

int T2GeometryTool::nStripetaEMwid ( const int layer ) const

inlineoverridevirtual

Get number of cells in eta for Wide Sized Cluster in EM calorimeter.

Parameters

[in]

layer

integer to identify the EM sampling

Implements IT2GeometryTool.

Definition at line 495 of file T2GeometryTool.h.

        { return m_nStripetaEMwid[layer]; }

nStripetaHADnar()

int T2GeometryTool::nStripetaHADnar ( const int layer ) const

inlineoverridevirtual

Get number of cells in eta for Narrow Sized Cluster in HAD calorimeter.

Parameters

[in]

layer

integer to identify the HAD sampling

Implements IT2GeometryTool.

Definition at line 505 of file T2GeometryTool.h.

        { return m_nStripetaHADnar[layer]; }

nStripetaHADnor()

int T2GeometryTool::nStripetaHADnor ( const int layer ) const

inlineoverridevirtual

Get number of cells in eta for Normal Sized Cluster in HAD calorimeter.

Parameters

[in]

layer

integer to identify the HAD sampling

Implements IT2GeometryTool.

Definition at line 515 of file T2GeometryTool.h.

        { return m_nStripetaHADnor[layer]; }

nStripetaHADwid()

int T2GeometryTool::nStripetaHADwid ( const int layer ) const

inlineoverridevirtual

Get number of cells in eta for Wide Sized Cluster in HADcalorimeter.

Parameters

[in]

layer

integer to identify the HAD sampling

Implements IT2GeometryTool.

Definition at line 510 of file T2GeometryTool.h.

        { return m_nStripetaHADwid[layer]; }

nStripphiEMnar()

int T2GeometryTool::nStripphiEMnar ( const int layer ) const

inlineoverridevirtual

Get number of cells in phi for Narrow Sized Cluster in EM calorimeter.

Parameters

[in]

layer

integer to identify the EM sampling

Implements IT2GeometryTool.

Definition at line 492 of file T2GeometryTool.h.

        { return m_nStripphiEMnar[layer]; }

nStripphiEMnor()

int T2GeometryTool::nStripphiEMnor ( const int layer ) const

inlineoverridevirtual

Get number of cells in phi for Normal Sized Cluster in EM calorimeter.

Parameters

[in]

layer

integer to identify the EM sampling

Implements IT2GeometryTool.

Definition at line 502 of file T2GeometryTool.h.

        { return m_nStripphiEMnor[layer]; }

nStripphiEMwid()

int T2GeometryTool::nStripphiEMwid ( const int layer ) const

inlineoverridevirtual

Get number of cells in phi for Wide Sized Cluster in EM calorimeter.

Parameters

[in]

layer

integer to identify the EM sampling

Implements IT2GeometryTool.

Definition at line 497 of file T2GeometryTool.h.

        { return m_nStripphiEMwid[layer]; }

nStripphiHADnar()

int T2GeometryTool::nStripphiHADnar ( const int layer ) const

inlineoverridevirtual

Get number of cells in phi for Narrow Sized Cluster in HAD calorimeter.

Parameters

[in]

layer

integer to identify the HAD sampling

Implements IT2GeometryTool.

Definition at line 507 of file T2GeometryTool.h.

        { return m_nStripphiHADnar[layer]; }

nStripphiHADnor()

int T2GeometryTool::nStripphiHADnor ( const int layer ) const

inlineoverridevirtual

Get number of cells in phi for Normal Sized Cluster in HAD calorimeter.

Parameters

[in]

layer

integer to identify the HAD sampling

Implements IT2GeometryTool.

Definition at line 517 of file T2GeometryTool.h.

        { return m_nStripphiHADnor[layer]; }

nStripphiHADwid()

int T2GeometryTool::nStripphiHADwid ( const int layer ) const

inlineoverridevirtual

Get number of cells in phi for Wide Sized Cluster in HAD calorimeter.

Parameters

[in]

layer

integer to identify the HAD sampling

Implements IT2GeometryTool.

Definition at line 512 of file T2GeometryTool.h.

        { return m_nStripphiHADwid[layer]; }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_CellEtaNarMax

double T2GeometryTool::m_CellEtaNarMax[7] {}

private

Definition at line 161 of file T2GeometryTool.h.

{}, m_CellEtaNarMax[7]{};

m_CellEtaNarMin

double T2GeometryTool::m_CellEtaNarMin[7] {}

private

Narrow cluster variables.

Definition at line 161 of file T2GeometryTool.h.

{}, m_CellEtaNarMax[7]{};

m_CellEtaNorMax

double T2GeometryTool::m_CellEtaNorMax[7] {}

private

Definition at line 158 of file T2GeometryTool.h.

{}, m_CellEtaNorMax[7]{};

m_CellEtaNorMin

double T2GeometryTool::m_CellEtaNorMin[7] {}

private

Normal cluster variables.

Definition at line 158 of file T2GeometryTool.h.

{}, m_CellEtaNorMax[7]{};

m_CellEtaWidMax

double T2GeometryTool::m_CellEtaWidMax[7] {}

private

Definition at line 164 of file T2GeometryTool.h.

{}, m_CellEtaWidMax[7]{};

m_CellEtaWidMin

double T2GeometryTool::m_CellEtaWidMin[7] {}

private

Wider cluster variables.

Definition at line 164 of file T2GeometryTool.h.

{}, m_CellEtaWidMax[7]{};

m_CellPhiNarMax

double T2GeometryTool::m_CellPhiNarMax[7] {}

private

Definition at line 162 of file T2GeometryTool.h.

{}, m_CellPhiNarMax[7]{};

m_CellPhiNarMin

double T2GeometryTool::m_CellPhiNarMin[7] {}

private

Definition at line 162 of file T2GeometryTool.h.

{}, m_CellPhiNarMax[7]{};

m_CellPhiNorMax

double T2GeometryTool::m_CellPhiNorMax[7] {}

private

Definition at line 159 of file T2GeometryTool.h.

{}, m_CellPhiNorMax[7]{};

m_CellPhiNorMin

double T2GeometryTool::m_CellPhiNorMin[7] {}

private

Definition at line 159 of file T2GeometryTool.h.

{}, m_CellPhiNorMax[7]{};

m_CellPhiWidMax

double T2GeometryTool::m_CellPhiWidMax[7] {}

private

Definition at line 165 of file T2GeometryTool.h.

{}, m_CellPhiWidMax[7]{};

m_CellPhiWidMin

double T2GeometryTool::m_CellPhiWidMin[7] {}

private

Definition at line 165 of file T2GeometryTool.h.

{}, m_CellPhiWidMax[7]{};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_etaGran

const double T2GeometryTool::m_etaGran

staticprivate

Initial value:

= {
      {0.025   ,0.024,0.025   ,0.0      ,0.0    ,0.0  ,0.0},
      {0.003125,0.025,0.003125,0.0041667,0.00625,0.025,0.1},
      {0.025   ,0.05 ,0.025   ,0.1      ,0.0    ,0.0  ,0.0},
      {0.05    ,0.001,0.05    ,0.0      ,0.0    ,0.0  ,0.0},
      {0.1     ,0.2  ,0.0     ,0.0      ,0.0    ,0.0  ,0.0},
      {0.1     ,0.2  ,0.0     ,0.0      ,0.0    ,0.0  ,0.0},
      {0.2     ,0.1  ,0.2     ,0.0      ,0.0    ,0.0  ,0.0}}

Definition at line 85 of file T2GeometryTool.h.

m_etareg

const double T2GeometryTool::m_etareg

staticprivate

Initial value:

= {
       {0.0,1.4 ,1.5  ,1.8,0.0,0.0,0.0,0.0},
       {0.0,1.4 ,1.5  ,1.8,2.0,2.4,2.5,3.2},
       {0.0,1.4 ,1.425,2.5,3.2,0.0,0.0,0.0},
       {0.0,1.35,1.5  ,2.5,0.0,0.0,0.0,0.0},
       {0.0,2.5 ,3.2  ,0.0,0.0,0.0,0.0,0.0},
       {0.0,2.5 ,3.2  ,0.0,0.0,0.0,0.0,0.0},
       {0.0,1.5 ,2.5  ,3.2,0.0,0.0,0.0,0.0}}

Definition at line 77 of file T2GeometryTool.h.

                    : integer to identify the EM sampling
  */
 
        virtual
    int nStripetaEMnar(const int layer) const override;

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_netareg

const int T2GeometryTool::m_netareg = {3,7,4,3,2,2,3}

staticprivate

Granularity and eta ranges relevant for above energy sums.

Definition at line 75 of file T2GeometryTool.h.

m_nphireg

const int T2GeometryTool::m_nphireg = {1,3,2,1,2,2,2}

staticprivate

Definition at line 95 of file T2GeometryTool.h.

m_nStripetaEMnar

int T2GeometryTool::m_nStripetaEMnar[4] {}

private

Reference changed to vectors.

Definition at line 174 of file T2GeometryTool.h.

m_nStripetaEMnor

int T2GeometryTool::m_nStripetaEMnor[4] {}

private

Definition at line 180 of file T2GeometryTool.h.

m_nStripetaEMwid

int T2GeometryTool::m_nStripetaEMwid[4] {}

private

Definition at line 177 of file T2GeometryTool.h.

m_nStripetaHADnar

int T2GeometryTool::m_nStripetaHADnar[3] {}

private

Definition at line 183 of file T2GeometryTool.h.

m_nStripetaHADnor

int T2GeometryTool::m_nStripetaHADnor[3] {}

private

Definition at line 189 of file T2GeometryTool.h.

m_nStripetaHADwid

int T2GeometryTool::m_nStripetaHADwid[3] {}

private

Definition at line 186 of file T2GeometryTool.h.

m_nStripphiEMnar

int T2GeometryTool::m_nStripphiEMnar[4] {}

private

Definition at line 175 of file T2GeometryTool.h.

m_nStripphiEMnor

int T2GeometryTool::m_nStripphiEMnor[4] {}

private

Definition at line 181 of file T2GeometryTool.h.

m_nStripphiEMwid

int T2GeometryTool::m_nStripphiEMwid[4] {}

private

Definition at line 178 of file T2GeometryTool.h.

m_nStripphiHADnar

int T2GeometryTool::m_nStripphiHADnar[3] {}

private

Definition at line 184 of file T2GeometryTool.h.

m_nStripphiHADnor

int T2GeometryTool::m_nStripphiHADnor[3] {}

private

Definition at line 190 of file T2GeometryTool.h.

m_nStripphiHADwid

int T2GeometryTool::m_nStripphiHADwid[3] {}

private

Definition at line 187 of file T2GeometryTool.h.

m_phiGran

const double T2GeometryTool::m_phiGran

staticprivate

Initial value:

Definition at line 106 of file T2GeometryTool.h.

m_phireg

const double T2GeometryTool::m_phireg

staticprivate

Initial value:

Definition at line 97 of file T2GeometryTool.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• T2GeometryTool.h
• T2GeometryTool.cxx

---

Generated on for ATLAS Offline Software by  1.14.0