InDet::SiCombinatorialTrackFinder_xk is algorithm which produce track-finding in the road of InDetDD::SiDetectorElement* of SCT and Pixels sorted in propagation order. More...

#include <SiCombinatorialTrackFinder_xk.h>

Inheritance diagram for InDet::SiCombinatorialTrackFinder_xk:

Collaboration diagram for InDet::SiCombinatorialTrackFinder_xk:

Public Member Functions
Standard tool methods
	SiCombinatorialTrackFinder_xk (const std::string &, const std::string &, const IInterface *)
virtual	~SiCombinatorialTrackFinder_xk ()=default
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
Main methods for local track finding
virtual const std::list< Trk::Track * > &	getTracks (SiCombinatorialTrackFinderData_xk &data, const Trk::TrackParameters &, const std::vector< const Trk::SpacePoint * > &, const std::vector< Amg::Vector3D > &, std::vector< const InDetDD::SiDetectorElement * > &, const TrackQualityCuts &, const EventContext &ctx) const override
virtual const std::list< Trk::Track * > &	getTracks (SiCombinatorialTrackFinderData_xk &data, const Trk::TrackParameters &, const std::vector< const Trk::SpacePoint * > &, const std::vector< Amg::Vector3D > &, std::vector< const InDetDD::SiDetectorElement * > &, std::multimap< const Trk::PrepRawData , const Trk::Track > &, const EventContext &ctx) const override
virtual const std::list< Trk::Track * > &	getTracksWithBrem (SiCombinatorialTrackFinderData_xk &data, const Trk::TrackParameters &, const std::vector< const Trk::SpacePoint * > &, const std::vector< Amg::Vector3D > &, std::vector< const InDetDD::SiDetectorElement * > &, std::multimap< const Trk::PrepRawData , const Trk::Track > &, bool, const EventContext &ctx) const override
virtual double	pTseed (SiCombinatorialTrackFinderData_xk &data, const Trk::TrackParameters &, const std::vector< const Trk::SpacePoint * > &, const EventContext &) const override
virtual void	newEvent (const EventContext &ctx, SiCombinatorialTrackFinderData_xk &data) const override
virtual void	newEvent (const EventContext &ctx, SiCombinatorialTrackFinderData_xk &data, Trk::TrackInfo, const TrackQualityCuts &) const override
virtual void	endEvent (SiCombinatorialTrackFinderData_xk &data) const override
Print internal tool parameters and status
MsgStream &	dump (SiCombinatorialTrackFinderData_xk &data, MsgStream &out) const override

Private Attributes
Tool handles
ToolHandle< IInDetConditionsTool >	m_pixelCondSummaryTool
ToolHandle< IInDetConditionsTool >	m_sctCondSummaryTool
ToolHandle< Trk::IBoundaryCheckTool >	m_boundaryCheckTool
PublicToolHandle< Trk::IPatternParametersPropagator >	m_proptool
PublicToolHandle< Trk::IPatternParametersUpdator >	m_updatortool
PublicToolHandle< Trk::IRIO_OnTrackCreator >	m_riocreator
Data handles
SG::ReadHandleKey< InDet::PixelClusterContainer >	m_pixcontainerkey
SG::ReadHandleKey< InDet::SCT_ClusterContainer >	m_sctcontainerkey
SG::ReadCondHandleKey< InDet::SiDetElementBoundaryLinks_xk >	m_boundaryPixelKey
SG::ReadCondHandleKey< InDet::SiDetElementBoundaryLinks_xk >	m_boundarySCTKey
SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"}
SG::ReadHandleKey< InDet::SiDetectorElementStatus >	m_pixelDetElStatus {this, "PixelDetElStatus", "", "Key of SiDetectorElementStatus for Pixel"}
	Optional read handle to get status data to test whether a pixel detector element is good.
SG::ReadHandleKey< InDet::SiDetectorElementStatus >	m_sctDetElStatus {this, "SCTDetElStatus", "", "Key of SiDetectorElementStatus for SCT"}
	Optional read handle to get status data to test whether a SCT detector element is good.
Properties
BooleanProperty	m_usePIX {this, "usePixel", true}
BooleanProperty	m_useSCT {this, "useSCT", true}
BooleanProperty	m_ITkGeometry {this, "ITkGeometry", false}
BooleanProperty	m_doFastTracking {this, "doFastTracking", false}
StringProperty	m_fieldmode {this, "MagneticFieldMode", "MapSolenoid", "Mode of magnetic field"}
DoubleProperty	m_qualityCut {this, "TrackQualityCut", 9.3, "Simple track quality cut"}
FloatProperty	m_minPtCut {this, "MinFinalPtCut", 100, "Cut on the pt of the final track. Must be >0 to avoid NANs when computing eta."}
float	m_minPt2Cut =0
BooleanProperty	m_writeHolesFromPattern {this, "writeHolesFromPattern", false,"Flag to activate writing hole info from the pattern recognition"}

Friends
class	SiTrajectory_xk

Data members, which are updated in only initialize
enum	EStat_t { TwoCluster , WrongRoad , WrongInit , CantFindTrk , NotNewTrk , BremAttempt , NumberOfStats , Success }
	Array entries for data.statistic counter. More...
int	m_outputlevel {0}
Trk::MagneticFieldProperties	m_fieldprop
	Magnetic field properties.
EStat_t	findTrack (SiCombinatorialTrackFinderData_xk &data, const Trk::TrackParameters &, const std::vector< const Trk::SpacePoint * > &, const std::vector< Amg::Vector3D > &, std::vector< const InDetDD::SiDetectorElement * > &, std::multimap< const Trk::PrepRawData , const Trk::Track > &, const EventContext &) const
Trk::Track *	convertToTrack (SiCombinatorialTrackFinderData_xk &data, const EventContext &ctx) const
Trk::Track *	convertToNextTrack (SiCombinatorialTrackFinderData_xk &data) const
void	magneticFieldInit ()
void	detectorElementLinks (std::vector< const InDetDD::SiDetectorElement * > &, std::vector< const InDet::SiDetElementBoundaryLink_xk * > &, const EventContext &ctx) const
MsgStream &	dumpconditions (MsgStream &out) const
void	initializeCombinatorialData (const EventContext &ctx, SiCombinatorialTrackFinderData_xk &data) const
virtual void	fillStatistic (SiCombinatorialTrackFinderData_xk &data, std::array< bool, NumberOfStats > &information) const override
static void	getTrackQualityCuts (SiCombinatorialTrackFinderData_xk &data, const TrackQualityCuts &)
static bool	spacePointsToClusters (const std::vector< const Trk::SpacePoint * > &, std::vector< const InDet::SiCluster * > &, std::optional< std::reference_wrapper< std::vector< const InDetDD::SiDetectorElement * > > >=std::nullopt)
static MsgStream &	dumpevent (SiCombinatorialTrackFinderData_xk &data, MsgStream &out)

Detailed Description

InDet::SiCombinatorialTrackFinder_xk is algorithm which produce track-finding in the road of InDetDD::SiDetectorElement* of SCT and Pixels sorted in propagation order.

In AthenaMT, event dependent cache inside SiCombinatorialTrackFinder_xk is not preferred. SiCombinatorialTrackFinderData_xk class holds event dependent data for SiCombinatorialTrackFinder_xk. Its object is instantiated in SiSPSeededTrackFinder::execute through SiTrackMakerEventData_xk.

Author: Igor..nosp@m.Gavr.nosp@m.ilenk.nosp@m.o@ce.nosp@m.rn.ch

Definition at line 60 of file SiCombinatorialTrackFinder_xk.h.

Member Enumeration Documentation

◆ EStat_t

enum InDet::SiCombinatorialTrackFinder_xk::EStat_t

private

Array entries for data.statistic counter.

Enumerator
TwoCluster
WrongRoad
WrongInit
CantFindTrk
NotNewTrk
BremAttempt
NumberOfStats
Success

Definition at line 206 of file SiCombinatorialTrackFinder_xk.h.

                   {
        TwoCluster,
        WrongRoad,
        WrongInit,
        CantFindTrk,
        NotNewTrk,
        BremAttempt,
        NumberOfStats,
        Success,
      } EStat_t;

Constructor & Destructor Documentation

◆ SiCombinatorialTrackFinder_xk()

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

Definition at line 34 of file SiCombinatorialTrackFinder_xk.cxx.

  : base_class(t, n, p)
{
}

◆ ~SiCombinatorialTrackFinder_xk()

virtual InDet::SiCombinatorialTrackFinder_xk::~SiCombinatorialTrackFinder_xk ( )

virtualdefault

Member Function Documentation

◆ convertToNextTrack()

Trk::Track * InDet::SiCombinatorialTrackFinder_xk::convertToNextTrack ( SiCombinatorialTrackFinderData_xk & data ) const

private

Definition at line 885 of file SiCombinatorialTrackFinder_xk.cxx.

{
  auto tsos = std::make_unique<Trk::TrackStates>(
      data.trajectory().convertToNextTrackStateOnSurface());
  if (tsos->empty()){
     return nullptr;
  }
 
  // verify first track parameters
  const Trk::TrackParameters *param = nullptr;
  for (const Trk::TrackStateOnSurface *a_tsos : *tsos) {
     const Trk::TrackParameters *param = a_tsos->trackParameters();
     if (param) {
        break;
     }
  }
 
  if (param) {
     auto momentum = param->momentum();
     const auto  pt2 = momentum.perp2();
     // reject tracks with small pT
     // The cut should be large enough otherwise eta computation of such tracks may yield NANs.
     if (pt2 < m_minPt2Cut) {
        ATH_MSG_WARNING( "Reject low pT track (pT = " << sqrt(pt2) << " < " << m_minPtCut.value() << ")");
        return nullptr;
     }
  }
  return new Trk::Track(data.trackinfo(),
                        std::move(tsos),
                        data.trajectory().convertToFitQuality());
}

◆ convertToTrack()

Trk::Track * InDet::SiCombinatorialTrackFinder_xk::convertToTrack	(	SiCombinatorialTrackFinderData_xk &	data,
		const EventContext &	ctx ) const

private

Definition at line 839 of file SiCombinatorialTrackFinder_xk.cxx.

{
  const Trk::PatternTrackParameters *param = data.trajectory().firstParameters();
  if (param) {
     double pt = param->transverseMomentum();
     // reject tracks with small pT
     // The cut should be large enough otherwise eta computation of such tracks may yield NANs.
     if (pt < m_minPtCut.value()) {
        ATH_MSG_DEBUG( "Reject low pT track (pT = " << pt << " < " << m_minPtCut.value() << ")");
        return nullptr;
     }
  }
  if (!data.simpleTrack()) {
     return new Trk::Track(
         data.trackinfo(),
         std::make_unique<Trk::TrackStates>(
             data.trajectory().convertToTrackStateOnSurface(
                 data.cosmicTrack())),
         data.trajectory().convertToFitQuality());
  }
 
  Trk::TrackInfo info = data.trackinfo();
  info.setPatternRecognitionInfo(Trk::TrackInfo::SiSPSeededFinderSimple);
  info.setParticleHypothesis(Trk::pion);
  if (!data.flagToReturnFailedTrack()) {
     return new Trk::Track(
         info,
         std::make_unique<Trk::TrackStates>(
             data.trajectory().convertToSimpleTrackStateOnSurface(
                 data.cosmicTrack(), ctx)),
         data.trajectory().convertToFitQuality());
  } else {
     return new Trk::Track(
         info,
         std::make_unique<Trk::TrackStates>(
             data.trajectory()
                 .convertToSimpleTrackStateOnSurfaceForDisTrackTrigger(
                     data.cosmicTrack(), ctx)),
         data.trajectory().convertToFitQuality());
  }
}

◆ detectorElementLinks()

void InDet::SiCombinatorialTrackFinder_xk::detectorElementLinks	(	std::vector< const InDetDD::SiDetectorElement * > &	DE,
		std::vector< const InDet::SiDetElementBoundaryLink_xk * > &	DEL,
		const EventContext &	ctx ) const

private

Definition at line 989 of file SiCombinatorialTrackFinder_xk.cxx.

{
  const InDet::SiDetElementBoundaryLinks_xk* boundaryPixel{nullptr};
  const InDet::SiDetElementBoundaryLinks_xk* boundarySCT{nullptr};
  if (m_usePIX) {
    SG::ReadCondHandle<InDet::SiDetElementBoundaryLinks_xk> boundaryPixelHandle(m_boundaryPixelKey,ctx);
    boundaryPixel = *boundaryPixelHandle;
    if (boundaryPixel==nullptr) {
      ATH_MSG_FATAL(m_boundaryPixelKey.fullKey() << " returns null pointer");
    }
  }
  if (m_useSCT) {
    SG::ReadCondHandle<InDet::SiDetElementBoundaryLinks_xk> boundarySCTHandle(m_boundarySCTKey,ctx);
    boundarySCT = *boundarySCTHandle;
    if (boundarySCT==nullptr) {
      ATH_MSG_FATAL(m_boundarySCTKey.fullKey() << " returns null pointer");
    }
  }
 
  DEL.reserve(DE.size());
  for (const InDetDD::SiDetectorElement* d: DE) {
    IdentifierHash id = d->identifyHash();
    if (d->isPixel() && boundaryPixel && id < boundaryPixel->size()) DEL.push_back(&(*boundaryPixel)[id]);
    else if (d->isSCT() && boundarySCT && id < boundarySCT->size()) DEL.push_back(&(*boundarySCT)[id]);
  }
}

◆ dump()

MsgStream & InDet::SiCombinatorialTrackFinder_xk::dump	(	SiCombinatorialTrackFinderData_xk &	data,
		MsgStream &	out ) const

override

Definition at line 119 of file SiCombinatorialTrackFinder_xk.cxx.

{
  if (not data.isInitialized()) initializeCombinatorialData(Gaudi::Hive::currentContext(), data);
 
  out<<std::endl;
  if (data.nprint()) return dumpevent(data, out);
  return dumpconditions(out);
}

◆ dumpconditions()

MsgStream & InDet::SiCombinatorialTrackFinder_xk::dumpconditions ( MsgStream & out ) const

private

Definition at line 132 of file SiCombinatorialTrackFinder_xk.cxx.

{
  int n = 62-m_proptool.type().size();
  std::string s1;
  for (int i=0; i<n; ++i) s1.append(" ");
  s1.append("|");
 
  std::string fieldmode[9] ={"NoField"       , "ConstantField", "SolenoidalField",
                 "ToroidalField" , "Grid3DField"  , "RealisticField" ,
                 "UndefinedField", "AthenaField"  , "?????"         };
 
  int mode = m_fieldprop.magneticFieldMode();
  if (mode<0 || mode>8 ) mode = 8;
 
  n     = 62-fieldmode[mode].size();
  std::string s3;
  for (int i=0; i<n; ++i) s3.append(" ");
  s3.append("|");
 
  n     = 62-m_updatortool.type().size();
  std::string s4;
  for (int i=0; i<n; ++i) s4.append(" ");
  s4.append("|");
 
  n     = 62-m_riocreator.type().size();
  std::string s5;
  for (int i=0; i<n; ++i) s5.append(" ");
  s5.append("|");
 
  n     = 62-m_pixcontainerkey.key().size();
  std::string s7;
  for (int i=0; i<n; ++i) s7.append(" ");
  s7.append("|");
 
  n     = 62-m_sctcontainerkey.key().size();
  std::string s8;
  for (int i=0; i<n; ++i) s8.append(" ");
  s8.append("|");
 
  out<<"|----------------------------------------------------------------------"
     <<"-------------------|"
     <<std::endl;
  if (m_usePIX) {
    out<<"| Pixel clusters location | "<<m_pixcontainerkey.key() <<s7<<std::endl;
  }
  if (m_useSCT) {
    out<<"| SCT   clusters location | "<<m_sctcontainerkey.key() <<s8<<std::endl;
  }
  out<<"| Tool for propagation    | "<<m_proptool   .type()<<s1<<std::endl;
  out<<"| Tool for updator        | "<<m_updatortool.type()<<s4<<std::endl;
  out<<"| Tool for rio  on track  | "<<m_riocreator .type()<<s5<<std::endl;
  out<<"| Magnetic field mode     | "<<fieldmode[mode]     <<s3<<std::endl;
  out<<"|----------------------------------------------------------------------"
     <<"-------------------|"
     <<std::endl;
  return out;
}

◆ dumpevent()

MsgStream & InDet::SiCombinatorialTrackFinder_xk::dumpevent	(	SiCombinatorialTrackFinderData_xk &	data,
		MsgStream &	out )

staticprivate

Definition at line 194 of file SiCombinatorialTrackFinder_xk.cxx.

{
  out<<"|---------------------------------------------------------------------|"
     <<std::endl;
  out<<"| Min pT of track (MeV)   | "<<std::setw(12)<<std::setprecision(5)<<data.pTmin()
     <<"                              |"<<std::endl;
  out<<"| Max Xi2 for cluster     | "<<std::setw(12)<<std::setprecision(5)<<data.xi2max()
     <<"                              |"<<std::endl;
  out<<"| Max Xi2 for outlayer    | "<<std::setw(12)<<std::setprecision(5)<<data.xi2maxNoAdd()
     <<"                              |"<<std::endl;
  out<<"| Max Xi2 for link        | "<<std::setw(12)<<std::setprecision(5)<<data.xi2maxlink()
     <<"                              |"<<std::endl;
  out<<"| Min number of clusters  | "<<std::setw(12)<<data.nclusmin()
     <<"                              |"<<std::endl;
  out<<"| Min number of wclusters | "<<std::setw(12)<<data.nwclusmin()
     <<"                              |"<<std::endl;
  out<<"| Max number holes        | "<<std::setw(12)<<data.nholesmax()
     <<"                              |"<<std::endl;
  out<<"| Max holes  gap          | "<<std::setw(12)<<data.dholesmax()
     <<"                              |"<<std::endl;
  out<<"| Use PRD to track assoc.?| "<<std::setw(12)<<(data.PRDtoTrackMap() ? "yes" : "no ")
     <<"                              |"<<std::endl;
  out<<"|---------------------------------------------------------------------|"
     <<std::endl;
  out<<"| Number input     seeds  | "<<std::setw(12)<<data.inputseeds()
     <<"                              |"<<std::endl;
  out<<"| Number accepted  seeds  | "<<std::setw(12)<<data.goodseeds()
     <<"                              |"<<std::endl;
  out<<"| Number initial  tracks  | "<<std::setw(12)<<data.inittracks()
     <<"                              |"<<std::endl;
  out<<"| Number wrong DE  roads  | "<<std::setw(12)<<data.roadbug()
     <<"                              |"<<std::endl;
  out<<"| Number output   tracks  | "<<std::setw(12)<<data.findtracks()
     <<"                              |"<<std::endl;
  out<<"|---------------------------------------------------------------------|"
     <<std::endl;
  return out;
}

◆ endEvent()

void InDet::SiCombinatorialTrackFinder_xk::endEvent ( SiCombinatorialTrackFinderData_xk & data ) const

overridevirtual

Definition at line 303 of file SiCombinatorialTrackFinder_xk.cxx.

{
  if (not data.isInitialized()) initializeCombinatorialData(Gaudi::Hive::currentContext(), data);
 
  // Print event information
  //
  if (m_outputlevel<=0) {
    data.nprint() = 1;
    dump(data, msg(MSG::DEBUG));
  }
}

◆ fillStatistic()

void InDet::SiCombinatorialTrackFinder_xk::fillStatistic	(	SiCombinatorialTrackFinderData_xk &	data,
		std::array< bool, NumberOfStats > &	information ) const

overrideprivatevirtual

Definition at line 1126 of file SiCombinatorialTrackFinder_xk.cxx.

{
   for(int i=0; i!=NumberOfStats; ++i) information[i] = data.statistic()[i];
}

◆ finalize()

StatusCode InDet::SiCombinatorialTrackFinder_xk::finalize ( )

overridevirtual

Definition at line 110 of file SiCombinatorialTrackFinder_xk.cxx.

{
  return AlgTool::finalize();
}

◆ findTrack()

InDet::SiCombinatorialTrackFinder_xk::EStat_t InDet::SiCombinatorialTrackFinder_xk::findTrack	(	SiCombinatorialTrackFinderData_xk &	data,
		const Trk::TrackParameters &	Tp,
		const std::vector< const Trk::SpacePoint * > &	Sp,
		const std::vector< Amg::Vector3D > &	Gp,
		std::vector< const InDetDD::SiDetectorElement * > &	DE,
		std::multimap< const Trk::PrepRawData , const Trk::Track > &	PT,
		const EventContext &	ctx ) const

private

init event data

populate a list of boundary links for the detector elements on our search road

Retrieve cached pointers to SG collections, or create the cache

Cluster list preparationn

in inside-out track finding, Sp.size() is typically 3 for TRT-seeded backtracking, it is 2 both applications go into this branch

returns false if two clusters are on the same detector element

use case if we have a set of global positions rather than space points to start from

use case if we have neither space-points nor global posittions, but track parameters to start from

Build initial trajectory

This will initialize the trajectory using the clusters we have and the parameter estimate

if the initialisation fails (indicating this is probably a bad attempt) and we are running with global positions instead of seeding

reset our cluster list

try again using the clusters from the track parameters only

if it worked now, set the quality flag to true

this can never happen?!

this can never happen either?!

if the last cluster on track is in the pixels, this is assumed to come from a pixel seed

max #iterations

Track finding

Strategy for pixel seeds

refine if needed

check if we found enough clusters

case of a strip seed or mixed PPS

first application of hit cut

backward smooting

apply hit cut again following smoothing step

refine if needed

quality cut

refine the hole cut

Definition at line 583 of file SiCombinatorialTrackFinder_xk.cxx.

{
  if (not data.isInitialized()) initializeCombinatorialData(ctx, data);
  std::vector<const InDet::SiDetElementBoundaryLink_xk*> DEL;
  detectorElementLinks(DE, DEL,ctx);
  const InDet::PixelClusterContainer* p_pixcontainer = data.pixContainer();
  if (m_usePIX && !p_pixcontainer) {
    SG::ReadHandle<InDet::PixelClusterContainer> pixcontainer(m_pixcontainerkey,ctx);
    if(!pixcontainer.isValid()){
      //the pixel container is not valid
      //we are supposed to use pixel so stop here
      REPORT_ERROR (StatusCode::FAILURE)<< "Unable to dereference pix container, handle not valid";
      return WrongInit;
    }
 
    p_pixcontainer = pixcontainer.ptr();
    data.setPixContainer(p_pixcontainer);
  }
  const InDet::SCT_ClusterContainer* p_sctcontainer = data.sctContainer();
  if (m_useSCT && !p_sctcontainer) {
    SG::ReadHandle<InDet::SCT_ClusterContainer> sctcontainer(m_sctcontainerkey,ctx);
 
    if(!sctcontainer.isValid()){
      //the sct container is not valid
      //we are supposed to use sct so stop here
      REPORT_ERROR (StatusCode::FAILURE)<< "Unable to dereference sct container, handle not valid";
      return WrongInit;
    }
 
    p_sctcontainer = sctcontainer.ptr();
    data.setSctContainer(p_sctcontainer);
  }
  std::vector<const InDet::SiCluster*> Cl;
  bool isTwoPointSeed = false;
  if (Sp.size() > 1) {
    if (!spacePointsToClusters(Sp,Cl)) {
      return TwoCluster;
    }
    if (Sp.size()==2) isTwoPointSeed = true;
  }
  else if (Gp.size() > 2) {
    if (!data.trajectory().globalPositionsToClusters(p_pixcontainer, p_sctcontainer, Gp, DEL, PT, Cl)) return TwoCluster;
  } else {
    if (!data.trajectory().trackParametersToClusters(p_pixcontainer, p_sctcontainer, Tp, DEL, PT, Cl, ctx)) return TwoCluster;
  }
  ++data.goodseeds();
  bool Qr;
  bool Q = data.trajectory().initialize(m_usePIX, m_useSCT, p_pixcontainer, p_sctcontainer, Tp, Cl, DEL, Qr,ctx);
  if (!Q && Sp.size() < 2 && Gp.size() > 3) {
    Cl.clear();
    if (!data.trajectory().trackParametersToClusters(p_pixcontainer, p_sctcontainer, Tp, DEL, PT, Cl,ctx)) return TwoCluster;
 
    if (!data.trajectory().initialize(m_usePIX, m_useSCT, p_pixcontainer, p_sctcontainer, Tp, Cl, DEL, Qr,ctx)) return TwoCluster;
 
    Q = Qr = true;
  }
  if (!Qr){
    ++data.roadbug();
    return WrongRoad;
  }
  if (!Q) return WrongInit;
 
  ++data.inittracks();
  bool pixseed = data.trajectory().isLastPixel();
  int itmax    = 30;
  if (!data.useFastTracking() and data.simpleTrack()) itmax = 10;
  if (data.heavyIon()) itmax = 50;
 
  //
  bool toReturnFailedTrack = data.flagToReturnFailedTrack();
  if( toReturnFailedTrack ) data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::Unrecoverable);
  if (pixseed) {      
    if (!data.trajectory().forwardExtension (false,itmax,ctx)) return CantFindTrk;
    if (!data.trajectory().backwardSmoother (false,ctx)      ) return CantFindTrk;
    if (!data.trajectory().backwardExtension(itmax,ctx)      ) return CantFindTrk;
    if (data.isITkGeometry() &&
    (data.trajectory().nclusters() < data.nclusmin() ||
     data.trajectory().ndf() < data.nwclusmin()) ) return CantFindTrk;
    if(!data.useFastTracking() && data.trajectory().difference() > 0){
      if (!data.trajectory().forwardFilter(ctx)) {
    if( toReturnFailedTrack ) {
      data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::PixSeedDiffKFFwd);
    }
    else {
      return CantFindTrk;
    }
      }
 
      if (!data.trajectory().backwardSmoother (false,ctx) ) {
    if( toReturnFailedTrack ) {
      data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::PixSeedDiffKFBwd);
    }
    else {
      return CantFindTrk;
    }
      }
    }
 
    int na = data.trajectory().nclustersNoAdd();
    if (data.trajectory().nclusters()+na < data.nclusmin() ||
    data.trajectory().ndf() < data.nwclusmin()) {
       if( toReturnFailedTrack ) {
      data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::PixSeedNCluster);
       }
       else {
      return CantFindTrk;
       }
    }
  }
  else {      // Strategy for mixed seeds
    if (!data.trajectory().backwardSmoother(isTwoPointSeed,ctx)       ) return CantFindTrk;
    if (!data.trajectory().backwardExtension(itmax,ctx)    ) return CantFindTrk;
    if (!data.trajectory().forwardExtension(true,itmax,ctx)) return CantFindTrk;
    int na = data.trajectory().nclustersNoAdd();
    if (data.trajectory().nclusters()+na < data.nclusmin() ||
    data.trajectory().ndf() < data.nwclusmin()) return CantFindTrk;
    if (!data.trajectory().backwardSmoother(false,ctx)    ) return CantFindTrk;
    na     = data.trajectory().nclustersNoAdd();
    if (data.trajectory().nclusters()+na < data.nclusmin() ||
    data.trajectory().ndf() < data.nwclusmin()) {
       if( toReturnFailedTrack ) {
      data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::MixSeedNCluster);
       }
       else {
      return CantFindTrk;
       }
    }
    if (data.trajectory().difference() > 0) {
      if (!data.trajectory().forwardFilter(ctx)) {
     if( toReturnFailedTrack ) {
        data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::MixSeedDiffKFFwd);
     }
     else {
        return CantFindTrk;
     }
      }
      if (!data.trajectory().backwardSmoother (false, ctx)) {
     if( toReturnFailedTrack ) {
        data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::MixSeedDiffKFBwd);
     }
     else {
        return CantFindTrk;
     }
      }
    }
  }
  if (data.trajectory().qualityOptimization() < (m_qualityCut*data.nclusmin())) {
     if( toReturnFailedTrack ) {
    data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::Quality);
     }
     else {
    return CantFindTrk;
     }
  }
 
  if (data.trajectory().pTfirst  () < data.pTmin   () &&
      data.trajectory().nclusters() < data.nclusmin()) {
     if( toReturnFailedTrack ) {
    data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::Pt);
     }
     else {
    return CantFindTrk;
     }
  }
 
  if (data.trajectory().nclusters() < data.nclusminb() ||
      data.trajectory().ndf      () < data.nwclusmin()) {
     if( toReturnFailedTrack ) {
    data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::NCluster);
     }
     else {
    return CantFindTrk;
     }
  }
  if (m_writeHolesFromPattern){
    data.trajectory().updateHoleSearchResult();
    if (!data.trajectory().getHoleSearchResult().passPatternHoleCut) {
       if( toReturnFailedTrack ) {
      data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::HoleCut);
       }
       else {
      return CantFindTrk;
       }
    }
  }
 
  if( toReturnFailedTrack ) {
     if(  data.resultCode() != SiCombinatorialTrackFinderData_xk::ResultCode::Unrecoverable ) {
    return CantFindTrk;
     }
     else {
    data.setResultCode(SiCombinatorialTrackFinderData_xk::ResultCode::Success);
     }
  }
 
  return Success;
}

◆ getTrackQualityCuts()

void InDet::SiCombinatorialTrackFinder_xk::getTrackQualityCuts	(	SiCombinatorialTrackFinderData_xk &	data,
		const TrackQualityCuts &	Cuts )

staticprivate

Definition at line 1023 of file SiCombinatorialTrackFinder_xk.cxx.

{
  // Integer cuts
  //
  int intCut = 0;
 
  if (!Cuts.getIntCut("CosmicTrack"         , intCut)) intCut = 0;
  data.setCosmicTrack(intCut);
 
  if (!Cuts.getIntCut("MinNumberOfClusters" , intCut)) intCut = 7;
  data.setNclusmin(intCut);
 
  data.setNclusminb(std::max(3, data.nclusmin()-1));
 
  if (!Cuts.getIntCut("MinNumberOfWClusters", intCut)) intCut = 7;
  data.setNwclusmin(intCut);
 
  if (!Cuts.getIntCut("MaxNumberOfHoles"    , intCut)) intCut = 2;
  if(!data.cosmicTrack() && intCut>2) intCut = 2;
  data.setNholesmax(intCut);
 
  if (!Cuts.getIntCut("MaxHolesGap"         , intCut)) intCut = 2;
  if(!data.cosmicTrack() && intCut>2) intCut = 2;
  if(intCut > data.nholesmax())       intCut = data.nholesmax();
  data.setDholesmax(intCut);
 
  data.tools().setHolesClusters(data.nholesmax(), data.dholesmax(),
                data.nclusmin());
 
  if (!Cuts.getIntCut("UseAssociationTool"  , intCut)) intCut = 0;
  data.tools().setAssociation(intCut);
 
  if (!Cuts.getIntCut("SimpleTrack"         , intCut)) intCut = 0;
  data.setSimpleTrack(bool(intCut));
 
  // Double cuts
  //
  double doubleCut = 0.;
 
  if (!Cuts.getDoubleCut("pTmin"             , doubleCut)) doubleCut = 500.;
  data.setPTmin(doubleCut);
 
  if (!Cuts.getDoubleCut("pTminBrem"         , doubleCut)) doubleCut = 1000.;
  data.setPTminBrem(doubleCut);
 
  if (!Cuts.getDoubleCut("MaxXi2forCluster"  , doubleCut)) doubleCut = 9.;
  data.setXi2max(doubleCut);
 
  if (!Cuts.getDoubleCut("MaxXi2forOutlier"  , doubleCut)) doubleCut = 25.;
  if (!data.cosmicTrack() && doubleCut > 25.) doubleCut = 25.;
  if (doubleCut <= data.xi2max()) doubleCut = data.xi2max()+5.;
  data.setXi2maxNoAdd(doubleCut);
 
  if (!Cuts.getDoubleCut("MaxXi2forSearch"    , doubleCut)) doubleCut = 100.;
  data.setXi2maxlink(doubleCut);
 
  data.tools().setXi2pTmin(data.xi2max(), data.xi2maxNoAdd(),
               data.xi2maxlink(), data.pTmin());
 
  if (!Cuts.getIntCut   ("doMultiTracksProd", intCut)) intCut = 0;
  if (!Cuts.getDoubleCut("MaxXi2MultiTracks", doubleCut)) doubleCut = 7.;
  if (!data.cosmicTrack() && doubleCut > 7.) doubleCut = 7.;
  data.tools().setMultiTracks(intCut, doubleCut);
 
  data.trajectory().setParameters();
}

◆ getTracks() [1/2]

const std::list< Trk::Track * > & InDet::SiCombinatorialTrackFinder_xk::getTracks	(	SiCombinatorialTrackFinderData_xk &	data,
		const Trk::TrackParameters &	Tp,
		const std::vector< const Trk::SpacePoint * > &	Sp,
		const std::vector< Amg::Vector3D > &	Gp,
		std::vector< const InDetDD::SiDetectorElement * > &	DE,
		const TrackQualityCuts &	Cuts,
		const EventContext &	ctx ) const

overridevirtual

turn off brem fit & electron flags

remove existing tracks

try to find the tracks

if we didn't find anything, bail out

sort in step order

Definition at line 319 of file SiCombinatorialTrackFinder_xk.cxx.

{
 
  if (not data.isInitialized()) initializeCombinatorialData(ctx, data);
 
  data.statistic().fill(false);
  data.tools().setBremNoise(false, false);
  data.tracks().erase(data.tracks().begin(), data.tracks().end());
 
  ++data.inputseeds();
  if (!m_usePIX && !m_useSCT) {
    return data.tracks();
  }
 
  // Get track quality cuts information and write them into the event data... again?
  getTrackQualityCuts(data, Cuts);
  std::multimap<const Trk::PrepRawData*, const Trk::Track*> PT;
  EStat_t FT = findTrack(data, Tp, Sp, Gp, DE, PT, ctx);
  if(FT!=Success) {
    data.statistic()[FT] = true;
    if( ! data.flagToReturnFailedTrack() ) return data.tracks();
  }
  data.trajectory().sortStep();
 
 
  Trk::Track* t = convertToTrack(data,ctx);
  if (t) {
     ++data.findtracks();
     if (m_writeHolesFromPattern) data.addPatternHoleSearchOutcome(t,data.trajectory().getHoleSearchResult());
     data.tracks().push_back(t);
  }
 
  if (!data.tools().multiTrack() ||
      data.simpleTrack() ||
      Sp.size()<=2 ||
      data.cosmicTrack() ||
      data.trajectory().pTfirst() < data.tools().pTmin()) return data.tracks();
 
  while ((t=convertToNextTrack(data))) {
    ++data.findtracks();
    if (m_writeHolesFromPattern) data.addPatternHoleSearchOutcome(t,data.trajectory().getHoleSearchResult());
    data.tracks().push_back(t);
  }
  return data.tracks();
}

◆ getTracks() [2/2]

const std::list< Trk::Track * > & InDet::SiCombinatorialTrackFinder_xk::getTracks	(	SiCombinatorialTrackFinderData_xk &	data,
		const Trk::TrackParameters &	Tp,
		const std::vector< const Trk::SpacePoint * > &	Sp,
		const std::vector< Amg::Vector3D > &	Gp,
		std::vector< const InDetDD::SiDetectorElement * > &	DE,
		std::multimap< const Trk::PrepRawData , const Trk::Track > &	PT,
		const EventContext &	ctx ) const

overridevirtual

Definition at line 385 of file SiCombinatorialTrackFinder_xk.cxx.

{
 
  if (not data.isInitialized()) initializeCombinatorialData(ctx, data);
 
  data.tools().setBremNoise(false, false);
  data.tracks().erase(data.tracks().begin(), data.tracks().end());
 
  data.statistic().fill(false);
  ++data.inputseeds();
  if (!m_usePIX && !m_useSCT) {
    return data.tracks();
  }
 
  EStat_t FT = findTrack(data, Tp, Sp, Gp, DE, PT,ctx);
  if(FT!=Success) {
    data.statistic()[FT] = true;
    if( ! data.flagToReturnFailedTrack() || std::abs(data.resultCode()) < SiCombinatorialTrackFinderData_xk::ResultCodeThreshold::RecoverableForDisTrk ) return data.tracks();
  }
  if (!data.trajectory().isNewTrack(PT))
  {
     data.statistic()[NotNewTrk] = true;
     return data.tracks();
  }
  data.trajectory().sortStep();
 
  if(data.useFastTracking()) {
    if(!data.trajectory().filterWithPreciseClustersError(ctx)) {
      data.statistic()[CantFindTrk] = true;
      return data.tracks();
    }
  }
 
  // Trk::Track production
  //
  Trk::Track* t = convertToTrack(data,ctx);
  if (t==nullptr) return data.tracks(); // @TODO should one check if convertToNextTrack would yield anything ?
  if (m_writeHolesFromPattern) data.addPatternHoleSearchOutcome(t,data.trajectory().getHoleSearchResult());
 
  ++data.findtracks();
  data.tracks().push_back(t);
 
  if (!data.tools().multiTrack() ||
      data.simpleTrack() ||
      Sp.size()<=2 ||
      data.cosmicTrack() ||
      data.trajectory().pTfirst() < data.tools().pTmin()) return data.tracks();
 
  while ((t=convertToNextTrack(data))) {
    if (m_writeHolesFromPattern) data.addPatternHoleSearchOutcome(t,data.trajectory().getHoleSearchResult());
    ++data.findtracks();
    data.tracks().push_back(t);
  }
  return data.tracks();
}

◆ getTracksWithBrem()

const std::list< Trk::Track * > & InDet::SiCombinatorialTrackFinder_xk::getTracksWithBrem	(	SiCombinatorialTrackFinderData_xk &	data,
		const Trk::TrackParameters &	Tp,
		const std::vector< const Trk::SpacePoint * > &	Sp,
		const std::vector< Amg::Vector3D > &	Gp,
		std::vector< const InDetDD::SiDetectorElement * > &	DE,
		std::multimap< const Trk::PrepRawData , const Trk::Track > &	PT,
		bool	isCaloCompatible,
		const EventContext &	ctx ) const

overridevirtual

Definition at line 453 of file SiCombinatorialTrackFinder_xk.cxx.

{
 
  if (not data.isInitialized()) initializeCombinatorialData(ctx, data);
 
  data.statistic().fill(false);
 
  // Old information
  //
  int mult = 0;
  if (data.tools().multiTrack()) mult = 1;
  double Xi2m = data.tools().xi2multi();
 
  data.tools().setBremNoise(false, true);
  data.tracks().erase(data.tracks().begin(), data.tracks().end());
 
  ++data.inputseeds();
  if (!m_usePIX && !m_useSCT) {
    return data.tracks();
  }
 
  EStat_t FT = findTrack(data,Tp,Sp,Gp,DE,PT,ctx);
 
  bool Q = (FT==Success);
  if (Q){
    Q = data.trajectory().isNewTrack(PT);
  }
 
  int na = 0;
  if (Q) {
    data.trajectory().sortStep();
 
    // Trk::Track production
    //
    Trk::TrackInfo oldinfo = data.trackinfo();
    if (isCaloCompatible) data.trackinfo().setPatternRecognitionInfo(Trk::TrackInfo::TrackInCaloROI);
 
    data.tools().setMultiTracks(0, Xi2m);
    Trk::Track* t = convertToTrack(data,ctx);
    data.trackinfo() = oldinfo;
    data.tools().setMultiTracks(mult,Xi2m);
 
    if (t==nullptr) return data.tracks(); // @TODO should one check whether the next findTrack call would yield something ?
    if (m_writeHolesFromPattern) data.addPatternHoleSearchOutcome(t,data.trajectory().getHoleSearchResult());
    ++data.findtracks();
    data.tracks().push_back(t);
    na = data.trajectory().nclusters();
    if (na >=12 && !data.trajectory().nclustersNoAdd()) return data.tracks();
 
    if (data.trajectory().pTfirst() < data.pTminBrem()) return data.tracks();
  }
  if ((*Sp.begin())->clusterList().second) return data.tracks();
 
  // Repeat track finding using electron noise model
  //
  data.statistic()[BremAttempt] = true;
 
  data.tools().setBremNoise(true,true);
  FT = findTrack(data, Tp, Sp, Gp, DE, PT,ctx);
 
  if (FT!=Success) {
    data.statistic()[FT] = true;
    return data.tracks();
  }
 
  if (!data.trajectory().isNewTrack(PT)) { data.statistic()[NotNewTrk] = true; return data.tracks(); }
 
  int nb = data.trajectory().nclusters();
  if (nb <= na) return data.tracks();
 
  data.trajectory().sortStep();
 
  if(data.useFastTracking()) {
    if(!data.trajectory().filterWithPreciseClustersError(ctx)) {
      data.statistic()[CantFindTrk] = true;
      return data.tracks();
    }
  }
 
  // Trk::Track production
  //
  Trk::TrackInfo oldinfo = data.trackinfo();
  data.trackinfo().setTrackProperties(Trk::TrackInfo::BremFit          );
  data.trackinfo().setTrackProperties(Trk::TrackInfo::BremFitSuccessful);
  if (isCaloCompatible) data.trackinfo().setPatternRecognitionInfo(Trk::TrackInfo::TrackInCaloROI);
 
  data.tools().setMultiTracks(0, Xi2m);
  Trk::Track* t = convertToTrack(data, ctx);
  data.trackinfo() = oldinfo;
  data.tools().setMultiTracks(mult, Xi2m);
 
  if (t==nullptr) return data.tracks();
  if (m_writeHolesFromPattern) data.addPatternHoleSearchOutcome(t,data.trajectory().getHoleSearchResult());
 
  ++data.findtracks();
  data.tracks().push_back(t);
  return data.tracks();
}

◆ initialize()

StatusCode InDet::SiCombinatorialTrackFinder_xk::initialize ( )

overridevirtual

Definition at line 44 of file SiCombinatorialTrackFinder_xk.cxx.

{
  // Get RungeKutta propagator tool
  //
  if ( m_proptool.retrieve().isFailure() ) {
    ATH_MSG_FATAL("Failed to retrieve tool " << m_proptool);
    return StatusCode::FAILURE;
  }
  ATH_MSG_DEBUG("Retrieved tool " << m_proptool);
 
  // Get updator tool
  //
  if ( m_updatortool.retrieve().isFailure() ) {
    ATH_MSG_FATAL("Failed to retrieve tool " << m_updatortool);
    return StatusCode::FAILURE;
  }
  ATH_MSG_DEBUG("Retrieved tool " << m_updatortool);
 
  // Get RIO_OnTrack creator
  //
  if ( m_riocreator.retrieve().isFailure() ) {
    ATH_MSG_FATAL("Failed to retrieve tool " << m_riocreator);
    return StatusCode::FAILURE;
  }
  ATH_MSG_DEBUG("Retrieved tool " << m_riocreator);
 
  // disable pixel/SCT conditions summary tool: pixel/SCT are not used, the status event data is used and not being validated
  ATH_CHECK( m_pixelCondSummaryTool.retrieve( DisableTool{!m_usePIX || (!m_pixelDetElStatus.empty() && !VALIDATE_STATUS_ARRAY_ACTIVATED)} ) );
  ATH_CHECK( m_sctCondSummaryTool.retrieve(   DisableTool{!m_useSCT || (!m_sctDetElStatus.empty()   && !VALIDATE_STATUS_ARRAY_ACTIVATED)} ) );
  //
  // Get InDetBoundaryCheckTool
  if ( m_boundaryCheckTool.retrieve().isFailure() ) {
    ATH_MSG_FATAL("Failed to retrieve tool " << m_boundaryCheckTool);
    return StatusCode::FAILURE;
  }
  ATH_MSG_DEBUG("Retrieved tool " << m_boundaryCheckTool);
 
 
 
  // Setup callback for magnetic field
  //
  magneticFieldInit();
 
  // Get output print level
  //
  m_outputlevel = msg().level()-MSG::DEBUG;
 
  ATH_CHECK( m_pixcontainerkey.initialize (m_usePIX) );
  ATH_CHECK( m_boundaryPixelKey.initialize (m_usePIX) );
 
  ATH_CHECK( m_sctcontainerkey.initialize (m_useSCT) );
  ATH_CHECK( m_boundarySCTKey.initialize (m_useSCT) );
 
  // initialize conditions object key for field cache
  //
  ATH_CHECK( m_fieldCondObjInputKey.initialize() );
  ATH_CHECK( m_pixelDetElStatus.initialize( !m_pixelDetElStatus.empty() && m_usePIX) );
  ATH_CHECK( m_sctDetElStatus.initialize( !m_sctDetElStatus.empty() && m_useSCT) );
  m_minPt2Cut = std::pow(m_minPtCut.value(),2);
  return StatusCode::SUCCESS;
}

◆ initializeCombinatorialData()

void InDet::SiCombinatorialTrackFinder_xk::initializeCombinatorialData	(	const EventContext &	ctx,
		SiCombinatorialTrackFinderData_xk &	data ) const

private

Add conditions object to SiCombinatorialTrackFinderData to be able to access the field cache for each new event Get conditions object for field cache

Must have set fieldCondObj BEFORE calling setTools because fieldCondObj is used there

Definition at line 1091 of file SiCombinatorialTrackFinder_xk.cxx.

                                                                                                                                           {
  SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCondObjInputKey, ctx};
  const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
  if (fieldCondObj == nullptr) {
      std::string msg = "InDet::SiCombinatorialTrackFinder_xk::initializeCombinatorialData: Failed to retrieve AtlasFieldCacheCondObj with key " + m_fieldCondObjInputKey.key();
      throw(std::runtime_error(msg));
  }
  data.setFieldCondObj(fieldCondObj);
  data.setTools(&*m_proptool,
                &*m_updatortool,
                &*m_riocreator,
                (m_usePIX && (m_pixelDetElStatus.empty() || VALIDATE_STATUS_ARRAY_ACTIVATED)? &*m_pixelCondSummaryTool : nullptr),
                (m_useSCT && (m_sctDetElStatus.empty()   || VALIDATE_STATUS_ARRAY_ACTIVATED)? &*m_sctCondSummaryTool   : nullptr),
                &m_fieldprop,
                &*m_boundaryCheckTool);
  if (!m_pixelDetElStatus.empty()) {
     SG::ReadHandle<InDet::SiDetectorElementStatus> pixelDetElStatus( m_pixelDetElStatus, ctx);
     data.setPixelDetectorElementStatus( pixelDetElStatus.cptr() );
  }
  if (!m_sctDetElStatus.empty()) {
     SG::ReadHandle<InDet::SiDetectorElementStatus> sctDetElStatus( m_sctDetElStatus, ctx);
     data.setSCTDetectorElementStatus( sctDetElStatus.cptr() );
  }
 
  // Set the ITk Geometry setup
  data.setITkGeometry(m_ITkGeometry);
  // Set the ITk Fast Tracking setup
  data.setFastTracking(m_doFastTracking);
}

◆ magneticFieldInit()

void InDet::SiCombinatorialTrackFinder_xk::magneticFieldInit ( )

private

Definition at line 921 of file SiCombinatorialTrackFinder_xk.cxx.

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

◆ newEvent() [1/2]

void InDet::SiCombinatorialTrackFinder_xk::newEvent	(	const EventContext &	ctx,
		SiCombinatorialTrackFinderData_xk &	data ) const

overridevirtual

Definition at line 237 of file SiCombinatorialTrackFinder_xk.cxx.

{
  if (not data.isInitialized()) initializeCombinatorialData(ctx, data);
 
  // Erase statistic information
  //
  data.inputseeds() = 0;
  data.goodseeds()  = 0;
  data.inittracks() = 0;
  data.findtracks() = 0;
  data.roadbug()    = 0;
 
  // Set track info
  //
  data.trackinfo().setPatternRecognitionInfo(Trk::TrackInfo::SiSPSeededFinder);
  data.setCosmicTrack(0);
 
  // Add conditions object to SiCombinatorialTrackFinderData to be able to access the field cache for each new event
  // Get conditions object for field cache
  SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCondObjInputKey, ctx};
  const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
  if (fieldCondObj == nullptr) {
      std::string msg = "InDet::SiCombinatorialTrackFinder_xk::newEvent: Failed to retrieve AtlasFieldCacheCondObj with key " + m_fieldCondObjInputKey.key();
      throw(std::runtime_error(msg));
  }
  data.setFieldCondObj(fieldCondObj);
}

◆ newEvent() [2/2]

void InDet::SiCombinatorialTrackFinder_xk::newEvent	(	const EventContext &	ctx,
		SiCombinatorialTrackFinderData_xk &	data,
		Trk::TrackInfo	info,
		const TrackQualityCuts &	Cuts ) const

overridevirtual

Get track quality cuts information from argument and write it into the event data

update pattern recognition flags in the event data based on track info arg

Definition at line 269 of file SiCombinatorialTrackFinder_xk.cxx.

{
 
  if (not data.isInitialized()) {
    //Check if to use PRDAssociation before initializing all the tools
    int useasso;
    if(!Cuts.getIntCut   ("UseAssociationTool"  ,useasso      )) useasso = 0;
    data.tools().setAssociation(useasso);
 
    initializeCombinatorialData(ctx, data);
  }
 
  newEvent(ctx, data);
  data.trackinfo() = info;
  getTrackQualityCuts(data, Cuts);
 
  data.setHeavyIon(false);
  data.setCosmicTrack(0);
  if (info.patternRecoInfo(Trk::TrackInfo::SiSpacePointsSeedMaker_Cosmic))
    data.setCosmicTrack(1);
  else if (info.patternRecoInfo(Trk::TrackInfo::SiSpacePointsSeedMaker_HeavyIon))
    data.setHeavyIon(true);
 
}

◆ pTseed()

double InDet::SiCombinatorialTrackFinder_xk::pTseed	(	SiCombinatorialTrackFinderData_xk &	data,
		const Trk::TrackParameters &	Tp,
		const std::vector< const Trk::SpacePoint * > &	Sp,
		const EventContext &	ctx ) const

overridevirtual

Definition at line 564 of file SiCombinatorialTrackFinder_xk.cxx.

{
  std::vector<const InDet::SiCluster*>           Cl;
  std::vector<const InDetDD::SiDetectorElement*> DE;
  if(!spacePointsToClusters(Sp,Cl,DE)) return 0.;
 
  std::vector<const InDet::SiDetElementBoundaryLink_xk*> DEL;
  detectorElementLinks(DE,DEL,ctx);
  return data.trajectory().pTseed(Tp,Cl,DEL,ctx);
}

◆ spacePointsToClusters()

bool InDet::SiCombinatorialTrackFinder_xk::spacePointsToClusters	(	const std::vector< const Trk::SpacePoint * > &	Sp,
		std::vector< const InDet::SiCluster * > &	Sc,
		std::optional< std::reference_wrapper< std::vector< const InDetDD::SiDetectorElement * > > >	DE = std::nullopt )

staticprivate

loop over all SP

get the first cluster on an SP

add to list

for strips, also make sure to pick up the second one!

Detector elments test

here we reject cases where two subsequent clusters are on the same detector element

Definition at line 934 of file SiCombinatorialTrackFinder_xk.cxx.

{
  Sc.reserve(Sp.size());
  for (const Trk::SpacePoint* s : Sp) {
     const Trk::PrepRawData* p = s->clusterList().first;
     if (p) {
        const InDet::SiCluster* c = static_cast<const InDet::SiCluster*>(p);
        if (c){
          Sc.push_back(c);
        }
     }
     p = s->clusterList().second;
     if (p) {
        const InDet::SiCluster* c = static_cast<const InDet::SiCluster*>(p);
        if (c){
          Sc.push_back(c);
        }
     }
  }
  std::vector<const InDet::SiCluster*>::iterator cluster = Sc.begin(), nextCluster, endClusters = Sc.end();
  if (DE) {
    DE->get().reserve(Sc.size());
  }
  for (; cluster != endClusters; ++cluster) {
 
     const InDetDD::SiDetectorElement* de = (*cluster)->detectorElement();
 
     nextCluster = cluster;
     ++nextCluster;
     for (; nextCluster != endClusters; ++nextCluster) {
        if (de == (*nextCluster)->detectorElement()){
          return false;
        }
     }
     if (DE) {
       DE->get().push_back(de);
     }
  }
  return true;
}

◆ SiTrajectory_xk

friend class SiTrajectory_xk

friend

Definition at line 64 of file SiCombinatorialTrackFinder_xk.h.

Member Data Documentation

◆ m_boundaryCheckTool

ToolHandle<Trk::IBoundaryCheckTool> InDet::SiCombinatorialTrackFinder_xk::m_boundaryCheckTool

private

Initial value:

{this, "BoundaryCheckTool",

"InDet::InDetBoundaryCheckTool", "Boundary checking tool for detector sensitivities"}

Definition at line 149 of file SiCombinatorialTrackFinder_xk.h.

149 {this, "BoundaryCheckTool",

150 "InDet::InDetBoundaryCheckTool", "Boundary checking tool for detector sensitivities"};

◆ m_boundaryPixelKey

SG::ReadCondHandleKey<InDet::SiDetElementBoundaryLinks_xk> InDet::SiCombinatorialTrackFinder_xk::m_boundaryPixelKey

private

Initial value:

{this, "PixelDetElementBoundaryLinks_xk",

"PixelDetElementBoundaryLinks_xk", "Key of InDet::SiDetElementBoundaryLinks_xk for Pixel"}

Definition at line 166 of file SiCombinatorialTrackFinder_xk.h.

166 {this, "PixelDetElementBoundaryLinks_xk",

167 "PixelDetElementBoundaryLinks_xk", "Key of InDet::SiDetElementBoundaryLinks_xk for Pixel"};

◆ m_boundarySCTKey

SG::ReadCondHandleKey<InDet::SiDetElementBoundaryLinks_xk> InDet::SiCombinatorialTrackFinder_xk::m_boundarySCTKey

private

Initial value:

{this, "SCT_DetElementBoundaryLinks_xk",

"SCT_DetElementBoundaryLinks_xk", "Key of InDet::SiDetElementBoundaryLinks_xk for SCT"}

Definition at line 168 of file SiCombinatorialTrackFinder_xk.h.

168 {this, "SCT_DetElementBoundaryLinks_xk",

169 "SCT_DetElementBoundaryLinks_xk", "Key of InDet::SiDetElementBoundaryLinks_xk for SCT"};

◆ m_doFastTracking

BooleanProperty InDet::SiCombinatorialTrackFinder_xk::m_doFastTracking {this, "doFastTracking", false}

private

Definition at line 191 of file SiCombinatorialTrackFinder_xk.h.

191{this, "doFastTracking", false};

◆ m_fieldCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> InDet::SiCombinatorialTrackFinder_xk::m_fieldCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"}

private

Definition at line 171 of file SiCombinatorialTrackFinder_xk.h.

171{this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"};

◆ m_fieldmode

StringProperty InDet::SiCombinatorialTrackFinder_xk::m_fieldmode {this, "MagneticFieldMode", "MapSolenoid", "Mode of magnetic field"}

private

Definition at line 192 of file SiCombinatorialTrackFinder_xk.h.

192{this, "MagneticFieldMode", "MapSolenoid", "Mode of magnetic field"};

◆ m_fieldprop

Trk::MagneticFieldProperties InDet::SiCombinatorialTrackFinder_xk::m_fieldprop

private

Magnetic field properties.

Definition at line 202 of file SiCombinatorialTrackFinder_xk.h.

◆ m_ITkGeometry

BooleanProperty InDet::SiCombinatorialTrackFinder_xk::m_ITkGeometry {this, "ITkGeometry", false}

private

Definition at line 190 of file SiCombinatorialTrackFinder_xk.h.

190{this, "ITkGeometry", false};

◆ m_minPt2Cut

float InDet::SiCombinatorialTrackFinder_xk::m_minPt2Cut =0

private

Definition at line 195 of file SiCombinatorialTrackFinder_xk.h.

◆ m_minPtCut

FloatProperty InDet::SiCombinatorialTrackFinder_xk::m_minPtCut {this, "MinFinalPtCut", 100, "Cut on the pt of the final track. Must be >0 to avoid NANs when computing eta."}

private

Definition at line 194 of file SiCombinatorialTrackFinder_xk.h.

194{this, "MinFinalPtCut", 100, "Cut on the pt of the final track. Must be >0 to avoid NANs when computing eta."};

◆ m_outputlevel

int InDet::SiCombinatorialTrackFinder_xk::m_outputlevel {0}

private

Definition at line 201 of file SiCombinatorialTrackFinder_xk.h.

201{0};

◆ m_pixcontainerkey

SG::ReadHandleKey<InDet::PixelClusterContainer> InDet::SiCombinatorialTrackFinder_xk::m_pixcontainerkey

private

Initial value:

{this, "PixelClusterContainer",

"PixelClusters"}

Definition at line 162 of file SiCombinatorialTrackFinder_xk.h.

162 {this, "PixelClusterContainer",

163 "PixelClusters"};

◆ m_pixelCondSummaryTool

ToolHandle<IInDetConditionsTool> InDet::SiCombinatorialTrackFinder_xk::m_pixelCondSummaryTool

private

Initial value:

{this, "PixelSummaryTool",

"PixelConditionsSummaryTool"}

Definition at line 145 of file SiCombinatorialTrackFinder_xk.h.

145 {this, "PixelSummaryTool",

146 "PixelConditionsSummaryTool"};

◆ m_pixelDetElStatus

SG::ReadHandleKey<InDet::SiDetectorElementStatus> InDet::SiCombinatorialTrackFinder_xk::m_pixelDetElStatus {this, "PixelDetElStatus", "", "Key of SiDetectorElementStatus for Pixel"}

private

Optional read handle to get status data to test whether a pixel detector element is good.

If set to e.g. PixelDetectorElementStatus the event data will be used instead of the pixel conditions summary tool.

Definition at line 177 of file SiCombinatorialTrackFinder_xk.h.

178{this, "PixelDetElStatus", "", "Key of SiDetectorElementStatus for Pixel"};

◆ m_proptool

PublicToolHandle<Trk::IPatternParametersPropagator> InDet::SiCombinatorialTrackFinder_xk::m_proptool

private

Initial value:

{this, "PropagatorTool",

"Trk::RungeKuttaPropagator/InDetPropagator"}

Definition at line 152 of file SiCombinatorialTrackFinder_xk.h.

152 {this, "PropagatorTool",

153 "Trk::RungeKuttaPropagator/InDetPropagator"};

◆ m_qualityCut

DoubleProperty InDet::SiCombinatorialTrackFinder_xk::m_qualityCut {this, "TrackQualityCut", 9.3, "Simple track quality cut"}

private

Definition at line 193 of file SiCombinatorialTrackFinder_xk.h.

193{this, "TrackQualityCut", 9.3, "Simple track quality cut"};

◆ m_riocreator

PublicToolHandle<Trk::IRIO_OnTrackCreator> InDet::SiCombinatorialTrackFinder_xk::m_riocreator

private

Initial value:

{this, "RIOonTrackTool",

"Trk::RIO_OnTrackCreator/RIO_OnTrackCreator"}

Definition at line 156 of file SiCombinatorialTrackFinder_xk.h.

156 {this, "RIOonTrackTool",

157 "Trk::RIO_OnTrackCreator/RIO_OnTrackCreator"};

◆ m_sctCondSummaryTool

ToolHandle<IInDetConditionsTool> InDet::SiCombinatorialTrackFinder_xk::m_sctCondSummaryTool

private

Initial value:

{this, "SctSummaryTool",

"InDetSCT_ConditionsSummaryTool/SCT_ConditionsSummaryTool", "Tool to retrieve SCT Conditions summary"}

Definition at line 147 of file SiCombinatorialTrackFinder_xk.h.

147 {this, "SctSummaryTool",

148 "InDetSCT_ConditionsSummaryTool/SCT_ConditionsSummaryTool", "Tool to retrieve SCT Conditions summary"};

◆ m_sctcontainerkey

SG::ReadHandleKey<InDet::SCT_ClusterContainer> InDet::SiCombinatorialTrackFinder_xk::m_sctcontainerkey

private

Initial value:

{this, "SCT_ClusterContainer",

"SCT_Clusters"}

Definition at line 164 of file SiCombinatorialTrackFinder_xk.h.

164 {this, "SCT_ClusterContainer",

165 "SCT_Clusters"};

◆ m_sctDetElStatus

SG::ReadHandleKey<InDet::SiDetectorElementStatus> InDet::SiCombinatorialTrackFinder_xk::m_sctDetElStatus {this, "SCTDetElStatus", "", "Key of SiDetectorElementStatus for SCT"}

private

Optional read handle to get status data to test whether a SCT detector element is good.

If set to e.g. SCTDetectorElementStatus the event data will be used instead of the SCT conditions summary tool.

Definition at line 183 of file SiCombinatorialTrackFinder_xk.h.

184{this, "SCTDetElStatus", "", "Key of SiDetectorElementStatus for SCT"};

◆ m_updatortool

PublicToolHandle<Trk::IPatternParametersUpdator> InDet::SiCombinatorialTrackFinder_xk::m_updatortool

private

Initial value:

{this, "UpdatorTool",

"Trk::KalmanUpdator_xk/InDetPatternUpdator"}

Definition at line 154 of file SiCombinatorialTrackFinder_xk.h.

154 {this, "UpdatorTool",

155 "Trk::KalmanUpdator_xk/InDetPatternUpdator"};

◆ m_usePIX

BooleanProperty InDet::SiCombinatorialTrackFinder_xk::m_usePIX {this, "usePixel", true}

private

Definition at line 188 of file SiCombinatorialTrackFinder_xk.h.

188{this, "usePixel", true};

◆ m_useSCT

BooleanProperty InDet::SiCombinatorialTrackFinder_xk::m_useSCT {this, "useSCT", true}

private

Definition at line 189 of file SiCombinatorialTrackFinder_xk.h.

189{this, "useSCT", true};

◆ m_writeHolesFromPattern

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

private

Definition at line 196 of file SiCombinatorialTrackFinder_xk.h.

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

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

Public Member Functions

Private Attributes

Friends

Data members, which are updated in only initialize

Detailed Description

Member Enumeration Documentation

◆ EStat_t

Constructor & Destructor Documentation

◆ SiCombinatorialTrackFinder_xk()

◆ ~SiCombinatorialTrackFinder_xk()

Member Function Documentation

◆ convertToNextTrack()

◆ convertToTrack()

◆ detectorElementLinks()

◆ dump()

◆ dumpconditions()

◆ dumpevent()

◆ endEvent()

◆ fillStatistic()

◆ finalize()

◆ findTrack()

◆ getTrackQualityCuts()

◆ getTracks() [1/2]

◆ getTracks() [2/2]

◆ getTracksWithBrem()

◆ initialize()

◆ initializeCombinatorialData()

◆ magneticFieldInit()

◆ newEvent() [1/2]

◆ newEvent() [2/2]

◆ pTseed()

◆ spacePointsToClusters()

Friends And Related Symbol Documentation

◆ SiTrajectory_xk

Member Data Documentation

◆ m_boundaryCheckTool

◆ m_boundaryPixelKey

◆ m_boundarySCTKey

◆ m_doFastTracking

◆ m_fieldCondObjInputKey

◆ m_fieldmode

◆ m_fieldprop

◆ m_ITkGeometry

◆ m_minPt2Cut

◆ m_minPtCut

◆ m_outputlevel

◆ m_pixcontainerkey

◆ m_pixelCondSummaryTool

◆ m_pixelDetElStatus

◆ m_proptool

◆ m_qualityCut

◆ m_riocreator

◆ m_sctCondSummaryTool

◆ m_sctcontainerkey

◆ m_sctDetElStatus

◆ m_updatortool

◆ m_usePIX

◆ m_useSCT

◆ m_writeHolesFromPattern