InDet::InDetConversionTrackSelectorTool Class Reference

#include <InDetConversionTrackSelectorTool.h>

Inheritance diagram for InDet::InDetConversionTrackSelectorTool:

Collaboration diagram for InDet::InDetConversionTrackSelectorTool:

Public Member Functions
virtual StatusCode	initialize () override
	InDetConversionTrackSelectorTool (const std::string &t, const std::string &n, const IInterface *p)
virtual	~InDetConversionTrackSelectorTool ()=default
virtual bool	decision (const Trk::Track &track, const Trk::Vertex *vertex) const override final
	Select a Trk::Track.
virtual bool	decision (const Trk::TrackParticleBase &track, const Trk::Vertex *vertex) const override final
	Select a Trk::TrackParticleBase.
virtual bool	decision (const xAOD::TrackParticle &, const xAOD::Vertex *) const override final
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
unsigned int	getEtaBin (const Trk::Perigee &perigee) const
Amg::Vector3D	getPosOrBeamSpot (const EventContext &ctx, const xAOD::Vertex *) const
Trk::Vertex *	getBeamSpot (const EventContext &ctx) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ToolHandle< Trk::IExtrapolator >	m_extrapolator
SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfo_key {this, "EventInfo", "EventInfo", "Input event information"}
Gaudi::Property< bool >	m_useEventInfoBs {this,"UseEventInfoBS",false}
DoubleProperty	m_maxSiD0 {this, "maxSiD0", 35., "Maximal d0 at (0,0,0) for tracks with Si hits"}
	Properties for track selection:all cuts are ANDed.
DoubleProperty	m_maxTrtD0 {this, "maxTrtD0", 100., "Maximal d0 at (0,0,0) for standalone TRT tracks"}
DoubleProperty	m_maxSiZ0 {this, "maxSiZ0", 200., "Maximal z0 at (0,0,0)"}
DoubleProperty	m_maxTrtZ0 {this, "maxTrtZ0", 1200., "Maximal z0 at (0,0,0) for standalone TRT tracks"}
DoubleProperty	m_minPt {this, "minPt", 500., "Minimum Pt of tracks"}
DoubleProperty	m_trRatio1 {this, "RatioCut1", 0.5, "TR ratio for tracks with 15 or less TRT hits"}
DoubleProperty	m_trRatio2 {this, "RatioCut2", 0.1, "TR ratio for tracks with 16 to 25 TRT hits"}
DoubleProperty	m_trRatio3 {this, "RatioCut3", 0.05, "TR ratio for tracks with 26 or more TRT hits"}
DoubleProperty	m_trRatioTRT {this, "RatioTRT", 0.1, "TR ratio for all TRT only tracks"}
DoubleArrayProperty	m_TRTTrksEtaBins
DoubleArrayProperty	m_TRTTrksBinnedRatioTRT
DoubleProperty	m_trRatioV0 {this, "RatioV0", 1., "TR ratio for pion selection during V0 reconstruction"}
DoubleProperty	m_sD0_Si {this, "significanceD0_Si", 2., "Cut on D0 significance of Si tracks"}
DoubleProperty	m_sD0_Trt {this, "significanceD0_Trt", 0.5, "Cut on D0 significance of TRT tracks"}
DoubleProperty	m_sZ0_Trt {this, "significanceZ0_Trt", 3., "Cut on Z0 significance of TRT tracks"}
BooleanProperty	m_isConv {this, "IsConversion", true, "Conversion flag"}
BooleanProperty	m_PIDonlyForXe {this, "PIDonlyForXe", false, "Only check TRT PID if all hits are Xe hits"}
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 32 of file InDetConversionTrackSelectorTool.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

InDetConversionTrackSelectorTool()

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

Definition at line 21 of file InDetConversionTrackSelectorTool.cxx.

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

~InDetConversionTrackSelectorTool()

virtual InDet::InDetConversionTrackSelectorTool::~InDetConversionTrackSelectorTool ( )

virtualdefault

Member Function Documentation

decision() [1/3]

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

finaloverridevirtual

Select a Trk::Track.

Implements Trk::ITrackSelectorTool.

Definition at line 78 of file InDetConversionTrackSelectorTool.cxx.

 {
   const EventContext& ctx = Gaudi::Hive::currentContext();
   bool pass = false;
   const Trk::Perigee* perigee=dynamic_cast<const Trk::Perigee*>(track.perigeeParameters());
   const bool vertexSuppliedByUser{vx!=nullptr};
   const Trk::Vertex* myVertex=vx;
   //in case no Vertex is provided by the user, beam position will be used if available
   if (not vertexSuppliedByUser) {
     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 (not vertexSuppliedByUser) delete myVertex;
       return false;
     }
   }
 
   const Trk::TrackParameters* extrapolatedParameters =
     m_extrapolator->extrapolate(ctx,
                                 *firstmeaspar,
                                 perigeeSurface,
                                 Trk::anyDirection,
                                 true,
                                 track.info().particleHypothesis()).release();
   perigee = extrapolatedParameters
               ? dynamic_cast<const Trk::Perigee*>(extrapolatedParameters)
               : nullptr;
   if (perigee==nullptr || !perigee->covariance() ) {
     ATH_MSG_WARNING( "Track Selector failed to extrapolate track to the vertex: " << myVertex->position() );
     if (extrapolatedParameters!=nullptr) {
       ATH_MSG_WARNING( "The return object of the extrapolator was not a perigee even if a perigeeSurface was used!");
       delete extrapolatedParameters;
       if (not vertexSuppliedByUser) delete myVertex;
       return false;
     }
     if (not vertexSuppliedByUser) delete myVertex;
     return false;
   }
 
   double qOverP = perigee->parameters()[Trk::qOverP];
   double pt = std::fabs(1./qOverP)*std::sin(perigee->parameters()[Trk::theta]);
   double d0 = perigee->parameters()[Trk::d0];
   double z0 = perigee->parameters()[Trk::z0];
   const Trk::TrackSummary* tSum = track.trackSummary();
   if(tSum){
     double ratioTrk = 1.0;
     int nclus  = tSum->get(Trk::numberOfPixelHits) + tSum->get(Trk::numberOfSCTHits);
     bool isSilicon = (nclus > 0);
     int nTrtHits       = tSum->get(Trk::numberOfTRTHits);
     int nTrtOutliers   = tSum->get(Trk::numberOfTRTOutliers);
     int ntrt           = nTrtHits + nTrtOutliers;
     int nTrtXenonHits  = tSum->get(Trk::numberOfTRTXenonHits);
     if(m_isConv) {
       if(ntrt > 0 && (!m_PIDonlyForXe || nTrtXenonHits==ntrt) ){
         // only check TRT PID if m_PIDonlyForXe is false or all TRT hits are Xenon hits
          ATH_MSG_FATAL( "eProbabilityHT not available for Trk::Tracks only xAOD::TrackParticle objects" );
       }
 
       // Start of track cuts
       if ( pt >= m_minPt ) {
 
         // Silicon track cuts
         if ( isSilicon && (std::fabs(d0)<=m_maxSiD0) && (fabs(z0)<=m_maxSiZ0) ) {
           if((ntrt<=15 && ratioTrk>=m_trRatio1) ||
              (ntrt>15 && ntrt<=25 && ratioTrk>=m_trRatio2) ||
              (ntrt>25 && ratioTrk>=m_trRatio3)) pass = true;
         }
 
         // TRT-only track cuts
         if ( (not isSilicon) && (std::fabs(d0)<=m_maxTrtD0) && (std::fabs(z0)<=m_maxTrtZ0) ) {
 
           unsigned int eta_bin = getEtaBin(*perigee);
           // TRT-only Tracks: eProbabilityHT cut below
           // See InDet Trt Track Scoring Tool for nTRT cut on TRT-only tracks
           double trRatioTRT = std::max( m_trRatioTRT.value(), m_TRTTrksBinnedRatioTRT[eta_bin] );
 
           if ( ratioTrk >= trRatioTRT ) pass = true; // TRT Track cuts
         }
       } // end of track cuts for isConv
 
     } else {
       //The cuts below are necessary for the V0 track selection
       const AmgSymMatrix(5)& err = *perigee->covariance();
       double sd0sq = err(0,0);
       double sd0 = (sd0sq>0.)?std::sqrt(sd0sq):0.;
       double sz0sq = err(1,1);
       double sz0 = (sz0sq>0.)?std::sqrt(sz0sq):0.;
       if(nclus == 0){
           if(std::fabs(d0)>=m_sD0_Trt*sd0 && std::fabs(d0)<=m_maxTrtD0 && std::fabs(z0)<=m_sZ0_Trt*sz0 && pt>=m_minPt) pass = true;
       }else{
           if(std::fabs(d0)>=m_sD0_Si*sd0 && std::fabs(z0)<=m_maxSiZ0 && pt>=m_minPt) pass = true;
       }
       ratioTrk = 1.0;
       if(ntrt > 0 && (!m_PIDonlyForXe || nTrtXenonHits==ntrt) ) { // only check TRT PID if m_PIDonlyForXe is false or all TRT hits are Xenon hits
          ATH_MSG_FATAL( "eProbabilityHT not available for Trk::Tracks only xAOD::TrackParticle objects" );
       }
       if(ratioTrk>m_trRatioV0) pass = false;
     }
 
   } else pass = false;
   if (not vertexSuppliedByUser) delete myVertex;
   if (perigee!=track.perigeeParameters()) delete perigee;
   return pass;
 }

decision() [2/3]

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

finaloverridevirtual

Select a Trk::TrackParticleBase.

Implements Trk::ITrackSelectorTool.

Definition at line 197 of file InDetConversionTrackSelectorTool.cxx.

 {
   const EventContext& ctx = Gaudi::Hive::currentContext();
   bool pass = false;
   const Trk::TrackParameters* definintParameters=&(track.definingParameters());
   const Trk::Perigee* perigee=dynamic_cast<const Trk::Perigee*>(definintParameters);
   const Trk::Vertex* myVertex=vx;
   const bool vertexSuppliedByUser{vx!=nullptr};
   //in case no Vertex is provided by the user, beam position will be used if available
   if (not vertexSuppliedByUser) {
     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) {
     //using perigee instead of firstmeasurement, since first measurement was not found...
     firstmeaspar=&(track.definingParameters());
     if (!firstmeaspar){
       ATH_MSG_WARNING( " Track Paraemters at first measurement not found. Perigee not found. Cannot do TrackSelection..." );
       if (not vertexSuppliedByUser) delete myVertex;
       return false;
     }
   }
   const Trk::TrackParameters* extrapolatedParameters =
     m_extrapolator->extrapolate(
       ctx, *firstmeaspar, perigeeSurface, Trk::anyDirection, true, Trk::pion).release();
   perigee = extrapolatedParameters
               ? dynamic_cast<const Trk::Perigee*>(extrapolatedParameters)
               : nullptr;
   if (perigee == nullptr || !perigee->covariance()) {
     ATH_MSG_WARNING( "Track Selector failed to extrapolate track to the vertex: " << myVertex->position() );
     if (extrapolatedParameters!=nullptr) {
       ATH_MSG_WARNING( "The return object of the extrapolator was not a perigee even if a perigeeSurface was used!" );
       delete extrapolatedParameters;
       if (not vertexSuppliedByUser) delete myVertex;
       return false;
     }
     if (not vertexSuppliedByUser) delete myVertex;
     return false;
   }
 
   double qOverP = perigee->parameters()[Trk::qOverP];
   double pt = std::fabs(1./qOverP)*std::sin(perigee->parameters()[Trk::theta]);
   double d0 = perigee->parameters()[Trk::d0];
   double z0 = perigee->parameters()[Trk::z0];
   const Trk::TrackSummary* tSum = track.trackSummary();
   if(tSum){
     double ratioTrk = 1.0;
     int nclus  = tSum->get(Trk::numberOfPixelHits) + tSum->get(Trk::numberOfSCTHits);
     bool isSilicon = (nclus > 0);
     int nTrtHits       = tSum->get(Trk::numberOfTRTHits);
     int nTrtOutliers   = tSum->get(Trk::numberOfTRTOutliers);
     int ntrt           = nTrtHits + nTrtOutliers;
     int nTrtXenonHits  = tSum->get(Trk::numberOfTRTXenonHits);
 
     if(m_isConv){
       if(ntrt > 0 && (!m_PIDonlyForXe || nTrtXenonHits==ntrt) ){
         // only check TRT PID if m_PIDonlyForXe is false or all TRT hits are Xenon hits
         ATH_MSG_FATAL( "eProbabilityHT not available for Trk::TrackParticleBase only xAOD::TrackParticle objects" );
       }
 
       // Start of track cuts
       if ( pt >= m_minPt ) {
 
         // Silicon track cuts
         if ( isSilicon && (std::fabs(d0)<=m_maxSiD0) && (fabs(z0)<=m_maxSiZ0) ) {
           if((ntrt<=15 && ratioTrk>=m_trRatio1) ||
              (ntrt>15 && ntrt<=25 && ratioTrk>=m_trRatio2) ||
              (ntrt>25 && ratioTrk>=m_trRatio3)) pass = true;
         }
 
         // TRT-only track cuts
         if ( (not isSilicon) && (std::fabs(d0)<=m_maxTrtD0) && (std::fabs(z0)<=m_maxTrtZ0) ) {
 
           unsigned int eta_bin = getEtaBin(*perigee);
           // TRT-only Tracks: eProbabilityHT cuts below
           // See InDet Trt Track Scoring Tool for nTRT cut on TRT-only tracks
           double trRatioTRT = std::max( m_trRatioTRT.value(), m_TRTTrksBinnedRatioTRT[eta_bin] );
 
           if ( ratioTrk >= trRatioTRT ) pass = true; // TRT Track cuts
         }
       } // end of track cuts for isConv
 
     } else {
       //The cuts below are necessary for the V0 track selection
       const AmgSymMatrix(5)& err = *perigee->covariance();
       double sd0sq = err(0,0);
       double sd0 = (sd0sq>0.)?std::sqrt(sd0sq):0.;
       double sz0sq = err(1,1);
       double sz0 = (sz0sq>0.)?std::sqrt(sz0sq):0.;
       if(nclus == 0){
           if(std::fabs(d0)>=m_sD0_Trt*sd0 && std::fabs(d0)<= m_maxTrtD0 && std::fabs(z0)<=m_sZ0_Trt*sz0 && pt>=m_minPt) pass = true;
       }else{
           if(std::fabs(d0)>=m_sD0_Si*sd0 && std::fabs(z0)<=m_maxSiZ0 && pt>=m_minPt) pass = true;
       }
 
       ratioTrk = 1.0;
       if(ntrt > 0 && (!m_PIDonlyForXe || nTrtXenonHits==ntrt) ) {// only check TRT PID if m_PIDonlyForXe is false or all TRT hits are Xenon hits
          ATH_MSG_FATAL( "eProbabilityHT not available for Trk::TrackParticleBase only xAOD::TrackParticle objects" );
       }
       if(ratioTrk>m_trRatioV0) pass = false;
     }
   } else pass = false;
   if (not vertexSuppliedByUser) delete myVertex;
   if (perigee!=&(track.definingParameters())) delete perigee;
 
   return pass;
 }

decision() [3/3]

bool InDet::InDetConversionTrackSelectorTool::decision ( const xAOD::TrackParticle & tp, const xAOD::Vertex * vertex ) const

finaloverridevirtual

Implements Trk::ITrackSelectorTool.

Definition at line 338 of file InDetConversionTrackSelectorTool.cxx.

  {
    const EventContext& ctx = Gaudi::Hive::currentContext();
    bool pass = false;
    const Trk::Perigee& perigee=tp.perigeeParameters();
    // in case no Vertex is provided by the user, beam position will be used if
    // available
    Trk::PerigeeSurface perigeeSurface(getPosOrBeamSpot(ctx, vertex));
    const Trk::TrackParameters* extrapolatedParameters =
      m_extrapolator->extrapolate(
        ctx, perigee, perigeeSurface, Trk::anyDirection, false, Trk::pion).release();
    if (extrapolatedParameters == nullptr) {
      ATH_MSG_WARNING("Extrapolation to the vertex failed: " << perigeeSurface
                                                             << "\n"
                                                             << perigee);
      return false;
    }
    double qOverP = perigee.parameters()[Trk::qOverP];
    double pt = std::fabs(1./qOverP)*std::sin(perigee.parameters()[Trk::theta]);
    double d0 = extrapolatedParameters->parameters()[Trk::d0];
    double z0 = extrapolatedParameters->parameters()[Trk::z0];
 
    double ratioTrk = 1.0;
    int nclus  = getCount(tp,xAOD::numberOfPixelHits) + getCount(tp,xAOD::numberOfSCTHits);
    bool isSilicon = (nclus > 0);
    int nTrtHits       = getCount(tp,xAOD::numberOfTRTHits);
    int nTrtOutliers   = getCount(tp,xAOD::numberOfTRTOutliers);
    int ntrt           = nTrtHits + nTrtOutliers;
    int nTrtXenonHits  = getCount(tp,xAOD::numberOfTRTXenonHits);
 
    if(m_isConv){
     float temp(0);
     if(ntrt > 0 && (!m_PIDonlyForXe || nTrtXenonHits==ntrt) ){
       // only check TRT PID if m_PIDonlyForXe is false or all TRT hits are Xenon hits
       ratioTrk = tp.summaryValue(temp,xAOD::eProbabilityHT)  ?  temp: 0 ;
     }
 
     // Start of track cuts
     if ( pt >= m_minPt ) {
 
       // Silicon track cuts
       if ( isSilicon && (std::fabs(d0)<=m_maxSiD0) && (fabs(z0)<=m_maxSiZ0) ) {
         if((ntrt<=15 && ratioTrk>=m_trRatio1) ||
            (ntrt>15 && ntrt<=25 && ratioTrk>=m_trRatio2) ||
            (ntrt>25 && ratioTrk>=m_trRatio3)) pass = true;
       }
 
       // TRT-only track cuts
       if ( (not isSilicon) && (std::fabs(d0)<=m_maxTrtD0) && (std::fabs(z0)<=m_maxTrtZ0) ) {
 
         unsigned int eta_bin = getEtaBin(perigee);
         // TRT-only Tracks: eProbabilityHT cuts below
         // See InDet Trt Track Scoring Tool for nTRT cut on TRT-only tracks
         double trRatioTRT = std::max( m_trRatioTRT.value(), m_TRTTrksBinnedRatioTRT[eta_bin] );
 
         if ( ratioTrk >= trRatioTRT ) pass = true; // TRT Track cuts
       }
     } // end of track cuts for isConv
 
    } else {
     //The cuts below are necessary for the V0 track selection
     const AmgSymMatrix(5)& err = *perigee.covariance();
     double sd0sq = err(0,0);
     double sd0 = (sd0sq>0.)?std::sqrt(sd0sq):0.;
     double sz0sq = err(1,1);
     double sz0 = (sz0sq>0.)?std::sqrt(sz0sq):0.;
     if(nclus == 0){
       if(std::fabs(d0)>=m_sD0_Trt*sd0 && std::fabs(d0)<= m_maxTrtD0 && std::fabs(z0)<=m_sZ0_Trt*sz0 && pt>=m_minPt) pass = true;
     }else{
       if(std::fabs(d0)>=m_sD0_Si*sd0 && std::fabs(z0)<=m_maxSiZ0 && pt>=m_minPt) pass = true;
     }
     ratioTrk = 1.0;
     float temp(0);
     if(ntrt > 0 && (!m_PIDonlyForXe || nTrtXenonHits==ntrt) ) // only check TRT PID if m_PIDonlyForXe is false or all TRT hits are Xenon hits
      ratioTrk = tp.summaryValue(temp,xAOD::eProbabilityHT)  ?  temp: 0 ;
     if(ratioTrk>m_trRatioV0) pass = false;
    }
 
 
    delete extrapolatedParameters;
    return pass;
  }

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

getBeamSpot()

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

private

Definition at line 37 of file InDetConversionTrackSelectorTool.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::InDetConversionTrackSelectorTool::getCount ( const xAOD::TrackParticle & tp, xAOD::SummaryType type ) const

inlineprivate

Definition at line 58 of file InDetConversionTrackSelectorTool.h.

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

getEtaBin()

unsigned int InDet::InDetConversionTrackSelectorTool::getEtaBin ( const Trk::Perigee & perigee ) const

private

Definition at line 62 of file InDetConversionTrackSelectorTool.cxx.

 {
   // Find the correct bin for applying eta-dependent cuts
 
   double tanThetaOver2 = std::tan( perigee.parameters()[Trk::theta] / 2.);
   double abs_eta = (tanThetaOver2 == 0) ? 999.0 : std::fabs( std::log(tanThetaOver2) );
 
   for (unsigned int i=0;i<m_TRTTrksEtaBins.size();++i) {
     if (abs_eta < m_TRTTrksEtaBins[i]) {
       return i;
     }
   }
   return m_TRTTrksEtaBins.size()-1;
 }

getPosOrBeamSpot()

Amg::Vector3D InDet::InDetConversionTrackSelectorTool::getPosOrBeamSpot ( const EventContext & ctx, const xAOD::Vertex * vertex ) const

private

Definition at line 314 of file InDetConversionTrackSelectorTool.cxx.

 {
   if (vertex) {
     return vertex->position();
   }
   if(m_useEventInfoBs){
      SG::ReadHandle<xAOD::EventInfo> evt(m_eventInfo_key, ctx);
      if (evt.isValid()) {
         return Amg::Vector3D(evt->beamPosX(), evt->beamPosY(), evt->beamPosZ());
      }else{
         return Amg::Vector3D(0, 0, 0);
      }
   }
   SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle{ m_beamSpotKey, ctx };
   if (beamSpotHandle.isValid()) {
     return beamSpotHandle->beamVtx().position();
   } else {
     return Amg::Vector3D(0, 0, 0);
   }
 }

initialize()

StatusCode InDet::InDetConversionTrackSelectorTool::initialize ( )

overridevirtual

Definition at line 28 of file InDetConversionTrackSelectorTool.cxx.

 {
   ATH_CHECK( m_extrapolator.retrieve() );
   ATH_CHECK(m_beamSpotKey.initialize(!m_useEventInfoBs));
   ATH_CHECK(m_eventInfo_key.initialize(m_useEventInfoBs));
 
   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.

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_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> InDet::InDetConversionTrackSelectorTool::m_beamSpotKey

private

Initial value:

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

Definition at line 76 of file InDetConversionTrackSelectorTool.h.

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

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_eventInfo_key

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

private

Definition at line 82 of file InDetConversionTrackSelectorTool.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::InDetConversionTrackSelectorTool::m_extrapolator

private

Initial value:

{ this,
                                                 "Extrapolator",
                                                 "",
                                                 "Extrapolator Tool" }

Definition at line 72 of file InDetConversionTrackSelectorTool.h.

                                             { this,
                                                 "Extrapolator",
                                                 "",
                                                 "Extrapolator Tool" };

m_isConv

BooleanProperty InDet::InDetConversionTrackSelectorTool::m_isConv {this, "IsConversion", true, "Conversion flag"}

private

Definition at line 120 of file InDetConversionTrackSelectorTool.h.

{this, "IsConversion", true, "Conversion flag"};

m_maxSiD0

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_maxSiD0 {this, "maxSiD0", 35., "Maximal d0 at (0,0,0) for tracks with Si hits"}

private

Properties for track selection:all cuts are ANDed.

Definition at line 86 of file InDetConversionTrackSelectorTool.h.

{this, "maxSiD0", 35., "Maximal d0 at (0,0,0) for tracks with Si hits"};

m_maxSiZ0

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_maxSiZ0 {this, "maxSiZ0", 200., "Maximal z0 at (0,0,0)"}

private

Definition at line 90 of file InDetConversionTrackSelectorTool.h.

{this, "maxSiZ0", 200., "Maximal z0 at (0,0,0)"};

m_maxTrtD0

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_maxTrtD0 {this, "maxTrtD0", 100., "Maximal d0 at (0,0,0) for standalone TRT tracks"}

private

Definition at line 88 of file InDetConversionTrackSelectorTool.h.

{this, "maxTrtD0", 100., "Maximal d0 at (0,0,0) for standalone TRT tracks"};

m_maxTrtZ0

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_maxTrtZ0 {this, "maxTrtZ0", 1200., "Maximal z0 at (0,0,0) for standalone TRT tracks"}

private

Definition at line 91 of file InDetConversionTrackSelectorTool.h.

{this, "maxTrtZ0", 1200., "Maximal z0 at (0,0,0) for standalone TRT tracks"};

m_minPt

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_minPt {this, "minPt", 500., "Minimum Pt of tracks"}

private

Definition at line 93 of file InDetConversionTrackSelectorTool.h.

{this, "minPt", 500., "Minimum Pt of tracks"};

m_PIDonlyForXe

BooleanProperty InDet::InDetConversionTrackSelectorTool::m_PIDonlyForXe {this, "PIDonlyForXe", false, "Only check TRT PID if all hits are Xe hits"}

private

Definition at line 121 of file InDetConversionTrackSelectorTool.h.

{this, "PIDonlyForXe", false, "Only check TRT PID if all hits are Xe hits"};

m_sD0_Si

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_sD0_Si {this, "significanceD0_Si", 2., "Cut on D0 significance of Si tracks"}

private

Definition at line 113 of file InDetConversionTrackSelectorTool.h.

{this, "significanceD0_Si", 2., "Cut on D0 significance of Si tracks"};

m_sD0_Trt

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_sD0_Trt {this, "significanceD0_Trt", 0.5, "Cut on D0 significance of TRT tracks"}

private

Definition at line 115 of file InDetConversionTrackSelectorTool.h.

{this, "significanceD0_Trt", 0.5, "Cut on D0 significance of TRT tracks"};

m_sZ0_Trt

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_sZ0_Trt {this, "significanceZ0_Trt", 3., "Cut on Z0 significance of TRT tracks"}

private

Definition at line 117 of file InDetConversionTrackSelectorTool.h.

{this, "significanceZ0_Trt", 3., "Cut on Z0 significance of TRT tracks"};

m_trRatio1

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_trRatio1 {this, "RatioCut1", 0.5, "TR ratio for tracks with 15 or less TRT hits"}

private

Definition at line 94 of file InDetConversionTrackSelectorTool.h.

{this, "RatioCut1", 0.5, "TR ratio for tracks with 15 or less TRT hits"};

m_trRatio2

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_trRatio2 {this, "RatioCut2", 0.1, "TR ratio for tracks with 16 to 25 TRT hits"}

private

Definition at line 96 of file InDetConversionTrackSelectorTool.h.

{this, "RatioCut2", 0.1, "TR ratio for tracks with 16 to 25 TRT hits"};

m_trRatio3

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_trRatio3 {this, "RatioCut3", 0.05, "TR ratio for tracks with 26 or more TRT hits"}

private

Definition at line 98 of file InDetConversionTrackSelectorTool.h.

{this, "RatioCut3", 0.05, "TR ratio for tracks with 26 or more TRT hits"};

m_trRatioTRT

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_trRatioTRT {this, "RatioTRT", 0.1, "TR ratio for all TRT only tracks"}

private

Definition at line 100 of file InDetConversionTrackSelectorTool.h.

{this, "RatioTRT", 0.1, "TR ratio for all TRT only tracks"};

m_trRatioV0

DoubleProperty InDet::InDetConversionTrackSelectorTool::m_trRatioV0 {this, "RatioV0", 1., "TR ratio for pion selection during V0 reconstruction"}

private

Definition at line 111 of file InDetConversionTrackSelectorTool.h.

{this, "RatioV0", 1., "TR ratio for pion selection during V0 reconstruction"};

m_TRTTrksBinnedRatioTRT

DoubleArrayProperty InDet::InDetConversionTrackSelectorTool::m_TRTTrksBinnedRatioTRT

private

Initial value:

{this, "TRTTrksBinnedRatioTRT", {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.},
     "Eta-binned eProbabilityHT for TRT-only track cuts"}

Definition at line 107 of file InDetConversionTrackSelectorTool.h.

    {this, "TRTTrksBinnedRatioTRT", {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.},
     "Eta-binned eProbabilityHT for TRT-only track cuts"};

m_TRTTrksEtaBins

DoubleArrayProperty InDet::InDetConversionTrackSelectorTool::m_TRTTrksEtaBins

private

Initial value:

{this, "TRTTrksEtaBins",
     {999., 999., 999., 999., 999., 999., 999., 999., 999., 999.},
     "Eta bins (10 expected) for TRT-only track cuts"}

Definition at line 103 of file InDetConversionTrackSelectorTool.h.

    {this, "TRTTrksEtaBins",
     {999., 999., 999., 999., 999., 999., 999., 999., 999., 999.},
     "Eta bins (10 expected) for TRT-only track cuts"};

m_useEventInfoBs

Gaudi::Property<bool> InDet::InDetConversionTrackSelectorTool::m_useEventInfoBs {this,"UseEventInfoBS",false}

private

Definition at line 83 of file InDetConversionTrackSelectorTool.h.

{this,"UseEventInfoBS",false};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• InDetConversionTrackSelectorTool.h
• InDetConversionTrackSelectorTool.cxx