Concrete implementation of the ITrackScoringTool pABC. More...

#include <InDetNNScoringTool.h>

Inheritance diagram for InDet::InDetNNScoringTool:

Collaboration diagram for InDet::InDetNNScoringTool:

Public Member Functions
	InDetNNScoringTool (const std::string &, const std::string &, const IInterface *)

virtual	~InDetNNScoringTool ()=default

virtual StatusCode	initialize () override

virtual Trk::TrackScore	score (const Trk::Track &track) const override
	create a score based on how good the passed track is More...

virtual Trk::TrackScore	simpleScore (const Trk::Track &track, const Trk::TrackSummary &trackSum) const override
	create a score based on how good the passed TrackSummary is More...

Trk::TrackScore	ambigScore (const Trk::Track &track, const Trk::TrackSummary &trackSum) const

Trk::TrackScore	calcNnScore (const Trk::Track &track, const Trk::TrackSummary &trackSum, const Trk::Perigee *extrapolatedPerigee) const

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

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

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	setupScoreModifiers ()

bool	isEmCaloCompatible (const Trk::Track &track, const EventContext &ctx) const
	Check if the cluster is compatible with a EM cluster. More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
std::unique_ptr< lwt::LightweightGraph >	m_graph

std::map< std::string, double >	m_DNN_inputValues

int	m_maxDblHoles

int	m_maxPixHoles

int	m_maxSCT_Holes

int	m_maxHits

int	m_maxSigmaChi2

int	m_maxTrtRatio

int	m_maxTrtFittedRatio

int	m_maxB_LayerHits

int	m_maxPixelHits

int	m_maxPixLay

int	m_maxGangedFakes

std::vector< double >	m_factorDblHoles

std::vector< double >	m_factorPixHoles

std::vector< double >	m_factorSCT_Holes

std::vector< double >	m_factorHits

std::vector< double >	m_factorSigmaChi2

std::vector< double >	m_factorB_LayerHits

std::vector< double >	m_factorPixelHits

std::vector< double >	m_factorPixLay

std::vector< double >	m_factorGangedFakes

std::vector< double >	m_boundsSigmaChi2

std::vector< double >	m_boundsTrtRatio

std::vector< double >	m_factorTrtRatio

std::vector< double >	m_boundsTrtFittedRatio

std::vector< double >	m_factorTrtFittedRatio

ToolHandle< ITrtDriftCircleCutTool >	m_selectortool
	Returns minimum number of expected TRT drift circles depending on eta. More...

std::vector< Trk::TrackScore >	m_summaryTypeScore
	holds the scores assigned to each Trk::SummaryType from the track's Trk::TrackSummary More...

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

ToolHandle< Trk::IExtrapolator >	m_extrapolator

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

bool	m_useAmbigFcn
	use the scoring tuned to Ambiguity processing or not More...

bool	m_useTRT_AmbigFcn

bool	m_useSigmaChi2

std::string	m_nnCutConfig

double	m_nnCutThreshold

bool	m_usePixel

bool	m_useSCT

double	m_minPt
	cuts for selecting good tracks More...

double	m_maxEta
	maximal Eta cut More...

double	m_maxRPhiImp
	maximal RPhi impact parameter cut More...

double	m_maxZImp
	maximal z impact parameter cut More...

int	m_minSiClusters
	minimal number of Si clusters More...

int	m_maxDoubleHoles

int	m_maxSiHoles
	maximal number of holes (Pixel+SCT) More...

int	m_maxPixelHoles
	maximal number of Pixel holes More...

int	m_maxSctHoles
	maximal number of SCT holes More...

int	m_minTRTonTrk
	minimum number of TRT hits More...

double	m_minTRTprecision
	minimum fraction of TRT precision hits More...

int	m_minPixel
	minimum number of pixel clusters More...

double	m_maxRPhiImpEM
	maximal RPhi impact parameter cut track that match EM clusters More...

bool	m_useEmClusSeed

float	m_phiWidthEm

float	m_etaWidthEm

SG::ReadHandleKey< ROIPhiRZContainer >	m_caloClusterROIKey {this, "EMROIPhiRZContainer", "", "Name of the calo cluster ROIs in Phi,R,Z parameterization"}

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Concrete implementation of the ITrackScoringTool pABC.

Definition at line 47 of file InDetNNScoringTool.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

◆ InDetNNScoringTool()

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

Definition at line 24 of file InDetNNScoringTool.cxx.

                                      :
   AthAlgTool(t,n,p),
   // Initialization of ScoreModifiers variables
   m_maxDblHoles(-1),
   m_maxPixHoles(-1),
   m_maxSCT_Holes(-1),
   m_maxHits(-1),
   m_maxSigmaChi2(-1),
   m_maxTrtRatio(-1),
   m_maxTrtFittedRatio(-1),
   m_maxB_LayerHits(-1),
   m_maxPixelHits(-1),
   m_maxPixLay(-1),
   m_maxGangedFakes(-1),
   m_selectortool("InDet::InDetTrtDriftCircleCutTool", this),
   m_summaryTypeScore(Trk::numberOfTrackSummaryTypes),
   m_extrapolator("Trk::Extrapolator", this)
 {
   declareInterface<Trk::ITrackScoringTool>(this);
   
   // declare properties
   declareProperty("minPt",             m_minPt              = 500.);
   declareProperty("maxEta",            m_maxEta             = 2.7);
   declareProperty("maxRPhiImp",        m_maxRPhiImp         = 10.);
   declareProperty("maxZImp",           m_maxZImp            = 250.);
  
   // changed cuts on the selction of tracks
   declareProperty("minSiClusters",     m_minSiClusters         = 7);
   declareProperty("maxDoubleHoles",    m_maxDoubleHoles        = 2);
   declareProperty("maxSiHoles",        m_maxSiHoles            = 5, "max number of Silicon (Pixel+SCT) holes");
   declareProperty("maxPixelHoles",     m_maxPixelHoles         = 5, "max number of Pixel holes");
   declareProperty("maxSCTHoles",       m_maxSctHoles           = 5, "max number of SCT holes");
   declareProperty("minTRTonTrk",       m_minTRTonTrk           = 9);
   declareProperty("minTRTPrecisionFraction", m_minTRTprecision = 0.5);
   declareProperty("minPixel",          m_minPixel              = 0);
   declareProperty("nnCutConfig",       m_nnCutConfig           = "");
   declareProperty("nnCutThreshold",    m_nnCutThreshold        = -0.1);
  
   // which detectors to use
   declareProperty("usePixel",          m_usePixel           = true);
   declareProperty("useSCT",            m_useSCT             = true);
  
   // switches and tools
   declareProperty("useAmbigFcn",       m_useAmbigFcn        = true);
   declareProperty("useTRT_AmbigFcn",   m_useTRT_AmbigFcn    = false);
   declareProperty("useSigmaChi2",      m_useSigmaChi2       = false);
  
   // tools
   declareProperty("Extrapolator",      m_extrapolator);
   declareProperty("DriftCircleCutTool",m_selectortool );
  
   declareProperty("maxRPhiImpEM",      m_maxRPhiImpEM  = 50.  );
   declareProperty("doEmCaloSeed",      m_useEmClusSeed = true );
   declareProperty("phiWidthEM",        m_phiWidthEm   = 0.075 );
   declareProperty("etaWidthEM",        m_etaWidthEm   = 0.05  );
  
   //set values for scores
   m_summaryTypeScore[Trk::numberOfPixelHits]            =  20;
   m_summaryTypeScore[Trk::numberOfPixelHoles]           = -10;  // a hole is bad
   m_summaryTypeScore[Trk::numberOfInnermostPixelLayerHits]           =  10;  // addition for being b-layer
   m_summaryTypeScore[Trk::numberOfGangedPixels]         =  -5;  // decrease for being ganged
   m_summaryTypeScore[Trk::numberOfGangedFlaggedFakes]   = -10;  // decrease for being ganged fake
   m_summaryTypeScore[Trk::numberOfSCTHits]              =  10;  // half of a pixel, since only 1dim
   m_summaryTypeScore[Trk::numberOfSCTHoles]             =  -5;  // a hole is bad !
   m_summaryTypeScore[Trk::numberOfSCTDoubleHoles]       = -15;  // a double hole is even worse !
   m_summaryTypeScore[Trk::numberOfTRTHits]              =   1;  // 10 straws ~ 1 SCT
   m_summaryTypeScore[Trk::numberOfTRTHighThresholdHits] =   0;  // addition for being TR
   m_summaryTypeScore[Trk::numberOfOutliersOnTrack]      =  -1;  // -ME- TRT oulier should not kill 5 TRT on track (was -5)
   
   // scoring for Muons not needed
   m_summaryTypeScore[Trk::numberOfMdtHits]    = 0;   
   m_summaryTypeScore[Trk::numberOfTgcPhiHits] = 0; 
   m_summaryTypeScore[Trk::numberOfTgcEtaHits] = 0;
   m_summaryTypeScore[Trk::numberOfCscPhiHits] = 0;     
   m_summaryTypeScore[Trk::numberOfCscEtaHits] = 0;
   m_summaryTypeScore[Trk::numberOfRpcPhiHits] = 0;
   m_summaryTypeScore[Trk::numberOfRpcEtaHits] = 0;
   
 }

◆ ~InDetNNScoringTool()

virtual InDet::InDetNNScoringTool::~InDetNNScoringTool ( )

virtualdefault

Member Function Documentation

◆ ambigScore()

Trk::TrackScore InDet::InDetNNScoringTool::ambigScore	(	const Trk::Track &	track,
		const Trk::TrackSummary &	trackSum
	)		const

Definition at line 375 of file InDetNNScoringTool.cxx.

 {
   //
   // --- start with bonus for high pt tracks
   //
   // double prob = 1.;
   double pt = fabs(track.trackParameters()->front()->pT());
   double prob = log10( pt ) - 1.; // 100 MeV is min and gets score 1
   ATH_MSG_DEBUG ("Modifier for pt = " << pt / 1000. << " GeV is: "<< prob);
  
   //
   // --- prob and cuts on holes
   //
   if (m_usePixel) {
     // --- Pixel holes
     int iPixHoles = trackSummary.get(Trk::numberOfPixelHoles);
     if ( iPixHoles > -1 && m_maxPixHoles > 0) {
       if (iPixHoles > m_maxPixHoles) {
         prob /= (iPixHoles - m_maxPixHoles + 1); // holes are bad ! 
         iPixHoles = m_maxPixHoles;
       }
       prob *= m_factorPixHoles[iPixHoles];
       ATH_MSG_DEBUG ("Modifier for " << iPixHoles << " Pixel holes: "<<m_factorPixHoles[iPixHoles]
          << "  New score now: " << prob);
     }  
   }
  
   if (m_useSCT) {
     // --- SCT holes
     int iSCT_Holes = trackSummary.get(Trk::numberOfSCTHoles);
     if (iSCT_Holes > -1 && m_maxSCT_Holes > 0) {
       if (iSCT_Holes > m_maxSCT_Holes) {
         prob /= (iSCT_Holes - m_maxSCT_Holes + 1); // holes are bad ! 
         iSCT_Holes = m_maxSCT_Holes;
       }
       prob *= m_factorSCT_Holes[iSCT_Holes];
       ATH_MSG_DEBUG ("Modifier for " << iSCT_Holes << " SCT holes: "<<m_factorSCT_Holes[iSCT_Holes]
          << "  New score now: " << prob);
     }
     // --- SCT double holes
     int iDblHoles = trackSummary.get(Trk::numberOfSCTDoubleHoles);
     if (iDblHoles > -1 && m_maxDblHoles > 0) {
       if (iDblHoles > m_maxDblHoles) {
         prob /= (iDblHoles - m_maxDblHoles + 1); // holes are bad !
         iDblHoles = m_maxDblHoles;
       }
       prob *= m_factorDblHoles[iDblHoles];
       ATH_MSG_DEBUG ("Modifier for " << iDblHoles << " double holes: "<<m_factorDblHoles[iDblHoles]
          << "  New score now: " << prob);
     }
   }
   //
   // --- prob for other counters
   //
   if (m_usePixel) {
  
     // ME: this if statement needs to be removed...
     // --- count layers only if holes are not searched for
     // if (trackSummary.get(Trk::numberOfPixelHoles) == -1) {
  
     // --- Pixel layers
     int iPixLay = trackSummary.get(Trk::numberOfContribPixelLayers);
     if (iPixLay > -1 && m_maxPixLay > 0) {
       if (iPixLay > m_maxPixLay) {
         prob *= (iPixLay - m_maxPixLay + 1); // layers are good !
         iPixLay = m_maxPixLay;
       }
       prob *= m_factorPixLay[iPixLay];
       ATH_MSG_DEBUG ("Modifier for " << iPixLay << " Pixel layers: "<<m_factorPixLay[iPixLay]
          << "  New score now: " << prob);
     }
  
     // --- Pixel hits
     int pixelHits = trackSummary.get(Trk::numberOfPixelHits);
     if (pixelHits > -1 && m_maxPixelHits > 0) {
       if (pixelHits > m_maxPixelHits) {
         prob *= (pixelHits - m_maxPixelHits + 1); // hits are good ! 
         pixelHits = m_maxPixelHits;
       }
       prob *= m_factorPixelHits[pixelHits];
       ATH_MSG_DEBUG ("Modifier for " << pixelHits << " Pixel hits: "<<m_factorPixelHits[pixelHits]
          << "  New score now: " << prob);
     }
     // --- Pixel blayer hits
     int bLayerHits = trackSummary.get(Trk::numberOfInnermostPixelLayerHits);
     if (bLayerHits > -1 && m_maxB_LayerHits > 0) {
       if (bLayerHits > m_maxB_LayerHits) {
         prob *= (bLayerHits - m_maxB_LayerHits + 1); // hits are good !
         bLayerHits = m_maxB_LayerHits;
       }
       prob *= m_factorB_LayerHits[bLayerHits];
       ATH_MSG_DEBUG ("Modifier for " << bLayerHits << " b-layer hits: "<<m_factorB_LayerHits[bLayerHits]
          << "  New score now: " << prob);
     }
     // --- Pixel Ganged fakes
     int pixelGangedFakes = trackSummary.get(Trk::numberOfGangedFlaggedFakes);
     if (pixelGangedFakes > -1 && m_maxGangedFakes > 0) {
       if (pixelGangedFakes > m_maxGangedFakes) {
         prob /= (pixelGangedFakes - m_maxGangedFakes + 1); // ganged are bad !
         pixelGangedFakes = m_maxGangedFakes;
       }
       prob *= m_factorGangedFakes[pixelGangedFakes];
       ATH_MSG_DEBUG ("Modifier for " << pixelGangedFakes << " ganged fakes hits: "<<m_factorGangedFakes[pixelGangedFakes]
          << "  New score now: " << prob);
     }
   }
  
   int iHits = 0;
   iHits += m_usePixel ? trackSummary.get(Trk::numberOfPixelHits) : 3;  // if Pixel off, do not count as inefficient
   iHits += m_useSCT   ? trackSummary.get(Trk::numberOfSCTHits)   : 8;  // if SCT off, do not count as inefficient
   if (iHits > -1 && m_maxHits > 0) {
     if (iHits > m_maxHits) {
       prob *= (iHits - m_maxHits + 1); // hits are good !
       iHits = m_maxHits; 
     }
     prob *= m_factorHits[iHits]; 
     ATH_MSG_DEBUG ("Modifier for " << iHits << " Sihits: "<<m_factorHits[iHits]
        << "  New score now: " << prob);
   }
  
   //
   // --- special treatment for TRT hits
   //
   int iTRT_Hits     = trackSummary.get(Trk::numberOfTRTHits);
   int iTRT_Outliers = trackSummary.get(Trk::numberOfTRTOutliers);
   //
   if ( iTRT_Hits > 0 && m_maxTrtRatio > 0) {
     // get expected number of TRT hits
     double nTrtExpected = 30.;
     assert( m_selectortool.isEnabled() );
     nTrtExpected = m_selectortool->minNumberDCs(track.trackParameters()->front());
     ATH_MSG_DEBUG ("Expected number of TRT hits: " << nTrtExpected << " for eta: "
        << fabs(track.trackParameters()->front()->eta()));
     double ratio = (nTrtExpected != 0) ? iTRT_Hits / nTrtExpected : 0;
     if (ratio > m_boundsTrtRatio[m_maxTrtRatio]) ratio = m_boundsTrtRatio[m_maxTrtRatio];
     for (int i=0; i<m_maxTrtRatio; ++i) {
       if ( m_boundsTrtRatio[i] < ratio && ratio <= m_boundsTrtRatio[i+1]) {
         prob*= m_factorTrtRatio[i];
         ATH_MSG_DEBUG ("Modifier for " << iTRT_Hits << " TRT hits (ratio " << ratio
                  << ") is : "<< m_factorTrtRatio[i] << "  New score now: " << prob);
         break;
       }
     } 
   }
   //
   if ( iTRT_Hits > 0 && iTRT_Outliers >= 0 && m_maxTrtFittedRatio > 0) {
     double fitted = double(iTRT_Hits) / double(iTRT_Hits + iTRT_Outliers);
     if (fitted > m_boundsTrtFittedRatio[m_maxTrtFittedRatio]) fitted = m_boundsTrtFittedRatio[m_maxTrtFittedRatio];
     for (int i=0; i<m_maxTrtFittedRatio; ++i) {
       if (fitted <= m_boundsTrtFittedRatio[i+1]) {
         prob*= m_factorTrtFittedRatio[i];
         ATH_MSG_DEBUG ("Modifier for TRT fitted ratio of " << fitted
                  << " is : "<< m_factorTrtFittedRatio[i] << "  New score now: " << prob);
         break;
       }
     } 
   }
  
   // is this a track from the pattern or a fitted track ?
   bool ispatterntrack = (track.info().trackFitter()==Trk::TrackInfo::Unknown);
  
   // 
   // --- non binned Chi2
   //
   if (!ispatterntrack) {
     if (track.fitQuality()!=nullptr && track.fitQuality()->chiSquared()>0 && track.fitQuality()->numberDoF()>0 ) {
       int    indf  = track.fitQuality()->numberDoF();
       double chi2  = track.fitQuality()->chiSquared();
       double fac   = 1. / log10 (10. + 10. * chi2 / indf); // very soft chi2 
       prob        *= fac;
       ATH_MSG_DEBUG ("Modifier for chi2 = " << chi2 << " and NDF = " << indf
          << " is : "<< fac << "  New score now: " << prob);
  
     }
   }
   //
   // --- fit quality prob
   //
   if ( !ispatterntrack && (m_useSigmaChi2) && track.fitQuality() ) {
  
     int    ndf  = track.fitQuality()->numberDoF();
     double chi2 = track.fitQuality()->chiSquared();
     if (ndf > 0) {
       //
       // --- special variable for bad chi2 distribution
       //
       if (m_useSigmaChi2) {
         int sigmaChi2times100 = trackSummary.get(Trk::standardDeviationOfChi2OS);
         if (sigmaChi2times100 > 0) {
           double  testvar = double(sigmaChi2times100)/100. - sqrt(2.*chi2/ndf);
           if (msgLvl(MSG::VERBOSE)) msg(MSG::VERBOSE) << "sigma chi2 = " << testvar << endmsg;
           if ( testvar< m_boundsSigmaChi2[0] ) {
             prob *= m_factorSigmaChi2[0];
             ATH_MSG_DEBUG ("Modifier for " << testvar << " sigma chi2: "<< m_factorSigmaChi2[0]
                << "  New score now: " << prob);
           } else if (m_boundsSigmaChi2[m_maxSigmaChi2] <= testvar) {
             prob *= m_factorSigmaChi2[m_maxSigmaChi2-1];
             ATH_MSG_DEBUG ("Modifier for " << testvar << " sigma chi2: "<< m_factorSigmaChi2[m_maxSigmaChi2-1]
                << "  New score now: " << prob);
           } else {
             for (int i = 0 ; i<m_maxSigmaChi2 ; ++i ) {
               if ( m_boundsSigmaChi2[i]<=testvar && testvar<m_boundsSigmaChi2[i+1] ) {
           prob *= m_factorSigmaChi2[i];
           ATH_MSG_DEBUG ("Modifier for " << testvar << " sigma chi2: "<< m_factorSigmaChi2[i]
                    << "  New score now: " << prob);
           break;
               }
             }
           }
         }
       }
     }
   }
  
   return Trk::TrackScore(prob); 
 }

◆ calcNnScore()

Trk::TrackScore InDet::InDetNNScoringTool::calcNnScore	(	const Trk::Track &	track,
		const Trk::TrackSummary &	trackSum,
		const Trk::Perigee *	extrapolatedPerigee
	)		const

Definition at line 848 of file InDetNNScoringTool.cxx.

 {
   ATH_MSG_DEBUG("Using NN Score Function");
   // initialize with dummy score
   double DNNscore(-1.0);
  
   // This calculates a variant of the delta-eta variable used in large-d0 seeding
   double d0 = extrapolatedPerigee->parameters()[Trk::d0];
   double z0 = extrapolatedPerigee->parameters()[Trk::z0];
   double deltaEta = std::abs(std::atan2(fabs(d0), z0) - 2 * std::atan(std::exp(-track.trackParameters()->front()->eta())));
  
   // Build dictionary of inputs for lwtnn to use
   // It is ok to fill this with more variables than the model uses
   // as long as no variables are missing
   std::map<std::string, double> trackInputs{
           {"pT",               track.trackParameters()->front()->pT()},
           {"eta",              track.trackParameters()->front()->eta()},
           {"numberOfSCTHoles", (double) trackSummary.get(Trk::numberOfSCTHoles)},
           {"numberOfSCTHits",  (double) trackSummary.get(Trk::numberOfSCTHits)},
           {"numberDoF",        (double) track.fitQuality()->numberDoF()},
           {"Sihits",           (double) (trackSummary.get(Trk::numberOfPixelHits) + trackSummary.get(Trk::numberOfSCTHits))},
           {"d0",               d0},
           {"z0",               z0},
           {"deltaEta",         deltaEta}
   };
  
   // Set up the nodes used for inputs
   std::map<std::string, std::map<std::string, double> > inputs{
           {"trackInputs", trackInputs}
   };
  
   // Evaluate the network
   lwt::ValueMap output = m_graph->compute(inputs);
  
   // Obtain the discriminant associated with the single output node
   DNNscore = output["nnScore"];
   // Return the discriminant as the score
   ATH_MSG_DEBUG("DNNscore: " << DNNscore);
   if (DNNscore < m_nnCutThreshold)
     ATH_MSG_DEBUG("DNNscore is below threshold of " << m_nnCutThreshold << ", rejecting track. ");
   if (DNNscore < 0)
     ATH_MSG_ERROR("DNNscore should be 0 or greater.");
   return Trk::TrackScore(DNNscore);
 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                  {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

◆ initialize()

StatusCode InDet::InDetNNScoringTool::initialize ( )

overridevirtual

Definition at line 108 of file InDetNNScoringTool.cxx.

 {
   StatusCode sc = AlgTool::initialize();
   if (sc.isFailure()) return sc;
   
   sc = m_extrapolator.retrieve();
   if (sc.isFailure()) {
     msg(MSG::FATAL) << "Failed to retrieve tool " << m_extrapolator << endmsg;
     return StatusCode::FAILURE;
   } else 
     msg(MSG::DEBUG) << "Retrieved tool " << m_extrapolator << endmsg;
  
   ATH_CHECK (m_selectortool.retrieve( DisableTool{ m_selectortool.empty() } ));
   // Get segment selector tool
   //
  
   ATH_CHECK(m_beamSpotKey.initialize());
  
   if (m_useAmbigFcn && m_useTRT_AmbigFcn) {
     msg(MSG::FATAL) << "Both on, normal ambi funciton and the one for back tracking, configuration problem, not recoverable" << endmsg;
     return StatusCode::FAILURE;
   }
   
   // Read handle for AtlasFieldCacheCondObj
   ATH_CHECK( m_fieldCacheCondObjInputKey.initialize() );
   
   if (m_useAmbigFcn || m_useTRT_AmbigFcn) setupScoreModifiers();
  
   // lwtnn initialization
   if (m_nnCutThreshold > 0.0) { // Load NN only if it will be used
       // Locate configuration file  
       std::string nnCutConfigPath = PathResolverFindCalibFile(m_nnCutConfig); //returns "" if file not found
       if (nnCutConfigPath.empty()){
         ATH_MSG_FATAL ( "Failed to find configuration file: " << m_nnCutConfig);
         return StatusCode::FAILURE;
       }
       ATH_MSG_DEBUG ("Loading configuration file:" << nnCutConfigPath);
       // Load NN weights
       std::ifstream inFile(nnCutConfigPath);
       lwt::GraphConfig config(lwt::parse_json_graph(inFile));
       m_graph = std::make_unique<lwt::LightweightGraph>(config);
       ATH_MSG_DEBUG ("NN configuration loaded");
       inFile.close();
   }
  
   ATH_CHECK( m_caloClusterROIKey.initialize(m_useEmClusSeed) );
   
   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::ITrackScoringTool::interfaceID ( )

inlinestaticinherited

Definition at line 36 of file ITrackScoringTool.h.

 {
   return IID_ITrackScoringTool;
 }

◆ isEmCaloCompatible()

bool InDet::InDetNNScoringTool::isEmCaloCompatible	(	const Trk::Track &	track,
		const EventContext &	ctx
	)		const

private

Check if the cluster is compatible with a EM cluster.

Definition at line 820 of file InDetNNScoringTool.cxx.

 {
   const Trk::TrackParameters * Tp = track.trackParameters()->front();
  
   //Switch to the track parameters of the first measurment instead of the perigee parameters 
   ATH_MSG_VERBOSE ("--> Looping over TSOS's");
   for (const auto *tsos : *track.trackStateOnSurfaces() ) {
     // get measurment from TSOS
     const auto *meas = tsos->measurementOnTrack();
     const auto *tp = tsos->trackParameters();
     // if we do not have a measurement, we should just mark it
     if (!meas || !tp) {
       continue;
     } else {
       Tp = tp;
       break;
     }
   }
  
   double F = Tp->momentum().phi();
   double E = Tp->momentum().eta();
   double R = Tp->position().perp();
   double Z = Tp->position().z();
  
   SG::ReadHandle<ROIPhiRZContainer> calo(m_caloClusterROIKey, ctx);
   return calo->hasMatchingROI(F, E,  R, Z, m_phiWidthEm, m_etaWidthEm);
 }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ score()

Trk::TrackScore InDet::InDetNNScoringTool::score ( const Trk::Track & track ) const

overridevirtual

create a score based on how good the passed track is

Implements Trk::ITrackScoringTool.

Definition at line 160 of file InDetNNScoringTool.cxx.

 {
    if (!track.trackSummary()) {
       ATH_MSG_FATAL("Track without a summary");
    }
    ATH_MSG_VERBOSE ("Track has TrackSummary "<<*track.trackSummary());
    Trk::TrackScore score = Trk::TrackScore( simpleScore(track, *track.trackSummary()) );
    ATH_MSG_DEBUG ("Track has Score: "<<score);
  
    return score;
 }

◆ setupScoreModifiers()

void InDet::InDetNNScoringTool::setupScoreModifiers ( )

private

Definition at line 594 of file InDetNNScoringTool.cxx.

 {
   //
   // --- number of Pixel holes
   //
   // --- NewTracking and BackTracking
   const int maxPixHoles = 2; // there is an explicit cut
   const double goodPixHoles[maxPixHoles+1] = {1.0 , 0.04 , 0.004 };
   const double fakePixHoles[maxPixHoles+1] = {1.0 , 0.30 , 0.200 };
   // put it into the private members
   m_maxPixHoles = maxPixHoles;
   for (int i=0; i<=m_maxPixHoles; ++i) m_factorPixHoles.push_back(goodPixHoles[i]/fakePixHoles[i]);
  
   //
   // --- number of SCT holes
   //
   if (!m_useTRT_AmbigFcn) {
     // --- NewTracking
     const int maxSCT_Holes = 5; // moved from 3 -> 5 , there is an explicit cut anyway
     const double goodSCT_Holes[maxSCT_Holes+1] = { 1.0 , 0.06 , 0.010 , 0.0007, 0.0005, 0.0003  };
     const double fakeSCT_Holes[maxSCT_Holes+1] = { 1.0 , 0.15 , 0.100 , 0.0100, 0.0100, 0.0100 };
     // put it into the private members
     m_maxSCT_Holes = maxSCT_Holes;
     for (int i=0; i<=m_maxSCT_Holes; ++i) m_factorSCT_Holes.push_back(goodSCT_Holes[i]/fakeSCT_Holes[i]);
   } else {
     // --- BackTracking
     const int maxSCT_Holes = 6;
     const double goodSCT_Holes[maxSCT_Holes+1] = {0.910, 0.074, 0.014, 0.001, 0.001, 0.00001, 0.00001};
     const double fakeSCT_Holes[maxSCT_Holes+1] = {0.910, 0.192, 0.229, 0.061, 0.065, 0.016  , 0.025};
     // put it into the private members
     m_maxSCT_Holes = maxSCT_Holes;
     for (int i=0; i<=m_maxSCT_Holes; ++i) m_factorSCT_Holes.push_back(goodSCT_Holes[i]/fakeSCT_Holes[i]);
   }
  
   //
   // --- number of SCT double holes
   //
   // --- NewTracking and BackTracking
   const int maxDblHoles = 3; // there is a cut on this anyway !
   const double goodDblHoles[maxDblHoles+1] = { 1.  , 0.03 , 0.007 , 0.0003  };
   const double fakeDblHoles[maxDblHoles+1] = { 1.  , 0.09 , 0.09  , 0.008  };
   // put it into the private members
   m_maxDblHoles = maxDblHoles;
   for (int i=0; i<=m_maxDblHoles; ++i) m_factorDblHoles.push_back(goodDblHoles[i]/fakeDblHoles[i]);
  
   //
   // --- number of Blayer hits
   //
   if (!m_useTRT_AmbigFcn) {
     // --- NewTracking
     const int maxB_LayerHits = 3;
     const double goodB_LayerHits[maxB_LayerHits+1] = {0.203, 0.732, 0.081, 0.010};
     const double fakeB_LayerHits[maxB_LayerHits+1] = {0.808, 0.174, 0.018, 0.002};
     // put it into the private members
     m_maxB_LayerHits = maxB_LayerHits;
     for (int i=0; i<=m_maxB_LayerHits; ++i) m_factorB_LayerHits.push_back(goodB_LayerHits[i]/fakeB_LayerHits[i]);
   } else {
     // --- BackTracking
     const int maxB_LayerHits = 3;
     const double goodB_LayerHits[maxB_LayerHits+1] = {0.605, 0.349, 0.044, 0.010};
     const double fakeB_LayerHits[maxB_LayerHits+1] = {0.865, 0.124, 0.011, 0.002};
     // put it into the private members
     m_maxB_LayerHits = maxB_LayerHits;
     for (int i=0; i<=m_maxB_LayerHits; ++i) m_factorB_LayerHits.push_back(goodB_LayerHits[i]/fakeB_LayerHits[i]);
   }
  
   //
   // --- number of Pixel hits without Blayer
   //
   if (!m_useTRT_AmbigFcn) {
     // --- NewTracking
     const int maxPixelHits = 8; // we see up to 8 with IBL (was 6)
     const double goodPixelHits[maxPixelHits+1] = {0.095, 0.031, 0.118, 0.615, 0.137, 0.011, 0.01   , 0.011, 0.012};
     const double fakePixelHits[maxPixelHits+1] = {0.658, 0.100, 0.091, 0.124, 0.026, 0.002, 0.001  , 0.001, 0.001};
     m_maxPixelHits = maxPixelHits;
     for (int i=0; i<=m_maxPixelHits; ++i) m_factorPixelHits.push_back(goodPixelHits[i]/fakePixelHits[i]);
   } else {
     // --- BackTracking
     const int maxPixelHits = 8; // we see up to 8 with IBL (was 6)
     const double goodPixelHits[maxPixelHits+1] = {0.401, 0.079, 0.140, 0.291, 0.011, 0.078, 0.01   , 0.011, 0.012};
     const double fakePixelHits[maxPixelHits+1] = {0.673, 0.138, 0.113, 0.057, 0.002, 0.011, 0.001  , 0.001, 0.001}; 
     m_maxPixelHits = maxPixelHits;
     for (int i=0; i<=m_maxPixelHits; ++i) m_factorPixelHits.push_back(goodPixelHits[i]/fakePixelHits[i]);
   }
  
   //
   // --- number of Pixel layers
   //
   if (!m_useTRT_AmbigFcn) {
     // --- NewTracking
     const int maxPixLay = 7; // 3 barrel, 3 endcap, IBL, in practice one should see maybe 5 (was 3)
     const double goodPixLay[maxPixLay+1] = {0.095, 0.033, 0.131, 0.740,   0.840, 0.940, 1.040,1.140};
     const double fakePixLay[maxPixLay+1] = {0.658, 0.106, 0.092, 0.144,   0.144, 0.144, 0.144,0.144};
     // put it into the private members
     m_maxPixLay = maxPixLay;
     for (int i=0; i<=m_maxPixLay; ++i) m_factorPixLay.push_back(goodPixLay[i]/fakePixLay[i]);
   } else {
     // --- BackTracking
     const int maxPixLay = 7; // 3 barrel, 3 endcap, IBL, in practice one should see maybe 5 (was 5)
     const double goodPixLay[maxPixLay+1] = {0.401, 0.088, 0.152, 0.355,   0.405, 0.455, 0.505, 0.555};
     const double fakePixLay[maxPixLay+1] = {0.673, 0.146, 0.115, 0.065,   0.065, 0.065, 0.065, 0.065};
     // put it into the private members
     m_maxPixLay = maxPixLay;
     for (int i=0; i<=m_maxPixLay; ++i) m_factorPixLay.push_back(goodPixLay[i]/fakePixLay[i]);
   }
   
   //
   // --- number of Pixel Ganged Fakes
   //
   // --- NewTracking and BackTracking
   const int maxGangedFakes = 2; // there is an explicit cut
   const double goodGangedFakes[maxGangedFakes+1] = {0.62 , 0.23 , 0.15 };
   const double fakeGangedFakes[maxGangedFakes+1] = {0.12 , 0.41 , 0.47 };
   // put it into the private members
   m_maxGangedFakes = maxGangedFakes;
   for (int i=0; i<=m_maxGangedFakes; ++i) m_factorGangedFakes.push_back(goodGangedFakes[i]/fakeGangedFakes[i]);
  
   //
   // --- total number of SCT+Pixel hits
   //
   // --- NewTracking and BackTracking
   const int maxHits = 19; // there is a min cut on this anyway !
   const double goodSiHits[maxHits+1] = { 0.001 , 0.002 , 0.003 , 0.004 , 0.01 , 0.01 , 0.01 ,
            0.015 , 0.02 , 0.06 , 0.1  , 0.3  , 0.2   , 0.50 , 0.055,
            0.03  , 0.015  , 0.010  , 0.002  , 0.0005 };
   const double fakeSiHits[maxHits+1] = { 1.0   , 1.0   , 1.0   , 1.0   , 0.5 , 0.25  , 0.15 ,
            0.20  , 0.1  , 0.2  , 0.08 , 0.07 , 0.035 , 0.08 , 0.008,
            0.004 , 0.0015 , 0.0008 , 0.0001 , 0.00001 };
   // put it into the private members
   m_maxHits = maxHits;
   for (int i=0; i<=m_maxHits; ++i) m_factorHits.push_back(goodSiHits[i]/fakeSiHits[i]);
  
   //
   // --- ratio of TRT hits over expected
   //
   // --- NewTracking and BackTracking
   const int maxTrtRatio = 7 ;
   const double TrtRatioBounds[maxTrtRatio+1] = {  0,   0.2,  0.4,  0.6,  0.8,  1.0,  1.2, 2.4};
   // this needs tuning !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   const double goodTrtRatio[maxTrtRatio]     = {      0.05, 0.11, 0.12, 0.15, 0.20, 0.16, 0.17};
   const double fakeTrtRatio[maxTrtRatio]     = {      0.6 , 0.08, 0.06, 0.05, 0.04, 0.03, 0.03};
   // put it into the private members
   m_maxTrtRatio = m_selectortool.isEnabled() ? maxTrtRatio : 0;
   for (int i=0; i<m_maxTrtRatio; ++i) m_factorTrtRatio.push_back(goodTrtRatio[i]/fakeTrtRatio[i]);
   for (int i=0; i<=m_maxTrtRatio; ++i) m_boundsTrtRatio.push_back(TrtRatioBounds[i]);
  
   //
   // --- ratio of TRT fitted to (fitted+outliers)
   //
   // --- NewTracking and BackTracking
   const int maxTrtFittedRatio = 4;
   const double TrtFittedRatioBounds[maxTrtFittedRatio+1] = {  0,   0.3, 0.6, 0.9, 1.0};
   // this needs tuning !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   const double goodTrtFittedRatio[maxTrtFittedRatio]     = {       0.1, 0.2, 0.3, 0.5};
   const double fakeTrtFittedRatio[maxTrtFittedRatio]     = {       0.6, 0.1, 0.1, 0.1};
   // put it into the private members
   m_maxTrtFittedRatio = maxTrtFittedRatio;
   for (int i=0; i<m_maxTrtFittedRatio; ++i) m_factorTrtFittedRatio.push_back(goodTrtFittedRatio[i]/fakeTrtFittedRatio[i]);
   for (int i=0; i<=m_maxTrtFittedRatio; ++i) m_boundsTrtFittedRatio.push_back(TrtFittedRatioBounds[i]);
  
  
   //
   // --- sigma chi2
   //
   if (!m_useSigmaChi2) {
     m_maxSigmaChi2 = -1 ;
   } else {
     // --- NewTracking and BackTracking
     const int maxSigmaChi2 = 13 ;
     const double SigmaChi2Bounds[maxSigmaChi2+1] = { -5.,-4.,-3.,-2.,-1.,0.,1.,2.,3.,4.,5.,6.,7.,8.};
     const double goodSigmaChi2[maxSigmaChi2] = {0.00004, 0.0004, 0.002, 0.15, 0.8, 0.015, 0.01 , 0.009, 0.008, 0.0007 , 0.0006 , 0.0005, 0.00004};
     const double fakeSigmaChi2[maxSigmaChi2] = {0.0008 , 0.005 , 0.02 , 0.2 , 0.3, 0.1  , 0.1  , 0.1  , 0.1  , 0.01   , 0.01   , 0.01  , 0.001};
     // put it into the private members
     m_maxSigmaChi2 = maxSigmaChi2;
     for (int i=0; i<m_maxSigmaChi2; ++i) m_factorSigmaChi2.push_back(goodSigmaChi2[i]/fakeSigmaChi2[i]);
     for (int i=0; i<=m_maxSigmaChi2; ++i) m_boundsSigmaChi2.push_back(SigmaChi2Bounds[i]);
   }
  
   //
   // --- debug output
   //
   if (msgLvl(MSG::VERBOSE)) { 
     
     for (int i=0; i<=m_maxPixHoles; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " Pixel holes: " << m_factorPixHoles[i] <<endmsg;
     
     for (int i=0; i<=m_maxSCT_Holes; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " SCT holes: " << m_factorSCT_Holes[i] <<endmsg;
     
     for (int i=0; i<=m_maxDblHoles; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " double SCT holes: " << m_factorDblHoles[i] <<endmsg;
     
     for (int i=0; i<=m_maxPixLay; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " Pixel layers: " << m_factorPixLay[i] <<endmsg;
     
     for (int i=0; i<=m_maxB_LayerHits; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " b-layer hits: " << m_factorB_LayerHits[i] <<endmsg;
     
     for (int i=0; i<=m_maxPixelHits; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " Pixel hits: " << m_factorPixelHits[i] <<endmsg;
     
     for (int i=0; i<=m_maxGangedFakes; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " ganged fakes: " << m_factorGangedFakes[i] <<endmsg;
     
     for (int i=0; i<=m_maxHits; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << i << " Si hits: " << m_factorHits[i] <<endmsg;
  
     if (m_selectortool.isEnabled()) {
        for (int i=0; i<m_maxTrtRatio; ++i)
           msg(MSG::VERBOSE) << "Modifier for " << m_boundsTrtRatio[i] << " < TRT ratio  < "
                             << m_boundsTrtRatio[i+1] <<"  : " <<m_factorTrtRatio[i] <<endmsg;
     }
  
     for (int i=0; i<m_maxTrtFittedRatio; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << m_boundsTrtFittedRatio[i] << " < TRT fitted ratio  < "
       << m_boundsTrtFittedRatio[i+1] <<"  : " <<m_factorTrtFittedRatio[i] <<endmsg;
     
     
     // only if used
     for (int i=0; i<m_maxSigmaChi2; ++i)
       msg(MSG::VERBOSE) << "Modifier for " << m_boundsSigmaChi2[i] << " < sigma(chi2) - sqrt(2chi2)  < " << m_boundsSigmaChi2[i+1]
       <<"  : " <<m_factorSigmaChi2[i] <<endmsg;
   }
  
 }

◆ simpleScore()

Trk::TrackScore InDet::InDetNNScoringTool::simpleScore	(	const Trk::Track &	track,
		const Trk::TrackSummary &	trackSum
	)		const

overridevirtual

create a score based on how good the passed TrackSummary is

Implements Trk::ITrackScoringTool.

Definition at line 174 of file InDetNNScoringTool.cxx.

 {
   int numPixel          = trackSummary.get(Trk::numberOfPixelHits);
   int numSCT            = trackSummary.get(Trk::numberOfSCTHits);
   int numTRT            = trackSummary.get(Trk::numberOfTRTHits);
   int numTRTTube        = trackSummary.get(Trk::numberOfTRTTubeHits);
   int numPixelHoles     = trackSummary.get(Trk::numberOfPixelHoles);
   int numSCTHoles       = trackSummary.get(Trk::numberOfSCTHoles);
   int numSCTDoubleHoles = trackSummary.get(Trk::numberOfSCTDoubleHoles);
   int numPixelDead      = trackSummary.get(Trk::numberOfPixelDeadSensors);  
   int numSCTDead        = trackSummary.get(Trk::numberOfSCTDeadSensors);  
  
   if (numPixelDead<0) numPixelDead = 0;
   if (numSCTDead<0)   numSCTDead   = 0;
   
   // is this a track from the pattern or a fitted track ?
   bool ispatterntrack = (track.info().trackFitter()==Trk::TrackInfo::Unknown);
   //
   // --- reject bad tracks
   //
   //
   // --- the following cuts we only apply to fitted tracks
   //
   if (ispatterntrack) {
     ATH_MSG_DEBUG ("==> this is a pattern track, no hit cuts !");
   } else {
     ATH_MSG_DEBUG ("==> this is a refitted track, so we can use the chi2 and other stuff ! ");
  
     // NdF cut :
     if (track.fitQuality()) {
       ATH_MSG_DEBUG ("numberDoF = "<<track.fitQuality()->numberDoF());
     }
     // Number of double Holes
     if (numSCTDoubleHoles>=0) {
       if (numSCTDoubleHoles > m_maxDoubleHoles ) {
         ATH_MSG_DEBUG ("Track has "<< numSCTDoubleHoles <<" double holes, reject it!");
         return Trk::TrackScore(0);
       }
     }
     // Number of Si (Pixel+SCT) Holes
     if (numSCTHoles>=0 && numPixelHoles>=0) {
       if (numPixelHoles+numSCTHoles > m_maxSiHoles ) {
         ATH_MSG_DEBUG ("Track has "<< numPixelHoles <<" Pixel and " << numSCTHoles << " SCT holes, reject it!");
         return Trk::TrackScore(0);
       }
     }
     // Number of Pixel Holes
     if ( numPixelHoles>=0 ) {
       if (numPixelHoles > m_maxPixelHoles ) {
         ATH_MSG_DEBUG ("Track has "<< numPixelHoles <<" Pixel  holes, reject it!");
         return Trk::TrackScore(0);
       }
     }
     // Number of SCT Holes
     if ( numSCTHoles>=0 ) {
       if ( numSCTHoles > m_maxSctHoles ) {
         ATH_MSG_DEBUG ("Track has "<< numSCTHoles << " SCT holes, reject it!");
         return Trk::TrackScore(0);
       }
     }
     // TRT cut
     if (numTRT > 0 && numTRT < m_minTRTonTrk) {
       ATH_MSG_DEBUG ("Track has " << numTRT << " TRT hits,  reject it");
       return Trk::TrackScore(0);
     }
     // TRT precision hits cut  
     if (numTRT >= 15 && ((double)(numTRT-numTRTTube))/numTRT < m_minTRTprecision) {  
       ATH_MSG_DEBUG ("Track has " << ((double)numTRTTube)/numTRT << " TRT tube hit fraction,  reject it");
       return Trk::TrackScore(0);  
     }  
     // Number of Si Clusters
     if ( numSCT>=0 && numPixel>=0) {
       if (numPixel+numSCT+numPixelDead+numSCTDead < m_minSiClusters) {
         ATH_MSG_DEBUG ("Track has " << numPixel+numSCT << " Si clusters and " << numPixelDead+numSCTDead << " dead sensors, reject it");
         return Trk::TrackScore(0);
       }
     }
     // Number of pixel clusters
     if (numPixel>=0 && numPixel < m_minPixel) {
       ATH_MSG_DEBUG ("Track has " << numPixel << " pixel hits, reject it");
       return Trk::TrackScore(0);
     }
   }
   //
   // --- kinematic selection (done as well on input ?)
   //
   // --- beam spot position 
   Amg::Vector3D beamSpotPosition(0,0,0);
   SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey };
   if (beamSpotHandle.isValid()) beamSpotPosition = beamSpotHandle->beamVtx().position();
   // --- create surface
   Trk::PerigeeSurface perigeeSurface(beamSpotPosition);
  
   const Trk::TrackParameters* input = track.trackParameters()->front();
  
   // cuts on parameters
   const EventContext& ctx = Gaudi::Hive::currentContext();
   SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCacheCondObjInputKey, ctx};
   const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
   if (fieldCondObj == nullptr) {
       ATH_MSG_ERROR("simpleScore: Failed to retrieve AtlasFieldCacheCondObj with key " << m_fieldCacheCondObjInputKey.key());
       return Trk::TrackScore(0);
   }
   MagField::AtlasFieldCache fieldCache;
   fieldCondObj->getInitializedCache (fieldCache);
  
   if (fieldCache.solenoidOn()){ 
     if (fabs(input->pT()) < m_minPt) {
       ATH_MSG_DEBUG ("Track pt < "<<m_minPt<<", reject it");
       return Trk::TrackScore(0);
     } 
   }
   if (fabs(input->eta()) > m_maxEta) {
     ATH_MSG_DEBUG ("Track eta > "<<m_maxEta<<", reject it");
     return Trk::TrackScore(0);
   }
  
   // uses perigee on track or extrapolates, no material in any case, we cut on impacts
   // add back extrapolation without errors
   {
   std::unique_ptr<const Trk::TrackParameters> parm( m_extrapolator->extrapolateDirectly(ctx, *input, perigeeSurface) );
  
   const Trk::Perigee*extrapolatedPerigee = dynamic_cast<const Trk::Perigee*> (parm.get());
   if (!extrapolatedPerigee) {
      ATH_MSG_WARNING( "Extrapolation of perigee failed, this should never happen" );
      return Trk::TrackScore(0);
   }
  
   ATH_MSG_VERBOSE ("extrapolated perigee: "<<*extrapolatedPerigee);
   if (fabs(extrapolatedPerigee->parameters()[Trk::z0]) > m_maxZImp) {
     ATH_MSG_DEBUG ("Track Z impact > "<<m_maxZImp<<", reject it");
     return Trk::TrackScore(0);
   }
  
   double maxD0 = m_maxRPhiImp;
   if(m_useEmClusSeed) {
      if (isEmCaloCompatible( track, ctx ) ) {
         maxD0 = m_maxRPhiImpEM;
      }
   }
   if (fabs(extrapolatedPerigee->parameters()[Trk::d0]) > maxD0) {
     ATH_MSG_DEBUG ("Track Rphi impact > "<<maxD0<<", reject it");
     return Trk::TrackScore(0);
   }
  
   if (m_useAmbigFcn) {
     if (m_nnCutThreshold > 0.0) { // check if nn cut is enabled
       if (track.fitQuality()) { // check if track has been fitted
         if (calcNnScore(track, trackSummary, extrapolatedPerigee) < m_nnCutThreshold) {
           ATH_MSG_DEBUG ("Rejecting track for falling below nn threshold");
           return Trk::TrackScore(0);  // scores of 0 are rejected by ambisolver
         }
       }
     }
   }
   }
  
   //
   // --- now start scoring
   //
  
   if ( m_useAmbigFcn || m_useTRT_AmbigFcn) {
     
     //
     // --- use scoring function
     //
     
     return ambigScore(track,trackSummary);
  
   } else {
   
     //
     // use classical score
     //
     // score of 100 per track
     Trk::TrackScore score (100);
     // --- summary score analysis
     for (int i=0; i<Trk::numberOfTrackSummaryTypes; ++i) {
       int value = trackSummary.get(static_cast<Trk::SummaryType>(i));
       //value is -1 if undefined.
       if (value>0) { 
         score+=m_summaryTypeScore[i]*value; 
         ATH_MSG_DEBUG ("\tType ["<<i<<"], value \t= "<<value<<"], score \t="<<score);
       }
     }
     // --- prob(chi2,NDF), protect for chi2 <= 0
     if (track.fitQuality()!=nullptr && track.fitQuality()->chiSquared()>0 && track.fitQuality()->numberDoF()>0 ) {
       double p = 1.0-Genfun::CumulativeChiSquare(track.fitQuality()->numberDoF())(track.fitQuality()->chiSquared());
       if ( p > 0 )
         score += log10( p );
       else
         score -= 50;
     }    
     
     return score;
   }
  
 }

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

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

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_beamSpotKey

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

private

Definition at line 89 of file InDetNNScoringTool.h.

◆ m_boundsSigmaChi2

std::vector<double> InDet::InDetNNScoringTool::m_boundsSigmaChi2

private

Definition at line 80 of file InDetNNScoringTool.h.

◆ m_boundsTrtFittedRatio

std::vector<double> InDet::InDetNNScoringTool::m_boundsTrtFittedRatio

private

Definition at line 81 of file InDetNNScoringTool.h.

◆ m_boundsTrtRatio

std::vector<double> InDet::InDetNNScoringTool::m_boundsTrtRatio

private

Definition at line 81 of file InDetNNScoringTool.h.

◆ m_caloClusterROIKey

SG::ReadHandleKey<ROIPhiRZContainer> InDet::InDetNNScoringTool::m_caloClusterROIKey {this, "EMROIPhiRZContainer", "", "Name of the calo cluster ROIs in Phi,R,Z parameterization"}

private

Definition at line 126 of file InDetNNScoringTool.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_DNN_inputValues

std::map<std::string, double> InDet::InDetNNScoringTool::m_DNN_inputValues

private

Definition at line 69 of file InDetNNScoringTool.h.

◆ m_etaWidthEm

float InDet::InDetNNScoringTool::m_etaWidthEm

private

Definition at line 124 of file InDetNNScoringTool.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extrapolator

ToolHandle<Trk::IExtrapolator> InDet::InDetNNScoringTool::m_extrapolator

private

Definition at line 91 of file InDetNNScoringTool.h.

◆ m_factorB_LayerHits

std::vector<double> InDet::InDetNNScoringTool::m_factorB_LayerHits

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_factorDblHoles

std::vector<double> InDet::InDetNNScoringTool::m_factorDblHoles

private

Definition at line 78 of file InDetNNScoringTool.h.

◆ m_factorGangedFakes

std::vector<double> InDet::InDetNNScoringTool::m_factorGangedFakes

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_factorHits

std::vector<double> InDet::InDetNNScoringTool::m_factorHits

private

Definition at line 78 of file InDetNNScoringTool.h.

◆ m_factorPixelHits

std::vector<double> InDet::InDetNNScoringTool::m_factorPixelHits

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_factorPixHoles

std::vector<double> InDet::InDetNNScoringTool::m_factorPixHoles

private

Definition at line 78 of file InDetNNScoringTool.h.

◆ m_factorPixLay

std::vector<double> InDet::InDetNNScoringTool::m_factorPixLay

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_factorSCT_Holes

std::vector<double> InDet::InDetNNScoringTool::m_factorSCT_Holes

private

Definition at line 78 of file InDetNNScoringTool.h.

◆ m_factorSigmaChi2

std::vector<double> InDet::InDetNNScoringTool::m_factorSigmaChi2

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_factorTrtFittedRatio

std::vector<double> InDet::InDetNNScoringTool::m_factorTrtFittedRatio

private

Definition at line 81 of file InDetNNScoringTool.h.

◆ m_factorTrtRatio

std::vector<double> InDet::InDetNNScoringTool::m_factorTrtRatio

private

Definition at line 81 of file InDetNNScoringTool.h.

◆ m_fieldCacheCondObjInputKey

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

private

Definition at line 94 of file InDetNNScoringTool.h.

◆ m_graph

std::unique_ptr<lwt::LightweightGraph> InDet::InDetNNScoringTool::m_graph

private

Definition at line 68 of file InDetNNScoringTool.h.

◆ m_maxB_LayerHits

int InDet::InDetNNScoringTool::m_maxB_LayerHits

private

Definition at line 77 of file InDetNNScoringTool.h.

◆ m_maxDblHoles

int InDet::InDetNNScoringTool::m_maxDblHoles

private

Definition at line 76 of file InDetNNScoringTool.h.

◆ m_maxDoubleHoles

int InDet::InDetNNScoringTool::m_maxDoubleHoles

private

Definition at line 113 of file InDetNNScoringTool.h.

◆ m_maxEta

double InDet::InDetNNScoringTool::m_maxEta

private

maximal Eta cut

Definition at line 109 of file InDetNNScoringTool.h.

◆ m_maxGangedFakes

int InDet::InDetNNScoringTool::m_maxGangedFakes

private

Definition at line 77 of file InDetNNScoringTool.h.

◆ m_maxHits

int InDet::InDetNNScoringTool::m_maxHits

private

Definition at line 76 of file InDetNNScoringTool.h.

◆ m_maxPixelHits

int InDet::InDetNNScoringTool::m_maxPixelHits

private

Definition at line 77 of file InDetNNScoringTool.h.

◆ m_maxPixelHoles

int InDet::InDetNNScoringTool::m_maxPixelHoles

private

maximal number of Pixel holes

Definition at line 115 of file InDetNNScoringTool.h.

◆ m_maxPixHoles

int InDet::InDetNNScoringTool::m_maxPixHoles

private

Definition at line 76 of file InDetNNScoringTool.h.

◆ m_maxPixLay

int InDet::InDetNNScoringTool::m_maxPixLay

private

Definition at line 77 of file InDetNNScoringTool.h.

◆ m_maxRPhiImp

double InDet::InDetNNScoringTool::m_maxRPhiImp

private

maximal RPhi impact parameter cut

Definition at line 110 of file InDetNNScoringTool.h.

◆ m_maxRPhiImpEM

double InDet::InDetNNScoringTool::m_maxRPhiImpEM

private

maximal RPhi impact parameter cut track that match EM clusters

Definition at line 120 of file InDetNNScoringTool.h.

◆ m_maxSCT_Holes

int InDet::InDetNNScoringTool::m_maxSCT_Holes

private

Definition at line 76 of file InDetNNScoringTool.h.

◆ m_maxSctHoles

int InDet::InDetNNScoringTool::m_maxSctHoles

private

maximal number of SCT holes

Definition at line 116 of file InDetNNScoringTool.h.

◆ m_maxSigmaChi2

int InDet::InDetNNScoringTool::m_maxSigmaChi2

private

Definition at line 76 of file InDetNNScoringTool.h.

◆ m_maxSiHoles

int InDet::InDetNNScoringTool::m_maxSiHoles

private

maximal number of holes (Pixel+SCT)

Definition at line 114 of file InDetNNScoringTool.h.

◆ m_maxTrtFittedRatio

int InDet::InDetNNScoringTool::m_maxTrtFittedRatio

private

Definition at line 76 of file InDetNNScoringTool.h.

◆ m_maxTrtRatio

int InDet::InDetNNScoringTool::m_maxTrtRatio

private

Definition at line 76 of file InDetNNScoringTool.h.

◆ m_maxZImp

double InDet::InDetNNScoringTool::m_maxZImp

private

maximal z impact parameter cut

Definition at line 111 of file InDetNNScoringTool.h.

◆ m_minPixel

int InDet::InDetNNScoringTool::m_minPixel

private

minimum number of pixel clusters

Definition at line 119 of file InDetNNScoringTool.h.

◆ m_minPt

double InDet::InDetNNScoringTool::m_minPt

private

cuts for selecting good tracks

minimal Pt cut

Definition at line 108 of file InDetNNScoringTool.h.

◆ m_minSiClusters

int InDet::InDetNNScoringTool::m_minSiClusters

private

minimal number of Si clusters

Definition at line 112 of file InDetNNScoringTool.h.

◆ m_minTRTonTrk

int InDet::InDetNNScoringTool::m_minTRTonTrk

private

minimum number of TRT hits

Definition at line 117 of file InDetNNScoringTool.h.

◆ m_minTRTprecision

double InDet::InDetNNScoringTool::m_minTRTprecision

private

minimum fraction of TRT precision hits

Definition at line 118 of file InDetNNScoringTool.h.

◆ m_nnCutConfig

std::string InDet::InDetNNScoringTool::m_nnCutConfig

private

Definition at line 101 of file InDetNNScoringTool.h.

◆ m_nnCutThreshold

double InDet::InDetNNScoringTool::m_nnCutThreshold

private

Definition at line 102 of file InDetNNScoringTool.h.

◆ m_phiWidthEm

float InDet::InDetNNScoringTool::m_phiWidthEm

private

Definition at line 123 of file InDetNNScoringTool.h.

◆ m_selectortool

ToolHandle<ITrtDriftCircleCutTool> InDet::InDetNNScoringTool::m_selectortool

private

Returns minimum number of expected TRT drift circles depending on eta.

Definition at line 84 of file InDetNNScoringTool.h.

◆ m_summaryTypeScore

std::vector<Trk::TrackScore> InDet::InDetNNScoringTool::m_summaryTypeScore

private

holds the scores assigned to each Trk::SummaryType from the track's Trk::TrackSummary

Definition at line 87 of file InDetNNScoringTool.h.

◆ m_useAmbigFcn

bool InDet::InDetNNScoringTool::m_useAmbigFcn

private

use the scoring tuned to Ambiguity processing or not

Definition at line 98 of file InDetNNScoringTool.h.

◆ m_useEmClusSeed

bool InDet::InDetNNScoringTool::m_useEmClusSeed

private

Definition at line 122 of file InDetNNScoringTool.h.

◆ m_usePixel

bool InDet::InDetNNScoringTool::m_usePixel

private

Definition at line 104 of file InDetNNScoringTool.h.

◆ m_useSCT

bool InDet::InDetNNScoringTool::m_useSCT

private

Definition at line 105 of file InDetNNScoringTool.h.

◆ m_useSigmaChi2

bool InDet::InDetNNScoringTool::m_useSigmaChi2

private

Definition at line 100 of file InDetNNScoringTool.h.

◆ m_useTRT_AmbigFcn

bool InDet::InDetNNScoringTool::m_useTRT_AmbigFcn

private

Definition at line 99 of file InDetNNScoringTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ InDetNNScoringTool()

◆ ~InDetNNScoringTool()

Member Function Documentation

◆ ambigScore()

◆ calcNnScore()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ isEmCaloCompatible()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ score()

◆ setupScoreModifiers()

◆ simpleScore()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_beamSpotKey

◆ m_boundsSigmaChi2

◆ m_boundsTrtFittedRatio

◆ m_boundsTrtRatio

◆ m_caloClusterROIKey

◆ m_detStore

◆ m_DNN_inputValues

◆ m_etaWidthEm

◆ m_evtStore

◆ m_extrapolator

◆ m_factorB_LayerHits

◆ m_factorDblHoles

◆ m_factorGangedFakes

◆ m_factorHits

◆ m_factorPixelHits

◆ m_factorPixHoles

◆ m_factorPixLay

◆ m_factorSCT_Holes

◆ m_factorSigmaChi2

◆ m_factorTrtFittedRatio

◆ m_factorTrtRatio

◆ m_fieldCacheCondObjInputKey

◆ m_graph

◆ m_maxB_LayerHits

◆ m_maxDblHoles

◆ m_maxDoubleHoles

◆ m_maxEta

◆ m_maxGangedFakes

◆ m_maxHits

◆ m_maxPixelHits

◆ m_maxPixelHoles

◆ m_maxPixHoles

◆ m_maxPixLay

◆ m_maxRPhiImp