#include <SiSPSeededTrackFinder.h>

Inheritance diagram for InDet::SiSPSeededTrackFinder:

Collaboration diagram for InDet::SiSPSeededTrackFinder:

Classes
class	Counter_t

Public Member Functions
virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...

virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...

virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm. More...

virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...

virtual bool	filterPassed (const EventContext &ctx) const

virtual void	setFilterPassed (bool state, const EventContext &ctx) const

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

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

Print internal tool parameters and status
MsgStream &	dump (MSG::Level lvl, const SiSPSeededTrackFinder::Counter_t *) const

MsgStream &	dumptools (MsgStream &out) const

MsgStream &	dumpevent (MsgStream &out, const SiSPSeededTrackFinder::Counter_t &counter) const

bool	isGoodEvent (const EventContext &ctx) const
	EventContext is used to specify which event. More...

void	filterSharedTracks (std::multimap< double, Trk::Track * > &scoredTracks) const
	cleans up the collection of quality filtered tracks. More...

void	filterSharedTracksFast (std::multimap< double, Trk::Track * > &scoredTracks) const

void	fillZHistogram (const Trk::Track *Tr, const Trk::PerigeeSurface &beamlinePerigee, std::vector< int > &numberWeightedhistogram, std::vector< double > &zWeightedHistogram, std::vector< double > &ptWeightedHistogram) const
	fills three z0 histograms (non-weighted, weighted by z, and weighted by pt) with the track z at the beam line estimated using the innermost measurement. More...

void	findZvertex (std::list< Trk::Vertex > &vertexList, std::pair< double, double > &zBoundaries, const std::vector< int > &numberWeightedhistogram, const std::vector< double > &zWeightedHistogram, const std::vector< double > &ptWeightedHistogram) const
	estimates a set of vertex positions and a z interval for the second track finding pass using the input histograms populated using fillZHistogram. More...

StatusCode	oldStrategy (const EventContext &ctx) const
	this method performs the track finding using the old strategy More...

StatusCode	newStrategy (const EventContext &ctx) const
	this method performs the track finding using the new strategy More...

StatusCode	itkFastTrackingStrategy (const EventContext &ctx) const
	EventContext is used to specify which event. More...

StatusCode	itkConvStrategy (const EventContext &ctx) const
	EventContext is used to specify which event. More...

void	magneticFieldInit ()

bool	passEtaDepCuts (const Trk::Track *track, int nClusters, int nFreeClusters, int nPixels) const

Protected Attributes
Flags to configure SiSPSeededTrackFinder
BooleanProperty	m_useZvertexTool {this, "useZvertexTool", true}

BooleanProperty	m_writeHolesFromPattern {this, "writeHolesFromPattern", false,"Flag to activate writing hole info from the pattern recognition"}

BooleanProperty	m_useMBTS {this, "useMBTSTimeDiff", false}

BooleanProperty	m_useNewStrategy {this, "useNewStrategy", false}

BooleanProperty	m_useZBoundaryFinding {this, "useZBoundFinding", false}

BooleanProperty	m_ITKGeometry {this, "ITKGeometry", false, "Flag to specify if it is ITK geometry"}

BooleanProperty	m_useITkConvSeeded {this, "useITkConvSeeded", false, "ITk EM-seeded conversion reco"}

BooleanProperty	m_doFastTracking {this, "doFastTracking", false, "ITk fast tracking reco"}

IntegerProperty	m_maxNumberSeeds {this, "maxNumberSeeds", 3000000, "Max. number used seeds"}

BooleanProperty	m_alwaysProtectAgainstBusyEvent {this, "alwaysProtectAgainstBusyEvents", false, "Flag that enables checking always maxNumberPIXsp and maxNumberSCTsp"}

IntegerProperty	m_maxPIXsp {this, "maxNumberPIXsp", 150000, "Max. number pixels space points"}

IntegerProperty	m_maxSCTsp {this, "maxNumberSCTsp", 500000, "Max. number sct space points"}

IntegerProperty	m_nfreeCut {this, "FreeClustersCut", 1, "Min number free clusters"}

IntegerProperty	m_histsize {this, "HistSize", 1400}

IntegerProperty	m_nvertex {this, "maxVertices", 4}

DoubleProperty	m_zcut {this, "Zcut", 350.}

DoubleProperty	m_deltaEta {this, "dEtaCaloRoI", .1}

DoubleProperty	m_deltaPhi {this, "dPhiCaloRoI", .25}

DoubleProperty	m_deltaZ {this, "dZCaloRoI", 300.}

StringProperty	m_fieldmode {this, "MagneticFieldMode", "MapSolenoid"}

Data members for new strategy reconstruction
double	m_pTcut {500}

double	m_imcut {2.}

double	m_zstep {0.}

Data handles for StoreGate access in AthenaMT
SG::ReadHandleKey< xAOD::EventInfo >	m_evtKey {this, "EventInfoKey", "EventInfo"}

SG::ReadHandleKey< SpacePointContainer >	m_SpacePointsSCTKey {this, "SpacePointsSCTName", "SCT_SpacePoints"}

SG::ReadHandleKey< SpacePointContainer >	m_SpacePointsPixelKey {this, "SpacePointsPixelName", "PixelSpacePoints"}

SG::ReadHandleKey< Trk::PRDtoTrackMap >	m_prdToTrackMap {this,"PRDtoTrackMap",""}

SG::ReadHandleKey< ROIPhiRZContainer >	m_caloClusterROIKey {this, "EMROIPhiRZContainer", "InDetCaloClusterROIPhiRZ15GeVUnordered", "Name of the calo cluster ROIs in Phi,R,Z parameterization"}

SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"}

SG::ReadDecorHandleKey< xAOD::EventInfo >	m_mbtsKey {this, "MBTSKey", "EventInfo.MBTSBackground"}

SG::WriteHandleKey< TrackCollection >	m_outputTracksKey {this, "TracksLocation", "SiSPSeededTracks"}

Tool handles
ToolHandle< ISiSpacePointsSeedMaker >	m_seedsmaker {this, "SeedsTool", "InDet::SiSpacePointsSeedMaker_ATLxk/InDetSpSeedsMaker", "Space poins seed maker"}

ToolHandle< ISiZvertexMaker >	m_zvertexmaker {this, "ZvertexTool", "InDet::SiZvertexMaker_xk/InDetSiZvertexMaker", "Space poins z-vertex maker"}

ToolHandle< ISiTrackMaker >	m_trackmaker {this, "TrackTool", "InDet::SiTrackMaker_xk/InDetSiTrackMaker", "Track maker"}

PublicToolHandle< Trk::IPatternParametersPropagator >	m_proptool {this, "PropagatorTool", "Trk::RungeKuttaPropagator/InDetPropagator"}

ToolHandle< Trk::IExtendedTrackSummaryTool >	m_trackSummaryTool {this, "TrackSummaryTool", "InDetTrackSummaryToolNoHoleSearch"}

ToolHandle< IRegSelTool >	m_regsel_strip { this, "RegSelTool_Strip", "RegSelTool/RegSelTool_SCT" }

Magnetic field properties
Trk::MagneticFieldProperties	m_fieldprop

Service handles
ServiceHandle< IInDetEtaDependentCutsSvc >	m_etaDependentCutsSvc {this, "InDetEtaDependentCutsSvc", "InDetEtaDependentCutsSvc"}
	service to get cut values depending on different variable More...

Counters
Counter_t m_counterTotal	ATLAS_THREAD_SAFE {}

std::atomic_int	m_neventsTotal {0}
	Number events. More...

std::atomic_int	m_neventsTotalV {0}
	Number events. More...

std::atomic_int	m_problemsTotal {0}
	Number events with number seeds > maxNumber. More...

std::atomic_int	m_problemsTotalV {0}
	Number events with number seeds > maxNumber. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
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
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. 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

Standard Algotithm methods
enum	ECounter { kNSeeds, kNTracks, kNCounter }
	enums for Counter_t More...

	SiSPSeededTrackFinder (const std::string &name, ISvcLocator *pSvcLocator)
	Constructor. More...

virtual	~SiSPSeededTrackFinder ()=default

virtual StatusCode	initialize () override
	Initialisation. More...

virtual StatusCode	execute (const EventContext &ctx) const override
	Execute. More...

virtual StatusCode	finalize () override

virtual bool	isClonable () const override
	Make this algorithm clonable. More...

Detailed Description

Class-algorithm for track finding in Pixels and SCT initiated by space points seeds

Definition at line 48 of file SiSPSeededTrackFinder.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ ECounter

enum InDet::SiSPSeededTrackFinder::ECounter

protected

enums for Counter_t

Enumerator
kNSeeds
kNTracks
kNCounter

Definition at line 78 of file SiSPSeededTrackFinder.h.

78 {kNSeeds, kNTracks, kNCounter};

Constructor & Destructor Documentation

◆ SiSPSeededTrackFinder()

InDet::SiSPSeededTrackFinder::SiSPSeededTrackFinder	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Constructor.

Definition at line 64 of file SiSPSeededTrackFinder.cxx.

                                                  : AthReentrantAlgorithm(name, pSvcLocator)
 {
 }

◆ ~SiSPSeededTrackFinder()

virtual InDet::SiSPSeededTrackFinder::~SiSPSeededTrackFinder ( )

virtualdefault

Member Function Documentation

◆ cardinality()

unsigned int AthReentrantAlgorithm::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 55 of file AthReentrantAlgorithm.cxx.

 {
   return 0;
 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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; }

◆ dump()

MsgStream & InDet::SiSPSeededTrackFinder::dump	(	MSG::Level	lvl,
		const SiSPSeededTrackFinder::Counter_t *	counter
	)		const

protected

Definition at line 724 of file SiSPSeededTrackFinder.cxx.

 {
   msg(assign_level) <<std::endl;
   MsgStream& out_msg=msg();
   if (counter) dumpevent(out_msg ,*counter);
   else dumptools(out_msg);
   out_msg << endmsg;
   return out_msg;
 }

◆ dumpevent()

MsgStream & InDet::SiSPSeededTrackFinder::dumpevent	(	MsgStream &	out,
		const SiSPSeededTrackFinder::Counter_t &	counter
	)		const

protected

Definition at line 778 of file SiSPSeededTrackFinder.cxx.

 {
   out<<"|-------------------------------------------------------------------";
   out<<"---------------------------------|"
      <<std::endl;
   out<<"|  Investigated "
      <<std::setw(9)<<counter[kNSeeds]<<" space points seeds and found ";
   out<<std::setw(9)<<counter[kNTracks];
   if      (m_ITKGeometry        ) out<<" tracks using new strategy for ITK   |"<<std::endl; 
   else if (m_useNewStrategy     ) out<<" tracks using new strategy ("<<std::setw(2)<<  m_nvertex <<")      |"<<std::endl;
   else if (m_useZBoundaryFinding) out<<" tracks using old strategy with Zb   |"<<std::endl;
   else                            out<<" tracks using old strategy           |"<<std::endl;
  
   out<<"|-------------------------------------------------------------------";
   out<<"---------------------------------|"
      <<std::endl;
   if (m_problemsTotal or m_problemsTotalV) {
     out<<"|  Events       "
        <<std::setw(7)<<m_neventsTotal   <<"            without Z-vertz   "
        <<std::setw(7)<<m_neventsTotalV  <<"            with    Z-vertex         |"
        <<std::endl;
     out<<"|  Problems     "
        <<std::setw(7)<<m_problemsTotal  <<"            without Z-vertz   "
        <<std::setw(7)<<m_problemsTotalV <<"            with    Z-vertex         |"
        <<std::endl;
     out<<"|-------------------------------------------------------------------";
     out<<"-----------------------------|"
        <<std::endl;
   }
   return out;
 }

◆ dumptools()

MsgStream & InDet::SiSPSeededTrackFinder::dumptools ( MsgStream & out ) const

protected

Definition at line 738 of file SiSPSeededTrackFinder.cxx.

 {
   int n = 65-m_zvertexmaker.type().size();
   std::string s1; for (int i=0; i<n; ++i) s1.append(" "); s1.append("|");
   n     = 65-m_seedsmaker.type().size();
   std::string s2; for (int i=0; i<n; ++i) s2.append(" "); s2.append("|");
   n     = 65-m_trackmaker.type().size();
   std::string s3; for (int i=0; i<n; ++i) s3.append(" "); s3.append("|");
   n     = 65-m_outputTracksKey.key().size();
   std::string s4; for (int i=0; i<n; ++i) s4.append(" "); s4.append("|");
  
   std::string s5; 
   if (m_useZvertexTool) s5= "Yes"; else s5 = "No";
   n     = 65-s5.size(); for (int i=0; i<n; ++i) s5.append(" "); s5.append("|");
   
   out<<"|----------------------------------------------------------------"
      <<"----------------------------------------------------|"
      <<std::endl;
   out<<"| Use primary vertices z-coordinates      finding?| "<<s5
      <<std::endl;
   if (m_useZvertexTool) {
     out<<"| Tool for primary vertices z-coordinates finding | "<<m_zvertexmaker.type()<<s1
        <<std::endl;
   }
   out<<"| Tool for space points seeds             finding | "<<m_seedsmaker.type()<<s2
      <<std::endl;
   out<<"| Tool for space points seeded track      finding | "<<m_trackmaker.type()<<s3
      <<std::endl;
   out<<"| Location of output tracks                       | "<<m_outputTracksKey.key()<<s4
      <<std::endl;
   out<<"|----------------------------------------------------------------"
      <<"----------------------------------------------------|"
      <<std::endl;
   return out;
 }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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; }

◆ execute()

StatusCode InDet::SiSPSeededTrackFinder::execute ( const EventContext & ctx ) const

overridevirtual

Execute.

Note that !m_useNewStrategy alone is NOT sufficient to trigger the oldStrategy call here! For example, run-3 central offline Si tracking has m_useNewStrategy=false, but m_useZBoundaryFinding true --> newStrategy

Definition at line 141 of file SiSPSeededTrackFinder.cxx.

 { 
   if (m_ITKGeometry and m_doFastTracking) return itkFastTrackingStrategy(ctx);
   else if (m_useITkConvSeeded) return itkConvStrategy(ctx);
   else if (not m_useNewStrategy and not m_useZBoundaryFinding and not m_ITKGeometry) {
     return oldStrategy(ctx);
   }
   return newStrategy(ctx);
 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthReentrantAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 79 of file AthReentrantAlgorithm.cxx.

 {
   // If we didn't find any symlinks to add, just return the collection
   // from the base class.  Otherwise, return the extended collection.
   if (!m_extendedExtraObjects.empty()) {
     return m_extendedExtraObjects;
   }
   return Algorithm::extraOutputDeps();
 }

◆ fillZHistogram()

void InDet::SiSPSeededTrackFinder::fillZHistogram	(	const Trk::Track *	Tr,
		const Trk::PerigeeSurface &	beamlinePerigee,
		std::vector< int > &	numberWeightedhistogram,
		std::vector< double > &	zWeightedHistogram,
		std::vector< double > &	ptWeightedHistogram
	)		const

protected

fills three z0 histograms (non-weighted, weighted by z, and weighted by pt) with the track z at the beam line estimated using the innermost measurement.

the first two parameters are input, the other three output.

Parameters

[in]	Tr	Track candidate to fill
[in]	beamlinePerigee	Perigee surface corresponding to the beam spot
[out]	numberWeightedhistogram	vector representing a histogram in z, counting the tracks per bin
[out]	zWeightedHistogram	vector representing a histogram in z, counting the tracks per bin weighted by their z values
[out]	ptWeightedHistogram	vector representing a histogram in z, counting the tracks per bin weighted by their pt values

only take into accounts tracks with a hit inside r < 60mm

propagate from innermost hit to beam spot

determine bin number - m_zstep is the inverse bin width, where the histo axis extends from -m_zcut to +m_zcut

fill histograms if we are not in the over/underflow

simple z histogram, counting tracks per z

z weighted histogram binned in z - used for vertex z calculation

pt weighted histogram binned in z - used for vertex sumpt calculation

Definition at line 973 of file SiSPSeededTrackFinder.cxx.

 {
   
   if (Tr->measurementsOnTrack()->size() < 10) return;
     
   const Trk::TrackParameters* paramsAtFirstSurface = Tr->trackStateOnSurfaces()->front()->trackParameters();
   Amg::Vector3D              position = paramsAtFirstSurface->position()          ;
   Amg::Vector3D              momentum = paramsAtFirstSurface->momentum()          ;
   
   constexpr double rSquare_max_forZHisto = 60.*60.; 
   if (position.x()*position.x()+position.y()*position.y() >= rSquare_max_forZHisto) return;
  
   double pT = sqrt(momentum.x()*momentum.x()+momentum.y()*momentum.y());
   if (pT < m_pTcut) return;
  
   Trk::PatternTrackParameters TP;
   if (not TP.production(paramsAtFirstSurface)) return;
       
   double step;
   if (not m_proptool->propagate(Gaudi::Hive::currentContext(),
                                 TP, beamPosPerigee, TP, Trk::anyDirection, m_fieldprop, step, Trk::pion)) return;
           
   const AmgVector(5)& parsAtBeamSpot = TP.parameters();
   if (std::abs(parsAtBeamSpot[0]) > m_imcut) return;
   int z = static_cast<int>((parsAtBeamSpot[1]+m_zcut)*m_zstep);
   if (z >=0 and z < m_histsize) {
     ++numberHistogram[z];
     zWeightedHistogram[z] += parsAtBeamSpot[1];
     ptWeightedHistogram[z] += pT;
   }
   
 }

◆ filterPassed()

virtual bool AthReentrantAlgorithm::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 135 of file AthReentrantAlgorithm.h.

                                                            {
     return execState( ctx ).filterPassed();
   }

◆ filterSharedTracks()

void InDet::SiSPSeededTrackFinder::filterSharedTracks ( std::multimap< double, Trk::Track * > & scoredTracks ) const

protected

cleans up the collection of quality filtered tracks.

Candidates which share most of their hits (steered by m_freeCut) with higher quality candidates are erased from the multimap

Parameters

[in,out] scoredTracks Track candidates, sorted by by score, best scored first (implemented by assigning negative sign to scores)

loop over all track candidates, sorted by quality

loop over measurements on the track candidate

get the PRD from the measurement

increase cluster count

and check if the cluster was already used in a previous ( = higher quality) track

if not, record as a free (not prevously used) cluster

check if the track has the minimum number of free clusters or if it has no shared clusters

if this is fulfilled, we keep the candidate add the free clusters to our cluster set

if we do not keep the track, clean up candidate

Definition at line 865 of file SiSPSeededTrackFinder.cxx.

 {
   std::set<const Trk::PrepRawData*> clusters;
   
   std::vector<const Trk::PrepRawData*> freeClusters;
   freeClusters.reserve(15);    
   
   std::multimap<double, Trk::Track*>::iterator it_qualityAndTrack = qualitySortedTracks.begin();
  
   while (it_qualityAndTrack!=qualitySortedTracks.end()) {
     freeClusters.clear();
  
     std::set<const Trk::PrepRawData*>::iterator it_clustersEnd = clusters.end();
  
     int nClusters = 0; 
     for (const Trk::MeasurementBase* m: *((*it_qualityAndTrack).second->measurementsOnTrack())) {
  
       const Trk::PrepRawData* pr = (static_cast<const Trk::RIO_OnTrack*>(m))->prepRawData();
       if (pr) {
         ++nClusters;
         if (clusters.find(pr)==it_clustersEnd) {
           freeClusters.push_back(pr); 
         }
       }
     }
  
     int nFreeClusters = static_cast<int>(freeClusters.size()); 
     if (nFreeClusters >= m_nfreeCut || nFreeClusters==nClusters) {
       clusters.insert(freeClusters.begin(), freeClusters.end());
       ++it_qualityAndTrack;
     } else {
       delete (*it_qualityAndTrack).second;
       qualitySortedTracks.erase(it_qualityAndTrack++);
     }
   }
 }

◆ filterSharedTracksFast()

void InDet::SiSPSeededTrackFinder::filterSharedTracksFast ( std::multimap< double, Trk::Track * > & scoredTracks ) const

protected

loop over all track candidates, sorted by quality

loop over track states on surface of the track candidate

get the PRD from the measurement

increase cluster count

and check if the cluster was already used in a previous ( = higher quality) track

if not, record as a free (not prevously used) cluster

add the free clusters to our cluster set

if this is fulfilled, we keep the candidate

if we do not keep the track, clean up candidate

Definition at line 913 of file SiSPSeededTrackFinder.cxx.

 {
   std::set<const Trk::PrepRawData*> clusters;
  
   std::vector<const Trk::PrepRawData*> freeClusters;
   freeClusters.reserve(15);
  
   std::multimap<double, Trk::Track*>::iterator it_qualityAndTrack = qualitySortedTracks.begin();
  
   while (it_qualityAndTrack!=qualitySortedTracks.end()) {
     freeClusters.clear();
  
     std::set<const Trk::PrepRawData*>::iterator it_clustersEnd = clusters.end();
  
     int nClusters = 0;
     int nPixels = 0;
     for (const Trk::TrackStateOnSurface* tsos: *((*it_qualityAndTrack).second->trackStateOnSurfaces())) {
  
       if(!tsos->type(Trk::TrackStateOnSurface::Measurement)) continue;
       const Trk::FitQualityOnSurface fq =  tsos->fitQualityOnSurface();
       if(!fq) continue;
       if(fq.numberDoF() == 2) ++nPixels;
  
       const Trk::MeasurementBase* mb = tsos->measurementOnTrack();
       const Trk::RIO_OnTrack*     ri = dynamic_cast<const Trk::RIO_OnTrack*>(mb);
       if(!ri) continue;
       const Trk::PrepRawData* pr = ri->prepRawData();
       if (pr) {
         ++nClusters;
         if (clusters.find(pr)==it_clustersEnd) {
           freeClusters.push_back(pr);
         }
       }
     }
  
     clusters.insert(freeClusters.begin(), freeClusters.end());
  
     int nFreeClusters = static_cast<int>(freeClusters.size());
     if( passEtaDepCuts( (*it_qualityAndTrack).second, nClusters, nFreeClusters, nPixels) ){
       ++it_qualityAndTrack;
     } else {
       delete (*it_qualityAndTrack).second;
       qualitySortedTracks.erase(it_qualityAndTrack++);
     }
   }
 }

◆ finalize()

StatusCode InDet::SiSPSeededTrackFinder::finalize ( )

overridevirtual

Definition at line 712 of file SiSPSeededTrackFinder.cxx.

 {
   
     dump(MSG::INFO, &m_counterTotal);
   
   return StatusCode::SUCCESS;
 }

◆ findZvertex()

void InDet::SiSPSeededTrackFinder::findZvertex	(	std::list< Trk::Vertex > &	vertexList,
		std::pair< double, double > &	zBoundaries,
		const std::vector< int > &	numberWeightedhistogram,
		const std::vector< double > &	zWeightedHistogram,
		const std::vector< double > &	ptWeightedHistogram
	)		const

protected

estimates a set of vertex positions and a z interval for the second track finding pass using the input histograms populated using fillZHistogram.

the first two arguments serve as output, while the three histograms are input.

Parameters

[out]	vertexList	- will be populated with vertex candidates if m_useNewStrategy is set
[out]	zBoundaries	- will be populated with edges of a z interval corresponding to the estimated vertex locations
[in]	numberWeightedhistogram	vector representing a histogram in z, counting the tracks per bin
[in]	zWeightedHistogram	vector representing a histogram in z, counting the tracks per bin weighted by their z values
[in]	ptWeightedHistogram	vector representing a histogram in z, counting the tracks per bin weighted by their pt values

loop over the bins in z. Start with the second, end with the second to last, as we investigate triplets of bins

total entries in the three neighbouring bins

if we have at least 3 tracks in this z0 interval with the peak of the local triplet in the central bin, add a vertex candidate

vertex z estimate as the mean z0 of the tracks in the three bins

if this vertex is to the left of the lower edge or the right of the upper edge, broaden the z interval accordingly to include this vertex.

also get the sum(pt) of the vertex candidate

stop when we have collected the requested number of leading vertices

insert one vertex from the ntracks-sorted list and one from the pt-sorted list. If both are the same, only one insertion is performed

if we did not find any PV candidates, assume a very wide interval

otherwise, add a reasonable bit of space on each side of the first / last vertex candidate

Definition at line 1021 of file SiSPSeededTrackFinder.cxx.

 {
   zBoundaries = {1000., -1000};
  
   std::multimap<int   ,double> vertexZ_sortedByNtracks;
   std::multimap<double,double> vertexZ_sortedBySumPt;
  
   int lastBin = m_histsize-1;
   int minBinContentSum = 3;
  
   for (int binIndex=1; binIndex<lastBin; ++binIndex) { 
  
     int vertexNtracks = numberHistogram.at(binIndex-1)+numberHistogram.at(binIndex)+numberHistogram.at(binIndex+1);
  
     if (vertexNtracks>=minBinContentSum and (numberHistogram.at(binIndex) >= numberHistogram.at(binIndex-1) and numberHistogram.at(binIndex) >= numberHistogram.at(binIndex+1))) {
       double vertexZestimate = (zWeightedHistogram.at(binIndex-1)+zWeightedHistogram.at(binIndex)+zWeightedHistogram.at(binIndex+1))/static_cast<double>(vertexNtracks);
  
       if (vertexZestimate < zBoundaries.first) zBoundaries.first = vertexZestimate;
       if (vertexZestimate > zBoundaries.second) zBoundaries.second = vertexZestimate;
  
       if (m_useNewStrategy) { 
         double vertexSumPt = ptWeightedHistogram.at(binIndex-1)+ptWeightedHistogram.at(binIndex)+ptWeightedHistogram.at(binIndex+1);
         vertexZ_sortedByNtracks.insert(std::make_pair(-vertexNtracks, vertexZestimate));
         vertexZ_sortedBySumPt.insert(std::make_pair(-vertexSumPt, vertexZestimate));
       }
     }
   }
  
   if (m_useNewStrategy) {
  
     std::set<double> leadingVertices;
     int n = 0;
     std::multimap<double, double>::iterator vertex_pt_and_z = vertexZ_sortedBySumPt.begin();
     for (std::pair<int, double> nTrackAndZ: vertexZ_sortedByNtracks) {
       if (n++ >= m_nvertex) break; 
       leadingVertices.insert(nTrackAndZ.second);
       leadingVertices.insert((*vertex_pt_and_z++).second);
     }
  
     for (double v: leadingVertices) {
       vertexZList.emplace_back(Amg::Vector3D{0.,0.,v});
     }
   }
   if (zBoundaries.first > zBoundaries.second) {
     zBoundaries.first = -1000.;
     zBoundaries.second = +1000.;
   } else {
     zBoundaries.first -= 20.;
     zBoundaries.second += 20.;
   }
 } 

◆ initialize()

StatusCode InDet::SiSPSeededTrackFinder::initialize ( )

overridevirtual

Initialisation.

optional PRD to track association map

Get tool for space points seed maker

Get track-finding tool

Get RungeKutta propagator tool

Setup for magnetic field

Setup for Z-histograms

Get output print level

Definition at line 73 of file SiSPSeededTrackFinder.cxx.

 {
   ATH_CHECK(m_evtKey.initialize());
   ATH_CHECK(m_mbtsKey.initialize(m_useMBTS));
   ATH_CHECK(m_SpacePointsPixelKey.initialize(SG::AllowEmpty));
   ATH_CHECK(m_SpacePointsSCTKey.initialize(SG::AllowEmpty));
   ATH_CHECK(m_outputTracksKey.initialize());
  
   ATH_CHECK( m_prdToTrackMap.initialize( !m_prdToTrackMap.key().empty() ) );
  
   ATH_CHECK( m_caloClusterROIKey.initialize(m_useITkConvSeeded) );
  
   ATH_CHECK( m_seedsmaker.retrieve() );
   ATH_CHECK( m_zvertexmaker.retrieve( DisableTool{ not m_useZvertexTool } ));
  
   ATH_CHECK( m_trackmaker.retrieve());
  
   ATH_CHECK( m_regsel_strip.retrieve( DisableTool{ not m_useITkConvSeeded } ) );
  
   if (m_ITKGeometry and m_doFastTracking) ATH_CHECK(m_etaDependentCutsSvc.retrieve());
  
   ATH_CHECK( m_trackSummaryTool.retrieve( DisableTool{ m_trackSummaryTool.name().empty()} ));
  
   if (m_useNewStrategy and m_beamSpotKey.key().empty()) {
     m_useNewStrategy = false;
     m_useZBoundaryFinding = false;
   }
  
   ATH_CHECK(m_beamSpotKey.initialize( (m_useNewStrategy or m_useZBoundaryFinding or m_ITKGeometry or m_useITkConvSeeded) && not m_beamSpotKey.key().empty() ));
   if (m_useNewStrategy or m_useZBoundaryFinding or m_ITKGeometry or m_useITkConvSeeded) {
  
     if (not m_beamSpotKey.key().empty()) {
       ATH_CHECK( m_proptool.retrieve() );
  
       magneticFieldInit();
       
       if (m_histsize < 100) m_histsize = 100;
       m_zstep = static_cast<double>(m_histsize)/(2.*m_zcut);
     } else {
       m_proptool.disable();
       m_useNewStrategy = false;
       m_useZBoundaryFinding = false;
     }
   } else {
     m_proptool.disable();
   }
  
   if (msgLvl(MSG::DEBUG)) {
     dump(MSG::DEBUG, nullptr);
   }
   m_neventsTotal   = 0;
   m_neventsTotalV  = 0;
   m_problemsTotal  = 0;
   m_problemsTotalV = 0;
   return StatusCode::SUCCESS;
 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ isClonable()

virtual bool InDet::SiSPSeededTrackFinder::isClonable ( ) const

inlineoverridevirtual

Make this algorithm clonable.

Reimplemented from AthReentrantAlgorithm.

Definition at line 68 of file SiSPSeededTrackFinder.h.

68 { return true; };

◆ isGoodEvent()

bool InDet::SiSPSeededTrackFinder::isGoodEvent ( const EventContext & ctx ) const

protected

EventContext is used to specify which event.

Definition at line 814 of file SiSPSeededTrackFinder.cxx.

                                                                           {
  
   if ( not m_alwaysProtectAgainstBusyEvent ) { // if not enabled the protection is only applied to specific events
     // Test MBTS information from calorimeter
     //
     if (not m_useMBTS) return true;
  
     SG::ReadHandle<xAOD::EventInfo> eventInfo{m_evtKey, ctx};
     if (not eventInfo->isEventFlagBitSet(xAOD::EventInfo::Background, xAOD::EventInfo::MBTSTimeDiffHalo) ) {
       return true;
     }
   }
  
   // Test total number pixels space points
   //
   unsigned int nsp = 0;
   if (not m_SpacePointsPixelKey.empty()) {
     SG::ReadHandle<SpacePointContainer> spacePointsPixel{m_SpacePointsPixelKey, ctx};
     if (spacePointsPixel.isValid()) {
       for (const SpacePointCollection* spc: *spacePointsPixel) {
         nsp += spc->size();
       }
       if (static_cast<int>(nsp) > m_maxPIXsp) {
         ATH_MSG_WARNING("Found more than "<<m_maxPIXsp<<" pixels space points in background event. Skip track finding");
         return false;
       }
     }
   }
   // Test total number sct space points
   //
   nsp = 0;
   if (not m_SpacePointsSCTKey.empty()) {
     SG::ReadHandle<SpacePointContainer> spacePointsSCT{m_SpacePointsSCTKey, ctx};
     if (spacePointsSCT.isValid()) {
       for (const SpacePointCollection* spc: *spacePointsSCT) {
         nsp += spc->size();
       }
       if (static_cast<int>(nsp) > m_maxSCTsp) {
         ATH_MSG_WARNING("Found more than "<<m_maxSCTsp<<" sct space points in background event. Skip track finding");
         return false;
       }
     }
   }
  
   return true;
 }

◆ itkConvStrategy()

StatusCode InDet::SiSPSeededTrackFinder::itkConvStrategy ( const EventContext & ctx ) const

protected

EventContext is used to specify which event.

For HI events we can use MBTS information from calorimeter

set up the seed maker for first pass

and run seeding - starting with an empty list of vertices for the first pass

set up the track maker

prepare a collection for the quality-sorted track canddiates

Get the value of the seed maker validation ntuple writing switch

Loop through all seeds from the first pass and attempt to form track candidates

copy all the tracks into trackList

record track candidates found, using combinatorial track finding, from the given seed

Call the ntuple writing method

Remove shared tracks with worse quality

Save good tracks in track collection

Note that for run-3 the tool here is configured to not perform a hole search, regardless of the 'false' argument below

Check if we have a hole search result for this guy

If yes: Write this information into the track summary.

Definition at line 562 of file SiSPSeededTrackFinder.cxx.

 {
   ATH_MSG_DEBUG("Executing " << name() << "::itkConvStrategy");
  
   SG::WriteHandle<TrackCollection> outputTracks{m_outputTracksKey, ctx};
   ATH_CHECK(outputTracks.record(std::make_unique<TrackCollection>()));
   if (not isGoodEvent(ctx)) {
     return StatusCode::SUCCESS;
   }
  
   SiSpacePointsSeedMakerEventData seedEventData;
  
   SG::ReadHandle<ROIPhiRZContainer> calo_rois(m_caloClusterROIKey, ctx);
   std::unique_ptr<RoiDescriptor> roiComp = std::make_unique<RoiDescriptor>(true);
  
   if(calo_rois.isValid()) {
     RoiDescriptor * roi =nullptr;
     SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle{m_beamSpotKey, ctx};
     double beamZ = beamSpotHandle->beamVtx().position().z();
     roiComp->clear();
     roiComp->setComposite();
  
     const ROIPhiRZContainer &calo_rois_ref=*calo_rois;
     for (const ROIPhiRZ &calo_roi : calo_rois_ref) {
         double phi = calo_roi.phi();
         if (std::abs(phi)>=M_PI && phi!=-M_PI) continue; // skip duplicates < -pi and >pi
         double eta = calo_roi.eta();
         double z = beamZ;
         double roiPhiMin = phi - m_deltaPhi;
         double roiPhiMax = phi + m_deltaPhi;
         double roiEtaMin = eta - m_deltaEta;
         double roiEtaMax = eta + m_deltaEta;
         double roiZMin = beamZ - m_deltaZ;
         double roiZMax = beamZ + m_deltaZ;
         roi = new RoiDescriptor( eta, roiEtaMin, roiEtaMax,phi, roiPhiMin ,roiPhiMax,z,roiZMin,roiZMax);
         roiComp->push_back(roi);
     }
   }
   else {
     ATH_MSG_ERROR("Calo RoI is not valid: " << m_caloClusterROIKey.key());
     return StatusCode::FAILURE;
   }
  
   std::vector<IdentifierHash> listOfStripIds;
   std::vector<IdentifierHash> listOfPixIds;
  
   m_regsel_strip->lookup(ctx)->HashIDList( *roiComp, listOfStripIds );
   
   m_seedsmaker->newRegion(ctx, seedEventData, listOfPixIds, listOfStripIds);
   std::list<Trk::Vertex> vertexList;
   m_seedsmaker->find3Sp(ctx, seedEventData, vertexList);
  
   const bool PIX = true ;
   const bool STRIP = true ;
   InDet::ExtendedSiTrackMakerEventData_xk trackEventData(m_prdToTrackMap);
   m_trackmaker->newEvent(ctx, trackEventData, PIX, STRIP);
  
   bool ERR = false;
   Counter_t counter{};
   const InDet::SiSpacePointsSeed* seed = nullptr;
  
   std::multimap<double, Trk::Track*> qualitySortedTrackCandidates;
  
   bool doWriteNtuple = m_seedsmaker->getWriteNtupleBoolProperty();
   long EvNumber = 0.;            //Event number variable to be used for the validation ntuple
  
   if (doWriteNtuple) {
     SG::ReadHandle<xAOD::EventInfo> eventInfo(m_evtKey,ctx);
     if(!eventInfo.isValid()) {EvNumber = -1.0;} else {EvNumber = eventInfo->eventNumber();}
   }
  
   while ((seed = m_seedsmaker->next(ctx, seedEventData))) {
  
     ++counter[kNSeeds];
  
     std::list<Trk::Track*> trackList = m_trackmaker->getTracks(ctx, trackEventData, seed->spacePoints());
     for (Trk::Track* t: trackList) {
       qualitySortedTrackCandidates.insert(std::make_pair(-trackQuality(t), t));
     }
  
     if(doWriteNtuple) { m_seedsmaker->writeNtuple(seed, !trackList.empty() ? trackList.front() : nullptr, ISiSpacePointsSeedMaker::StripSeed, EvNumber) ; }
  
     if (counter[kNSeeds] >= m_maxNumberSeeds) {
       ERR = true;
       ++m_problemsTotal;
       break;
     }
   }
  
   m_trackmaker->endEvent(trackEventData);
  
   filterSharedTracks(qualitySortedTrackCandidates);
  
   for (const std::pair<const double, Trk::Track*> & qualityAndTrack: qualitySortedTrackCandidates) {
     ++counter[kNTracks];
     if (m_trackSummaryTool.isEnabled()) {
       m_trackSummaryTool->computeAndReplaceTrackSummary(*qualityAndTrack.second,
                             false /* DO NOT suppress hole search*/);
       InDet::PatternHoleSearchOutcome theOutcome;
       if (m_writeHolesFromPattern && trackEventData.combinatorialData().findPatternHoleSearchOutcome(qualityAndTrack.second,theOutcome)){
     qualityAndTrack.second->trackSummary()->update(Trk::numberOfPixelHoles, theOutcome.nPixelHoles);
     qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTHoles, theOutcome.nSCTHoles);
     qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTDoubleHoles, theOutcome.nSCTDoubleHoles);
     qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTDeadSensors, theOutcome.nSCTDeads);
     qualityAndTrack.second->trackSummary()->update(Trk::numberOfPixelDeadSensors, theOutcome.nPixelDeads);
       }
     }
     outputTracks->push_back(qualityAndTrack.second);
   }
  
   m_counterTotal[kNSeeds] += counter[kNSeeds] ;
  
   ++m_neventsTotal;
  
   if (ERR) {
     outputTracks->clear();
   } else {
     m_counterTotal[kNTracks] += counter[kNTracks];
   }
  
   // Print common event information
   //
   if (msgLvl(MSG::DEBUG)) {
     dump(MSG::DEBUG, &counter);
   }
   return StatusCode::SUCCESS;
 }

◆ itkFastTrackingStrategy()

StatusCode InDet::SiSPSeededTrackFinder::itkFastTrackingStrategy ( const EventContext & ctx ) const

protected

EventContext is used to specify which event.

set up the track maker

We run a single passe of seeding & track finding with PPP seeds

Set up the first pass (pixel seeds), and prepare to obtain a vertex Z estimate from the candidates we find

set up the seed maker for first pass

and run seeding - starting with an empty list of vertices for the first pass

prepare a collection for the quality-sorted track canddiates

Get the value of the seed maker validation ntuple writing switch

Loop through all seeds from the first pass and attempt to form track candidates

copy all the tracks into trackList

record track candidates found, using combinatorial track finding, from the given seed

Call the ntuple writing method

Remove shared tracks with worse quality

Save good tracks in track collection

Check if we have a hole search result for this guy

If yes: Write this information into the track summary.

Definition at line 446 of file SiSPSeededTrackFinder.cxx.

 {
   SG::WriteHandle<TrackCollection> outputTracks{m_outputTracksKey, ctx};
   ATH_CHECK(outputTracks.record(std::make_unique<TrackCollection>()));
  
   const bool PIX = true ;
   const bool STRIP = true ;
   InDet::ExtendedSiTrackMakerEventData_xk trackEventData(m_prdToTrackMap);
   m_trackmaker->newTrigEvent(ctx, trackEventData, PIX, STRIP);
  
   SiSpacePointsSeedMakerEventData seedEventData;
  
   m_seedsmaker->newEvent(ctx, seedEventData, 0);
   std::list<Trk::Vertex> vertexList;
   m_seedsmaker->find3Sp(ctx, seedEventData, vertexList);
  
   bool ERR = false;
   Counter_t counter{};
   const InDet::SiSpacePointsSeed* seed = nullptr;
  
   std::multimap<double, Trk::Track*> qualitySortedTrackCandidates;
  
   bool doWriteNtuple = m_seedsmaker->getWriteNtupleBoolProperty();
   long EvNumber = 0.;            //Event number variable to be used for the validation ntuple
  
   if (doWriteNtuple) {
     SG::ReadHandle<xAOD::EventInfo> eventInfo(m_evtKey,ctx);
     if(!eventInfo.isValid()) {EvNumber = -1.0;} else {EvNumber = eventInfo->eventNumber();}
   }
  
   while ((seed = m_seedsmaker->next(ctx, seedEventData))) {
  
     ++counter[kNSeeds];
  
     std::list<Trk::Track*> trackList = m_trackmaker->getTracks(ctx, trackEventData, seed->spacePoints());
     for (Trk::Track* t: trackList) {
  
       qualitySortedTrackCandidates.insert(std::make_pair(-trackQuality(t), t));
  
     }
     if(doWriteNtuple) { m_seedsmaker->writeNtuple(seed, !trackList.empty() ? trackList.front() : nullptr, ISiSpacePointsSeedMaker::PixelSeed, EvNumber) ; }
  
     if (counter[kNSeeds] >= m_maxNumberSeeds) {
       ERR = true;
       ++m_problemsTotal;
       break;
     }
   }
  
   m_trackmaker->endEvent(trackEventData);
  
   filterSharedTracksFast(qualitySortedTrackCandidates);
  
   for (const std::pair<const double, Trk::Track*> & qualityAndTrack: qualitySortedTrackCandidates) {
     ++counter[kNTracks];
  
     if (m_trackSummaryTool.isEnabled()) {
       m_trackSummaryTool->computeAndReplaceTrackSummary(*qualityAndTrack.second);
       InDet::PatternHoleSearchOutcome theOutcome;
       if (m_writeHolesFromPattern && trackEventData.combinatorialData().findPatternHoleSearchOutcome(qualityAndTrack.second,theOutcome)){
         qualityAndTrack.second->trackSummary()->update(Trk::numberOfPixelHoles, theOutcome.nPixelHoles);
         qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTHoles, theOutcome.nSCTHoles);
         qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTDoubleHoles, theOutcome.nSCTDoubleHoles);
         qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTDeadSensors, theOutcome.nSCTDeads);
         qualityAndTrack.second->trackSummary()->update(Trk::numberOfPixelDeadSensors, theOutcome.nPixelDeads);
       }
     }
  
     outputTracks->push_back(qualityAndTrack.second);
   }
  
   m_counterTotal[kNSeeds] += counter[kNSeeds] ;
  
   ++m_neventsTotal;
  
   if (ERR) {
     outputTracks->clear();
   } else {
     m_counterTotal[kNTracks] += counter[kNTracks];
   }
  
   // Print common event information
   //
   if (msgLvl(MSG::DEBUG)) {
     dump(MSG::DEBUG, &counter);
   }
   return StatusCode::SUCCESS;
 }

◆ magneticFieldInit()

void InDet::SiSPSeededTrackFinder::magneticFieldInit ( )

protected

Definition at line 1095 of file SiSPSeededTrackFinder.cxx.

 {
   // Build MagneticFieldProperties 
   //
   if(m_fieldmode == "NoField") {
     m_fieldprop = Trk::MagneticFieldProperties(Trk::NoField);
   } else {
     m_fieldprop = Trk::MagneticFieldProperties(Trk::FastField);
   }
 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ newStrategy()

StatusCode InDet::SiSPSeededTrackFinder::newStrategy ( const EventContext & ctx ) const

protected

this method performs the track finding using the new strategy

EventContext is used to specify which event

For HI events we can use MBTS information from calorimeter

Get beam information and preparation for z -histogramming

We run two passes of seeding & track finding.

First, we use strip seeds, perform finding with them and obtain an estimate of the beam spot region from the found tracks.

Then, we run PPP-seeded tracking within the resulting z interval.

Set up the first pass (strip seeds), and prepare to obtain a vertex Z estimate from the candidates we find

set up the seed maker for first pass

and run seeding - starting with an empty list of vertices for the first pass

set up the track maker

initialize empty histograms for the vertex estimate

prepare a collection for the quality-sorted track canddiates

Get the value of the seed maker validation ntuple writing switch

Loop through all seeds from the first pass and attempt to form track candidates

we only want to fill the Z histo with the first candidate for each seed.

copy all the tracks into trackList

record track candidates found, using combinatorial track finding, from the given seed

For the first (highest quality) track from each seed, populate the vertex finding histograms

Call the ntuple writing method

Now set up the second seeding pass, using pixel seeds The seed maker will internally reconfigure itself based on the "1" argument for the "iteration" argument in this call.

perform vertex Z estimation and run second seeding pass

Estimate a Z vertex interval and, if running the new strategy, also a list of the HS candidates

pass the Z boundary pair c-array-style to satisfy existing interfaces of the seeds maker family. Trigger second seed finding pass (PPP)

Again, loop over the newly found seeds and attempt to form track candidates

Remove shared tracks with worse quality

Save good tracks in track collection

Note that for run-3 the tool here is configured to not perform a hole search, regardless of the 'false' argument below

Check if we have a hole search result for this guy

If yes: Write this information into the track summary.

Definition at line 260 of file SiSPSeededTrackFinder.cxx.

 {
   SG::WriteHandle<TrackCollection> outputTracks{m_outputTracksKey, ctx};
   ATH_CHECK(outputTracks.record(std::make_unique<TrackCollection>()));
   if (not isGoodEvent(ctx)) {
     return StatusCode::SUCCESS;
   }
  
   SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle{m_beamSpotKey, ctx};
   Trk::PerigeeSurface beamPosPerigee(beamSpotHandle->beamPos());
  
   SiSpacePointsSeedMakerEventData seedEventData;
   
   m_seedsmaker->newEvent(ctx, seedEventData, 0);
   std::list<Trk::Vertex> vertexList;
   m_seedsmaker->find3Sp(ctx, seedEventData, vertexList);
  
   const bool PIX = true ;
   const bool SCT = true ;
   InDet::ExtendedSiTrackMakerEventData_xk trackEventData(m_prdToTrackMap);
   m_trackmaker->newEvent(ctx, trackEventData, PIX, SCT);
  
   std::vector<int> numberHistogram(m_histsize, 0);
   std::vector<double> zWeightedHistogram(m_histsize, 0.);
   std::vector<double> ptWeightedHistogram(m_histsize, 0.);
  
   bool ERR = false;
   Counter_t counter{};
   const InDet::SiSpacePointsSeed* seed = nullptr;
  
   std::multimap<double, Trk::Track*> qualitySortedTrackCandidates;
  
   bool doWriteNtuple = m_seedsmaker->getWriteNtupleBoolProperty();
   long EvNumber = 0.;            //Event number variable to be used for the validation ntuple 
  
   if (doWriteNtuple) {
     SG::ReadHandle<xAOD::EventInfo> eventInfo(m_evtKey,ctx);
     if(!eventInfo.isValid()) {EvNumber = -1.0;} else {EvNumber = eventInfo->eventNumber();}
   }
  
   while ((seed = m_seedsmaker->next(ctx, seedEventData))) {
  
     ++counter[kNSeeds];
     bool firstTrack{true};
  
       std::list<Trk::Track*> trackList = m_trackmaker->getTracks(ctx, trackEventData, seed->spacePoints());
       for (Trk::Track* t: trackList) {
  
         qualitySortedTrackCandidates.insert(std::make_pair(-trackQuality(t), t));
  
         if (firstTrack and not m_ITKGeometry) {
           fillZHistogram(t, beamPosPerigee, numberHistogram, zWeightedHistogram, ptWeightedHistogram);
         }
         firstTrack = false;
       }  
       if(doWriteNtuple) { m_seedsmaker->writeNtuple(seed, !trackList.empty() ? trackList.front() : nullptr, ISiSpacePointsSeedMaker::StripSeed, EvNumber) ; } 
         
  
     if (counter[kNSeeds] >= m_maxNumberSeeds) {
       ERR = true;
       ++m_problemsTotal;
       break;
     }
   }
  
   if(not m_SpacePointsPixelKey.empty()) {
     m_seedsmaker->newEvent(ctx, seedEventData, 1);
   } else {
     ATH_MSG_WARNING("SpacePointsPixelKey is empty. Skipping the second seeding pass that uses pixel seeds.");
   }
  
   std::pair<double,double> zBoundaries;
   if (not m_ITKGeometry) {
     findZvertex(vertexList, zBoundaries, numberHistogram, zWeightedHistogram, ptWeightedHistogram);
     m_seedsmaker->find3Sp(ctx, seedEventData, vertexList, &(zBoundaries.first));
   } else {
     m_seedsmaker->find3Sp(ctx, seedEventData, vertexList);
   }
  
   while ((seed = m_seedsmaker->next(ctx, seedEventData))) {
  
     ++counter[kNSeeds];
  
     std::list<Trk::Track*> trackList = m_trackmaker->getTracks(ctx, trackEventData, seed->spacePoints());
  
     for (Trk::Track* t: trackList) {
       qualitySortedTrackCandidates.insert(std::make_pair(-trackQuality(t), t));
     }
  
     if(doWriteNtuple) { m_seedsmaker->writeNtuple(seed, !trackList.empty() ? trackList.front() : nullptr, ISiSpacePointsSeedMaker::PixelSeed, EvNumber); }
  
     if (counter[kNSeeds] >= m_maxNumberSeeds) {
       ERR = true;
       ++m_problemsTotal;
       break;
     }
   }
  
   m_trackmaker->endEvent(trackEventData);
  
   filterSharedTracks(qualitySortedTrackCandidates);
  
   for (const std::pair<const double, Trk::Track*> & qualityAndTrack: qualitySortedTrackCandidates) {
     ++counter[kNTracks];
     if (m_trackSummaryTool.isEnabled()) {
        m_trackSummaryTool->computeAndReplaceTrackSummary(*qualityAndTrack.second,
                                                          false /* DO NOT suppress hole search*/);
        InDet::PatternHoleSearchOutcome theOutcome; 
        if (m_writeHolesFromPattern && trackEventData.combinatorialData().findPatternHoleSearchOutcome(qualityAndTrack.second,theOutcome)){
          qualityAndTrack.second->trackSummary()->update(Trk::numberOfPixelHoles, theOutcome.nPixelHoles); 
          qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTHoles, theOutcome.nSCTHoles); 
          qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTDoubleHoles, theOutcome.nSCTDoubleHoles); 
          qualityAndTrack.second->trackSummary()->update(Trk::numberOfSCTDeadSensors, theOutcome.nSCTDeads); 
          qualityAndTrack.second->trackSummary()->update(Trk::numberOfPixelDeadSensors, theOutcome.nPixelDeads); 
        }
     }
     outputTracks->push_back(qualityAndTrack.second);
   }
  
   m_counterTotal[kNSeeds] += counter[kNSeeds] ;
  
   ++m_neventsTotal;
  
   if (ERR) {
     outputTracks->clear();
   } else {
     m_counterTotal[kNTracks] += counter[kNTracks];
   }
  
   // Print common event information
   //
   if (msgLvl(MSG::DEBUG)) {
     dump(MSG::DEBUG, &counter);
   }
   return StatusCode::SUCCESS;
 }

◆ oldStrategy()

StatusCode InDet::SiSPSeededTrackFinder::oldStrategy ( const EventContext & ctx ) const

protected

this method performs the track finding using the old strategy

EventContext is used to specify which event

For HI events we can use MBTS information from calorimeter

Definition at line 175 of file SiSPSeededTrackFinder.cxx.

 {
   SG::WriteHandle<TrackCollection> outputTracks{m_outputTracksKey, ctx};
   ATH_CHECK(outputTracks.record(std::make_unique<TrackCollection>()));
   if (not isGoodEvent(ctx)) {
     return StatusCode::SUCCESS;
   }
  
   SiSpacePointsSeedMakerEventData seedEventData;
   bool ZVE = false;
   if (m_useZvertexTool) {
     std::list<Trk::Vertex> vertices = m_zvertexmaker->newEvent(ctx, seedEventData);
     if (not vertices.empty()) ZVE = true;
     m_seedsmaker->find3Sp(ctx, seedEventData, vertices);
   } else {
     m_seedsmaker->newEvent(ctx, seedEventData, -1);
     std::list<Trk::Vertex> vertexList;
     m_seedsmaker->find3Sp(ctx, seedEventData, vertexList);
   }
  
   const bool PIX = true;
   const bool SCT = true;
   InDet::ExtendedSiTrackMakerEventData_xk trackEventData(m_prdToTrackMap);
   m_trackmaker->newEvent(ctx, trackEventData, PIX, SCT);
  
   bool ERR = false;
   Counter_t counter{};
   const InDet::SiSpacePointsSeed* seed = nullptr;
   std::multimap<double, Trk::Track*> qualitySortedTrackCandidates;
   // Loop through all seed and reconsrtucted tracks collection preparation
   //
   while ((seed = m_seedsmaker->next(ctx, seedEventData))) {
     ++counter[kNSeeds];
     std::list<Trk::Track*> trackList = m_trackmaker->getTracks(ctx, trackEventData, seed->spacePoints());
     for (Trk::Track* t: trackList) {
       qualitySortedTrackCandidates.insert(std::make_pair(-trackQuality(t), t));
     }
     if (not ZVE and (counter[kNSeeds] >= m_maxNumberSeeds)) {
       ERR = true;
       ++m_problemsTotal;
       break;
     }
   }
   m_trackmaker->endEvent(trackEventData);
  
   // Remove shared tracks with worse quality
   //
   filterSharedTracks(qualitySortedTrackCandidates);
  
   // Save good tracks in track collection
   //
   for (const std::pair<const double, Trk::Track*> & qualityAndTrack: qualitySortedTrackCandidates) {
     ++counter[kNTracks];
     if (m_trackSummaryTool.isEnabled()) {
        m_trackSummaryTool->computeAndReplaceTrackSummary(*(qualityAndTrack.second),
                                                          false /* DO NOT suppress hole search*/);
     }
     outputTracks->push_back(qualityAndTrack.second);
   }
  
   m_counterTotal[kNSeeds] += counter[kNSeeds];
  
   if (ZVE) ++m_neventsTotalV;
   else     ++m_neventsTotal;
  
   if (ERR) {
     outputTracks->clear();
   } else {
     m_counterTotal[kNTracks] += counter[kNTracks];
   }
  
   // Print common event information
   //
   if (msgLvl(MSG::DEBUG)) {
     dump(MSG::DEBUG, &counter);
   }
  
   return StatusCode::SUCCESS;
 }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ passEtaDepCuts()

bool InDet::SiSPSeededTrackFinder::passEtaDepCuts	(	const Trk::Track *	track,
		int	nClusters,
		int	nFreeClusters,
		int	nPixels
	)		const

protected

Definition at line 1111 of file SiSPSeededTrackFinder.cxx.

 {
   Trk::TrackStates::const_iterator  m = track->trackStateOnSurfaces()->begin();
   const Trk::TrackParameters* par = (*m)->trackParameters();
   if(!par) return false;
  
   double eta = std::abs(par->eta());
   if(nClusters               < m_etaDependentCutsSvc->getMinSiHitsAtEta(eta)) return false;
   if(nFreeClusters           < m_etaDependentCutsSvc->getMinSiNotSharedAtEta(eta)) return false;
   if(nClusters-nFreeClusters > m_etaDependentCutsSvc->getMaxSharedAtEta(eta)) return false;
   if(nPixels                 < m_etaDependentCutsSvc->getMinPixelHitsAtEta(eta)) return false;
  
   if(par->pT() < m_etaDependentCutsSvc->getMinPtAtEta(eta)) return false;
   if(!(*m)->type(Trk::TrackStateOnSurface::Perigee)) return true ;
   if(std::abs(par->localPosition()[0]) > m_etaDependentCutsSvc->getMaxPrimaryImpactAtEta(eta)) return false;
   return true;
 }

◆ 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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ setFilterPassed()

virtual void AthReentrantAlgorithm::setFilterPassed	(	bool	state,
		const EventContext &	ctx
	)		const

inlinevirtualinherited

Definition at line 139 of file AthReentrantAlgorithm.h.

                                                                             {
     execState( ctx ).setFilterPassed( state );
   }

◆ sysExecute()

StatusCode AthReentrantAlgorithm::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 67 of file AthReentrantAlgorithm.cxx.

 {
   return Gaudi::Algorithm::sysExecute (ctx);
 }

◆ sysInitialize()

StatusCode AthReentrantAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in InputMakerBase, and HypoBase.

Definition at line 96 of file AthReentrantAlgorithm.cxx.

                                                 {
   StatusCode sc=AthCommonDataStore<AthCommonMsg<Gaudi::Algorithm>>::sysInitialize();
  
   if (sc.isFailure()) {
     return sc;
   }
   
   ServiceHandle<ICondSvc> cs("CondSvc",name());
   for (auto h : outputHandles()) {
     if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
       // do this inside the loop so we don't create the CondSvc until needed
       if ( cs.retrieve().isFailure() ) {
         ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
         return StatusCode::SUCCESS;
       }      
       if (cs->regHandle(this,*h).isFailure()) {
         sc = StatusCode::FAILURE;
         ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                       << " with CondSvc");
       }
     }
   }
   return sc;  
 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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

◆ ATLAS_THREAD_SAFE

Counter_t m_counterTotal InDet::SiSPSeededTrackFinder::ATLAS_THREAD_SAFE {}

mutableprotected

Definition at line 162 of file SiSPSeededTrackFinder.h.

◆ m_alwaysProtectAgainstBusyEvent

BooleanProperty InDet::SiSPSeededTrackFinder::m_alwaysProtectAgainstBusyEvent {this, "alwaysProtectAgainstBusyEvents", false, "Flag that enables checking always maxNumberPIXsp and maxNumberSCTsp"}

protected

Definition at line 104 of file SiSPSeededTrackFinder.h.

◆ m_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> InDet::SiSPSeededTrackFinder::m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"}

protected

Definition at line 132 of file SiSPSeededTrackFinder.h.

◆ m_caloClusterROIKey

SG::ReadHandleKey<ROIPhiRZContainer> InDet::SiSPSeededTrackFinder::m_caloClusterROIKey {this, "EMROIPhiRZContainer", "InDetCaloClusterROIPhiRZ15GeVUnordered", "Name of the calo cluster ROIs in Phi,R,Z parameterization"}

protected

Definition at line 130 of file SiSPSeededTrackFinder.h.

◆ m_deltaEta

DoubleProperty InDet::SiSPSeededTrackFinder::m_deltaEta {this, "dEtaCaloRoI", .1}

protected

Definition at line 111 of file SiSPSeededTrackFinder.h.

◆ m_deltaPhi

DoubleProperty InDet::SiSPSeededTrackFinder::m_deltaPhi {this, "dPhiCaloRoI", .25}

protected

Definition at line 112 of file SiSPSeededTrackFinder.h.

◆ m_deltaZ

DoubleProperty InDet::SiSPSeededTrackFinder::m_deltaZ {this, "dZCaloRoI", 300.}

protected

Definition at line 113 of file SiSPSeededTrackFinder.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_doFastTracking

BooleanProperty InDet::SiSPSeededTrackFinder::m_doFastTracking {this, "doFastTracking", false, "ITk fast tracking reco"}

protected

Definition at line 102 of file SiSPSeededTrackFinder.h.

◆ m_etaDependentCutsSvc

ServiceHandle<IInDetEtaDependentCutsSvc> InDet::SiSPSeededTrackFinder::m_etaDependentCutsSvc {this, "InDetEtaDependentCutsSvc", "InDetEtaDependentCutsSvc"}

protected

service to get cut values depending on different variable

Definition at line 157 of file SiSPSeededTrackFinder.h.

◆ m_evtKey

SG::ReadHandleKey<xAOD::EventInfo> InDet::SiSPSeededTrackFinder::m_evtKey {this, "EventInfoKey", "EventInfo"}

protected

Definition at line 126 of file SiSPSeededTrackFinder.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthReentrantAlgorithm::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 153 of file AthReentrantAlgorithm.h.

◆ m_fieldmode

StringProperty InDet::SiSPSeededTrackFinder::m_fieldmode {this, "MagneticFieldMode", "MapSolenoid"}

protected

Definition at line 114 of file SiSPSeededTrackFinder.h.

◆ m_fieldprop

Trk::MagneticFieldProperties InDet::SiSPSeededTrackFinder::m_fieldprop

protected

Definition at line 151 of file SiSPSeededTrackFinder.h.

◆ m_histsize

IntegerProperty InDet::SiSPSeededTrackFinder::m_histsize {this, "HistSize", 1400}

protected

Definition at line 108 of file SiSPSeededTrackFinder.h.

◆ m_imcut

double InDet::SiSPSeededTrackFinder::m_imcut {2.}

protected

Definition at line 120 of file SiSPSeededTrackFinder.h.

◆ m_ITKGeometry

BooleanProperty InDet::SiSPSeededTrackFinder::m_ITKGeometry {this, "ITKGeometry", false, "Flag to specify if it is ITK geometry"}

protected

Definition at line 100 of file SiSPSeededTrackFinder.h.

◆ m_maxNumberSeeds

IntegerProperty InDet::SiSPSeededTrackFinder::m_maxNumberSeeds {this, "maxNumberSeeds", 3000000, "Max. number used seeds"}

protected

Definition at line 103 of file SiSPSeededTrackFinder.h.

◆ m_maxPIXsp

IntegerProperty InDet::SiSPSeededTrackFinder::m_maxPIXsp {this, "maxNumberPIXsp", 150000, "Max. number pixels space points"}

protected

Definition at line 105 of file SiSPSeededTrackFinder.h.

◆ m_maxSCTsp

IntegerProperty InDet::SiSPSeededTrackFinder::m_maxSCTsp {this, "maxNumberSCTsp", 500000, "Max. number sct space points"}

protected

Definition at line 106 of file SiSPSeededTrackFinder.h.

◆ m_mbtsKey

SG::ReadDecorHandleKey<xAOD::EventInfo> InDet::SiSPSeededTrackFinder::m_mbtsKey {this, "MBTSKey", "EventInfo.MBTSBackground"}

protected

Definition at line 134 of file SiSPSeededTrackFinder.h.

◆ m_neventsTotal

std::atomic_int InDet::SiSPSeededTrackFinder::m_neventsTotal {0}

mutableprotected

Number events.

Definition at line 164 of file SiSPSeededTrackFinder.h.

◆ m_neventsTotalV

std::atomic_int InDet::SiSPSeededTrackFinder::m_neventsTotalV {0}

mutableprotected

Number events.

Definition at line 165 of file SiSPSeededTrackFinder.h.

◆ m_nfreeCut

IntegerProperty InDet::SiSPSeededTrackFinder::m_nfreeCut {this, "FreeClustersCut", 1, "Min number free clusters"}

protected

Definition at line 107 of file SiSPSeededTrackFinder.h.

◆ m_nvertex

IntegerProperty InDet::SiSPSeededTrackFinder::m_nvertex {this, "maxVertices", 4}

protected

Definition at line 109 of file SiSPSeededTrackFinder.h.

◆ m_outputTracksKey

SG::WriteHandleKey<TrackCollection> InDet::SiSPSeededTrackFinder::m_outputTracksKey {this, "TracksLocation", "SiSPSeededTracks"}

protected

Definition at line 136 of file SiSPSeededTrackFinder.h.

◆ m_prdToTrackMap

SG::ReadHandleKey<Trk::PRDtoTrackMap> InDet::SiSPSeededTrackFinder::m_prdToTrackMap {this,"PRDtoTrackMap",""}

protected

Definition at line 129 of file SiSPSeededTrackFinder.h.

◆ m_problemsTotal

std::atomic_int InDet::SiSPSeededTrackFinder::m_problemsTotal {0}

mutableprotected

Number events with number seeds > maxNumber.

Definition at line 166 of file SiSPSeededTrackFinder.h.

◆ m_problemsTotalV

std::atomic_int InDet::SiSPSeededTrackFinder::m_problemsTotalV {0}

mutableprotected

Number events with number seeds > maxNumber.

Definition at line 167 of file SiSPSeededTrackFinder.h.

◆ m_proptool

PublicToolHandle<Trk::IPatternParametersPropagator> InDet::SiSPSeededTrackFinder::m_proptool {this, "PropagatorTool", "Trk::RungeKuttaPropagator/InDetPropagator"}

protected

Definition at line 144 of file SiSPSeededTrackFinder.h.

◆ m_pTcut

double InDet::SiSPSeededTrackFinder::m_pTcut {500}

protected

Definition at line 119 of file SiSPSeededTrackFinder.h.

◆ m_regsel_strip

ToolHandle<IRegSelTool> InDet::SiSPSeededTrackFinder::m_regsel_strip { this, "RegSelTool_Strip", "RegSelTool/RegSelTool_SCT" }

protected

Definition at line 146 of file SiSPSeededTrackFinder.h.

◆ m_seedsmaker

ToolHandle<ISiSpacePointsSeedMaker> InDet::SiSPSeededTrackFinder::m_seedsmaker {this, "SeedsTool", "InDet::SiSpacePointsSeedMaker_ATLxk/InDetSpSeedsMaker", "Space poins seed maker"}

protected

Definition at line 141 of file SiSPSeededTrackFinder.h.

◆ m_SpacePointsPixelKey

SG::ReadHandleKey<SpacePointContainer> InDet::SiSPSeededTrackFinder::m_SpacePointsPixelKey {this, "SpacePointsPixelName", "PixelSpacePoints"}

protected

Definition at line 128 of file SiSPSeededTrackFinder.h.

◆ m_SpacePointsSCTKey

SG::ReadHandleKey<SpacePointContainer> InDet::SiSPSeededTrackFinder::m_SpacePointsSCTKey {this, "SpacePointsSCTName", "SCT_SpacePoints"}

protected

Definition at line 127 of file SiSPSeededTrackFinder.h.

◆ m_trackmaker

ToolHandle<ISiTrackMaker> InDet::SiSPSeededTrackFinder::m_trackmaker {this, "TrackTool", "InDet::SiTrackMaker_xk/InDetSiTrackMaker", "Track maker"}

protected

Definition at line 143 of file SiSPSeededTrackFinder.h.

◆ m_trackSummaryTool

ToolHandle<Trk::IExtendedTrackSummaryTool> InDet::SiSPSeededTrackFinder::m_trackSummaryTool {this, "TrackSummaryTool", "InDetTrackSummaryToolNoHoleSearch"}

protected

Definition at line 145 of file SiSPSeededTrackFinder.h.

◆ m_useITkConvSeeded

BooleanProperty InDet::SiSPSeededTrackFinder::m_useITkConvSeeded {this, "useITkConvSeeded", false, "ITk EM-seeded conversion reco"}

protected

Definition at line 101 of file SiSPSeededTrackFinder.h.

◆ m_useMBTS

BooleanProperty InDet::SiSPSeededTrackFinder::m_useMBTS {this, "useMBTSTimeDiff", false}

protected

Definition at line 97 of file SiSPSeededTrackFinder.h.

◆ m_useNewStrategy

BooleanProperty InDet::SiSPSeededTrackFinder::m_useNewStrategy {this, "useNewStrategy", false}

protected

Definition at line 98 of file SiSPSeededTrackFinder.h.

◆ m_useZBoundaryFinding

BooleanProperty InDet::SiSPSeededTrackFinder::m_useZBoundaryFinding {this, "useZBoundFinding", false}

protected

Definition at line 99 of file SiSPSeededTrackFinder.h.

◆ m_useZvertexTool

BooleanProperty InDet::SiSPSeededTrackFinder::m_useZvertexTool {this, "useZvertexTool", true}

protected

Definition at line 95 of file SiSPSeededTrackFinder.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_writeHolesFromPattern

BooleanProperty InDet::SiSPSeededTrackFinder::m_writeHolesFromPattern {this, "writeHolesFromPattern", false,"Flag to activate writing hole info from the pattern recognition"}

protected

Definition at line 96 of file SiSPSeededTrackFinder.h.

◆ m_zcut

DoubleProperty InDet::SiSPSeededTrackFinder::m_zcut {this, "Zcut", 350.}

protected

Definition at line 110 of file SiSPSeededTrackFinder.h.

◆ m_zstep

double InDet::SiSPSeededTrackFinder::m_zstep {0.}

protected

Definition at line 121 of file SiSPSeededTrackFinder.h.

◆ m_zvertexmaker

ToolHandle<ISiZvertexMaker> InDet::SiSPSeededTrackFinder::m_zvertexmaker {this, "ZvertexTool", "InDet::SiZvertexMaker_xk/InDetSiZvertexMaker", "Space poins z-vertex maker"}

protected

Definition at line 142 of file SiSPSeededTrackFinder.h.

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

Classes

Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Standard Algotithm methods

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Member Enumeration Documentation

◆ ECounter

Constructor & Destructor Documentation

◆ SiSPSeededTrackFinder()

◆ ~SiSPSeededTrackFinder()

Member Function Documentation

◆ cardinality()

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

◆ dump()

◆ dumpevent()

◆ dumptools()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ fillZHistogram()

◆ filterPassed()

◆ filterSharedTracks()

◆ filterSharedTracksFast()

◆ finalize()

◆ findZvertex()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ isGoodEvent()

◆ itkConvStrategy()

◆ itkFastTrackingStrategy()

◆ magneticFieldInit()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ newStrategy()

◆ oldStrategy()

◆ outputHandles()

◆ passEtaDepCuts()

◆ renounce()

◆ renounceArray()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_alwaysProtectAgainstBusyEvent

◆ m_beamSpotKey

◆ m_caloClusterROIKey

◆ m_deltaEta

◆ m_deltaPhi

◆ m_deltaZ

◆ m_detStore

◆ m_doFastTracking

◆ m_etaDependentCutsSvc

◆ m_evtKey

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_fieldmode

◆ m_fieldprop