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
virtual Trk::TrackScore	score (const Trk::Track &track, bool checkBasicSel) const override
	create a score based on how good the passed track is
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
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`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	setupScoreModifiers ()
bool	isEmCaloCompatible (const Trk::Track &track, const EventContext &ctx) const
	Check if the cluster is compatible with a EM cluster.
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

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.
std::vector< Trk::TrackScore >	m_summaryTypeScore
	holds the scores assigned to each Trk::SummaryType from the track's Trk::TrackSummary
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
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
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)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
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()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ 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()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

◆ 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 asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

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

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_beamSpotKey

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

private

Definition at line 95 of file InDetNNScoringTool.h.

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

◆ 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.

147{this, "EMROIPhiRZContainer", "", "Name of the calo cluster ROIs in Phi,R,Z parameterization"};

◆ 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.

144{this, "etaWidthEM", 0.05};

◆ 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.

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

◆ 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.

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

◆ 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.

126{this, "maxDoubleHoles", 2, "maximum number of SCT double holes"};

◆ m_maxEta

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

private

Definition at line 117 of file InDetNNScoringTool.h.

117{this, "maxEta", 2.7, "maximal Eta cut"};

◆ 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.

130{this, "maxPixelHoles", 5, "max number of Pixel holes"};

◆ 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.

119{this, "maxRPhiImp", 10., "maximal RPhi impact parameter cut"};

◆ 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.

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

◆ 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.

132{this, "maxSCTHoles", 5, "max number of SCT holes"};

◆ 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.

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

◆ 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.

121{this, "maxZImp", 250., "maximal z impact parameter cut"};

◆ m_minPixel

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

private

Definition at line 137 of file InDetNNScoringTool.h.

138{this, "minPixel", 0, "minimum number of pixel clusters"};

◆ 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.

116{this, "minPt", 500., "minimal Pt cut"};

◆ m_minSiClusters

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

private

Definition at line 123 of file InDetNNScoringTool.h.

124{this, "minSiClusters", 7, "minimal number of Si clusters"};

◆ m_minTRTonTrk

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

private

Definition at line 133 of file InDetNNScoringTool.h.

134{this, "minTRTonTrk", 9, "minimum number of TRT hits"};

◆ 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.

136{this, "minTRTPrecisionFraction", 0.5, "minimum fraction of TRT precision hits"};

◆ m_nnCutConfig

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

private

Definition at line 109 of file InDetNNScoringTool.h.

109{this, "nnCutConfig", ""};

◆ m_nnCutThreshold

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

private

Definition at line 110 of file InDetNNScoringTool.h.

110{this, "nnCutThreshold", -0.1};

◆ m_phiWidthEm

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

private

Definition at line 143 of file InDetNNScoringTool.h.

143{this, "phiWidthEM", 0.075};

◆ 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.

90{this, "DriftCircleCutTool", "InDet::InDetTrtDriftCircleCutTool"};

◆ 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.

105{this, "useAmbigFcn", true};

◆ m_useEmClusSeed

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

private

Definition at line 142 of file InDetNNScoringTool.h.

142{this, "doEmCaloSeed", true};

◆ m_usePixel

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

private

Definition at line 112 of file InDetNNScoringTool.h.

112{this, "usePixel", true};

◆ m_useSCT

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

private

Definition at line 113 of file InDetNNScoringTool.h.

113{this, "useSCT", true};

◆ m_useSigmaChi2

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

private

Definition at line 107 of file InDetNNScoringTool.h.

107{this, "useSigmaChi2", false};

◆ m_useTRT_AmbigFcn

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

private

Definition at line 106 of file InDetNNScoringTool.h.

106{this, "useTRT_AmbigFcn", false};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

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()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ extraDeps_update_handler()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ isEmCaloCompatible()

◆ msg()

◆ 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

◆ m_maxRPhiImp

◆ m_maxRPhiImpEM

◆ m_maxSCT_Holes

◆ m_maxSctHoles

◆ m_maxSigmaChi2

◆ m_maxSiHoles

◆ m_maxTrtFittedRatio

◆ m_maxTrtRatio

◆ m_maxZImp

◆ m_minPixel