InDet::SimpleTRT_SeededSpacePointFinder_ATL Class Reference

#include <SimpleTRT_SeededSpacePointFinder_ATL.h>

Inheritance diagram for InDet::SimpleTRT_SeededSpacePointFinder_ATL:

Collaboration diagram for InDet::SimpleTRT_SeededSpacePointFinder_ATL:

Public Types
typedef std::list< std::pair< int, int > >	modulLookupTable

Public Member Functions
	SimpleTRT_SeededSpacePointFinder_ATL (const std::string &, const std::string &, const IInterface *)
virtual	~SimpleTRT_SeededSpacePointFinder_ATL ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual std::unique_ptr< InDet::ITRT_SeededSpacePointFinder::IEventData >	newEvent () const override
	obsolete, do not do anything.
virtual std::unique_ptr< InDet::ITRT_SeededSpacePointFinder::IEventData >	newRegion (const std::vector< IdentifierHash > &, const std::vector< IdentifierHash > &) const override
const SiSpacePointsSeed *	next (InDet::ITRT_SeededSpacePointFinder::IEventData &event_data) const override
std::list< std::pair< const Trk::SpacePoint *, const Trk::SpacePoint * > >	find2Sp (const EventContext &ctx, const Trk::TrackParameters &, InDet::ITRT_SeededSpacePointFinder::IEventData &event_data) const override
	main method, calls the private methods and returns a pointer to the list of SpacePointpairs.
MsgStream &	dump (MsgStream &out) const override
std::ostream &	dump (std::ostream &out) const override
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () 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
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.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	getSpacePointsInROI (const EventContext &ctx, std::set< IdentifierHash > &setOfSCT_Hashes, int modulTRT, std::multimap< int, const Trk::SpacePoint * > &relevantSpacePoints) const
	retrieves SP Collections of modules in the ROI and sorts them by SCT layer
void	getHashesInROI (const EventContext &ctx, const Trk::TrackParameters &directionTRT, std::set< IdentifierHash > &setOfSCT_Hashes) const
	obtains the hashes of modules in the ROI
void	combineSpacePoints (const std::multimap< int, const Trk::SpacePoint * > &relevantSpacePoints, const Trk::TrackParameters &directionTRT, int modulTRT, std::list< std::pair< const Trk::SpacePoint *, const Trk::SpacePoint * > > &listOfSpacePointPairsBuffer) const
	builds pairs of SP according to the kook-up table
bool	pairIsOk (const Trk::SpacePoint *sp1, const Trk::SpacePoint *sp2, const Trk::TrackParameters &directionTRT) const
	applies rough cuts on the quality of a SP pair
void	setupLookUpTable ()
	the name says it
void	printLookupTable () const
int	TRT_Module (const Trk::TrackParameters &directionTRT) const
	returns the number of the TRT wheel or barrel where the TP belongs to
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Static Private Member Functions
static void	getSearchRange (double &deltaPhi, double &deltaEta)
	List with SP pairs as seed for the Si part of the back-track.

Private Attributes
BooleanProperty	m_useROI {this, "RestrictROI", true}
	Controls, if SP have to be checked with the AssociationTool of the forward tracking and to avoid double use of measurements.
IntegerProperty	m_maxHoles {this, "MaxHoles", 1}
	controls how many not considered SCT layers are allowed between two SP in order to form a seed pair
DoubleProperty	m_perigeeCut {this, "PerigeeCut", 200.}
	rough cuts on the quality of the suggested SP pair
DoubleProperty	m_directionPhiCut {this, "DirectionPhiCut", 0.05}
SG::ReadHandleKey< SpacePointContainer >	m_spacepointsSCTname {this,"SpacePointsSCTName","SCT_SpacePoints","RHK to retrieve SCT SpacePointContainer"}
SG::ReadHandleKey< SpacePointOverlapCollection >	m_spacepointsOverlapname {this,"SpacePointsSCTName","OverlapSpacePoints","RHK to retrieve OverlapCollection"}
SG::ReadHandleKey< Trk::PRDtoTrackMap >	m_prdToTrackMap {this,"PRDtoTrackMap",""}
modulLookupTable	m_modulLookupTable [2 *SIMPLE_TRT_INDEX_OFFSET+1]
	Lookup table that contains the SCT Layers to be considered to provide SP for the pairing in dependence on the TRT part the track parameter is on.
std::set< int >	m_modulLookupTableIndex [2 *SIMPLE_TRT_INDEX_OFFSET+1]
ToolHandle< IRegSelTool >	m_pRegionSelector { this, "RegSelTool_SCT", "RegSelTool/RegSelTool_SCT" }
	Region Selector.
const SCT_ID *	m_sctId = nullptr
	ID SCT helper.
const TRT_ID *	m_trtId = nullptr
	ID TRT helper.
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 38 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

Member Typedef Documentation

modulLookupTable

typedef std::list<std::pair< int,int > > InDet::SimpleTRT_SeededSpacePointFinder_ATL::modulLookupTable

Definition at line 47 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

SimpleTRT_SeededSpacePointFinder_ATL()

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

Definition at line 31 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

  : AthAlgTool(t,n,p)
{
  declareInterface<ITRT_SeededSpacePointFinder>(this);
}

~SimpleTRT_SeededSpacePointFinder_ATL()

InDet::SimpleTRT_SeededSpacePointFinder_ATL::~SimpleTRT_SeededSpacePointFinder_ATL ( )

virtualdefault

Member Function Documentation

combineSpacePoints()

void InDet::SimpleTRT_SeededSpacePointFinder_ATL::combineSpacePoints ( const std::multimap< int, const Trk::SpacePoint * > & relevantSpacePoints, const Trk::TrackParameters & directionTRT, int modulTRT, std::list< std::pair< const Trk::SpacePoint *, const Trk::SpacePoint * > > & listOfSpacePointPairsBuffer ) const

private

builds pairs of SP according to the kook-up table

This method combines SpacePoints contained in the first argument to pairs which are then stored in the internal list of SpacePoint pairs.

Only SpacePoints which are on two different SCT layers or wheels specified as pair in the look-up table are combined. The pairing is only accepted, if if satisfies quality cuts specified in the method pairIsOk.

The second argument specifies the TRT track parameter which is used for quality cuts on the SpacePoint pairs. The third argument specifies the part of the TRT the track parameter comes from. This number is used to inquire the relevant pairings of SCT layers from the look-up table.

Definition at line 417 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
  
  if ( modulNumber+SIMPLE_TRT_INDEX_OFFSET < 29 && -1 < modulNumber+SIMPLE_TRT_INDEX_OFFSET)
    {
      modulLookupTable::const_iterator  itTable = m_modulLookupTable[modulNumber+SIMPLE_TRT_INDEX_OFFSET].begin();
      modulLookupTable::const_iterator endTable = m_modulLookupTable[modulNumber+SIMPLE_TRT_INDEX_OFFSET].end();
      for ( ; itTable != endTable ; ++itTable)
    {
      msg(MSG::VERBOSE) << " Combining Space Poinst from modules " << itTable->first << " and " << itTable->second << endmsg; 
      // SPs from SCT layer 1
      std::pair< std::multimap<int,const Trk::SpacePoint*>::const_iterator,
        std::multimap<int,const Trk::SpacePoint*>::const_iterator  >  
        range1 = relevantSpacePoints.equal_range(itTable->first);
      //      if (range1.second != relevantSpacePoints.end() ) ++(range1.second);  // transforming the 2nd pointer to an end of range
 
      //SPs from SCT layer 2
      std::pair< std::multimap<int,const Trk::SpacePoint*>::const_iterator,
        std::multimap<int,const Trk::SpacePoint*>::const_iterator  >  
        range2 = relevantSpacePoints.equal_range(itTable->second);
      //      if (range2.second != relevantSpacePoints.end() ) ++(range2.second);  // transforming the 2nd pointer to an end of range
 
      // add the SP pairs
      for ( std::multimap<int,const Trk::SpacePoint*>::const_iterator it1 = range1.first; 
        it1 != range1.second ; ++it1 )        
        for ( std::multimap<int,const Trk::SpacePoint*>::const_iterator it2 = range2.first; 
          it2 != range2.second ; ++it2 )          
          if ( pairIsOk(it1->second,it2->second,directionTRT) )
        listOfSpacePointPairsBuffer.emplace_back(it2->second,it1->second);
 
    }
    }
  else
    msg(MSG::WARNING) << "TRT module not in look-up table --> no SP pairs formed" << endmsg;
 
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & 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());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & 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.

{ return m_detStore; }

dump() [1/2]

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

overridevirtual

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 758 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
  out << "to be implemented soon..." << std::endl;
 
  return out;
}

dump() [2/2]

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

overridevirtual

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 766 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
  out << "to be implemented soon..." << std::endl;
  return out;
}

evtStore()

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.

{ return m_evtStore; }

extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

finalize()

StatusCode InDet::SimpleTRT_SeededSpacePointFinder_ATL::finalize ( )

overridevirtual

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 91 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
   StatusCode sc = AthAlgTool::finalize(); 
   return sc;
}

find2Sp()

std::list< std::pair< const Trk::SpacePoint *, const Trk::SpacePoint * > > InDet::SimpleTRT_SeededSpacePointFinder_ATL::find2Sp ( const EventContext & ctx, const Trk::TrackParameters & directionTRT, InDet::ITRT_SeededSpacePointFinder::IEventData & event_data ) const

overridevirtual

main method, calls the private methods and returns a pointer to the list of SpacePointpairs.

The list itself is a private member of this class

The main method. This method is called from outside to form pairs of SCT SpacePoints which could seed Si tracks linked to the TRT TrackParameter given as input. iModus is a dummy variable to satisfy the interface.

This method calls several private methods to perform specific tasks. These are:

getHashesInROI : If the restriction of the Region Of Intereset is used, this method returns a set of IdHashes of SCT modules within the ROI

getSpacePointsInROI : This Method retrieves the SpacePointCollections from StoreGate and fills a map with the SpacePoints. The map key is a number which identifies the layer or disc the SpacePoint is on.

combineSpacePoints : The retrieved and sorted SpacePoints are combined to pairs. Pairs are built only of SpacePoints which are on the layer or disc pairs specified in the internal look-up table.

The method returns a pointer to an internal list of SpacePoint pairs, the list may very well be empty nevertheless.

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 109 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{               
 
 
  msg(MSG::VERBOSE) << "Enter getListOfSpacePointPairs, TrackParameter given is: " << endmsg;
  msg(MSG::VERBOSE) << directionTRT << endmsg;
 
  // clear output buffer
 
  std::list<std::pair<const Trk::SpacePoint*, const Trk::SpacePoint*> > listOfSpacePointPairsBuffer;
 
  // IdHash in ROI
  std::set<IdentifierHash> setOfSCT_Hashes;
 
  if (m_useROI)
    {
      // fill IdHashes in ROI
      getHashesInROI(ctx,directionTRT,setOfSCT_Hashes);
      msg(MSG::VERBOSE) << "Retrieved " << setOfSCT_Hashes.size() << " potentially interesting detector elements." << endmsg;
    }
 
  if ( !m_useROI || !setOfSCT_Hashes.empty())
    {
      // map of < SCT layer number, SP* >
      std::multimap<int,const Trk::SpacePoint*> relevantSpacePoints; 
 
      // fill the map of relevant SP as defined by hashes from ROI
      int modulTRT = TRT_Module(directionTRT);
      getSpacePointsInROI(ctx, setOfSCT_Hashes, modulTRT, relevantSpacePoints);
 
      msg(MSG::VERBOSE) << "Retrieved " << relevantSpacePoints.size() << " potentially interesting SpacePoints" << endmsg;
 
      // build pairs of the relevant SP according to the look-up table
      combineSpacePoints(relevantSpacePoints, directionTRT, modulTRT, listOfSpacePointPairsBuffer);
 
      /* output for debug purposes, deactivated now. Once development is finished, it will be removed.
       */
      Amg::Vector3D r0 = directionTRT.position();
      const Amg::Vector3D& v0(directionTRT.momentum());
      msg(MSG::VERBOSE) << "------------------------------------------------------------------------------------------" << endmsg;
      msg(MSG::VERBOSE) << "Final SpacePoint pairs: " << listOfSpacePointPairsBuffer.size() << endmsg;
      msg(MSG::VERBOSE) << "   Position of initial vector:  ( " << r0.x() << ", " << r0.y() << ", "<< r0.z() << " ) " << endmsg;
      msg(MSG::VERBOSE) << "   Direction of initial vector:  ( " << v0.unit().x() << ", " << v0.unit().y() << ", "<< v0.unit().z() << " ) , phi = " 
        << v0.phi() <<  "   theta = " << v0.theta() << "    eta = "<< v0.eta() << endmsg;
      msg(MSG::VERBOSE) << "------------------------------------------------------------------------------------------" << endmsg;
      msg(MSG::VERBOSE) << "   Direction of space Point vectors: "<< endmsg;
      msg(MSG::VERBOSE) << ". . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . " << endmsg;
      for (auto & it : listOfSpacePointPairsBuffer)
    {
      Amg::Vector3D s1 = it.first->globalPosition();
      Amg::Vector3D s2 = it.second->globalPosition();
      Amg::Vector3D s1s2 = s2-s1;  // vector from s1 to s2
      msg(MSG::VERBOSE) << "   Positions:  ( " << s1.x() << " , "<< s1.y() << " , "<< s1.z() << " ) " << endmsg;
      msg(MSG::VERBOSE) << "               ( " << s2.x() << " , "<< s2.y() << " , "<< s2.z() << " ) " << endmsg;
      msg(MSG::VERBOSE) << "   direction:  ( " 
        << s1s2.unit().x() << ", " << s1s2.unit().y() << ", " << s1s2.unit().z() <<" ) , phi = "
        << s1s2.phi() <<  "   theta = " << s1s2.theta() << "    eta = "<< s1s2.eta() << endmsg;
      msg(MSG::VERBOSE) << ". . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . " << endmsg;
    }
      msg(MSG::VERBOSE) << "------------------------------------------------------------------------------------------" << endmsg;
    }
 
  return listOfSpacePointPairsBuffer;
}

getHashesInROI()

void InDet::SimpleTRT_SeededSpacePointFinder_ATL::getHashesInROI ( const EventContext & ctx, const Trk::TrackParameters & directionTRT, std::set< IdentifierHash > & setOfSCT_Hashes ) const

private

obtains the hashes of modules in the ROI

This Method finds the IdHashes of SCT modules inside a Region Of Interest around the TRT track parameter. First the allowed tolerance in phi and eta is inquired, then a vector with the relevant hashes is retrieved from the RegionSelectorService. In order to be able to search conveniently for hashes later on, they are then copied into a STL set, which is then returned via call by reference.

Definition at line 196 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
 
  //double eta = directionTRT.eta();
  //double phi = directionTRT.position().phi();
  double phi = directionTRT.parameters()[Trk::phi];
  double eta = directionTRT.position().eta();
  double deltaPhi = 0.;
  double deltaEta = 0.;
 
  // define tolerance for ROI
  //getSearchRange(directionTRT, deltaPhi, deltaEta);
  getSearchRange(deltaPhi, deltaEta);
 
  // retrieve SCT hashes in Region of interest  
  std::vector<IdentifierHash> listOfSCT_Hashes;

  RoiDescriptor roi(  eta-deltaEta, eta+deltaEta, phi-deltaPhi, phi+deltaPhi);
 
  m_pRegionSelector->lookup(ctx)->HashIDList( roi, listOfSCT_Hashes );
 
  // copy Hashes into Set to be able to search them
  for (auto listOfSCT_Hashe : listOfSCT_Hashes)
    setOfSCT_Hashes.insert(listOfSCT_Hashe);
}

getSearchRange()

void InDet::SimpleTRT_SeededSpacePointFinder_ATL::getSearchRange ( double & deltaPhi, double & deltaEta )

staticprivate

List with SP pairs as seed for the Si part of the back-track.

It is cleared and filled in every call of find2Sp, the pointer to it is returned. determines the width of the considered ROI around the TP. deltaPhi and deltaEta are return values

This method defines the tolerance for the ROI around the TRT track parameter. At the moment it is set constant, but it should be considered to define it in dependence on the track parameter. In order to suppress warnings for unused variables, the method signature including the track parameter as an argument is replaced for the time being by a signature without the track parameter.

Definition at line 403 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
  // first guess: the ROI tolerance is constant
  deltaPhi = 0.05;
  deltaEta = .2;
}

getSpacePointsInROI()

void InDet::SimpleTRT_SeededSpacePointFinder_ATL::getSpacePointsInROI ( const EventContext & ctx, std::set< IdentifierHash > & setOfSCT_Hashes, int modulTRT, std::multimap< int, const Trk::SpacePoint * > & relevantSpacePoints ) const

private

retrieves SP Collections of modules in the ROI and sorts them by SCT layer

This method retrieves the SpacePoints in a given region of interest. The region of interest is defined by the hashes in the first argument of the method.

The second argument specifies the part of the TRT in which the TRT track parameter is situated. This number is used to inquire the index of the look-up table, which layers of the SCT are relevant for the pair building process. Only SpacePointCollections from relevant layers or discs are looped over in order to save time.

If the restriction of the ROI is used, it is inquired for every SpacePointCollection, if the corresponding SCT module is listed in the set of IdHashes which defines the ROI. In this case, only Collections which satisfy this are looped over.

Relevant SpacePoints are then stored in a map. Key of the map is a number which identifies the SCT layer the SpacePoints are on.

In a last step, the SpacePointOverlapCollection is retrieved and looped over. SpacePoints in the ROI (if use is specified) are added to the list of SpacePoints

If the AssociationTool is used, only SpacePoints are stored, which do not contain a PRD that is not yet registered in the AssociationTool, i.e. which is not already used by the inside-out tracking.

Definition at line 230 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
 
  const std::set<IdentifierHash>::const_iterator endSCT_Hashes = setOfSCT_Hashes.end();
  
  SG::ReadHandle<Trk::PRDtoTrackMap>  prd_to_track_map;
  const Trk::PRDtoTrackMap *prd_to_track_map_cptr = nullptr;
  if (!m_prdToTrackMap.key().empty()) {
    prd_to_track_map=SG::ReadHandle<Trk::PRDtoTrackMap>(m_prdToTrackMap, ctx);
    if (!prd_to_track_map.isValid()) {
      ATH_MSG_ERROR("Failed to read PRD to track association map: " << m_prdToTrackMap.key());
    }
    prd_to_track_map_cptr = prd_to_track_map.cptr();
  }
 
  // retrieve SP Container
  SG::ReadHandle<SpacePointContainer> spacepointsSCT(m_spacepointsSCTname, ctx);
  if(spacepointsSCT.isValid()) 
    {
      // loop over SP collections in SP container
      SpacePointContainer::const_iterator itCont  =  spacepointsSCT->begin();
      SpacePointContainer::const_iterator endCont =  spacepointsSCT->end  ();
      for(; itCont != endCont; ++itCont) 
    {
      bool acceptCollection = true;
 
      if (m_useROI)
        {
          // Check if IdHash of the Collection is in the list of ROI hashes
          IdentifierHash idHash = (*itCont)->identifyHash();
          acceptCollection = ( idHash.is_valid() && setOfSCT_Hashes.find( idHash ) != endSCT_Hashes );
        }
 
      int SCT_LayerNumber = -99;
      if (acceptCollection)
        {
          // check if the module is on a relevant layer
          Identifier id  = (*itCont)->identify();
          int detRegionIndex = m_sctId->barrel_ec(id) ;
          if ( detRegionIndex == -2 ||  detRegionIndex == 0 ||  detRegionIndex == 2 )
        {
          // conversion factors for layer numbering
          int detLayerOffset = 1;         
          int detLayerSign   = 1;
          if (detRegionIndex == 0) detLayerOffset = 10;
          if (detRegionIndex == -2) detLayerSign = -1;
          
          // retrieve layer number from IdHelper
          int detLayerIndex = m_sctId->layer_disk(id);
          SCT_LayerNumber =  detLayerSign*(detLayerIndex+detLayerOffset);
        }
          else
        {
          msg(MSG::VERBOSE) << "detRegionIndex " << detRegionIndex << " recieved. " << endmsg;
          acceptCollection = false;
        }
        }
 
      if (acceptCollection)
        {
          if (modulTRT > -14 && modulTRT < 15 )
        {
          std::set<int>::const_iterator pos = m_modulLookupTableIndex[modulTRT+SIMPLE_TRT_INDEX_OFFSET].find(SCT_LayerNumber);
          if ( pos == m_modulLookupTableIndex[modulTRT+SIMPLE_TRT_INDEX_OFFSET].end() ) acceptCollection = false;
        }
          else
        {
          msg(MSG::WARNING) << "Received TRTmodul number " << modulTRT << endmsg;
          acceptCollection = false;
        }
        }
 
      if (acceptCollection)
        {         
          // Loop over SP Collection and add SP if they are not yet used (or if this does not matter). 
          SpacePointCollection::const_iterator itColl  = (*itCont)->begin();
          SpacePointCollection::const_iterator endColl = (*itCont)->end  ();
          for(; itColl != endColl; ++itColl) 
        if (    !prd_to_track_map_cptr
            || !(    prd_to_track_map_cptr->isUsed(*((*itColl)->clusterList().first))
                 || prd_to_track_map_cptr->isUsed(*((*itColl)->clusterList().second)) ) )
          {
            relevantSpacePoints.insert(std::make_pair( SCT_LayerNumber   ,*itColl));
            msg(MSG::VERBOSE) << "Added SpacePoint for layer " << SCT_LayerNumber << " at  ( " 
              << (*itColl)->globalPosition().x() << " , " 
              << (*itColl)->globalPosition().y() << " , " 
              << (*itColl)->globalPosition().z() << " ) " << endmsg; 
          }
        }
    }
    }
 
  // retrieve the overlap collection
  SG::ReadHandle<SpacePointOverlapCollection> spacepointsOverlap(m_spacepointsOverlapname, ctx);
  if(spacepointsOverlap.isValid()) 
    {
 
      // Loop over Overlap SP
      SpacePointOverlapCollection::const_iterator  itColl  = spacepointsOverlap->begin();
      SpacePointOverlapCollection::const_iterator endColl = spacepointsOverlap->end  ();
      for (; itColl != endColl; ++itColl) 
    {
      // check if SP is in ROI
      std::pair<IdentifierHash, IdentifierHash> idHashPair = (*itColl)->elementIdList();
      if (    !m_useROI  
           || (idHashPair.first.is_valid()  && setOfSCT_Hashes.find( idHashPair.first )  != endSCT_Hashes)
           || (idHashPair.second.is_valid() && setOfSCT_Hashes.find( idHashPair.second ) != endSCT_Hashes) )
      {
 
        // find out if one of the Clusters has already been used, if relevant
        if(prd_to_track_map_cptr)
          {
        bool u1=false; 
        bool u2=false;
        const Trk::PrepRawData* p1=(*itColl)->clusterList().first; 
        u1=prd_to_track_map->isUsed(*p1);
        const Trk::PrepRawData* p2=(*itColl)->clusterList().second;
        u2=prd_to_track_map->isUsed(*p2);
        if(u1 || u2) continue;
          }
        
        // retrieve identifier and fill SP in SP buffer.
        Identifier id = (*itColl)->associatedSurface().associatedDetectorElementIdentifier();
        if ( id.is_valid() )
          {
        int detRegionIndex = m_sctId->barrel_ec(id) ;          
        if ( detRegionIndex == -2 ||  detRegionIndex == 0 ||  detRegionIndex == 2 )
          {
            int detLayerOffset = 1;       
            int detLayerSign   = 1;
            if (detRegionIndex == 0) detLayerOffset = 10;
            if (detRegionIndex == -2) detLayerSign = -1;
            int detLayerIndex = m_sctId->layer_disk(id);
            
            relevantSpacePoints.insert(std::make_pair(detLayerSign*(detLayerIndex+detLayerOffset),*itColl));
            msg(MSG::VERBOSE) << "Added OverlapSpacePoint for layer " << detLayerSign*(detLayerIndex+detLayerOffset) << " at  ( " 
              << (*itColl)->globalPosition().x() << " , " 
              << (*itColl)->globalPosition().y() << " , " 
              << (*itColl)->globalPosition().z() << " ) " << endmsg; 
          }
        else
          msg(MSG::VERBOSE) << "detRegionIndex " << detRegionIndex << " recieved. " << endmsg;
          }
        else
          msg(MSG::DEBUG) << "Invalid Id from OverlapCollection recieved. " << endmsg;
      }
    }
    }
}

initialize()

StatusCode InDet::SimpleTRT_SeededSpacePointFinder_ATL::initialize ( )

overridevirtual

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 49 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
  // msg().setLevel(outputLevel());
 
  // PRD-to-track association (optional)
  ATH_CHECK( m_prdToTrackMap.initialize( !m_prdToTrackMap.key().empty()));
 
  // get region selector
  StatusCode sc = m_pRegionSelector.retrieve();
  if( sc.isFailure() ) 
    {
      msg(MSG::FATAL) << "Failed to retrieve RegionSelector Service";
      return sc;
    }
 
  sc = detStore()->retrieve(m_sctId, "SCT_ID");
  if (sc.isFailure())
    {
      msg(MSG::ERROR) << "Could not get SCT_ID helper !" << endmsg;
      return StatusCode::FAILURE;
    }
 
  sc = detStore()->retrieve(m_trtId, "TRT_ID");
  if (sc.isFailure())
    {
      msg(MSG::ERROR) << "Could not get TRT_ID helper !" << endmsg;
      return StatusCode::FAILURE;
    }
 
  setupLookUpTable();
 
  ATH_CHECK(m_spacepointsSCTname.initialize());
  ATH_CHECK(m_spacepointsOverlapname.initialize());
 
  return sc;
}

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 & InDet::ITRT_SeededSpacePointFinder::interfaceID ( )

inlinestaticinherited

Definition at line 108 of file ITRT_SeededSpacePointFinder.h.

    {
      return IID_ITRT_SeededSpacePointFinder;
    }

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

newEvent()

std::unique_ptr< InDet::ITRT_SeededSpacePointFinder::IEventData > InDet::SimpleTRT_SeededSpacePointFinder_ATL::newEvent ( ) const

inlineoverridevirtual

obsolete, do not do anything.

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 178 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

                                                                                                                         {
     return std::unique_ptr<InDet::ITRT_SeededSpacePointFinder::IEventData>();
   }

newRegion()

std::unique_ptr< InDet::ITRT_SeededSpacePointFinder::IEventData > InDet::SimpleTRT_SeededSpacePointFinder_ATL::newRegion ( const std::vector< IdentifierHash > & , const std::vector< IdentifierHash > &  ) const

overridevirtual

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 102 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
   return {};
}

next()

const SiSpacePointsSeed * InDet::SimpleTRT_SeededSpacePointFinder_ATL::next ( InDet::ITRT_SeededSpacePointFinder::IEventData & event_data ) const

inlineoverridevirtual

Implements InDet::ITRT_SeededSpacePointFinder.

Definition at line 172 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

    {
      return nullptr;
    }

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.

pairIsOk()

bool InDet::SimpleTRT_SeededSpacePointFinder_ATL::pairIsOk ( const Trk::SpacePoint * sp1, const Trk::SpacePoint * sp2, const Trk::TrackParameters & directionTRT ) const

private

applies rough cuts on the quality of a SP pair

This method applies quality cuts on a given pairing of SpacePoints. The quality cuts at the moment are:

• the closest approach of the line connecting the two SpacePoints to the beam axis in (x,y) has to be smaller than a given limit.
• the direction in phi of the line connecting the two SpacePoints must not deviate from the phi direction of the TRT track parameter by more than a given tolerance.

If the cuts are satisfied, true is returned, false otherwise.

Definition at line 468 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
 
  Amg::Vector3D s1 = sp1->globalPosition();
  Amg::Vector3D s2 = sp2->globalPosition();
  //Amg::Vector3D r0 = directionTRT.position();
  const Amg::Vector3D& v0(directionTRT.momentum());
  Amg::Vector3D s1s2 = s2-s1;  // vector from s1 to s2
  
  msg(MSG::VERBOSE) << "Checking Space Point Pair at ( "<< s1.x() << ", " << s1.y() << ", " << s1.z()
    << " )  and  ( "<< s2.x() << ", " << s2.y() << ", " << s2.z() << " )" << endmsg;
 
  // Cut on closest approach to z-axis
  double t = ( s1.x()-s2.x() )*( s1.x()-s2.x() ) + ( s1.y()-s2.y() )*( s1.y()-s2.y() ) ;
  if (t)
    {
      t =  ( s1.x()*(s2.x()-s1.x()) + s1.y()*(s2.y()-s1.y()) ) / t;
      Amg::Vector3D perigee = s1 + t*(s1-s2);
      msg(MSG::VERBOSE) << " closest approach to beam-pipe at ( "<< perigee.x() << ", " << perigee.y() << ", " << perigee.z() << " )    transversal: " << perigee.perp()<<endmsg;
      if (perigee.perp() > m_perigeeCut ) return false;
    }
  msg(MSG::VERBOSE) << " Passed cut on r-phi impact parameter" <<endmsg;
 
 
  // Cut on angle between s1-s2 and parameter direction
  double diffPhi= std::abs(s1s2.phi()-v0.phi());
 
  if (diffPhi > CLHEP::pi ) diffPhi = 2.*CLHEP::pi - diffPhi;
  if (diffPhi > CLHEP::pi/2. ) diffPhi = CLHEP::pi - diffPhi;
  msg(MSG::VERBOSE) << "Phi directions differ by "<< diffPhi << endmsg;
  if ( diffPhi > m_directionPhiCut ) return false;
 
  msg(MSG::VERBOSE) << " Passed cut on direction phi deviation" <<endmsg;
 
  return true;
}

printLookupTable()

void InDet::SimpleTRT_SeededSpacePointFinder_ATL::printLookupTable ( ) const

private

This method provides a printout of the look-up table for debugging purposes. It is called in the setup.

Definition at line 775 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
 
  msg(MSG::VERBOSE) << "=====================================================================================================================" << endmsg;
  msg(MSG::VERBOSE) << "  Module Lookup table: " << endmsg;
  msg(MSG::VERBOSE) << "---------------------------------------------------------------------------------------------------------------------" << endmsg;
  for (int i = -14; i<15 ; ++i)
    {
      msg(MSG::VERBOSE) << "Module "<< i << "  : " << endmsg;
      std::list<std::pair<int,int> >::const_iterator  itList = m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].begin();
      std::list<std::pair<int,int> >::const_iterator endList = m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].end();
      for ( ; itList != endList ; ++itList )
    msg(MSG::VERBOSE) << " (" << itList->first << ", " << itList->second << ") " ; 
      msg(MSG::VERBOSE) << "  " << endmsg;
      std::set<int>::const_iterator  itIndex = m_modulLookupTableIndex[i+SIMPLE_TRT_INDEX_OFFSET].begin();
      std::set<int>::const_iterator endIndex = m_modulLookupTableIndex[i+SIMPLE_TRT_INDEX_OFFSET].end();
      msg(MSG::VERBOSE) << "Needs Layers: { " ;
      for ( ; itIndex != endIndex ; ++ itIndex )
    msg(MSG::VERBOSE) << *itIndex << ", " ;
      msg(MSG::VERBOSE) << " } " << endmsg;
      msg(MSG::VERBOSE) << "...................................................................................................................." << endmsg;
 
    }
  msg(MSG::VERBOSE) << "=====================================================================================================================" << endmsg;
 
}

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

setupLookUpTable()

void InDet::SimpleTRT_SeededSpacePointFinder_ATL::setupLookUpTable ( )

private

the name says it

This method sets up the look-up table. This table is an array of lists of integer pairs. The index of the array corresponds to a number which identifies the part of the TRT the TRT track parameter comes from, the two integer numbers in the pairs identify the SCT layers or wheels from which SpacePoints have to be paired if the TRT track parameter comes from the corresponding TRT part.

The numbering scheme for the TRT detectors: 0 -> barrel straws 1 - 14 -> forward wheels, counted from inside out -1 - -14 -> backward wheels, -1 the innermost

The numbering scheme for the SCT detector (first number ATLAS, second number within this Tool):

        ____backw_________________________     __barrel__    ________forw_____________

Index = 8 7 6 5 4 3 2 1 0 3 (13) 0 1 2 3 4 5 6 7 8 -9 -8 -7 -6 -5 -4 -3 -2 -1 2 (12) 1 2 3 4 5 6 7 8 9 1 (11) 0 (10)

The numbering scheme used here deviates from the numbering used in the ATLAS Detector description for two reasons:

• In order to be fast, the identification has to be done by one id number only. ATLAS id is based on two information: the barrel-endcap flag and the modul number. As the modul numbers alone are ambiguous (forward, backward, barrel may feature the same module numbers), an unambiguous id number is derived from both informations.

As anyway a specific numbering scheme has to be employed here, it is designed to be anti-symmetric with respect to forward and backward direction. This allows a convenient automatic set-up of the backward part of the look-up table, once the forward part is defined and thus reduces the likelihood of mistakes. In order to compensate for negative module numbers, an offset is added when they are used as index to the look-up table, because C++ does not allow negative indices (as FORTRAN would do ... :-) )

Depending on the number of allowed holes, i.e. basically the number of SCT layers or wheels omitted between the SpacePoints of a pair, additional pairings of SCT layers are added to the look-up table.

Definition at line 513 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
  m_modulLookupTable[0+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,13);
  if (m_maxHoles > 0 )
    {
      m_modulLookupTable[0+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,12);
      m_modulLookupTable[0+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,13);
      m_modulLookupTable[0+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,12);
    }
  if (m_maxHoles > 1 )
    {
      m_modulLookupTable[0+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,13);
      m_modulLookupTable[0+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,11);
    }
 
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,13);
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,12);
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,11);
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,1);
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,1);
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,1);
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,1);
  m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,2);
  if (m_maxHoles > 0 )
    {
      m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,13);
      m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,12);
      m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,2);
      m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,2);
      m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,2);
      m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,2);
    }
  if (m_maxHoles > 1 )
    {
      m_modulLookupTable[1+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,13);
    }
 
  //  TRT wheels 2 and 3
  for (int i = 2; i<4 ; ++ i)
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,13);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,3);
      if (m_maxHoles > 0 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,12);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,3);
    }
      if (m_maxHoles > 1 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,11);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,12);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,3);
    }
    }
 
  // TRT wheels 4 - 6
  for (int i = 4 ; i < 7 ; ++i )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,13);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,1);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,5);
      if (m_maxHoles > 0 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,13);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,12);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,5);
    }
      if (m_maxHoles > 1 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(11,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(10,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,5);
    }
    }
 
 
  // TRT wheels 7 - 9
  for ( int i = 7 ;  i < 10 ; ++i)
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(5,6);  
      if (m_maxHoles > 0 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,2);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,6);
    }
      if (m_maxHoles > 1 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(13,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(12,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(1,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,6);
    }
    }
 
 
  // TRT wheels 10 - 12
  for (int i = 10 ; i < 13 ; ++i )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(5,6);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(6,7);
      if (m_maxHoles > 0 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,6);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(5,7);
    }
      if (m_maxHoles > 1 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,3);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(2,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,6);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,7);
    }
    }
 
 
  // TRT - wheels 13 and 14
  for (int i = 13 ; i < 15 ; ++i )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(5,6);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(6,7);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(7,8);
      if (m_maxHoles > 0 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,6);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(5,7);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(6,8);
    }
      if (m_maxHoles > 0 )
    {
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,4);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,5);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(3,6);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(4,7);
      m_modulLookupTable[i+SIMPLE_TRT_INDEX_OFFSET].emplace_back(5,8);
    }
    }
 
 
  // mirror the endcap information
  for ( int iTable = 1 ; iTable < 15 ; ++iTable )
    {
      std::list<std::pair<int,int> >::const_iterator  itList = m_modulLookupTable[iTable+SIMPLE_TRT_INDEX_OFFSET].begin();
      std::list<std::pair<int,int> >::const_iterator endList = m_modulLookupTable[iTable+SIMPLE_TRT_INDEX_OFFSET].end();
      for ( ; itList != endList ; ++itList )
    {
      int i = itList->first;
      int j = itList->second;
      if (i<10) i*=-1;   // don't change the sign of the barrel modules...
      if (j<10) j*=-1;
      m_modulLookupTable[-1*iTable+SIMPLE_TRT_INDEX_OFFSET].emplace_back(i,j);
    }
    }     
 
  // fill the short index for the lookup table
  for ( int iTable = -14 ; iTable < 15 ; ++iTable )
    {
      std::list<std::pair<int,int> >::const_iterator  itList = m_modulLookupTable[iTable+SIMPLE_TRT_INDEX_OFFSET].begin();
      std::list<std::pair<int,int> >::const_iterator endList = m_modulLookupTable[iTable+SIMPLE_TRT_INDEX_OFFSET].end();
      for ( ; itList != endList ; ++itList )
    {
      m_modulLookupTableIndex[iTable+SIMPLE_TRT_INDEX_OFFSET].insert(itList->first);
      m_modulLookupTableIndex[iTable+SIMPLE_TRT_INDEX_OFFSET].insert(itList->second);
    }
    }
 
  printLookupTable();
 
}

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 asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

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.

TRT_Module()

int InDet::SimpleTRT_SeededSpacePointFinder_ATL::TRT_Module ( const Trk::TrackParameters & directionTRT ) const

private

returns the number of the TRT wheel or barrel where the TP belongs to

This method returns the id number of the TRT part, the given trackparameter comes from following the numbering scheme applied in this Tool (see documentation to setupLookupTable).

Definition at line 806 of file SimpleTRT_SeededSpacePointFinder_ATL.cxx.

{
 
  // find out from which SCT part the segment is built
  Identifier id  ;
 
    id = directionTRT.associatedSurface().associatedDetectorElementIdentifier();
 
  if (!id.is_valid())
    {
      msg(MSG::WARNING) << " Id not valid "<<endmsg;
      return -999;
    }
  int modulNumber;
  int detRegionIndex = m_trtId->barrel_ec(id) ;  
  if ( detRegionIndex == -1 || detRegionIndex == 1 )
    modulNumber = 0;  // central TRT
  else if ( detRegionIndex == -2 || detRegionIndex == 2 )
    {
      modulNumber = m_trtId->layer_or_wheel(id) + 1  ; // forward/bachward wheels
      modulNumber *= (detRegionIndex/2) ;  // add sign to distinguish forward and backward
    }
  else
    {
      msg(MSG::WARNING) << "TRT barrel-endcap id not in {-2;-1;1;2}" << endmsg;
      return -999;
    }
 
  return modulNumber;
 
}

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_directionPhiCut

DoubleProperty InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_directionPhiCut {this, "DirectionPhiCut", 0.05}

private

Definition at line 101 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{this, "DirectionPhiCut", 0.05};

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_maxHoles

IntegerProperty InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_maxHoles {this, "MaxHoles", 1}

private

controls how many not considered SCT layers are allowed between two SP in order to form a seed pair

Definition at line 97 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{this, "MaxHoles", 1};

m_modulLookupTable

modulLookupTable InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_modulLookupTable[2 *SIMPLE_TRT_INDEX_OFFSET+1]

private

Lookup table that contains the SCT Layers to be considered to provide SP for the pairing in dependence on the TRT part the track parameter is on.

Definition at line 116 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

m_modulLookupTableIndex

std::set<int> InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_modulLookupTableIndex[2 *SIMPLE_TRT_INDEX_OFFSET+1]

private

Definition at line 117 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

m_perigeeCut

DoubleProperty InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_perigeeCut {this, "PerigeeCut", 200.}

private

rough cuts on the quality of the suggested SP pair

Definition at line 100 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{this, "PerigeeCut", 200.};

m_prdToTrackMap

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

private

Definition at line 111 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{this,"PRDtoTrackMap",""};

m_pRegionSelector

ToolHandle<IRegSelTool> InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_pRegionSelector { this, "RegSelTool_SCT", "RegSelTool/RegSelTool_SCT" }

private

Region Selector.

Definition at line 121 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{ this, "RegSelTool_SCT", "RegSelTool/RegSelTool_SCT" };

m_sctId

const SCT_ID* InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_sctId = nullptr

private

ID SCT helper.

Definition at line 124 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

m_spacepointsOverlapname

SG::ReadHandleKey<SpacePointOverlapCollection> InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_spacepointsOverlapname {this,"SpacePointsSCTName","OverlapSpacePoints","RHK to retrieve OverlapCollection"}

private

Definition at line 110 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{this,"SpacePointsSCTName","OverlapSpacePoints","RHK to retrieve OverlapCollection"}        ;

m_spacepointsSCTname

SG::ReadHandleKey<SpacePointContainer> InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_spacepointsSCTname {this,"SpacePointsSCTName","SCT_SpacePoints","RHK to retrieve SCT SpacePointContainer"}

private

Definition at line 109 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{this,"SpacePointsSCTName","SCT_SpacePoints","RHK to retrieve SCT SpacePointContainer"}            ;

m_trtId

const TRT_ID* InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_trtId = nullptr

private

ID TRT helper.

Definition at line 127 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

m_useROI

BooleanProperty InDet::SimpleTRT_SeededSpacePointFinder_ATL::m_useROI {this, "RestrictROI", true}

private

Controls, if SP have to be checked with the AssociationTool of the forward tracking and to avoid double use of measurements.

It is set true, if an AssociationTool is specified in the setup controls if the ROI is restricted around the TRT trackparameters

Definition at line 90 of file SimpleTRT_SeededSpacePointFinder_ATL.h.

{this, "RestrictROI", true};

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:

• SimpleTRT_SeededSpacePointFinder_ATL.h
• SimpleTRT_SeededSpacePointFinder_ATL.cxx