InDetNNScoringTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/

// Markus Elsing
 
#ifndef INDETNNSCORINGTOOL_H
#define INDETNNSCORINGTOOL_H
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ToolHandle.h"
#include "GaudiKernel/ServiceHandle.h"
#include "StoreGate/ReadHandleKey.h"
#include "TrkCaloClusterROI/ROIPhiRZContainer.h"
#include "TrkEventPrimitives/TrackScore.h"
#include "TrkToolInterfaces/ITrackScoringTool.h"
#include "lwtnn/lightweight_network_config.hh"
#include "lwtnn/LightweightNeuralNetwork.hh"
#include "lwtnn/LightweightGraph.hh"
#include "lwtnn/parse_json.hh"
#include "PathResolver/PathResolver.h"
#include "TrkParameters/TrackParameters.h"
#include "AthenaKernel/CLASS_DEF.h"
#include "AthenaKernel/IOVSvcDefs.h"
// MagField cache
#include "MagFieldConditions/AtlasFieldCacheCondObj.h"
#include "MagFieldElements/AtlasFieldCache.h"
#include <fstream>
#include <vector>
#include <string>
#include "BeamSpotConditionsData/BeamSpotData.h"
 
namespace Trk {
  class IExtrapolator;
  class Track;
  class TrackSummary;
 
}
class EventContext;
 
namespace InDet {
class ITrtDriftCircleCutTool;

class InDetNNScoringTool : virtual public Trk::ITrackScoringTool,  
                             public AthAlgTool
{
 
 public:
  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;

  virtual Trk::TrackScore score( const Trk::Track& track, bool checkBasicSel ) const override;
  
  virtual Trk::TrackScore simpleScore( const Trk::Track& track, const Trk::TrackSummary& trackSum ) const override;
  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;
  
 
 private:
  void setupScoreModifiers();
  
    // DNN tools
    std::unique_ptr<lwt::LightweightGraph> m_graph;
    std::map<std::string, double> m_DNN_inputValues;
  
  bool isEmCaloCompatible(const Trk::Track& track, const EventContext& ctx) const;
  
  
  //these are used for ScoreModifiers 
  int m_maxDblHoles = -1, m_maxPixHoles = -1, m_maxSCT_Holes = -1,
    m_maxHits = -1, m_maxSigmaChi2 = -1, m_maxTrtRatio = -1,
    m_maxTrtFittedRatio = -1, m_maxB_LayerHits = -1, m_maxPixelHits = -1,
    m_maxPixLay = -1,  m_maxGangedFakes = -1;
  std::vector<double> m_factorDblHoles, m_factorPixHoles, m_factorSCT_Holes,  m_factorHits,
    m_factorSigmaChi2, m_factorB_LayerHits, m_factorPixelHits, m_factorPixLay, m_factorGangedFakes;
  std::vector<double> m_boundsSigmaChi2,
    m_boundsTrtRatio, m_factorTrtRatio, m_boundsTrtFittedRatio, m_factorTrtFittedRatio;

  ToolHandle<ITrtDriftCircleCutTool> m_selectortool
    {this, "DriftCircleCutTool", "InDet::InDetTrtDriftCircleCutTool"};
  
  std::vector<Trk::TrackScore>           m_summaryTypeScore;
  
  SG::ReadCondHandleKey<InDet::BeamSpotData> m_beamSpotKey { this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" };
  
  ToolHandle<Trk::IExtrapolator> m_extrapolator
    {this, "Extrapolator", "Trk::Extrapolator"};
 
  // Read handle for conditions object to get the field cache
  SG::ReadCondHandleKey<AtlasFieldCacheCondObj> m_fieldCacheCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"};
 
  
  BooleanProperty m_useAmbigFcn{this, "useAmbigFcn", true};
  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"};
  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"};
};
 
 
} // namespace InDet
 
 
#endif