InDet::SCTGeometryManagerTool Class Reference

#include <SCTGeometryManagerTool.h>

Inheritance diagram for InDet::SCTGeometryManagerTool:

Collaboration diagram for InDet::SCTGeometryManagerTool:

Public Member Functions
	SCTGeometryManagerTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~SCTGeometryManagerTool ()
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 SCT 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 ()
	creates L0 AlignModules for SCT
void	buildL1Barrel ()
	creates L1 AlignModules for SCT barrel
void	buildL1Endcaps ()
	creates L1 AlignModules for SCT endcaps
void	buildL2Barrel ()
	creates L2 AlignModules for SCT barrel
void	buildL27Barrel ()
	creates L27 AlignModules for SCT barrel
void	buildL2Endcaps ()
	creates L2 AlignModules for SCT endcaps
void	buildL25Endcaps ()
	creates L25 AlignModules for SCT endcaps
void	buildL2Old ()
	creates L2 AlignModules for the full Pixel in the ordering equal to the old code
void	buildL3Barrel ()
	creates L3 AlignModules for SCT barrel
void	buildL3Endcaps ()
	creates L3 AlignModules for SCT endcaps
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
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::SCT_DetectorManager *	m_detManager
	pointer to SCT detector manager
const SCT_ID *	m_idHelper
	pointer to SCT detector manager
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_doModuleSelection
std::vector< unsigned long long >	m_moduleSelection
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
bool	m_useOldL2
bool	m_dumpGeometry
Gaudi::Property< std::string >	m_stripDetManagerName {this,"StripDetectorManager","SCT","Name of the Strip Detector Manager to attempt to retrieve"}
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 43 of file SCTGeometryManagerTool.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

SCTGeometryManagerTool()

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

Definition at line 30 of file SCTGeometryManagerTool.cxx.

    : AthAlgTool(type,name,parent)
    , m_detManager(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);
 
    // Pixel joboptions
    declareProperty("AlignBarrel",    m_alignBarrel      = true);
    declareProperty("AlignEndcaps",   m_alignEndcaps     = true);
 
    declareProperty("AlignBarrelX",    m_alignBarrelX    = true);
    declareProperty("AlignBarrelY",    m_alignBarrelY    = true);
    declareProperty("AlignBarrelZ",    m_alignBarrelZ    = true);
    declareProperty("AlignBarrelRotX", m_alignBarrelRotX = true);
    declareProperty("AlignBarrelRotY", m_alignBarrelRotY = true);
    declareProperty("AlignBarrelRotZ", m_alignBarrelRotZ = 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.);
 
    // defines alignment level
    declareProperty("AlignmentLevel",        m_alignLevel        =  1);
    declareProperty("AlignmentLevelBarrel",  m_alignLevelBarrel  = -1);
    declareProperty("AlignmentLevelEndcaps", m_alignLevelEndcaps = -1);
 
    declareProperty("doModuleSelection",  m_doModuleSelection = false);
    declareProperty("ModuleSelection",    m_moduleSelection);
 
    declareProperty("UseOldL2",           m_useOldL2 = false);
    declareProperty("DumpGeometry",       m_dumpGeometry = true);
 
    m_hashCounter = 0;
    m_logStream = nullptr;
  }

~SCTGeometryManagerTool()

InDet::SCTGeometryManagerTool::~SCTGeometryManagerTool ( )

virtual

Definition at line 106 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::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 1234 of file SCTGeometryManagerTool.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 AlimodID = module->identify();
      Identifier modID = ((module->detElementCollection(Trk::AlignModule::SCT))->at(0))->identify();
 
      if(m_idHelper->is_barrel(modID)) {
        ATH_MSG_DEBUG("SCT barrel module with id "<<AlimodID);
        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("SCT endcap module with id "<<AlimodID);
        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");
      }
    }
 
    // now add parameters to the list
    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::SCTGeometryManagerTool::buildGeometry ( )

virtual

builds geometry for SCT alignment

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 286 of file SCTGeometryManagerTool.cxx.

  {
    ATH_MSG_INFO("Preparing the SCT geometry");
 
    if(m_alignLevel == 0)
      buildL0();
    else if (m_alignLevel == 2 && m_useOldL2)
      buildL2Old();
    else {
 
      if(m_alignBarrel)
        switch(m_alignLevelBarrel) {
          case 1:
            buildL1Barrel();
            break;
          case 2:
            buildL2Barrel();
            break;
          case 27:
            buildL27Barrel();
            break;
          case 3:
            buildL3Barrel();
            break;
          default:
            break;
        }
 
      if(m_alignEndcaps)
        switch(m_alignLevelEndcaps) {
          case 1:
            buildL1Endcaps();
            break;
          case 2:
            buildL2Endcaps();
            break;
          case 25:
            buildL25Endcaps();
            break;
          case 3:
            buildL3Endcaps();
            break;
          default:
            break;
        }
 
    }
  }

buildL0()

void InDet::SCTGeometryManagerTool::buildL0 ( )

private

creates L0 AlignModules for SCT

Definition at line 336 of file SCTGeometryManagerTool.cxx.

  {
    ATH_MSG_INFO("Preparing the SCT geometry for L0: the whole SCT");
    // ===================================================
    // Level 0 is just one module containing the whole SCT
 
    Trk::AlignModule * sct = new Trk::AlignModule(this);
    sct->setIdHash(getNextIDHash());
    sct->setName("SCT");
 
    // use the identifier of the SCT barrel for the whole SCT at L0
    sct->setIdentifier(m_idHelper->wafer_id(0,0,0,0,0));
 
    ATH_MSG_DEBUG("Created module "<<sct->name()<<"  idHash: "<<sct->identifyHash()<<"  identifier: "<<sct->identify());
 
    // for L0 alignment the alignment frame is equal to the global frame
    // and since the SCT detector element positions are also stored
    // in the global frame in DB, transform is identity
    const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
    // get maximum number of elements from the helper
    unsigned int sctmaxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT: "<<sctmaxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(sctmaxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    // ================================================================
    // loop over SCT elements and add them to respective alignModules
    // ================================================================
    for (unsigned int index = 0; index < sctmaxHash; index++) {
      IdentifierHash idHash = index;
      Identifier id = m_idHelper->wafer_id(idHash);
 
      ATH_MSG_DEBUG(" SCT DetectorElement idhash: "<<index);
      ATH_MSG_DEBUG(" DetectorElement id: "<<id);
 
      // get the element via hash
      const SiDetectorElement * element = m_detManager->getDetectorElement(id);
      if (element) {
        // add element to respective AlignModule
        // add to the sct barrel structure
        if(msgLvl(MSG::DEBUG)) {
          if (m_idHelper->is_barrel(id))
            msg(MSG::DEBUG)<<"... SCT barrel element"<<endmsg;
          else
            msg(MSG::DEBUG)<<"... SCT endcap element"<<endmsg;
        }
        sct->addDetElement(Trk::AlignModule::SCT,element,transform);
 
        // and fill the corresponding map
        (*sctIdHashMap)[idHash] = sct;
      }
      else
        ATH_MSG_DEBUG("No SCT detector with id: "<<id);
    }
 
    // add created module to the geometry
    m_alignModuleListPtr->push_back(sct);
 
    ATH_MSG_DEBUG("SCT L0 module successfully added to the geometry");
  }

buildL1Barrel()

void InDet::SCTGeometryManagerTool::buildL1Barrel ( )

private

creates L1 AlignModules for SCT barrel

Definition at line 401 of file SCTGeometryManagerTool.cxx.

  {
    // ===========================
    // BARREL
    // ===========================
    ATH_MSG_INFO("Preparing the SCT Barrel geometry for L1: barrel + 2 endcaps");
 
    // ========================================
    // Level 1 is whole SCT barrel
 
    // get maximum number of elements from the helper
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    // use the 'real' identifier of the Pixel
    // build barrel module
    Trk::AlignModule * mod = new Trk::AlignModule(this);
    mod->setIdHash(getNextIDHash());
    mod->setName("SCT Barrel");
    mod->setIdentifier(m_idHelper->wafer_id(0,0,0,0,0));
 
    // check if selected
    if(!m_alignBarrel || !moduleSelected(mod)) {
      ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
      delete mod;
      return;
    }
 
    // fill module with detector elements
    ATH_MSG_DEBUG("Building module "<<mod->name());
 
    const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
    for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
      if (!m_detManager->numerology().useLayer(iLayer))
        ATH_MSG_INFO("  Layer "<<iLayer<<" not present");
 
      for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
        ATH_MSG_DEBUG("iPhi "<<iPhi);
        for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
 
          if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
            // iEta=0 is not defined for the SCT
            continue;
 
          for (int is = 0; is < 2; is++) { // module side
            const SiDetectorElement* element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
 
            if (element) {
              // get element location for debugging
              // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
              // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
              // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
              // add element to the AlignModule
              mod->addDetElement(Trk::AlignModule::SCT,element,transform);
              // and fill the corresponding map
              (*sctIdHashMap)[element->identifyHash()] = mod;
            }
          }
        }
      }
    }
 
    m_alignModuleListPtr->push_back(mod);
 
    ATH_MSG_DEBUG("SCT Barrel successfully added to the geometry");
  }

buildL1Endcaps()

void InDet::SCTGeometryManagerTool::buildL1Endcaps ( )

private

creates L1 AlignModules for SCT endcaps

Definition at line 476 of file SCTGeometryManagerTool.cxx.

  {
    // ===========================
    // ENDCAPs
    // ===========================
    ATH_MSG_INFO("Preparing the SCT Endcap geometry for L1: 2 endcaps");
 
    // ========================================
    // Level 1 is 2 SCT endcaps
 
    // get maximum number of elements from the helper
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
      int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
 
      Trk::AlignModule * mod = new Trk::AlignModule(this);
      mod->setIdHash(getNextIDHash());
      mod->setIdentifier(m_idHelper->wafer_id(iSide,0,0,0,0));
      if(iEndcapIndex == 0)
        mod->setName("SCT EndCap A");
      else
        mod->setName("SCT EndCap C");
 
      if(!m_alignEndcaps || !moduleSelected(mod)) {
        ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
        delete mod;
        continue;
      }
 
      ATH_MSG_DEBUG("Building module "<<mod->name());
 
      const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
      for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
        if (!m_detManager->numerology().useDisk(iWheel))
          ATH_MSG_INFO("  Disk "<<iWheel<<" not present");
 
        for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
          for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
            ATH_MSG_DEBUG("iPhi "<<iPhi);
            for(int is=0; is<2; is++) { // module side
 
              const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
              if (element) {
                // get element location for debugging
                // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
                // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
                // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
                // add element to the AlignModule
                mod->addDetElement(Trk::AlignModule::SCT,element,transform);
                // and fill the corresponding map
                (*sctIdHashMap)[element->identifyHash()] =  mod;
              }
            }
          }
        }
      }
 
      // for endcap we move the CoG
      Amg::Translation3D translation(mod->centerOfGravity());
      ATH_MSG_DEBUG("Endcap:" << iSide << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
      Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
 
      mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
 
      m_alignModuleListPtr->push_back(mod);
    }
 
    ATH_MSG_DEBUG("SCT End-caps successfully added to the geometry");
  }

buildL25Endcaps()

void InDet::SCTGeometryManagerTool::buildL25Endcaps ( )

private

creates L25 AlignModules for SCT endcaps

Definition at line 805 of file SCTGeometryManagerTool.cxx.

  {
    // ========================================
    // END CAPS
    ATH_MSG_INFO("Preparing the SCT Endcap geometry for L25: 2 x 22 rings");
 
    if(!m_alignEndcaps) {
      ATH_MSG_DEBUG("Not aligning endcaps");
      return;
    }
 
    // ========================================
    // get all modules for Level 2 alignment of the Endcaps
    // 2 x 22 SCT Endcap rings
 
    const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
      int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
      for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
        if (!m_detManager->numerology().useDisk(iWheel))
          ATH_MSG_INFO("  Disk "<<iWheel<<" not present");
        for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
 
          Trk::AlignModule * mod = new Trk::AlignModule(this);
          mod->setIdHash(getNextIDHash());
          mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,iEta,0,0));
 
          std::stringstream name;
          if(iEndcapIndex == 0)
            name <<"SCT/EndcapA/Disk_"<<iWheel<<"/Ring_"<<iEta;
          else
            name <<"SCT/EndcapC/Disk_"<<iWheel<<"/Ring_"<<iEta;
          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 (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
            ATH_MSG_DEBUG("iPhi "<<iPhi);
            for(int is=0; is<2; is++) { // module side
 
              const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
              if (element) {
                // get element location for debugging
                // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
                // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
                // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
                // add element to the AlignModule
                mod->addDetElement(Trk::AlignModule::SCT,element,transform);
                // and fill the corresponding map
                (*sctIdHashMap)[element->identifyHash()] =  mod;
              }
            }
          }
 
          // for endcap we move the CoG
 
 
 
          Amg::Translation3D translation(mod->centerOfGravity());
          ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
          Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
          mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
 
          m_alignModuleListPtr->push_back(mod);
        }
      }
    }
  }

buildL27Barrel()

void InDet::SCTGeometryManagerTool::buildL27Barrel ( )

private

creates L27 AlignModules for SCT barrel

Definition at line 626 of file SCTGeometryManagerTool.cxx.

  {
    // ========================================
    // BARREL
    ATH_MSG_INFO("Preparing the SCT Barrel geometry for L27: 176 staves (ladders)");
 
    if(!m_alignBarrel) {
      ATH_MSG_DEBUG("Not aligning barrel");
      return;
    }
 
    // ========================================
    // get all modules for Level 27 alignment of the Barrel
    // 176 staves (ladders)
 
    const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
      if (!m_detManager->numerology().useLayer(iLayer))
        ATH_MSG_INFO("  Layer "<<iLayer<<" not present");
      for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
        ATH_MSG_DEBUG("iPhi "<<iPhi);
 
        // create the AlignModule
        Trk::AlignModule * mod = new Trk::AlignModule(this);
        mod->setIdHash(getNextIDHash());
 
        // even though there is no iEta=0 module for SCT barrel, the
        // Identifier is still valid so we use it for the stave
        mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, iPhi, 0, 0));
 
        std::stringstream name;
        name <<"SCT/Barrel/Layer_"<<iLayer<<"/PhiStave_"<<iPhi;
        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());
 
        // rotation of the stave alignment frame
        Amg::RotationMatrix3D rotation; rotation.setIdentity();
 
        for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
          if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
            // iEta=0 is not defined for the SCT
            continue;
          for (int is = 0; is < 2; is++) { // module side
 
            const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
            if (element) {
              // get element location for debugging
              // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
              // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
              // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
              // add element to the AlignModule
              mod->addDetElement(Trk::AlignModule::SCT,element,transform);
              // and fill the corresponding map
              (*sctIdHashMap)[element->identifyHash()] = mod;
            }
 
            // for the stave alignment frame rotation we use the one of the iEta=1
            // non-stereo side (which is the module local frame)
            if(iEta==1 && !element->isStereo())
              rotation = element->moduleTransform().rotation();
          }
        }
 
        // we set the alignment frame to be the CoG of the stave
        // with rotation being the one of the iEta=1 module set above
        Amg::Translation3D translation(mod->centerOfGravity());
        Amg::Transform3D localToGlobal = translation * rotation;
 
 
        ATH_MSG_DEBUG("Prepared local to global transform :");
        ATH_MSG_DEBUG(" - translation: "<<localToGlobal.translation().x()<<"  "<<localToGlobal.translation().y()<<"  "<<localToGlobal.translation().z());
        ATH_MSG_DEBUG(" - rotation:");
        ATH_MSG_DEBUG("      "<<localToGlobal.rotation()(0,0)<<"  "<<localToGlobal.rotation()(0,1)<<"  "<<localToGlobal.rotation()(0,2));
        ATH_MSG_DEBUG("      "<<localToGlobal.rotation()(1,0)<<"  "<<localToGlobal.rotation()(1,1)<<"  "<<localToGlobal.rotation()(1,2));
        ATH_MSG_DEBUG("      "<<localToGlobal.rotation()(2,0)<<"  "<<localToGlobal.rotation()(2,1)<<"  "<<localToGlobal.rotation()(2,2));
 
        mod->setGlobalFrameToAlignFrameTransform(localToGlobal.inverse());
 
        m_alignModuleListPtr->push_back(mod);
      }
    }
  }

buildL2Barrel()

void InDet::SCTGeometryManagerTool::buildL2Barrel ( )

private

creates L2 AlignModules for SCT barrel

Definition at line 556 of file SCTGeometryManagerTool.cxx.

  {
    // ========================================
    // BARREL
    ATH_MSG_INFO("Preparing the SCT Barrel geometry for L2: 4 layers");
 
    // ========================================
    // get all modules for Level 2 alignment of the Barrel
    // Level 2 is 4 SCT barrel layers
 
    const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
      if (!m_detManager->numerology().useLayer(iLayer))
        ATH_MSG_INFO("  Layer "<<iLayer<<" not present");
 
      // create the AlignModule
      Trk::AlignModule * mod = new Trk::AlignModule(this);
      mod->setIdHash(getNextIDHash());
      mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, 0, 0, 0));
 
      std::stringstream name;
      name <<"SCT/Barrel/Layer_"<<iLayer;
      mod->setName(name.str());
 
      if(!m_alignBarrel || !moduleSelected(mod)) {
        ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
        delete mod;
        continue;
      }
 
      ATH_MSG_DEBUG("Building module "<<mod->name());
 
      for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
        ATH_MSG_DEBUG("iPhi "<<iPhi);
        for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
          if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
            // iEta=0 is not defined for the SCT
            continue;
          for (int is = 0; is < 2; is++) { // module side
 
            const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
            if (element) {
              // get element location for debugging
              // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
              // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
              // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
              // add element to the AlignModule
              mod->addDetElement(Trk::AlignModule::SCT,element,transform);
              // and fill the corresponding map
              (*sctIdHashMap)[element->identifyHash()] = mod;
            }
          }
        }
      }
 
      m_alignModuleListPtr->push_back(mod);
    }
  }

buildL2Endcaps()

void InDet::SCTGeometryManagerTool::buildL2Endcaps ( )

private

creates L2 AlignModules for SCT endcaps

Definition at line 727 of file SCTGeometryManagerTool.cxx.

  {
    // ========================================
    // END CAPS
    ATH_MSG_INFO("Preparing the SCT Endcap geometry for L2: 2 x 9 disks");
 
    // ========================================
    // get all modules for Level 2 alignment of the Endcaps
    // Level 2 is 2 x 9 SCT Endcap disks
 
    const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
      int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
      for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
        if (!m_detManager->numerology().useDisk(iWheel))
          ATH_MSG_INFO("  Disk "<<iWheel<<" not present");
 
        Trk::AlignModule * mod = new Trk::AlignModule(this);
        mod->setIdHash(getNextIDHash());
        mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,0,0,0));
 
        std::stringstream name;
        if(iEndcapIndex == 0)
          name <<"SCT/EndcapA/Disk_" << iWheel;
        else
          name <<"SCT/EndcapC/Disk_" << iWheel;
        mod->setName(name.str());
 
        if(!m_alignEndcaps || !moduleSelected(mod)) {
          ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
          delete mod;
          continue;
        }
 
        ATH_MSG_DEBUG("Building module "<<mod->name());
 
        for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
          for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
            ATH_MSG_DEBUG("iPhi "<<iPhi);
            for(int is=0; is<2; is++) { // module side
 
              const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
              if (element) {
                // get element location for debugging
                // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
                // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
                // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
                // add element to the AlignModule
                mod->addDetElement(Trk::AlignModule::SCT,element,transform);
                // and fill the corresponding map
                (*sctIdHashMap)[element->identifyHash()] =  mod;
              }
            }
          }
        }
 
        // for endcap we move the CoG
        Amg::Translation3D translation(mod->centerOfGravity());
        ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
        Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
        mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
 
        m_alignModuleListPtr->push_back(mod);
      }
    }
  }

buildL2Old()

void InDet::SCTGeometryManagerTool::buildL2Old ( )

private

creates L2 AlignModules for the full Pixel in the ordering equal to the old code

Definition at line 891 of file SCTGeometryManagerTool.cxx.

  {
    ATH_MSG_INFO("Preparing the old SCT Barrel geometry for L2: 4 barrel layers + 2 x 9 endcap disks");
 
    const Amg::Transform3D transform = Amg::Transform3D::Identity();
 
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    // first endcaap C
    int iEndcapIndex = 1;
    int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
    for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
      if (!m_detManager->numerology().useDisk(iWheel))
        ATH_MSG_INFO("  Disk "<<iWheel<<" not present");
 
      Trk::AlignModule * mod = new Trk::AlignModule(this);
      mod->setIdHash(getNextIDHash());
      mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,0,0,0));
 
      std::stringstream name;
      name <<"SCT/EndcapC/Disk_" << iWheel;
      mod->setName(name.str());
 
      if(!m_alignEndcaps || !moduleSelected(mod)) {
        ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
        delete mod;
        continue;
      }
 
      ATH_MSG_DEBUG("Building module "<<mod->name());
 
      for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
        ATH_MSG_DEBUG("iEta "<<iEta);
        for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
          ATH_MSG_DEBUG("iPhi "<<iPhi);
          for(int is=0; is<2; is++) { // module side
            const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
            if (element) {
              // get element location for debugging
              // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
              // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
              // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
              // add element to the AlignModule
              mod->addDetElement(Trk::AlignModule::SCT,element,transform);
              // and fill the corresponding map
              (*sctIdHashMap)[element->identifyHash()] =  mod;
            }
          }
        }
      }
 
      // for endcap we move the CoG
      Amg::Translation3D translation(mod->centerOfGravity());
      ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
      Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
      mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
 
      m_alignModuleListPtr->push_back(mod);
    }
 
    // now barrel
    for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
      if (!m_detManager->numerology().useLayer(iLayer))
        ATH_MSG_INFO("  Layer "<<iLayer<<" not present");
 
      // create the AlignModule
      Trk::AlignModule * mod = new Trk::AlignModule(this);
      mod->setIdHash(getNextIDHash());
      mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, 0, 0, 0));
 
      std::stringstream name;
      name <<"SCT/Barrel/Layer_"<<iLayer;
      mod->setName(name.str());
 
      if(!m_alignBarrel || !moduleSelected(mod)) {
        ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
        delete mod;
        continue;
      }
 
      ATH_MSG_DEBUG("Building module "<<mod->name());
 
      for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
        ATH_MSG_DEBUG("iPhi "<<iPhi);
        for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
          if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
            // iEta=0 is not defined for the SCT
            continue;
          for (int is = 0; is < 2; is++) { // module side
 
            const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
            if (element) {
              // get element location for debugging
              // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
              // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
              // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
              // add element to the AlignModule
              mod->addDetElement(Trk::AlignModule::SCT,element,transform);
              // and fill the corresponding map
              (*sctIdHashMap)[element->identifyHash()] = mod;
            }
          }
        }
      }
 
      m_alignModuleListPtr->push_back(mod);
    }
 
    // and endcaap A
    iEndcapIndex = 0;
    iSide = m_detManager->numerology().endcapId(iEndcapIndex);
    for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
      if (!m_detManager->numerology().useDisk(iWheel))
        ATH_MSG_INFO("  Disk "<<iWheel<<" not present");
 
      Trk::AlignModule * mod = new Trk::AlignModule(this);
      mod->setIdHash(getNextIDHash());
      mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,0,0,0));
 
      std::stringstream name;
      name <<"SCT/EndcapA/Disk_" << iWheel;
      mod->setName(name.str());
 
      if(!m_alignEndcaps || !moduleSelected(mod)) {
        ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
        delete mod;
        continue;
      }
 
      ATH_MSG_DEBUG("Building module "<<mod->name());
 
      for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
        ATH_MSG_DEBUG("iEta "<<iEta);
        for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
          ATH_MSG_DEBUG("iPhi "<<iPhi);
          for(int is=0; is<2; is++) { // module side
            const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
            if (element) {
              // get element location for debugging
              // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
              // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
              // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
              // add element to the AlignModule
              mod->addDetElement(Trk::AlignModule::SCT,element,transform);
              // and fill the corresponding map
              (*sctIdHashMap)[element->identifyHash()] =  mod;
            }
          }
        }
      }
 
      // for endcap we move the CoG
      Amg::Translation3D translation(mod->centerOfGravity());
      ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
 
      Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
 
      mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
 
      m_alignModuleListPtr->push_back(mod);
    }
  }

buildL3Barrel()

void InDet::SCTGeometryManagerTool::buildL3Barrel ( )

private

creates L3 AlignModules for SCT barrel

Definition at line 1064 of file SCTGeometryManagerTool.cxx.

  {
    // ========================================
    // BARREL
    ATH_MSG_INFO("Preparing the SCT Barrel geometry for L3: module (bi-wafer) level");
 
    // =========================================
    // get all modules for SCT Barrel Level 3 alignment
    // 2112 modules
 
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
      if (!m_detManager->numerology().useLayer(iLayer))
        ATH_MSG_INFO("  Layer "<<iLayer<<" not present");
      for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
        ATH_MSG_DEBUG("iPhi "<<iPhi);
        for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
          if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
            // iEta=0 is not defined for the SCT
            continue;
 
          // create the AlignModule
          Trk::AlignModule * mod = new Trk::AlignModule(this);
          mod->setIdHash(getNextIDHash());
          mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, iPhi, iEta, 0));
 
          std::stringstream name;
          name <<"SCT/Barrel/Layer_"<<iLayer<<"/Phi_"<<iPhi<<"/Eta_"<<iEta;
          mod->setName(name.str());
 
          if(!m_alignBarrel || !moduleSelected(mod)) {
            ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
            delete mod;
            continue;
          }
 
          ATH_MSG_DEBUG("Building module "<<mod->name());
 
          for(int is=0;is<2;is++) { // module side
 
            const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
 
            if (element) {
              // get element location for debugging
              // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
              // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
              // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
              // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
              // should this be identity? if yes, it means that the
              // db frame constants are in the global frame. Or not?
              const Amg::Transform3D transform = Amg::Transform3D::Identity();
              // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
              // add element to the AlignModule
              mod->addDetElement(Trk::AlignModule::SCT,element,transform);
              // and fill the corresponding map
              (*sctIdHashMap)[element->identifyHash()] = mod;
 
              // the local frame for the SCT barrel modules is always that
              // of the Rphi element, NOT stereo !
              // in any case the correct frame can be obtained using moduleTransform() call
              // regardless of the side, but we only need it once so we retrieve it
              // for the non-stereo side
              if(!element->isStereo())
                mod->setGlobalFrameToAlignFrameTransform(element->moduleTransform().inverse());
            }
          }
 
          m_alignModuleListPtr->push_back(mod);
        }
      }
 
    }
  }

buildL3Endcaps()

void InDet::SCTGeometryManagerTool::buildL3Endcaps ( )

private

creates L3 AlignModules for SCT endcaps

Definition at line 1148 of file SCTGeometryManagerTool.cxx.

  {
    // ========================================
    // ENDCAPs
    ATH_MSG_INFO("Preparing the SCT Endcap geometry for L3: module (bi-wafer) level");
 
    // =========================================
    // get all modules for SCT Endcap Level 3 alignment
    // 2 * 988
 
    unsigned int maxHash = m_idHelper->wafer_hash_max();
    ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
 
    if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
      m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
    Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
 
    for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
      int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
      for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
        if (!m_detManager->numerology().useDisk(iWheel))
          ATH_MSG_INFO("  Disk "<<iWheel<<" not present");
        for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
          ATH_MSG_DEBUG("iEta "<<iEta);
          for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
            ATH_MSG_DEBUG("iPhi "<<iPhi);
 
            Trk::AlignModule * mod = new Trk::AlignModule(this);
            mod->setIdHash(getNextIDHash());
            mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,iPhi,iEta,0));
 
            std::stringstream name;
            if(iEndcapIndex == 0)
              name<<"SCT/EndcapA";
            else
              name<<"SCT/EndcapC";
            name<<"/Disk_"<<iWheel<<"/Phi_"<<iPhi<<"/Eta_"<<iEta;
            mod->setName(name.str());
 
            if(!m_alignEndcaps || !moduleSelected(mod)) {
              ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
              delete mod;
              continue;
            }
 
            ATH_MSG_DEBUG("Building module "<<mod->name());
 
            for(int is=0; is<2; is++) { // module side
 
              const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
              if (element) {
                // get element location for debugging
                // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
                // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
                // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
 
                // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                // should this be identity? if yes, it means that the
                // db frame constants are in the global frame. Or not?
                const Amg::Transform3D transform = Amg::Transform3D::Identity();
                // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
                // add element to the AlignModule
                mod->addDetElement(Trk::AlignModule::SCT,element,transform);
                // and fill the corresponding map
                (*sctIdHashMap)[element->identifyHash()] =  mod;
 
                // the local frame for the SCT endcap modules is always that
                // of the Rphi element, NOT stereo !
                // in any case the correct frame can be obtained using moduleTransform() call
                // regardless of the side, but we only need it once so we retrieve it
                // for the non-stereo side
                if(!element->isStereo())
                  mod->setGlobalFrameToAlignFrameTransform(element->moduleTransform().inverse());
              }
            }
 
            m_alignModuleListPtr->push_back(mod);
          }
        }
      }
    }
 
     }

checkAlignLevel()

bool InDet::SCTGeometryManagerTool::checkAlignLevel ( )

virtual

check whether the alignment level is correct

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 178 of file SCTGeometryManagerTool.cxx.

  {
    if(m_alignLevel == 0) {
      ATH_MSG_INFO("Alignment level for SCT is "<<m_alignLevel);
      return true;
    }
 
    if(m_alignLevel == 2 && m_useOldL2) {
      ATH_MSG_INFO("Setting up OLD L2 alignment level for SCT");
      m_alignLevelBarrel = m_alignLevel;
      m_alignLevelEndcaps = m_alignLevel;
      return true;
    }
 
    if(m_alignLevelBarrel == -1)
      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::SCTGeometryManagerTool::checkAlignLevelBarrel ( )

virtual

check whether the Barrel alignment level is correct

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 209 of file SCTGeometryManagerTool.cxx.

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

checkAlignLevelEndcaps()

bool InDet::SCTGeometryManagerTool::checkAlignLevelEndcaps ( )

virtual

check whether the Endcap alignment level is correct

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 226 of file SCTGeometryManagerTool.cxx.

  {
    if(!m_alignEndcaps)
      return true;
 
    // check whether geometry level is allowed
    switch(m_alignLevelEndcaps) {
      case 1: case 2: case 25: case 3:
        ATH_MSG_INFO("Alignment level for SCT Endcaps is "<<m_alignLevelBarrel);
        return true;
      default:
        msg(MSG::FATAL)<<"Alignment level "<<m_alignLevelEndcaps<<" does not exist for SCT 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::SCTGeometryManagerTool::dumpGeometry ( )

private

print basic geometry info to screen

Definition at line 1378 of file SCTGeometryManagerTool.cxx.

  {
    ATH_MSG_INFO("---------------------------------------------------");
    ATH_MSG_INFO("Summary of the alignment geometry");
    ATH_MSG_INFO("Number of alignable objects: "<< m_alignModuleList.size());
    for(unsigned int i=0;i<m_alignModuleList.size();i++) {
      const Trk::AlignModule* module = m_alignModuleList.at(i);
      ATH_MSG_INFO(i<<". "<< module->name());
      ATH_MSG_INFO("   - identifier: "<<module->identify());
      ATH_MSG_INFO("   - has "<<module->detElementCollection(Trk::AlignModule::SCT)->size()<<" SCT 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::SCTGeometryManagerTool::finalize ( )

virtual

Implements Trk::IGeometryManagerTool.

Definition at line 163 of file SCTGeometryManagerTool.cxx.

  {
    ATH_MSG_DEBUG("finalize() of SCTGeometryManagerTool");
 
    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::SCTGeometryManagerTool::initialize ( )

virtual

Implements Trk::IGeometryManagerTool.

Definition at line 120 of file SCTGeometryManagerTool.cxx.

  {
    ATH_MSG_DEBUG("initialize() of SCTGeometryManagerTool");
 
    // 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 SCT helper
    if ( detStore()->retrieve(m_idHelper).isFailure() ) {
      msg(MSG::FATAL) << " Cannot retrieve SCT Helper " << endmsg;
      return StatusCode::FAILURE;
    }
    else
      ATH_MSG_INFO("retrieved Silicon SCT Helper");
 
    // retrieve SCT detector manager
    if ( detStore()->retrieve(m_detManager,m_stripDetManagerName).isFailure() ) {
      msg(MSG::FATAL) << " Cannot retrieve SCT Detector Manager " << endmsg;
      return StatusCode::FAILURE;
    }
    else
      ATH_MSG_INFO("retrieved SCT 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::SCTGeometryManagerTool::moduleSelected ( Trk::AlignModule * mod )

private

check wheather module is selected for module pointer

Definition at line 1372 of file SCTGeometryManagerTool.cxx.

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

moduleSelected() [2/2]

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

private

check wheather module is selected for module id

Definition at line 1358 of file SCTGeometryManagerTool.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;
  }

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::SCTGeometryManagerTool::PrintModPositions ( )

inlinevirtual

Implements Trk::IGeometryManagerTool.

Definition at line 59 of file SCTGeometryManagerTool.h.

{};

ReadGeometry()

int InDet::SCTGeometryManagerTool::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 243 of file SCTGeometryManagerTool.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
    ATH_MSG_DEBUG("Adding module parameters to 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::SCTGeometryManagerTool::setAlignLevel ( int level )

virtual

set alignment level

Reimplemented from Trk::IGeometryManagerTool.

Definition at line 171 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::setSigmas ( Trk::AlignModule * mod, DataVector< Trk::AlignPar > * modPars )

private

sets sigmas for modules

Definition at line 1291 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::StoreModPositions ( std::string )

inlinevirtual

Implements Trk::IGeometryManagerTool.

Definition at line 60 of file SCTGeometryManagerTool.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 >, 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_alignBarrel

bool InDet::SCTGeometryManagerTool::m_alignBarrel

private

Definition at line 142 of file SCTGeometryManagerTool.h.

m_alignBarrelRotX

bool InDet::SCTGeometryManagerTool::m_alignBarrelRotX

private

Definition at line 148 of file SCTGeometryManagerTool.h.

m_alignBarrelRotY

bool InDet::SCTGeometryManagerTool::m_alignBarrelRotY

private

Definition at line 149 of file SCTGeometryManagerTool.h.

m_alignBarrelRotZ

bool InDet::SCTGeometryManagerTool::m_alignBarrelRotZ

private

Definition at line 150 of file SCTGeometryManagerTool.h.

m_alignBarrelX

bool InDet::SCTGeometryManagerTool::m_alignBarrelX

private

Definition at line 145 of file SCTGeometryManagerTool.h.

m_alignBarrelY

bool InDet::SCTGeometryManagerTool::m_alignBarrelY

private

Definition at line 146 of file SCTGeometryManagerTool.h.

m_alignBarrelZ

bool InDet::SCTGeometryManagerTool::m_alignBarrelZ

private

Definition at line 147 of file SCTGeometryManagerTool.h.

m_alignEndcapRotX

bool InDet::SCTGeometryManagerTool::m_alignEndcapRotX

private

Definition at line 155 of file SCTGeometryManagerTool.h.

m_alignEndcapRotY

bool InDet::SCTGeometryManagerTool::m_alignEndcapRotY

private

Definition at line 156 of file SCTGeometryManagerTool.h.

m_alignEndcapRotZ

bool InDet::SCTGeometryManagerTool::m_alignEndcapRotZ

private

Definition at line 157 of file SCTGeometryManagerTool.h.

m_alignEndcaps

bool InDet::SCTGeometryManagerTool::m_alignEndcaps

private

Definition at line 143 of file SCTGeometryManagerTool.h.

m_alignEndcapX

bool InDet::SCTGeometryManagerTool::m_alignEndcapX

private

Definition at line 152 of file SCTGeometryManagerTool.h.

m_alignEndcapY

bool InDet::SCTGeometryManagerTool::m_alignEndcapY

private

Definition at line 153 of file SCTGeometryManagerTool.h.

m_alignEndcapZ

bool InDet::SCTGeometryManagerTool::m_alignEndcapZ

private

Definition at line 154 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::m_alignModuleList

private

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

private

pointer to AlignModuleTool

Definition at line 129 of file SCTGeometryManagerTool.h.

m_alignParList

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

private

Definition at line 136 of file SCTGeometryManagerTool.h.

m_detManager

const InDetDD::SCT_DetectorManager* InDet::SCTGeometryManagerTool::m_detManager

private

pointer to SCT detector manager

Definition at line 125 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::m_doModuleSelection

private

Definition at line 139 of file SCTGeometryManagerTool.h.

m_dumpGeometry

bool InDet::SCTGeometryManagerTool::m_dumpGeometry

private

Definition at line 188 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::m_fullAlignParList

private

Definition at line 137 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::m_idHashToAlignModuleMap

private

Definition at line 132 of file SCTGeometryManagerTool.h.

m_idHashToAlignModuleMaps

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

private

Definition at line 134 of file SCTGeometryManagerTool.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 SCT_ID* InDet::SCTGeometryManagerTool::m_idHelper

private

pointer to SCT detector manager

Definition at line 127 of file SCTGeometryManagerTool.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::SCTGeometryManagerTool::m_moduleSelection

private

Definition at line 140 of file SCTGeometryManagerTool.h.

m_sigmaBarrelRotX

double InDet::SCTGeometryManagerTool::m_sigmaBarrelRotX

private

Definition at line 162 of file SCTGeometryManagerTool.h.

m_sigmaBarrelRotY

double InDet::SCTGeometryManagerTool::m_sigmaBarrelRotY

private

Definition at line 163 of file SCTGeometryManagerTool.h.

m_sigmaBarrelRotZ

double InDet::SCTGeometryManagerTool::m_sigmaBarrelRotZ

private

Definition at line 164 of file SCTGeometryManagerTool.h.

m_sigmaBarrelX

double InDet::SCTGeometryManagerTool::m_sigmaBarrelX

private

Definition at line 159 of file SCTGeometryManagerTool.h.

m_sigmaBarrelY

double InDet::SCTGeometryManagerTool::m_sigmaBarrelY

private

Definition at line 160 of file SCTGeometryManagerTool.h.

m_sigmaBarrelZ

double InDet::SCTGeometryManagerTool::m_sigmaBarrelZ

private

Definition at line 161 of file SCTGeometryManagerTool.h.

m_sigmaEndcapRotX

double InDet::SCTGeometryManagerTool::m_sigmaEndcapRotX

private

Definition at line 169 of file SCTGeometryManagerTool.h.

m_sigmaEndcapRotY

double InDet::SCTGeometryManagerTool::m_sigmaEndcapRotY

private

Definition at line 170 of file SCTGeometryManagerTool.h.

m_sigmaEndcapRotZ

double InDet::SCTGeometryManagerTool::m_sigmaEndcapRotZ

private

Definition at line 171 of file SCTGeometryManagerTool.h.

m_sigmaEndcapX

double InDet::SCTGeometryManagerTool::m_sigmaEndcapX

private

Definition at line 166 of file SCTGeometryManagerTool.h.

m_sigmaEndcapY

double InDet::SCTGeometryManagerTool::m_sigmaEndcapY

private

Definition at line 167 of file SCTGeometryManagerTool.h.

m_sigmaEndcapZ

double InDet::SCTGeometryManagerTool::m_sigmaEndcapZ

private

Definition at line 168 of file SCTGeometryManagerTool.h.

m_softcutBarrelRotX

double InDet::SCTGeometryManagerTool::m_softcutBarrelRotX

private

Definition at line 176 of file SCTGeometryManagerTool.h.

m_softcutBarrelRotY

double InDet::SCTGeometryManagerTool::m_softcutBarrelRotY

private

Definition at line 177 of file SCTGeometryManagerTool.h.

m_softcutBarrelRotZ

double InDet::SCTGeometryManagerTool::m_softcutBarrelRotZ

private

Definition at line 178 of file SCTGeometryManagerTool.h.

m_softcutBarrelX

double InDet::SCTGeometryManagerTool::m_softcutBarrelX

private

Definition at line 173 of file SCTGeometryManagerTool.h.

m_softcutBarrelY

double InDet::SCTGeometryManagerTool::m_softcutBarrelY

private

Definition at line 174 of file SCTGeometryManagerTool.h.

m_softcutBarrelZ

double InDet::SCTGeometryManagerTool::m_softcutBarrelZ

private

Definition at line 175 of file SCTGeometryManagerTool.h.

m_softcutEndcapRotX

double InDet::SCTGeometryManagerTool::m_softcutEndcapRotX

private

Definition at line 183 of file SCTGeometryManagerTool.h.

m_softcutEndcapRotY

double InDet::SCTGeometryManagerTool::m_softcutEndcapRotY

private

Definition at line 184 of file SCTGeometryManagerTool.h.

m_softcutEndcapRotZ

double InDet::SCTGeometryManagerTool::m_softcutEndcapRotZ

private

Definition at line 185 of file SCTGeometryManagerTool.h.

m_softcutEndcapX

double InDet::SCTGeometryManagerTool::m_softcutEndcapX

private

Definition at line 180 of file SCTGeometryManagerTool.h.

m_softcutEndcapY

double InDet::SCTGeometryManagerTool::m_softcutEndcapY

private

Definition at line 181 of file SCTGeometryManagerTool.h.

m_softcutEndcapZ

double InDet::SCTGeometryManagerTool::m_softcutEndcapZ

private

Definition at line 182 of file SCTGeometryManagerTool.h.

m_stripDetManagerName

Gaudi::Property<std::string> InDet::SCTGeometryManagerTool::m_stripDetManagerName {this,"StripDetectorManager","SCT","Name of the Strip Detector Manager to attempt to retrieve"}

private

Definition at line 190 of file SCTGeometryManagerTool.h.

{this,"StripDetectorManager","SCT","Name of the Strip Detector Manager to attempt to retrieve"};

m_useOldL2

bool InDet::SCTGeometryManagerTool::m_useOldL2

private

Definition at line 187 of file SCTGeometryManagerTool.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:

• SCTGeometryManagerTool.h
• SCTGeometryManagerTool.cxx