#include <TRTTrackHoleSearchTool.h>

Inheritance diagram for TRTTrackHoleSearchTool:

Collaboration diagram for TRTTrackHoleSearchTool:

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

StatusCode	initialize ()

StatusCode	finalize ()

void	countHoles (const Trk::Track &track, std::vector< int > &information, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, partHyp Output: Changes in information This method fills the fields relevant to the hole counts in the vector information. More...

const Trk::TrackStates *	getHolesOnTrack (const Trk::Track &track, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, parthyp Return: A DataVector containing pointers to TrackStateOnSurfaces which each represent an identified hole on the track. More...

const Trk::Track *	getTrackWithHoles (const Trk::Track &track, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes The parthyp argument is relevant for the extrapolation steps in the hole search. More...

const Trk::Track *	getTrackWithHolesAndOutliers (const Trk::Track &track, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes or outliers The parthyp argument is relevant for the extrapolation steps in the hole search. 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...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
int	extrapolateBetweenHits (const Trk::TrackParameters start_parameters, const Trk::Surface &end_surf, Trk::TrackStates holes, const Trk::ParticleHypothesis partHyp=Trk::pion) const

void	dump_bad_straw_log () const

Trk::TrackStates::const_iterator	find_first_trt_hit (const Trk::TrackStates &track_states) const

Trk::TrackStates::const_iterator	find_last_hit_before_trt (const Trk::TrackStates &track_states) const

const Trk::Track *	addHolesToTrack (const Trk::Track &track, const Trk::TrackStates *holes) const

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< Trk::IExtrapolator >	m_extrapolator

ServiceHandle< IInDetConditionsSvc >	m_conditions_svc

bool	m_use_conditions_svc

float	m_outer_radius

float	m_max_z

int	m_max_trailing_holes

bool	m_begin_at_first_trt_hit

bool	m_end_at_last_trt_hit

bool	m_bcheck

bool	m_do_dump_bad_straw_log

float	m_locR_cut

float	m_locR_sigma_cut

float	m_locZ_cut

std::atomic_bool	m_has_been_called

const TRT_ID *	m_TRT_ID

Trk::CylinderSurface *	m_trt_outer_surf

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 35 of file TRTTrackHoleSearchTool.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

◆ TRTTrackHoleSearchTool()

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

Definition at line 39 of file TRTTrackHoleSearchTool.cxx.

     : AthAlgTool(type, name, parent),
       m_extrapolator("Trk::Extrapolator"),
       m_conditions_svc("TRT_ConditionsSummarySvc", name),
       m_has_been_called(false),
       m_TRT_ID(nullptr),
       m_trt_outer_surf(nullptr)
 {
     declareInterface<ITrackHoleSearchTool>(this);
  
     declareProperty("extrapolator",             m_extrapolator );
     declareProperty("conditions_svc",           m_conditions_svc);
     declareProperty("use_conditions_svc",       m_use_conditions_svc =  true);
     declareProperty("outer_radius",             m_outer_radius =        1075.0*CLHEP::mm );  // barrel 1075.0, EC 1010.0
     declareProperty("max_z",                    m_max_z =               2715.0*CLHEP::mm );  // barrel 715.0, EC 2715.0
     declareProperty("max_trailing_holes",       m_max_trailing_holes =  1 );   // only used if not end_at_last_trt_hit
     declareProperty("begin_at_first_trt_hit",   m_begin_at_first_trt_hit = false ); // if not, extrapolate from last Si hit
     declareProperty("end_at_last_trt_hit",      m_end_at_last_trt_hit = false ); // if not, continue hole search to the edge of the TRT
     declareProperty("bcheck",                   m_bcheck =              false );
     declareProperty("do_dump_bad_straw_log",    m_do_dump_bad_straw_log = false );
     declareProperty("locR_cut",                 m_locR_cut =            -1.0 );   // 1.4*CLHEP::mm // negative means no cut
     declareProperty("locR_sigma_cut",           m_locR_sigma_cut =      -1.0 );
     declareProperty("locZ_cut",                 m_locZ_cut =            5.0*CLHEP::mm ); // 5.0*CLHEP::mm );
 }

Member Function Documentation

◆ addHolesToTrack()

const Trk::Track * TRTTrackHoleSearchTool::addHolesToTrack	(	const Trk::Track &	track,
		const Trk::TrackStates *	holes
	)		const

private

Definition at line 440 of file TRTTrackHoleSearchTool.cxx.

 {
   ATH_MSG_DEBUG("TRTTrackHoleSearchTool::addHolesToTrack");
   /*
     This method was basically coppied from here:
     http://alxr.usatlas.bnl.gov/lxr-stb4/source/atlas/InnerDetector/InDetRecTools/InDetTrackHoleSearch/src/InDetTrackHoleSearchTool.cxx#931
   */
  
   // get states from track
   auto tsos = std::make_unique<Trk::TrackStates>();
   for (const auto *it : *track.trackStateOnSurfaces()) {
     // veto old holes
     if (!it->type(Trk::TrackStateOnSurface::Hole)) {
       tsos->push_back(new Trk::TrackStateOnSurface(*it));
     }
   }
  
   // if we have no holes on the old track and no holes found by search, then we
   // just copy the track
   if (track.trackStateOnSurfaces()->size() == tsos->size() && holes->empty()) {
     // create copy of track
     const Trk::Track* new_track = new Trk::Track(
       track.info(),
       std::move(tsos),
       track.fitQuality() ? track.fitQuality()->uniqueClone() : nullptr);
     return new_track;
   }
  
   // add new holes
   tsos->insert(tsos->end(), holes->begin(), holes->end());
  
   // sort
   const Trk::TrackParameters* perigee = track.perigeeParameters();
   if (!perigee)
     perigee = (*(track.trackStateOnSurfaces()->begin()))->trackParameters();
  
   if (perigee) {
     Trk::TrackStateOnSurfaceComparisonFunction CompFunc(perigee->momentum());
  
     if (fabs(perigee->parameters()[Trk::qOverP]) > 0.002) {
       /* invest n*(logN)**2 sorting time for lowPt, coping with a possibly
          not 100% transitive comparison functor. */
       if (msgLvl(MSG::DEBUG)) {
         msg() << "sorting vector with stable_sort" << endmsg;
       }
       std::stable_sort(tsos->begin(), tsos->end(), CompFunc);
     } else {
       tsos->sort(CompFunc); // respects DV object ownership
     }
   }
  
   // create copy of track
   const Trk::Track* new_track =
     new Trk::Track(track.info(),
                    std::move(tsos),
                    track.fitQuality() ? track.fitQuality()->uniqueClone() : nullptr);
  
   return new_track;
 }

◆ countHoles()

void TRTTrackHoleSearchTool::countHoles	(	const Trk::Track &	track,
		std::vector< int > &	information,
		const Trk::ParticleHypothesis	partHyp = `Trk::pion`
	)		const

virtual

Input : track, partHyp Output: Changes in information This method fills the fields relevant to the hole counts in the vector information.

These fields should be initialised to 0. The relevant indices are specified by the enumeration in Tracking/TrkEvent/TrkTrackSummary. If problems occur, the information counters are reset to -1 flagging them as not set. The parthyp argument is relevant for the extrapolation steps in the hole search.

Implements Trk::ITrackHoleSearchTool.

Definition at line 95 of file TRTTrackHoleSearchTool.cxx.

                                                                                    {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::countHoles" );
     const Trk::TrackStates* holes = getHolesOnTrack(track, partHyp);
     if (holes) {
         information[Trk::numberOfTRTHoles] = holes->size();
         delete holes;
     }   else {
         information[Trk::numberOfTRTHoles] = -1;
     }
 }

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

◆ dump_bad_straw_log()

void TRTTrackHoleSearchTool::dump_bad_straw_log ( ) const

private

Definition at line 372 of file TRTTrackHoleSearchTool.cxx.

                                                       {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::dump_bad_straw_log" );
     std::ofstream out("TRT_ConditionsSummarySvc_bad_straws.log");
     out << "# id  barrel_ec  phi_module  layer_or_wheel  straw_layer  straw" << std::endl;
     for(std::vector<Identifier>::const_iterator it = m_TRT_ID->straw_layer_begin(); it != m_TRT_ID->straw_layer_end(); ++it) {
         for(int i=0; i<= m_TRT_ID->straw_max(*it); ++i) {
             Identifier id = m_TRT_ID->straw_id(*it, i);
             if(!m_conditions_svc->isGood(id)) {
                 out << id.getString()
                     << std::setw(3) << m_TRT_ID->barrel_ec(id)
                     << std::setw(4) << m_TRT_ID->phi_module(id)
                     << std::setw(4) << m_TRT_ID->layer_or_wheel(id)
                     << std::setw(4) << m_TRT_ID->straw_layer(id)
                     << std::setw(4) << m_TRT_ID->straw(id)
                     << std::endl;
             }
         }
     }
     out.close();
 }

◆ 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

◆ extrapolateBetweenHits()

int TRTTrackHoleSearchTool::extrapolateBetweenHits	(	const Trk::TrackParameters *	start_parameters,
		const Trk::Surface &	end_surf,
		Trk::TrackStates *	holes,
		const Trk::ParticleHypothesis	partHyp = `Trk::pion`
	)		const

private

Definition at line 245 of file TRTTrackHoleSearchTool.cxx.

                                                                                                {
     int hole_count = 0;
     // initialize previous id
     const Trk::Surface& start_surf = start_parameters->associatedSurface();
  
     Identifier previous_id = start_surf.associatedDetectorElementIdentifier();
  
     // get end id
     const Identifier end_id = end_surf.associatedDetectorElementIdentifier();
  
     // look for holes
     std::vector<std::unique_ptr<Trk::TrackParameters> > steps =
           m_extrapolator->extrapolateStepwise(Gaudi::Hive::currentContext(),
                                               *start_parameters, 
                                               end_surf, 
                                               Trk::alongMomentum, m_bcheck, partHyp);
  
     if(steps.empty()) {
         ATH_MSG_DEBUG("extrapolateBetweenHits: extrapolateStepwise returned null");
     } else {
         // loop over parameters from extrapolation
         // note: the last element in the vector is always the track parameters at the destination surface
                 for(std::vector<std::unique_ptr<Trk::TrackParameters> >::const_iterator step = steps.begin(); step != steps.end()-1; ++step) {
             // check for surface
             const Trk::Surface& surf = (*step)->associatedSurface();
             
             // get id
             Identifier id = surf.associatedDetectorElementIdentifier();
  
             // check that the surface is a TRT straw
             if(!m_TRT_ID->is_trt(id)) {
                 ATH_MSG_DEBUG("extrapolateBetweenHits: surf is not a TRT straw.  Skipping.");
                 continue;
             }
  
             // check that we are not still on the previous id
             if(id == previous_id) {
                 ATH_MSG_DEBUG("extrapolateBetweenHits: id == previous_id");
                 continue;
             }
  
             // check that we haven't hit the destination surface
             if(id == end_id) {
                 ATH_MSG_DEBUG("extrapolateBetweenHits: id == end_id");
                 continue;
             }
  
             // ignore id 0xffffffff from m_trt_outer_surf
             if(id == 0xffffffff) {
                 ATH_MSG_DEBUG("extrapolateBetweenHits: id == 0xffffffff. Skipping.");
                 continue;
             }
  
             // don't count bad straws as holes
             if( (m_use_conditions_svc)&&(!m_conditions_svc->isGood(id)) ) {
                 ATH_MSG_DEBUG("extrapolateBetweenHits: ConditionsSvc says this straw is bad. Skipping.");
                 continue;
             }
  
             // check fiducial |locR| by sigma
             const float locR = (*step)->parameters()[Trk::locR];
             if(m_locR_sigma_cut > 0.0) {
                 //const Trk::MeasuredTrackParameters* meas = dynamic_cast< const Trk::MeasuredTrackParameters* >(*step);
                 //if(meas)
                 //{
  
                 const AmgSymMatrix(5)* merr = (*step)->covariance();
  
                 if(merr){
                     //const float locR_error = (*merr)(Trk::locR,Trk::locR);
                     
                     const float locR_error = Amg::error(*merr,Trk::locR);
  
                     // hole must be within the straw by some tolerance
                     if( (2.0*CLHEP::mm - fabs(locR))/locR_error < m_locR_sigma_cut ) {
                         continue;
                     }
                 } else {
                     ATH_MSG_WARNING("extrapolateBetweenHits: Track parameters failed to dynamic_cast< const Trk::MeasuredTrackParameters* >.");
                     continue;
                 }
             }
             // check fiducial |locR| by value
             if(m_locR_cut > 0.0) {
                 if( fabs(locR) > m_locR_cut )
                     continue;
             }
             // check fiducial |locZ| by value
             if(m_locZ_cut > 0.0) {
                 const float locZ = (*step)->parameters()[Trk::locZ];
                 const Trk::CylinderBounds* cylb = dynamic_cast<const Trk::CylinderBounds *>(&surf.bounds());
                 if(cylb) {
                     const float halfz=cylb->halflengthZ();
                     if( fabs(locZ) > halfz - m_locZ_cut ) continue;
                 } else {
                     ATH_MSG_WARNING("extrapolateBetweenHits: Surface failed to dynamic_cast to Trk::CylinderBounds. Skipping.");
                     continue;
                 }
             }
  
             // if we've gotten here, it is a legitimate hole
             ATH_MSG_DEBUG("      HOLE Found! Identifier = " << id.getString() << " indicates:\n"
                           << "            is_pixel =       " << m_TRT_ID->is_pixel(id) << "\n"
                           << "            is_sct =         " << m_TRT_ID->is_sct(id) << "\n"
                           << "            is_trt =         " << m_TRT_ID->is_trt(id) << "\n"
                           << "            barrel_ec =      " << m_TRT_ID->barrel_ec(id) << "\n"
                           << "            phi_module =     " << m_TRT_ID->phi_module(id) << "\n"
                           << "            layer_or_wheel = " << m_TRT_ID->layer_or_wheel(id) << "\n"
                           << "            straw_layer =    " << m_TRT_ID->straw_layer(id) << "\n"
                           << "            straw =          " << m_TRT_ID->straw(id) << "\n"
                           << "            rho = " << (*step)->position().perp() );
             
             std::bitset<Trk::TrackStateOnSurface::NumberOfTrackStateOnSurfaceTypes> typePattern;
             typePattern.set(Trk::TrackStateOnSurface::Hole);
             holes->push_back( new Trk::TrackStateOnSurface(nullptr, (*step)->uniqueClone(), nullptr, typePattern) );
             hole_count++;
             previous_id = id;
         } // end loop over parameters from extrapolation
     }
     return hole_count;
 }

◆ finalize()

StatusCode TRTTrackHoleSearchTool::finalize ( )

Definition at line 87 of file TRTTrackHoleSearchTool.cxx.

                                             {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::finalize()" );
     delete m_trt_outer_surf;
     return StatusCode::SUCCESS;
 }

◆ find_first_trt_hit()

Trk::TrackStates::const_iterator TRTTrackHoleSearchTool::find_first_trt_hit ( const Trk::TrackStates & track_states ) const

private

Definition at line 396 of file TRTTrackHoleSearchTool.cxx.

                                                                                    {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::find_first_trt_hit" );
     // finds the first TRT hit or returns end
     Trk::TrackStates::const_iterator track_state = track_states.begin();
     for(; track_state != track_states.end(); ++track_state) {
         // skip track states that are not measurements
         if(!(*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
             ATH_MSG_VERBOSE( "    TrackStateOnSurface is not of type Trk::TrackStateOnSurface::Measurement." );
             continue;
         }
         if( dynamic_cast<const InDet::TRT_DriftCircleOnTrack*>( (*track_state)->measurementOnTrack() )  ) {
             ATH_MSG_DEBUG("    Found first TRT hit.");
             break;
         }
     }
     return track_state;
 }

◆ find_last_hit_before_trt()

Trk::TrackStates::const_iterator TRTTrackHoleSearchTool::find_last_hit_before_trt ( const Trk::TrackStates & track_states ) const

private

Definition at line 417 of file TRTTrackHoleSearchTool.cxx.

                                                                                          {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::find_last_hit_before_trt" );
     // start at first TRT hit
     Trk::TrackStates::const_iterator track_state = find_first_trt_hit(track_states);
     // if there is no room to step backwards in the hits, return end
     if(track_states.size() < 2 || track_state == track_states.begin()) {
         return track_states.end();
     }
     --track_state; // step back and look for last measurement before the TRT hit
     for(; track_state != track_states.begin(); --track_state) {
         if((*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
             break;
         }
     }
     // if the loop ended without finding a measurement, return end
     if(!(*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
         return track_states.end();
     }
     return track_state;
 }

◆ getHolesOnTrack()

const Trk::TrackStates * TRTTrackHoleSearchTool::getHolesOnTrack	(	const Trk::Track &	track,
		const Trk::ParticleHypothesis	partHyp = `Trk::pion`
	)		const

virtual

Input : track, parthyp Return: A DataVector containing pointers to TrackStateOnSurfaces which each represent an identified hole on the track.

The parthyp argument is relevant for the extrapolation steps in the hole search. Attention: This is a factory, ownership of the return vector is passed to the calling method.

Implements Trk::ITrackHoleSearchTool.

Definition at line 110 of file TRTTrackHoleSearchTool.cxx.

                                                                    {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::getHolesOnTrack" );
  
     // write out list of bad straws for debugging purposes.
     // only dump log on first call of this function.
     if (m_do_dump_bad_straw_log && !m_has_been_called && m_use_conditions_svc) {
         dump_bad_straw_log();
     }
     m_has_been_called = true;
  
     if (track.perigeeParameters()) {
         ATH_MSG_VERBOSE( "  This track has perigee: \n"
                          << "    pT =  " << track.perigeeParameters()->pT()/CLHEP::GeV << " CLHEP::GeV\n"
                          << "    eta = " << track.perigeeParameters()->eta() << "\n"
                          << "    phi = " << track.perigeeParameters()->parameters()[Trk::phi0] << "\n"
                          << "    d0 =  " << track.perigeeParameters()->parameters()[Trk::d0] << "\n"
                          << "    z0 =  " << track.perigeeParameters()->parameters()[Trk::z0] );
     }
  
     // get TrackStateOnSurfaces
     const Trk::TrackStates* track_states = track.trackStateOnSurfaces();
     if (track_states) {
         ATH_MSG_DEBUG( "  This track has " << track_states->size() << " track states on surface." );
     } else {
         ATH_MSG_WARNING( "  This track has null track states on surface. Returning 0." );
         return nullptr;
     }
  
     if (track_states->size() < 2) {
         ATH_MSG_WARNING( "  Fewer than 2 TrackStatesOnSurface. Returning 0." );
         return nullptr;
     }
  
     // set beginning point of extrapolation
     Trk::TrackStates::const_iterator beginning_track_state;
     if (m_begin_at_first_trt_hit) {
         // begin at first TRT hit
         beginning_track_state = find_first_trt_hit(*track_states);
     } else {
         // begin at last Si hit
         beginning_track_state = find_last_hit_before_trt(*track_states);
         if (beginning_track_state == track_states->end()) {
             ATH_MSG_WARNING( "  beginning_track_state == track_states->end().  There must be no Si hits.\n"
                              << "  Will try to begin at the first TRT hit." );
             beginning_track_state = find_first_trt_hit(*track_states);
         }
     }
  
     if (beginning_track_state == track_states->end()) {
         ATH_MSG_WARNING( "  beginning_track_state == track_states->end(). No where to extrapolate to. Returning 0." );
         return nullptr;
     }
  
     // to be returned:
     Trk::TrackStates* holes = new Trk::TrackStates;
  
     Trk::TrackStates::const_iterator track_state = beginning_track_state;
     const Trk::TrackParameters* start_parameters = (*track_state)->trackParameters();
     ++track_state;
     // loop over TrackStateOnSurfaces (destination surfaces for extrapolation)
     for (; track_state != track_states->end(); ++track_state) {
         // skip track states that are not measurements
         if (!(*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
             ATH_MSG_VERBOSE( "    TrackStateOnSurface is not of type Trk::TrackStateOnSurface::Measurement." );
             continue;
         }
  
         const Trk::TrackParameters* end_parameters = (*track_state)->trackParameters();
         if(end_parameters) {
             ATH_MSG_VERBOSE( "    TrackStateOnSurface: (x, y, z) = (" 
                              << end_parameters->position().x() << ", "
                              << end_parameters->position().y() << ", "
                              << end_parameters->position().z() << "); (rho, eta, phi) = (" 
                              << end_parameters->position().perp() << ", "
                              << end_parameters->position().eta() << ", "
                              << end_parameters->position().phi() << ")");
         } else {
             ATH_MSG_WARNING( "    TrackStateOnSurface has no TrackParameters." );
             continue;
         }
  
         const Trk::Surface& end_surf = end_parameters->associatedSurface();
         if (!start_parameters) {
             start_parameters = end_parameters;
         } else {         
             extrapolateBetweenHits(start_parameters, end_surf, holes, partHyp);
             start_parameters = end_parameters;
         }
     } // end loop over TrackStateOnSurfaces
  
     if( !m_end_at_last_trt_hit ) {
         // final extrapolation to the edge of TRT to check for holes in outermost straws
         int trailing_hole_count = extrapolateBetweenHits(start_parameters, *m_trt_outer_surf, holes, partHyp);
  
         // remove trailing holes
         if(trailing_hole_count > m_max_trailing_holes) {
             ATH_MSG_DEBUG("There are " << trailing_hole_count << " trailing holes removed.");
             for(int i=0; i < trailing_hole_count; i++) {
                 holes->pop_back();
                 /*
                   Note: DataVector::pop_back() deletes pointers if it owns them.
                 */
             }
         }
     }
  
     return holes;
 }

◆ getTrackWithHoles()

const Trk::Track * TRTTrackHoleSearchTool::getTrackWithHoles	(	const Trk::Track &	track,
		const Trk::ParticleHypothesis	partHyp = `Trk::pion`
	)		const

virtual

Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes The parthyp argument is relevant for the extrapolation steps in the hole search.

Attention: This is a factory, ownership of the return track is passed to the calling method.

Implements Trk::ITrackHoleSearchTool.

Definition at line 223 of file TRTTrackHoleSearchTool.cxx.

                                                                                                        {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::getTrackWithHoles" );
     const Trk::TrackStates* holes = getHolesOnTrack(track, partHyp);
     const Trk::Track* new_track = addHolesToTrack(track, holes);
     delete holes;
     return new_track;
 }

◆ getTrackWithHolesAndOutliers()

const Trk::Track * TRTTrackHoleSearchTool::getTrackWithHolesAndOutliers	(	const Trk::Track &	track,
		const Trk::ParticleHypothesis	partHyp = `Trk::pion`
	)		const

virtual

Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes or outliers The parthyp argument is relevant for the extrapolation steps in the hole search.

Attention: This is a factory, ownership of the return track is passed to the calling method.

Implements Trk::ITrackHoleSearchTool.

Definition at line 234 of file TRTTrackHoleSearchTool.cxx.

                                                                                                                   {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::getTrackWithHolesAndOutliers" );
     return getTrackWithHoles(track, partHyp);
 }

◆ initialize()

StatusCode TRTTrackHoleSearchTool::initialize ( )

Definition at line 66 of file TRTTrackHoleSearchTool.cxx.

                                               {
     ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::initialize()" );
  
     // retrieve extrapolator
     ATH_CHECK(m_extrapolator.retrieve());
  
     // retrieve TRT_ID
     ATH_CHECK(detStore()->retrieve(m_TRT_ID, "TRT_ID"));
  
     // retrieve ConditionsSummarySvc
     ATH_CHECK(m_conditions_svc.retrieve());
  
     m_trt_outer_surf = new Trk::CylinderSurface(Amg::Transform3D(Amg::Transform3D::Identity()),
      m_outer_radius, m_max_z);
     // note: HepGeom::Translate3D is deleted by Trk::Surface destructor
  
     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 & Trk::ITrackHoleSearchTool::interfaceID ( )

inlinestaticinherited

Definition at line 78 of file ITrackHoleSearchTool.h.

     { 
       return IID_ITrackHoleSearchTool; 
     }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

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

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_bcheck

bool TRTTrackHoleSearchTool::m_bcheck

private

Definition at line 91 of file TRTTrackHoleSearchTool.h.

◆ m_begin_at_first_trt_hit

bool TRTTrackHoleSearchTool::m_begin_at_first_trt_hit

private

Definition at line 89 of file TRTTrackHoleSearchTool.h.

◆ m_conditions_svc

ServiceHandle<IInDetConditionsSvc> TRTTrackHoleSearchTool::m_conditions_svc

private

Definition at line 84 of file TRTTrackHoleSearchTool.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_do_dump_bad_straw_log

bool TRTTrackHoleSearchTool::m_do_dump_bad_straw_log

private

Definition at line 92 of file TRTTrackHoleSearchTool.h.

◆ m_end_at_last_trt_hit

bool TRTTrackHoleSearchTool::m_end_at_last_trt_hit

private

Definition at line 90 of file TRTTrackHoleSearchTool.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_extrapolator

ToolHandle<Trk::IExtrapolator> TRTTrackHoleSearchTool::m_extrapolator

private

Definition at line 83 of file TRTTrackHoleSearchTool.h.

◆ m_has_been_called

std::atomic_bool TRTTrackHoleSearchTool::m_has_been_called

mutableprivate

Definition at line 99 of file TRTTrackHoleSearchTool.h.

◆ m_locR_cut

float TRTTrackHoleSearchTool::m_locR_cut

private

Definition at line 93 of file TRTTrackHoleSearchTool.h.

◆ m_locR_sigma_cut

float TRTTrackHoleSearchTool::m_locR_sigma_cut

private

Definition at line 94 of file TRTTrackHoleSearchTool.h.

◆ m_locZ_cut

float TRTTrackHoleSearchTool::m_locZ_cut

private

Definition at line 95 of file TRTTrackHoleSearchTool.h.

◆ m_max_trailing_holes

int TRTTrackHoleSearchTool::m_max_trailing_holes

private

Definition at line 88 of file TRTTrackHoleSearchTool.h.

◆ m_max_z

float TRTTrackHoleSearchTool::m_max_z

private

Definition at line 87 of file TRTTrackHoleSearchTool.h.

◆ m_outer_radius

float TRTTrackHoleSearchTool::m_outer_radius

private

Definition at line 86 of file TRTTrackHoleSearchTool.h.

◆ m_TRT_ID

const TRT_ID* TRTTrackHoleSearchTool::m_TRT_ID

private

Definition at line 100 of file TRTTrackHoleSearchTool.h.

◆ m_trt_outer_surf

Trk::CylinderSurface* TRTTrackHoleSearchTool::m_trt_outer_surf

private

Definition at line 101 of file TRTTrackHoleSearchTool.h.

◆ m_use_conditions_svc

bool TRTTrackHoleSearchTool::m_use_conditions_svc

private

Definition at line 85 of file TRTTrackHoleSearchTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ TRTTrackHoleSearchTool()

Member Function Documentation

◆ addHolesToTrack()

◆ countHoles()

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

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ extrapolateBetweenHits()

◆ finalize()

◆ find_first_trt_hit()

◆ find_last_hit_before_trt()

◆ getHolesOnTrack()

◆ getTrackWithHoles()

◆ getTrackWithHolesAndOutliers()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_bcheck

◆ m_begin_at_first_trt_hit

◆ m_conditions_svc

◆ m_detStore

◆ m_do_dump_bad_straw_log

◆ m_end_at_last_trt_hit

◆ m_evtStore

◆ m_extrapolator

◆ m_has_been_called

◆ m_locR_cut

◆ m_locR_sigma_cut

◆ m_locZ_cut

◆ m_max_trailing_holes

◆ m_max_z

◆ m_outer_radius

◆ m_TRT_ID

◆ m_trt_outer_surf

◆ m_use_conditions_svc

◆ m_varHandleArraysDeclared

◆ m_vhka