InDetDD::PixelDetectorManager Class Reference

Dedicated detector manager extending the functionality of the SiDetectorManager with dedicated pixel information, access. More...

#include <PixelDetectorManager.h>

Inheritance diagram for InDetDD::PixelDetectorManager:

Collaboration diagram for InDetDD::PixelDetectorManager:

Public Member Functions
	PixelDetectorManager (StoreGateSvc *detStore)
	Constructor.
	PixelDetectorManager (StoreGateSvc *detStore, const std::string &name, const std::string &pixelIDName)
	Constructor with name.
virtual unsigned int	getNumTreeTops () const override
	Access to raw geometry:
virtual PVConstLink	getTreeTop (unsigned int i) const override
void	addTreeTop (const PVConstLink &vol)
	Add a Tree top:
virtual const SiDetectorElement *	getDetectorElement (const Identifier &id) const override
	access to individual elements : via Identifier
virtual const SiDetectorElement *	getDetectorElement (const IdentifierHash &idHash) const override
	access to individual elements : via IdentifierHash
const SiDetectorElement *	getDetectorElement (int barrel_endcap, int layer_wheel, int phi_module, int eta_module) const
	access to individual elements : via element identification
virtual const SiDetectorElementCollection *	getDetectorElementCollection () const override
	access to whole collection via Iterators
virtual SiDetectorElementCollection::const_iterator	getDetectorElementBegin () const override
virtual SiDetectorElementCollection::const_iterator	getDetectorElementEnd () const override
virtual SiDetectorElementCollection::iterator	getDetectorElementBegin () override
virtual SiDetectorElementCollection::iterator	getDetectorElementEnd () override
virtual void	addDetectorElement (SiDetectorElement *element) override
	Add elememts.
virtual void	addAlignableTransform (int level, const Identifier &id, GeoAlignableTransform *xf, const GeoVFullPhysVol *child)
	Add alignable transforms.
virtual void	addAlignableTransform (int level, const Identifier &id, GeoAlignableTransform *xf, const GeoVPhysVol *child)
	As above but does a dynamic_cast to GeoVFullPhysVol.
virtual void	addAlignableTransform (int, const Identifier &, GeoAlignableTransform *) override
	Add alignable transforms.
virtual void	initNeighbours () override
	Initialize the neighbours.
virtual bool	isPixel () const override
	Methods to query which manager we have.
virtual bool	identifierBelongs (const Identifier &id) const override
	Check identifier is for this detector.
unsigned int	getBarrelLayers () const
	Layers built.
unsigned int	getEndcapLayers () const
	Layers built.
virtual const PixelModuleDesign *	getPixelDesign (int i) const
	Access to module design, Casts to PixelModuleDesign.
bool	processSpecialAlignment (const std::string &, InDetDD::AlignFolderType) override
	Process new IBLDist DB folder.
bool	processSpecialAlignment (const std::string &key, const CondAttrListCollection *obj=nullptr, GeoVAlignmentStore *alignStore=nullptr) const override
virtual bool	processGlobalAlignment (const std::string &, int level, FrameType frame, const CondAttrListCollection *obj, GeoVAlignmentStore *alignStore) const override
	Process new global DB folders for L1 and L2.
const std::string &	tag () const
	Get tag used in dictionary.
bool	isSCT () const
virtual void	invalidateAll ()
	Invalidate cache for all detector elements.
virtual void	updateAll () const
	Update all caches.
const SiDetectorDesign *	addDesign (std::unique_ptr< const SiDetectorDesign > &&)
	Access to module design; returns an observer pointer.
int	numDesigns () const
const SiDetectorDesign *	getDesign (int i) const
const SiNumerology &	numerology () const
	Access Numerology.
SiNumerology &	numerology ()
void	setCommonItems (std::unique_ptr< const SiCommonItems > &&commonItems)
	Set SiCommonItems.
const Version &	getVersion () const
	Get version information.
const std::string &	getLayout () const
void	setVersion (const Version &version)
void	addChannel (const std::string &key, int level, FrameType frame)
	Alignment access.
void	addFolder (const std::string &key)
void	addSpecialFolder (const std::string &key)
void	addGlobalFolder (const std::string &key)
void	addAlignFolderType (const AlignFolderType alignfolder)
StatusCode	align (IOVSVC_CALLBACK_ARGS)
StatusCode	align (const RawAlignmentObjects &alignObjects, GeoVAlignmentStore *alignStore) const
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

Static Public Member Functions
static bool	setAlignableTransformGlobalDelta (ExtendedAlignableTransform *extXF, const Amg::Transform3D &delta, GeoVAlignmentStore *alignStore=nullptr)
	Helper method to set delta transform from a global delta - Amg interface.
static bool	setAlignableTransformLocalDelta (ExtendedAlignableTransform *extXF, const Amg::Transform3D &localToGlobalXF, const Amg::Transform3D &delta, GeoVAlignmentStore *alignStore=nullptr)
	Helper method to set delta transform from a local delta - Amg interface.

Public Attributes
AlignFolderType	m_alignfoldertype

Protected Attributes
StoreGateSvc *	m_detStore
std::atomic_bool	m_suppressWarnings

Private Types
typedef std::map< Identifier, std::unique_ptr< ExtendedAlignableTransform > >	AlignableTransformMap

Private Member Functions
virtual bool	setAlignableTransformDelta (int level, const Identifier &id, const Amg::Transform3D &delta, FrameType frame, GeoVAlignmentStore *alignStore) const override
	implements the main alignment update for delta transforms in different frames, it translates into the LocalDelta or GlobalDelta function of SiDetectorManager
const PixelDetectorManager &	operator= (const PixelDetectorManager &right)
	prevent copy and assignment
	PixelDetectorManager (const PixelDetectorManager &right)
virtual const PixelID *	getIdHelper () const override
void	printTransform (MSG::Level level, const Amg::Transform3D &tr) const
const LevelInfo &	getLevel (const std::string &key) const
	Retrieve level information.
bool	processAlignmentContainer (const std::string &key) const
	return align folder string to use
bool	processAlignmentContainer (const AlignableTransformContainer *container, GeoVAlignmentStore *alignStore) const
bool	processKey (const std::string &key, const AlignableTransform *transformCollection, GeoVAlignmentStore *alignStore=nullptr) const
	Called by processAlignmentContainer, applies only one key on the transform Collections.
bool	processGlobalAlignmentContainer (const std::string &key, const CondAttrListCollection *obj=nullptr, GeoVAlignmentStore *alignStore=nullptr) const
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
std::vector< PVConstLink >	m_volume
SiDetectorElementCollection	m_elementCollection
std::vector< AlignableTransformMap >	m_higherAlignableTransforms
std::vector< std::unique_ptr< ExtendedAlignableTransform > >	m_alignableTransforms
const PixelID *	m_idHelper
bool	m_isLogical
	This variable switches the how the local alignment corrections are applied If true they will be calcualted on top of all of other corrections but in the default reference frame If false they will be calcualted on top of all of other corrections but in the globally aligned reference frame.
std::string	m_tag
SiNumerology	m_numerology
std::vector< std::unique_ptr< const SiDetectorDesign > >	m_designs
std::unique_ptr< const SiCommonItems >	m_commonItems
Version	m_version
std::map< std::string, LevelInfo >	m_keys
std::set< std::string >	m_folders
std::set< std::string >	m_specialFolders
std::set< std::string >	m_globalFolders
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels)
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging)

Static Private Attributes
static const LevelInfo	s_invalidLevel

Detailed Description

Dedicated detector manager extending the functionality of the SiDetectorManager with dedicated pixel information, access.

Author

: Grant Gorfine

• modified and maintained by Nick Styles & Andreas Salzburger

Definition at line 47 of file PixelDetectorManager.h.

Member Typedef Documentation

AlignableTransformMap

typedef std::map<Identifier, std::unique_ptr<ExtendedAlignableTransform> > InDetDD::PixelDetectorManager::AlignableTransformMap

private

Definition at line 164 of file PixelDetectorManager.h.

Constructor & Destructor Documentation

PixelDetectorManager() [1/3]

InDetDD::PixelDetectorManager::PixelDetectorManager

(

StoreGateSvc *

detStore

)

Constructor.

Definition at line 63 of file PixelDetectorManager.cxx.

     : PixelDetectorManager(detStore, "Pixel", "PixelID"){ }

PixelDetectorManager() [2/3]

InDetDD::PixelDetectorManager::PixelDetectorManager	(	StoreGateSvc *	detStore,
		const std::string &	name,
		const std::string &	pixelIDName )

Constructor with name.

Definition at line 25 of file PixelDetectorManager.cxx.

    : SiDetectorManager(detStore, name),
      m_idHelper(nullptr),
      m_isLogical(false) // Change to true to change the definition of local module corrections
  {
    //
    // Initialized the Identifier helper.
    //
    if (pixelIDName == "PLR_ID") {
      const PLR_ID* plr_idHelper;
      StatusCode sc = detStore->retrieve(plr_idHelper, pixelIDName);
      if (sc.isFailure()) {
        ATH_MSG_ERROR("Could not retrieve PLR_ID helper");
      }
      m_idHelper = plr_idHelper;
      // make the symlink
      sc = detStore->symLink(plr_idHelper, m_idHelper);
      if (sc.isFailure()) {
        ATH_MSG_ERROR("Could not make PLR_ID symlink");
      }
    } else {
      StatusCode sc = detStore->retrieve(m_idHelper, pixelIDName);
      if (sc.isFailure()) {
        ATH_MSG_ERROR("Could not retrieve PixelID helper");
      }
    }
 
    // Initialize the collections.
    if (m_idHelper) {
      m_elementCollection.resize(m_idHelper->wafer_hash_max());
      m_alignableTransforms.resize(m_idHelper->wafer_hash_max());
    }
    //https://its.cern.ch/jira/browse/ATLASRECTS-4886
    InDetDetectorManager::m_suppressWarnings = true;
  }

PixelDetectorManager() [3/3]

InDetDD::PixelDetectorManager::PixelDetectorManager ( const PixelDetectorManager & right )

private

Member Function Documentation

addAlignableTransform() [1/3]

void InDetDD::PixelDetectorManager::addAlignableTransform ( int level, const Identifier & id, GeoAlignableTransform * xf, const GeoVFullPhysVol * child )

virtual

Add alignable transforms.

No access to these, they will be changed by manager:

Definition at line 274 of file PixelDetectorManager.cxx.

  {
    GeoNodePtr<GeoAlignableTransform> tmp_transform(transform);
    if (m_idHelper) {
      if (level == 0) {
        IdentifierHash idHash = m_idHelper->wafer_hash(id);
        if (idHash.is_valid()) {
          m_alignableTransforms[idHash]= std::make_unique<ExtendedAlignableTransform>(transform, child);
        }
      } else {
        // Higher levels are saved in a map. NB the index is level-1 as level=0 is treated above.
        int index = level - FIRST_HIGHER_LEVEL;
        if (index >= static_cast<int>(m_higherAlignableTransforms.size())) m_higherAlignableTransforms.resize(index+1);
        m_higherAlignableTransforms[index][id] = std::make_unique<ExtendedAlignableTransform>(transform, child);
      }
    }
  }

addAlignableTransform() [2/3]

void InDetDD::PixelDetectorManager::addAlignableTransform ( int level, const Identifier & id, GeoAlignableTransform * xf, const GeoVPhysVol * child )

virtual

As above but does a dynamic_cast to GeoVFullPhysVol.

Definition at line 257 of file PixelDetectorManager.cxx.

  {
    GeoNodePtr<GeoAlignableTransform> tmp_transform(transform);
    if (m_idHelper) {
 
      const GeoVFullPhysVol * childFPV = dynamic_cast<const GeoVFullPhysVol *>(child);
      if (!childFPV) {
        ATH_MSG_ERROR("Child of alignable transform is not a full physical volume");
      } else {
        addAlignableTransform (level, id, transform, childFPV);
      }
    }
  }

addAlignableTransform() [3/3]

virtual void InDetDD::PixelDetectorManager::addAlignableTransform ( int level, const Identifier & id, GeoAlignableTransform * xf )

inlineoverridevirtual

Add alignable transforms.

No access to these, they will be changed by manager:

Implements InDetDD::SiDetectorManager.

Definition at line 101 of file PixelDetectorManager.h.

                                                                            {}; // For backward compatibility

addAlignFolderType()

void InDetDD::InDetDetectorManager::addAlignFolderType ( const AlignFolderType alignfolder )

inherited

Definition at line 81 of file InDetDetectorManager.cxx.

    {
        m_alignfoldertype = alignfolder;
    }

addChannel()

void InDetDD::InDetDetectorManager::addChannel ( const std::string & key, int level, FrameType frame )

inherited

Alignment access.

Definition at line 56 of file InDetDetectorManager.cxx.

    {
        std::string frameStr = "other";
        if (frame == InDetDD::global) frameStr = "global";
        if (frame == InDetDD::local) frameStr  = "local";
        ATH_MSG_INFO("Registering alignment channel with key " << key << ", level " << level
                     << ", with frame " << frameStr << ".");
        m_keys[key] = LevelInfo(level, frame);
    }

addDesign()

const SiDetectorDesign * InDetDD::SiDetectorManager::addDesign ( std::unique_ptr< const SiDetectorDesign > && design )

inherited

Access to module design; returns an observer pointer.

Definition at line 134 of file SiDetectorManager.cxx.

    {
        m_designs.push_back(std::move(design));
        return m_designs.back().get();
    }

addDetectorElement()

void InDetDD::PixelDetectorManager::addDetectorElement ( SiDetectorElement * element )

overridevirtual

Add elememts.

Implements InDetDD::SiDetectorManager.

Definition at line 128 of file PixelDetectorManager.cxx.

  {
    IdentifierHash idHash = element->identifyHash();
    if (idHash >=  m_elementCollection.size())
      throw std::runtime_error("PixelDetectorManager: Error adding detector element.");
    m_elementCollection[idHash] = element;
  }

addFolder()

void InDetDD::InDetDetectorManager::addFolder ( const std::string & key )

inherited

Definition at line 66 of file InDetDetectorManager.cxx.

    {
        m_folders.insert(key);
    }

addGlobalFolder()

void InDetDD::InDetDetectorManager::addGlobalFolder ( const std::string & key )

inherited

Definition at line 76 of file InDetDetectorManager.cxx.

    {
        m_globalFolders.insert(key);
    }

addSpecialFolder()

void InDetDD::InDetDetectorManager::addSpecialFolder ( const std::string & key )

inherited

Definition at line 71 of file InDetDetectorManager.cxx.

    {
        m_specialFolders.insert(key);
    }

addTreeTop()

void InDetDD::PixelDetectorManager::addTreeTop

(

const PVConstLink &

vol

)

Add a Tree top:

Definition at line 76 of file PixelDetectorManager.cxx.

                                                             {
    m_volume.push_back(vol);
  }

align() [1/2]

StatusCode InDetDD::InDetDetectorManager::align ( const RawAlignmentObjects & alignObjects, GeoVAlignmentStore * alignStore ) const

inherited

Definition at line 209 of file InDetDetectorManager.cxx.

    {
 
        ATH_MSG_DEBUG("align() called from an alignment CondAlg");
        if (!getIdHelper()) return StatusCode::SUCCESS; // To Do: is it really a success?
 
        bool alignmentChange = false;
        //      const AlignInfo &aligninfo = AlignInfo(m_alignfoldertype);
 
        for(const auto& alignObj : alignObjects) {
            const std::string& key = alignObj.first;
 
            ATH_MSG_DEBUG(" Processing folder  " << key);
 
            if(m_globalFolders.find(key)!=m_globalFolders.end()) {
                try {
                    // New global alignemnts
                    const CondAttrListCollection* obj = static_cast<const CondAttrListCollection*>(alignObj.second);
                    bool status = processGlobalAlignmentContainer(key,obj,alignStore);
                    alignmentChange = (alignmentChange || status);
                } catch(std::runtime_error& err) {
                    // alignments should always exist so we return fatal if we could not process the alignment for this key
                    ATH_MSG_FATAL(err.what());
                    return StatusCode::FAILURE;
                }
            }
            else if(m_folders.find(key)!=m_folders.end()) {
                try {
                    // Regular alignemnts
                    const AlignableTransformContainer* container = static_cast<const AlignableTransformContainer*>(alignObj.second);
                    bool status = processAlignmentContainer(container,alignStore);
                    alignmentChange = (alignmentChange || status);
                } catch(std::runtime_error& err) {
                    // alignments should always exist so we return fatal if we could not process the alignment for this key
                    ATH_MSG_FATAL(err.what());
                    return StatusCode::FAILURE;
                }
            }
            else if(m_specialFolders.find(key)!=m_specialFolders.end()) {
                try {
                  // Detector specific alignments
                  const CondAttrListCollection *obj =
                    static_cast<const CondAttrListCollection*>(alignObj.second);
                  bool status = processSpecialAlignment(key, obj, alignStore);
                  alignmentChange = (alignmentChange || status);
                }
                catch(std::runtime_error& err) {
                  // Should always exist if the folder was requested so we return fatal if
                  // we could not process the alignment for this key
                  ATH_MSG_FATAL(err.what());
                  return StatusCode::FAILURE;
                }
            }
            else {
                // Should not be any other keys specified in raw alignment object.
                ATH_MSG_ERROR("Unrecognized folder name "<<key<<". Expected names are:");
                for (const std::string& out:m_globalFolders) ATH_MSG_ERROR("--"<<out);
                for (const std::string& out:m_folders) ATH_MSG_ERROR("--"<<out);
                for (const std::string&  out:m_specialFolders) ATH_MSG_ERROR("--"<<out);
 
                return StatusCode::RECOVERABLE;
            }
        }
        return StatusCode::SUCCESS;
    }

align() [2/2]

StatusCode InDetDD::InDetDetectorManager::align ( IOVSVC_CALLBACK_ARGS )

inherited

getBarrelLayers()

unsigned int InDetDD::PixelDetectorManager::getBarrelLayers

(

)

const

Layers built.

Deprecated. Recommend to access via numerology (see SiDetectorManager) Warning: decodes in "human-readable", e.g. using 3 out of 3 pixel layers yiels 111, the laset two out of 3 in 110, etc.

Definition at line 172 of file PixelDetectorManager.cxx.

  {
    unsigned int barrelLayers = 0;
    for (int i = numerology().numLayers() - 1; i >= 0; --i) {
      barrelLayers *= 10;
      barrelLayers += numerology().useLayer(i);
    }
    return barrelLayers;
  }

getDesign()

const SiDetectorDesign * InDetDD::SiDetectorManager::getDesign ( int i ) const

inherited

Definition at line 146 of file SiDetectorManager.cxx.

    {
        return m_designs[i].get();
    }

getDetectorElement() [1/3]

const SiDetectorElement * InDetDD::PixelDetectorManager::getDetectorElement ( const Identifier & id ) const

overridevirtual

access to individual elements : via Identifier

Implements InDetDD::SiDetectorManager.

Definition at line 80 of file PixelDetectorManager.cxx.

  {
  // NB the id helpers implementation for getting a hash is not optimal.
  // Essentially does a binary search.
  // Make sure it is a wafer Id
    Identifier waferId = m_idHelper->wafer_id(id);
    IdentifierHash idHash = m_idHelper->wafer_hash(waferId);
    if (idHash.is_valid()) {
      return m_elementCollection[idHash];
    } else {
      return nullptr;
    }
  }

getDetectorElement() [2/3]

const SiDetectorElement * InDetDD::PixelDetectorManager::getDetectorElement ( const IdentifierHash & idHash ) const

overridevirtual

access to individual elements : via IdentifierHash

Implements InDetDD::SiDetectorManager.

Definition at line 94 of file PixelDetectorManager.cxx.

  {
    return m_elementCollection[idHash];
  }

getDetectorElement() [3/3]

const SiDetectorElement * InDetDD::PixelDetectorManager::getDetectorElement	(	int	barrel_endcap,
		int	layer_wheel,
		int	phi_module,
		int	eta_module ) const

access to individual elements : via element identification

Definition at line 99 of file PixelDetectorManager.cxx.

  {
    return getDetectorElement(m_idHelper->wafer_id(barrel_endcap, layer_wheel, phi_module, eta_module));
  }

getDetectorElementBegin() [1/2]

SiDetectorElementCollection::const_iterator InDetDD::PixelDetectorManager::getDetectorElementBegin ( ) const

overridevirtual

Implements InDetDD::SiDetectorManager.

Definition at line 109 of file PixelDetectorManager.cxx.

  {
    return m_elementCollection.begin();
  }

getDetectorElementBegin() [2/2]

SiDetectorElementCollection::iterator InDetDD::PixelDetectorManager::getDetectorElementBegin ( )

overridevirtual

Implements InDetDD::SiDetectorManager.

Definition at line 118 of file PixelDetectorManager.cxx.

  {
    return m_elementCollection.begin();
  }

getDetectorElementCollection()

const SiDetectorElementCollection * InDetDD::PixelDetectorManager::getDetectorElementCollection ( ) const

overridevirtual

access to whole collection via Iterators

Implements InDetDD::SiDetectorManager.

Definition at line 104 of file PixelDetectorManager.cxx.

  {
    return &m_elementCollection;
  }

getDetectorElementEnd() [1/2]

SiDetectorElementCollection::const_iterator InDetDD::PixelDetectorManager::getDetectorElementEnd ( ) const

overridevirtual

Implements InDetDD::SiDetectorManager.

Definition at line 114 of file PixelDetectorManager.cxx.

  {
    return m_elementCollection.end();
  }

getDetectorElementEnd() [2/2]

SiDetectorElementCollection::iterator InDetDD::PixelDetectorManager::getDetectorElementEnd ( )

overridevirtual

Implements InDetDD::SiDetectorManager.

Definition at line 123 of file PixelDetectorManager.cxx.

  {
    return m_elementCollection.end();
  }

getEndcapLayers()

unsigned int InDetDD::PixelDetectorManager::getEndcapLayers

(

)

const

Layers built.

Deprecated. Recommend to access via numerology (see SiDetectorManager) Warning: decodes in "human-readable", e.g. using 3 out of 3 pixel layers yiels 111, the laset two out of 3 in 110, etc.

Definition at line 182 of file PixelDetectorManager.cxx.

  {
    unsigned int endcapLayers = 0;
    for (int i = numerology().numDisks() - 1; i >= 0; --i) {
      endcapLayers *= 10;
      endcapLayers += numerology().useDisk(i);
    }
    return endcapLayers;
  }

getIdHelper()

const PixelID * InDetDD::PixelDetectorManager::getIdHelper ( ) const

overrideprivatevirtual

Implements InDetDD::InDetDetectorManager.

Definition at line 167 of file PixelDetectorManager.cxx.

  {
    return m_idHelper;
  }

getLayout()

const std::string & InDetDD::InDetDetectorManager::getLayout ( ) const

inherited

Definition at line 38 of file InDetDetectorManager.cxx.

    {
        return m_version.layout();
    }

getLevel()

const InDetDetectorManager::LevelInfo & InDetDD::InDetDetectorManager::getLevel ( const std::string & key ) const

privateinherited

Retrieve level information.

Definition at line 87 of file InDetDetectorManager.cxx.

    {
        std::map<std::string, LevelInfo>::const_iterator iter;
        iter = m_keys.find(key);
        if (iter == m_keys.end()) return s_invalidLevel;
        return iter->second;
    }

getNumTreeTops()

unsigned int InDetDD::PixelDetectorManager::getNumTreeTops ( ) const

overridevirtual

Access to raw geometry:

Definition at line 66 of file PixelDetectorManager.cxx.

  {
    return m_volume.size();
  }

getPixelDesign()

const PixelModuleDesign * InDetDD::PixelDetectorManager::getPixelDesign ( int i ) const

virtual

Access to module design, Casts to PixelModuleDesign.

Definition at line 302 of file PixelDetectorManager.cxx.

  {
    return dynamic_cast<const PixelModuleDesign *>(getDesign(i));
  }

getTreeTop()

PVConstLink InDetDD::PixelDetectorManager::getTreeTop ( unsigned int i ) const

overridevirtual

Definition at line 71 of file PixelDetectorManager.cxx.

  {
    return m_volume[i];
  }

getVersion()

const Version & InDetDD::InDetDetectorManager::getVersion ( ) const

inherited

Get version information.

Definition at line 33 of file InDetDetectorManager.cxx.

    {
        return m_version;
    }

identifierBelongs()

bool InDetDD::PixelDetectorManager::identifierBelongs ( const Identifier & id ) const

overridevirtual

Check identifier is for this detector.

Implements InDetDD::InDetDetectorManager.

Definition at line 295 of file PixelDetectorManager.cxx.

  {
    return getIdHelper()->is_pixel(id);
  }

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

initNeighbours()

void InDetDD::PixelDetectorManager::initNeighbours ( )

overridevirtual

Initialize the neighbours.

This can only be done when all elements are built.

Implements InDetDD::SiDetectorManager.

Definition at line 136 of file PixelDetectorManager.cxx.

  {
    SiDetectorElementCollection::iterator iter;
 
    // Loop over all elements and set the neighbours
    for (iter = m_elementCollection.begin(); iter != m_elementCollection.end(); ++iter){
 
      SiDetectorElement * element = *iter;
      if (element) {
 
        IdentifierHash idHash = element->identifyHash();
        IdentifierHash idHashOther;
 
        int result;
        // If no neighbour, result != 0 in which case we leave neighbour as null
        result = m_idHelper->get_next_in_eta(idHash, idHashOther);
        if (result==0) element->setNextInEta(m_elementCollection[idHashOther]);
 
        result = m_idHelper->get_prev_in_eta(idHash, idHashOther);
        if (result==0) element->setPrevInEta(m_elementCollection[idHashOther]);
 
        result = m_idHelper->get_next_in_phi(idHash, idHashOther);
        if (result==0) element->setNextInPhi(m_elementCollection[idHashOther]);
 
        result = m_idHelper->get_prev_in_phi(idHash, idHashOther);
        if (result==0) element->setPrevInPhi(m_elementCollection[idHashOther]);
      }
    }
  }

invalidateAll()

void InDetDD::SiDetectorManager::invalidateAll ( )

virtualinherited

Invalidate cache for all detector elements.

Implements InDetDD::InDetDetectorManager.

Definition at line 39 of file SiDetectorManager.cxx.

    {
        for (SiDetectorElementCollection::iterator element_iter = getDetectorElementBegin();
        element_iter != getDetectorElementEnd();
        ++element_iter) {
 
            if (*element_iter) {
                (*element_iter)->invalidate();
            }
        }
    }

isPixel()

virtual bool InDetDD::PixelDetectorManager::isPixel ( ) const

inlineoverridevirtual

Methods to query which manager we have.

Implements InDetDD::SiDetectorManager.

Definition at line 110 of file PixelDetectorManager.h.

{return true;}

isSCT()

bool InDetDD::SiDetectorManager::isSCT ( ) const

inlineinherited

Definition at line 100 of file SiDetectorManager.h.

{return !isPixel();}

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 163 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 178 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

numDesigns()

int InDetDD::SiDetectorManager::numDesigns ( ) const

inherited

Definition at line 140 of file SiDetectorManager.cxx.

    {
        return m_designs.size();
    }

numerology() [1/2]

SiNumerology & InDetDD::SiDetectorManager::numerology ( )

inlineinherited

Definition at line 129 of file SiDetectorManager.h.

{return m_numerology;}

numerology() [2/2]

const SiNumerology & InDetDD::SiDetectorManager::numerology ( ) const

inlineinherited

Access Numerology.

Definition at line 128 of file SiDetectorManager.h.

{return m_numerology;}

operator=()

const PixelDetectorManager & InDetDD::PixelDetectorManager::operator= ( const PixelDetectorManager & right )

private

prevent copy and assignment

printTransform()

void InDetDD::PixelDetectorManager::printTransform ( MSG::Level level, const Amg::Transform3D & tr ) const

private

Definition at line 663 of file PixelDetectorManager.cxx.

  {
    ATH_MSG_LVL_NOCHK(level, " - translation: " << tr.translation().x() << "  " <<tr.translation().y() << "  " <<tr.translation().z());
    ATH_MSG_LVL_NOCHK(level, " - rotation:");
    ATH_MSG_LVL_NOCHK(level, "    " << tr(0,0) << "  " << tr(0,1) << "  " << tr(0,2));
    ATH_MSG_LVL_NOCHK(level, "    " << tr(1,0) << "  " << tr(1,1) << "  " << tr(1,2));
    ATH_MSG_LVL_NOCHK(level, "    " << tr(2,0) << "  " << tr(2,1) << "  " << tr(2,2));
 }

processAlignmentContainer() [1/2]

bool InDetDD::InDetDetectorManager::processAlignmentContainer ( const AlignableTransformContainer * container, GeoVAlignmentStore * alignStore ) const

privateinherited

Definition at line 303 of file InDetDetectorManager.cxx.

    {
        bool alignmentChange = false;
 
        // Check if container is empty - this can occur if it is an invalid IOV.
        if (container->empty()) {
            ATH_MSG_ERROR("AlignableTransformContainer "
                          << " is empty. Probably due to out of range IOV"); // To Do: add key to this printout for making it more informative
            // This should not occur in normal situations so we force job to abort.
            throw std::runtime_error("Unable to apply Inner Detector alignments.");
        }
        // loop over all the AlignableTransform objects in the collection
        // use only the last ones.
        // /Indet/AlignL3/SCTEA9 appear repeatedly in tags of the /Indet/AlignL3 folder
        std::map<const std::string, const AlignableTransform*> stringToTransform;
        for (const auto *pat : *container) {
            stringToTransform[pat->tag()] = pat;
        }
        for (const std::pair<const std::string, const AlignableTransform*>& value: stringToTransform) {
            bool status = processKey(value.first, value.second, alignStore);
            alignmentChange = (alignmentChange || status);
        }
        return alignmentChange;
    }

processAlignmentContainer() [2/2]

bool InDetDD::InDetDetectorManager::processAlignmentContainer ( const std::string & key ) const

privateinherited

return align folder string to use

Process the alignment container, calls processKey

Definition at line 275 of file InDetDetectorManager.cxx.

    {
        bool alignmentChange = false;
 
        ATH_MSG_DEBUG("Dealing with key as container");
        const AlignableTransformContainer* container;
        if (StatusCode::SUCCESS!=m_detStore->retrieve(container, key)) {
            ATH_MSG_ERROR("Cannot find AlignableTransformContainer for key "
                          << key << " - no misalignment");
            // This should not occur in normal situations so we force job to abort.
            throw std::runtime_error("Unable to apply Inner Detector alignments");
        }
        // Check if container is empty - this can occur if it is an invalid IOV.
        if (container->empty()) {
            ATH_MSG_ERROR("AlignableTransformContainer for key "
                          << key << " is empty. Probably due to out of range IOV");
            // This should not occur in normal situations so we force job to abort.
            throw std::runtime_error("Unable to apply Inner Detector alignments.");
        }
        // loop over all the AlignableTransform objects in the collection
        for (const auto *pat : *container) {
 
            bool status = processKey(pat->tag(),pat);
            alignmentChange = (alignmentChange || status);
        }
        return alignmentChange;
    }

processGlobalAlignment()

bool InDetDD::PixelDetectorManager::processGlobalAlignment ( const std::string & key, int level, FrameType frame, const CondAttrListCollection * obj, GeoVAlignmentStore * alignStore ) const

overridevirtual

Process new global DB folders for L1 and L2.

Reimplemented from InDetDD::InDetDetectorManager.

Definition at line 577 of file PixelDetectorManager.cxx.

  {
 
    bool alignmentChange = false;
 
    ATH_MSG_INFO("Processing new global alignment containers with key " << key
                 << " in the " << frame << " frame at level " << level);
 
    Identifier ident=Identifier();
    const CondAttrListCollection* atrlistcol=obj;
    //cppcheck-suppress nullPointerRedundantCheck
    if (!atrlistcol) {
       ATH_MSG_INFO("Read alignment from detector store with key " << key);
       if (StatusCode::SUCCESS!=m_detStore->retrieve(atrlistcol,key)) {
          ATH_MSG_WARNING("Cannot find new global align Container for key "
                          << key << " - no new global alignment");
          return alignmentChange;
       }
    }
    // Avoid cppcheck warning.
    if (!atrlistcol) {
      return alignmentChange;
    }
    {
      // loop over objects in collection
      //cppcheck-suppress nullPointerRedundantCheck
      for (CondAttrListCollection::const_iterator citr=atrlistcol->begin(); citr!=atrlistcol->end();++citr) {
        const coral::AttributeList& atrlist=citr->second;
        // We are in the Pixel manager, therefore ignore all that is not Pixel Identifier
        if (atrlist["det"].data<int>()!=1) continue;
 
        ident = getIdHelper()->wafer_id(atrlist["bec"].data<int>(),
                atrlist["layer"].data<int>(),
                atrlist["ring"].data<int>(),
                atrlist["sector"].data<int>());
 
        // Make sure we have valid HashID (This ONLY works here as the 0-modules exist)
        // Precaution which does not work for e.g. SCT
        if (!(getIdHelper()->wafer_hash(ident)).is_valid()){
          ATH_MSG_WARNING("Invalid HashID for identifier " << getIdHelper()->show_to_string(ident));
          ATH_MSG_WARNING("No global alignment corrections can be applied for invalid HashID's - exiting");
          return false;
        }
 
        // construct new transform
        // Order of rotations is defined as around z, then y, then x.
        Amg::Translation3D  newtranslation(atrlist["Tx"].data<float>(),atrlist["Ty"].data<float>(),atrlist["Tz"].data<float>());
        Amg::Transform3D newtrans = newtranslation * Amg::RotationMatrix3D::Identity();
        newtrans *= Amg::AngleAxis3D(atrlist["Rz"].data<float>()*CLHEP::mrad, Amg::Vector3D(0.,0.,1.));
        newtrans *= Amg::AngleAxis3D(atrlist["Ry"].data<float>()*CLHEP::mrad, Amg::Vector3D(0.,1.,0.));
        newtrans *= Amg::AngleAxis3D(atrlist["Rx"].data<float>()*CLHEP::mrad, Amg::Vector3D(1.,0.,0.));
 
        ATH_MSG_DEBUG("New global DB -- channel: " << citr->first
                      << " ,det: "    << atrlist["det"].data<int>()
                      << " ,bec: "    << atrlist["bec"].data<int>()
                      << " ,layer: "  << atrlist["layer"].data<int>()
                      << " ,ring: "   << atrlist["ring"].data<int>()
                      << " ,sector: " << atrlist["sector"].data<int>()
                      << " ,Tx: "     << atrlist["Tx"].data<float>()
                      << " ,Ty: "     << atrlist["Ty"].data<float>()
                      << " ,Tz: "     << atrlist["Tz"].data<float>()
                      << " ,Rx: "     << atrlist["Rx"].data<float>()
                      << " ,Ry: "     << atrlist["Ry"].data<float>()
                      << " ,Rz: "     << atrlist["Rz"].data<float>());
 
        // Set the new transform; Will replace existing one with updated transform
        bool status = setAlignableTransformDelta(level,
                                                 ident,
                                                 newtrans,
                                                 frame,
                                                 alignStore);
 
        if (!status) {
          ATH_MSG_DEBUG("Cannot set AlignableTransform for identifier."
                        << getIdHelper()->show_to_string(ident)
                        << " at level " << level << " for new global DB");
        }
 
        alignmentChange = (alignmentChange || status);
      }
    }
    return alignmentChange;
  }

processGlobalAlignmentContainer()

bool InDetDD::InDetDetectorManager::processGlobalAlignmentContainer ( const std::string & key, const CondAttrListCollection * obj = nullptr, GeoVAlignmentStore * alignStore = nullptr ) const

privateinherited

Definition at line 391 of file InDetDetectorManager.cxx.

    {
      bool alignmentChange = false;
 
      ATH_MSG_DEBUG("processing GlobalAlignmentContainer with key:  " << key);
      // From the key determine what level in hierarchy we are dealing with.
      // returns -1 if unrecognized.
      const LevelInfo & levelInfo = getLevel(key);
      if (levelInfo.isValid()) {
          ATH_MSG_VERBOSE("Processing channel: " << key);
      } else {
          ATH_MSG_DEBUG("Channel " << key << " not registered in this manager");
      }
      // return silently if unrecognised - this can happen in container mode
      // when a single container holds transforms for both pixel and SCT
      if (!levelInfo.isValid() ) return false;
 
      // Within detector specific code
      bool status = processGlobalAlignment(key, levelInfo.level(), levelInfo.frame(), obj, alignStore);
 
      alignmentChange = (alignmentChange || status);
 
      return alignmentChange;
 
    }

processKey()

bool InDetDD::InDetDetectorManager::processKey ( const std::string & key, const AlignableTransform * transformCollection, GeoVAlignmentStore * alignStore = nullptr ) const

privateinherited

Called by processAlignmentContainer, applies only one key on the transform Collections.

Definition at line 328 of file InDetDetectorManager.cxx.

    {
        bool alignmentChange = false;
 
        // From the key determine what level in hierarchy we are dealing with.
        // returns -1 if unrecognized.
        const LevelInfo & levelInfo = getLevel(key);
        if (levelInfo.isValid()) {
            ATH_MSG_VERBOSE("Processing channel: " << key);
        } else {
            ATH_MSG_DEBUG("Channel " << key << " not registered in this manager");
        }
        // return silently if unrecognised - this can happen in container mode
        // when a single container holds transforms for both pixel and SCT
        if (!levelInfo.isValid() ) return false;
 
        //Loop over the effected nodes.
        for (AlignableTransform::AlignTransMem_citr trans_iter = transformCollection->begin();
             trans_iter != transformCollection->end();
             ++trans_iter) {
            ATH_MSG_DEBUG( "Get alignment for identifier "
                           << getIdHelper()->show_to_string(trans_iter->identify())
                           << " at level " << levelInfo.level());
 
            // The delta in the conditions DB is not necessarily the same as what is needed in the
            // alignable transform. At the moment we support global frame, local frame or an alternative frame
            // The setAlignableTransformDelta method takes care of this correction - this is CLHEP <--> Amg interfaced
            bool status = setAlignableTransformDelta(levelInfo.level(),
                                                     trans_iter->identify(),
                                                     Amg::CLHEPTransformToEigen(trans_iter->transform()),
                                                     levelInfo.frame(),
                                                     alignStore);
 
            alignmentChange = (alignmentChange || status);
 
            if (!status) {
                if (!identifierBelongs(trans_iter->identify())) {
                    // Its probably OK. Eg /Indet/Align/ID contains alse pixel and sct ids.
                    ATH_MSG_DEBUG("Cannot set AlignableTransform for identifier."
                                  << " Probably OK if its /Indet/Align/ID folder. "
                                  << getIdHelper()->show_to_string(trans_iter->identify())
                                  << " at level " << levelInfo.level());
                } else {
                    if (m_suppressWarnings) {
                        ATH_MSG_DEBUG("WARNING: Cannot set AlignableTransform for identifier  "
                                      << getIdHelper()->show_to_string(trans_iter->identify())
                                      << " at level " << levelInfo.level());
                    } else {
                        ATH_MSG_WARNING("Cannot set AlignableTransform for identifier  "
                                        << getIdHelper()->show_to_string(trans_iter->identify())
                                        << " at level " << levelInfo.level());
                        ATH_MSG_WARNING("Subsequent WARNINGS will be printed at DEBUG level.");
                        m_suppressWarnings = true;
                    }
                }
            }
        }
        return alignmentChange;
    }

processSpecialAlignment() [1/2]

bool InDetDD::PixelDetectorManager::processSpecialAlignment ( const std::string & key, InDetDD::AlignFolderType alignfolderType )

overridevirtual

Process new IBLDist DB folder.

Matthias D. (Oct. 2016): The idea of this first implementation is to get the AlignTransforms for IBL modules from their folder; Calculate from the new DB entries the Tx displacement for this module; Make a simple transfrom in local frame; Add the Tx to the aligntransform; Apply this new transform with the setAlignableTransformDelta() function; This is non-optimal as it simply overrides the previous Delta (these are included in total Delta); Possibly we could extend this in future to add tweak() functionality to GeoAlignableTransform?

Verbose level debug section for transforms

End of verbose level debug section

Implements InDetDD::InDetDetectorManager.

Definition at line 311 of file PixelDetectorManager.cxx.

  {
 
  bool alignmentChange = false;
 
  std::string alignfolder;
  if (alignfolderType == InDetDD::static_run1) alignfolder = "/Indet/Align";
  if (alignfolderType == InDetDD::timedependent_run2) alignfolder = "/Indet/AlignL3";
 
  ATH_MSG_INFO("Processing IBLDist alignment container with key (" << key << ")"
               << "and alignment folder pointing to " << alignfolder);
 
  int nstaves = 0;
  if (numerology().numPhiModulesForLayer(0)<14) nstaves = 14;
  else nstaves = numerology().numPhiModulesForLayer(0);
 
  std::vector<float> ibldist;
  std::vector<float> iblbaseline;
  ibldist.resize(nstaves);
  iblbaseline.resize(nstaves);
 
  const CondAttrListCollection* atrlistcol=nullptr;
  if (StatusCode::SUCCESS==m_detStore->retrieve(atrlistcol,key)) {
    // loop over objects in collection
    for (CondAttrListCollection::const_iterator citr=atrlistcol->begin(); citr!=atrlistcol->end();++citr) {
      const coral::AttributeList& atrlist=citr->second;
      ibldist[atrlist["stave"].data<int>()] = atrlist["mag"].data<float>();
      iblbaseline[atrlist["stave"].data<int>()] = atrlist["base"].data<float>();
 
      ATH_MSG_VERBOSE("IBLDist DB -- channel: " << citr->first
                      << " ,stave: " << atrlist["stave"].data<int>()
                      << " ,mag: " << atrlist["mag"].data<float>()
                      << " ,base: " << atrlist["base"].data<float>());
    }
  }
  else {
    ATH_MSG_WARNING("Cannot find IBLDist Container for key "
                    << key << " - no IBL bowing alignment");
    return alignmentChange;
  }

 
  const AlignableTransformContainer* container;
  if (StatusCode::SUCCESS!=m_detStore->retrieve(container, alignfolder)) {
    ATH_MSG_ERROR("Cannot find AlignableTransformContainer for key "
                  << key << " - no misalignment");
    // This should not occur in normal situations so we force job to abort.
    throw std::runtime_error("Unable to apply Inner Detector alignments");
  }
  // Check if container is empty - this can occur if it is an invalid IOV.
  if (container->empty()) {
    ATH_MSG_ERROR("AlignableTransformContainer for key "
                  << key << " is empty. Probably due to out of range IOV");
    // This should not occur in normal situations so we force job to abort.
    throw std::runtime_error("Unable to apply Inner Detector alignments.");
  }
  // loop over all the AlignableTransform objects in the collection
  std::string IBLalignfolder = alignfolder;
  IBLalignfolder.append("/PIXB1");// "/Indet/Align/PIXB1"
  for (const auto *pat : *container)
  {
    if (!( pat->tag()==IBLalignfolder &&
     numerology().numPhiModulesForLayer(0)==14 &&
     numerology().numLayers()==4) ){  // hard-coded to IBL for now; no other geometry should really apply this!
      ATH_MSG_DEBUG("IBLDist; ignoring collections " << pat->tag());
    }
    else{
      const AlignableTransform* transformCollection = pat;
      for (AlignableTransform::AlignTransMem_citr trans_iter = transformCollection->begin();
         trans_iter != transformCollection->end();
         ++trans_iter)
      {
        ATH_MSG_DEBUG("IBLDist alignment for identifier "
                      << getIdHelper()->show_to_string(trans_iter->identify()));
 
        IdentifierHash idHash = getIdHelper()->wafer_hash(trans_iter->identify());
        if (!idHash.is_valid()){
          ATH_MSG_WARNING("Invalid HashID for identifier "
                          << getIdHelper()->show_to_string(trans_iter->identify()));
          ATH_MSG_WARNING("No IBLDist corrections can be applied for invalid HashID's - exiting");
          return false;
        }
        const SiDetectorElement * sielem = m_elementCollection[idHash];
        //This should work as Bowing is in L3 frame, i.e. local module frame
        Amg::Vector3D center = sielem->defTransform() * Amg::Vector3D{0, 0, 0};
        double z = center[2];
        const double  y0y0  = 366.5*366.5;
 
        double bowx = ibldist[getIdHelper()->phi_module(trans_iter->identify())] * ( z*z - y0y0 ) / y0y0;
        double basex= iblbaseline[getIdHelper()->phi_module(trans_iter->identify())];
        // This is in the module frame, as bowing corrections are directly L3
 
        ATH_MSG_DEBUG("Total IBL-module Tx shift (baseline+bowing): " << basex+bowx);
        if ( (basex+bowx)==0 ) continue; // make sure we ignore NULL corrections
 
        Amg::Transform3D shift = Amg::Translation3D(basex+bowx,0,0) * Amg::RotationMatrix3D::Identity();
 
        const AlignableTransform* cpat = pat;
        AlignableTransform::AlignTransMem_citr this_trans=cpat->findIdent(trans_iter->identify());
        HepGeom::Transform3D newtrans = Amg::EigenTransformToCLHEP(shift)*this_trans->transform();

        if (msgLvl(MSG::VERBOSE)) {
          ATH_MSG_LVL_NOCHK(MSG::VERBOSE, "Bowing Transformation only:");
          printTransform(MSG::VERBOSE, shift);
          ATH_MSG_LVL_NOCHK(MSG::VERBOSE, "Original alignable Transformation from StoreGate:");
          printTransform(MSG::VERBOSE, Amg::CLHEPTransformToEigen(this_trans->transform()));
          ATH_MSG_LVL_NOCHK(MSG::VERBOSE, "Final mModified Transformation:");
          printTransform(MSG::VERBOSE, Amg::CLHEPTransformToEigen(newtrans));
        }
 
        // Set the new transform; Will replace existing one with updated transform
        bool status = setAlignableTransformDelta(0,
                                                 trans_iter->identify(),
                                                 Amg::CLHEPTransformToEigen(newtrans),
                                                 InDetDD::local,
                                                 nullptr);
 
        if (!status) {
          ATH_MSG_DEBUG("Cannot set AlignableTransform for identifier."
                        << getIdHelper()->show_to_string(trans_iter->identify())
                        << " at level 0 for IBLDist bowing deformation");
        }
        alignmentChange = (alignmentChange || status);
      }
    }
  }
 
  return alignmentChange;
  }

processSpecialAlignment() [2/2]

bool InDetDD::PixelDetectorManager::processSpecialAlignment ( const std::string & key, const CondAttrListCollection * obj = nullptr, GeoVAlignmentStore * alignStore = nullptr ) const

overridevirtual

Paul Gessinger (Jun 2019): For MT, we won't retrieve the underlying actual alignable transform container, since that's hard to come by. We loop over all detector elements, check if they match the detector elements we have IBL dist info for, retrieve the actual transform from the provided GeoAlignmentStore, and re-set it to the corrected value including IBL bowing.

The calculation of the bowing is taken from the non-MT implementation above.

Implements InDetDD::InDetDetectorManager.

Definition at line 453 of file PixelDetectorManager.cxx.

                                                                                           {
    bool alignmentChange = false;
 
    ATH_MSG_INFO("Processing IBLDist alignment container with key " << key);
    if(numerology().numLayers() != 4) {
      // this doesn't appear to be Run 2, i.e. the IBL isn't there. Bailing
      return alignmentChange;
    }
 
    int nstaves = 0;
    if (numerology().numPhiModulesForLayer(0) < 14)
      nstaves = 14;
    else
      nstaves = numerology().numPhiModulesForLayer(0);
 
    std::vector<float> ibldist;
    std::vector<float> iblbaseline;
    ibldist.resize(nstaves);
    iblbaseline.resize(nstaves);
 
    // loop over objects in collection
    for (CondAttrListCollection::const_iterator citr = obj->begin();
        citr != obj->end(); ++citr) {
 
      const coral::AttributeList &atrlist = citr->second;
      ibldist[atrlist["stave"].data<int>()] = atrlist["mag"].data<float>();
      iblbaseline[atrlist["stave"].data<int>()] = atrlist["base"].data<float>();
 
      ATH_MSG_VERBOSE("IBLDist DB -- channel: " << citr->first
                      << ", stave: " << atrlist["stave"].data<int>()
                      << ", mag: " << atrlist["mag"].data<float>()
                      << ", base: " << atrlist["base"].data<float>());
    }

 
 
    for(const auto* detElem : *getDetectorElementCollection()) {
      if(!detElem->isInnermostPixelLayer()) {
        // skip non-IBL elements. This only works if the innermost pixel layer is in fact
        // the IBL. That should be the case for Run2 until replacement of the ID.
        continue;
      }
      std::string repr = getIdHelper()->show_to_string(detElem->identify());
      ATH_MSG_DEBUG("IBLDist alignment for identifier " << repr);
 
 
      IdentifierHash idHash = getIdHelper()->wafer_hash(detElem->identify());
      if (!idHash.is_valid()) {
        ATH_MSG_WARNING("Invalid HashID for identifier " << repr);
        ATH_MSG_WARNING("No IBLDist corrections can be applied for "
                        "invalid HashID's - exiting ");
        return false;
      }
 
      // extract the stave number
      int stave = getIdHelper()->phi_module(detElem->identify());
 
      Amg::Vector3D center = detElem->defTransform() * Amg::Vector3D{0, 0, 0};
      double z = center[2];
      const double y0y0 = 366.5 * 366.5;
 
      double bowx = ibldist[stave] * (z * z - y0y0) / y0y0;
      double basex = iblbaseline[stave];
      // This is in the module frame, as bowing corrections are directly L3
 
      ATH_MSG_DEBUG("Total IBL-module Tx shift (baseline+bowing): " << basex + bowx);
 
      if ((basex + bowx) == 0) {
        continue; // make sure we ignore NULL corrections
      }
 
      Amg::Transform3D shift = Amg::Translation3D(basex + bowx, 0, 0) *
                               Amg::RotationMatrix3D::Identity();
 
      // now we need to get the original alignment delta to apply this additional
      // shift to
      ExtendedAlignableTransform* eat = m_alignableTransforms[idHash].get();
      const GeoTrf::Transform3D* currentDelta = alignStore->getDelta(eat->alignableTransform());
      if (currentDelta == nullptr) {
        ATH_MSG_ERROR("Have IBL Dist for element which does not have an alignment delta."
                      << " This indicates inconsistent alignment data");
        return false;
      }
 
      HepGeom::Transform3D recoToHitTransform = detElem->recoToHitTransform();
      Amg::Transform3D r2h = Amg::CLHEPTransformToEigen(recoToHitTransform);
 
      ATH_MSG_VERBOSE("Previous delta for " << repr << ":\n" << currentDelta->matrix());
      ATH_MSG_VERBOSE("Bowing-only delta for " << repr << ":\n" << shift.matrix());
      Amg::Transform3D newDelta = shift * r2h.inverse() * (*currentDelta) * r2h;
 
      // We can probably just write it back to the geo alignment store.
      // The IBL bowing is always at the module level, and that's the delta that
      // we retrieved from the geo alignment store.
 
      bool status = setAlignableTransformDelta(0,
                                               detElem->identify(),
                                               newDelta,
                                               InDetDD::local,
                                               alignStore);
 
 
      ATH_MSG_VERBOSE("New delta for " << repr << ":\n"
                      << alignStore->getDelta(eat->alignableTransform())->matrix());
 
      alignmentChange |= status;
    }
 
    return alignmentChange;
  }

setAlignableTransformDelta()

bool InDetDD::PixelDetectorManager::setAlignableTransformDelta ( int level, const Identifier & id, const Amg::Transform3D & delta, FrameType frame, GeoVAlignmentStore * alignStore ) const

overrideprivatevirtual

implements the main alignment update for delta transforms in different frames, it translates into the LocalDelta or GlobalDelta function of SiDetectorManager

Implements InDetDD::SiDetectorManager.

Definition at line 193 of file PixelDetectorManager.cxx.

  {
 
    if (level == 0) { // At the element level - local shift
 
      // We retrieve it via a hashId.
      IdentifierHash idHash = m_idHelper->wafer_hash(id);
      if (!idHash.is_valid()) return false;
 
      if (frame == InDetDD::global) {
 
        return setAlignableTransformGlobalDelta(m_alignableTransforms[idHash].get(), delta, alignStore);
 
      } else if (frame == InDetDD::local) {
 
        const SiDetectorElement * element =  m_elementCollection[idHash];
        if (!element) return false;
 
        // Its a local transform
        //See header file for definition of m_isLogical
        if( m_isLogical ){
          // Ensure cache is up to date and use the alignment corrected local to
          // global transform
          element->updateCache();
          return setAlignableTransformLocalDelta(
              m_alignableTransforms[idHash].get(), element->transform(), delta,
              alignStore);
        } else
          //Use default local to global transform
          return setAlignableTransformLocalDelta(
              m_alignableTransforms[idHash].get(), element->defTransform(),
              delta, alignStore);
      } else {
        // other not supported
        ATH_MSG_WARNING("Frames other than global or local are not supported.");
        return false;
      }
 
    } else { // higher level
 
      if (frame != InDetDD::global) {
        ATH_MSG_WARNING("Non global shift at higher levels is not possible.");
        return false;
      }
 
      int index = level - FIRST_HIGHER_LEVEL; // level 0 is treated separately.
      if (index >=  static_cast<int>(m_higherAlignableTransforms.size())) return false;
 
      // We retrieve it from a map.
      AlignableTransformMap::const_iterator iter;
      iter = m_higherAlignableTransforms[index].find(id);
      if (iter == m_higherAlignableTransforms[index].end()) return false;
 
      // Its a global transform
      return setAlignableTransformGlobalDelta((iter->second).get(), delta, alignStore);
    }
 
  }

setAlignableTransformGlobalDelta()

bool InDetDD::SiDetectorManager::setAlignableTransformGlobalDelta ( ExtendedAlignableTransform * extXF, const Amg::Transform3D & delta, GeoVAlignmentStore * alignStore = nullptr )

staticinherited

Helper method to set delta transform from a global delta - Amg interface.

Definition at line 101 of file SiDetectorManager.cxx.

                                                                                             {
        // ATTENTION -------------------------------------------------------- (A.S.)
        // CLHEP < -- > AMG interface method
 
        // Sets the alignable transform delta when the supplied delta is in the global frame.
 
        // If the default transform down to the alignable transform is
        // T = A*B*C
        // and the alignable transform is C with delta c and the delta in the global frame is g, then
        // A*B*C*c = g*A*B*C
        // T*c = g*T
        //  c = T.inverse() * g * T
 
        // To get default transform up and including the alignable transform,
        // we assume the next volume down is a fullphys volume and so its
        // transform is the transform we want (ie T=A*B*C in the above example).
 
 
        if (!extXF) return false;
 
        const GeoVFullPhysVol * child = extXF->child();
        if (child && extXF->alignableTransform()) {
            // do the calculation in CLHEP
            const GeoTrf::Transform3D& transform = child->getDefAbsoluteTransform(alignStore);
            extXF->alignableTransform()->setDelta(transform.inverse() * delta * transform, alignStore);
            return true;
        } else {
            return false;
        }
    }

setAlignableTransformLocalDelta()

bool InDetDD::SiDetectorManager::setAlignableTransformLocalDelta ( ExtendedAlignableTransform * extXF, const Amg::Transform3D & localToGlobalXF, const Amg::Transform3D & delta, GeoVAlignmentStore * alignStore = nullptr )

staticinherited

Helper method to set delta transform from a local delta - Amg interface.

Definition at line 62 of file SiDetectorManager.cxx.

    {
        // ATTENTION -------------------------------------------------------- (A.S.)
        // CLHEP < -- > AMG interface method
 
        // Sets the alignable transform delta when the supplied delta is in the local
        // reconstruction frame
 
        // If the default transform to the local recostruction frame is
        // T = A*B*C*D*E
        // and the alignable transform is C with delta c and the delat in the local frame is l, then
        // A*B*C*c*D*E = A*B*C*D*E*l
        //  c = (D*E) * l * (D*E).inverse()
        //  c = (A*B*C).inverse * T * l * T.inverse() * (A*B*C)
 
        // To get default transform up and including the alignable transform,
        // we assume the next volume down is a fullphys volume and so its
        // transform is the transform we want (ie A*B*C in the above example).
 
        if (!extXF) return false;
 
        const GeoVFullPhysVol* child = extXF->child();
        if (child && extXF->alignableTransform()) {
            // the definitiv absolut transform is in CLHEP -> do the calculation in CLHEP
            const GeoTrf::Transform3D& transform = child->getDefAbsoluteTransform(alignStore);
            // calucluate the corrected delta according to the formula above
            GeoTrf::Transform3D correctedDelta = transform.inverse()*localToGlobalXF      // (A*B*C).inverse() * T
                                             * delta                                  // l
                                                 * localToGlobalXF.inverse() * transform; // T.inverse() * (A*B*C)
            extXF->alignableTransform()->setDelta(correctedDelta, alignStore);
            return true;
        } else {
            return false;
        }
    }

setCommonItems()

void InDetDD::SiDetectorManager::setCommonItems ( std::unique_ptr< const SiCommonItems > && commonItems )

inherited

Set SiCommonItems.

Definition at line 151 of file SiDetectorManager.cxx.

    {
        m_commonItems = std::move(commonItems);
    }

setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

setVersion()

void InDetDD::InDetDetectorManager::setVersion ( const Version & version )

inherited

Definition at line 43 of file InDetDetectorManager.cxx.

    {
        m_version = version;
 
    // Since default alignments are for final layout, Pixel Rome-Initial
    // layout will result in several (harmless) WARNING message. We suppress these.
    // Also the SR1 layout produce warnings due to missing parts. We suppress these also.
        m_suppressWarnings = ( (getName() == "Pixel" &&
            (version.tag() == "Pixel-01" || version.tag() == "Pixel-DC2-Initial-00"))
            || version.layout() == "SR1" || version.layout() == "SR1-EndcapC");
 
    }

tag()

const std::string & InDetDD::SiDetectorManager::tag ( ) const

inherited

Get tag used in dictionary.

Definition at line 34 of file SiDetectorManager.cxx.

    {
        return m_tag;
    }

updateAll()

void InDetDD::SiDetectorManager::updateAll ( ) const

virtualinherited

Update all caches.

Implements InDetDD::InDetDetectorManager.

Definition at line 51 of file SiDetectorManager.cxx.

    {
        for (SiDetectorElementCollection::const_iterator element_iter = getDetectorElementBegin();
        element_iter != getDetectorElementEnd();
        ++element_iter) {
            if (*element_iter) {
                (*element_iter)->updateCache();
            }
        }
    }

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_alignableTransforms

std::vector< std::unique_ptr<ExtendedAlignableTransform> > InDetDD::PixelDetectorManager::m_alignableTransforms

private

Definition at line 166 of file PixelDetectorManager.h.

m_alignfoldertype

AlignFolderType InDetDD::InDetDetectorManager::m_alignfoldertype

inherited

Definition at line 96 of file InDetDetectorManager.h.

m_commonItems

std::unique_ptr<const SiCommonItems> InDetDD::SiDetectorManager::m_commonItems

privateinherited

Definition at line 151 of file SiDetectorManager.h.

m_designs

std::vector<std::unique_ptr<const SiDetectorDesign> > InDetDD::SiDetectorManager::m_designs

privateinherited

Definition at line 150 of file SiDetectorManager.h.

m_detStore

StoreGateSvc* InDetDD::InDetDetectorManager::m_detStore

protectedinherited

Definition at line 99 of file InDetDetectorManager.h.

m_elementCollection

SiDetectorElementCollection InDetDD::PixelDetectorManager::m_elementCollection

private

Definition at line 163 of file PixelDetectorManager.h.

m_folders

std::set<std::string> InDetDD::InDetDetectorManager::m_folders

privateinherited

Definition at line 182 of file InDetDetectorManager.h.

m_globalFolders

std::set<std::string> InDetDD::InDetDetectorManager::m_globalFolders

privateinherited

Definition at line 184 of file InDetDetectorManager.h.

m_higherAlignableTransforms

std::vector< AlignableTransformMap > InDetDD::PixelDetectorManager::m_higherAlignableTransforms

private

Definition at line 165 of file PixelDetectorManager.h.

m_idHelper

const PixelID* InDetDD::PixelDetectorManager::m_idHelper

private

Definition at line 167 of file PixelDetectorManager.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_isLogical

bool InDetDD::PixelDetectorManager::m_isLogical

private

This variable switches the how the local alignment corrections are applied If true they will be calcualted on top of all of other corrections but in the default reference frame If false they will be calcualted on top of all of other corrections but in the globally aligned reference frame.

Definition at line 178 of file PixelDetectorManager.h.

m_keys

std::map<std::string, LevelInfo> InDetDD::InDetDetectorManager::m_keys

privateinherited

Definition at line 181 of file InDetDetectorManager.h.

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_numerology

SiNumerology InDetDD::SiDetectorManager::m_numerology

privateinherited

Definition at line 149 of file SiDetectorManager.h.

m_specialFolders

std::set<std::string> InDetDD::InDetDetectorManager::m_specialFolders

privateinherited

Definition at line 183 of file InDetDetectorManager.h.

m_suppressWarnings

std::atomic_bool InDetDD::InDetDetectorManager::m_suppressWarnings

mutableprotectedinherited

Definition at line 100 of file InDetDetectorManager.h.

m_tag

std::string InDetDD::SiDetectorManager::m_tag

privateinherited

Definition at line 148 of file SiDetectorManager.h.

m_version

Version InDetDD::InDetDetectorManager::m_version

privateinherited

Definition at line 180 of file InDetDetectorManager.h.

m_volume

std::vector<PVConstLink> InDetDD::PixelDetectorManager::m_volume

private

Definition at line 162 of file PixelDetectorManager.h.

s_invalidLevel

const InDetDetectorManager::LevelInfo InDetDD::InDetDetectorManager::s_invalidLevel

staticprivateinherited

Definition at line 186 of file InDetDetectorManager.h.

---

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

• PixelDetectorManager.h
• PixelDetectorManager.cxx