#include <TRT_SeededTrackFinder.h>

Inheritance diagram for InDet::TRT_SeededTrackFinder:

Collaboration diagram for InDet::TRT_SeededTrackFinder:

Classes
struct	Stat_t
	Global Counters for final algorithm statistics. More...

Public Member Functions
	TRT_SeededTrackFinder (const std::string &name, ISvcLocator *pSvcLocator)
	Standard Algorithm methods More...

virtual	~TRT_SeededTrackFinder ()

StatusCode	initialize () override

StatusCode	execute (const EventContext &ctx) const override

StatusCode	finalize () override

MsgStream &	dump (MsgStream &out) const
	Print internal tool parameters and status More...

std::ostream &	dump (std::ostream &out) const

virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...

virtual bool	isClonable () const override
	Specify if the algorithm is clonable. 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
Trk::Track *	mergeSegments (const Trk::Track &, const Trk::TrackSegment &) const
	Protected methods More...

Trk::Track *	mergeExtension (const Trk::Track &, std::vector< const Trk::MeasurementBase * > &) const
	Merge a TRT track extension and a Si track component into one global ID track. More...

Trk::Track *	segToTrack (const EventContext &, const Trk::TrackSegment &) const
	Transform a TRT track segment into a track More...

void	Analyze (TrackCollection *) const
	Do some statistics analysis at the end of each event. More...

MsgStream &	dumptools (MsgStream &out) const

MsgStream &	dumpevent (MsgStream &out, const InDet::TRT_SeededTrackFinder::Stat_t &stat) const

void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

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

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Protected Attributes
bool	m_doRefit

bool	m_doExtension
	Do final careful refit of tracks. More...

bool	m_rejectShortExten
	Find the TRT extension of the Si track segment. More...

bool	m_doStat
	use extension only if better than original track More...

bool	m_saveTRT
	Statistics of final tracks. More...

int	m_MaxSegNum
	Output stand-alone TRT segments. More...

unsigned int	m_minTRTonSegment
	Maximum number of segments to be handled. More...

unsigned int	m_minTRTonly
	Minimum number of TRT hits on segment. More...

ToolHandle< ITRT_SeededTrackFinder >	m_trackmaker
	Minimum number of TRT hits on TRT only. More...

ToolHandle< Trk::ITrackFitter >	m_fitterTool
	Track maker tool. More...

ToolHandle< ITRT_TrackExtensionTool >	m_trtExtension
	Refitting tool. More...

SG::ReadHandleKey< Trk::SegmentCollection >	m_SegmentsKey

SG::WriteHandleKey< TrackCollection >	m_outTracksKey
	TRT segments to use. More...

SG::ReadHandleKey< Trk::PRDtoTrackMap >	m_prdToTrackMap

ToolHandle< Trk::IExtendedTrackSummaryTool >	m_trackSummaryTool

ToolHandle< Trk::IExtrapolator >	m_extrapolator
	the extrapoator More...

SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey

bool	m_SiExtensionCuts
	enable cuts after Si segment finding More...

double	m_minPt
	minimal Pt cut More...

double	m_maxEta
	maximal Eta cut More...

double	m_maxRPhiImp
	maximal RPhi impact parameter cut More...

double	m_maxZImp
	maximal z impact parameter cut More...

bool	m_caloSeededRoI

SG::ReadHandleKey< ROIPhiRZContainer >	m_caloClusterROIKey

ToolHandle< IRegSelTool >	m_regionSelector

float	m_deltaEta
	delta Eta used for RoI creation More...

float	m_deltaPhi
	delta Phi used for RoI creation More...

float	m_deltaZ
	delta Z used for RoI creation More...

std::mutex m_statMutex	ATLAS_THREAD_SAFE

Stat_t m_totalStat	ATLAS_THREAD_SAFE

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

Detailed Description

InDet::TRT_SeededTrackFinde is an algorithm which produces tracks moving outside-in in the Inner Detector.

Definition at line 54 of file TRT_SeededTrackFinder.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TRT_SeededTrackFinder()

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

Standard Algorithm methods

Definition at line 38 of file TRT_SeededTrackFinder.cxx.

   : AthReentrantAlgorithm(name, pSvcLocator),
     m_trackmaker("InDet::TRT_SeededTrackFinderTool"),
     m_fitterTool("Trk::GlobalChi2Fitter/InDetTrackFitter"),
     m_SegmentsKey("TRTSegments"),
     m_outTracksKey("TRTSeededTracks")
 {
   m_doRefit          = false                                ;       //Do a final careful refit of tracks
   m_doExtension      = false                                ;       //Find the track TRT extension
   m_rejectShortExten = false                                ;
   m_doStat           = false                                ;       //Statistics of final tracks
   m_saveTRT          = true                                 ;       //Save stand-alone TRT segments
   m_MaxSegNum        = 5000                                 ;
  
   declareProperty("TrackTool"                  ,m_trackmaker       ); //Back tracking tool
   declareProperty("RefitterTool"               ,m_fitterTool       ); //Track refit tool
   declareProperty("InputSegmentsLocation"      ,m_SegmentsKey      ); //Input track collection
   declareProperty("OutputTracksLocation"       ,m_outTracksKey     ); //Output track collection
   declareProperty("FinalRefit"                 ,m_doRefit          ); //Do a final careful refit of tracks
   declareProperty("TrtExtension"               ,m_doExtension      ); //Find the TRT extension of the track
   declareProperty("RejectShortExtension"       ,m_rejectShortExten ); //Reject short extensions
   declareProperty("FinalStatistics"            ,m_doStat           ); //Statistics of final tracks
   declareProperty("OutputSegments"             ,m_saveTRT          ); //Save stand-alone TRT segments
   declareProperty("MaxNumberSegments"          ,m_MaxSegNum        ); //Maximum Number of segments to be handled
   declareProperty("MinTRTonSegment"            ,m_minTRTonSegment = 10); //Minimum Number of Hits on segment
   declareProperty("MinTRTonly"                 ,m_minTRTonly      = 15); //Minimum Number of Hits on TRT only
   // --- selection cuts after SI extension
   declareProperty("SiExtensionCuts",   m_SiExtensionCuts    = false);
   declareProperty("minPt"             ,m_minPt              = 500. );
   declareProperty("maxEta"            ,m_maxEta             = 2.7  );
   declareProperty("maxRPhiImp"        ,m_maxRPhiImp         = 10.  );
   declareProperty("maxZImp"           ,m_maxZImp            = 250. );
   declareProperty("Extrapolator"      ,m_extrapolator              );
   declareProperty("CaloSeededRoI"     ,m_caloSeededRoI=false       );
   declareProperty("dEtaCaloRoI"       ,m_deltaEta=0.1              );
   declareProperty("dPhiCaloRoI"       ,m_deltaPhi=0.25             );
   declareProperty("dZCaloRoI"         ,m_deltaZ=300                );
 }

◆ ~TRT_SeededTrackFinder()

virtual InDet::TRT_SeededTrackFinder::~TRT_SeededTrackFinder ( )

inlinevirtual

Definition at line 67 of file TRT_SeededTrackFinder.h.

67 {}

Member Function Documentation

◆ Analyze()

void InDet::TRT_SeededTrackFinder::Analyze ( TrackCollection * ) const

protected

Do some statistics analysis at the end of each event.

Loop over tracks in track collection

Concentrate on the Si component of the track

Definition at line 701 of file TRT_SeededTrackFinder.cxx.

                                                             {
   if(msgLvl(MSG::DEBUG)) {
     ATH_MSG_DEBUG( "Analyzing tracks..." );
     ATH_MSG_DEBUG( "Number of back tracks " << (tC->size()) );
     int nsct1{}, nsct2{}, nsct3{}, nsct4{}; //SCT layer counters
     int nsctTot1{}, nsctTot2{}, nsctTot3{}, nsctTot4{}; //SCT layer counters
     int npix1{}, npix2{}, npix3{}; //Pixel layer counters
     int npixTot1{}, npixTot2{}, npixTot3{}; //Pixel layer counters
     TrackCollection::const_iterator r    = tC->begin();
     TrackCollection::const_iterator re = tC->end();
     for (; r != re ; ++r){
       nsct1=nsct2=nsct3=nsct4=0; npix1=npix2=npix3=0; 
       const Trk::TrackStates* newtsos = (*r)->trackStateOnSurfaces();
       if(!newtsos) continue;
       Trk::TrackStates::const_iterator itp, itpe=newtsos->end();
       for(itp=newtsos->begin(); itp!=itpe; ++itp){
         const InDet::SiClusterOnTrack* clus = dynamic_cast<const InDet::SiClusterOnTrack*>((*itp)->measurementOnTrack());
         if(clus && ((*itp)->type(Trk::TrackStateOnSurface::Measurement))){  //Count the number of hits used in the track
           double rc = clus->globalPosition().perp();
           if((40.<=rc)&&(rc<80.)){npix1++;}    //1st pixel layer
           if((80.<=rc)&&(rc<100.)){npix2++;}   //2nd pixel layer
           if((100.<=rc)&&(rc<150.)){npix3++;}  //3rd pixel layer
           if((280.<=rc)&&(rc<340.)){nsct1++;}  //1st SCT layer
           if((340.<=rc)&&(rc<390.)){nsct2++;}  //2nd SCT layer
           if((390.<=rc)&&(rc<460.)){nsct3++;}  //3rd SCT layer
           if((460.<=rc)&&(rc<550.)){nsct4++;}  //4th SCT layer
         }
       }
       nsctTot1+=nsct1; nsctTot2+=nsct2; nsctTot3+=nsct3; nsctTot4+=nsct4;
       npixTot1+=npix1; npixTot2+=npix2; npixTot3+=npix3;
     }
     ATH_MSG_DEBUG("Total hits on 1st SCT: "<<nsctTot1<<" 2nd SCT: "<<nsctTot2<<" 3rd SCT: "<<nsctTot3<<" 4th SCT: "<<nsctTot4);
     ATH_MSG_DEBUG("Total hits on 1st Pixel: "<<npixTot1<<" 2nd Pixel: "<<npixTot2<<" 3rd Pixel: "<<npixTot3);
   }
 }

◆ 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() [1/2]

MsgStream& InDet::TRT_SeededTrackFinder::dump ( MsgStream & out ) const

Print internal tool parameters and status

◆ dump() [2/2]

std::ostream& InDet::TRT_SeededTrackFinder::dump ( std::ostream & out ) const

◆ dumpevent()

MsgStream & InDet::TRT_SeededTrackFinder::dumpevent	(	MsgStream &	out,
		const InDet::TRT_SeededTrackFinder::Stat_t &	stat
	)		const

protected

Definition at line 494 of file TRT_SeededTrackFinder.cxx.

 {
   out<<"|-------------------------------------------------------------------"<<std::endl;
   out<<"|  Investigated :"<<std::endl
      <<"| "<<std::setw(7)<<stat.m_counter[Stat_t::kNTrtSeg]       <<" TRT Segments on input"<<std::endl
      <<"| "<<std::setw(7)<<stat.m_counter[Stat_t::kNTrtFailSel]   <<" TRT Segments fail selection on input"<<std::endl
      <<"| "<<std::setw(7)<<stat.m_counter[Stat_t::Stat_t::kNTrtSegGood]   <<" TRT Segments after selection"<<std::endl;
   if (stat.m_counter[Stat_t::kNTrtLimit]>0) {
     out<<"|-------------------------------------------------------------------"<<std::endl;
     out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::kNTrtLimit]   <<" TRT segments lost because of processing limit"<<std::endl;
   }
   out<<"|-------------------------------------------------------------------"<<std::endl;
   out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::kNTrtNoSiExt]   <<" TRT segments without Si extension"<<std::endl;
   if (m_SiExtensionCuts) {
     out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::kNExtCut]     <<" number of Si extension failing cuts"<<std::endl;
   }
   if (stat.m_counter[Stat_t::Stat_t::kNBckTrkTrt]>0) {
     out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::Stat_t::kNBckTrkTrt]  <<" number ot TRT only tracks created"<<std::endl;
   }
   if (m_doExtension) {
     out<<"|-------------------------------------------------------------------"<<std::endl;
     out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::Stat_t::kNTrtExtCalls]<<" number of times TRT extension is called"<<std::endl
        <<"| "<<std::setw(7)<<stat.m_counter[Stat_t::kNTrtExt]     <<" number of good TRT extension"<<std::endl;
     if (m_rejectShortExten) {
       out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::Stat_t::kNTrtExtBad]<<" number of bad TRT extension"<<std::endl;
     }
     out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::Stat_t::kNTrtExtFail]<<" number of failed TRT extension and fallback"<<std::endl;
   }
   out<<"|-------------------------------------------------------------------"<<std::endl;
   out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::Stat_t::kNBckTrkSi]    <<" TRT+Si tracks created of output"<<std::endl;
   if (stat.m_counter[Stat_t::kNBckTrkSi] != stat.m_counter[Stat_t::kNBckTrk]) {
     out<<"| "<<std::setw(7)<<stat.m_counter[Stat_t::kNBckTrk]    <<" total number of tracks on output"<<std::endl;
   }
   out<<"|-------------------------------------------------------------------";
   return out;
 }

◆ dumptools()

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

protected

Definition at line 463 of file TRT_SeededTrackFinder.cxx.

 {
   int n = 65-m_trackmaker.type().size();
   std::string s1; for(int i=0; i<n; ++i) s1.append(" "); s1.append("|");
   n     = 65-m_fitterTool.type().size();
   std::string s2; for(int i=0; i<n; ++i) s2.append(" "); s2.append("|");
   n     = 65-m_trtExtension.type().size();
   std::string s3; for(int i=0; i<n; ++i) s3.append(" "); s3.append("|");
   n     = 65-m_SegmentsKey.key().size();
   std::string s4; for(int i=0; i<n; ++i) s4.append(" "); s4.append("|");
   n     = 65-m_outTracksKey.key().size();
   std::string s5; for(int i=0; i<n; ++i) s5.append(" "); s5.append("|");
  
   out<<"|----------------------------------------------------------------------"
      <<"-------------------|"
        <<std::endl;
   out<<"| Tool for TRT seeded track finding | "<<m_trackmaker.type()    <<s1<<std::endl;
   out<<"| Tool for final track refitting    | "<<m_fitterTool.type()    <<s2<<std::endl;
   out<<"| Tool for TRT trac extension       | "<<m_trtExtension.type()  <<s3<<std::endl;
   out<<"| Location of input tracks          | "<<m_SegmentsKey.key()      <<s4<<std::endl;
   out<<"| Location of output tracks         | "<<m_outTracksKey.key()     <<s5<<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::TRT_SeededTrackFinder::execute ( const EventContext & ctx ) const

override

Transform the original TRT segment into a track

Definition at line 153 of file TRT_SeededTrackFinder.cxx.

                                                                 {
   //Counters. See the include file for definitions
   Stat_t ev_stat;
   // counter
   int nTrtSegCur = 0;
   SG::ReadHandle<Trk::SegmentCollection>  segments(m_SegmentsKey,ctx);
   if(!segments.isValid()){
     ATH_MSG_FATAL ("No segment with name " << segments.name() << " found in StoreGate!");
     return StatusCode::FAILURE;
   } else {
     ATH_MSG_DEBUG ("Found segments collection " << segments.name() << " in StoreGate!");
   }
   // number of segments + statistics
   ev_stat.m_counter[Stat_t::kNTrtSeg] = int(segments->size());
   ATH_MSG_DEBUG ("TRT track container size " << ev_stat.m_counter[Stat_t::kNTrtSeg]);
   if(ev_stat.m_counter[Stat_t::kNTrtSeg]>m_MaxSegNum) {
     ATH_MSG_DEBUG ("TRT track container size huge; will process event partially if number of max segments reached !!!");
   }
   // Event dependent data of SiCombinatorialTrackFinder_xk
   InDet::ExtendedSiCombinatorialTrackFinderData_xk combinatorialData(m_prdToTrackMap);
   std::unique_ptr<InDet::ITRT_SeededTrackFinder::IEventData> event_data_p;
   if(m_caloSeededRoI ) {
     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();
       for (const ROIPhiRZ &the_roi : *calo_rois) {
           double eta = the_roi.eta();
           double phi = the_roi.phi();
           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, beamZ, roiZMin,roiZMax);
           roiComp->push_back(roi);
       }
  
     }
     std::vector<IdentifierHash> listOfSCTIds;
     std::vector<IdentifierHash> listOfPixIds;
     m_regionSelector->HashIDList(*roiComp, listOfSCTIds );
     event_data_p = m_trackmaker->newRegion(ctx, combinatorialData, listOfPixIds, listOfSCTIds);
   } else {
     event_data_p = m_trackmaker->newEvent(ctx, combinatorialData);
   }
  
   std::unique_ptr<InDet::ITRT_TrackExtensionTool::IEventData> ext_event_data_p( m_trtExtension->newEvent(ctx) );
   std::unique_ptr<TrackCollection> outTracks = std::make_unique<TrackCollection>();
   std::vector<Trk::Track*> tempTracks;                           //Temporary track collection
   tempTracks.reserve(128);
   // loop over event
   ATH_MSG_DEBUG ("Begin looping over all TRT segments in the event");
   Trk::SegmentCollection::const_iterator iseg    = segments->begin();
   Trk::SegmentCollection::const_iterator isegEnd = segments->end();
   InDet::TRT_DetElementLink_xk::TRT_DetElemUsedMap map;
   for(; iseg != isegEnd; ++ iseg) {
     // Get the track segment
     const Trk::TrackSegment *trackTRT = dynamic_cast<const Trk::TrackSegment*>(*iseg);
     if(!trackTRT){
       ATH_MSG_ERROR ("No pointer to segment !");
       continue;
     } else {
       // the segment finder is applying a TRT(eta) cut and a pt preselection, so we don't do that here
       //Ask for at least 10 TRT hits in order to process
       if(trackTRT->numberOfMeasurementBases() < m_minTRTonSegment) {
         ATH_MSG_DEBUG ("TRT segment fails nTRT hit cut, reject.");
         // statistics
         ev_stat.m_counter[Stat_t::kNTrtFailSel]++;
       } else {
         // do we continue to process ?
         nTrtSegCur++;
         if(nTrtSegCur>=m_MaxSegNum) {
           ATH_MSG_DEBUG ("====> Reached maximal number of segments in event, stop !!!");
           // statistics
           ev_stat.m_counter[Stat_t::kNTrtLimit]++;
           break;
         }
         // Get the number of the TRT track segment ROTs
         ATH_MSG_DEBUG ("=> New segment to process, number Of TRT ROTs : " << (trackTRT->numberOfMeasurementBases()));
         // statistics
         ev_stat.m_counter[Stat_t::Stat_t::kNTrtSegGood]++;
         // ok, call track maker and get list of possible track candidates
         std::list<Trk::Track*> trackSi = m_trackmaker->getTrack(ctx, *event_data_p, *trackTRT); //Get the possible Si extensions
         if (trackSi.empty()) {
           ATH_MSG_DEBUG ("No Si track candidates associated to the TRT track ");
           // statistics
           ev_stat.m_counter[Stat_t::kNTrtNoSiExt]++;
           // obsolete backup of TRT only
           if(m_saveTRT && trackTRT->numberOfMeasurementBases() > m_minTRTonly){
             Trk::Track* trtSeg = nullptr;trtSeg = segToTrack(ctx, *trackTRT);
             if(!trtSeg) {
               ATH_MSG_DEBUG ("Failed to make a track out of the TRT segment!");
               continue;
             }
             // statistics
             ev_stat.m_counter[Stat_t::kNBckTrk]++; ev_stat.m_counter[Stat_t::Stat_t::kNBckTrkTrt]++;
             // add track to output list
             outTracks->push_back(trtSeg);
           }
           continue;
         } else {
           // Found useful extensions
           ATH_MSG_DEBUG ("Found " << (trackSi.size()) << " Si tracks associated to the TRT track ");
           // Merge the resolved Si extensions with the original TRT track segment
           std::list<Trk::Track*>::const_iterator itt    = trackSi.begin();
           std::list<Trk::Track*>::const_iterator ittEnd = trackSi.end();
           for (; itt != ittEnd ; ++itt){
             tempTracks.push_back(*itt);
             // get list of TSOS
             const Trk::TrackStates* temptsos = (*itt)->trackStateOnSurfaces();
             if (!temptsos) {
               ATH_MSG_DEBUG ("Silicon extension empty ???");
               continue;
             }
             // Add the track to the list of tracks in the event
             ATH_MSG_DEBUG ("Silicon extension found has length of : " << temptsos->size());
             // do we do a preselection ?
             if (m_SiExtensionCuts) {
               // get parameters without errors
               auto input = (*itt)->trackParameters()->front()->uniqueClone();
               // cuts on parameters
               if (std::abs(input->pT()) < m_minPt) {
                 ATH_MSG_DEBUG ("Track pt < "<<m_minPt<<", reject it");
                 // statistics
                 ev_stat.m_counter[Stat_t::kNExtCut]++;
                 continue;
               }
               if (std::abs(input->eta()) > m_maxEta) {
                 ATH_MSG_DEBUG ("Track eta > "<<m_maxEta<<", reject it");
                 // statistics
                 ev_stat.m_counter[Stat_t::kNExtCut]++;
                 continue;
               }
               // --- beam spot position
               Amg::Vector3D beamSpotPosition(0,0,0);
               if (m_SiExtensionCuts){
                        SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey, ctx };
                        beamSpotPosition = beamSpotHandle->beamVtx().position();
               }
               // --- create surface
               Trk::PerigeeSurface perigeeSurface(beamSpotPosition);
  
               // uses perigee on track or extrapolates, no material in any case, we cut on impacts
               std::unique_ptr<const Trk::TrackParameters> parm =
                 m_extrapolator->extrapolateDirectly(
                   ctx, *input, perigeeSurface);
               std::unique_ptr<const Trk::Perigee> extrapolatedPerigee = nullptr;
               if (parm && parm->associatedSurface().type() == Trk::SurfaceType::Perigee) {
                 extrapolatedPerigee.reset(static_cast<const Trk::Perigee*>(parm.release()));
               }
               if (!extrapolatedPerigee) {
                 ATH_MSG_WARNING("Extrapolation of perigee failed, this should never happen" );
                 // statistics
                 ev_stat.m_counter[Stat_t::kNExtCut]++;
                 continue;
               }
  
               ATH_MSG_VERBOSE ("extrapolated perigee: "<<*extrapolatedPerigee);
               if (std::abs(extrapolatedPerigee->parameters()[Trk::d0]) > m_maxRPhiImp) {
                 ATH_MSG_DEBUG ("Track Rphi impact > "<<m_maxRPhiImp<<", reject it");
                 // statistics
                 ev_stat.m_counter[Stat_t::kNExtCut]++;
                 continue;
               }
               if (std::abs(extrapolatedPerigee->parameters()[Trk::z0]) > m_maxZImp) {
                 ATH_MSG_DEBUG ("Track Z impact > "<<m_maxZImp<<", reject it");
                 // statistics
                 ev_stat.m_counter[Stat_t::kNExtCut]++;
                 continue;
               }
             }
  
             // do re run a Track extension into TRT ?
             Trk::Track* globalTrackNew = nullptr;
             // do we have 4 and extension is enabled ?
             if(int(temptsos->size())>=4 && m_doExtension){
               // Add the track to the list of tracks in the event
               ATH_MSG_DEBUG ("Try to improve TRT calling extension tool.");
               // statistics
               ev_stat.m_counter[Stat_t::Stat_t::kNTrtExtCalls]++;
               // call extension tool
               std::vector<const Trk::MeasurementBase*>& tn =
                   m_trtExtension->extendTrack(ctx, *(*itt), *ext_event_data_p, map);
               if (tn.empty()) {
                 // Fallback if extension failed
                 ATH_MSG_DEBUG ("No new segment found, use input segment as fallback.");
                 // statistics
                 ev_stat.m_counter[Stat_t::Stat_t::kNTrtExtFail]++;
                 // merge Si with input track segments
                 globalTrackNew = mergeSegments(**itt,*trackTRT);
               } else if (!m_rejectShortExten ||
                          tn.size() >= trackTRT->numberOfMeasurementBases()) {
                 // Use the extension to instead of the segment
                 ATH_MSG_DEBUG ("Successful extension, number of TRT hits : " << tn.size() << " was : " << (trackTRT->numberOfMeasurementBases()));
                 // merge the extension with the Si track
                 globalTrackNew = mergeExtension(**itt,tn);
                 // Add the track to the list of tracks in the event
                 ATH_MSG_DEBUG ("Merged extension with Si segment");
                 // statistics
                 ev_stat.m_counter[Stat_t::kNTrtExt]++;
                 // clean up
                 std::vector<const Trk::MeasurementBase*>::const_iterator iv, ive=tn.end();
                 for(iv=tn.begin(); iv!=ive; ++iv) delete (*iv);
               } else {
                 // Extension is shorter, let's fall back onto the original
                 ATH_MSG_DEBUG ("Extension too short, number of TRT hits : " << tn.size() << " was : " << (trackTRT->numberOfMeasurementBases()) << ". Use Segement !");
                 // merge segments
                 globalTrackNew = mergeSegments(**itt,*trackTRT);
                 // Add the track to the list of tracks in the event
                 ATH_MSG_DEBUG ("Merged TRT segment with Si segment");
                 // statistics
                 ev_stat.m_counter[Stat_t::Stat_t::kNTrtExtBad]++;
                 // clean up
                 std::vector<const Trk::MeasurementBase*>::const_iterator iv, ive=tn.end();
                 for(iv=tn.begin(); iv!=ive; ++iv) delete (*iv);
               }
             } else {
               // no extension tool, jsut add the two
               ATH_MSG_DEBUG ("Do not try to extend Si track, merging it with input TRT.");
               // merge segments
               globalTrackNew = mergeSegments(**itt,*trackTRT);
             }
             // do we have an track candidate ?
             if(!globalTrackNew){
               ATH_MSG_DEBUG ("Failed to merge TRT+Si track segment !");
               if(m_saveTRT && trackTRT->numberOfMeasurementBases() > m_minTRTonly) {
                 Trk::Track* trtSeg = nullptr;trtSeg = segToTrack(ctx, *trackTRT);
                 if(!trtSeg){
                   ATH_MSG_DEBUG ("Failed to make a track out of the  TRT segment!");
                   continue;
                 }
                 ATH_MSG_DEBUG ("Add TRT only to output list");
                 // statistis
                 ev_stat.m_counter[Stat_t::kNBckTrk]++; ev_stat.m_counter[Stat_t::Stat_t::kNBckTrkTrt]++;
                 // add it to output list
                 if (m_trackSummaryTool.isEnabled()) {
                   m_trackSummaryTool->computeAndReplaceTrackSummary(*trtSeg,
                                       false /* DO NOT suppress hole search*/);
                 }
                 outTracks->push_back(trtSeg);
               }
             } else {
               ATH_MSG_DEBUG ("Save merged TRT+Si track segment!");
               // statistics
               ev_stat.m_counter[Stat_t::kNBckTrk]++; ev_stat.m_counter[Stat_t::Stat_t::kNBckTrkSi]++;
               // add it to output list
               if (m_trackSummaryTool.isEnabled()) {
                 m_trackSummaryTool->computeAndReplaceTrackSummary(*globalTrackNew,
                                     false /* DO NOT suppress hole search*/);
               }
               outTracks->push_back(globalTrackNew);
             }
           }
         }
       }
     }
   }
   // further debugging of results
   if(m_doStat){
     Analyze(outTracks.get());
   }
   if (SG::WriteHandle<TrackCollection>(m_outTracksKey,ctx).record(std::move(outTracks)).isFailure()){
      ATH_MSG_ERROR("Failed to record " << m_outTracksKey.key());
      return StatusCode::FAILURE;
   }
   // Update the total counters
   {
      std::lock_guard<std::mutex> lock(m_statMutex);
      m_totalStat += ev_stat;
   }
  
   for (auto *p : tempTracks){
     delete p;
   }
   m_trackmaker->endEvent(*event_data_p);
  
   //Print common event information
   if(msgLvl(MSG::DEBUG)){
      dumpevent(msg(MSG::DEBUG), ev_stat);
   }
  
   return StatusCode::SUCCESS;
 }

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

◆ filterPassed()

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

inlinevirtualinherited

Definition at line 135 of file AthReentrantAlgorithm.h.

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

◆ finalize()

StatusCode InDet::TRT_SeededTrackFinder::finalize ( )

override

Definition at line 448 of file TRT_SeededTrackFinder.cxx.

 {
   if(msgLvl(MSG::INFO)){
      msg(MSG::INFO) << "\n";
      dumpevent(msg(MSG::INFO), m_totalStat);
      msg(MSG::INFO) << endmsg;
   }
   return StatusCode::SUCCESS;
 }

◆ initialize()

StatusCode InDet::TRT_SeededTrackFinder::initialize ( )

override

TRT segments to use

Definition at line 82 of file TRT_SeededTrackFinder.cxx.

 {
  
   //Get the TRT seeded track maker tool
   //
   ATH_CHECK(m_trackmaker.retrieve());
  
   //Get the refitting tool
   //
   ATH_CHECK( m_fitterTool.retrieve(      DisableTool{ !m_doRefit } ));
   ATH_CHECK( m_extrapolator.retrieve(    DisableTool{ !m_SiExtensionCuts } ));
   ATH_CHECK( m_beamSpotKey.initialize(   m_SiExtensionCuts) );
  
   // optional PRD to track association map
   ATH_CHECK( m_prdToTrackMap.initialize( !m_prdToTrackMap.key().empty() ) );
  
   ATH_CHECK( m_trackSummaryTool.retrieve( DisableTool{ m_trackSummaryTool.name().empty() } ) );
  
   // Get tool for track extension to TRT
   //
   ATH_CHECK( m_trtExtension.retrieve(    DisableTool{ !m_doExtension} ));
  
  
   ATH_CHECK(  m_SegmentsKey.initialize()) ;  
   ATH_CHECK( m_outTracksKey.initialize());
  
   ATH_CHECK( m_caloClusterROIKey.initialize(m_caloSeededRoI) );
  
   if(m_caloSeededRoI){
     ATH_CHECK( m_regionSelector.retrieve());
   } else {
     m_regionSelector.disable();
   }
  
   // Get output print level
   //
   if(msgLvl(MSG::DEBUG)) {
      dumptools(msg(MSG::DEBUG));
   }
    //Global counters. See the include file for definitions
   m_totalStat = Stat_t();
   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()

bool AthReentrantAlgorithm::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in Simulation::BeamEffectsAlg, InDet::SiTrackerSpacePointFinder, InDet::SCT_Clusterization, InDet::SiSPSeededTrackFinder, SCTRawDataProvider, InDet::GNNSeedingTrackMaker, SCT_PrepDataToxAOD, RoIBResultToxAOD, SCT_CablingCondAlgFromCoraCool, SCT_ReadCalibDataTestAlg, SCT_CablingCondAlgFromText, InDet::SiSPGNNTrackMaker, SCT_ReadCalibChipDataTestAlg, SCT_TestCablingAlg, SCT_ConfigurationConditionsTestAlg, ITkPixelCablingAlg, SCTEventFlagWriter, SCT_ConditionsSummaryTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_TdaqEnabledTestAlg, SCT_SiliconConditionsTestAlg, SCTSiLorentzAngleTestAlg, SCT_ByteStreamErrorsTestAlg, SCT_ConditionsParameterTestAlg, SCT_FlaggedConditionTestAlg, SCT_StripVetoTestAlg, SCT_RawDataToxAOD, and SCTSiPropertiesTestAlg.

Definition at line 44 of file AthReentrantAlgorithm.cxx.

 {
   // Reentrant algorithms are clonable.
   return true;
 }

◆ mergeExtension()

Trk::Track * InDet::TRT_SeededTrackFinder::mergeExtension	(	const Trk::Track &	tT,
		std::vector< const Trk::MeasurementBase * > &	tS
	)		const

protected

Merge a TRT track extension and a Si track component into one global ID track.

Construct the new track

Definition at line 652 of file TRT_SeededTrackFinder.cxx.

                                                                                    {
   // TSOS from the track
   const Trk::TrackStates* stsos = tT.trackStateOnSurfaces();
   // fitQuality from track
   auto fq = tT.fitQuality()->uniqueClone();
   // output datavector of TSOS
   auto ntsos = std::make_unique<Trk::TrackStates>();
   // copy track Si states into track
   Trk::TrackStates::const_iterator p_stsos;
   for (p_stsos = stsos->begin(); p_stsos != stsos->end(); ++p_stsos) {
     ntsos->push_back((*p_stsos)->clone());
   }
   // loop over TRT track extension
   for (auto & it : tS) {
     std::bitset<Trk::TrackStateOnSurface::NumberOfTrackStateOnSurfaceTypes> typePattern;
     typePattern.set(Trk::TrackStateOnSurface::Measurement);
     const Trk::TrackStateOnSurface* seg_tsos = new Trk::TrackStateOnSurface(it->uniqueClone(), nullptr, nullptr, typePattern);
     ntsos->push_back(seg_tsos);
   }
   Trk::TrackInfo info;
   info.setPatternRecognitionInfo(Trk::TrackInfo::TRTSeededTrackFinder);
   std::unique_ptr<Trk::Track> newTrack( std::make_unique<Trk::Track>(info, std::move(ntsos), std::move(fq)) );
   //Careful refitting at the end
   if (m_doRefit) {
     newTrack = (m_fitterTool->fit(Gaudi::Hive::currentContext(),*newTrack, false, Trk::pion) ) ;
     if (!newTrack) {
       ATH_MSG_DEBUG ("Refit of TRT+Si track segment failed!");
       return nullptr;
     }
     //Protect for tracks that have no really defined locz and theta parameters
     const Trk::Perigee* perTrack=newTrack->perigeeParameters();
     if (perTrack) {
       const AmgSymMatrix(5)* CM = perTrack->covariance();
       if (!CM || std::sqrt((*CM)(1,1)) == 0. || std::sqrt((*CM)(3,3)) == 0.) {
         return nullptr;
       }
     }
   }
  
   return newTrack.release();
 }

◆ mergeSegments()

Trk::Track * InDet::TRT_SeededTrackFinder::mergeSegments	(	const Trk::Track &	tT,
		const Trk::TrackSegment &	tS
	)		const

protected

Protected methods

Merge a TRT track segment and a Si track component into one global ID track

Construct the new track

Definition at line 536 of file TRT_SeededTrackFinder.cxx.

                                                                                                        {
   // TSOS from the track
   const Trk::TrackStates* stsos = tT.trackStateOnSurfaces();
   // fitQuality from track
   auto fq = tT.fitQuality()->uniqueClone();
   // output datavector of TSOS
   auto ntsos = std::make_unique<Trk::TrackStates>();
   int siHits = 0;
   // copy track Si states into track
   Trk::TrackStates::const_iterator p_stsos;
   for (p_stsos=stsos->begin(); p_stsos != stsos->end(); ++p_stsos) {
     ntsos->push_back( (*p_stsos)->clone() );
     if ((*p_stsos)->type(Trk::TrackStateOnSurface::Measurement)) siHits++;
   }
   // loop over segment
   for (int it = 0; it < int(tS.numberOfMeasurementBases()); it++) {
     //test if it is a pseudo measurement
     if ( dynamic_cast<const Trk::PseudoMeasurementOnTrack*>(tS.measurement(it)) ) {
       if (siHits < 4) {
         ATH_MSG_DEBUG ("Too few Si hits.Will keep pseudomeasurement...");
         const Trk::TrackStateOnSurface* seg_tsos = new Trk::TrackStateOnSurface(tS.measurement(it)->uniqueClone(), nullptr);
         ntsos->push_back(seg_tsos);
       }
     } else {
       std::bitset<Trk::TrackStateOnSurface::NumberOfTrackStateOnSurfaceTypes> typePattern;
       typePattern.set(Trk::TrackStateOnSurface::Measurement);
       const Trk::TrackStateOnSurface* seg_tsos = new Trk::TrackStateOnSurface(tS.measurement(it)->uniqueClone(), nullptr, nullptr, typePattern);
       ntsos->push_back(seg_tsos);
     }
   }
  
   Trk::TrackInfo info;
   info.setPatternRecognitionInfo(Trk::TrackInfo::TRTSeededTrackFinder);
   std::unique_ptr<Trk::Track> newTrack(std::make_unique<Trk::Track>(info, std::move(ntsos), std::move(fq)));
  
   //Careful refitting at the end
   if (m_doRefit) {
     newTrack=m_fitterTool->fit(Gaudi::Hive::currentContext(),*newTrack, false, Trk::pion);
     if (!newTrack) {
       ATH_MSG_DEBUG ("Refit of TRT+Si track segment failed!");
       return nullptr;
     }
     const Trk::Perigee* perTrack=newTrack->perigeeParameters();
     if (perTrack) {
       const AmgSymMatrix(5)* CM = perTrack->covariance();
       if (!CM || std::sqrt((*CM)(1,1)) == 0. || std::sqrt((*CM)(3,3)) == 0.) {
         return nullptr;
       }
     }
   }
   return newTrack.release();
 }

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

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

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

◆ segToTrack()

Trk::Track * InDet::TRT_SeededTrackFinder::segToTrack	(	const EventContext &	ctx,
		const Trk::TrackSegment &	tS
	)		const

protected

Transform a TRT track segment into a track

Definition at line 594 of file TRT_SeededTrackFinder.cxx.

                                                                                                          {
   ATH_MSG_DEBUG ("Transforming the TRT segment into a track...");
  
   //Get the track segment information and build the initial track parameters
   const Trk::StraightLineSurface* surf = dynamic_cast<const Trk::StraightLineSurface*>(&(tS.associatedSurface()));
   if (!surf) {
     throw std::logic_error("Unhandled surface.");
   }
   const AmgVector(5)& p = tS.localParameters();
   AmgSymMatrix(5) ep = AmgSymMatrix(5)(tS.localCovariance());
   auto ntsos = std::make_unique<Trk::TrackStates>();
   std::unique_ptr<Trk::TrackParameters> segPar =
     surf->createUniqueParameters<5, Trk::Charged>(
       p(0), p(1), p(2), p(3), p(4), std::move(ep));
   if (segPar) {
     ATH_MSG_DEBUG( "Initial TRT Segment Parameters for refitting " << (*segPar) );
   } else {
     ATH_MSG_DEBUG( "Could not get initial TRT segment parameters! " );
     return nullptr;
   }
   for (int it = 0; it < int(tS.numberOfMeasurementBases()); it++) {
     // on first measurement add parameters
     const Trk::TrackStateOnSurface* seg_tsos = nullptr;
     std::bitset<Trk::TrackStateOnSurface::NumberOfTrackStateOnSurfaceTypes> typePattern;
     typePattern.set(Trk::TrackStateOnSurface::Measurement);
     if (it == 0){
       seg_tsos = new Trk::TrackStateOnSurface(tS.measurement(it)->uniqueClone(), std::move(segPar), nullptr, typePattern);
     } else {
       seg_tsos = new Trk::TrackStateOnSurface(tS.measurement(it)->uniqueClone(), nullptr, nullptr, typePattern);
     }
     ntsos->push_back(seg_tsos);
   }
   Trk::TrackInfo info;
   info.setPatternRecognitionInfo(Trk::TrackInfo::TRTSeededTrackFinder);
   std::unique_ptr<Trk::Track> newTrack = std::make_unique<Trk::Track>(info, std::move(ntsos), nullptr);
   // Careful refitting of the TRT stand alone track
   if (m_doRefit) {
     newTrack = m_fitterTool->fit(ctx,*newTrack, false, Trk::pion);
     if (!newTrack) {
       ATH_MSG_DEBUG ("Refit of TRT track segment failed!");
       return nullptr;
     }
     //Protect for tracks that have no really defined locz and theta parameters
     const Trk::Perigee* perTrack=newTrack->perigeeParameters();
     if (perTrack) {
       const AmgSymMatrix(5)* CM = perTrack->covariance();
       if (!CM || std::sqrt((*CM)(1,1)) == 0. || std::sqrt((*CM)(3,3)) == 0.) {
          return nullptr;
       }
     }
   }
   return newTrack.release();
 }

◆ 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 [1/2]

std::mutex m_statMutex InDet::TRT_SeededTrackFinder::ATLAS_THREAD_SAFE

mutableprotected

Definition at line 180 of file TRT_SeededTrackFinder.h.

◆ ATLAS_THREAD_SAFE [2/2]

Stat_t m_totalStat InDet::TRT_SeededTrackFinder::ATLAS_THREAD_SAFE

mutableprotected

Definition at line 181 of file TRT_SeededTrackFinder.h.

◆ m_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> InDet::TRT_SeededTrackFinder::m_beamSpotKey

protected

Initial value:

{
    this,
    "BeamSpotKey",
    "BeamSpotData",
    "SG key for beam spot"
  }

Definition at line 116 of file TRT_SeededTrackFinder.h.

◆ m_caloClusterROIKey

SG::ReadHandleKey<ROIPhiRZContainer> InDet::TRT_SeededTrackFinder::m_caloClusterROIKey

protected

Initial value:

{
    this,
    "EMROIPhiRZContainer",
    "",
    "Name of the calo cluster ROIs in Phi,R,Z parameterization"
  }

Definition at line 129 of file TRT_SeededTrackFinder.h.

◆ m_caloSeededRoI

bool InDet::TRT_SeededTrackFinder::m_caloSeededRoI

protected

Definition at line 128 of file TRT_SeededTrackFinder.h.

◆ m_deltaEta

float InDet::TRT_SeededTrackFinder::m_deltaEta

protected

delta Eta used for RoI creation

Definition at line 143 of file TRT_SeededTrackFinder.h.

◆ m_deltaPhi

float InDet::TRT_SeededTrackFinder::m_deltaPhi

protected

delta Phi used for RoI creation

Definition at line 144 of file TRT_SeededTrackFinder.h.

◆ m_deltaZ

float InDet::TRT_SeededTrackFinder::m_deltaZ

protected

delta Z used for RoI creation

Definition at line 145 of file TRT_SeededTrackFinder.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_doExtension

bool InDet::TRT_SeededTrackFinder::m_doExtension

protected

Do final careful refit of tracks.

Definition at line 85 of file TRT_SeededTrackFinder.h.

◆ m_doRefit

bool InDet::TRT_SeededTrackFinder::m_doRefit

protected

Definition at line 84 of file TRT_SeededTrackFinder.h.

◆ m_doStat

bool InDet::TRT_SeededTrackFinder::m_doStat

protected

use extension only if better than original track

Definition at line 87 of file TRT_SeededTrackFinder.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_extrapolator

ToolHandle<Trk::IExtrapolator> InDet::TRT_SeededTrackFinder::m_extrapolator

protected

the extrapoator

Definition at line 115 of file TRT_SeededTrackFinder.h.

◆ m_fitterTool

ToolHandle<Trk::ITrackFitter> InDet::TRT_SeededTrackFinder::m_fitterTool

protected

Track maker tool.

Definition at line 94 of file TRT_SeededTrackFinder.h.

◆ m_maxEta

double InDet::TRT_SeededTrackFinder::m_maxEta

protected

maximal Eta cut

Definition at line 124 of file TRT_SeededTrackFinder.h.

◆ m_maxRPhiImp

double InDet::TRT_SeededTrackFinder::m_maxRPhiImp

protected

maximal RPhi impact parameter cut

Definition at line 125 of file TRT_SeededTrackFinder.h.

◆ m_MaxSegNum

int InDet::TRT_SeededTrackFinder::m_MaxSegNum

protected

Output stand-alone TRT segments.

Definition at line 89 of file TRT_SeededTrackFinder.h.

◆ m_maxZImp

double InDet::TRT_SeededTrackFinder::m_maxZImp

protected

maximal z impact parameter cut

Definition at line 126 of file TRT_SeededTrackFinder.h.

◆ m_minPt

double InDet::TRT_SeededTrackFinder::m_minPt

protected

minimal Pt cut

Definition at line 123 of file TRT_SeededTrackFinder.h.

◆ m_minTRTonly

unsigned int InDet::TRT_SeededTrackFinder::m_minTRTonly

protected

Minimum number of TRT hits on segment.

Definition at line 91 of file TRT_SeededTrackFinder.h.

◆ m_minTRTonSegment

unsigned int InDet::TRT_SeededTrackFinder::m_minTRTonSegment

protected

Maximum number of segments to be handled.

Definition at line 90 of file TRT_SeededTrackFinder.h.

◆ m_outTracksKey

SG::WriteHandleKey<TrackCollection> InDet::TRT_SeededTrackFinder::m_outTracksKey

protected

TRT segments to use.

Definition at line 104 of file TRT_SeededTrackFinder.h.

◆ m_prdToTrackMap

SG::ReadHandleKey<Trk::PRDtoTrackMap> InDet::TRT_SeededTrackFinder::m_prdToTrackMap

protected

Initial value:

{ this,
                                                         "PRDtoTrackMap",
                                                         "" }

Definition at line 106 of file TRT_SeededTrackFinder.h.

◆ m_regionSelector

ToolHandle<IRegSelTool> InDet::TRT_SeededTrackFinder::m_regionSelector

protected

Initial value:

{
    this,
    "RegSelTool",
    "RegSelTool/RegSel_SCT",
    "Region selector service instance"
  }

Definition at line 136 of file TRT_SeededTrackFinder.h.

◆ m_rejectShortExten

bool InDet::TRT_SeededTrackFinder::m_rejectShortExten

protected

Find the TRT extension of the Si track segment.

Definition at line 86 of file TRT_SeededTrackFinder.h.

◆ m_saveTRT

bool InDet::TRT_SeededTrackFinder::m_saveTRT

protected

Statistics of final tracks.

Definition at line 88 of file TRT_SeededTrackFinder.h.

◆ m_SegmentsKey

SG::ReadHandleKey<Trk::SegmentCollection> InDet::TRT_SeededTrackFinder::m_SegmentsKey

protected

Definition at line 103 of file TRT_SeededTrackFinder.h.

◆ m_SiExtensionCuts

bool InDet::TRT_SeededTrackFinder::m_SiExtensionCuts

protected

enable cuts after Si segment finding

Definition at line 122 of file TRT_SeededTrackFinder.h.

◆ m_trackmaker

ToolHandle<ITRT_SeededTrackFinder> InDet::TRT_SeededTrackFinder::m_trackmaker

protected

Minimum number of TRT hits on TRT only.

Definition at line 93 of file TRT_SeededTrackFinder.h.

◆ m_trackSummaryTool

ToolHandle<Trk::IExtendedTrackSummaryTool> InDet::TRT_SeededTrackFinder::m_trackSummaryTool

protected

Initial value:

{
    this,
    "TrackSummaryTool",
    "InDetTrackSummaryToolNoHoleSearch"
  }

Definition at line 109 of file TRT_SeededTrackFinder.h.

◆ m_trtExtension

ToolHandle<ITRT_TrackExtensionTool> InDet::TRT_SeededTrackFinder::m_trtExtension

protected

Initial value:

{
    this,
    "TrackExtensionTool",
    "InDet::TRT_TrackExtensionTool_xk",
    "TRT track extension tool "
  }

Refitting tool.

Definition at line 95 of file TRT_SeededTrackFinder.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.

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

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ TRT_SeededTrackFinder()

◆ ~TRT_SeededTrackFinder()

Member Function Documentation

◆ Analyze()

◆ 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() [1/2]

◆ dump() [2/2]

◆ dumpevent()

◆ dumptools()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ filterPassed()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ mergeExtension()

◆ mergeSegments()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ segToTrack()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/2]

◆ ATLAS_THREAD_SAFE [2/2]

◆ m_beamSpotKey

◆ m_caloClusterROIKey

◆ m_caloSeededRoI

◆ m_deltaEta

◆ m_deltaPhi

◆ m_deltaZ

◆ m_detStore

◆ m_doExtension

◆ m_doRefit

◆ m_doStat

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_extrapolator

◆ m_fitterTool

◆ m_maxEta

◆ m_maxRPhiImp

◆ m_MaxSegNum

◆ m_maxZImp

◆ m_minPt

◆ m_minTRTonly

◆ m_minTRTonSegment

◆ m_outTracksKey