InDet::TRTGeometryManagerTool Class Reference

#include <TRTGeometryManagerTool.h>

Inheritance diagram for InDet::TRTGeometryManagerTool:

Collaboration diagram for InDet::TRTGeometryManagerTool:

Public Member Functions
	TRTGeometryManagerTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~TRTGeometryManagerTool ()
StatusCode	initialize ()
StatusCode	finalize ()
int	ReadGeometry (int solveLevel)
	read the geometry Method is called from the main AlignAlg to build the geometry based on the requested level of alignment.
void	PrintModPositions ()
void	StoreModPositions (std::string)
void	buildGeometry ()
	builds geometry for TRT alignment
void	addModuleParameters (Trk::AlignModule *module, DataVector< DataVector< Trk::AlignPar > > *allFullModPars, DataVector< DataVector< Trk::AlignPar > > *allActiveModPars)
	adds alignment parameters for the module checks for active parameters and calls setSigmas()
bool	checkAlignLevel ()
	check whether the alignment level is correct
bool	checkAlignLevelBarrel ()
	check whether the Barrel alignment level is correct
bool	checkAlignLevelEndcaps ()
	check whether the Endcap alignment level is correct
void	setAlignLevel (int level)
	set alignment level
virtual void	setNtuple (TFile *)
	Sets output ntuple file, to be implemented if needed for detector-specific implementation.
virtual void	setLogStream (std::ostream *os)
	sets the output stream for the logfile
virtual int	getNextIDHash ()
	get next free IDHash usable for new AlignModule
virtual void	setFirstIDHash (unsigned int idhash)
	set first free IDHash usable for new AlignModule
virtual void	setModuleList (Trk::AlignModuleList *list)
	set pointer to module list to which the modules should be added to
virtual void	setHashMaps (std::vector< Trk::AlignModuleList * > *maps)
	set pointer to vector of hashMaps to which the elements should be added to
virtual int	alignLevel ()
	get alignment level
virtual int	alignLevelBarrel ()
	get alignment level
virtual void	setAlignLevelBarrel (int level)
	set alignment level
virtual int	alignLevelEndcaps ()
	get alignment level
virtual void	setAlignLevelEndcaps (int level)
	set alignment level
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

Static Public Member Functions
static const InterfaceID &	interfaceID ()

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.

Protected Attributes
Trk::AlignModuleList *	m_alignModuleListPtr = nullptr
	pointer to module list to which the modules are added
std::vector< Trk::AlignModuleList * > *	m_idHashToAlignModuleMapsPtr = nullptr
	pointer to vector of hashMaps to which the elements are added
int	m_hashCounter = 0
	variable for setting the idHash of the AlignModules
std::ostream *	m_logStream = nullptr
	logfile output stream
int	m_alignLevel = 0
	alignment level
int	m_alignLevelBarrel = 0
	alignment level
int	m_alignLevelEndcaps = 0
	alignment level

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	buildL0 ()
	builds geometry for L1 TRT alignment Creates 1 module - the whole TRT
void	buildL1Barrel ()
	builds geometry for L1 TRT Barrel alignment Creates single module for barrel
void	buildL1Endcaps ()
	builds geometry for L1 TRT Endcap alignment Creates 2 modules, one for each endcap
void	buildL2Barrel ()
	builds geometry for L2 TRT Barrel alignment Creates 96 modules - 3 layers with 32 modules per layer
void	buildL2Endcaps ()
	builds geometry for L2 TRT Endcapalignment Creates 80 modules - 40 rings per endcap
void	buildL2EndcapsOLD ()
	builds geometry for L2 TRT Endcapalignment Creates 28 modules - 14 wheels per endcap
void	buildL3Barrel ()
	builds geometry for L3 TRT Barrel alignment Creates many modules - one module per straw
void	buildL3Endcap ()
	builds geometry for L3 TRT Endcap alignment Creates many modules - one module per straw
void	setSigmas (Trk::AlignModule *mod, DataVector< Trk::AlignPar > *modPars)
	sets sigmas for modules
bool	moduleSelected (Trk::AlignModule *mod)
	check wheather module is selected for module pointer
bool	moduleSelected (unsigned long long modId)
	check wheather module is selected for module id
bool	moduleSelectedBarrel (unsigned int iLayer, unsigned int iPhi)
	check wheather barrel module is selected by layer and phi
bool	moduleSelectedEndcap (unsigned int iwheel, unsigned int iphi)
	check wheather we are aligning the different endcaps and sectors
void	dumpGeometry ()
	print basic geometry info to screen
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
const InDetDD::TRT_DetectorManager *	m_trtDetManager
	pointer to TRT detector manager
const TRT_ID *	m_idHelper
ToolHandle< Trk::IAlignModuleTool >	m_alignModuleTool
	pointer to AlignModuleTool
Trk::AlignModuleList	m_alignModuleList
Trk::AlignModuleList	m_idHashToAlignModuleMap
std::vector< Trk::AlignModuleList * >	m_idHashToAlignModuleMaps
DataVector< DataVector< Trk::AlignPar > > *	m_alignParList
DataVector< DataVector< Trk::AlignPar > > *	m_fullAlignParList
bool	m_alignBarrel
bool	m_alignEndcaps
bool	m_alignBarrelX
bool	m_alignBarrelY
bool	m_alignBarrelZ
bool	m_alignBarrelRotX
bool	m_alignBarrelRotY
bool	m_alignBarrelRotZ
bool	m_alignEndcapX
bool	m_alignEndcapY
bool	m_alignEndcapZ
bool	m_alignEndcapRotX
bool	m_alignEndcapRotY
bool	m_alignEndcapRotZ
double	m_sigmaBarrelX
double	m_sigmaBarrelY
double	m_sigmaBarrelZ
double	m_sigmaBarrelRotX
double	m_sigmaBarrelRotY
double	m_sigmaBarrelRotZ
double	m_sigmaEndcapX
double	m_sigmaEndcapY
double	m_sigmaEndcapZ
double	m_sigmaEndcapRotX
double	m_sigmaEndcapRotY
double	m_sigmaEndcapRotZ
double	m_softcutBarrelX
double	m_softcutBarrelY
double	m_softcutBarrelZ
double	m_softcutBarrelRotX
double	m_softcutBarrelRotY
double	m_softcutBarrelRotZ
double	m_softcutEndcapX
double	m_softcutEndcapY
double	m_softcutEndcapZ
double	m_softcutEndcapRotX
double	m_softcutEndcapRotY
double	m_softcutEndcapRotZ
std::vector< int >	m_barrelPhiSectors
	list of barrel phi sectors to align (if empty, align everything)
std::vector< int >	m_barrelLayers
	list of barrel layers to align (if empty, align everything)
int	m_endcap
	endcap_id of the endcap to align (-2 or 2 for endcap A or C, 0 for both)
std::vector< int >	m_endcapPhiSectors
	list of endcap phi sectors to align (if empty, align everything)
std::vector< int >	m_endcapWheels
	list of endcap wheels to align (if empty, align everything)
bool	m_doModuleSelection
std::vector< unsigned long long >	m_moduleSelection
bool	m_dumpGeometry
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

Detailed Description

Definition at line 42 of file TRTGeometryManagerTool.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

TRTGeometryManagerTool()

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

Definition at line 31 of file TRTGeometryManagerTool.cxx.

  : AthAlgTool(type,name,parent)
      , m_trtDetManager(nullptr)
      , m_idHelper()
      , m_alignModuleTool("Trk::AlignModuleTool/AlignModuleTool")
      , m_idHashToAlignModuleMaps(Trk::AlignModule::NDetectorTypes,(Trk::AlignModuleList*)nullptr)
      , m_alignParList(nullptr)
      , m_fullAlignParList(nullptr)
{
   declareInterface<IGeometryManagerTool>(this);
   declareProperty("AlignModuleTool",    m_alignModuleTool);
 
   declareProperty("AlignBarrel",        m_alignBarrel        = true);
   declareProperty("AlignBarrelX",       m_alignBarrelX       = true);
   declareProperty("AlignBarrelY",       m_alignBarrelY       = true);
   declareProperty("AlignBarrelZ",       m_alignBarrelZ       = false);
   declareProperty("AlignBarrelRotX",    m_alignBarrelRotX    = true);
   declareProperty("AlignBarrelRotY",    m_alignBarrelRotY    = true);
   declareProperty("AlignBarrelRotZ",    m_alignBarrelRotZ    = true);
 
   declareProperty("AlignEndcaps",       m_alignEndcaps       = true);
   declareProperty("AlignEndcapX",       m_alignEndcapX       = true);
   declareProperty("AlignEndcapY",       m_alignEndcapY       = true);
   declareProperty("AlignEndcapZ",       m_alignEndcapZ       = true);
   declareProperty("AlignEndcapRotX",    m_alignEndcapRotX    = true);
   declareProperty("AlignEndcapRotY",    m_alignEndcapRotY    = true);
   declareProperty("AlignEndcapRotZ",    m_alignEndcapRotZ    = true);
 
   declareProperty("SetSigmaBarrelX",    m_sigmaBarrelX       = 1.);
   declareProperty("SetSigmaBarrelY",    m_sigmaBarrelY       = 1.);
   declareProperty("SetSigmaBarrelZ",    m_sigmaBarrelZ       = 1.);
   declareProperty("SetSigmaBarrelRotX", m_sigmaBarrelRotX    = 0.001);
   declareProperty("SetSigmaBarrelRotY", m_sigmaBarrelRotY    = 0.001);
   declareProperty("SetSigmaBarrelRotZ", m_sigmaBarrelRotZ    = 0.001);
 
   declareProperty("SetSigmaEndcapX",    m_sigmaEndcapX       = 1.);
   declareProperty("SetSigmaEndcapY",    m_sigmaEndcapY       = 1.);
   declareProperty("SetSigmaEndcapZ",    m_sigmaEndcapZ       = 1.);
   declareProperty("SetSigmaEndcapRotX", m_sigmaEndcapRotX    = 0.001);
   declareProperty("SetSigmaEndcapRotY", m_sigmaEndcapRotY    = 0.001);
   declareProperty("SetSigmaEndcapRotZ", m_sigmaEndcapRotZ    = 0.001);
 
   declareProperty("SetSoftCutBarrelX",    m_softcutBarrelX       = 1.);
   declareProperty("SetSoftCutBarrelY",    m_softcutBarrelY       = 1.);
   declareProperty("SetSoftCutBarrelZ",    m_softcutBarrelZ       = 1.);
   declareProperty("SetSoftCutBarrelRotX", m_softcutBarrelRotX    = 1.);
   declareProperty("SetSoftCutBarrelRotY", m_softcutBarrelRotY    = 1.);
   declareProperty("SetSoftCutBarrelRotZ", m_softcutBarrelRotZ    = 1.);
 
   declareProperty("SetSoftCutEndcapX",    m_softcutEndcapX       = 1.);
   declareProperty("SetSoftCutEndcapY",    m_softcutEndcapY       = 1.);
   declareProperty("SetSoftCutEndcapZ",    m_softcutEndcapZ       = 1.);
   declareProperty("SetSoftCutEndcapRotX", m_softcutEndcapRotX    = 1.);
   declareProperty("SetSoftCutEndcapRotY", m_softcutEndcapRotY    = 1.);
   declareProperty("SetSoftCutEndcapRotZ", m_softcutEndcapRotZ    = 1.);
 
   declareProperty("AlignmentLevel",        m_alignLevel        = 1);
   declareProperty("AlignmentLevelBarrel",  m_alignLevelBarrel  = -1);
   declareProperty("AlignmentLevelEndcaps", m_alignLevelEndcaps = -1);
 
   declareProperty("AlignBarrelPhiSectors", m_barrelPhiSectors);
   declareProperty("AlignBarrelLayers",     m_barrelLayers);
 
   declareProperty("AlignSingleEndcap",     m_endcap = 0);
   declareProperty("AlignEndcapWheels",     m_endcapWheels);
   declareProperty("AlignEndcapPhiSectors", m_endcapPhiSectors);
 
   declareProperty("doModuleSelection",  m_doModuleSelection  = false);
   declareProperty("ModuleSelection",    m_moduleSelection);
 
   declareProperty("DumpGeometry",       m_dumpGeometry = true);
 
   m_hashCounter = 0;
   m_logStream = nullptr;
}

~TRTGeometryManagerTool()

InDet::TRTGeometryManagerTool::~TRTGeometryManagerTool ( )

virtual

Definition at line 110 of file TRTGeometryManagerTool.cxx.

{
   ATH_MSG_DEBUG("deleting alignModuleList");
   for (int i=0;i<(int)m_alignModuleList.size();i++)
      delete m_alignModuleList[i];
   m_alignModuleList.clear();
 
   ATH_MSG_DEBUG("deleting fullAlignParList");
   delete m_fullAlignParList;
   ATH_MSG_DEBUG("deleting alignParList");
   delete m_alignParList;
}

Member Function Documentation

addModuleParameters()

void InDet::TRTGeometryManagerTool::addModuleParameters ( Trk::AlignModule * module, DataVector< DataVector< Trk::AlignPar > > * allFullModPars, DataVector< DataVector< Trk::AlignPar > > * allActiveModPars )

virtual

adds alignment parameters for the module checks for active parameters and calls setSigmas()

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 1168 of file TRTGeometryManagerTool.cxx.

{
   // prepare all parameters
   DataVector<Trk::AlignPar> * fullModPars = new DataVector<Trk::AlignPar>(SG::OWN_ELEMENTS);
   fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::TransX));
   fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::TransY));
   fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::TransZ));
   fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::RotX));
   fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::RotY));
   fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::RotZ));
 
   // set sigmas
   setSigmas(module,fullModPars);
 
   // select active parameters based on jobOption properties
   DataVector<Trk::AlignPar> * activeModPars = new DataVector<Trk::AlignPar>(SG::VIEW_ELEMENTS);
   for(unsigned int ipar=0;ipar<fullModPars->size();++ipar) {
 
      Identifier modID = module->identify();
 
      if(m_idHelper->is_barrel(modID)) {
         ATH_MSG_DEBUG("Barrel module with id "<<modID);
         if( (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransX && m_alignBarrelX)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransY && m_alignBarrelY)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransZ && m_alignBarrelZ)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotX && m_alignBarrelRotX)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotY && m_alignBarrelRotY)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotZ && m_alignBarrelRotZ) ) {
            ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is now active");
            activeModPars->push_back(fullModPars->at(ipar));
         }
         else
            ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is NOT active");
      }
      else {
         ATH_MSG_DEBUG("Endcap module with id "<<modID);
         if( (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransX && m_alignEndcapX)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransY && m_alignEndcapY)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransZ && m_alignEndcapZ)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotX && m_alignEndcapRotX)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotY && m_alignEndcapRotY)
            || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotZ && m_alignEndcapRotZ) ) {
            ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is now active");
            activeModPars->push_back(fullModPars->at(ipar));
         }
         else
            ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is NOT active");
      }
   }
 
   // add to parameter lists
   allFullModPars->push_back(fullModPars);
   allActiveModPars->push_back(activeModPars);
}

alignLevel()

virtual int Trk::IGeometryManagerTool::alignLevel ( )

inlinevirtualinherited

get alignment level

Definition at line 80 of file IGeometryManagerTool.h.

{ return m_alignLevel; }

alignLevelBarrel()

virtual int Trk::IGeometryManagerTool::alignLevelBarrel ( )

inlinevirtualinherited

get alignment level

Definition at line 86 of file IGeometryManagerTool.h.

{ return m_alignLevelBarrel; }

alignLevelEndcaps()

virtual int Trk::IGeometryManagerTool::alignLevelEndcaps ( )

inlinevirtualinherited

get alignment level

Definition at line 95 of file IGeometryManagerTool.h.

{ return m_alignLevelEndcaps; }

buildGeometry()

void InDet::TRTGeometryManagerTool::buildGeometry ( )

virtual

builds geometry for TRT alignment

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 286 of file TRTGeometryManagerTool.cxx.

{
   ATH_MSG_INFO("Preparing the TRT geometry");
 
   if(m_alignLevel == 0)
      buildL0();
   else {
 
      if(m_alignBarrel)
         switch(m_alignLevelBarrel) {
            case 1:
               buildL1Barrel();
               break;
            case 2:
               buildL2Barrel();
               break;
            case 3:
               buildL3Barrel();
               m_alignModuleTool->setSubDetElement(Trk::AlignModule::TRT,true);
               break;
            default:
               break;
         }
 
      if(m_alignEndcaps)
         switch(m_alignLevelEndcaps) {
            case 1:
               buildL1Endcaps();
               break;
            case 2:
               buildL2Endcaps();
               break;
            case 200:
               buildL2EndcapsOLD();
               break;
            case 3:
               buildL3Endcap();
               m_alignModuleTool->setSubDetElement(Trk::AlignModule::TRT,true);
               break;
            default:
               break;
         }
   }
}

buildL0()

void InDet::TRTGeometryManagerTool::buildL0 ( )

private

builds geometry for L1 TRT alignment Creates 1 module - the whole TRT

Definition at line 332 of file TRTGeometryManagerTool.cxx.

{
   ATH_MSG_INFO("Preparing the TRT geometry for L0: whole TRT");
 
   // ========================================
   // get all modules for Level 0 alignment
   // Level 0 is the whole TRT, i.e., 1 AlignModule in total
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // for L1 alignment the alignment frame is equal to the global frame
   // and since the TRT detector element positions are also stored in
   // the global frame in DB, transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // use the 'real' identifier of the TRT barrel for TRT L0 module
   Trk::AlignModule * trt = new Trk::AlignModule(this);
   trt->setIdHash(getNextIDHash());
   trt->setName("TRT");
   trt->setIdentifier(m_idHelper->barrel_ec_id(-1));
 
   if(!moduleSelected(trt)) {
      ATH_MSG_DEBUG("Module "<<trt->name()<<" NOT selected");
      delete trt;
      return;
   }
 
   // get maximum number of elements from the helper
   unsigned int trtmaxHash = m_idHelper->straw_layer_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<trtmaxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)(trtmaxHash),nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
   // ================================================================
   // loop over TRT elements and add them to respective alignModules
   // ================================================================
   for (unsigned int index = 0; index < trtmaxHash; index++) {
      IdentifierHash idHash = index;
      Identifier id = m_idHelper->layer_id(idHash);
 
      // get the element via hash
      const TRT_BaseElement * element = m_trtDetManager->getElement(idHash);
      if (element) {
         // get element location for debugging
//         HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
//         ATH_MSG_DEBUG(" TRTDetectorElement id: " << id << " with center = " << center);
 
         // add element to respective AlignModule
 
         // add to the Inner Detector structure
         if(msgLvl(MSG::DEBUG)) {
            if (m_idHelper->is_barrel(id))
               msg(MSG::DEBUG)<<"... TRT barrel element"<<endmsg;
            else
               msg(MSG::DEBUG)<<"... TRT endcap element"<<endmsg;
         }
         trt->addDetElement(Trk::AlignModule::TRT,element,transform);
 
         // and fill the corresponding map
         (*trtIdHashMap)[idHash] = trt;
      }
      else
         ATH_MSG_DEBUG("No TRTDetectorElement with id: "<<id);
   }
 
   // add created module to the geometry
   m_alignModuleListPtr->push_back(trt);
 
   ATH_MSG_DEBUG("TRT L0 module successfully added to the list");
}

buildL1Barrel()

void InDet::TRTGeometryManagerTool::buildL1Barrel ( )

private

builds geometry for L1 TRT Barrel alignment Creates single module for barrel

Definition at line 409 of file TRTGeometryManagerTool.cxx.

{
   // ===========================
   // BARREL
   // ===========================
   ATH_MSG_INFO("Preparing the TRT Barrelgeometry for L1: whole barrel");
 
   // ========================================
   // get all modules for barrel Level 1 alignment
   // Level 1 is barrel as a single AlignModule
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // get maximum number of elements from the helper
   unsigned int maxHash = 0;
   if(m_alignLevelEndcaps==3)
      maxHash = m_idHelper->straw_hash_max();
   else
      maxHash = m_idHelper->straw_layer_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<maxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)maxHash,nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
   // for L1 alignment the alignment frame is equal to the global frame
   // and since the TRT detector element positions are also stored in
   // the global frame in DB, transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // BARREL
   // use the 'real' identifier of the TRT barrel that can be
   // used by idHelper->barrel_ec()
   // TRT barrel has two sides so both -1 and +1 are ok
   // but -1 is used in the DB so we use it here as well
   Trk::AlignModule * trtBarrel = new Trk::AlignModule(this);
   trtBarrel->setIdHash(getNextIDHash());
   trtBarrel->setName("TRT Barrel");
   trtBarrel->setIdentifier(m_idHelper->barrel_ec_id(-1));
 
   if(!moduleSelected(trtBarrel)) {
      ATH_MSG_DEBUG("Module "<<trtBarrel->name()<<" NOT selected");
      delete trtBarrel;
   }
   else {
      // loop over detector elements in the barrel and add them to respective alignModules
      for (unsigned int iLayer = 0; iLayer < m_trtDetManager->getNumerology()->getNBarrelRings(); iLayer++) // barrel layers
         for (unsigned int iPhiModule = 0; iPhiModule < m_trtDetManager->getNumerology()->getNBarrelPhi(); iPhiModule++) // modules in phi
            for (unsigned int iStrawLayer = 0; iStrawLayer < m_trtDetManager->getNumerology()->getNBarrelLayers(iLayer); iStrawLayer++) // straw layers inside one phi module
               for (int iSide = 0; iSide < 2; iSide++) { // side of the straw
                  const TRT_BarrelElement * element = m_trtDetManager->getBarrelElement(iSide, iLayer, iPhiModule, iStrawLayer);
                  if (element) {
                     // get element location for debugging
//                     HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
//                     ATH_MSG_DEBUG(" TRTDetectorElement id: " << element->identify() << " with center = " << center);
 
                     // add to the barrel
                     trtBarrel->addDetElement(Trk::AlignModule::TRT,element,transform);
                     // and fill the corresponding map
                     (*trtIdHashMap)[element->identifyHash()] = trtBarrel;
                  }
                  else
                     ATH_MSG_DEBUG("No TRT_BarrelElement with side-layer-phi-strawLayer:" <<iSide<<"-"<<iLayer<<"-"<<iPhiModule<<"-"<<iStrawLayer);
               }
 
      // add
      m_alignModuleListPtr->push_back(trtBarrel);
   }
}

buildL1Endcaps()

void InDet::TRTGeometryManagerTool::buildL1Endcaps ( )

private

builds geometry for L1 TRT Endcap alignment Creates 2 modules, one for each endcap

Definition at line 484 of file TRTGeometryManagerTool.cxx.

{
   // ===========================
   // ENDCAPs
   // ===========================
   ATH_MSG_INFO("Preparing the TRT Endcap geometry for L1: 2 endcaps");
 
   // ========================================
   // get all modules for Endcap Level 1 alignment
   // Level 1 is endcap AlignModules
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // get maximum number of elements from the helper
   unsigned int maxHash = 0;
   if(m_alignLevelBarrel==3)
      maxHash = m_idHelper->straw_hash_max();
   else
      maxHash = m_idHelper->straw_layer_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<maxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)maxHash,nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
   // for L1 alignment the alignment frame is equal to the global frame
   // and since the TRT detector element positions are also stored in
   // the global frame in DB, transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // ENDCAPs
   // identifiers for endcaps are set in analogy to barrel
   // above to -2 and +2 so they can be used by
   // idHelper->barrel_ec()
   Trk::AlignModule * trtEndCap[2];
   for (int iSide=0;iSide<2;iSide++) {
      trtEndCap[iSide] = new Trk::AlignModule(this);
      trtEndCap[iSide] -> setIdHash(getNextIDHash());
 
      if(iSide)
         trtEndCap[iSide]->setName("TRT EndCap C");
      else
         trtEndCap[iSide]->setName("TRT EndCap A");
 
      int idSide = iSide ? 2 : -2;
      trtEndCap[iSide]->setIdentifier(m_idHelper->barrel_ec_id(idSide));
 
      if( !moduleSelected(trtEndCap[iSide]) || (m_endcap!=idSide && m_endcap!=0 ) ) {
         ATH_MSG_DEBUG("Module "<<trtEndCap[iSide]->name()<<" NOT selected");
         delete trtEndCap[iSide];
      }
      else {
         for (unsigned int iWheel = 0; iWheel < m_trtDetManager->getNumerology()->getNEndcapWheels(); iWheel++)
            for (unsigned int iStrawLayer = 0; iStrawLayer < m_trtDetManager->getNumerology()->getNEndcapLayers(iWheel); iStrawLayer++)
               for(unsigned int iPhi = 0; iPhi < m_trtDetManager->getNumerology()->getNEndcapPhi(); iPhi++) {
                  const TRT_EndcapElement * element = m_trtDetManager->getEndcapElement(iSide, iWheel, iStrawLayer, iPhi);
                  if (element) {
                     // get element location for debugging
//                     HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
//                     ATH_MSG_DEBUG(" TRTDetectorElement id: " << id << " with center = " << center);
 
                     // add element to respective AlignModule
                     trtEndCap[iSide]->addDetElement(Trk::AlignModule::TRT,element,transform);
                     // and fill the corresponding map
                     (*trtIdHashMap)[element->identifyHash()] = trtEndCap[iSide];
                  }
                  else
                     ATH_MSG_DEBUG("No TRT_EndcapElement with side-wheel-strawLayer-phi:" <<iSide<<"-"<<iWheel<<"-"<<iStrawLayer<<"-"<<iPhi);
               }
 
         // for endcap we move the CoG
         Amg::Translation3D translation(trtEndCap[iSide]->centerOfGravity());
         Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
         trtEndCap[iSide]->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
 
         // add
         m_alignModuleListPtr->push_back(trtEndCap[iSide]);
      }
   }
}

buildL2Barrel()

void InDet::TRTGeometryManagerTool::buildL2Barrel ( )

private

builds geometry for L2 TRT Barrel alignment Creates 96 modules - 3 layers with 32 modules per layer

Definition at line 570 of file TRTGeometryManagerTool.cxx.

{
   // ===========================
   // BARREL
   // ===========================
   ATH_MSG_INFO("Preparing the TRT Barrel geometry for L2: 3 layers with 32 modules each (96 modules)");
 
   // ========================================
   // get all modules for TRT Barrel Level 2 alignment
   // Level 2 in the barrel is 3 layers with 32 phi modules each
   // i.e. 96 AlignModules in total
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // TRT detector element positions are stored in the DB in global frame
   // so the DB frame to global frame transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // get maximum number of elements from the helper
   unsigned int maxHash = 0;
   if(m_alignLevelEndcaps==3)
      maxHash = m_idHelper->straw_hash_max();
   else
      maxHash = m_idHelper->straw_layer_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<maxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)maxHash,nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
   // for the alignment frame we would like to use something like
   // the 'center-of-gravity' (cog) of the module and the rotation
   // of the strawlayer planes. the main problem here is that the
   // planes are not entirely parallel. for the definition of the
   // alignment of the module this is all irrelevant, but we would
   // like to give the twist some meaning, like 'the rotation
   // around the most central straw': since strawlayer planes are
   // even and straws as well, there is no such thing as a really
   // central straw. we'll do something that brings us as close as
   // possible: for the central plane take strawlayer that is at
   // (N+1)/2, where N is the maximum strawlayer. (strawlayers
   // start counting at 0). The reason to pick the outer of the two
   // central strawlayers is that there are more straws per layer
   // in the outer layers and that the outer layers are also
   // further apart, so this is closer to the 'real' cog.
 
   for (unsigned int iPhiModule = 0; iPhiModule < m_trtDetManager->getNumerology()->getNBarrelPhi(); iPhiModule++) { // modules in phi
      for (unsigned int iLayer = 0; iLayer < m_trtDetManager->getNumerology()->getNBarrelRings(); iLayer++) {// barrel layers
         // create the AlignModule
         Trk::AlignModule * mod = new Trk::AlignModule(this);
         mod->setIdHash(getNextIDHash());
         mod->setIdentifier(m_idHelper->module_id(-1, iPhiModule, iLayer));
 
         std::stringstream name;
         name<<"TRT/Barrel/Module_-1_"<<iPhiModule<<"_"<<iLayer;
         mod->setName(name.str());
 
         if(!moduleSelected(mod) || !moduleSelectedBarrel(iLayer,iPhiModule)) {
            ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
            delete mod;
            continue;
         }
 
         ATH_MSG_DEBUG("Building module "<<mod->name());
 
         // create the collection of central layers for cog and transform calculation
         std::vector<const Trk::TrkDetElementBase *> centralLayers;
 
         int maxstrawlayer = m_trtDetManager->getNumerology()->getNBarrelLayers(iLayer);
         ATH_MSG_DEBUG("maxstrawlayer "<<maxstrawlayer);
         for (int iStrawLayer = 0; iStrawLayer < maxstrawlayer; iStrawLayer++) { // straw layers inside one phi module
            ATH_MSG_DEBUG("iStrawLayer "<<iStrawLayer);
            for (int iSide = 0; iSide < 2; iSide++) { // side of the straw
               ATH_MSG_DEBUG("iSide "<<iSide);
               const TRT_BarrelElement * element = m_trtDetManager->getBarrelElement(iSide, iLayer, iPhiModule, iStrawLayer);
               if (element) {
                  // get element location for debugging
//                  HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
//                  ATH_MSG_DEBUG(" TRTDetectorElement id: " << element->identify() << " with center = " << center);
 
                  // update maximum straw layer
                  maxstrawlayer = std::max(maxstrawlayer,m_idHelper->straw_layer(element->identify()));
 
                  // add element to respective AlignModule
                  mod->addDetElement(Trk::AlignModule::TRT,element,transform);
 
                  // and fill the corresponding map
                  (*trtIdHashMap)[element->identifyHash()] = mod;
 
                  // add central layers to centralLayers collection
                  // maxstrawlayer is probably odd (if there is an even number
                  // of strawlayers), but I could be wrong, so let's make this always work:
                  // if maxstrawlayer==23, we want to pick strawlayer 12
                  // if maxstrawlayer==24, we still want to pick strawlayer 12
                  if( iStrawLayer == maxstrawlayer/2 ) {
                     //Amg::Vector3D center = element->transform() * Amg::Vector3D(0,0,0);
                     centralLayers.push_back(element);
                  }
               }
               else
                  ATH_MSG_DEBUG("No TRT_BarrelElement with side-layer-phi-strawLayer:" <<iSide<<"-"<<iLayer<<"-"<<iPhiModule<<"-"<<iStrawLayer);
            }
         }
 
         // This is how the frames line up for a bec=-1 module at global
         // position (0,R,-360):
         //
         // global       straw   det-element
         //     x          +x       -x
         //     y          +y       +z
         //     z          +z       +y
         //
         // For other bec=-1 modules you just need to rotate around
         // global y. For bec=+1 modules it looks like this:
         //
         // global       straw   det-element
         //     x          -x       -x
         //     y          +y       +z
         //     z          -z       +y
         // So that just involves one rotation around local-y for the strawreference frame.
         //
         std::vector<std::vector<const Trk::TrkDetElementBase *> *> centralLayersVec;
         centralLayersVec.push_back(&centralLayers);
 
         Amg::Translation3D translation(Trk::AlignModule::centerOfGravity(centralLayersVec));
         // for the rotation, we choose the element on the '-1' side.
         Amg::RotationMatrix3D rotation = ( m_idHelper->barrel_ec(centralLayers.front()->identify()) == -1 ) ?
            centralLayers.front()->transform().rotation() : centralLayers.back()->transform().rotation();
         Amg::Transform3D localToGlobal = translation * rotation;
 
 
 
         ATH_MSG_DEBUG("centerOfGravity "<< mod->centerOfGravity(centralLayersVec));
         ATH_MSG_DEBUG("Prepared local to global transform : "<<std::setprecision(12)<<localToGlobal.translation()<<" "<<localToGlobal.rotation());
 
         mod->setGlobalFrameToAlignFrameTransform(localToGlobal.inverse());
 
         // add AlignModule to the geometry
         m_alignModuleListPtr->push_back(mod);
      }
   }
 
   }

buildL2Endcaps()

void InDet::TRTGeometryManagerTool::buildL2Endcaps ( )

private

builds geometry for L2 TRT Endcapalignment Creates 80 modules - 40 rings per endcap

Definition at line 719 of file TRTGeometryManagerTool.cxx.

{
   // ===========================
   // ENDCAPs
   // ===========================
   ATH_MSG_INFO("Preparing the TRT Endcap geometry for L2: 2 x 40 rings");
 
   // ========================================
   // get all modules for TRT Endcap Level 2 alignment
   // Level 2 in the endcaps is 40 rings per endcap
   // i.e. 40 AlignModules in total
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // TRT detector element positions are stored in the DB in global frame
   // so the DB frame to global frame transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // get maximum number of elements from the helper
   unsigned int maxHash = 0;
   if(m_alignLevelBarrel==3)
      maxHash = m_idHelper->straw_hash_max();
   else
      maxHash = m_idHelper->straw_layer_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<maxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)maxHash,nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
   for (int iSide = 0; iSide < 2; iSide++) {
      int idSide = iSide ? 2 : -2;
      ATH_MSG_DEBUG("Side : "<<iSide<<" , idSide : "<<idSide);
      if( m_endcap!=idSide && m_endcap!=0 ) {
         ATH_MSG_DEBUG("Endcap "<<idSide<<" NOT aligned");
         continue;
      }
 
      int iRing(-1);
      std::unique_ptr<Trk::AlignModule> mod;
      for (unsigned int iWheel = 0; iWheel < m_trtDetManager->getNumerology()->getNEndcapWheels(); iWheel++) {
         ATH_MSG_DEBUG("Wheel : "<<iWheel);
         for (unsigned int iStrawLayer = 0; iStrawLayer < m_trtDetManager->getNumerology()->getNEndcapLayers(iWheel); iStrawLayer++) {
            ATH_MSG_DEBUG("StrawLayer : "<<iStrawLayer);
 
            if(iStrawLayer%4==0) {
               if(iRing >= 0 && mod) {
                  // before creating module for new ring we set the alignment
                  // frame for the previous one and add it to the list of modules
                  ATH_MSG_DEBUG("Setting frame for ring : "<<iRing);
 
                  // for endcap we move the CoG with no additional rotation
                  Amg::Translation3D translation(mod->centerOfGravity());
            
                  Amg::Transform3D localToGlobal = translation * Amg::RotationMatrix3D::Identity();
 
                  mod->setGlobalFrameToAlignFrameTransform(localToGlobal.inverse());
 
                  // add AlignModule to the geometry
                  m_alignModuleListPtr->push_back(mod.release());
               }
 
               // new ring
               iRing++;
 
               ATH_MSG_DEBUG("Ring : "<<iRing);
               // create the AlignModule
               mod.reset(new Trk::AlignModule(this));
               mod->setIdHash(getNextIDHash());
               // Identifier for a ring is the Identifier for the first Straw layer in that ring
               // one ring has 4 straw layers
               int strawLayerForRing = (iStrawLayer/4) * 4;
               mod->setIdentifier(m_idHelper->layer_id((iSide ? 2:-2),0,iWheel,strawLayerForRing));
 
               std::stringstream name;
               name<<"TRT/Endcap/Module_"<<(iSide ? 2:-2)<<"_"<<iRing;
               mod->setName(name.str());
 
               if(!moduleSelected(mod.get())) {
                  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
                  mod.reset();
                  continue;
               }
 
               ATH_MSG_DEBUG("Building module "<<mod->name());
            }
            if (!mod) continue;
            for(unsigned int iPhi = 0; iPhi < m_trtDetManager->getNumerology()->getNEndcapPhi(); iPhi++) {
               const TRT_EndcapElement * element = m_trtDetManager->getEndcapElement(iSide, iWheel, iStrawLayer, iPhi);
               if (element) {
                  // get element location for debugging
//                  HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
//                  ATH_MSG_DEBUG(" TRTDetectorElement id: " << id << " with center = " << center);
 
                  // add element to respective AlignModule
                  mod->addDetElement(Trk::AlignModule::TRT,element,transform);
 
                  // and fill the corresponding map
                  (*trtIdHashMap)[element->identifyHash()] = mod.get();
               }
               else
                  ATH_MSG_DEBUG("No TRT_EndcapElement with side-wheel-strawLayer-phi:" <<iSide<<"-"<<iWheel<<"-"<<iStrawLayer<<"-"<<iPhi);
            }
         }
      }
 
      if (mod) {
        // for the last ring we have to explicitly set the alignment
        // frame add it to the list of modules here, at the end of
        // loop over wheels
        ATH_MSG_DEBUG("Setting frame for ring : "<<iRing);
        
        // for endcap we move the CoG with no additional rotation
        Amg::Translation3D translation(mod->centerOfGravity());      
        Amg::Transform3D localToGlobal = translation * Amg::RotationMatrix3D::Identity();
  
        mod->setGlobalFrameToAlignFrameTransform(localToGlobal.inverse());
  
        // add AlignModule to the geometry
        m_alignModuleListPtr->push_back(mod.release());
      }
   }
}

buildL2EndcapsOLD()

void InDet::TRTGeometryManagerTool::buildL2EndcapsOLD ( )

private

builds geometry for L2 TRT Endcapalignment Creates 28 modules - 14 wheels per endcap

Definition at line 848 of file TRTGeometryManagerTool.cxx.

{
   // !!!!   WARNING   !!!!
   // !!!!   This method sets up the TRT Endcap level 2 geometry which
   // !!!!   is now obsolete. The method buildL2Endcap() should be used
   // !!!!   instead !!!
 
   // ===========================
   // ENDCAPs
   // ===========================
   ATH_MSG_INFO("Preparing the _OLD_ TRT Endcap geometry for L2: 2 x 14 wheels");
 
   // ========================================
   // get all modules for TRT Endcap Level 2 alignment
   // Level 2 in the endcaps is 14 wheels per endcap
   // i.e. 28 AlignModules in total
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // TRT detector element positions are stored in the DB in global frame
   // so the DB frame to global frame transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // get maximum number of elements from the helper
   unsigned int maxHash = 0;
   if(m_alignLevelBarrel==3)
      maxHash = m_idHelper->straw_hash_max();
   else
      maxHash = m_idHelper->straw_layer_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<maxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)maxHash,nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
   for (int iSide = 0; iSide < 2; iSide++) {
      int idSide = iSide ? 2 : -2;
      ATH_MSG_DEBUG("Side : "<<iSide<<" , idSide : "<<idSide);
      if( m_endcap!=idSide && m_endcap!=0 ) {
         ATH_MSG_DEBUG("Endcap "<<idSide<<" NOT aligned");
         continue;
      }
 
      for (unsigned int iWheel = 0; iWheel < m_trtDetManager->getNumerology()->getNEndcapWheels(); iWheel++) {
         // create the AlignModule
         Trk::AlignModule * mod = new Trk::AlignModule(this);
         mod->setIdHash(getNextIDHash());
         mod->setIdentifier(m_idHelper->module_id(idSide,0,iWheel));
 
         std::stringstream name;
         name<<"TRT/Endcap/Module_"<<idSide<<"_"<<iWheel;
         mod->setName(name.str());
 
         if(!moduleSelected(mod)) {
            ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
            delete mod;
            continue;
         }
 
         ATH_MSG_DEBUG("Building module "<<mod->name());
 
         for (unsigned int iStrawLayer = 0; iStrawLayer < m_trtDetManager->getNumerology()->getNEndcapLayers(iWheel); iStrawLayer++) {
            for(unsigned int iPhi = 0; iPhi < m_trtDetManager->getNumerology()->getNEndcapPhi(); iPhi++) {
               const TRT_EndcapElement * element = m_trtDetManager->getEndcapElement(iSide, iWheel, iStrawLayer, iPhi);
               if (element) {
                  // get element location for debugging
//                  HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
//                  ATH_MSG_DEBUG(" TRTDetectorElement id: " << id << " with center = " << center);
 
                  // add element to respective AlignModule
                  mod->addDetElement(Trk::AlignModule::TRT,element,transform);
 
                  // and fill the corresponding map
                  (*trtIdHashMap)[element->identifyHash()] = mod;
               }
               else
                  ATH_MSG_DEBUG("No TRT_EndcapElement with side-wheel-strawLayer-phi:" <<iSide<<"-"<<iWheel<<"-"<<iStrawLayer<<"-"<<iPhi);
            }
         }
 
         // for endcap we move the CoG with no additional rotation
         Amg::Translation3D translation(mod->centerOfGravity());
         Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
 
         mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
 
         // add AlignModule to the geometry
         m_alignModuleListPtr->push_back(mod);
      }
   }
 
   }

buildL3Barrel()

void InDet::TRTGeometryManagerTool::buildL3Barrel ( )

private

builds geometry for L3 TRT Barrel alignment Creates many modules - one module per straw

Definition at line 946 of file TRTGeometryManagerTool.cxx.

{
   // ===========================
   // BARREL Level 3
   // ===========================
   ATH_MSG_INFO("Preparing the TRT Barrel geometry for L3");
 
   // ========================================
   // get all modules for TRT Barrel Level 3 alignment
   // Level 3 is one AlignModule per straw
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // TRT detector element positions are stored in the DB in global frame
   // so the DB frame to global frame transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // get maximum number of elements from the helper
   unsigned int maxHash = m_idHelper->straw_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<maxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)maxHash,nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
   for (unsigned int iPhiModule = 0; iPhiModule < m_trtDetManager->getNumerology()->getNBarrelPhi(); iPhiModule++) { // modules in phi
      ATH_MSG_DEBUG("iPhiModule "<<iPhiModule);
      for (unsigned int iLayer = 0; iLayer < m_trtDetManager->getNumerology()->getNBarrelRings(); iLayer++) {// barrel layers
         ATH_MSG_DEBUG("iLayer "<<iLayer);
 
         if(!moduleSelectedBarrel(iLayer,iPhiModule)) {
            ATH_MSG_INFO("Skipping barrel layer "<<iLayer<<"  phi sector "<<iPhiModule);
            continue;
         }
 
         for (unsigned int iStrawLayer = 0; iStrawLayer < m_trtDetManager->getNumerology()->getNBarrelLayers(iLayer); iStrawLayer++) { // straw layers inside one phi module
            ATH_MSG_DEBUG("iStrawLayer "<<iStrawLayer);
            for (int iSide = 0; iSide < 2; iSide++) { // side of the straw
               ATH_MSG_DEBUG("iSide "<<iSide);
               int idSide = iSide ? 1 : -1;
 
               const TRT_BarrelElement * element = m_trtDetManager->getBarrelElement(iSide, iLayer, iPhiModule, iStrawLayer);
               if(element) {
                  for (unsigned int iStraw = 0; iStraw< element->nStraws(); iStraw++)
                  {
                     ATH_MSG_DEBUG("iStraw "<<iStraw);
 
                     // create the AlignModule
                     Trk::AlignModule * mod = new Trk::AlignModule(this);
                     mod->setIdHash(getNextIDHash());
                     Identifier moduleId = m_idHelper->straw_id(idSide,iPhiModule,iLayer,iStrawLayer,iStraw);
                     mod->setIdentifier(moduleId);
 
                     std::stringstream name;
                     name<<"TRT/Barrel/Straw_"<<idSide<<"_"<<iPhiModule<<"_"<<iLayer<<"_"<<iStrawLayer<<"_"<<iStraw;
                     mod->setName(name.str());
 
                    if(!moduleSelected(mod)) {
                       ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
                       delete mod;
                       mod=nullptr;
                       continue;
                    }
 
                    ATH_MSG_DEBUG("Building module "<<mod->name());
 
                    // add element to respective AlignModule
                    mod->addDetElement(Trk::AlignModule::TRT, element, transform, moduleId);
 
                    // and fill the corresponding map
                    (*trtIdHashMap)[m_idHelper->straw_hash(moduleId)] = mod;
 
                    // set the align frame transform
                    Amg::Transform3D localToGlobal = element->strawTransform(iStraw);
                    Amg::RotationMatrix3D lgrot = localToGlobal.rotation();
                    ATH_MSG_DEBUG("Prepared local to global transform : "<<std::setprecision(12)
                                  <<localToGlobal.translation()<<" "<<localToGlobal.rotation());
                    ATH_MSG_DEBUG("  rotation = (" << lgrot(2,1) << "," << lgrot(0,2) << "," << lgrot(1,0) << ")");
 
                    mod->setGlobalFrameToAlignFrameTransform(localToGlobal.inverse());
 
                    Amg::Transform3D aatr = element->transform();
                    Amg::RotationMatrix3D aarot = aatr.rotation();
 
                    ATH_MSG_DEBUG("Straw layer transform :              "<<std::setprecision(12)
                          <<aatr.translation()<<" "<<aatr.rotation());
                    ATH_MSG_DEBUG("  rotation = (" << aarot(2,1) << "," << aarot(0,2) << "," << aarot(1,0) << ")");
 
                    // add AlignModule to the geometry
                    m_alignModuleListPtr->push_back(mod);
                  }
               }
               else
                  ATH_MSG_DEBUG("No TRT_BarrelElement with side-layer-phi-strawLayer:" <<iSide<<"-"<<iLayer<<"-"<<iPhiModule<<"-"<<iStrawLayer);
            }
         }
      }
   }
 
   }

buildL3Endcap()

void InDet::TRTGeometryManagerTool::buildL3Endcap ( )

private

builds geometry for L3 TRT Endcap alignment Creates many modules - one module per straw

Definition at line 1051 of file TRTGeometryManagerTool.cxx.

{
   // ===========================
   // ENDCAP Level 3
   // ===========================
   ATH_MSG_INFO("Preparing the TRT Endcap geometry for L3");
 
   // ========================================
   // get all modules for TRT Endcap Level 3 alignment
   // Level 3 is one AlignModule per straw
 
   // define modules
   // we'll stick with the base AlignModule here since we don't want to
   // do shifting for the moment but we might need to extend AlignModule
   // in the future
 
   // TRT detector element positions are stored in the DB in the local frame
   // so the DB frame to local frame transform is identity
   Amg::Transform3D transform = Amg::Transform3D::Identity();
 
   // get maximum number of elements from the helper
   unsigned int maxHash = m_idHelper->straw_hash_max();
   ATH_MSG_DEBUG("maxHash for the TRT "<<maxHash);
 
   if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT))
     m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT) = new Trk::AlignModuleList((size_t)maxHash,nullptr);
   Trk::AlignModuleList * trtIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::TRT);
 
 
   for (int iSide = 0; iSide < 2; iSide++) { // endcap
     int idSide = iSide ? 2 : -2;
     ATH_MSG_DEBUG("Side : "<<iSide<<" , idSide : "<<idSide);
     if( m_endcap!=idSide && m_endcap!=0 ) {
        ATH_MSG_DEBUG("Endcap "<<idSide<<" NOT aligned");
        continue;
     }
 
     // wheels in Z
     for (unsigned int iWheel = 0; iWheel < m_trtDetManager->getNumerology()->getNEndcapWheels(); iWheel++) {
       ATH_MSG_DEBUG("iWheel "<<iWheel);
 
       // sectors in Phi
       for (unsigned int iPhiModule = 0; iPhiModule < m_trtDetManager->getNumerology()->getNEndcapPhi(); iPhiModule++) {
 
         ATH_MSG_DEBUG("iPhiModule "<<iPhiModule);
 
         if(!moduleSelectedEndcap(iWheel,iPhiModule)) {
           ATH_MSG_INFO("Skipping endcap wheel "<<iWheel<<"  phi sector "<<iPhiModule);
           continue;
         }
 
         // endcap straw layers in z (with in a wheel)
         for (unsigned int iStrawLayer = 0; iStrawLayer < m_trtDetManager->getNumerology()->getNEndcapLayers(iWheel); iStrawLayer++) {
 
           const TRT_EndcapElement * element = m_trtDetManager->getEndcapElement(iSide, iWheel, iStrawLayer, iPhiModule);
 
           if(element) {
              for (unsigned int iStraw = 0; iStraw< element->nStraws(); iStraw++) {
                 ATH_MSG_DEBUG("iStraw "<<iStraw);
 
                 // create the AlignModule
                 Trk::AlignModule * mod = new Trk::AlignModule(this);
                 mod->setIdHash(getNextIDHash());
                 Identifier moduleId = m_idHelper->straw_id(idSide,iPhiModule,iWheel,iStrawLayer,iStraw);
                 mod->setIdentifier(moduleId);
 
                 std::stringstream name;
                 name<<"TRT/Endcap/Straw_"<<idSide<<"_"<<iPhiModule<<"_"<<iWheel<<"_"<<iStrawLayer<<"_"<<iStraw;
                 mod->setName(name.str());
 
                 if(!moduleSelected(mod)) {
                    ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
                    delete mod;
                    mod=nullptr;
                    continue;
                 }
 
                 ATH_MSG_DEBUG("Building module "<<mod->name());
 
                 // add element to respective AlignModule
                 mod->addDetElement(Trk::AlignModule::TRT, element, transform, moduleId);
 
                 // and fill the corresponding map
                 (*trtIdHashMap)[m_idHelper->straw_hash(moduleId)] = mod;
 
                 // set the align frame transform
                 Amg::Transform3D localToGlobal = element->strawTransform(iStraw);
                 Amg::RotationMatrix3D lgrot = localToGlobal.rotation();
 
                 ATH_MSG_DEBUG("Prepared local to global transform : "<<std::setprecision(12)
                               <<localToGlobal.translation()<<" "<<localToGlobal.rotation());
 
                 ATH_MSG_DEBUG("  rotation = (" << lgrot(2,1) << "," << lgrot(0,2) << "," << lgrot(1,0) << ")");
 
                 mod->setGlobalFrameToAlignFrameTransform(localToGlobal.inverse());
 
                 Amg::Transform3D aatr = element->transform();
                 Amg::RotationMatrix3D aarot = aatr.rotation();
 
                 ATH_MSG_DEBUG("Straw layer transform :              "<<std::setprecision(12)
                               <<aatr.translation()<<" "<<aatr.rotation());
                 ATH_MSG_DEBUG("  rotation = (" << aarot(2,1) << "," << aarot(0,2) << "," << aarot(1,0) << ")");
 
                 // add AlignModule to the geometry
                 m_alignModuleListPtr->push_back(mod);
              }// Straws
           }// if the EndcapElement Exists
           else
             ATH_MSG_DEBUG("No TRT_EndcapElement with endcap-wheel-strawLayer-phi:" <<iSide<<"-"<<iWheel<<"-"<<iStrawLayer<<"-"<<iPhiModule);
         }// StrawLayer
       }// PhiModule
     }// Wheel
   }// Endcap
 
   }

checkAlignLevel()

bool InDet::TRTGeometryManagerTool::checkAlignLevel ( )

virtual

check whether the alignment level is correct

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 182 of file TRTGeometryManagerTool.cxx.

{
   if(m_alignLevel == 0) {
      ATH_MSG_INFO("Alignment level for TRT is "<<m_alignLevel);
      return true;
   }
 
   if(m_alignLevelBarrel == -1) {
      if(m_alignLevel == 200)
         m_alignLevelBarrel = 2;
      else
         m_alignLevelBarrel = m_alignLevel;
   }
 
   if(m_alignLevelEndcaps == -1)
      m_alignLevelEndcaps = m_alignLevel;
 
   bool ok = checkAlignLevelBarrel() && checkAlignLevelEndcaps();
 
   // if the alignment level for barrel and EC is set correctly
   // reset the global alignment level
   if(ok)
      m_alignLevel = -1;
 
   return ok;
}

checkAlignLevelBarrel()

bool InDet::TRTGeometryManagerTool::checkAlignLevelBarrel ( )

virtual

check whether the Barrel alignment level is correct

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 210 of file TRTGeometryManagerTool.cxx.

{
   if(!m_alignBarrel)
      return true;
 
   // check whether geometry level is allowed
   switch(m_alignLevelBarrel) {
      case 1: case 2: case 3:
         ATH_MSG_INFO("Alignment level for TRT Barrel is "<<m_alignLevelBarrel);
         return true;
      default:
         msg(MSG::FATAL)<<"Alignment level "<<m_alignLevelBarrel<<" does not exist for TRT Barrel"<<endmsg;
         return false;
   }
}

checkAlignLevelEndcaps()

bool InDet::TRTGeometryManagerTool::checkAlignLevelEndcaps ( )

virtual

check whether the Endcap alignment level is correct

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 227 of file TRTGeometryManagerTool.cxx.

{
   if(!m_alignEndcaps)
      return true;
 
   // check whether geometry level is allowed
   switch(m_alignLevelEndcaps) {
      case 1: case 2: case 200: case 3:
         ATH_MSG_INFO("Alignment level for TRT Endcaps is "<<m_alignLevelEndcaps);
         return true;
      default:
         msg(MSG::FATAL)<<"Alignment level "<<m_alignLevelEndcaps<<" does not exist for TRT Endcaps"<<endmsg;
         return false;
   }
}

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());
 
  }

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; }

dumpGeometry()

void InDet::TRTGeometryManagerTool::dumpGeometry ( )

private

print basic geometry info to screen

Definition at line 1376 of file TRTGeometryManagerTool.cxx.

{
   ATH_MSG_INFO("---------------------------------------------------");
   ATH_MSG_INFO("Summary of the alignment geometry");
   ATH_MSG_INFO("Number of alignable objects: "<<m_alignModuleList.size());
   ATH_MSG_INFO("Total number of aligment DoFs: "<<m_alignModuleTool->nAlignParameters());
 
   for(unsigned int i=0;i<m_alignModuleList.size();i++) {
      const Trk::AlignModule * module = m_alignModuleList.at(i);
      //ATH_MSG_DEBUG(i<<". "<< module->name());
      //ATH_MSG_DEBUG("   - identifier: "<<module->identify());
      //ATH_MSG_DEBUG("   - has "<<module->detElementCollection(Trk::AlignModule::TRT)->size()<<" TRT modules");
 
      Amg::Transform3D localtoglobal = (module->globalFrameToAlignFrame()).inverse();
      ATH_MSG_DEBUG("   - local to global : "<<std::setprecision(12)<<localtoglobal.translation()<<" "<<localtoglobal.rotation());
 
      DataVector<Trk::AlignPar> * pars = m_alignModuleTool->getAlignPars(module);
      int npars = pars->size();
      ATH_MSG_DEBUG("   - number of active transform parameters: "<<npars);
      for(int j=0;j<npars;j++)
         ATH_MSG_DEBUG("      * par "<<j<<" \'"<<(*pars)[j]->dumpType()<<"\' : sigma = "<<(*pars)[j]->sigma()<<" , softCut = "<<(*pars)[j]->softCut());
   }
   ATH_MSG_INFO("---------------------------------------------------");
}

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

finalize()

StatusCode InDet::TRTGeometryManagerTool::finalize ( )

virtual

Implements Trk::IGeometryManagerTool.

Definition at line 167 of file TRTGeometryManagerTool.cxx.

{
   ATH_MSG_DEBUG("finalize() of TRTGeometryManagerTool");
 
   return StatusCode::SUCCESS;
}

getNextIDHash()

virtual int Trk::IGeometryManagerTool::getNextIDHash ( )

inlinevirtualinherited

get next free IDHash usable for new AlignModule

Definition at line 65 of file IGeometryManagerTool.h.

{ return m_hashCounter++; };

initialize()

StatusCode InDet::TRTGeometryManagerTool::initialize ( )

virtual

Implements Trk::IGeometryManagerTool.

Definition at line 124 of file TRTGeometryManagerTool.cxx.

{
   ATH_MSG_DEBUG("initialize() of TRTGeometryManagerTool");
 
   // retrieve AlignModuleTool
   if ( m_alignModuleTool.retrieve().isFailure() ) {
      msg(MSG::FATAL)<<"Could not get " << m_alignModuleTool << endmsg;
      return StatusCode::FAILURE;
   }
   else
      ATH_MSG_INFO("Retrieved " << m_alignModuleTool);
 
   // retrieve TRT ID helper
   if ( detStore()->retrieve(m_idHelper, "TRT_ID").isFailure() ) {
      msg(MSG::FATAL) << " Cannot retrieve TRT ID Helper " << endmsg;
      return StatusCode::FAILURE;
   }
   else
      ATH_MSG_INFO("retrieved TRT ID Helper");
 
   // retrieve TRT detector manager
   if ( detStore()->retrieve(m_trtDetManager).isFailure() ) {
      msg(MSG::FATAL) << " Cannot retrieve TRT Detector Manager " << endmsg;
      return StatusCode::FAILURE;
   }
   else
      ATH_MSG_INFO("retrieved TRT Detector Manager");
 
   // dump module selection
   if(m_doModuleSelection && msgLvl(MSG::INFO)) {
      msg(MSG::INFO)<<"Creating geometry for selected "<<m_moduleSelection.size()<<" modules:"<<endmsg;
      for(unsigned int i=0;i<m_moduleSelection.size();i++)
         msg(MSG::INFO)<<"   "<<i<<".  "<<m_moduleSelection.at(i)<<endmsg;
   }
 
   // check the allowed geometry levels
   if(!checkAlignLevel())
      return StatusCode::FAILURE;
 
   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.

interfaceID()

const InterfaceID & IGeometryManagerTool::interfaceID ( )

inlinestaticinherited

Definition at line 121 of file IGeometryManagerTool.h.

                                                            {
  return   IID_TRKALIGNINTERFACES_IGeometryManagerTool;
}

moduleSelected() [1/2]

bool InDet::TRTGeometryManagerTool::moduleSelected ( Trk::AlignModule * mod )

private

check wheather module is selected for module pointer

Definition at line 1308 of file TRTGeometryManagerTool.cxx.

{
   return moduleSelected(mod->identify().get_compact());
}

moduleSelected() [2/2]

bool InDet::TRTGeometryManagerTool::moduleSelected ( unsigned long long modId )

private

check wheather module is selected for module id

Definition at line 1293 of file TRTGeometryManagerTool.cxx.

{
   if(!m_doModuleSelection)
      return true;
 
   int nsel = m_moduleSelection.size();
   for(int i=0;i<nsel;++i)
      if(m_moduleSelection.at(i) == id)
         return true;
 
   return false;
}

moduleSelectedBarrel()

bool InDet::TRTGeometryManagerTool::moduleSelectedBarrel ( unsigned int iLayer, unsigned int iPhi )

private

check wheather barrel module is selected by layer and phi

Definition at line 1315 of file TRTGeometryManagerTool.cxx.

{
   bool select = true;
   if(!m_barrelLayers.empty()) {
      select = false;
      for(unsigned int i=0;i<m_barrelLayers.size();i++) {
         if(m_barrelLayers[i] == (int)ilayer) {
            select = true;
            break;
         }
      }
   }
 
   if(!select)
      return false;
 
   if(!m_barrelPhiSectors.empty()) {
      select = false;
      for(unsigned int i=0;i<m_barrelPhiSectors.size();i++) {
         if(m_barrelPhiSectors[i] == (int)iphi) {
            select = true;
            break;
         }
      }
   }
 
   return select;
}

moduleSelectedEndcap()

bool InDet::TRTGeometryManagerTool::moduleSelectedEndcap ( unsigned int iwheel, unsigned int iphi )

private

check wheather we are aligning the different endcaps and sectors

Definition at line 1345 of file TRTGeometryManagerTool.cxx.

{
  bool select = true;
  if(!m_endcapWheels.empty()) {
    select = false;
    for(unsigned int i=0;i<m_endcapWheels.size();++i) {
      if(m_endcapWheels[i] == (int)iwheel) {
        select = true;
        break;
      }
    }
  }
 
  if(!select)
    return false;
 
  if(!m_endcapPhiSectors.empty()) {
    select = false;
    for(unsigned int i=0;i<m_endcapPhiSectors.size();i++) {
      if(m_endcapPhiSectors[i] == (int)iphi) {
        select = true;
        break;
      }
    }
  }
 
  return select;
}

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);
  }

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.

PrintModPositions()

void InDet::TRTGeometryManagerTool::PrintModPositions ( )

inlinevirtual

Implements Trk::IGeometryManagerTool.

Definition at line 58 of file TRTGeometryManagerTool.h.

{};

ReadGeometry()

int InDet::TRTGeometryManagerTool::ReadGeometry ( int solveLevel )

virtual

read the geometry Method is called from the main AlignAlg to build the geometry based on the requested level of alignment.

It calls the other build() methods based on that.

Implements Trk::IGeometryManagerTool.

Definition at line 244 of file TRTGeometryManagerTool.cxx.

{
   ATH_MSG_DEBUG("in ReadGeometry() solveLevel="<<solveLevel);
 
   // set pointers
   m_idHashToAlignModuleMapsPtr = &m_idHashToAlignModuleMaps;
   m_alignModuleListPtr = &m_alignModuleList;
 
   // build alignment geometry
   buildGeometry();
 
   // now set the alignment parameters
   // first prepare the parameter lists
   m_alignParList     = new DataVector< DataVector<Trk::AlignPar> >(SG::OWN_ELEMENTS);
   m_fullAlignParList = new DataVector< DataVector<Trk::AlignPar> >(SG::OWN_ELEMENTS);
   // loop over modules
   std::vector<Trk::AlignModule *>::const_iterator imod = m_alignModuleList.begin();
   std::vector<Trk::AlignModule *>::const_iterator imod_end = m_alignModuleList.end();
   for( ; imod!=imod_end; ++imod)
      addModuleParameters(*imod,m_fullAlignParList,m_alignParList);
 
   // set alignModuleList and hash table in the alignModuleTool
   m_alignModuleTool->setAlignModules(&m_alignModuleList, &m_idHashToAlignModuleMaps);
   ATH_MSG_DEBUG(" geometry set in m_alignModuleTool");
 
   // set alignPar lists in the alignModuleTool
   ATH_MSG_DEBUG(" alignParList = "<<m_alignParList);
   ATH_MSG_DEBUG(" fullAlignParList = "<<m_fullAlignParList);
   m_alignModuleTool->setAlignParLists(m_fullAlignParList, m_alignParList);
   ATH_MSG_DEBUG(" AlignParLists set in m_alignModuleTool");
 
   // dump summary about the geometry setup
   if (m_dumpGeometry)
      dumpGeometry();
 
   int nDoF= m_alignModuleTool->nAlignParameters();
   ATH_MSG_INFO("Total number of degrees of freedom: "<<nDoF);
 
   return nDoF;
}

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();
  }

setAlignLevel()

void InDet::TRTGeometryManagerTool::setAlignLevel ( int level )

virtual

set alignment level

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 175 of file TRTGeometryManagerTool.cxx.

{
   m_alignLevelBarrel = m_alignLevelEndcaps = -1;
   m_alignLevel = level;
}

setAlignLevelBarrel()

virtual void Trk::IGeometryManagerTool::setAlignLevelBarrel ( int level )

inlinevirtualinherited

set alignment level

Definition at line 89 of file IGeometryManagerTool.h.

{ m_alignLevelBarrel = level; }

setAlignLevelEndcaps()

virtual void Trk::IGeometryManagerTool::setAlignLevelEndcaps ( int level )

inlinevirtualinherited

set alignment level

Definition at line 98 of file IGeometryManagerTool.h.

{ m_alignLevelEndcaps = level; }

setFirstIDHash()

virtual void Trk::IGeometryManagerTool::setFirstIDHash ( unsigned int idhash )

inlinevirtualinherited

set first free IDHash usable for new AlignModule

Definition at line 68 of file IGeometryManagerTool.h.

{ m_hashCounter=idhash; };

setHashMaps()

virtual void Trk::IGeometryManagerTool::setHashMaps ( std::vector< Trk::AlignModuleList * > * maps )

inlinevirtualinherited

set pointer to vector of hashMaps to which the elements should be added to

Definition at line 74 of file IGeometryManagerTool.h.

{ m_idHashToAlignModuleMapsPtr=maps; }

setLogStream()

virtual void Trk::IGeometryManagerTool::setLogStream ( std::ostream * os )

inlinevirtualinherited

sets the output stream for the logfile

Reimplemented in InDet::InDetGeometryManagerTool, and InDet::SiGeometryManagerTool.

Definition at line 55 of file IGeometryManagerTool.h.

{ m_logStream = os; }

setModuleList()

virtual void Trk::IGeometryManagerTool::setModuleList ( Trk::AlignModuleList * list )

inlinevirtualinherited

set pointer to module list to which the modules should be added to

Definition at line 71 of file IGeometryManagerTool.h.

{ m_alignModuleListPtr=list; }

setNtuple()

virtual void Trk::IGeometryManagerTool::setNtuple ( TFile * )

inlinevirtualinherited

Sets output ntuple file, to be implemented if needed for detector-specific implementation.

Definition at line 52 of file IGeometryManagerTool.h.

{} 

setSigmas()

void InDet::TRTGeometryManagerTool::setSigmas ( Trk::AlignModule * mod, DataVector< Trk::AlignPar > * modPars )

private

sets sigmas for modules

Definition at line 1225 of file TRTGeometryManagerTool.cxx.

{
   for(unsigned int ipar=0;ipar<modPars->size();++ipar) {
      // barrel
      if(m_idHelper->is_barrel(module->identify()))
         switch(modPars->at(ipar)->paramType()) {
            case Trk::AlignModule::TransX:
               modPars->at(ipar)->setSigma(m_sigmaBarrelX);
               modPars->at(ipar)->setSoftCut(m_softcutBarrelX);
               break;
            case Trk::AlignModule::TransY:
               modPars->at(ipar)->setSigma(m_sigmaBarrelY);
               modPars->at(ipar)->setSoftCut(m_softcutBarrelY);
               break;
            case Trk::AlignModule::TransZ:
               modPars->at(ipar)->setSigma(m_sigmaBarrelZ);
               modPars->at(ipar)->setSoftCut(m_softcutBarrelZ);
               break;
            case Trk::AlignModule::RotX:
               modPars->at(ipar)->setSigma(m_sigmaBarrelRotX);
               modPars->at(ipar)->setSoftCut(m_softcutBarrelRotX);
               break;
            case Trk::AlignModule::RotY:
               modPars->at(ipar)->setSigma(m_sigmaBarrelRotY);
               modPars->at(ipar)->setSoftCut(m_softcutBarrelRotY);
               break;
            case Trk::AlignModule::RotZ:
               modPars->at(ipar)->setSigma(m_sigmaBarrelRotZ);
               modPars->at(ipar)->setSoftCut(m_softcutBarrelRotZ);
               break;
            default:
               break;
         }
      // end-caps
      else
         switch(modPars->at(ipar)->paramType()) {
            case Trk::AlignModule::TransX:
               modPars->at(ipar)->setSigma(m_sigmaEndcapX);
               modPars->at(ipar)->setSoftCut(m_softcutEndcapX);
               break;
            case Trk::AlignModule::TransY:
               modPars->at(ipar)->setSigma(m_sigmaEndcapY);
               modPars->at(ipar)->setSoftCut(m_softcutEndcapY);
               break;
            case Trk::AlignModule::TransZ:
               modPars->at(ipar)->setSigma(m_sigmaEndcapZ);
               modPars->at(ipar)->setSoftCut(m_softcutEndcapZ);
               break;
            case Trk::AlignModule::RotX:
               modPars->at(ipar)->setSigma(m_sigmaEndcapRotX);
               modPars->at(ipar)->setSoftCut(m_softcutEndcapRotX);
               break;
            case Trk::AlignModule::RotY:
               modPars->at(ipar)->setSigma(m_sigmaEndcapRotY);
               modPars->at(ipar)->setSoftCut(m_softcutEndcapRotY);
               break;
            case Trk::AlignModule::RotZ:
               modPars->at(ipar)->setSigma(m_sigmaEndcapRotZ);
               modPars->at(ipar)->setSoftCut(m_softcutEndcapRotZ);
               break;
            default:
               break;
         }
   }
}

StoreModPositions()

void InDet::TRTGeometryManagerTool::StoreModPositions ( std::string )

inlinevirtual

Implements Trk::IGeometryManagerTool.

Definition at line 59 of file TRTGeometryManagerTool.h.

{};

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 >, and AthCheckedComponent<::AthAlgTool >.

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_alignBarrel

bool InDet::TRTGeometryManagerTool::m_alignBarrel

private

Definition at line 145 of file TRTGeometryManagerTool.h.

m_alignBarrelRotX

bool InDet::TRTGeometryManagerTool::m_alignBarrelRotX

private

Definition at line 151 of file TRTGeometryManagerTool.h.

m_alignBarrelRotY

bool InDet::TRTGeometryManagerTool::m_alignBarrelRotY

private

Definition at line 152 of file TRTGeometryManagerTool.h.

m_alignBarrelRotZ

bool InDet::TRTGeometryManagerTool::m_alignBarrelRotZ

private

Definition at line 153 of file TRTGeometryManagerTool.h.

m_alignBarrelX

bool InDet::TRTGeometryManagerTool::m_alignBarrelX

private

Definition at line 148 of file TRTGeometryManagerTool.h.

m_alignBarrelY

bool InDet::TRTGeometryManagerTool::m_alignBarrelY

private

Definition at line 149 of file TRTGeometryManagerTool.h.

m_alignBarrelZ

bool InDet::TRTGeometryManagerTool::m_alignBarrelZ

private

Definition at line 150 of file TRTGeometryManagerTool.h.

m_alignEndcapRotX

bool InDet::TRTGeometryManagerTool::m_alignEndcapRotX

private

Definition at line 158 of file TRTGeometryManagerTool.h.

m_alignEndcapRotY

bool InDet::TRTGeometryManagerTool::m_alignEndcapRotY

private

Definition at line 159 of file TRTGeometryManagerTool.h.

m_alignEndcapRotZ

bool InDet::TRTGeometryManagerTool::m_alignEndcapRotZ

private

Definition at line 160 of file TRTGeometryManagerTool.h.

m_alignEndcaps

bool InDet::TRTGeometryManagerTool::m_alignEndcaps

private

Definition at line 146 of file TRTGeometryManagerTool.h.

m_alignEndcapX

bool InDet::TRTGeometryManagerTool::m_alignEndcapX

private

Definition at line 155 of file TRTGeometryManagerTool.h.

m_alignEndcapY

bool InDet::TRTGeometryManagerTool::m_alignEndcapY

private

Definition at line 156 of file TRTGeometryManagerTool.h.

m_alignEndcapZ

bool InDet::TRTGeometryManagerTool::m_alignEndcapZ

private

Definition at line 157 of file TRTGeometryManagerTool.h.

m_alignLevel

int Trk::IGeometryManagerTool::m_alignLevel = 0

protectedinherited

alignment level

Definition at line 115 of file IGeometryManagerTool.h.

m_alignLevelBarrel

int Trk::IGeometryManagerTool::m_alignLevelBarrel = 0

protectedinherited

alignment level

Definition at line 116 of file IGeometryManagerTool.h.

m_alignLevelEndcaps

int Trk::IGeometryManagerTool::m_alignLevelEndcaps = 0

protectedinherited

alignment level

Definition at line 117 of file IGeometryManagerTool.h.

m_alignModuleList

Trk::AlignModuleList InDet::TRTGeometryManagerTool::m_alignModuleList

private

Definition at line 137 of file TRTGeometryManagerTool.h.

m_alignModuleListPtr

Trk::AlignModuleList* Trk::IGeometryManagerTool::m_alignModuleListPtr = nullptr

protectedinherited

pointer to module list to which the modules are added

Definition at line 106 of file IGeometryManagerTool.h.

m_alignModuleTool

ToolHandle<Trk::IAlignModuleTool> InDet::TRTGeometryManagerTool::m_alignModuleTool

private

pointer to AlignModuleTool

Definition at line 135 of file TRTGeometryManagerTool.h.

m_alignParList

DataVector< DataVector<Trk::AlignPar> >* InDet::TRTGeometryManagerTool::m_alignParList

private

Definition at line 142 of file TRTGeometryManagerTool.h.

m_barrelLayers

std::vector<int> InDet::TRTGeometryManagerTool::m_barrelLayers

private

list of barrel layers to align (if empty, align everything)

Definition at line 191 of file TRTGeometryManagerTool.h.

m_barrelPhiSectors

std::vector<int> InDet::TRTGeometryManagerTool::m_barrelPhiSectors

private

list of barrel phi sectors to align (if empty, align everything)

Definition at line 190 of file TRTGeometryManagerTool.h.

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_doModuleSelection

bool InDet::TRTGeometryManagerTool::m_doModuleSelection

private

Definition at line 197 of file TRTGeometryManagerTool.h.

m_dumpGeometry

bool InDet::TRTGeometryManagerTool::m_dumpGeometry

private

Definition at line 200 of file TRTGeometryManagerTool.h.

m_endcap

int InDet::TRTGeometryManagerTool::m_endcap

private

endcap_id of the endcap to align (-2 or 2 for endcap A or C, 0 for both)

Definition at line 193 of file TRTGeometryManagerTool.h.

m_endcapPhiSectors

std::vector<int> InDet::TRTGeometryManagerTool::m_endcapPhiSectors

private

list of endcap phi sectors to align (if empty, align everything)

Definition at line 194 of file TRTGeometryManagerTool.h.

m_endcapWheels

std::vector<int> InDet::TRTGeometryManagerTool::m_endcapWheels

private

list of endcap wheels to align (if empty, align everything)

Definition at line 195 of file TRTGeometryManagerTool.h.

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_fullAlignParList

DataVector< DataVector<Trk::AlignPar> >* InDet::TRTGeometryManagerTool::m_fullAlignParList

private

Definition at line 143 of file TRTGeometryManagerTool.h.

m_hashCounter

int Trk::IGeometryManagerTool::m_hashCounter = 0

protectedinherited

variable for setting the idHash of the AlignModules

Definition at line 111 of file IGeometryManagerTool.h.

m_idHashToAlignModuleMap

Trk::AlignModuleList InDet::TRTGeometryManagerTool::m_idHashToAlignModuleMap

private

Definition at line 138 of file TRTGeometryManagerTool.h.

m_idHashToAlignModuleMaps

std::vector<Trk::AlignModuleList*> InDet::TRTGeometryManagerTool::m_idHashToAlignModuleMaps

private

Definition at line 140 of file TRTGeometryManagerTool.h.

m_idHashToAlignModuleMapsPtr

std::vector<Trk::AlignModuleList *>* Trk::IGeometryManagerTool::m_idHashToAlignModuleMapsPtr = nullptr

protectedinherited

pointer to vector of hashMaps to which the elements are added

Definition at line 109 of file IGeometryManagerTool.h.

m_idHelper

const TRT_ID* InDet::TRTGeometryManagerTool::m_idHelper

private

Definition at line 133 of file TRTGeometryManagerTool.h.

m_logStream

std::ostream* Trk::IGeometryManagerTool::m_logStream = nullptr

protectedinherited

logfile output stream

Definition at line 113 of file IGeometryManagerTool.h.

m_moduleSelection

std::vector<unsigned long long> InDet::TRTGeometryManagerTool::m_moduleSelection

private

Definition at line 198 of file TRTGeometryManagerTool.h.

m_sigmaBarrelRotX

double InDet::TRTGeometryManagerTool::m_sigmaBarrelRotX

private

Definition at line 165 of file TRTGeometryManagerTool.h.

m_sigmaBarrelRotY

double InDet::TRTGeometryManagerTool::m_sigmaBarrelRotY

private

Definition at line 166 of file TRTGeometryManagerTool.h.

m_sigmaBarrelRotZ

double InDet::TRTGeometryManagerTool::m_sigmaBarrelRotZ

private

Definition at line 167 of file TRTGeometryManagerTool.h.

m_sigmaBarrelX

double InDet::TRTGeometryManagerTool::m_sigmaBarrelX

private

Definition at line 162 of file TRTGeometryManagerTool.h.

m_sigmaBarrelY

double InDet::TRTGeometryManagerTool::m_sigmaBarrelY

private

Definition at line 163 of file TRTGeometryManagerTool.h.

m_sigmaBarrelZ

double InDet::TRTGeometryManagerTool::m_sigmaBarrelZ

private

Definition at line 164 of file TRTGeometryManagerTool.h.

m_sigmaEndcapRotX

double InDet::TRTGeometryManagerTool::m_sigmaEndcapRotX

private

Definition at line 172 of file TRTGeometryManagerTool.h.

m_sigmaEndcapRotY

double InDet::TRTGeometryManagerTool::m_sigmaEndcapRotY

private

Definition at line 173 of file TRTGeometryManagerTool.h.

m_sigmaEndcapRotZ

double InDet::TRTGeometryManagerTool::m_sigmaEndcapRotZ

private

Definition at line 174 of file TRTGeometryManagerTool.h.

m_sigmaEndcapX

double InDet::TRTGeometryManagerTool::m_sigmaEndcapX

private

Definition at line 169 of file TRTGeometryManagerTool.h.

m_sigmaEndcapY

double InDet::TRTGeometryManagerTool::m_sigmaEndcapY

private

Definition at line 170 of file TRTGeometryManagerTool.h.

m_sigmaEndcapZ

double InDet::TRTGeometryManagerTool::m_sigmaEndcapZ

private

Definition at line 171 of file TRTGeometryManagerTool.h.

m_softcutBarrelRotX

double InDet::TRTGeometryManagerTool::m_softcutBarrelRotX

private

Definition at line 179 of file TRTGeometryManagerTool.h.

m_softcutBarrelRotY

double InDet::TRTGeometryManagerTool::m_softcutBarrelRotY

private

Definition at line 180 of file TRTGeometryManagerTool.h.

m_softcutBarrelRotZ

double InDet::TRTGeometryManagerTool::m_softcutBarrelRotZ

private

Definition at line 181 of file TRTGeometryManagerTool.h.

m_softcutBarrelX

double InDet::TRTGeometryManagerTool::m_softcutBarrelX

private

Definition at line 176 of file TRTGeometryManagerTool.h.

m_softcutBarrelY

double InDet::TRTGeometryManagerTool::m_softcutBarrelY

private

Definition at line 177 of file TRTGeometryManagerTool.h.

m_softcutBarrelZ

double InDet::TRTGeometryManagerTool::m_softcutBarrelZ

private

Definition at line 178 of file TRTGeometryManagerTool.h.

m_softcutEndcapRotX

double InDet::TRTGeometryManagerTool::m_softcutEndcapRotX

private

Definition at line 186 of file TRTGeometryManagerTool.h.

m_softcutEndcapRotY

double InDet::TRTGeometryManagerTool::m_softcutEndcapRotY

private

Definition at line 187 of file TRTGeometryManagerTool.h.

m_softcutEndcapRotZ

double InDet::TRTGeometryManagerTool::m_softcutEndcapRotZ

private

Definition at line 188 of file TRTGeometryManagerTool.h.

m_softcutEndcapX

double InDet::TRTGeometryManagerTool::m_softcutEndcapX

private

Definition at line 183 of file TRTGeometryManagerTool.h.

m_softcutEndcapY

double InDet::TRTGeometryManagerTool::m_softcutEndcapY

private

Definition at line 184 of file TRTGeometryManagerTool.h.

m_softcutEndcapZ

double InDet::TRTGeometryManagerTool::m_softcutEndcapZ

private

Definition at line 185 of file TRTGeometryManagerTool.h.

m_trtDetManager

const InDetDD::TRT_DetectorManager* InDet::TRTGeometryManagerTool::m_trtDetManager

private

pointer to TRT detector manager

Definition at line 131 of file TRTGeometryManagerTool.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:

• TRTGeometryManagerTool.h
• TRTGeometryManagerTool.cxx