#include <SiTrkAlignDBTool.h>

Inheritance diagram for InDet::SiTrkAlignDBTool:

Collaboration diagram for InDet::SiTrkAlignDBTool:

Public Member Functions
	SiTrkAlignDBTool (const std::string &type, const std::string &name, const IInterface *parent)
	constructor More...

virtual	~SiTrkAlignDBTool ()
	destructor More...

virtual StatusCode	initialize ()
	initialize method More...

virtual StatusCode	finalize ()
	finalize method More...

void	writeAlignPar ()
	Writes alignment parameters. More...

void	setNtuple (TFile *file)
	sets ntuple to be used to write output More...

virtual void	setLogStream (std::ostream *os)
	sets the output stream for the logfile More...

virtual StatusCode	preSolve ()
	sets up initial alignment constants before solving More...

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) 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 More...

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. More...

Protected Attributes
std::ostream *	m_logStream = nullptr
	logfile output stream More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	updateDB ()
	updates alignment constants in the DB More...

bool	checkPixelLevel ()

bool	checkSCTLevel ()

bool	checkSiLevel ()

void	updateSiL0asL1 (Identifier idL0, const Amg::Transform3D &transform)

void	updateAsL3 (const Trk::AlignModule *module, const Amg::Transform3D &transform, double)

void	printTransform (const Amg::Transform3D &tr) const

void	updateAsL2 (const Trk::AlignModule *module, const Amg::Transform3D &transform, const int)

void	updateAsL16 (const Trk::AlignModule *module, double)

virtual Identifier	GetL1L2fromL3Identifier (const Identifier &ident, const int &level) const
	convert L3 module identifier to L1 or L2 More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
const PixelID *	m_pixHelper
	pointer to PIX detector manager More...

const SCT_ID *	m_sctHelper
	pointer to SCT detector manager More...

const AtlasDetectorID *	m_idHelper

bool	m_writeSQLFile
	flag wether to write out .db file or not (file is not needed in parallel jobs) More...

std::string	m_SQLiteTag
	tag name for the ConditionsDB More...

std::string	m_outputAlignFile
	filename for the textfile with final alignment constants More...

std::string	m_outIBLDistFile
	filename for the textfile with final IBLDist params More...

std::string	m_outGlobalFolderFile
	filename for the textfile with final GlobalFolder constants More...

bool	m_writeTextFile
	flag whether to write out the final constants into text file More...

bool	m_writeOldConstants
	flag whether to write out the initial constants into text file More...

std::string	m_oldIBLDistFile
	filename for the textfile with initial IBLDist params More...

std::string	m_oldAlignFile
	filename for the textfile with initial alignment constants More...

std::string	m_oldGlobalFolderFile
	filename for the textfile with initial GlobalFolder constants More...

bool	m_updateConstants
	flag whether to update the constants More...

ToolHandle< IInDetAlignDBTool >	m_IDAlignDBTool

ToolHandle< Trk::IAlignModuleTool >	m_alignModuleTool

ToolHandle< Trk::IGeometryManagerTool >	m_pixelGeoManager

ToolHandle< Trk::IGeometryManagerTool >	m_sctGeoManager

ToolHandle< Trk::IGeometryManagerTool >	m_siGeoManager

TFile *	m_ntuple

int	m_pixelAlignLevel

int	m_pixelAlignLevelBarrel

int	m_pixelAlignLevelEndcaps

bool	m_pixelaAlignDBM

int	m_sctAlignLevel

int	m_sctAlignLevelBarrel

int	m_sctAlignLevelEndcaps

int	m_siAlignLevel

bool	m_doSi

bool	m_doPixel

bool	m_doSCT

bool	m_writeAsL3

bool	m_writeAsL2

bool	m_writeAsL16

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 44 of file SiTrkAlignDBTool.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

◆ SiTrkAlignDBTool()

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

constructor

Definition at line 36 of file SiTrkAlignDBTool.cxx.

    : AthAlgTool(type,name,parent)
    , m_pixHelper(nullptr)
    , m_sctHelper(nullptr)
    , m_idHelper(nullptr)
    , m_writeSQLFile(true)
    , m_SQLiteTag("test_tag")
    , m_outputAlignFile("OutputSiAlignment.txt")
    , m_outIBLDistFile("OutputIBLDist.txt")
    , m_outGlobalFolderFile("OutputSiGlobalFolder.txt")
    , m_writeTextFile(true)
    , m_writeOldConstants(false)
    , m_oldIBLDistFile("OldIBLDist.txt")
    , m_oldAlignFile("OldSiAlignment.txt")
    , m_oldGlobalFolderFile("OldSiGlobalFolder.txt")
    , m_updateConstants(true)
    , m_IDAlignDBTool("InDetAlignDBTool")
    , m_alignModuleTool("Trk::AlignModuleTool/AlignModuleTool")
    , m_pixelGeoManager("")
    , m_sctGeoManager("")
    , m_siGeoManager("")
    , m_ntuple{}
    , m_pixelAlignLevel(-1)
    , m_pixelAlignLevelBarrel(-1)
    , m_pixelAlignLevelEndcaps(-1)
    , m_pixelaAlignDBM(false)
    , m_sctAlignLevel(-1)
    , m_sctAlignLevelBarrel(-1)
    , m_sctAlignLevelEndcaps(-1)
    , m_siAlignLevel(-1)
    , m_doSi(false)
    , m_doPixel(false)
    , m_doSCT(false)
    , m_writeAsL3(false)
    , m_writeAsL2(false)
    , m_writeAsL16(false)
 {
   declareInterface<ITrkAlignDBTool>(this);
  
   declareProperty("WriteSQLFile",      m_writeSQLFile);
   declareProperty("SQLiteTag",         m_SQLiteTag);
   declareProperty("OutputTextFile",    m_outputAlignFile);
   declareProperty("OutputIBLDistFile", m_outIBLDistFile);
   declareProperty("OutputGlobalFolderFile",m_outGlobalFolderFile);
   declareProperty("WriteTextFile",     m_writeTextFile);
   declareProperty("WriteOldConstants", m_writeOldConstants);
   declareProperty("OldAlignFile",      m_oldAlignFile);
   declareProperty("OldIBLDistFile",    m_oldIBLDistFile);
   declareProperty("OldGlobalFolderFile",m_oldGlobalFolderFile);
   declareProperty("UpdateConstants",   m_updateConstants);
   declareProperty("WriteAsL3",         m_writeAsL3);
  
   declareProperty("IDAlignDBTool",     m_IDAlignDBTool);
   declareProperty("AlignModuleTool",   m_alignModuleTool, "Tool for handling of align module collections");
   declareProperty("AlignmentDBM",      m_pixelaAlignDBM);
  
   declareProperty("PixelGeometryManager", m_pixelGeoManager, "Tool for handling the Pixel alignment geometry");
   declareProperty("SCTGeometryManager",   m_sctGeoManager,   "Tool for handling the SCT alignment geometry");
   declareProperty("SiGeometryManager",    m_siGeoManager,    "Tool for handling the Silicon alignment geometry");
 }

◆ ~SiTrkAlignDBTool()

InDet::SiTrkAlignDBTool::~SiTrkAlignDBTool ( )

virtual

destructor

Definition at line 98 of file SiTrkAlignDBTool.cxx.

99 {}

Member Function Documentation

◆ checkPixelLevel()

bool InDet::SiTrkAlignDBTool::checkPixelLevel ( )

private

Definition at line 233 of file SiTrkAlignDBTool.cxx.

 {
   if(!m_pixelaAlignDBM && m_pixelAlignLevel==1)
     return true;
  
   if(!m_pixelaAlignDBM && m_pixelAlignLevel==11){// the new one with IBL (seems simplest solution)
     m_writeAsL2 = true;
     return true;
   }
  
   if(m_pixelaAlignDBM && (m_pixelAlignLevel==1 || m_pixelAlignLevel==2 )){// the new one with DBM (seems simplest solution)
     m_writeAsL2 = true;
     return true;
   }
   
   if(m_pixelaAlignDBM && m_pixelAlignLevel==3 ){// the new one with DBM (seems simplest solution)
     return true;
   }
   
   if(!m_pixelaAlignDBM && m_pixelAlignLevel==16){// the new IBL bowing level                   
     m_writeAsL16 = true; 
     // We need some fix here with respect to current CL implementation! There bowing is updated at L3!! (this would currently do both..)
     return true;
   }
  
   if (m_writeAsL3)
     return true;
  
   m_pixelAlignLevel = -1;
  
   bool ok = false;
  
   switch(m_pixelAlignLevelBarrel) {
     case 11: case 12: case 15: case 16: case 2: case 22: case 26: case 27: case 3:
       ok = true;
       break;
     default:
       ATH_MSG_ERROR(" Storing of alignment constants for level "<<m_pixelAlignLevelBarrel<<" for Pixel Barrel is not implemented.");
       m_pixelAlignLevelBarrel = -1;
       ok = false;
       break;
   }
  
   switch(m_pixelAlignLevelEndcaps) {
     case 11: case 16: case 2: case 3: case 12:
       ok = ok && true;
       break;
     default:
       ATH_MSG_ERROR(" Storing of alignment constants for level "<<m_pixelAlignLevelEndcaps<<" for Pixel Endcaps is not implemented.");
       m_pixelAlignLevelEndcaps = -1;
       ok = false;
       break;
   }
  
   return ok;
 }

◆ checkSCTLevel()

bool InDet::SiTrkAlignDBTool::checkSCTLevel ( )

private

Definition at line 291 of file SiTrkAlignDBTool.cxx.

 {
   if(m_sctAlignLevel==0)
     return true;
  
   if (m_writeAsL3)
     return true;
  
   m_sctAlignLevel = -1;
  
   bool ok = false;
  
   switch(m_sctAlignLevelBarrel) {
     case 1: case 2: case 27: case 3:
       ok = true;
       break;
     default:
       ATH_MSG_ERROR(" Storing of alignment constants for level "<<m_sctAlignLevelBarrel<<" for SCT Barrel is not implemented.");
       m_sctAlignLevelBarrel = -1;
       ok = false;
       break;
   }
  
   switch(m_sctAlignLevelEndcaps) {
     case 1: case 2: case 25: case 3:
       ok = ok && true;
       break;
     default:
       ATH_MSG_ERROR(" Storing of alignment constants for level "<<m_sctAlignLevelEndcaps<<" for SCT Endcaps is not implemented.");
       m_sctAlignLevelEndcaps = -1;
       ok = false;
       break;
   }
  
   return ok;
 }

◆ checkSiLevel()

bool InDet::SiTrkAlignDBTool::checkSiLevel ( )

private

Definition at line 329 of file SiTrkAlignDBTool.cxx.

 {
   if (m_writeAsL3)
     return true;
  
   switch(m_siAlignLevel) {
     case 0:
     case 1:
     case 2:
     case 3:
       return true;
  
     default:
       return false;
   }
 }

◆ declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	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());
  
   }

◆ declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & 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.

95 { return m_detStore; }

◆ evtStore() [1/2]

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.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { 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::SiTrkAlignDBTool::finalize ( )

virtual

finalize method

Implements Trk::ITrkAlignDBTool.

Definition at line 225 of file SiTrkAlignDBTool.cxx.

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

◆ GetL1L2fromL3Identifier()

Identifier InDet::SiTrkAlignDBTool::GetL1L2fromL3Identifier	(	const Identifier &	ident,
		const int &	level
	)		const

privatevirtual

convert L3 module identifier to L1 or L2

< no translation needed

check whether PIX

< Whole pixel det. at L1

check whether SCT

< barrel + 2 x EC at L1

< take care of the case where level != 1,2,3 or ident neither pix nor sct

Definition at line 866 of file SiTrkAlignDBTool.cxx.

                                                                                                      {
   if( level == 3 ) return ident ; 
  
   if( m_pixHelper->is_pixel(ident) ) {
     if( level == 1 ) {
       return m_pixHelper->wafer_id( 0, 0, 0, 0 ) ; 
     }
     if( level == 2 ) {
       int barrel_ec  = m_pixHelper->barrel_ec(  ident ) ;
       int layer_disk = m_pixHelper->layer_disk( ident ) ;
       int phi_sector = m_pixHelper->phi_module( ident ) ;
       if (abs(barrel_ec)==4) return m_pixHelper->wafer_id( barrel_ec, layer_disk, phi_sector, 0 ) ;
       else return m_pixHelper->wafer_id( barrel_ec, layer_disk, 0, 0 ) ;
     }
   }
   if( m_sctHelper->is_sct(ident) ) {
     if( level == 1 ) {
       int barrel_ec  = m_sctHelper->barrel_ec(  ident ) ;
       return m_sctHelper->wafer_id( barrel_ec, 0, 0, 0, 0 ) ; 
     }
     if( level == 2 ) {
       int barrel_ec  = m_sctHelper->barrel_ec(  ident ) ;
       int layer_disk = m_sctHelper->layer_disk( ident ) ;
       return m_sctHelper->wafer_id( barrel_ec, layer_disk, 0, 0, 0 ) ;
     }
   }
   return ident ; 
 }

◆ initialize()

StatusCode InDet::SiTrkAlignDBTool::initialize ( )

virtual

initialize method

Implements Trk::ITrkAlignDBTool.

Definition at line 102 of file SiTrkAlignDBTool.cxx.

                                         {
  
   ATH_MSG_DEBUG("initialize() of SiTrkAlignDBTool");
  
   if(!m_updateConstants) {
     ATH_MSG_WARNING("Not updating Silicon constants!");
     return StatusCode::SUCCESS;
   }
  
   // Si alignment database tool
   ATH_CHECK( m_IDAlignDBTool.retrieve() );
   // get AlignModuleTool
   ATH_CHECK( m_alignModuleTool.retrieve() );
   // retrieve Pixel helper
   ATH_CHECK( detStore()->retrieve(m_pixHelper) );
   // retrieve SCT helper
   ATH_CHECK( detStore()->retrieve(m_sctHelper) );
   ATH_CHECK(detStore()->retrieve(m_idHelper, "AtlasID"));
  
   // check which geometry managers are loaded and which are to be used
  
   // get SiGeometryManagerTool
   if( !m_siGeoManager.empty() ) {
     ATH_CHECK( m_siGeoManager.retrieve());
  
     ATH_MSG_DEBUG("retrieved " << m_siGeoManager);
     m_siAlignLevel = m_siGeoManager->alignLevel();
     ATH_MSG_INFO(" Requested update of Level "<<m_siAlignLevel<<" alignment constants for Silicon");
  
     // check alignment level but don't stop on error
     // if unknown level given, set to -1
     if(!checkSiLevel()) {
       ATH_MSG_WARNING(" Storing of alignment constants for level "<<m_siAlignLevel<<" for Silicon is not implemented.");
       m_siAlignLevel=-1;
     }
   }
  
   // for anything else then level 0 we need the other two geometry managers
   if(m_siAlignLevel!=0) {
  
     // get PixelGeometryManagerTool
     if( !m_pixelGeoManager.empty() ) {
       ATH_CHECK( m_pixelGeoManager.retrieve());
       ATH_MSG_DEBUG("retrieved " << m_pixelGeoManager);
       m_pixelAlignLevel = m_pixelGeoManager->alignLevel();
       m_pixelAlignLevelBarrel = m_pixelGeoManager->alignLevelBarrel();
       m_pixelAlignLevelEndcaps = m_pixelGeoManager->alignLevelEndcaps();
 //      m_pixelaAlignDBM = m_pixelGeoManager->alignDBM();
       ATH_MSG_DEBUG("Retrieved Pixel Global alignment level: "<<m_pixelAlignLevel);
       ATH_MSG_DEBUG("Retrieved Pixel Barrel alignment level: "<<m_pixelAlignLevelBarrel);
       ATH_MSG_DEBUG("Retrieved Pixel Endcap alignment level: "<<m_pixelAlignLevelEndcaps);
       ATH_MSG_DEBUG("Retrieved DBM alignment configuration: "<<m_pixelaAlignDBM);
       
       if(m_pixelAlignLevel == 1)
         ATH_MSG_INFO(" Requested update of Level "<<m_pixelAlignLevel<<" alignment constants for Pixel");
       else {
         ATH_MSG_INFO(" Requested update of Level "<<m_pixelAlignLevelBarrel<<" alignment constants for Pixel Barrel");
         ATH_MSG_INFO(" Requested update of Level "<<m_pixelAlignLevelEndcaps<<" alignment constants for Pixel Endcaps");
       }
  
        // check alignment level but don't stop on error
       if(!checkPixelLevel()) {
         ATH_MSG_FATAL(" Pixel alignment level not given.");
 //        return StatusCode::FAILURE;
       }
     }
  
     // get SCTGeometryManagerTool
     if( !m_sctGeoManager.empty() ) {
       ATH_CHECK( m_sctGeoManager.retrieve());
       ATH_MSG_DEBUG("retrieved " << m_sctGeoManager);
       m_sctAlignLevel = m_sctGeoManager->alignLevel();
       m_sctAlignLevelBarrel = m_sctGeoManager->alignLevelBarrel();
       m_sctAlignLevelEndcaps = m_sctGeoManager->alignLevelEndcaps();
       ATH_MSG_DEBUG("Retrieved SCT Global alignment level: "<<m_sctAlignLevel);
       ATH_MSG_DEBUG("Retrieved SCT Barrel alignment level: "<<m_sctAlignLevelBarrel);
       ATH_MSG_DEBUG("Retrieved SCT Endcap alignment level: "<<m_sctAlignLevelEndcaps);
  
         if(m_sctAlignLevel == 0)
         ATH_MSG_INFO(" Requested update of Level "<<m_sctAlignLevel<<" alignment constants for SCT");
       else {
         ATH_MSG_INFO(" Requested update of Level "<<m_sctAlignLevelBarrel<<" alignment constants for SCT Barrel");
         ATH_MSG_INFO(" Requested update of Level "<<m_sctAlignLevelEndcaps<<" alignment constants for SCT Endcaps");
       }
  
        // check alignment level but don't stop on error
       if(!checkSCTLevel()) {
         ATH_MSG_FATAL(" SCT alignment level not given.");
 //        return StatusCode::FAILURE;
       }
     }
  
     // if we're here it means that we need at least one the two managers
     if(m_pixelAlignLevel<0 && m_pixelAlignLevelBarrel<0 && m_pixelAlignLevelEndcaps<0
        && m_sctAlignLevel<0 && m_sctAlignLevelBarrel<0 && m_sctAlignLevelEndcaps<0) {
       ATH_MSG_FATAL( " No geometry manager available or alignment level not given. " );
       // we don't want to stop if there's errors
       // so we comment out the return statement
 //      return StatusCode::FAILURE;
     }
   }
  
   m_doSi = (m_siAlignLevel==0);
   m_doPixel = (checkPixelLevel() && !m_doSi);
   m_doSCT = (checkSCTLevel() && !m_doSi);
  
   if (m_writeAsL3)
     ATH_MSG_INFO(" Storing as level 3 constants.");
  
   if (m_writeAsL16)
     ATH_MSG_INFO(" Storing IBLDist as level 16 constants in separate DB");
  
   if(!m_doPixel && !m_doSCT && !m_doSi) {
     ATH_MSG_FATAL(" No geometry manager available or alignment level not given.");
     ATH_MSG_FATAL(" DB not updated.");
 //    return;
   }
   ATH_MSG_DEBUG("Made it till here..... "<<m_doSi<<" "<<m_doPixel<<" "<<m_doSCT);
   ATH_MSG_DEBUG("update constants..... "<<m_updateConstants);
   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 & ITrkAlignDBTool::interfaceID ( )

inlinestaticinherited

Definition at line 58 of file ITrkAlignDBTool.h.

59 { return IID_TRKALIGNINTERFACES_ITrkAlignDBTool; }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ 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.

◆ preSolve()

virtual StatusCode Trk::ITrkAlignDBTool::preSolve ( )

inlinevirtualinherited

sets up initial alignment constants before solving

Definition at line 52 of file ITrkAlignDBTool.h.

52 { return StatusCode::SUCCESS; }

◆ printTransform()

void InDet::SiTrkAlignDBTool::printTransform ( const Amg::Transform3D & tr ) const

private

Definition at line 854 of file SiTrkAlignDBTool.cxx.

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

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

◆ setLogStream()

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

inlinevirtualinherited

sets the output stream for the logfile

Definition at line 49 of file ITrkAlignDBTool.h.

49 { m_logStream = os; }

◆ setNtuple()

void InDet::SiTrkAlignDBTool::setNtuple ( TFile * file )

inlinevirtual

sets ntuple to be used to write output

Reimplemented from Trk::ITrkAlignDBTool.

Definition at line 64 of file SiTrkAlignDBTool.h.

65 { m_ntuple=file; }

◆ 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 DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ 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.

◆ updateAsL16()

void InDet::SiTrkAlignDBTool::updateAsL16	(	const Trk::AlignModule *	module,
		double	bowx
	)

private

Definition at line 807 of file SiTrkAlignDBTool.cxx.

 {
   ATH_MSG_INFO("Updating constants for module "<<module->name()<<" as level 16 IBLDist constants (new)");
  
   Trk::AlignModule::DetectorType det = Trk::AlignModule::Pixel;
   const std::vector<const Trk::TrkDetElementBase *> * elements = module->detElementCollection(det);
  
   if( !elements ) {
     ATH_MSG_FATAL("no elements of type "<<det);
     return;
   }
  
   ATH_MSG_DEBUG("looping over "<<elements->size()<<" elements");
   std::vector<int> stave_phis; // where we store the unique identifiers we want to update                                                                                
   std::vector<const Trk::TrkDetElementBase*>::const_iterator ielem     = elements->begin();
   std::vector<const Trk::TrkDetElementBase*>::const_iterator ielem_end = elements->end();
   for ( ; ielem != ielem_end; ++ielem) {
     const InDetDD::SiDetectorElement * sielem = dynamic_cast<const InDetDD::SiDetectorElement *>(*ielem);
     if (!sielem) {
       ATH_MSG_WARNING("Should be Silicon detector element but is not. Skipping.");
       continue;
     }
  
     Identifier elemID = sielem->identify();
     const int this_stave_phi = m_pixHelper->phi_module( elemID );
  
     ATH_MSG_DEBUG("module ID: "<<m_idHelper->show_to_string(elemID,nullptr,'/'));
     ATH_MSG_DEBUG("phi identifier: "<<m_pixHelper->phi_module( elemID )) ;    
  
     std::vector<int>::const_iterator ix = find(stave_phis.begin(),stave_phis.end(),this_stave_phi); // check whether it is unique
     if (ix==stave_phis.end()) {
       stave_phis.push_back(this_stave_phi);
  
       // tweak applies the transform onto already existing transform in the DB                                              
       if( !(m_IDAlignDBTool->tweakIBLDist(this_stave_phi,bowx)) ) {
           ATH_MSG_ERROR("Error tweaking IBLDist DB for stave "<<this_stave_phi);
       }
     }
     else {
       ATH_MSG_DEBUG("Skipping ModuleID: "<<m_idHelper->show_to_string(elemID,nullptr,'/')<<" --> not unique");
     }
     
   } // end loop over detElements                                                                                                                                                 
 }

◆ updateAsL2()

void InDet::SiTrkAlignDBTool::updateAsL2	(	const Trk::AlignModule *	module,
		const Amg::Transform3D &	transform,
		const int	level
	)

private

Definition at line 736 of file SiTrkAlignDBTool.cxx.

 {
   ATH_MSG_INFO("Updating constants for module "<<module->name()<<" as level 2 constants (new)");
     
   // check the level we store in
   ATH_MSG_INFO("DB aligning at level: "<<level);
   
   Trk::AlignModule::DetectorType det = Trk::AlignModule::Pixel;
   const std::vector<const Trk::TrkDetElementBase *> * elements = module->detElementCollection(det);
  
   if( !elements ) {
     ATH_MSG_FATAL("no elements of type "<<det);
     return;
   }
  
   ATH_MSG_DEBUG("looping over "<<elements->size()<<" elements");
   std::vector<std::string> level_mods; // where we store the unique identifiers we want to update
     
   std::vector<const Trk::TrkDetElementBase*>::const_iterator ielem     = elements->begin();
   std::vector<const Trk::TrkDetElementBase*>::const_iterator ielem_end = elements->end();
   for ( ; ielem != ielem_end; ++ielem) {
     const InDetDD::SiDetectorElement * sielem = dynamic_cast<const InDetDD::SiDetectorElement *>(*ielem);
     if (!sielem) {
       ATH_MSG_WARNING("Should be Silicon detector element but is not. Skipping.");
       continue;
     }
     
     Identifier elemID = sielem->identify();
     Identifier level_mod = GetL1L2fromL3Identifier(elemID,level);
     std::string key = m_IDAlignDBTool->dirkey(level_mod,level+1);
     
     ATH_MSG_DEBUG("TESTING value: "<<level_mod);
     ATH_MSG_DEBUG("TESTING ID: "<<m_idHelper->show_to_string(level_mod,nullptr,'/'));
     ATH_MSG_DEBUG("TESTING retrieved key: "<<key);
     std::vector<std::string>::const_iterator ix = find(level_mods.begin(),level_mods.end(),key); // check whether it is unique
     if (ix==level_mods.end()) {
       level_mods.push_back(key);
       
       ATH_MSG_DEBUG("-----------------------------------------------------");
       
       ATH_MSG_DEBUG("ModuleID value: "<<level_mod);
       ATH_MSG_DEBUG("ModuleID: "<<m_idHelper->show_to_string(level_mod,nullptr,'/'));
       
       // for levels 1 and 2 the DB frame equals to the global frame but the align
       // frame doesnot, so we have to apply an additional transform
       const Amg::Transform3D& dbFrameToAlignFrame = module->globalFrameToAlignFrame();
         
       ATH_MSG_DEBUG("DB to align");
       printTransform(dbFrameToAlignFrame);
         
       Amg::Transform3D dbtransform = dbFrameToAlignFrame.inverse() * transform * dbFrameToAlignFrame;
         
       ATH_MSG_DEBUG("DB transform");
       printTransform(dbtransform);
         
         
       // tweak applies the transform onto already existing transform in the DB
       if( !(m_IDAlignDBTool->tweakTrans(level_mod, level, dbtransform)) ) {
         // set replaces the transform in the DB
         if( !(m_IDAlignDBTool->setTrans(level_mod, level, dbtransform)) )
           ATH_MSG_ERROR("Error setting level "<<level<<" constants for element "<<level_mod<<" in module "<<module->name());
       }
     }
     else {
       ATH_MSG_DEBUG("Skipping ModuleID: "<<m_idHelper->show_to_string(elemID,nullptr,'/')<<" --> not unique");
     }
   } // end loop over detElements
 }

◆ updateAsL3()

void InDet::SiTrkAlignDBTool::updateAsL3	(	const Trk::AlignModule *	module,
		const Amg::Transform3D &	transform,
		double
	)

private

  if (apBowX!=0){
double z = sielem->center()[2];        
const double  y0y0  = 366.5*366.5;

ATH_MSG_DEBUG("Z position : "<< z); double bowx = apBowX * ( z*z - y0y0 ) / y0y0; This is in the module frame, as bowing corrections are directly L3 Amg::Translation3D translation_bow(bowx , 0, 0 ); dbtransform *= translation_bow; }

Definition at line 639 of file SiTrkAlignDBTool.cxx.

 {
   ATH_MSG_INFO("Updating constants for module "<<module->name()<<" as level 3 constants");
  
   // store as level 3 constants
   const int level = 3;
  
   // transform from global frame to align frame
   const Amg::Transform3D& globaltoalign = module->globalFrameToAlignFrame();
  
   // same code pplies to both Pixel and SCT
   for (unsigned int idet=0; idet<2; idet++) {
  
     Trk::AlignModule::DetectorType det = idet ? Trk::AlignModule::SCT : Trk::AlignModule::Pixel;
  
     ATH_MSG_DEBUG("getting elements of type "<<det<<" from the module");
  
     const std::vector<const Trk::TrkDetElementBase *> * elements = module->detElementCollection(det);
  
     if( !elements ) {
       ATH_MSG_DEBUG("no elements of type "<<det);
       continue;
     }
  
     ATH_MSG_DEBUG("looping over "<<elements->size()<<" elements");
  
     std::vector<const Trk::TrkDetElementBase*>::const_iterator ielem     = elements->begin();
     std::vector<const Trk::TrkDetElementBase*>::const_iterator ielem_end = elements->end();
     for ( ; ielem != ielem_end; ++ielem) {
       const InDetDD::SiDetectorElement * sielem = dynamic_cast<const InDetDD::SiDetectorElement *>(*ielem);
       if (!sielem) {
         ATH_MSG_WARNING("Should be Silicon detector element but is not. Skipping.");
         continue;
       }  
  
       Identifier elemID = sielem->identify();
  
       // for the SCT the 0 side Identifier is used in the DB
       if (m_idHelper->is_sct(elemID)) {
         if (m_sctHelper->side(elemID) > 0)
           continue;
       }
  
       ATH_MSG_DEBUG("-----------------------------------------------------");
  
       ATH_MSG_DEBUG("ModuleID value: "<<elemID);
       ATH_MSG_DEBUG("ModuleID: "<<m_idHelper->show_to_string(elemID,nullptr,'/'));
  
       // this should work for both Pixel and SCT
       // module transform is the transform between the db frame and global frame
       Amg::Transform3D dbtoglobal = sielem->moduleTransform();
       //      ATH_MSG_DEBUG("DB to Global");
       //      printTransform(dbtoglobal);
  
       //      ATH_MSG_DEBUG("Global to Align");
       //      printTransform(globaltoalign);
  
       // from DB frame to alignment frame transform
       Amg::Transform3D dbtoalign = globaltoalign * dbtoglobal;
       //      ATH_MSG_DEBUG("DB to Align");
       //      printTransform(dbtoalign);
  
       ATH_MSG_DEBUG("Alignment transform");
       printTransform(transform);
  
       // alignment transform in DB frame
       Amg::Transform3D dbtransform = dbtoalign.inverse() * transform * dbtoalign; 
  
       // This should work as Bowing is in DB frame, i.e. local module frame
       // This implementation is need in the CL as of Oct2015; 
       // For now, just comment for offline
       ATH_MSG_DEBUG("DB transform");
       printTransform(dbtransform);
  
       // tweak applies the transform onto already existing transform in the DB
       if( !(m_IDAlignDBTool->tweakTrans(elemID, level, dbtransform)) ) {
         // set replaces the transform in the DB
         if( !(m_IDAlignDBTool->setTrans(elemID, level, dbtransform)) )
           ATH_MSG_ERROR("Error setting level "<<level<<" constants for element "<<elemID<<" in module "<<module->name());
       }
     } // end loop over detElements
     ATH_MSG_DEBUG("-----------------------------------------------------");
   } // end loop over Pixel and SCT
 }

◆ updateDB()

void InDet::SiTrkAlignDBTool::updateDB ( )

private

updates alignment constants in the DB

Definition at line 380 of file SiTrkAlignDBTool.cxx.

 {
   ATH_MSG_DEBUG(" in SiTrkAlignDBTool::updateDB()");
  
   Amg::Transform3D dbtransform;
  
   ATH_MSG_INFO("-------------------------------------------------");
   // updating level 1, level 2 or level 3 constants in the DB
   if(m_doSi)
     ATH_MSG_INFO("Updating Level "<<m_siAlignLevel<<" alignment constants for Silicon");
  
   if(m_doPixel) {
     if(m_pixelAlignLevel==1)
       ATH_MSG_INFO("Updating Level "<<m_pixelAlignLevel<<" alignment constants for Pixel");
     else {
       ATH_MSG_INFO("Updating Level "<<m_pixelAlignLevelBarrel<<" alignment constants for Pixel Barrel");
       ATH_MSG_INFO("Updating Level "<<m_pixelAlignLevelEndcaps<<" alignment constants for Pixel Endcaps");
     }
   }
  
   if(m_doSCT) {
     if(m_sctAlignLevel==0)
       ATH_MSG_INFO("Updating Level "<<m_sctAlignLevel<<" alignment constants for SCT");
     else {
       ATH_MSG_INFO("Updating Level "<<m_sctAlignLevelBarrel<<" alignment constants for SCT Barrel");
       ATH_MSG_INFO("Updating Level "<<m_sctAlignLevelEndcaps<<" alignment constants for SCT Endcaps");
     }
   }
  
   if (m_writeAsL3)
     ATH_MSG_INFO("Updating Level 3 constants.");
  
   // loop over all AlignModules
   std::vector<Trk::AlignModule *>::const_iterator imod     = m_alignModuleTool->getAlignModuleCollection()->begin();
   std::vector<Trk::AlignModule *>::const_iterator imod_end = m_alignModuleTool->getAlignModuleCollection()->end();
   for( ; imod != imod_end; ++imod ) {
  
     Trk::AlignModule * module = *imod;
  
     ATH_MSG_DEBUG("-------------------------------------------------------");
     ATH_MSG_DEBUG("Alignment corrections for module: "<<module->name());
  
     // the the identifier of the module for the DB
     Identifier modID = module->identify();
  
     // The Si AlignModules of levels 0, 1, 2 and 3 use the Silicon Identifier.
     if(!m_idHelper->is_pixel(modID) && !m_idHelper->is_sct(modID)) {
       ATH_MSG_DEBUG("Module with ID: "<<modID<<" is not a Silicon module. Skipping.");
       continue;
     }
  
     // get active alignment parameters for the module
     DataVector<Trk::AlignPar> * alignPars =  m_alignModuleTool->getAlignPars(module);
     if(alignPars->size() == 0) {
       ATH_MSG_INFO("Alignment constants for module "<<module->name()<<" not available.");
       continue;
     }
  
     ATH_MSG_DEBUG("Retrieved alignment constants for module "<<module->name());
     ATH_MSG_DEBUG("Module ID : "<<modID);
     ATH_MSG_DEBUG("Expanded module ID : "<<m_idHelper->show_to_string(modID,nullptr,'/'));
  
     // get final alignment parameters
     DataVector<Trk::AlignPar> * fullAlignPars = m_alignModuleTool->getFullAlignPars(module);
     double apTraX = fullAlignPars->at(Trk::AlignModule::TransX)->par();
     double apTraY = fullAlignPars->at(Trk::AlignModule::TransY)->par();
     double apTraZ = fullAlignPars->at(Trk::AlignModule::TransZ)->par();
     double apRotX = fullAlignPars->at(Trk::AlignModule::RotX)->par();
     double apRotY = fullAlignPars->at(Trk::AlignModule::RotY)->par();
     double apRotZ = fullAlignPars->at(Trk::AlignModule::RotZ)->par();
     
     // Need to make sure not to call SCT modules - only defined for PIX 
     double apBowX = 0;
     // This should be sufficient as bowing-DoF is only enabled for IBL
     if(m_idHelper->is_pixel(modID) && m_pixHelper->is_barrel(modID) ){
       apBowX = fullAlignPars->at(Trk::AlignModule::BowX)->par();
     
       // The last check ensures that it also updates for L11 in case bowing DoF is enables.
       if(m_pixelAlignLevelBarrel==16 || m_writeAsL16 || apBowX!=0){     
     updateAsL16(module,apBowX);  
     // Can we get a baseline as well? include here in case it is possible!
     // Allow other DoF to be updated as well!
     // Initial CL (Oct. 2016) update at L3 --> commented for now
       }
     }
     
     Amg::Translation3D translation( apTraX, apTraY, apTraZ );
     Amg::Transform3D transform = translation * Amg::RotationMatrix3D::Identity();
     transform *= Amg::AngleAxis3D(apRotZ, Amg::Vector3D(0.,0.,1.));
     transform *= Amg::AngleAxis3D(apRotY, Amg::Vector3D(0.,1.,0.));
     transform *= Amg::AngleAxis3D(apRotX, Amg::Vector3D(1.,0.,0.));
     
     
     ATH_MSG_DEBUG(" - translation X = "<<apTraX);
     ATH_MSG_DEBUG(" - translation Y = "<<apTraY);
     ATH_MSG_DEBUG(" - translation Z = "<<apTraZ);
     ATH_MSG_DEBUG(" - rotation X    = "<<apRotX);
     ATH_MSG_DEBUG(" - rotation Y    = "<<apRotY);
     ATH_MSG_DEBUG(" - rotation Z    = "<<apRotZ);
  
     int level=-1;
  
     // this is really ugly and quite confusing !!!
     // but will do for now...
     if(m_idHelper->is_pixel(modID)) {
       if(!m_doPixel && !m_doSi) {
         ATH_MSG_INFO("--> Pixel module - NOT storing");
         continue;
       }
       if(m_doPixel) {
         if(m_pixelAlignLevel==1 || m_pixelAlignLevel==2 || m_pixelAlignLevel==3)
           level = m_pixelAlignLevel;
         else if(m_pixHelper->is_barrel(modID))
           level = m_pixelAlignLevelBarrel;
         else
           level = m_pixelAlignLevelEndcaps;
       }
       else // this makes sure that Si L0 is handled since pixel identifier is used for it
         level = m_siAlignLevel;
     }
     else if(m_idHelper->is_sct(modID)) {
       if(!m_doSCT) {
         ATH_MSG_INFO("--> SCT module - NOT storing");
         continue;
       }
       if(m_sctAlignLevel==0)
         level = m_sctAlignLevel;
       else if(m_sctHelper->is_barrel(modID))
         level = m_sctAlignLevelBarrel;
       else
         level = m_sctAlignLevelEndcaps;
     }
     // to better understand what is going on
     ATH_MSG_DEBUG("We are having level: "<<level);
  
     if (m_writeAsL3) {
       updateAsL3(module,transform,apBowX);
       continue;
     }
  
     if (m_idHelper->is_pixel(modID) && m_writeAsL2) {
       level=2;
       updateAsL2(module,transform,level);
       continue;
     }
     
     // calculate the DB transform
     switch(level) {
  
       case -1:
         // level not recognized, do nothing
         continue;
  
       case 0:
         updateSiL0asL1(modID,transform);
         continue;
  
       case 1: case 2: {
         // for levels 1 and 2 the DB frame equals to the global frame but the align
         // frame doesnot, so we have to apply an additional transform
         Amg::Transform3D dbFrameToAlignFrame = module->globalFrameToAlignFrame();
  
         ATH_MSG_DEBUG("DB to align");
         printTransform(dbFrameToAlignFrame);
  
         dbtransform = dbFrameToAlignFrame.inverse() * transform * dbFrameToAlignFrame;
  
         ATH_MSG_DEBUG("DB transform");
         printTransform(dbtransform);
         break;
       }
  
       case 3:
         // for level 3 alignment the DB frame and the align frame equals to the local frame
         // so we don't apply any additional transform
         dbtransform = transform;
  
         ATH_MSG_DEBUG("DB transform");
         printTransform(dbtransform);
         break;
  
       default:
         // otherwise we split the AlignModule into level 3 modules
         ATH_MSG_DEBUG("We are doing this one");
         updateAsL3(module,transform,apBowX);
         continue;
     }
  
     // tweak applies the transform onto already existing transform in the DB
     if( m_IDAlignDBTool->tweakTrans(modID, level, dbtransform) )
       ATH_MSG_INFO("Module \'"<<module->name()<<"\': Level "<<level<<" constants updated.");
     else {
       // set replaces the transform in the DB
       ATH_MSG_DEBUG("Module \'"<<module->name()<<"\': Call to tweak failed. Constants probably don't exist yet. Use set instead.");
       if( m_IDAlignDBTool->setTrans(modID, level, dbtransform) )
         ATH_MSG_INFO("Module \'"<<module->name()<<"\': Level "<<level<<" constants updated.");
       else
         ATH_MSG_ERROR("Error setting level "<<level<<" constants for module "<<module->name());
     }
  
   }
   ATH_MSG_INFO("-------------------------------------------------");
  
 }

◆ updateSiL0asL1()

void InDet::SiTrkAlignDBTool::updateSiL0asL1	(	Identifier	idL0,
		const Amg::Transform3D &	transform
	)

private

Definition at line 586 of file SiTrkAlignDBTool.cxx.

 {
   // splitting the whole Silicon into Pixel, SCT barrel and 2 SCT endcaps
   // the id has to be equal to L1 Pixel ID for Silicon level 0 update
   // or to SCT barrel L1 ID for SCT level 0 update
   if( idL0 != m_pixHelper->wafer_id(0,0,0,0) && idL0 != m_sctHelper->wafer_id(0,0,0,0,0) ) {
     ATH_MSG_ERROR("wrong identifier for L0 Silicon in updateSiL0asL1()");
     return;
   }
  
   // for both L0 and L1 the global frame is equal to the DB frame
   // for L0 it is also equal to align frame so no additional transform is needed
  
   if(idL0 == m_pixHelper->wafer_id(0,0,0,0)) {
     ATH_MSG_INFO("Splitting Silicon into L1 modules and updating in the DB");
  
     // Pixel
     Identifier idPixel = m_pixHelper->wafer_id(0,0,0,0);
     ATH_MSG_INFO("Updating L1 constants for Pixel");
     // tweak applies the transform onto already existing transform in the DB
     if( m_IDAlignDBTool->tweakTrans(idPixel, 1, transform) )
       ATH_MSG_INFO("Call to tweak OK. Constants updated.");
     else {
       // set replaces the transform in the DB
       ATH_MSG_INFO("Call to tweak failed. Constants probably don't exist yet. Use set instead.");
       if( !m_IDAlignDBTool->setTrans(idPixel, 1, transform) )
         ATH_MSG_ERROR("Error setting constants.");
     }
   }
   else
     ATH_MSG_INFO("Splitting SCT into L1 modules and updating in the DB");
  
   // SCT
   std::vector<Identifier> ids;
   ids.push_back(m_sctHelper->wafer_id(0,0,0,0,0));  // SCT barrel
   ids.push_back(m_sctHelper->wafer_id(-2,0,0,0,0)); // SCT endcap A
   ids.push_back(m_sctHelper->wafer_id(2,0,0,0,0));  // SCT endcap C
  
   ATH_MSG_INFO("Updating L1 constants for SCT");
   for(unsigned int i=0; i<ids.size(); ++i) {
     // tweak applies the transform onto already existing transform in the DB
     if( m_IDAlignDBTool->tweakTrans(ids[i], 1, transform) )
       ATH_MSG_INFO("Module "<<i<<". - id "<<ids[i]<<": Call to tweak OK. Constants updated.");
     else {
       // set replaces the transform in the DB
       ATH_MSG_INFO("Module "<<i<<". - id "<<ids[i]<<": Call to tweak failed. Constants probably don't exist yet. Use set instead.");
       if( !m_IDAlignDBTool->setTrans(ids[i], 1, transform) )
         ATH_MSG_ERROR("Error setting constants.");
     }
   }
 }

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

◆ writeAlignPar()

void InDet::SiTrkAlignDBTool::writeAlignPar ( )

virtual

Writes alignment parameters.

Reimplemented from Trk::ITrkAlignDBTool.

Definition at line 347 of file SiTrkAlignDBTool.cxx.

 {
   ATH_MSG_DEBUG("in SiTrkAlignDBTool::writeAlignPar()");
  
   // do we want to update the DB?
   if(!m_updateConstants) {
     ATH_MSG_INFO("Update of Silicon constants not requested!");
     return;
   }
  
   // write initial alignment constants
   if(m_writeOldConstants) {
     ATH_MSG_INFO("Writing old Silicon alignment constants to file "<<m_oldAlignFile);
     m_IDAlignDBTool->writeFile(false, m_oldAlignFile);
     m_IDAlignDBTool->writeIBLDistFile(m_oldIBLDistFile);
     m_IDAlignDBTool->writeGlobalFolderFile(m_oldGlobalFolderFile);
   }
  
   // now update the DB
   updateDB();
  
   // and write the final alignment constants
   if(m_writeTextFile) {
     ATH_MSG_INFO("Writing Silicon alignment constants to file "<<m_outputAlignFile);
     m_IDAlignDBTool->writeFile(false, m_outputAlignFile);
     m_IDAlignDBTool->writeIBLDistFile(m_outIBLDistFile);
     m_IDAlignDBTool->writeGlobalFolderFile(m_outGlobalFolderFile);
   }
  
   return;
 }

Member Data Documentation

◆ m_alignModuleTool

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

private

Definition at line 106 of file SiTrkAlignDBTool.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_doPixel

bool InDet::SiTrkAlignDBTool::m_doPixel

private

Definition at line 126 of file SiTrkAlignDBTool.h.

◆ m_doSCT

bool InDet::SiTrkAlignDBTool::m_doSCT

private

Definition at line 127 of file SiTrkAlignDBTool.h.

◆ m_doSi

bool InDet::SiTrkAlignDBTool::m_doSi

private

Definition at line 125 of file SiTrkAlignDBTool.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_IDAlignDBTool

ToolHandle<IInDetAlignDBTool> InDet::SiTrkAlignDBTool::m_IDAlignDBTool

private

Definition at line 105 of file SiTrkAlignDBTool.h.

◆ m_idHelper

const AtlasDetectorID* InDet::SiTrkAlignDBTool::m_idHelper

private

Definition at line 91 of file SiTrkAlignDBTool.h.

◆ m_logStream

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

protectedinherited

logfile output stream

Definition at line 55 of file ITrkAlignDBTool.h.

◆ m_ntuple

TFile* InDet::SiTrkAlignDBTool::m_ntuple

private

Definition at line 112 of file SiTrkAlignDBTool.h.

◆ m_oldAlignFile

std::string InDet::SiTrkAlignDBTool::m_oldAlignFile

private

filename for the textfile with initial alignment constants

Definition at line 101 of file SiTrkAlignDBTool.h.

◆ m_oldGlobalFolderFile

std::string InDet::SiTrkAlignDBTool::m_oldGlobalFolderFile

private

filename for the textfile with initial GlobalFolder constants

Definition at line 102 of file SiTrkAlignDBTool.h.

◆ m_oldIBLDistFile

std::string InDet::SiTrkAlignDBTool::m_oldIBLDistFile

private

filename for the textfile with initial IBLDist params

Definition at line 100 of file SiTrkAlignDBTool.h.

◆ m_outGlobalFolderFile

std::string InDet::SiTrkAlignDBTool::m_outGlobalFolderFile

private

filename for the textfile with final GlobalFolder constants

Definition at line 97 of file SiTrkAlignDBTool.h.

◆ m_outIBLDistFile

std::string InDet::SiTrkAlignDBTool::m_outIBLDistFile

private

filename for the textfile with final IBLDist params

Definition at line 96 of file SiTrkAlignDBTool.h.

◆ m_outputAlignFile

std::string InDet::SiTrkAlignDBTool::m_outputAlignFile

private

filename for the textfile with final alignment constants

Definition at line 95 of file SiTrkAlignDBTool.h.

◆ m_pixelaAlignDBM

bool InDet::SiTrkAlignDBTool::m_pixelaAlignDBM

private

Definition at line 117 of file SiTrkAlignDBTool.h.

◆ m_pixelAlignLevel

int InDet::SiTrkAlignDBTool::m_pixelAlignLevel

private

Definition at line 114 of file SiTrkAlignDBTool.h.

◆ m_pixelAlignLevelBarrel

int InDet::SiTrkAlignDBTool::m_pixelAlignLevelBarrel

private

Definition at line 115 of file SiTrkAlignDBTool.h.

◆ m_pixelAlignLevelEndcaps

int InDet::SiTrkAlignDBTool::m_pixelAlignLevelEndcaps

private

Definition at line 116 of file SiTrkAlignDBTool.h.

◆ m_pixelGeoManager

ToolHandle<Trk::IGeometryManagerTool> InDet::SiTrkAlignDBTool::m_pixelGeoManager

private

Definition at line 108 of file SiTrkAlignDBTool.h.

◆ m_pixHelper

const PixelID* InDet::SiTrkAlignDBTool::m_pixHelper

private

pointer to PIX detector manager

Definition at line 89 of file SiTrkAlignDBTool.h.

◆ m_sctAlignLevel

int InDet::SiTrkAlignDBTool::m_sctAlignLevel

private

Definition at line 119 of file SiTrkAlignDBTool.h.

◆ m_sctAlignLevelBarrel

int InDet::SiTrkAlignDBTool::m_sctAlignLevelBarrel

private

Definition at line 120 of file SiTrkAlignDBTool.h.

◆ m_sctAlignLevelEndcaps

int InDet::SiTrkAlignDBTool::m_sctAlignLevelEndcaps

private

Definition at line 121 of file SiTrkAlignDBTool.h.

◆ m_sctGeoManager

ToolHandle<Trk::IGeometryManagerTool> InDet::SiTrkAlignDBTool::m_sctGeoManager

private

Definition at line 109 of file SiTrkAlignDBTool.h.

◆ m_sctHelper

const SCT_ID* InDet::SiTrkAlignDBTool::m_sctHelper

private

pointer to SCT detector manager

Definition at line 90 of file SiTrkAlignDBTool.h.

◆ m_siAlignLevel

int InDet::SiTrkAlignDBTool::m_siAlignLevel

private

Definition at line 123 of file SiTrkAlignDBTool.h.

◆ m_siGeoManager

ToolHandle<Trk::IGeometryManagerTool> InDet::SiTrkAlignDBTool::m_siGeoManager

private

Definition at line 110 of file SiTrkAlignDBTool.h.

◆ m_SQLiteTag

std::string InDet::SiTrkAlignDBTool::m_SQLiteTag

private

tag name for the ConditionsDB

Definition at line 94 of file SiTrkAlignDBTool.h.

◆ m_updateConstants

bool InDet::SiTrkAlignDBTool::m_updateConstants

private

flag whether to update the constants

Definition at line 103 of file SiTrkAlignDBTool.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.

◆ m_writeAsL16

bool InDet::SiTrkAlignDBTool::m_writeAsL16

private

Definition at line 131 of file SiTrkAlignDBTool.h.

◆ m_writeAsL2

bool InDet::SiTrkAlignDBTool::m_writeAsL2

private

Definition at line 130 of file SiTrkAlignDBTool.h.

◆ m_writeAsL3

bool InDet::SiTrkAlignDBTool::m_writeAsL3

private

Definition at line 129 of file SiTrkAlignDBTool.h.

◆ m_writeOldConstants

bool InDet::SiTrkAlignDBTool::m_writeOldConstants

private

flag whether to write out the initial constants into text file

Definition at line 99 of file SiTrkAlignDBTool.h.

◆ m_writeSQLFile

bool InDet::SiTrkAlignDBTool::m_writeSQLFile

private

flag wether to write out .db file or not (file is not needed in parallel jobs)

Definition at line 93 of file SiTrkAlignDBTool.h.

◆ m_writeTextFile

bool InDet::SiTrkAlignDBTool::m_writeTextFile

private

flag whether to write out the final constants into text file

Definition at line 98 of file SiTrkAlignDBTool.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ SiTrkAlignDBTool()

◆ ~SiTrkAlignDBTool()

Member Function Documentation

◆ checkPixelLevel()

◆ checkSCTLevel()

◆ checkSiLevel()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ finalize()

◆ GetL1L2fromL3Identifier()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ preSolve()

◆ printTransform()

◆ renounce()

◆ renounceArray()

◆ setLogStream()

◆ setNtuple()

◆ sysInitialize()

◆ sysStart()

◆ updateAsL16()

◆ updateAsL2()

◆ updateAsL3()

◆ updateDB()

◆ updateSiL0asL1()

◆ updateVHKA()

◆ writeAlignPar()

Member Data Documentation

◆ m_alignModuleTool

◆ m_detStore

◆ m_doPixel

◆ m_doSCT

◆ m_doSi

◆ m_evtStore

◆ m_IDAlignDBTool

◆ m_idHelper

◆ m_logStream

◆ m_ntuple

◆ m_oldAlignFile

◆ m_oldGlobalFolderFile

◆ m_oldIBLDistFile

◆ m_outGlobalFolderFile

◆ m_outIBLDistFile

◆ m_outputAlignFile

◆ m_pixelaAlignDBM

◆ m_pixelAlignLevel

◆ m_pixelAlignLevelBarrel

◆ m_pixelAlignLevelEndcaps

◆ m_pixelGeoManager

◆ m_pixHelper

◆ m_sctAlignLevel

◆ m_sctAlignLevelBarrel