InDet::InDetDetailedTrackSelectorTool Class Reference

#include <InDetDetailedTrackSelectorTool.h>

Inheritance diagram for InDet::InDetDetailedTrackSelectorTool:

Collaboration diagram for InDet::InDetDetailedTrackSelectorTool:

Public Types
enum	Grade { Undefined , Good , Shared , nbGrades }

Public Member Functions
StatusCode	initialize ()
StatusCode	finalize ()
	InDetDetailedTrackSelectorTool (const std::string &t, const std::string &n, const IInterface *p)
	~InDetDetailedTrackSelectorTool ()
bool	decision (const Trk::Track &track, const Trk::Vertex *vertex) const
bool	decision (const Trk::TrackParticleBase &track, const Trk::Vertex *vertex) const
bool	decision (const xAOD::TrackParticle &track, const xAOD::Vertex *vertex) const
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
int	getCount (const xAOD::TrackParticle &tp, xAOD::SummaryType type) const
bool	decision (const Trk::Perigee *track, const AmgSymMatrix(3) *covariancePosition) const
bool	decision (const Trk::FitQuality *TrkQuality) const
bool	decision (double chi2, int ndf) const
bool	decision (const Trk::TrackSummary *summary, const xAOD::TrackParticle *tp, bool useSharedHitInfo, bool useTrtHitInfo, const Trk::Perigee *track, const int nHitTrt, const int nHitTrtPlusOutliers) const
bool	preselectionBeforeExtrapolation (const Trk::Perigee &myPerigee) const
Amg::Vector3D	getPosOrBeamSpot (const xAOD::Vertex *) const
Trk::Vertex *	getBeamSpot (const EventContext &) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
DoubleProperty	m_pTMin {this, "pTMin", 1.*CLHEP::GeV, "min. pT: |pT|>pTMin"}
DoubleProperty	m_pMin {this, "pMin", 0., "min. p = pT/cos(theta): |p| > pMin"}
DoubleProperty	m_IPd0Max {this, "IPd0Max", 2.*CLHEP::mm, "max. d0: |d0|<d0Max"}
DoubleProperty	m_IPz0Max {this, "IPz0Max", 1.5*CLHEP::mm, "max. z0: |z0*sin(theta)|<z0Max"}
DoubleProperty	m_z0Max {this, "z0Max", 9999.*CLHEP::mm, "max. z0: |z0|<z0Max"}
DoubleProperty	m_sigIPd0Max {this, "sigIPd0Max", 999.*CLHEP::mm, "max d0 error"}
DoubleProperty	m_sigIPz0Max {this, "sigIPz0Max", 999.*CLHEP::mm, "max (error only due to z0)*sin(theta)"}
DoubleProperty	m_d0significanceMax {this, "d0significanceMax", -1., "max IP significance d0 (-1 switches it off)"}
DoubleProperty	m_z0significanceMax {this, "z0significanceMax", -1., "max IP significance z0 (-1 switches it off)"}
DoubleProperty	m_etaMax {this, "etaMax", 9999., "max. pseudo-rapidity"}
BooleanProperty	m_useTrackSummaryInfo {this, "useTrackSummaryInfo", true}
IntegerProperty	m_nHitBLayer {this, "nHitBLayer", 1, "at least n hits in Blayer"}
IntegerProperty	m_nHitPix {this, "nHitPix", 2, "at least n hits in pixels"}
IntegerProperty	m_nHitSct {this, "nHitSct", 0, "at least n hits in SCT"}
IntegerProperty	m_nHitSi {this, "nHitSi", 7, "at least n hits in pixels+SCT"}
IntegerProperty	m_nHitPixPhysical {this, "nHitPixPhysical", 0, "at least n physical hits in pixel"}
IntegerProperty	m_nHitSiPhysical {this, "nHitSiPhysical", 3, "at least n physical hits in pixel+SCT"}
IntegerProperty	m_nHitTrt {this, "nHitTrt", 0, "at least n hits in TRT"}
IntegerProperty	m_nHitTrtPlusOutliers
IntegerProperty	m_nHitTrtHighE {this, "nHitTrtHighE", 0, "at least n high threshold hits in TRT"}
IntegerProperty	m_nHitTrtPlusOutliersHighE
DoubleProperty	m_nHitTrtHighEFraction {this, "nHitTrtHighEFractionMax", 999., "maximum x fraction of transition hits in TRT"}
DoubleProperty	m_nHitTrtHighEFractionWithOutliers
DoubleProperty	m_TrtMaxEtaAcceptance
BooleanProperty	m_useSharedHitInfo {this, "useSharedHitInfo", false}
IntegerProperty	m_nSharedBLayer {this, "nSharedBLayer", 0, "max. number of shared hits in B layer"}
IntegerProperty	m_nSharedPix {this, "nSharedPix", 0, "max. number of shared hits in pixels"}
IntegerProperty	m_nSharedSct {this, "nSharedSct", 1, "max. number of shared hits in SCT"}
IntegerProperty	m_nSharedSi {this, "nSharedSi", 999, "max. number of shared hits in pixels+SCT"}
IntegerProperty	m_nHoles {this, "nHoles", 999, "max. number of holes in pixel+SCT"}
IntegerProperty	m_nDoubleHoles {this, "nDoubleHoles", 999, "max number of double-holes in SCT"}
IntegerProperty	m_nHolesPix {this, "nHolesPixel", 999, "max. number of holes in pixels"}
IntegerProperty	m_nHolesSct {this, "nHolesSct", 999, "max. number of holes in SCT"}
BooleanProperty	m_useTrackQualityInfo {this, "useTrackQualityInfo", true}
DoubleProperty	m_fitChi2 {this, "fitChi2", 99999., "max. fit chi2"}
DoubleProperty	m_fitProb {this, "fitProb", -1., "min. fit chi2 probability"}
DoubleProperty	m_fitChi2OnNdfMax {this, "fitChi2OnNdfMax", 999., "max. fitchi2/ndf"}
DoubleProperty	m_scaleMinHitTrt
IntegerProperty	m_addToMinHitTrt
DoubleProperty	m_scaleMinHitTrtWithOutliers
IntegerProperty	m_addToMinHitTrtWithOutliers
BooleanProperty	m_usePreselectionCuts {this, "usePreselectionCuts", false}
DoubleProperty	m_d0MaxPreselection {this, "d0MaxPreselection", 10.}
BooleanProperty	m_useEtaDepententMinHitTrt {this, "useEtaDepententMinHitTrt", false}
BooleanProperty	m_useEtaDepententMinHitTrtWithOutliers {this, "useEtaDepententMinHitTrtWithOutliers", false}
ToolHandle< Trk::ITrackSummaryTool >	m_trackSumTool {this, "TrackSummaryTool", "Trk::TrackSummaryTool"}
ToolHandle< Trk::ITrackParticleCreatorTool >	m_particleCreator {this, "TrackParticleCreatorTool", ""}
ToolHandle< Trk::IExtrapolator >	m_extrapolator {this, "Extrapolator", "Trk::Extrapolator"}
SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"}
ToolHandle< ITrtDriftCircleCutTool >	m_trtDCTool
ToolHandle< InDet::IInDetTestPixelLayerTool >	m_inDetTestPixelLayerTool
SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCacheCondObjInputKey
bool	m_trackSumToolAvailable = true
bool	m_partCreatorToolAvailable = false
BooleanProperty	m_usePtDependentCuts {this, "UsePtDependentCuts", false}
FloatArrayProperty	m_ptBenchmarks {this, "PtBenchmarks", {}}
IntegerArrayProperty	m_nSCTValues {this, "SCTCutValues", {}}
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfo_key {this, "EventInfo", "EventInfo", "Input event information"}
BooleanProperty	m_useEventInfoBs {this, "UseEventInfoBS", false, "Access beamspot via the EvenInfo object"}
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 60 of file InDetDetailedTrackSelectorTool.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

Grade

enum InDet::InDetDetailedTrackSelectorTool::Grade

Enumerator
Undefined
Good
Shared
nbGrades

Definition at line 65 of file InDetDetailedTrackSelectorTool.h.

{ Undefined, Good, Shared, nbGrades };

Constructor & Destructor Documentation

InDetDetailedTrackSelectorTool()

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

Definition at line 32 of file InDetDetailedTrackSelectorTool.cxx.

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

~InDetDetailedTrackSelectorTool()

InDet::InDetDetailedTrackSelectorTool::~InDetDetailedTrackSelectorTool ( )

default

Member Function Documentation

decision() [1/7]

bool InDet::InDetDetailedTrackSelectorTool::decision ( const Trk::FitQuality * TrkQuality ) const

private

Definition at line 801 of file InDetDetailedTrackSelectorTool.cxx.

  {
    if(nullptr == trkQuality) {
      ATH_MSG_WARNING( "Null FitQuality pointer passed. No track Quality cut possible. Selection failed." );
      return false;
    }
    return decision(trkQuality->chiSquared(),trkQuality->numberDoF());
  }

decision() [2/7]

bool InDet::InDetDetailedTrackSelectorTool::decision ( const Trk::Perigee * track, const AmgSymMatrix(3) * covariancePosition ) const

private

Definition at line 660 of file InDetDetailedTrackSelectorTool.cxx.

                                                                                                                       {
 
    // checking pointer first
    if(nullptr==track || !track->covariance()) {
      ATH_MSG_WARNING( "Decision on measured perigee: Zero pointer to measured perigee passed. Selection failed." );
      return false;
    }
 
    const AmgVector(5)& perigeeParms = track->parameters();
 
    // only check pt if mag. field is on
    const EventContext& ctx = Gaudi::Hive::currentContext();
    SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCacheCondObjInputKey, ctx};
    const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
    if (fieldCondObj == nullptr) {
      ATH_MSG_ERROR("execute: Failed to retrieve AtlasFieldCacheCondObj with key " << m_fieldCacheCondObjInputKey.key());
      return false;
    }
    MagField::AtlasFieldCache fieldCache;
    fieldCondObj->getInitializedCache (fieldCache);
 
    if (fieldCache.solenoidOn()){//B field
      if (perigeeParms[Trk::qOverP] == 0.) {
    ATH_MSG_DEBUG("Track rejected because of qOverP == 0.");
    return false;
      }
      double p = std::fabs(1./perigeeParms[Trk::qOverP]);
      if (p<m_pMin) {
    ATH_MSG_DEBUG("Track rejected because of p " << p << " < " << m_pMin);
    return false;
      }
      double pt = p*std::sin(perigeeParms[Trk::theta]);
      if (pt<m_pTMin) {
    ATH_MSG_DEBUG("Track rejected because of pt " << pt << " < " << m_pTMin);
    return false;
      }
    }
 
    if (std::fabs(perigeeParms[Trk::d0]) > m_IPd0Max) {
      ATH_MSG_DEBUG("Track rejected because of fabs(d0) " << std::fabs(perigeeParms[Trk::d0]) << " > " << m_IPd0Max);
      return false;
    }
 
    if (std::fabs(perigeeParms[Trk::z0]*std::sin(perigeeParms[Trk::theta])) > m_IPz0Max) {
      ATH_MSG_DEBUG("Track rejected because of fabs(z0*sin(theta)) " << std::fabs(perigeeParms[Trk::z0]*std::sin(perigeeParms[Trk::theta])) << " > " << m_IPz0Max);
      return false;
    }
 
    if (std::fabs(perigeeParms[Trk::z0]) > m_z0Max) {
      ATH_MSG_DEBUG("Track rejected because of fabs(z0) " << std::fabs(perigeeParms[Trk::z0]) << " > " << m_z0Max);
      return false;
    }
 
    if (sqrt( (*track->covariance())(Trk::z0,Trk::z0) )*std::sin(perigeeParms[Trk::theta])>m_sigIPz0Max) {
      ATH_MSG_DEBUG("Track rejected because of err(z0)*sin(theta) " << sqrt( (*track->covariance())(Trk::z0,Trk::z0) )*std::sin(perigeeParms[Trk::theta]) << " > " << m_sigIPz0Max);
      return false;
    }
 
    if (sqrt( (*track->covariance())(Trk::d0,Trk::d0) )>m_sigIPd0Max) {
      ATH_MSG_DEBUG("Track rejected because of err(d0) " << sqrt( (*track->covariance())(Trk::d0,Trk::d0) ) << " > " << m_sigIPd0Max);
      return false;
    }
 
    if (m_d0significanceMax>0 || m_z0significanceMax>0) {
 
      double sinTheta = std::sin(perigeeParms[Trk::theta]);
      double cosTheta = std::cos(perigeeParms[Trk::theta]);
      double d0wrtPriVtx = perigeeParms[Trk::d0];
      double deltaZ = perigeeParms[Trk::z0];
      double z0wrtPriVtx = deltaZ*sinTheta;
      double testtrackSigD0 = sqrt( (*track->covariance())(Trk::d0,Trk::d0) );
      double testtrackSigZ0 = sqrt( (*track->covariance())(Trk::z0,Trk::z0) );
      double testtrackSigTh = sqrt( (*track->covariance())(Trk::theta,Trk::theta) );
      // error on IP:
      double trackPhi = perigeeParms[Trk::phi];
      double dIPdx = std::sin(trackPhi);
      double dIPdy = -std::cos(trackPhi);
      double DD0 = testtrackSigD0*testtrackSigD0;
      double newD0Err=0;
      if (covariancePosition) {
        double DXX = dIPdx*dIPdx* (*covariancePosition)(0,0);
        double DYY = dIPdy*dIPdy* (*covariancePosition)(1,1);
        double DXY = 2.*dIPdx*dIPdy* (*covariancePosition)(0,1);
        newD0Err = DD0 + DXX + DYY + DXY;
      } else {
    newD0Err = DD0;
      }
 
      double d0ErrwrtPriVtx = (newD0Err>0 ? sqrt(newD0Err) : -10e-9);
 
      if (d0ErrwrtPriVtx<0) {
    ATH_MSG_WARNING( " error on d0 is negative: numeric error... (not expected. please report!)" );
      }
 
      if (m_d0significanceMax>0) {
    if (std::fabs(d0wrtPriVtx/d0ErrwrtPriVtx)>m_d0significanceMax) {
      ATH_MSG_DEBUG("Track rejected because of fabs(d0wrtPriVtx/d0ErrwrtPriVtx) " << std::fabs(d0wrtPriVtx/d0ErrwrtPriVtx) << " > " << m_d0significanceMax);
      return false;
    }
      }
 
      if (m_z0significanceMax>0) {
 
    // error on zIP:
    double dZIPdTheta = deltaZ*cosTheta;
    double dZIPdz0 = sinTheta;
    double dZIPdzV = -sinTheta;
    double DTheta2 = dZIPdTheta*dZIPdTheta*testtrackSigTh*testtrackSigTh;
    double DZ02 = dZIPdz0*dZIPdz0*testtrackSigZ0*testtrackSigZ0;
    double DThetaZ0 = 2.*dZIPdTheta*dZIPdz0*(*track->covariance())(Trk::theta,Trk::z0);
    double newZ0Err(0);
    if (covariancePosition) {
      double DZV2 = dZIPdzV*dZIPdzV* (*covariancePosition)(2,2);
      newZ0Err = DTheta2 + DZ02 + DZV2 + DThetaZ0;
    } else {
      newZ0Err = DTheta2 + DZ02 + DThetaZ0;
    }
 
    double z0ErrwrtPriVtx = (newZ0Err>0 ? sqrt(newZ0Err) : -10e-9);
 
    if (z0ErrwrtPriVtx<0) {
      ATH_MSG_WARNING( " error on z0 is negative: numeric error... (not expected. please report!)" );
    }
 
    if (std::fabs(z0wrtPriVtx/z0ErrwrtPriVtx)>m_z0significanceMax) {
      ATH_MSG_DEBUG("Track rejected because of fabs(z0wrtPriVtx/z0ErrwrtPriVtx) " << std::fabs(z0wrtPriVtx/z0ErrwrtPriVtx) << " > " << m_z0significanceMax);
      return false;
    }
      }
 
    }
 
    if (std::fabs(track->momentum().eta())>m_etaMax) {
      ATH_MSG_DEBUG("Track rejected because of fabs(eta) " << std::fabs(track->momentum().eta()) << " > " << m_etaMax);
      return false;
    }
 
    return true;
  }

decision() [3/7]

bool InDet::InDetDetailedTrackSelectorTool::decision ( const Trk::Track & track, const Trk::Vertex * vertex ) const

virtual

Implements Trk::ITrackSelectorTool.

Definition at line 129 of file InDetDetailedTrackSelectorTool.cxx.

                                                                                             {
 
    const EventContext& ctx = Gaudi::Hive::currentContext();
 
    int nHitTrt = m_nHitTrt;
    int nHitTrtPlusOutliers = m_nHitTrtPlusOutliers;
    const Trk::Perigee* perigeeBeforeExtrapolation=dynamic_cast<const Trk::Perigee*>(track.perigeeParameters());
    if (perigeeBeforeExtrapolation && m_usePreselectionCuts){
        bool preselectionDecision=preselectionBeforeExtrapolation(*perigeeBeforeExtrapolation);
        if (!preselectionDecision) {
          ATH_MSG_DEBUG("Track rejected because of preselection decision!");
          return false;
        }
    } else if (m_usePreselectionCuts){
      ATH_MSG_DEBUG( " Preselection was requested but cannot be made since no Perigee in Track is available. This is not an error." );
    }
    const Trk::Vertex* myVertex=vertex;
    //in case no Vertex is provided by the user, beam position will be used if available
    if (myVertex==nullptr) {
      myVertex = getBeamSpot(ctx);
    }
    Trk::PerigeeSurface perigeeSurface(myVertex->position());
    const Trk::TrackParameters *firstmeaspar=nullptr;
    for (const auto *i : *track.trackParameters()){
      if ( i->covariance() && !dynamic_cast<const Trk::Perigee*>(i)) {
        firstmeaspar=i;
        break;
      }
    }
    if (!firstmeaspar) {
      //assumes perigeeParameters exist...
      //no track selection if firstmeas + perigee does not exist !
      firstmeaspar=track.perigeeParameters();
      if (!firstmeaspar){
          ATH_MSG_WARNING( " First measurment on track is missing. Using perigee Parameters, but they are missing: 0 pointer! Track selection failed " );
          //clean up vertex
          if (myVertex!=vertex) {
            delete myVertex;
            myVertex=nullptr;
          }
          return false;
        }
    }
    const Trk::TrackParameters* extrapolatedParameters= m_extrapolator->extrapolate(ctx,
                                                                                    *firstmeaspar,
                                                                                    perigeeSurface,
                                                                                    Trk::anyDirection,
                                                                                    true,
                                                                                    track.info().particleHypothesis() ).release();
    const Trk::Perigee* extrapolatedPerigee = extrapolatedParameters ? dynamic_cast<const Trk::Perigee*>(extrapolatedParameters) : nullptr;
    if (!extrapolatedPerigee || !extrapolatedPerigee->covariance() ) {
      ATH_MSG_WARNING( "Track Selector failed to extrapolate track to the vertex: " << myVertex->position() );
      if (extrapolatedParameters) {
        ATH_MSG_WARNING( "The return object of the extrapolator was not a perigee even if a perigeeSurface was used!" );
        delete extrapolatedParameters;
        extrapolatedParameters=nullptr;
      }
    }
 
    //decision based on the track parameters
    const Trk::RecVertex* recVertex = dynamic_cast<const Trk::RecVertex*>(myVertex);
    bool dec = decision(extrapolatedPerigee, recVertex ? &recVertex->covariancePosition() : nullptr );
    if (myVertex!=vertex) {
      delete myVertex;
      myVertex=nullptr;
    }
    bool isInTrtAcceptance=true;
    if (!extrapolatedPerigee || std::fabs(extrapolatedPerigee->momentum().eta())>m_TrtMaxEtaAcceptance) {
      isInTrtAcceptance=false;
    }
    if (extrapolatedPerigee!=track.perigeeParameters()) {
      delete extrapolatedPerigee;
      extrapolatedPerigee=nullptr;
    }
    if(!dec) {
      ATH_MSG_DEBUG("Track rejected because of perigee parameters!");
      return false;
    }
    if (m_useTrackQualityInfo) {
      const Trk::FitQuality*  TrkQuality=track.fitQuality();
      if (TrkQuality==nullptr) {
        ATH_MSG_WARNING( "Requested cut on track quality was not possible. Track has no FitQuality object attached. Selection failed." );
        return false;
      }
      if (!decision(TrkQuality)) {
        return false;
      }
    }
    if (m_useTrackSummaryInfo) {
      //number of hits, silicon hits, b-layer
      // first ask track for summary
      std::unique_ptr<Trk::TrackSummary> summaryUniquePtr;
      const Trk::TrackSummary* summary = track.trackSummary();
      if (m_trackSumToolAvailable && summary == nullptr) {
        summaryUniquePtr = m_trackSumTool->summary(ctx, track);
        summary = summaryUniquePtr.get();
      }
      if (nullptr==summary ) {
        ATH_MSG_FATAL( "Track preselection: cannot create a track summary (but useTrackSummary is true). Selection failed." );
        return false;
      }
 
      // Create xAOD::TrackParticle to retrieve shared hit info
      const xAOD::TrackParticle* tp = m_partCreatorToolAvailable ? m_particleCreator->createParticle(ctx, track) : nullptr;
      if(m_useSharedHitInfo && tp==nullptr){
        ATH_MSG_FATAL( "Track preselection: cannot create a track particle (but useSharedHitInfo is true). Selection failed." );
        return false;
      }
 
      // get the minimum nimber of TRT hits based on eta of the track
      if(m_useEtaDepententMinHitTrt) {
        nHitTrt = m_trtDCTool->minNumberDCs( (*track.trackParameters())[0] );
        if(m_addToMinHitTrt!=0){
          nHitTrt += m_addToMinHitTrt;
        }else{
          nHitTrt = (int)((double)nHitTrt*m_scaleMinHitTrt);
        }
      }
 
      // get the minimum nimber of TRT hits + outliers based on eta of the track
      if(m_useEtaDepententMinHitTrtWithOutliers) {
        nHitTrtPlusOutliers = m_trtDCTool->minNumberDCs( (*track.trackParameters())[0] );
        if(m_addToMinHitTrtWithOutliers!=0){
          nHitTrtPlusOutliers += m_addToMinHitTrtWithOutliers;
        }else{
          nHitTrtPlusOutliers = (int)((double)nHitTrtPlusOutliers*m_scaleMinHitTrtWithOutliers);
        }
      }
 
      if (!decision(summary, tp, m_useSharedHitInfo,isInTrtAcceptance, perigeeBeforeExtrapolation,
                    nHitTrt, nHitTrtPlusOutliers)) {
        return false;
      }
    }
    return true;
  }

decision() [4/7]

bool InDet::InDetDetailedTrackSelectorTool::decision ( const Trk::TrackParticleBase & track, const Trk::Vertex * vertex ) const

virtual

Implements Trk::ITrackSelectorTool.

Definition at line 268 of file InDetDetailedTrackSelectorTool.cxx.

                                                                                                         {
 
    const EventContext& ctx = Gaudi::Hive::currentContext();
 
    int nHitTrt = m_nHitTrt;
    int nHitTrtPlusOutliers = m_nHitTrtPlusOutliers;
    const Trk::TrackParameters* definintParameters=&(track.definingParameters());
    const Trk::Perigee* perigeeBeforeExtrapolation=dynamic_cast<const Trk::Perigee*>(definintParameters);
    if (perigeeBeforeExtrapolation && m_usePreselectionCuts) {
      bool preselectionDecision=preselectionBeforeExtrapolation(*perigeeBeforeExtrapolation);
      if (!preselectionDecision) {
        ATH_MSG_DEBUG("Track rejected because of preselection decision!");
        return false;
      }
    } else if (m_usePreselectionCuts) {
      ATH_MSG_WARNING( " Preselection was requested but cannot be made since the Perigee is not the defining Parameter of the TrackParticle. This is not an error." );
    }
    bool isInTrtAcceptance=true;
    if (!perigeeBeforeExtrapolation || std::fabs(perigeeBeforeExtrapolation->momentum().eta())>m_TrtMaxEtaAcceptance) {
      isInTrtAcceptance=false;
    }
    if (m_useTrackQualityInfo) {
      const Trk::FitQuality*  TrkQuality=track.fitQuality();
      if (TrkQuality==nullptr) {
        ATH_MSG_WARNING( "Requested cut on track quality was not possible. TrackParticleBase has no FitQuality object attached. Selection failed." );
        return false;
      }
      if (!decision(TrkQuality)) {
        return false;
      }
    }
    if (m_useTrackSummaryInfo) {
      //number of hits, silicon hits, b-layer
      const Trk::TrackSummary* summary = track.trackSummary();
      if (nullptr==summary ) {
        ATH_MSG_WARNING( "Track preselection: cannot create a track summary (but useTrackSummary is true). Selection failed." );
        return false;
      }
 
 
      if (m_useSharedHitInfo) {
    ATH_MSG_ERROR( "Use of InDetDetailedTrackSelectorTool with Trk::TrackParticleBase and useSharedHitInfo is not supported");
    return false;
      }
      const xAOD::TrackParticle* tp = nullptr;
 
      if(m_useEtaDepententMinHitTrt) {
        nHitTrt = m_trtDCTool->minNumberDCs( (track.trackParameters())[0] );
        if(m_addToMinHitTrt!=0){
          nHitTrt += m_addToMinHitTrt;
        }else{
          nHitTrt = (int)((double)nHitTrt*m_scaleMinHitTrt);
        }
      }
 
      if(m_useEtaDepententMinHitTrtWithOutliers) {
        nHitTrtPlusOutliers = m_trtDCTool->minNumberDCs( (track.trackParameters())[0] );
        if(m_addToMinHitTrtWithOutliers!=0){
          nHitTrtPlusOutliers += m_addToMinHitTrtWithOutliers;
        }else{
          nHitTrtPlusOutliers = (int)((double)nHitTrtPlusOutliers*m_scaleMinHitTrtWithOutliers);
        }
      }
 
      if ((!perigeeBeforeExtrapolation) or
          (!decision(summary, tp, m_useSharedHitInfo, isInTrtAcceptance, perigeeBeforeExtrapolation,
                     nHitTrt, nHitTrtPlusOutliers))) {
          return false;
      }
    }
    const Trk::Perigee* extrapolatedPerigee=dynamic_cast<const Trk::Perigee*>(definintParameters);
    const Trk::Vertex* myVertex=vertex;
    if (vertex==nullptr) {
      myVertex = getBeamSpot(ctx);
    }
    Trk::PerigeeSurface perigeeSurface(myVertex->position());
    const Trk::TrackParameters *firstmeaspar=nullptr;
    for (const auto *i : track.trackParameters()) {
      if (i->covariance() &&
         !dynamic_cast<const Trk::Perigee*>(i)) {
        firstmeaspar=i;
        break;
      }
    }
    if (!firstmeaspar) {
      if (!extrapolatedPerigee || !extrapolatedPerigee->covariance() ) {
          ATH_MSG_DEBUG( " Track Paraemters at first measurement not found. Perigee not found. Cannot do TrackSelection..." );
          if (myVertex!=vertex) {
            delete myVertex;
            myVertex=nullptr;
          }
          return false;
      }
      //using perigee instead of firstmeasurement, since first measurement was not found...
      firstmeaspar=&(track.definingParameters());
    }
 
    ATH_MSG_VERBOSE ("Input to extrapolation: "    << *firstmeaspar);
    ATH_MSG_VERBOSE ("Extrapolating to position: " << myVertex->position()[0] << " , " <<
             myVertex->position()[1] << " , " << myVertex->position()[2]);
    const Trk::TrackParameters* extrapolatedParameters= firstmeaspar ?
      m_extrapolator->extrapolate(ctx,
                                  *firstmeaspar,
                                  perigeeSurface,
                                  Trk::anyDirection,
                                  true,Trk::pion ).release() : nullptr;
    extrapolatedPerigee = extrapolatedParameters ? dynamic_cast<const Trk::Perigee*>(extrapolatedParameters) : nullptr;
    if (extrapolatedPerigee==nullptr || !extrapolatedPerigee->covariance()) {
      ATH_MSG_WARNING( "Track Selector failed to extrapolate track to the vertex: " << myVertex->position() );
      if (extrapolatedParameters) {
          ATH_MSG_WARNING( "The return object of the extrapolator was not a perigee even if a perigeeSurface was used!" );
          delete extrapolatedParameters;
          extrapolatedParameters = nullptr;
      }
    }
    if (extrapolatedParameters) ATH_MSG_VERBOSE ("Result: " << *extrapolatedParameters);
    const Trk::RecVertex* recVertex = dynamic_cast<const Trk::RecVertex*>(myVertex);
    bool dec = decision(extrapolatedPerigee, recVertex ? &recVertex->covariancePosition() : nullptr );
    if (myVertex!=vertex) {
      delete myVertex;
      myVertex=nullptr;
    }
    if (extrapolatedPerigee!=&(track.definingParameters())) {
      delete extrapolatedPerigee;
      extrapolatedPerigee=nullptr;
    }
    if(!dec) {
      ATH_MSG_DEBUG("Track rejected because of perigee parameters!");
      return false;
    }
    return true;
    }

decision() [5/7]

bool InDet::InDetDetailedTrackSelectorTool::decision ( const Trk::TrackSummary * summary, const xAOD::TrackParticle * tp, bool useSharedHitInfo, bool useTrtHitInfo, const Trk::Perigee * track, const int nHitTrt, const int nHitTrtPlusOutliers ) const

private

Definition at line 842 of file InDetDetailedTrackSelectorTool.cxx.

  {
    if (summary==nullptr) {
      ATH_MSG_WARNING( "Null TrackSummary pointer passed. Selection failed." );
      return false;
    }
 
    int nb = summary->get(Trk::numberOfInnermostPixelLayerHits);
 
    if(nb<0) nb=0;
 
    int np  = summary->get(Trk::numberOfPixelHits);
    if(np<0) np=0;
 
    int npd = summary->get(Trk::numberOfPixelDeadSensors);
    if(npd<0) npd=0;
 
    int ns = summary->get(Trk::numberOfSCTHits);
    if(ns<0) ns=0;
 
    int nhp = summary->get(Trk::numberOfPixelHoles);
    if (nhp < 0) nhp = 0;
 
    int nhs = summary->get(Trk::numberOfSCTHoles);
    if (nhs < 0) nhs = 0;
 
    int ndhs = summary->get(Trk::numberOfSCTDoubleHoles);
    if (ndhs < 0) ndhs = 0;
 
    //**-----------------------------------------------------------------------
 
    if(m_usePtDependentCuts) {
      if (!track) {
         return false;
      }
      const AmgVector(5)& perigeeParms = track->parameters();
      double p = std::fabs(1./perigeeParms[Trk::qOverP]);
      double pt = p*std::sin(perigeeParms[Trk::theta]);
 
      unsigned int it = 0;
      for(; it< m_ptBenchmarks.size()-1; ++it ) {
    if(pt>m_ptBenchmarks[it] && pt <=m_ptBenchmarks[it+1] && ns < m_nSCTValues[it]) {
      ATH_MSG_DEBUG("Track rejected because of Pt-Dependent SCT Hit cut (CAREFUL! Excludes dead modules)") ;
      return false;
    }
      }//end of  pt intervals loop
 
      //now cutting all the rest by the last value in the vector
      if(pt>m_ptBenchmarks[it+1] && ns < m_nSCTValues[it+1]) {
    ATH_MSG_DEBUG("Track rejected because of Pt-Dependent SCT Hit cut (CAREFUL! Excludes dead modules)") ;
    return false;
      }
 
    }
 
    //*--------------------------------------------------------------------------------
 
    //normal cuts in all their variety
 
    if(nb == 0 && nb < m_nHitBLayer) {
      ATH_MSG_DEBUG("Track rejected because of nHitBLayer "<<nb<<" < "<<m_nHitBLayer);
      if(m_inDetTestPixelLayerTool.empty()) {
    ATH_MSG_DEBUG("and no blayer tool configured, so will not try to recover track");
    return false;
      } else if (m_inDetTestPixelLayerTool->expectHitInInnermostPixelLayer(track)) {
    ATH_MSG_DEBUG("and track rejected because at least one hit is expected in the innermost pixel layer") ;
    return false;
      }else  ATH_MSG_DEBUG("recovered track as no b-layer expected") ;
    }//end of checking the b-layer
 
    if(np+npd < m_nHitPix) {
      ATH_MSG_DEBUG("Track rejected because of nHitPix "<<np+npd<<" < "<<m_nHitPix);
      return false;
    }
 
    if(np < m_nHitPixPhysical) {
      ATH_MSG_DEBUG("Track rejected because of nHitPixPhysical "<<np<<" < "<<m_nHitPixPhysical);
      return false;
    }
 
    int nsd = summary->get(Trk::numberOfSCTDeadSensors);
    if(nsd<0)
      nsd=0;
 
    if(ns+nsd < m_nHitSct)
      {
    ATH_MSG_DEBUG("Track rejected because of nHitSct "<<ns+nsd<<" < "<<m_nHitSct);
    return false;
      }
 
    if((np+ns+npd+nsd) < m_nHitSi)
      {
    ATH_MSG_DEBUG("Track rejected because of nHitSi "<<np+npd+ns+nsd<<" < "<<m_nHitSi);
    return false;
      }
 
    if((np+ns) < m_nHitSiPhysical)
      {
    ATH_MSG_DEBUG("Track rejected because of nHitSiPhysical "<<np+ns<<" < "<<m_nHitSiPhysical);
    return false;
      }
 
    // Cuts on number of Holes
 
    if ((nhp+nhs) > m_nHoles)
      {
    ATH_MSG_DEBUG("Track rejected because of nHolesPixPlusSCT "<<nhp+nhs<<" > "<<m_nHoles);
    return false;
      }
 
    if (ndhs > m_nDoubleHoles)
      {
    ATH_MSG_DEBUG("Track rejected because of nDoubleHolesSCT "<<ndhs<<" > "<<m_nDoubleHoles);
    return false;
      }
 
    if (nhp > m_nHolesPix)
      {
    ATH_MSG_DEBUG("Track rejected because of nHolesPix "<<nhp<<" > "<<m_nHolesPix);
    return false;
      }
 
    if (nhs > m_nHolesSct)
      {
    ATH_MSG_DEBUG("Track rejected because of nHolesSct "<<nhs<<" > "<<m_nHolesSct);
    return false;
      }
 
    if (useTrtHitInfo) {
 
      int nh = summary->get(Trk::numberOfTRTHits);
      if(nh<0) nh=0;
      if(nh < nHitTrt) {
    ATH_MSG_DEBUG("Track rejected because of nHitTrt "<<nh<<" < "<<nHitTrt);
    return false;
      }
 
      int nhh = summary->get( Trk::numberOfTRTHits ) + summary->get( Trk::numberOfTRTOutliers );
      if (nhh<0) nhh=0;
      if (nhh<nHitTrtPlusOutliers) {
    ATH_MSG_DEBUG("Track rejected because of nHitTrtPlusOutliers "<<nhh<<" < "<<nHitTrtPlusOutliers);
    return false;
      }
 
      int nhthits=summary->get(Trk::numberOfTRTHighThresholdHits);
      if (nhthits<0) nhthits=0;
      if (nhthits<m_nHitTrtHighE) {
    ATH_MSG_DEBUG("Track rejected because of nHitTrtHighE "<<nhthits<<" < "<<m_nHitTrtHighE);
    return false;
      }
 
      int nhthitsWithOutliers=summary->get(Trk::numberOfTRTHighThresholdHits) + summary->get(Trk::numberOfTRTHighThresholdOutliers);
      if (nhthitsWithOutliers<0) nhthitsWithOutliers=0;
      if (nhthitsWithOutliers<m_nHitTrtPlusOutliersHighE) {
    ATH_MSG_DEBUG("Track rejected because of nHitTrtPlusOutliersHighE "<<nhthitsWithOutliers<<" < "<<m_nHitTrtPlusOutliersHighE);
    return false;
      }
 
      if (summary->get( Trk :: numberOfTRTHits )>0) {
    double nhe = (double)summary->get(Trk::numberOfTRTHighThresholdHits) / (double)summary->get( Trk::numberOfTRTHits );
    if(nhe<0.) nhe=0.;
    if(nhe > m_nHitTrtHighEFraction ) {
      ATH_MSG_DEBUG("Track rejected because of nHitTrtHighEFraction "<<nhe<<" < "<<m_nHitTrtHighEFraction);
      return false;
    }
      }
 
      if ( summary->get( Trk :: numberOfTRTHits ) + summary->get( Trk :: numberOfTRTOutliers ) > 0 ) {
    double nheh = (double)(summary->get(Trk::numberOfTRTHighThresholdHits) + summary->get(Trk::numberOfTRTHighThresholdOutliers))/
      (double)(summary->get( Trk::numberOfTRTHits) + summary->get( Trk :: numberOfTRTOutliers ) );
    if(nheh<0.) nheh=0.;
    if (nheh>1.) nheh=1.;
    if(nheh > m_nHitTrtHighEFractionWithOutliers ) {
      ATH_MSG_DEBUG("Track rejected because of nHitTrtHighEFractionWithOutliers "<<nheh<<" < "<<m_nHitTrtHighEFractionWithOutliers);
      return false;
    }
      }
    }
 
    if (useSharedHitInfo) {
      if(!tp){
    ATH_MSG_DEBUG("Track rejected because xAOD::TrackParticle not available");
    return false;
      }
 
      int nbs = getCount(*tp,xAOD::numberOfInnermostPixelLayerSharedHits);
      if(nbs < 0) nbs = 0;
      if (nbs>1) nbs=1;
      if(nbs>m_nSharedBLayer) {
    ATH_MSG_DEBUG("Track rejected because of nSharedBLayer "<<nbs<<" < "<<m_nSharedBLayer);
    return false;
      }
 
      int nps = getCount(*tp,xAOD::numberOfPixelSharedHits);
      if(nps < 0) nps = 0;
      if(nps>m_nSharedPix) {
    ATH_MSG_DEBUG("Track rejected because of nSharedPix "<<nps<<" < "<<m_nSharedPix);
    return false;
      }
 
      int nss = getCount(*tp,xAOD::numberOfSCTSharedHits);
      if(nss < 0) nss = 0;
      if(nss > m_nSharedSct) {
    ATH_MSG_DEBUG("Track rejected because of nSharedSct "<<nss<<" < "<<m_nSharedSct);
    return false;
      }
 
      int nst = nps + nss;
      if(nst>m_nSharedSi) {
    ATH_MSG_DEBUG("Track rejected because of nSharedSi "<<nst<<" < "<<m_nSharedSi);
    return false;
      }
    }
 
    return true;
 
  }

decision() [6/7]

bool InDet::InDetDetailedTrackSelectorTool::decision ( const xAOD::TrackParticle & track, const xAOD::Vertex * vertex ) const

virtual

Implements Trk::ITrackSelectorTool.

Definition at line 428 of file InDetDetailedTrackSelectorTool.cxx.

  {
    const EventContext& ctx = Gaudi::Hive::currentContext();
 
    int nHitTrt = m_nHitTrt;
    int nHitTrtPlusOutliers = m_nHitTrtPlusOutliers;
 
    const Trk::Perigee& perigee=tp.perigeeParameters();
    if (m_usePreselectionCuts && !preselectionBeforeExtrapolation(perigee)) {
      ATH_MSG_DEBUG("Track rejected because of preselection decision!");
      return false;
    }
 
    if (m_useTrackQualityInfo && !decision(tp.chiSquared(),tp.numberDoF())) {
      ATH_MSG_DEBUG("Track rejected because of bad fit quality!");
      return false;
    }
 
    if (m_useTrackSummaryInfo) {
      //number of hits, silicon hits, b-layer
 
      if(m_useEtaDepententMinHitTrt) {
        nHitTrt = m_trtDCTool->minNumberDCs( &perigee );
        if(m_addToMinHitTrt!=0){
          nHitTrt += m_addToMinHitTrt;
        }else{
          nHitTrt = (int)((double)nHitTrt*m_scaleMinHitTrt);
        }
      }
      if(m_useEtaDepententMinHitTrtWithOutliers) {
        nHitTrtPlusOutliers = m_trtDCTool->minNumberDCs( &perigee );
        if(m_addToMinHitTrtWithOutliers!=0){
          nHitTrtPlusOutliers += m_addToMinHitTrtWithOutliers;
        }else{
          nHitTrtPlusOutliers = (int)((double)nHitTrtPlusOutliers*m_scaleMinHitTrtWithOutliers);
        }
      }
      int nb = getCount(tp,xAOD::numberOfInnermostPixelLayerHits );
      int np = getCount(tp,xAOD::numberOfPixelHits );
      int npd = getCount(tp,xAOD::numberOfPixelDeadSensors );
      int ns  = getCount(tp,xAOD::numberOfSCTHits );
      int nhp = getCount(tp,xAOD::numberOfPixelHoles );
      int nhs = getCount(tp,xAOD::numberOfSCTHoles );
      int ndhs = getCount(tp,xAOD::numberOfSCTDoubleHoles);
      bool eiph = (getCount(tp,xAOD::expectInnermostPixelLayerHit)==1);
 
      //**-----------------------------------------------------------------------
      if(m_usePtDependentCuts) {
          double pt = tp.pt();
          unsigned int it = 0;
        for(; it< m_ptBenchmarks.size()-1; ++it ) {
          if(pt>m_ptBenchmarks[it] && pt <=m_ptBenchmarks[it+1] && ns < m_nSCTValues[it]) {
            ATH_MSG_DEBUG("Track rejected because of Pt-Dependent SCT Hit cut (CAREFUL! Excludes dead modules)") ;
            return false;
          }
        }//end of  pt intervals loop
 
        //now cutting all the rest by the last value in the vector
        if(pt>m_ptBenchmarks[it+1] && ns < m_nSCTValues[it+1]) {
          ATH_MSG_DEBUG("Track rejected because of Pt-Dependent SCT Hit cut (CAREFUL! Excludes dead modules)") ;
          return false;
        }
      }
 
      //*--------------------------------------------------------------------------------
 
      //normal cuts in all their variety
 
      if(nb == 0 && nb < m_nHitBLayer) {
        ATH_MSG_DEBUG("Track rejected because of nHitBLayer "<<nb<<" < "<<m_nHitBLayer);
    if (eiph) {
          ATH_MSG_DEBUG("and track rejected because at least one hit is expected in the innermost pixel layer") ;
          return false;
        }else  ATH_MSG_DEBUG("recovered track as no b-layer expected") ;
      }//end of checking the b-layer
 
      if(np+npd < m_nHitPix) {
        ATH_MSG_DEBUG("Track rejected because of nHitPix "<<np+npd<<" < "<<m_nHitPix);
        return false;
      }
 
      if(np < m_nHitPixPhysical) {
        ATH_MSG_DEBUG("Track rejected because of nHitPixPhysical "<<np<<" < "<<m_nHitPixPhysical);
        return false;
      }
 
      int nsd = getCount(tp,xAOD::numberOfSCTDeadSensors);
      if(ns+nsd < m_nHitSct) {
        ATH_MSG_DEBUG("Track rejected because of nHitSct "<<ns+nsd<<" < "<<m_nHitSct);
        return false;
      }
 
      if(np+ns+npd+nsd < m_nHitSi) {
        ATH_MSG_DEBUG("Track rejected because of nHitSi "<<np+npd+ns+nsd<<" < "<<m_nHitSi);
        return false;
      }
 
      if(np+ns < m_nHitSiPhysical) {
        ATH_MSG_DEBUG("Track rejected because of nHitSiPhysical "<<np+ns<<" < "<<m_nHitSiPhysical);
        return false;
      }
 
      // Cuts on number of Holes
      if (nhp+nhs > m_nHoles){
    ATH_MSG_DEBUG("Track rejected because of nHolesPixPlusSCT "<<nhp+nhs<<" > "<<m_nHoles);
    return false;
      }
 
      if(ndhs > m_nDoubleHoles){
    ATH_MSG_DEBUG("Track rejected because of nDoubleHolesSCT "<<ndhs<<" > "<<m_nDoubleHoles);
    return false;
      }
 
      if(nhp > m_nHolesPix){
        ATH_MSG_DEBUG("Track rejected because of nHolesPix "<<nhp<<" > "<<m_nHolesPix);
        return false;
      }
 
      if (nhs > m_nHolesSct){
        ATH_MSG_DEBUG("Track rejected because of nHolesSct "<<nhs<<" > "<<m_nHolesSct);
        return false;
      }
 
      if (std::fabs(tp.eta())>m_TrtMaxEtaAcceptance) {
        int nh = getCount(tp,xAOD::numberOfTRTHits);
        if(nh < nHitTrt) {
          ATH_MSG_DEBUG("Track rejected because of nHitTrt "<<nh<<" < "<<nHitTrt);
          return false;
        }
 
      int nhh = getCount(tp, xAOD::numberOfTRTHits ) + getCount(tp, xAOD::numberOfTRTOutliers );
      if (nhh<nHitTrtPlusOutliers) {
        ATH_MSG_DEBUG("Track rejected because of nHitTrtPlusOutliers "<<nhh<<" < "<<nHitTrtPlusOutliers);
        return false;
      }
 
      int nhthits= getCount(tp,xAOD::numberOfTRTHighThresholdHits);
      if (nhthits<m_nHitTrtHighE) {
        ATH_MSG_DEBUG("Track rejected because of nHitTrtHighE "<<nhthits<<" < "<<m_nHitTrtHighE);
        return false;
      }
 
      int nhthitsWithOutliers= getCount(tp,xAOD::numberOfTRTHighThresholdHits) + getCount(tp,xAOD::numberOfTRTHighThresholdOutliers);
      if (nhthitsWithOutliers<m_nHitTrtPlusOutliersHighE) {
        ATH_MSG_DEBUG("Track rejected because of nHitTrtPlusOutliersHighE "<<nhthitsWithOutliers<<" < "<<m_nHitTrtPlusOutliersHighE);
        return false;
      }
 
      if ( getCount(tp, xAOD::numberOfTRTHits )>0) {
        double nhe = getCount(tp,xAOD::numberOfTRTHighThresholdHits);
        nhe /= getCount(tp, xAOD::numberOfTRTHits);
        if(nhe > m_nHitTrtHighEFraction ) {
          ATH_MSG_DEBUG("Track rejected because of nHitTrtHighEFraction "<<nhe<<" < "<<m_nHitTrtHighEFraction);
          return false;
        }
      }
 
      if ( getCount(tp, xAOD::numberOfTRTHits ) + getCount(tp, xAOD::numberOfTRTOutliers ) > 0 ) {
        double nheh = (double)(getCount(tp,xAOD::numberOfTRTHighThresholdHits) + getCount(tp,xAOD::numberOfTRTHighThresholdOutliers))/
                (double)(getCount(tp, xAOD::numberOfTRTHits) + getCount(tp, xAOD::numberOfTRTOutliers ) );
        if(nheh<0.) nheh=0.;
        if (nheh>1.) nheh=1.;
        if(nheh > m_nHitTrtHighEFractionWithOutliers ) {
          ATH_MSG_DEBUG("Track rejected because of nHitTrtHighEFractionWithOutliers "<<nheh<<" < "<<m_nHitTrtHighEFractionWithOutliers);
          return false;
        }
      }
      }
      if (m_useSharedHitInfo) {
        int nbs = getCount(tp,xAOD::numberOfInnermostPixelLayerSharedHits);
        if (nbs>1) nbs=1;
        if(nbs>m_nSharedBLayer) {
          ATH_MSG_DEBUG("Track rejected because of nSharedBLayer "<<nbs<<" < "<<m_nSharedBLayer);
          return false;
        }
 
    int nps = getCount(tp,xAOD::numberOfPixelSharedHits);
    if(nps>m_nSharedPix) {
      ATH_MSG_DEBUG("Track rejected because of nSharedPix "<<nps<<" < "<<m_nSharedPix);
      return false;
    }
 
    int nss = getCount(tp,xAOD::numberOfSCTSharedHits);
    if(nss > m_nSharedSct) {
      ATH_MSG_DEBUG("Track rejected because of nSharedSct "<<nss<<" < "<<m_nSharedSct);
      return false;
    }
 
    int nst = nps + nss;
    if(nst>m_nSharedSi) {
      ATH_MSG_DEBUG("Track rejected because of nSharedSi "<<nst<<" < "<<m_nSharedSi);
      return false;
    }
      }
    }
    Trk::PerigeeSurface perigeeSurface( getPosOrBeamSpot(vertex) );
 
    const Trk::TrackParameters* extrapolatedParameters= m_extrapolator->extrapolate(
      ctx,
      perigee,perigeeSurface,
      Trk::anyDirection,true,Trk::pion).release();
    const Trk::Perigee* extrapolatedPerigee = extrapolatedParameters ? dynamic_cast<const Trk::Perigee*>(extrapolatedParameters) : nullptr;
    if (extrapolatedPerigee==nullptr) {
      ATH_MSG_WARNING( "Extrapolation to the vertex failed: " << perigeeSurface << std::endl << perigee );
      if (extrapolatedParameters!=nullptr) {
        ATH_MSG_WARNING( "The return object of the extrapolator was not a perigee even if a perigeeSurface was used!" );
        delete extrapolatedParameters;
        extrapolatedParameters=nullptr;
      }
      return false;
    }
    bool dec = false;
    if( vertex ){
      // for now copy the position error
      const AmgSymMatrix(3)& vertexError = vertex->covariancePosition();
      dec = decision(extrapolatedPerigee,&vertexError);
    }else{
      dec = decision(extrapolatedPerigee,nullptr);
    }
 
    delete extrapolatedPerigee;
 
    if(!dec) {
      ATH_MSG_DEBUG("Track rejected because of perigee parameters!");
      return false;
    }
 
    return true;
  }

decision() [7/7]

bool InDet::InDetDetailedTrackSelectorTool::decision ( double chi2, int ndf ) const

private

Definition at line 810 of file InDetDetailedTrackSelectorTool.cxx.

                                                                          {
 
    double proba = 1.;
 
    if(ndf>0 && chi2>=0.) {
      Genfun::CumulativeChiSquare myCumulativeChiSquare(ndf);
      proba = 1.-myCumulativeChiSquare(chi2);
    }
 
    if(chi2>m_fitChi2) {
      ATH_MSG_DEBUG("Track rejected because of chi2 "<<chi2<<" > "<<m_fitChi2);
      return false;
    }
 
    if(proba<m_fitProb) {
      ATH_MSG_DEBUG("Track rejected because of fit probability "<<proba<<" > "<<m_fitProb);
      return false;
    }
    if(!ndf) {
      ATH_MSG_DEBUG("Track rejected because of ndof = "<<ndf);
      return false;
    }
    if(chi2/double(ndf)>m_fitChi2OnNdfMax) {
      ATH_MSG_DEBUG("Track rejected because of chi2/ndof "<<chi2/double(ndf)<<" > "<<m_fitChi2OnNdfMax);
      return false;
    }
 
    return true;
  }

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

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::InDetDetailedTrackSelectorTool::finalize

(

)

Definition at line 121 of file InDetDetailedTrackSelectorTool.cxx.

  {
    ATH_MSG_DEBUG( "Finalize successful" );
    return StatusCode::SUCCESS;
  }

getBeamSpot()

Trk::Vertex * InDet::InDetDetailedTrackSelectorTool::getBeamSpot ( const EventContext & ctx ) const

private

Definition at line 38 of file InDetDetailedTrackSelectorTool.cxx.

  {
    if(m_useEventInfoBs){
      SG::ReadHandle<xAOD::EventInfo> evt(m_eventInfo_key, ctx);
      if (evt.isValid()) {
        InDet::BeamSpotData temp(evt->beamStatus(), evt->beamPosX(), evt->beamPosY(), evt->beamPosZ(),
                                 evt->beamPosSigmaX(), evt->beamPosSigmaY(), evt->beamPosSigmaZ(),
                                 evt->beamTiltXZ(), evt->beamTiltYZ(), evt->beamPosSigmaXY());
        return new Trk::RecVertex(temp.beamVtx());
      } else {
        ATH_MSG_WARNING( " Cannot get beamSpot center from xAOD::EventInfo. Using (0,0,0)... " );
        return new Trk::Vertex(Amg::Vector3D(0,0,0));
      }
    }else{
      SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey, ctx };
      if (beamSpotHandle.isValid()) {
        return new Trk::RecVertex(beamSpotHandle->beamVtx());
      } else {
        ATH_MSG_WARNING( " Cannot get beamSpot center from BeamSpotData. Using (0,0,0)... " );
        return new Trk::Vertex(Amg::Vector3D(0,0,0));
      }
    }
  }

getCount()

int InDet::InDetDetailedTrackSelectorTool::getCount ( const xAOD::TrackParticle & tp, xAOD::SummaryType type ) const

inlineprivate

Definition at line 82 of file InDetDetailedTrackSelectorTool.h.

                                                                            {
      uint8_t val;
      if( !tp.summaryValue(val,type) ) return 0;
      return val > 0 ? val : 0;
    }

getPosOrBeamSpot()

Amg::Vector3D InDet::InDetDetailedTrackSelectorTool::getPosOrBeamSpot ( const xAOD::Vertex * vertex ) const

private

Definition at line 401 of file InDetDetailedTrackSelectorTool.cxx.

  {
    if(vertex) return vertex->position();
    if(m_useEventInfoBs){
      SG::ReadHandle<xAOD::EventInfo> evt(m_eventInfo_key);
      if (evt.isValid()) {
        InDet::BeamSpotData temp(evt->beamStatus(), evt->beamPosX(), evt->beamPosY(), evt->beamPosZ(),
                                 evt->beamPosSigmaX(), evt->beamPosSigmaY(), evt->beamPosSigmaZ(),
                                 evt->beamTiltXZ(), evt->beamTiltYZ(), evt->beamPosSigmaXY());
        return temp.beamVtx().position();
      } else {
        ATH_MSG_WARNING( " Cannot get beamSpot center from xAOD::EventInfo. Using (0,0,0)... " );
        return Amg::Vector3D(0,0,0);
      }
    }else{
      SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey };
      if (beamSpotHandle.isValid()) {
        return beamSpotHandle->beamVtx().position();
      } else {
        ATH_MSG_WARNING( " Cannot get beamSpot center from BeamSpotData. Using (0,0,0)... " );
        return Amg::Vector3D(0,0,0);
      }
    }
  }

initialize()

StatusCode InDet::InDetDetailedTrackSelectorTool::initialize

(

)

Definition at line 68 of file InDetDetailedTrackSelectorTool.cxx.

                                            {
    if(m_trackSumTool.empty()){
      ATH_MSG_DEBUG("No TrackSummaryTool set. OK if running on AOD.");
    }
    m_trackSumToolAvailable = !m_trackSumTool.empty();
    ATH_CHECK(m_trackSumTool.retrieve(DisableTool{!m_trackSumToolAvailable}));
 
    if(m_useSharedHitInfo && m_particleCreator.empty()){
      ATH_MSG_DEBUG("No TrackParticleCreatorTool set but shared hit selection used. OK if running on AOD.");
    }
    m_partCreatorToolAvailable = m_useSharedHitInfo && !m_particleCreator.empty();
    ATH_CHECK(m_particleCreator.retrieve(DisableTool{!m_partCreatorToolAvailable}));
 
    ATH_CHECK( m_extrapolator.retrieve() );
    ATH_CHECK(m_beamSpotKey.initialize(!m_useEventInfoBs));
    ATH_CHECK(m_eventInfo_key.initialize(m_useEventInfoBs));
    if (m_useEtaDepententMinHitTrt || m_useEtaDepententMinHitTrtWithOutliers){
      if(m_trtDCTool.empty()) {
    ATH_MSG_ERROR(" Eta dependent cut on number of TRT hits requested but TrtDCCutTool not specified. ");
    return StatusCode::FAILURE;
      } else if(m_trtDCTool.retrieve().isFailure()) {
    ATH_MSG_ERROR(" Unable to retrieve tool "<<m_trtDCTool);
    return StatusCode::FAILURE;
      }
      ATH_MSG_DEBUG("Retrieved tool "<<m_trtDCTool);
      if(m_useEtaDepententMinHitTrt){
    ATH_MSG_DEBUG("Using eta dependent cut on number of TRT hits.");
      }
      if(m_useEtaDepententMinHitTrtWithOutliers){
    ATH_MSG_DEBUG("Using eta dependent cut on number of TRT hits + outliers.");
      }
    }else{
      m_trtDCTool.disable();
    }
 
    // Read handle for AtlasFieldCacheCondObj
    ATH_CHECK( m_fieldCacheCondObjInputKey.initialize() );
 
    ATH_MSG_DEBUG("Using cuts on the number of Silicon hits");
    if(m_usePtDependentCuts) {
        //checking whether sizes of cuts and pt interval expressed in vectors match
        if( m_ptBenchmarks.size() != m_nSCTValues.size()){
          ATH_MSG_ERROR( "Number of cuts DOES NOT match the number of intervals to apply. Please check jobOptions. ");
          return StatusCode::FAILURE;
        } else if (m_ptBenchmarks.empty()){
          ATH_MSG_ERROR( "Zero vectors for number of cuts and pt intervals. Please check jobOptions. ");
          return StatusCode::FAILURE;
        }//end of vector size protection block
    }//end of memory protection
    return StatusCode::SUCCESS;
  }

inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

interfaceID()

const InterfaceID & Trk::ITrackSelectorTool::interfaceID ( )

inlinestaticinherited

Definition at line 38 of file ITrackSelectorTool.h.

    {
     return IID_ITrackSelectorTool;
    }

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

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.

preselectionBeforeExtrapolation()

bool InDet::InDetDetailedTrackSelectorTool::preselectionBeforeExtrapolation ( const Trk::Perigee & myPerigee ) const

private

Definition at line 1067 of file InDetDetailedTrackSelectorTool.cxx.

  {
    const AmgVector(5)& perigeeParms = myPerigee.parameters();
 
    // only check pt if mag. field is on
    const EventContext& ctx = Gaudi::Hive::currentContext();
    SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCacheCondObjInputKey, ctx};
    const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
    if (fieldCondObj == nullptr) {
      ATH_MSG_ERROR("execute: Failed to retrieve AtlasFieldCacheCondObj with key " << m_fieldCacheCondObjInputKey.key());
      return false;
    }
    MagField::AtlasFieldCache fieldCache;
    fieldCondObj->getInitializedCache (fieldCache);
 
    if (fieldCache.solenoidOn()){//B field
      if (perigeeParms[Trk::qOverP] == 0.) {
    ATH_MSG_DEBUG("Track rejected because of perigee qOverP == 0.");
    return false;
      }
      double p = std::fabs(1./perigeeParms[Trk::qOverP]);
      if (p<m_pMin) {
    ATH_MSG_DEBUG("Track rejected because of p " << p << " < " << m_pMin);
    return false;
      }
      double pt = p*std::sin(perigeeParms[Trk::theta]);
      if (pt<m_pTMin) {
    ATH_MSG_DEBUG("Track rejected because of pt " << pt << " < " << m_pTMin);
    return false;
      }
    }
 
    if (std::fabs(perigeeParms[Trk::d0]) > m_d0MaxPreselection) {
      ATH_MSG_DEBUG("Track rejected because of fabs(d0) "<<std::fabs(perigeeParms[Trk::d0])<<" < "<<m_d0MaxPreselection);
      return false;
    }
 
    return true;
  }

renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

Reimplemented in 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.

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_addToMinHitTrt

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_addToMinHitTrt

private

Initial value:

{this, "addToMinHitTrt", 0,
       "add to/subtract from eta dependent minimum nimber of TRT hits"}

Definition at line 177 of file InDetDetailedTrackSelectorTool.h.

      {this, "addToMinHitTrt", 0,
       "add to/subtract from eta dependent minimum nimber of TRT hits"};

m_addToMinHitTrtWithOutliers

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_addToMinHitTrtWithOutliers

private

Initial value:

{this, "addToMinHitTrtWithOutliers", 0,
       "add to/subtract from eta dependent minimum nimber of TRT hits + outliers"}

Definition at line 183 of file InDetDetailedTrackSelectorTool.h.

      {this, "addToMinHitTrtWithOutliers", 0,
       "add to/subtract from eta dependent minimum nimber of TRT hits + outliers"};

m_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> InDet::InDetDetailedTrackSelectorTool::m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"}

private

Definition at line 201 of file InDetDetailedTrackSelectorTool.h.

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

m_d0MaxPreselection

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_d0MaxPreselection {this, "d0MaxPreselection", 10.}

private

Definition at line 188 of file InDetDetailedTrackSelectorTool.h.

{this, "d0MaxPreselection", 10.};

m_d0significanceMax

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_d0significanceMax {this, "d0significanceMax", -1., "max IP significance d0 (-1 switches it off)"}

private

Definition at line 109 of file InDetDetailedTrackSelectorTool.h.

{this, "d0significanceMax", -1., "max IP significance d0 (-1 switches it off)"};

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_etaMax

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_etaMax {this, "etaMax", 9999., "max. pseudo-rapidity"}

private

Definition at line 113 of file InDetDetailedTrackSelectorTool.h.

{this, "etaMax", 9999., "max. pseudo-rapidity"};

m_eventInfo_key

SG::ReadHandleKey<xAOD::EventInfo> InDet::InDetDetailedTrackSelectorTool::m_eventInfo_key {this, "EventInfo", "EventInfo", "Input event information"}

private

Definition at line 225 of file InDetDetailedTrackSelectorTool.h.

{this, "EventInfo", "EventInfo", "Input event information"};

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_extrapolator

ToolHandle<Trk::IExtrapolator> InDet::InDetDetailedTrackSelectorTool::m_extrapolator {this, "Extrapolator", "Trk::Extrapolator"}

private

Definition at line 199 of file InDetDetailedTrackSelectorTool.h.

{this, "Extrapolator", "Trk::Extrapolator"};

m_fieldCacheCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> InDet::InDetDetailedTrackSelectorTool::m_fieldCacheCondObjInputKey

private

Initial value:

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

Definition at line 211 of file InDetDetailedTrackSelectorTool.h.

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

m_fitChi2

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_fitChi2 {this, "fitChi2", 99999., "max. fit chi2"}

private

Definition at line 170 of file InDetDetailedTrackSelectorTool.h.

{this, "fitChi2", 99999., "max. fit chi2"};

m_fitChi2OnNdfMax

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_fitChi2OnNdfMax {this, "fitChi2OnNdfMax", 999., "max. fitchi2/ndf"}

private

Definition at line 172 of file InDetDetailedTrackSelectorTool.h.

{this, "fitChi2OnNdfMax", 999., "max. fitchi2/ndf"};

m_fitProb

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_fitProb {this, "fitProb", -1., "min. fit chi2 probability"}

private

Definition at line 171 of file InDetDetailedTrackSelectorTool.h.

{this, "fitProb", -1., "min. fit chi2 probability"};

m_inDetTestPixelLayerTool

ToolHandle< InDet::IInDetTestPixelLayerTool > InDet::InDetDetailedTrackSelectorTool::m_inDetTestPixelLayerTool

private

Initial value:

{this, "InDetTestPixelLayerTool", "",
       "Tool to test if the track crosses a dead module on the b-layer"}

Definition at line 207 of file InDetDetailedTrackSelectorTool.h.

      {this, "InDetTestPixelLayerTool", "",
       "Tool to test if the track crosses a dead module on the b-layer"};

m_IPd0Max

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_IPd0Max {this, "IPd0Max", 2.*CLHEP::mm, "max. d0: |d0|<d0Max"}

private

Definition at line 103 of file InDetDetailedTrackSelectorTool.h.

{this, "IPd0Max", 2.*CLHEP::mm, "max. d0: |d0|<d0Max"};

m_IPz0Max

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_IPz0Max {this, "IPz0Max", 1.5*CLHEP::mm, "max. z0: |z0*sin(theta)|<z0Max"}

private

Definition at line 104 of file InDetDetailedTrackSelectorTool.h.

{this, "IPz0Max", 1.5*CLHEP::mm, "max. z0: |z0*sin(theta)|<z0Max"};

m_nDoubleHoles

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nDoubleHoles {this, "nDoubleHoles", 999, "max number of double-holes in SCT"}

private

Definition at line 161 of file InDetDetailedTrackSelectorTool.h.

{this, "nDoubleHoles", 999, "max number of double-holes in SCT"};

m_nHitBLayer

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitBLayer {this, "nHitBLayer", 1, "at least n hits in Blayer"}

private

Definition at line 117 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitBLayer", 1, "at least n hits in Blayer"};

m_nHitPix

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitPix {this, "nHitPix", 2, "at least n hits in pixels"}

private

Definition at line 119 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitPix", 2, "at least n hits in pixels"};

m_nHitPixPhysical

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitPixPhysical {this, "nHitPixPhysical", 0, "at least n physical hits in pixel"}

private

Definition at line 122 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitPixPhysical", 0, "at least n physical hits in pixel"};

m_nHitSct

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitSct {this, "nHitSct", 0, "at least n hits in SCT"}

private

Definition at line 120 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitSct", 0, "at least n hits in SCT"};

m_nHitSi

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitSi {this, "nHitSi", 7, "at least n hits in pixels+SCT"}

private

Definition at line 121 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitSi", 7, "at least n hits in pixels+SCT"};

m_nHitSiPhysical

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitSiPhysical {this, "nHitSiPhysical", 3, "at least n physical hits in pixel+SCT"}

private

Definition at line 124 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitSiPhysical", 3, "at least n physical hits in pixel+SCT"};

m_nHitTrt

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitTrt {this, "nHitTrt", 0, "at least n hits in TRT"}

private

Definition at line 126 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitTrt", 0, "at least n hits in TRT"};

m_nHitTrtHighE

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitTrtHighE {this, "nHitTrtHighE", 0, "at least n high threshold hits in TRT"}

private

Definition at line 132 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitTrtHighE", 0, "at least n high threshold hits in TRT"};

m_nHitTrtHighEFraction

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_nHitTrtHighEFraction {this, "nHitTrtHighEFractionMax", 999., "maximum x fraction of transition hits in TRT"}

private

Definition at line 138 of file InDetDetailedTrackSelectorTool.h.

{this, "nHitTrtHighEFractionMax", 999., "maximum x fraction of transition hits in TRT"};

m_nHitTrtHighEFractionWithOutliers

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_nHitTrtHighEFractionWithOutliers

private

Initial value:

{this, "nHitTrtHighEFractionWithOutliersMax", 999.,
       "maximum x fraction of transition hits in TRT (including outliers)"}

Definition at line 140 of file InDetDetailedTrackSelectorTool.h.

      {this, "nHitTrtHighEFractionWithOutliersMax", 999.,
       "maximum x fraction of transition hits in TRT (including outliers)"};

m_nHitTrtPlusOutliers

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitTrtPlusOutliers

private

Initial value:

{this, "nHitTrtPlusOutliers", 0,
       "at least n hits in TRT (including outliers)"}

Definition at line 127 of file InDetDetailedTrackSelectorTool.h.

      {this, "nHitTrtPlusOutliers", 0,
       "at least n hits in TRT (including outliers)"};

m_nHitTrtPlusOutliersHighE

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHitTrtPlusOutliersHighE

private

Initial value:

{this, "nHitTrtPlusOutliersHighE", 0,
       "at least n high threshold hits in TRT (including outliers)"}

Definition at line 134 of file InDetDetailedTrackSelectorTool.h.

      {this, "nHitTrtPlusOutliersHighE", 0,
       "at least n high threshold hits in TRT (including outliers)"};

m_nHoles

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHoles {this, "nHoles", 999, "max. number of holes in pixel+SCT"}

private

Definition at line 159 of file InDetDetailedTrackSelectorTool.h.

{this, "nHoles", 999, "max. number of holes in pixel+SCT"};

m_nHolesPix

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHolesPix {this, "nHolesPixel", 999, "max. number of holes in pixels"}

private

Definition at line 163 of file InDetDetailedTrackSelectorTool.h.

{this, "nHolesPixel", 999, "max. number of holes in pixels"};

m_nHolesSct

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nHolesSct {this, "nHolesSct", 999, "max. number of holes in SCT"}

private

Definition at line 165 of file InDetDetailedTrackSelectorTool.h.

{this, "nHolesSct", 999, "max. number of holes in SCT"};

m_nSCTValues

IntegerArrayProperty InDet::InDetDetailedTrackSelectorTool::m_nSCTValues {this, "SCTCutValues", {}}

private

Definition at line 223 of file InDetDetailedTrackSelectorTool.h.

{this, "SCTCutValues", {}};

m_nSharedBLayer

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nSharedBLayer {this, "nSharedBLayer", 0, "max. number of shared hits in B layer"}

private

Definition at line 150 of file InDetDetailedTrackSelectorTool.h.

{this, "nSharedBLayer", 0, "max. number of shared hits in B layer"};

m_nSharedPix

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nSharedPix {this, "nSharedPix", 0, "max. number of shared hits in pixels"}

private

Definition at line 152 of file InDetDetailedTrackSelectorTool.h.

{this, "nSharedPix", 0, "max. number of shared hits in pixels"};

m_nSharedSct

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nSharedSct {this, "nSharedSct", 1, "max. number of shared hits in SCT"}

private

Definition at line 154 of file InDetDetailedTrackSelectorTool.h.

{this, "nSharedSct", 1, "max. number of shared hits in SCT"};

m_nSharedSi

IntegerProperty InDet::InDetDetailedTrackSelectorTool::m_nSharedSi {this, "nSharedSi", 999, "max. number of shared hits in pixels+SCT"}

private

Definition at line 156 of file InDetDetailedTrackSelectorTool.h.

{this, "nSharedSi", 999, "max. number of shared hits in pixels+SCT"};

m_partCreatorToolAvailable

bool InDet::InDetDetailedTrackSelectorTool::m_partCreatorToolAvailable = false

private

Definition at line 216 of file InDetDetailedTrackSelectorTool.h.

m_particleCreator

ToolHandle<Trk::ITrackParticleCreatorTool> InDet::InDetDetailedTrackSelectorTool::m_particleCreator {this, "TrackParticleCreatorTool", ""}

private

Definition at line 197 of file InDetDetailedTrackSelectorTool.h.

{this, "TrackParticleCreatorTool", ""};

m_pMin

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_pMin {this, "pMin", 0., "min. p = pT/cos(theta): |p| > pMin"}

private

Definition at line 102 of file InDetDetailedTrackSelectorTool.h.

{this, "pMin", 0., "min. p = pT/cos(theta): |p| > pMin"};

m_ptBenchmarks

FloatArrayProperty InDet::InDetDetailedTrackSelectorTool::m_ptBenchmarks {this, "PtBenchmarks", {}}

private

Definition at line 221 of file InDetDetailedTrackSelectorTool.h.

{this, "PtBenchmarks", {}};

m_pTMin

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_pTMin {this, "pTMin", 1.*CLHEP::GeV, "min. pT: |pT|>pTMin"}

private

Definition at line 101 of file InDetDetailedTrackSelectorTool.h.

{this, "pTMin", 1.*CLHEP::GeV, "min. pT: |pT|>pTMin"};

m_scaleMinHitTrt

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_scaleMinHitTrt

private

Initial value:

{this, "scaleMinHitTrt", 1.,
       "scale the eta dependent minimum number of TRT hits; scaling is only applied if m_addToMinHitTrt==0"}

Definition at line 174 of file InDetDetailedTrackSelectorTool.h.

      {this, "scaleMinHitTrt", 1.,
       "scale the eta dependent minimum number of TRT hits; scaling is only applied if m_addToMinHitTrt==0"};

m_scaleMinHitTrtWithOutliers

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_scaleMinHitTrtWithOutliers

private

Initial value:

{this, "scaleMinHitTrtWithOutliers", 1.,
       "scale the eta dependent minimum number of TRT hits + outliers; scaling is only applied if m_addToMinHitTrtWithOutliers==0"}

Definition at line 180 of file InDetDetailedTrackSelectorTool.h.

      {this, "scaleMinHitTrtWithOutliers", 1.,
       "scale the eta dependent minimum number of TRT hits + outliers; scaling is only applied if m_addToMinHitTrtWithOutliers==0"};

m_sigIPd0Max

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_sigIPd0Max {this, "sigIPd0Max", 999.*CLHEP::mm, "max d0 error"}

private

Definition at line 106 of file InDetDetailedTrackSelectorTool.h.

{this, "sigIPd0Max", 999.*CLHEP::mm, "max d0 error"};

m_sigIPz0Max

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_sigIPz0Max {this, "sigIPz0Max", 999.*CLHEP::mm, "max (error only due to z0)*sin(theta)"}

private

Definition at line 107 of file InDetDetailedTrackSelectorTool.h.

{this, "sigIPz0Max", 999.*CLHEP::mm, "max (error only due to z0)*sin(theta)"};

m_trackSumTool

ToolHandle<Trk::ITrackSummaryTool> InDet::InDetDetailedTrackSelectorTool::m_trackSumTool {this, "TrackSummaryTool", "Trk::TrackSummaryTool"}

private

Definition at line 195 of file InDetDetailedTrackSelectorTool.h.

{this, "TrackSummaryTool", "Trk::TrackSummaryTool"};

m_trackSumToolAvailable

bool InDet::InDetDetailedTrackSelectorTool::m_trackSumToolAvailable = true

private

Definition at line 215 of file InDetDetailedTrackSelectorTool.h.

m_trtDCTool

ToolHandle<ITrtDriftCircleCutTool> InDet::InDetDetailedTrackSelectorTool::m_trtDCTool

private

Initial value:

{this, "TrtDCCutTool", "InDet::InDetTrtDriftCircleCutTool",
       "Tool to get eta dependent cut on number of TRT hits"}

Definition at line 203 of file InDetDetailedTrackSelectorTool.h.

      {this, "TrtDCCutTool", "InDet::InDetTrtDriftCircleCutTool",
       "Tool to get eta dependent cut on number of TRT hits"};

m_TrtMaxEtaAcceptance

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_TrtMaxEtaAcceptance

private

Initial value:

{this, "TrtMaxEtaAcceptance", 999.,
       "limit of eta regions where trt hits are expected"}

Definition at line 144 of file InDetDetailedTrackSelectorTool.h.

      {this, "TrtMaxEtaAcceptance", 999.,
       "limit of eta regions where trt hits are expected"};

m_useEtaDepententMinHitTrt

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_useEtaDepententMinHitTrt {this, "useEtaDepententMinHitTrt", false}

private

Definition at line 190 of file InDetDetailedTrackSelectorTool.h.

{this, "useEtaDepententMinHitTrt", false};

m_useEtaDepententMinHitTrtWithOutliers

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_useEtaDepententMinHitTrtWithOutliers {this, "useEtaDepententMinHitTrtWithOutliers", false}

private

Definition at line 192 of file InDetDetailedTrackSelectorTool.h.

{this, "useEtaDepententMinHitTrtWithOutliers", false};

m_useEventInfoBs

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_useEventInfoBs {this, "UseEventInfoBS", false, "Access beamspot via the EvenInfo object"}

private

Definition at line 227 of file InDetDetailedTrackSelectorTool.h.

{this, "UseEventInfoBS", false, "Access beamspot via the EvenInfo object"};

m_usePreselectionCuts

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_usePreselectionCuts {this, "usePreselectionCuts", false}

private

Definition at line 187 of file InDetDetailedTrackSelectorTool.h.

{this, "usePreselectionCuts", false};

m_usePtDependentCuts

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_usePtDependentCuts {this, "UsePtDependentCuts", false}

private

Definition at line 219 of file InDetDetailedTrackSelectorTool.h.

{this, "UsePtDependentCuts", false};

m_useSharedHitInfo

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_useSharedHitInfo {this, "useSharedHitInfo", false}

private

Definition at line 149 of file InDetDetailedTrackSelectorTool.h.

{this, "useSharedHitInfo", false};

m_useTrackQualityInfo

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_useTrackQualityInfo {this, "useTrackQualityInfo", true}

private

Definition at line 169 of file InDetDetailedTrackSelectorTool.h.

{this, "useTrackQualityInfo", true};

m_useTrackSummaryInfo

BooleanProperty InDet::InDetDetailedTrackSelectorTool::m_useTrackSummaryInfo {this, "useTrackSummaryInfo", true}

private

Definition at line 116 of file InDetDetailedTrackSelectorTool.h.

{this, "useTrackSummaryInfo", 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.

m_z0Max

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_z0Max {this, "z0Max", 9999.*CLHEP::mm, "max. z0: |z0|<z0Max"}

private

Definition at line 105 of file InDetDetailedTrackSelectorTool.h.

{this, "z0Max", 9999.*CLHEP::mm, "max. z0: |z0|<z0Max"};

m_z0significanceMax

DoubleProperty InDet::InDetDetailedTrackSelectorTool::m_z0significanceMax {this, "z0significanceMax", -1., "max IP significance z0 (-1 switches it off)"}

private

Definition at line 111 of file InDetDetailedTrackSelectorTool.h.

{this, "z0significanceMax", -1., "max IP significance z0 (-1 switches it off)"};

---

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

• InDetDetailedTrackSelectorTool.h
• InDetDetailedTrackSelectorTool.cxx