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 bool	passBasicSelections (const Trk::Track &track) const override
	check track selections independent from TrackSummary More...

virtual Trk::TrackScore	score (const Trk::Track &track, bool checkBasicSel) 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 = -1

int	m_maxPixHoles = -1

int	m_maxSCT_Holes = -1

int	m_maxHits = -1

int	m_maxSigmaChi2 = -1

int	m_maxTrtRatio = -1

int	m_maxTrtFittedRatio = -1

int	m_maxB_LayerHits = -1

int	m_maxPixelHits = -1

int	m_maxPixLay = -1

int	m_maxGangedFakes = -1

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 {this, "DriftCircleCutTool", "InDet::InDetTrtDriftCircleCutTool"}
	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 {this, "Extrapolator", "Trk::Extrapolator"}

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

BooleanProperty	m_useAmbigFcn {this, "useAmbigFcn", true}
	use the scoring tuned to Ambiguity processing or not More...

BooleanProperty	m_useTRT_AmbigFcn {this, "useTRT_AmbigFcn", false}

BooleanProperty	m_useSigmaChi2 {this, "useSigmaChi2", false}

StringProperty	m_nnCutConfig {this, "nnCutConfig", ""}

DoubleProperty	m_nnCutThreshold {this, "nnCutThreshold", -0.1}

BooleanProperty	m_usePixel {this, "usePixel", true}

BooleanProperty	m_useSCT {this, "useSCT", true}

DoubleProperty	m_minPt {this, "minPt", 500., "minimal Pt cut"}
	cuts for selecting good tracks More...

DoubleProperty	m_maxEta {this, "maxEta", 2.7, "maximal Eta cut"}

DoubleProperty	m_maxRPhiImp {this, "maxRPhiImp", 10., "maximal RPhi impact parameter cut"}

DoubleProperty	m_maxZImp {this, "maxZImp", 250., "maximal z impact parameter cut"}

IntegerProperty	m_minSiClusters {this, "minSiClusters", 7, "minimal number of Si clusters"}

IntegerProperty	m_maxDoubleHoles {this, "maxDoubleHoles", 2, "maximum number of SCT double holes"}

IntegerProperty	m_maxSiHoles {this, "maxSiHoles", 5, "max number of Silicon (Pixel+SCT) holes"}

IntegerProperty	m_maxPixelHoles {this, "maxPixelHoles", 5, "max number of Pixel holes"}

IntegerProperty	m_maxSctHoles {this, "maxSCTHoles", 5, "max number of SCT holes"}

IntegerProperty	m_minTRTonTrk {this, "minTRTonTrk", 9, "minimum number of TRT hits"}

DoubleProperty	m_minTRTprecision {this, "minTRTPrecisionFraction", 0.5, "minimum fraction of TRT precision hits"}

IntegerProperty	m_minPixel {this, "minPixel", 0, "minimum number of pixel clusters"}

DoubleProperty	m_maxRPhiImpEM {this, "maxRPhiImpEM", 50., "maximal RPhi impact parameter cut track that match EM clusters"}

BooleanProperty	m_useEmClusSeed {this, "doEmCaloSeed", true}

FloatProperty	m_phiWidthEm {this, "phiWidthEM", 0.075}

FloatProperty	m_etaWidthEm {this, "etaWidthEM", 0.05}

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_summaryTypeScore(Trk::numberOfTrackSummaryTypes)
 {
   declareInterface<Trk::ITrackScoringTool>(this);
  
   //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 342 of file InDetNNScoringTool.cxx.

 {
   //
   // --- start with bonus for high pt tracks
   //
   // double prob = 1.;
   double pt = std::abs(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: "
        << std::abs(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 815 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(std::abs(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 58 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 37 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 787 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.

◆ passBasicSelections()

bool InDet::InDetNNScoringTool::passBasicSelections ( const Trk::Track & track ) const

overridevirtual

check track selections independent from TrackSummary

Implements Trk::ITrackScoringTool.

Definition at line 109 of file InDetNNScoringTool.cxx.

 {
   //
   // --- 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 false;
   }
   MagField::AtlasFieldCache fieldCache;
   fieldCondObj->getInitializedCache (fieldCache);
  
   if (fieldCache.solenoidOn()){
     if (std::abs(input->pT()) < m_minPt) {
       ATH_MSG_DEBUG ("Track pt < "<<m_minPt<<", reject it");
       return false;
     }
   }
   if (std::abs(input->eta()) > m_maxEta) {
     ATH_MSG_DEBUG ("Track eta > "<<m_maxEta<<", reject it");
     return false;
   }
  
   // 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) );
  
   if (!parm || parm->surfaceType()!=Trk::SurfaceType::Perigee) {
      ATH_MSG_WARNING( "Extrapolation of perigee failed, this should never happen" );
      return false;
   }
  
   ATH_MSG_VERBOSE ("extrapolated perigee: "<<*parm);
   if (std::abs(parm->parameters()[Trk::z0]) > m_maxZImp) {
     ATH_MSG_DEBUG ("Track Z impact > "<<m_maxZImp<<", reject it");
     return false;
   }
  
   double maxD0 = m_maxRPhiImp;
   if(m_useEmClusSeed) {
      if (isEmCaloCompatible( track, ctx ) ) {
         maxD0 = m_maxRPhiImpEM;
      }
   }
   if (std::abs(parm->parameters()[Trk::d0]) > maxD0) {
     ATH_MSG_DEBUG ("Track Rphi impact > "<<maxD0<<", reject it");
     return false;
   }
  
   return true;
 }

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ score()

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

overridevirtual

create a score based on how good the passed track is

Implements Trk::ITrackScoringTool.

Definition at line 176 of file InDetNNScoringTool.cxx.

 {
    if(checkBasicSel && !passBasicSelections(track)){
       ATH_MSG_VERBOSE ("Track fail basic selections");
       return Trk::TrackScore(0);
    }
    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 561 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 193 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);
     }
   }
  
   // check if nn cut is enabled
   if (m_useAmbigFcn && m_nnCutThreshold > 0.0) {
     if (track.fitQuality()) { // check if track has been fitted
       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();
       const EventContext& ctx = Gaudi::Hive::currentContext();
       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 (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 95 of file InDetNNScoringTool.h.

◆ m_boundsSigmaChi2

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

private

Definition at line 85 of file InDetNNScoringTool.h.

◆ m_boundsTrtFittedRatio

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

private

Definition at line 86 of file InDetNNScoringTool.h.

◆ m_boundsTrtRatio

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

private

Definition at line 86 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 146 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 72 of file InDetNNScoringTool.h.

◆ m_etaWidthEm

FloatProperty InDet::InDetNNScoringTool::m_etaWidthEm {this, "etaWidthEM", 0.05}

private

Definition at line 144 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 {this, "Extrapolator", "Trk::Extrapolator"}

private

Definition at line 97 of file InDetNNScoringTool.h.

◆ m_factorB_LayerHits

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

private

Definition at line 84 of file InDetNNScoringTool.h.

◆ m_factorDblHoles

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

private

Definition at line 83 of file InDetNNScoringTool.h.

◆ m_factorGangedFakes

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

private

Definition at line 84 of file InDetNNScoringTool.h.

◆ m_factorHits

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

private

Definition at line 83 of file InDetNNScoringTool.h.

◆ m_factorPixelHits

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

private

Definition at line 84 of file InDetNNScoringTool.h.

◆ m_factorPixHoles

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

private

Definition at line 83 of file InDetNNScoringTool.h.

◆ m_factorPixLay

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

private

Definition at line 84 of file InDetNNScoringTool.h.

◆ m_factorSCT_Holes

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

private

Definition at line 83 of file InDetNNScoringTool.h.

◆ m_factorSigmaChi2

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

private

Definition at line 84 of file InDetNNScoringTool.h.

◆ m_factorTrtFittedRatio

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

private

Definition at line 86 of file InDetNNScoringTool.h.

◆ m_factorTrtRatio

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

private

Definition at line 86 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 101 of file InDetNNScoringTool.h.

◆ m_graph

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

private

Definition at line 71 of file InDetNNScoringTool.h.

◆ m_maxB_LayerHits

int InDet::InDetNNScoringTool::m_maxB_LayerHits = -1

private

Definition at line 81 of file InDetNNScoringTool.h.

◆ m_maxDblHoles

int InDet::InDetNNScoringTool::m_maxDblHoles = -1

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_maxDoubleHoles

IntegerProperty InDet::InDetNNScoringTool::m_maxDoubleHoles {this, "maxDoubleHoles", 2, "maximum number of SCT double holes"}

private

Definition at line 125 of file InDetNNScoringTool.h.

◆ m_maxEta

DoubleProperty InDet::InDetNNScoringTool::m_maxEta {this, "maxEta", 2.7, "maximal Eta cut"}

private

Definition at line 117 of file InDetNNScoringTool.h.

◆ m_maxGangedFakes

int InDet::InDetNNScoringTool::m_maxGangedFakes = -1

private

Definition at line 82 of file InDetNNScoringTool.h.

◆ m_maxHits

int InDet::InDetNNScoringTool::m_maxHits = -1

private

Definition at line 80 of file InDetNNScoringTool.h.

◆ m_maxPixelHits

int InDet::InDetNNScoringTool::m_maxPixelHits = -1

private

Definition at line 81 of file InDetNNScoringTool.h.

◆ m_maxPixelHoles

IntegerProperty InDet::InDetNNScoringTool::m_maxPixelHoles {this, "maxPixelHoles", 5, "max number of Pixel holes"}

private

Definition at line 129 of file InDetNNScoringTool.h.

◆ m_maxPixHoles

int InDet::InDetNNScoringTool::m_maxPixHoles = -1

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_maxPixLay

int InDet::InDetNNScoringTool::m_maxPixLay = -1

private

Definition at line 82 of file InDetNNScoringTool.h.

◆ m_maxRPhiImp

DoubleProperty InDet::InDetNNScoringTool::m_maxRPhiImp {this, "maxRPhiImp", 10., "maximal RPhi impact parameter cut"}

private

Definition at line 118 of file InDetNNScoringTool.h.

◆ m_maxRPhiImpEM

DoubleProperty InDet::InDetNNScoringTool::m_maxRPhiImpEM {this, "maxRPhiImpEM", 50., "maximal RPhi impact parameter cut track that match EM clusters"}

private

Definition at line 140 of file InDetNNScoringTool.h.

◆ m_maxSCT_Holes

int InDet::InDetNNScoringTool::m_maxSCT_Holes = -1

private

Definition at line 79 of file InDetNNScoringTool.h.

◆ m_maxSctHoles

IntegerProperty InDet::InDetNNScoringTool::m_maxSctHoles {this, "maxSCTHoles", 5, "max number of SCT holes"}

private

Definition at line 131 of file InDetNNScoringTool.h.

◆ m_maxSigmaChi2

int InDet::InDetNNScoringTool::m_maxSigmaChi2 = -1

private

Definition at line 80 of file InDetNNScoringTool.h.

◆ m_maxSiHoles

IntegerProperty InDet::InDetNNScoringTool::m_maxSiHoles {this, "maxSiHoles", 5, "max number of Silicon (Pixel+SCT) holes"}

private

Definition at line 127 of file InDetNNScoringTool.h.

◆ m_maxTrtFittedRatio

int InDet::InDetNNScoringTool::m_maxTrtFittedRatio = -1

private

Definition at line 81 of file InDetNNScoringTool.h.

◆ m_maxTrtRatio

int InDet::InDetNNScoringTool::m_maxTrtRatio = -1

private

Definition at line 80 of file InDetNNScoringTool.h.

◆ m_maxZImp

DoubleProperty InDet::InDetNNScoringTool::m_maxZImp {this, "maxZImp", 250., "maximal z impact parameter cut"}

private

Definition at line 120 of file InDetNNScoringTool.h.

◆ m_minPixel

IntegerProperty InDet::InDetNNScoringTool::m_minPixel {this, "minPixel", 0, "minimum number of pixel clusters"}

private

Definition at line 137 of file InDetNNScoringTool.h.

◆ m_minPt

DoubleProperty InDet::InDetNNScoringTool::m_minPt {this, "minPt", 500., "minimal Pt cut"}

private

cuts for selecting good tracks

Definition at line 116 of file InDetNNScoringTool.h.

◆ m_minSiClusters

IntegerProperty InDet::InDetNNScoringTool::m_minSiClusters {this, "minSiClusters", 7, "minimal number of Si clusters"}

private

Definition at line 123 of file InDetNNScoringTool.h.

◆ m_minTRTonTrk

IntegerProperty InDet::InDetNNScoringTool::m_minTRTonTrk {this, "minTRTonTrk", 9, "minimum number of TRT hits"}

private

Definition at line 133 of file InDetNNScoringTool.h.

◆ m_minTRTprecision

DoubleProperty InDet::InDetNNScoringTool::m_minTRTprecision {this, "minTRTPrecisionFraction", 0.5, "minimum fraction of TRT precision hits"}

private

Definition at line 135 of file InDetNNScoringTool.h.

◆ m_nnCutConfig

StringProperty InDet::InDetNNScoringTool::m_nnCutConfig {this, "nnCutConfig", ""}

private

Definition at line 109 of file InDetNNScoringTool.h.

◆ m_nnCutThreshold

DoubleProperty InDet::InDetNNScoringTool::m_nnCutThreshold {this, "nnCutThreshold", -0.1}

private

Definition at line 110 of file InDetNNScoringTool.h.

◆ m_phiWidthEm

FloatProperty InDet::InDetNNScoringTool::m_phiWidthEm {this, "phiWidthEM", 0.075}

private

Definition at line 143 of file InDetNNScoringTool.h.

◆ m_selectortool

ToolHandle<ITrtDriftCircleCutTool> InDet::InDetNNScoringTool::m_selectortool {this, "DriftCircleCutTool", "InDet::InDetTrtDriftCircleCutTool"}

private

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

Definition at line 89 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 93 of file InDetNNScoringTool.h.

◆ m_useAmbigFcn

BooleanProperty InDet::InDetNNScoringTool::m_useAmbigFcn {this, "useAmbigFcn", true}

private

use the scoring tuned to Ambiguity processing or not

Definition at line 105 of file InDetNNScoringTool.h.

◆ m_useEmClusSeed

BooleanProperty InDet::InDetNNScoringTool::m_useEmClusSeed {this, "doEmCaloSeed", true}

private

Definition at line 142 of file InDetNNScoringTool.h.

◆ m_usePixel

BooleanProperty InDet::InDetNNScoringTool::m_usePixel {this, "usePixel", true}

private

Definition at line 112 of file InDetNNScoringTool.h.

◆ m_useSCT

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

private

Definition at line 113 of file InDetNNScoringTool.h.

◆ m_useSigmaChi2

BooleanProperty InDet::InDetNNScoringTool::m_useSigmaChi2 {this, "useSigmaChi2", false}

private

Definition at line 107 of file InDetNNScoringTool.h.

◆ m_useTRT_AmbigFcn

BooleanProperty InDet::InDetNNScoringTool::m_useTRT_AmbigFcn {this, "useTRT_AmbigFcn", false}

private

Definition at line 106 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()

◆ passBasicSelections()

◆ 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