#include <AlignTrackPreProcessor.h>

Inheritance diagram for Trk::AlignTrackPreProcessor:

Collaboration diagram for Trk::AlignTrackPreProcessor:

Public Member Functions
	AlignTrackPreProcessor (const std::string &type, const std::string &name, const IInterface *parent)

virtual	~AlignTrackPreProcessor ()

StatusCode	initialize ()

StatusCode	finalize ()

DataVector< Track > *	processTrackCollection (const DataVector< Track > *trks)
	creates AlignTrack containing all TSOS on track More...

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

virtual void	setNtuple (TFile *)
	Sets output ntuple file, to be implemented if needed for detector-specific implementation. More...

virtual StatusCode	fillNtuple ()
	writes tree and histogrms to ntuple More...

virtual void	accumulateVTX (AlignTrack *)
	methods added for the full VTX fit: More...

virtual void	solveVTX ()

virtual void	printSummary ()
	Print processing summary to logfile. More...

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

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

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

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

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

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

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

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

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

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

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

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

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

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

MsgStream &	msg () const

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

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

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

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

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

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

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Track *	performSiliconHitSelection (const Track *, const ToolHandle< Trk::IGlobalTrackFitter > &)
	select silicon hits by quality. More...

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

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

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

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

Private Attributes
ToolHandle< IGlobalTrackFitter >	m_trackFitterTool

ToolHandle< IGlobalTrackFitter >	m_SLTrackFitterTool

ToolHandle< InDet::IInDetTrackSelectionTool >	m_trackSelectorTool

ToolHandle< IInDetAlignHitQualSelTool >	m_hitQualityTool

bool	m_refitTracks

bool	m_storeFitMatricesAfterRefit
	flag to refit tracks More...

bool	m_runOutlierRemoval
	flag to store derivative and covariance matrices after refit More...

int	m_particleHypothesis
	run outlier removal in track refit More...

bool	m_useSingleFitter
	particle hypothesis in track refit More...

bool	m_selectTracks
	do the track selection More...

bool	m_selectHits
	perform the hit InnerDetector selection More...

bool	m_fixMomentum
	Fix the momentum of the track so it is not refitted. 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

Definition at line 37 of file AlignTrackPreProcessor.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ AlignTrackPreProcessor()

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

Definition at line 30 of file AlignTrackPreProcessor.cxx.

     : AthAlgTool(type,name,parent)
     , m_trackFitterTool("Trk::GlobalChi2Fitter/InDetTrackFitter")
     , m_SLTrackFitterTool("")
     , m_trackSelectorTool("")
     , m_runOutlierRemoval(false)
     , m_particleHypothesis(Trk::nonInteracting)
     , m_useSingleFitter(false)
     , m_selectHits(false)
     , m_fixMomentum(false)
   {
     declareInterface<IAlignTrackPreProcessor>(this);
  
     declareProperty("RefitTracks", m_refitTracks = true);
  
     declareProperty("TrackFitterTool",   m_trackFitterTool);
     declareProperty("SLTrackFitterTool", m_SLTrackFitterTool);
     declareProperty("UseSingleFitter",   m_useSingleFitter = false);
  
     declareProperty("StoreFitMatricesAfterRefit", m_storeFitMatricesAfterRefit = true);
  
     declareProperty("TrackSelectorTool", m_trackSelectorTool);
     declareProperty("SelectTracks", m_selectTracks = false);
  
     declareProperty("ParticleHypothesis", m_particleHypothesis);
     declareProperty("RunOutlierRemoval",  m_runOutlierRemoval = false);
  
     declareProperty("HitQualityTool", m_hitQualityTool);
     declareProperty("SelectHits", m_selectHits);
     declareProperty("FixMomentum", m_fixMomentum);
  
     m_logStream = nullptr;
   }

◆ ~AlignTrackPreProcessor()

AlignTrackPreProcessor::~AlignTrackPreProcessor ( )

virtualdefault

Member Function Documentation

◆ accumulateVTX()

virtual void Trk::IAlignTrackPreProcessor::accumulateVTX ( AlignTrack * )

inlinevirtualinherited

methods added for the full VTX fit:

Reimplemented in Trk::BeamspotVertexPreProcessor.

Definition at line 57 of file IAlignTrackPreProcessor.h.

57 {}

◆ declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

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

inlineprivateinherited

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

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 245 of file AthCommonDataStore.h.

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

◆ declareProperty() [2/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 221 of file AthCommonDataStore.h.

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

◆ declareProperty() [3/6]

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

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

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

◆ declareProperty() [4/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 333 of file AthCommonDataStore.h.

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

◆ declareProperty() [5/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 352 of file AthCommonDataStore.h.

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

◆ declareProperty() [6/6]

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                  {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

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

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ fillNtuple()

virtual StatusCode Trk::IAlignTrackPreProcessor::fillNtuple ( )

inlinevirtualinherited

writes tree and histogrms to ntuple

Definition at line 54 of file IAlignTrackPreProcessor.h.

54 { return StatusCode::SUCCESS; }

◆ finalize()

StatusCode AlignTrackPreProcessor::finalize ( )

Definition at line 121 of file AlignTrackPreProcessor.cxx.

   {
     return StatusCode::SUCCESS;
   }

◆ initialize()

StatusCode AlignTrackPreProcessor::initialize ( )

Definition at line 71 of file AlignTrackPreProcessor.cxx.

   {
     // get track fitter tools
     if (m_trackFitterTool.retrieve().isSuccess())
       ATH_MSG_INFO("Retrieved " << m_trackFitterTool);
     else{
       msg(MSG::FATAL) << "Could not get " << m_trackFitterTool << endmsg;
       return StatusCode::FAILURE;
     }
  
     if (m_useSingleFitter) {
       
       if (m_SLTrackFitterTool.retrieve().isSuccess())
         ATH_MSG_INFO("Retrieved " << m_SLTrackFitterTool);
       else {
         msg(MSG::FATAL) << "Could not get " << m_SLTrackFitterTool << endmsg;
         return StatusCode::FAILURE;
       }
     }
  
     if(m_selectTracks) {
       if(m_trackSelectorTool.empty()) {
         msg(MSG::FATAL) << "TrackSelectorTool not specified : " << m_trackSelectorTool << endmsg;
         return StatusCode::FAILURE;
       }
       else if(m_trackSelectorTool.retrieve().isFailure())
       {
         msg(MSG::FATAL) << "Could not get " << m_trackSelectorTool << endmsg;
         return StatusCode::FAILURE;
       }
       ATH_MSG_INFO("Retrieved " << m_trackSelectorTool);
     }
  
     if (m_selectHits) {
       if(m_hitQualityTool.empty()) {
   msg(MSG::FATAL) << "HitQualityTool not specified : " << m_hitQualityTool << endmsg;
   return StatusCode::FAILURE;
       }
       else if(m_hitQualityTool.retrieve().isFailure())
   {
     msg(MSG::FATAL) << "Could not get " << m_trackSelectorTool << endmsg;
     return StatusCode::FAILURE;
   }
       ATH_MSG_INFO("Retrieved " << m_hitQualityTool);
     }
  
     return StatusCode::SUCCESS;
   }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

const InterfaceID & IAlignTrackPreProcessor::interfaceID ( )

inlinestaticinherited

Definition at line 70 of file IAlignTrackPreProcessor.h.

   {
     return IID_Trk_IAlignTrackPreProcessor;
   }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ performSiliconHitSelection()

Track * AlignTrackPreProcessor::performSiliconHitSelection	(	const Track *	inputTrack,
		const ToolHandle< Trk::IGlobalTrackFitter > &	fitter
	)

private

select silicon hits by quality.

keep all the rest

select silicon hits by quality. keep all the rest

Definition at line 239 of file AlignTrackPreProcessor.cxx.

   {
     ATH_MSG_DEBUG(" -- performSiliconHitSelection -- before removing bad Silicon hits, this track has "<< inputTrack->trackStateOnSurfaces()->size()<< " tsos");
     Track * newTrack;
  
     std::vector<const Trk::MeasurementBase*> selectedMeasurementSet;
  
     // loop on track hits
     int nhits = 0;
     for (const Trk::TrackStateOnSurface* tsos : *inputTrack->trackStateOnSurfaces())
     {
       nhits++;
       if (m_hitQualityTool->isGoodSiHit(tsos)) {
         selectedMeasurementSet.push_back( tsos->measurementOnTrack() );
       }
       else {
         ATH_MSG_DEBUG(" -- performSiliconHitSelection -- hit # "<< nhits << " status = BAD HIT ");
       }
     }
     ATH_MSG_DEBUG(" -- performSiliconHitSelection -- after removing bad Silicon hits, the selected measurement collection has "<< selectedMeasurementSet.size()<< " elements");
  
     newTrack = (fitter->fit(Gaudi::Hive::currentContext(),
                             selectedMeasurementSet,
                             *inputTrack->perigeeParameters(),
                             m_runOutlierRemoval,ParticleHypothesis(m_particleHypothesis))).release();
  
     return newTrack;
   }

◆ printSummary()

virtual void Trk::IAlignTrackPreProcessor::printSummary ( )

inlinevirtualinherited

Print processing summary to logfile.

Reimplemented in Trk::BeamspotVertexPreProcessor.

Definition at line 62 of file IAlignTrackPreProcessor.h.

62 {}

◆ processTrackCollection()

DataVector< Track > * AlignTrackPreProcessor::processTrackCollection ( const DataVector< Track > * trks )

virtual

creates AlignTrack containing all TSOS on track

select silicon hits by quality. keep all the rest

Implements Trk::IAlignTrackPreProcessor.

Definition at line 127 of file AlignTrackPreProcessor.cxx.

   {
     ATH_MSG_DEBUG("AlignTrackPreProcessor::processTrackCollection()");
  
     if (!tracks || tracks->empty())
       return nullptr;
  
     // the output collection of AlignTracks
     // we define it as collection of Tracks but fill AlignTracks inside
     DataVector<Track> * newTracks = new DataVector<Track>;
  
     int itrk(0);
     // loop over tracks and create AlignTracks
     TrackCollection::const_iterator it     = tracks->begin();
     TrackCollection::const_iterator it_end = tracks->end();
  
     for ( ; it != it_end ; ++it, ++itrk) {
  
       ATH_MSG_DEBUG(" ** processTrackCollection ** Processing track "<<itrk);
       const Track * origTrack = *it;
       const Track * newTrack = *it;
       AlignTrack* at;
  
       // check whether the original track passes selection
       if (m_selectTracks) {
         if(!m_trackSelectorTool->accept(*origTrack)) {
           ATH_MSG_VERBOSE(" ** processTrackCollection ** Original track did not pass the selection."<<itrk);
           continue;
         }
  
         ToolHandle<Trk::IGlobalTrackFitter> fitter=m_trackFitterTool;
         if (!m_useSingleFitter && AlignTrack::isSLTrack(origTrack) )
           fitter = m_SLTrackFitterTool;
  
         if (m_selectHits) {
           ATH_MSG_DEBUG(" ** processTrackCollection ** entering the silicon hit selection ");
           newTrack = performSiliconHitSelection(origTrack, fitter);
  
           if (!newTrack) {
             ATH_MSG_DEBUG(" ** processTrackCollection ** Track refit yielded no track. Skipping the track.");
             continue;
           }
  
           // check whether the track passes selection
           if(!m_trackSelectorTool->accept(*newTrack)) {
             ATH_MSG_DEBUG(" ** processTrackCollection ** Track did not pass the selection.");
             delete newTrack;
             continue;
           }
         }
  
         IGlobalTrackFitter::AlignmentCache alignCache;
         // refit track
         if (m_refitTracks &!m_selectHits) {
  
           newTrack=fitter->alignmentFit(alignCache,*origTrack,m_runOutlierRemoval,ParticleHypothesis(m_particleHypothesis));
           if (!newTrack) {
             ATH_MSG_DEBUG("Track refit yielded no track. Skipping the track.");
             continue;
           }
           // check that the refitted track satisfies the aligntrack selection
           if(!m_trackSelectorTool->accept(*newTrack)) {
             ATH_MSG_DEBUG("New track did not pass the selection.");
             delete newTrack;
             continue;
           }
         }
  
         at = new AlignTrack(*newTrack);
  
         if (msgLvl(MSG::DEBUG) && !msgLvl(MSG::VERBOSE)) {
           msg(MSG::DEBUG)<<"before refit: "<<endmsg;
           AlignTrack::dumpLessTrackInfo(*origTrack,msg(MSG::DEBUG));
           msg(MSG::DEBUG)<<"after refit: "<<endmsg;
           AlignTrack::dumpLessTrackInfo(*newTrack,msg(MSG::DEBUG));
           msg(MSG::DEBUG)<<endmsg;
         }
  
         // store fit matrices
         if (m_storeFitMatricesAfterRefit) {
           at->setFullCovarianceMatrix(alignCache.m_fullCovarianceMatrix);
           at->setDerivativeMatrix(alignCache.m_derivMatrix);
         }
         if (m_fixMomentum)
           {
              at->AlignTrack::setRefitQovP(false);
           }
         // delete newTrack since it's copied in AlignTrack
         delete newTrack;
       }
       else { // in case no selection is performed, keep all tracks
         at=new AlignTrack(*newTrack);
       }
     if (m_fixMomentum)
        {
          at->AlignTrack::setRefitQovP(false);
        }
  
       newTracks->push_back(at);
     }
  
  
     if (newTracks->empty()) {
       delete newTracks;
       return nullptr;
     }
  
     return newTracks;
   }

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ setLogStream()

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

inlinevirtualinherited

sets the output stream for the logfile

Definition at line 48 of file IAlignTrackPreProcessor.h.

48 { m_logStream = os; }

◆ setNtuple()

virtual void Trk::IAlignTrackPreProcessor::setNtuple ( TFile * )

inlinevirtualinherited

Sets output ntuple file, to be implemented if needed for detector-specific implementation.

Definition at line 51 of file IAlignTrackPreProcessor.h.

51 {}

◆ solveVTX()

virtual void Trk::IAlignTrackPreProcessor::solveVTX ( )

inlinevirtualinherited

Reimplemented in Trk::BeamspotVertexPreProcessor.

Definition at line 59 of file IAlignTrackPreProcessor.h.

59 {}

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_fixMomentum

bool Trk::AlignTrackPreProcessor::m_fixMomentum

private

Fix the momentum of the track so it is not refitted.

Definition at line 72 of file AlignTrackPreProcessor.h.

◆ m_hitQualityTool

ToolHandle<IInDetAlignHitQualSelTool> Trk::AlignTrackPreProcessor::m_hitQualityTool

private

Definition at line 57 of file AlignTrackPreProcessor.h.

◆ m_logStream

std::ostream* Trk::IAlignTrackPreProcessor::m_logStream

protectedinherited

logfile output stream

Definition at line 66 of file IAlignTrackPreProcessor.h.

◆ m_particleHypothesis

int Trk::AlignTrackPreProcessor::m_particleHypothesis

private

run outlier removal in track refit

Definition at line 66 of file AlignTrackPreProcessor.h.

◆ m_refitTracks

bool Trk::AlignTrackPreProcessor::m_refitTracks

private

Definition at line 62 of file AlignTrackPreProcessor.h.

◆ m_runOutlierRemoval

bool Trk::AlignTrackPreProcessor::m_runOutlierRemoval

private

flag to store derivative and covariance matrices after refit

Definition at line 65 of file AlignTrackPreProcessor.h.

◆ m_selectHits

bool Trk::AlignTrackPreProcessor::m_selectHits

private

perform the hit InnerDetector selection

Definition at line 71 of file AlignTrackPreProcessor.h.

◆ m_selectTracks

bool Trk::AlignTrackPreProcessor::m_selectTracks

private

do the track selection

Definition at line 70 of file AlignTrackPreProcessor.h.

◆ m_SLTrackFitterTool

ToolHandle<IGlobalTrackFitter> Trk::AlignTrackPreProcessor::m_SLTrackFitterTool

private

Definition at line 53 of file AlignTrackPreProcessor.h.

◆ m_storeFitMatricesAfterRefit

bool Trk::AlignTrackPreProcessor::m_storeFitMatricesAfterRefit

private

flag to refit tracks

Definition at line 63 of file AlignTrackPreProcessor.h.

◆ m_trackFitterTool

ToolHandle<IGlobalTrackFitter> Trk::AlignTrackPreProcessor::m_trackFitterTool

private

Definition at line 52 of file AlignTrackPreProcessor.h.

◆ m_trackSelectorTool

ToolHandle<InDet::IInDetTrackSelectionTool> Trk::AlignTrackPreProcessor::m_trackSelectorTool

private

Definition at line 56 of file AlignTrackPreProcessor.h.

◆ m_useSingleFitter

bool Trk::AlignTrackPreProcessor::m_useSingleFitter

private

particle hypothesis in track refit

only use 1 fitter for refitting track

Definition at line 68 of file AlignTrackPreProcessor.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ AlignTrackPreProcessor()

◆ ~AlignTrackPreProcessor()

Member Function Documentation

◆ accumulateVTX()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ fillNtuple()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ performSiliconHitSelection()

◆ printSummary()

◆ processTrackCollection()

◆ renounce()

◆ renounceArray()

◆ setLogStream()

◆ setNtuple()

◆ solveVTX()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_evtStore

◆ m_fixMomentum

◆ m_hitQualityTool

◆ m_logStream

◆ m_particleHypothesis

◆ m_refitTracks

◆ m_runOutlierRemoval

◆ m_selectHits

◆ m_selectTracks

◆ m_SLTrackFitterTool

◆ m_storeFitMatricesAfterRefit

◆ m_trackFitterTool

◆ m_trackSelectorTool

◆ m_useSingleFitter

◆ m_varHandleArraysDeclared

◆ m_vhka