The InDetAlignHitQualSelTool is to select good quality hits for alignment to build residuals with possible cuts on outlier hits, hits which are too large clusters, hits with large incidence angles, etc. More...

#include <InDetAlignHitQualSelTool.h>

Inheritance diagram for InDetAlignHitQualSelTool:

Collaboration diagram for InDetAlignHitQualSelTool:

Public Member Functions
	InDetAlignHitQualSelTool (const std::string &, const std::string &, const IInterface *)

virtual	~InDetAlignHitQualSelTool ()

virtual StatusCode	initialize ()

virtual StatusCode	finalize ()

const Trk::RIO_OnTrack *	getGoodHit (const Trk::TrackStateOnSurface *tsos) const
	main method: from a TrackStateOnSurface select a good hit cutting on outlier hits, hits with too many pixels/strips, hits with large incidence angles More...

bool	isGoodSiHit (const Trk::TrackStateOnSurface *tsos) const

bool	getGoodHole (const Trk::TrackStateOnSurface *tsos) const
	from a TrackStateOnSurface select a good hole in track cutting on large incidence angles only 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 ()
	for ToolHandle functionality More...

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

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
bool	isGangedPixel (const Trk::PrepRawData *prd) const
	check, whether cluster contains a ganged pixel More...

bool	isGoodClusterSize (const std::vector< Identifier > &idVec) const
	check, whether cluster size within limits of m_maxClusterSize More...

bool	isEdgeChannel (const std::vector< Identifier > &idVec) const
	check, whether the strip/pixel is an edge channel More...

double	incidAngle (const Trk::TrackParameters trkPar, const InDetDD::SiDetectorElement detEle) const
	calculate track incidence angle in local x-z frame More...

bool	isGoodAngle (const Trk::TrackParameters trkPar, const InDetDD::SiDetectorElement detEle) const
	check whether track incidence angle within limits of m_maxIncidAngle 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
bool	m_rejectOutliers
	reject hits labeled as outliers by the track fitter More...

int	m_maxClusterSize

bool	m_rejectEdgeChannels
	reject clusters containing edge channels More...

bool	m_rejectGangedPixels
	reject clusters containing ganged pixels More...

float	m_maxIncidAngle
	maximum incidence angle of a track (to which the hit belongs) on the Si-module. More...

bool	m_acceptIBLHits

bool	m_acceptPixelHits

bool	m_acceptSCTHits

const InDetDD::PixelDetectorManager *	m_PIXManager
	to get pixel phi and eta identifiers More...

const PixelID *	m_pixelid
	Pixel id helper. More...

const SCT_ID *	m_sctID
	Pixel id helper. More...

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

The InDetAlignHitQualSelTool is to select good quality hits for alignment to build residuals with possible cuts on outlier hits, hits which are too large clusters, hits with large incidence angles, etc.

The InDetAlignHitQualSelTool is to select good quality hits to build residuals for alignment. The following cuts are possible: outlier hits, hits which have too large clusters, hits with large incidence angles, ganged pixels, edge channels.

Author: Oleg Brandt http://consult.cern.ch/xwho

Definition at line 34 of file InDetAlignHitQualSelTool.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

◆ InDetAlignHitQualSelTool()

InDetAlignHitQualSelTool::InDetAlignHitQualSelTool	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 21 of file InDetAlignHitQualSelTool.cxx.

   : AthAlgTool(t,n,p)
   , m_rejectOutliers( true )
   , m_maxClusterSize( 5 )
   , m_rejectEdgeChannels( true )
   , m_rejectGangedPixels( false )
   , m_maxIncidAngle( 0.8 ) 
   , m_acceptIBLHits( true )
   , m_acceptPixelHits( true )
   , m_acceptSCTHits( true )
   , m_PIXManager{}
   , m_pixelid{}
   , m_sctID{} 
 {
   declareInterface<IInDetAlignHitQualSelTool>(this) ;
   declareProperty( "RejectOutliers",         m_rejectOutliers     ) ;
   declareProperty( "MaxClusterSize",         m_maxClusterSize     ) ;
   declareProperty( "RejectEdgeChannels",     m_rejectEdgeChannels ) ;
   declareProperty( "RejectGangedPixels",     m_rejectGangedPixels ) ;
   declareProperty( "MaxIncidAngle",          m_maxIncidAngle      ) ;
   declareProperty( "AcceptIBLHits",          m_acceptIBLHits    ) ;
   declareProperty( "AcceptPixelHits",        m_acceptPixelHits    ) ;
   declareProperty( "AcceptSCTHits",          m_acceptSCTHits    ) ;
 }

◆ ~InDetAlignHitQualSelTool()

InDetAlignHitQualSelTool::~InDetAlignHitQualSelTool ( )

virtual

Definition at line 50 of file InDetAlignHitQualSelTool.cxx.

50 {}

Member Function Documentation

◆ 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 InDetAlignHitQualSelTool::finalize ( )

virtual

Definition at line 66 of file InDetAlignHitQualSelTool.cxx.

                                               {
   ATH_MSG_DEBUG( "finalize() successful in " << name() ) ; 
   return AlgTool::finalize() ;
 }

◆ getGoodHit()

const Trk::RIO_OnTrack * InDetAlignHitQualSelTool::getGoodHit ( const Trk::TrackStateOnSurface * tsos ) const

virtual

main method: from a TrackStateOnSurface select a good hit cutting on outlier hits, hits with too many pixels/strips, hits with large incidence angles

Implements IInDetAlignHitQualSelTool.

Definition at line 71 of file InDetAlignHitQualSelTool.cxx.

                                                                                                      {
   ATH_MSG_DEBUG( "** getGoodHit ** Dealing with a new tsos ** START ** " ) ;
  
   if( tsos == nullptr ) {
     ATH_MSG_ERROR( "0 pointer passed for TSOS!" ) ;
     return nullptr ;
   }
   if( !tsos->type(Trk::TrackStateOnSurface::Measurement) ) {
     ATH_MSG_DEBUG( "not a hit, cast to MeasurementBase will fail, so reject" ) ;
     return nullptr ;
   }
   if( m_rejectOutliers && tsos->type(Trk::TrackStateOnSurface::Outlier) ) {
     ATH_MSG_DEBUG( "outlier, reject" ) ;
     return nullptr ;
   }
   const Trk::MeasurementBase* measBase = tsos->measurementOnTrack() ;
   if( measBase == nullptr) {
     ATH_MSG_DEBUG( "tsos->measurementOnTrack() returned 0 pointer" ) ;
     return nullptr ;
   }
   
   const Trk::RIO_OnTrack* hit = dynamic_cast <const Trk::RIO_OnTrack*>( measBase ) ;
   if( hit == nullptr ) {
     ATH_MSG_DEBUG( "dynamic_cast <const Trk::RIO_OnTrack*>( measBase ) returned 0 pointer" ) ;
     return nullptr ;
   }
  
   const Trk::PrepRawData* prd = hit->prepRawData() ;
   if( prd == nullptr ) {
     ATH_MSG_WARNING( "hit->prepRawData() method failed" ) ;
     return nullptr ;
   }
   
   
   if( m_rejectGangedPixels && isGangedPixel( prd ) ) return nullptr ;
   const vector<Identifier> &idVec = prd->rdoList() ;
  
   // cut on the cluster size
   if( m_maxClusterSize > 0 && !isGoodClusterSize( idVec ) ) return nullptr ;
  
   // cut on edge channels
   if( m_rejectEdgeChannels && isEdgeChannel( idVec ) ) return nullptr ;
   // cut on the track incidence angle alpha
   const Trk::TrackParameters* trkPar = tsos->trackParameters() ;
   if( trkPar == nullptr ) {
     ATH_MSG_WARNING( "tsos->trackParameters() returned 0 pointer" ) ;
     return nullptr ;
   }
   const InDetDD::SiDetectorElement *detEle
     = dynamic_cast<const InDetDD::SiDetectorElement*>( hit->detectorElement() ) ;
   if( detEle == nullptr ) {
     ATH_MSG_WARNING( "hit cast to SiDetectorElement returned 0 pointer" ) ;
     return nullptr ;
   }
   if( !isGoodAngle( trkPar, detEle ) ) return nullptr ;  
   return hit ;
 }

◆ getGoodHole()

bool InDetAlignHitQualSelTool::getGoodHole ( const Trk::TrackStateOnSurface * tsos ) const

virtual

from a TrackStateOnSurface select a good hole in track cutting on large incidence angles only

Implements IInDetAlignHitQualSelTool.

Definition at line 233 of file InDetAlignHitQualSelTool.cxx.

                                                                                      {
   if( tsos == nullptr ) {
     ATH_MSG_ERROR( "0 pointer passed for TSOS!" ) ;
     return false ;
   }
   if( !tsos->type(Trk::TrackStateOnSurface::Hole) ) {
     ATH_MSG_DEBUG( "This is not a hole, reject" ) ;
     return false ;
   }
   // for holes only cut on the track incidence angle alpha
   const Trk::TrackParameters* trkPar = tsos->trackParameters() ;
   if( trkPar == nullptr ) {
     ATH_MSG_WARNING( "tsos->trackParameters() returned 0 pointer" ) ;
     return false ;
   }
   const InDetDD::SiDetectorElement *detEle = dynamic_cast<const InDetDD::SiDetectorElement*>( 
     tsos->trackParameters()->associatedSurface().associatedDetectorElement() ) ;
   if( detEle == nullptr ) {
     ATH_MSG_WARNING( "hole cast to SiDetectorElement returned 0 pointer" ) ;
     return false ;
   }
   if( !isGoodAngle( trkPar, detEle ) ) return false ;
   return true;
 }

◆ incidAngle()

double InDetAlignHitQualSelTool::incidAngle	(	const Trk::TrackParameters *	trkPar,
		const InDetDD::SiDetectorElement *	detEle
	)		const

private

calculate track incidence angle in local x-z frame

< local x axis in global frame

< local z axis in global frame

< scalar product

Definition at line 323 of file InDetAlignHitQualSelTool.cxx.

                                                          {
   Amg::Vector3D trkDir       = trkPar->momentum() ;
   const Amg::Vector3D& detElePhi    = detEle->phiAxis() ; 
   const Amg::Vector3D& detEleNormal = detEle->normal() ; 
   double trkDotPhi = trkDir.dot( detElePhi ) ; 
   double trkDotNormal = trkDir.dot( detEleNormal ) ;
   double trkIncidAngle = atan( trkDotPhi/trkDotNormal ) ;
   ATH_MSG_DEBUG( "trkIncidAngle = " << trkIncidAngle ) ;
   return trkIncidAngle ;
 }

◆ initialize()

StatusCode InDetAlignHitQualSelTool::initialize ( )

virtual

Definition at line 53 of file InDetAlignHitQualSelTool.cxx.

                                                 {
   StatusCode sc = AlgTool::initialize() ;
   if( sc.isFailure() ) return sc ;
   // get DetectorStore service
   ATH_CHECK( detStore().retrieve()) ;
   ATH_CHECK(detStore()->retrieve(m_sctID, "SCT_ID"));
   ATH_CHECK(detStore()->retrieve(m_pixelid, "PixelID"));
   // get pixel manager
   ATH_CHECK(detStore()->retrieve( m_PIXManager, "Pixel" )) ;
   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 & IInDetAlignHitQualSelTool::interfaceID ( )

inlinestaticinherited

for ToolHandle functionality

Definition at line 42 of file IInDetAlignHitQualSelTool.h.

                                                                  {
   return IID_INDETALIGN_IInDetAlignHitQualSelTool ; 
 }

◆ isEdgeChannel()

bool InDetAlignHitQualSelTool::isEdgeChannel ( const std::vector< Identifier > & idVec ) const

private

check, whether the strip/pixel is an edge channel

Definition at line 284 of file InDetAlignHitQualSelTool.cxx.

                                                                                     {
   for( unsigned int i=0, i_max=idVec.size() ; i!=i_max ; ++i ) {
     if( m_sctID->is_sct(idVec[i]) ) {
       int stripId = m_sctID->strip(idVec[i]) ;
       if( stripId == 0 || stripId == 767 ) {
         ATH_MSG_DEBUG( " SCT strip " << i << " with id " << stripId << " is an edge channel " ) ;
         return true ;
       }
       if( stripId < 0 || stripId > 767 ) {
         ATH_MSG_FATAL( " WRONG DETECTOR INFORMATION " ) ;
       }
     }
     if( m_PIXManager->identifierBelongs(idVec[i]) ) {
       int pixelIdPhi = m_pixelid->phi_index(idVec[i]) ;
       int pixelIdEta = m_pixelid->eta_index(idVec[i]) ;
       if( pixelIdPhi == 0 || pixelIdPhi == 327 || pixelIdEta == 0 || pixelIdEta == 143 ) {
         ATH_MSG_DEBUG( " pixel hit " << i << " with idPhi " << pixelIdPhi << " and idEta " << pixelIdEta << " is an edge channel " ) ;
         return true ;
       }
       if( pixelIdPhi < 0 || pixelIdPhi > 327 || pixelIdEta < 0 || pixelIdEta > 143 ) {
         ATH_MSG_FATAL( " WRONG DETECTOR INFORMATION " ) ;
       }
     }
   }
   return false ;
 }

◆ isGangedPixel()

bool InDetAlignHitQualSelTool::isGangedPixel ( const Trk::PrepRawData * prd ) const

private

check, whether cluster contains a ganged pixel

< cut only if m_maxClusterSize set

Definition at line 259 of file InDetAlignHitQualSelTool.cxx.

                                                                               {
   const InDet::SiCluster* cluster = dynamic_cast<const InDet::SiCluster*>( prd ) ;
   if( cluster == nullptr ) {
     ATH_MSG_WARNING( "dynamic_cast<const InDet::SiCluster*>( prd ) failed!" ) ;
     return false ;
   }
   if( cluster->gangedPixel() ) { 
     ATH_MSG_DEBUG( "cluster contains a ganged pixel, reject" ) ;
     return true ;
   }
   return false ;
 }

◆ isGoodAngle()

bool InDetAlignHitQualSelTool::isGoodAngle	(	const Trk::TrackParameters *	trkPar,
		const InDetDD::SiDetectorElement *	detEle
	)		const

private

check whether track incidence angle within limits of m_maxIncidAngle

Definition at line 311 of file InDetAlignHitQualSelTool.cxx.

                                                         {
   const double trkIncidAngle = incidAngle( trkPar, detEle ) ;
   if( std::abs(trkIncidAngle) > m_maxIncidAngle ) {
     ATH_MSG_DEBUG( "trkIncidAngle = |" << trkIncidAngle << "| > " << m_maxIncidAngle << ", reject" ) ;
     return false ;
   }
   return true;
 }

◆ isGoodClusterSize()

bool InDetAlignHitQualSelTool::isGoodClusterSize ( const std::vector< Identifier > & idVec ) const

private

check, whether cluster size within limits of m_maxClusterSize

< cut only if m_maxClusterSize set

Definition at line 273 of file InDetAlignHitQualSelTool.cxx.

                                                                                            {
   int clusterSize = idVec.size() ;
   ATH_MSG_DEBUG( "clusterSize = " << clusterSize ) ;
   if( clusterSize > m_maxClusterSize ) { 
     ATH_MSG_DEBUG( "clusterSize = " << clusterSize << " > " << m_maxClusterSize << ", reject" ) ;
     return false ;
   }
   return true ;
 }

◆ isGoodSiHit()

bool InDetAlignHitQualSelTool::isGoodSiHit ( const Trk::TrackStateOnSurface * tsos ) const

virtual

Implements IInDetAlignHitQualSelTool.

Definition at line 129 of file InDetAlignHitQualSelTool.cxx.

                                                                                      {
   bool isSiliconHit = false;
   bool isPixelHit = false;
   bool isIBLHit = false;
  
   if( tsos == nullptr ) {
     ATH_MSG_ERROR( "0 pointer passed for TSOS!" ) ;
     return false ;
   }
   if( !tsos->type(Trk::TrackStateOnSurface::Measurement) ) {
     ATH_MSG_DEBUG( "not a hit, cast to MeasurementBase will fail --> keep it anyway" ) ;
     return false ;
   }
  
   if( m_rejectOutliers && tsos->type(Trk::TrackStateOnSurface::Outlier) ) {
     ATH_MSG_DEBUG( "outlier --> keep it" ) ;
     return false ;
   }
   const Trk::MeasurementBase* measBase = tsos->measurementOnTrack() ;
   if( measBase == nullptr) {
     ATH_MSG_DEBUG( "tsos->measurementOnTrack() returned 0 pointer" ) ;
     return false ;
   }
   
   const Trk::RIO_OnTrack* hit = dynamic_cast <const Trk::RIO_OnTrack*>( measBase ) ;
   if( hit == nullptr ) {
     ATH_MSG_DEBUG( "dynamic_cast <const Trk::RIO_OnTrack*>( measBase ) returned 0 pointer" ) ;
     return false ;
   }
  
   const Trk::PrepRawData* prd = hit->prepRawData() ;
   if( prd == nullptr ) {
     ATH_MSG_WARNING( "hit->prepRawData() method failed" ) ;
     return false ;
   }
   
   
   const Identifier & hitId = hit->identify();
   if (m_sctID->is_sct(hitId)) {
     ATH_MSG_DEBUG( " this is a SCT hit - SCT - SCT - ");
     isSiliconHit = true;
   }
   if (m_pixelid->is_pixel(hitId)) {
     ATH_MSG_DEBUG( " this is a PIX hit - PIX - PIX - ");
     isSiliconHit = true;
     isPixelHit = true; // assume that is pixel hit
     // but check if it is IBL
     if (m_pixelid->layer_disk(hitId) == 0 && m_pixelid->barrel_ec(hitId) == 0 ) {isIBLHit = true; isPixelHit = false;}
   }
   
   if (!isSiliconHit) {
     ATH_MSG_DEBUG( "This is not a silicon hit. Keep it as good" ) ;
     return true;
   }
  
   // accept IBL hits ?
   if (!m_acceptIBLHits && isIBLHit) {
     ATH_MSG_INFO( "this is an IBL hit --> user wants to drop it" ) ;
     return false;
   }
  
   // accept pixel hits ?
   if (!m_acceptPixelHits && isPixelHit) {
     ATH_MSG_INFO( "this is a pixel hit --> user wants to drop it" ) ;
     return false;
   }
  
   // accept SCT hits ?
   if (!m_acceptSCTHits) {
     if (isSiliconHit && m_sctID->is_sct(hitId)) {
       ATH_MSG_INFO( "this is a SCT hit --> user wants to drop it" ) ;
       return false;
     }
   }
  
   if( m_rejectGangedPixels && isGangedPixel( prd ) ) return false ;
   const vector<Identifier> &idVec = prd->rdoList() ;
  
   // cut on the cluster size
   if( m_maxClusterSize > 0 && !isGoodClusterSize( idVec ) ) return false ;
  
   // cut on edge channels
   if( m_rejectEdgeChannels && isEdgeChannel( idVec ) ) return false ;
  
   // cut on the track incidence angle alpha
   const Trk::TrackParameters* trkPar = tsos->trackParameters() ;
   if( trkPar == nullptr ) {
     ATH_MSG_WARNING( "tsos->trackParameters() returned 0 pointer" ) ;
     return false ;
   }
   // incidence angle
   const InDetDD::SiDetectorElement *detEle
     = dynamic_cast<const InDetDD::SiDetectorElement*>( hit->detectorElement() ) ;
   if( detEle == nullptr ) {
     ATH_MSG_WARNING( "hit cast to SiDetectorElement returned 0 pointer" ) ;
     return false ;
   }
   if( !isGoodAngle( trkPar, detEle ) ) return false ;  
  
   return true ;
 }

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

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

Reimplemented in 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.

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_acceptIBLHits

bool InDetAlignHitQualSelTool::m_acceptIBLHits

private

Definition at line 85 of file InDetAlignHitQualSelTool.h.

◆ m_acceptPixelHits

bool InDetAlignHitQualSelTool::m_acceptPixelHits

private

Definition at line 89 of file InDetAlignHitQualSelTool.h.

◆ m_acceptSCTHits

bool InDetAlignHitQualSelTool::m_acceptSCTHits

private

Definition at line 93 of file InDetAlignHitQualSelTool.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_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_maxClusterSize

int InDetAlignHitQualSelTool::m_maxClusterSize

private

Definition at line 72 of file InDetAlignHitQualSelTool.h.

◆ m_maxIncidAngle

float InDetAlignHitQualSelTool::m_maxIncidAngle

private

maximum incidence angle of a track (to which the hit belongs) on the Si-module.

It is caculated in the local xz frame to ensure that the only the angle component perpendicular to the strips is considered. It is defined w/r/t the local z-axis.

Definition at line 81 of file InDetAlignHitQualSelTool.h.

◆ m_pixelid

const PixelID* InDetAlignHitQualSelTool::m_pixelid

private

Pixel id helper.

Definition at line 98 of file InDetAlignHitQualSelTool.h.

◆ m_PIXManager

const InDetDD::PixelDetectorManager* InDetAlignHitQualSelTool::m_PIXManager

private

to get pixel phi and eta identifiers

Definition at line 97 of file InDetAlignHitQualSelTool.h.

◆ m_rejectEdgeChannels

bool InDetAlignHitQualSelTool::m_rejectEdgeChannels

private

reject clusters containing edge channels

Definition at line 74 of file InDetAlignHitQualSelTool.h.

◆ m_rejectGangedPixels

bool InDetAlignHitQualSelTool::m_rejectGangedPixels

private

reject clusters containing ganged pixels

Definition at line 76 of file InDetAlignHitQualSelTool.h.

◆ m_rejectOutliers

bool InDetAlignHitQualSelTool::m_rejectOutliers

private

reject hits labeled as outliers by the track fitter

Definition at line 69 of file InDetAlignHitQualSelTool.h.

◆ m_sctID

const SCT_ID* InDetAlignHitQualSelTool::m_sctID

private

Pixel id helper.

Definition at line 99 of file InDetAlignHitQualSelTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ InDetAlignHitQualSelTool()

◆ ~InDetAlignHitQualSelTool()

Member Function Documentation

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

◆ getGoodHit()

◆ getGoodHole()

◆ incidAngle()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ isEdgeChannel()

◆ isGangedPixel()

◆ isGoodAngle()

◆ isGoodClusterSize()

◆ isGoodSiHit()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_acceptIBLHits

◆ m_acceptPixelHits

◆ m_acceptSCTHits

◆ m_detStore

◆ m_evtStore

◆ m_maxClusterSize

◆ m_maxIncidAngle

◆ m_pixelid

◆ m_PIXManager

◆ m_rejectEdgeChannels

◆ m_rejectGangedPixels

◆ m_rejectOutliers

◆ m_sctID

◆ m_varHandleArraysDeclared

◆ m_vhka