'Use EF tracking cuts in SiTrackMaker' (TrigFastTrackFinder-00-00-98)

Trigger/TrigAlgorithms/TrigFastTrackFinder/TrigFastTrackFinder/TrigFastTrackFinder.h

+/*
+  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+////////////////////////////////////////////////////////////////////////////////
+//
+// filename: TrigFastTrackFinder.h
+// 
+// Description: a part of L2+EF HLT ID tracking
+// 
+// date: 16/04/2013
+// 
+// -------------------------------
+// ATLAS Collaboration
+////////////////////////////////////////////////////////////////////////////////
+
+#ifndef __TRIG_FAST_TRACK_FINDER_H__
+#define __TRIG_FAST_TRACK_FINDER_H__
+
+#include<string>
+#include<vector>
+#include<map>
+//#include<algorithm>
+
+#include "GaudiKernel/ToolHandle.h"
+#include "TrigInterfaces/FexAlgo.h"
+
+#include "TrkEventPrimitives/ParticleHypothesis.h"
+
+#include "TrigInDetEvent/TrigInDetTrackCollection.h"
+#include "TrigInDetPattRecoEvent/TrigInDetPattRecoEvent/TrigInDetRoad.h"
+#include "TrigInDetPattRecoTools/TrigCombinatorialSettings.h"
+
+class IFTK_DataProviderTool;
+class ITrigL2LayerNumberTool;
+class ITrigL2LayerSetPredictorTool;
+class ITrigSpacePointConversionTool;
+class ITrigInDetRoadMakerTool;
+class ITrigL2SpacePointTruthTool;
+class ITrigInDetTrackFitter;
+class IRegSelSvc;
+class TrigRoiDescriptor;
+class Identifier;
+namespace InDet { 
+  class SiSpacePointsSeed;
+  class ISiTrackMaker; 
+}
+
+namespace MagField {	
+  class IMagFieldSvc;
+}
+
+namespace Trk {
+  class ITrackSummaryTool;
+}
+
+class TrigL2LayerSetLUT;
+class TrigSpacePointStorage;
+class TrigInDetDoublet;
+class TrigInDetTriplet;
+class TrigInDetTripletCluster;
+class TrigCombinatorialTrackFinding;
+class IBeamCondSvc;
+//class EventID;
+class PixelID;
+class SCT_ID;
+class AtlasDetectorID;
+
+class TrigFastTrackFinder : public HLT::FexAlgo {
+
+ public:
+  
+  TrigFastTrackFinder(const std::string& name, ISvcLocator* pSvcLocator);
+  ~TrigFastTrackFinder();
+  HLT::ErrorCode hltInitialize();
+  HLT::ErrorCode hltFinalize();
+  HLT::ErrorCode hltBeginRun();
+
+  HLT::ErrorCode hltExecute(const HLT::TriggerElement* inputTE,
+			    HLT::TriggerElement* outputTE);
+
+  bool isInRoad(const TrigSiSpacePointBase*, const TrigInDetRoad*);
+
+  void convertToTrigInDetTrack(const TrackCollection& offlineTracks, TrigInDetTrackCollection& trigInDetTracks);
+  double trackQuality(const Trk::Track* Tr);
+  void filterSharedTracks(std::vector<std::tuple<bool, double, Trk::Track*>>& QT);
+  void createOfflineSeeds(const std::vector<std::shared_ptr<TrigInDetTriplet>>& input, std::vector<InDet::SiSpacePointsSeed>& output);
+
+protected: 
+
+  void updateClusterMap(long int, const Trk::Track*, std::map<Identifier, std::vector<long int> >&);
+  void extractClusterIds(const Trk::SpacePoint*, std::vector<Identifier>&);
+  bool usedByAnyTrack(const std::vector<Identifier>&, std::map<Identifier, std::vector<long int> >&);
+
+  int findBarCodeInData(int, const std::vector<TrigSiSpacePointBase>&);
+  void showBarCodeInData(int, const std::vector<TrigSiSpacePointBase>&);
+  int findBarCodeInTriplets(int, const std::vector<std::shared_ptr<TrigInDetTriplet>>&);
+  int findBarCodeInDoublets(int, const std::vector<TrigInDetDoublet*>&);
+  void assignTripletBarCodes(const std::vector<std::shared_ptr<TrigInDetTriplet>>&, std::vector<int>&);
+  void assignTripletBarCodes(const std::vector<TrigInDetTriplet*>&, std::vector<int>&);
+  
+ private:
+
+  // AlgTools and Services
+
+  ToolHandle<ITrigL2LayerNumberTool> m_numberingTool;
+  ToolHandle<ITrigL2LayerSetPredictorTool> m_predictorTool;
+  ToolHandle<ITrigSpacePointConversionTool> m_spacePointTool;
+  ToolHandle<ITrigInDetRoadMakerTool> m_roadMakerTool;
+  ToolHandle<ITrigL2SpacePointTruthTool> m_TrigL2SpacePointTruthTool;
+  ToolHandle<InDet::ISiTrackMaker> m_trackMaker;   // Track maker 
+  ToolHandle<ITrigInDetTrackFitter> m_trigInDetTrackFitter;   // Track maker 
+  ToolHandle< Trk::ITrackSummaryTool > m_trackSummaryTool;
+
+  // ToolHandle<IFTK_DataProviderTool> m_ftkReader;
+ 
+  ServiceHandle<IRegSelSvc>     m_regionSelector;      //!< region selector service
+
+  ServiceHandle<MagField::IMagFieldSvc> m_MagFieldSvc;
+
+  double m_shift_x, m_shift_y;
+
+  // Control flags
+
+  bool m_useNewSeeding;//to activate roadless seed generator
+
+  bool m_ftkMode;
+  bool m_useBeamSpot; 
+  bool m_vertexSeededMode;
+
+  HLT::ErrorCode makeTripletClusters(const TrigSpacePointStorage& spacePointStorage, 
+                                     TrigCombinatorialTrackFinding& combinatorial, 
+                                     std::vector<TrigInDetTripletCluster*>& tripletClusters,  
+                                     const TrigL2LayerSetLUT* pLUT,
+                                     std::map<int,int>& nGoodDoublets,
+                                     const std::vector<TrigSiSpacePointBase>& convertedSpacePoints);
+
+  // Cuts and settings
+  TrigCombinatorialSettings m_tcs;
+
+  int  m_minHits;
+
+  int                     m_nSeeds;       //!< Number seeds 
+  int                     m_nTracksNew;      //!< Number found tracks 
+  int                     m_nfreeCut;     // Min number free clusters 
+
+  int                     m_nUsedLayers;
+
+  bool m_retrieveBarCodes;
+
+  float m_tripletMinPtFrac;
+  float m_pTmin;
+
+  bool m_checkSeedRedundancy;
+
+  // Roads
+
+  std::vector<TrigInDetRoad> m_roads;
+
+  // Names of IDCs with input data 
+  //
+  //
+  // Vertex (from beamspot)
+  Amg::Vector3D m_vertex;
+
+ 
+  // Reconstructed tracks 
+
+  TrigInDetTrackCollection* m_recoTracks;
+
+  IBeamCondSvc* m_iBeamCondSvc;
+  
+  // Data members for monitoring
+
+  int m_ntracks;
+  int m_nPixSPsInRoI;  // Total number of (filtered) pixel SPs in the RoI
+  int m_nSCTSPsInRoI;  // Total number of (filtered) SCT SPs in the RoI
+  int m_currentStage;  // The last stage reached during the processing of a given RoI
+
+  int m_roi_nSPs;
+
+  double m_roiPhi, m_roiEta;  
+  double m_roiPhiWidth, m_roiEtaWidth;  
+  double m_timePattReco;
+
+  int m_nZvertices; 
+  std::vector<float> m_zVertices; 
+
+  std::vector<double> m_a0beam;
+  std::vector<double> m_trkdPhi0, m_trkdEta;
+  //std::vector<double> m_sp_x, m_sp_y, m_sp_z, m_sp_r;//Spacepoint coordinates
+
+  std::vector<double> m_pixResPhiBarrel;
+  std::vector<double> m_pixResEtaBarrel;
+  std::vector<double> m_pixPullPhiBarrel;
+  std::vector<double> m_pixPullEtaBarrel;
+  std::vector<double> m_sctResBarrel;
+  std::vector<double> m_sctPullBarrel;
+  std::vector<double> m_pixResPhiEC;
+  std::vector<double> m_pixResEtaEC;
+  std::vector<double> m_pixPullPhiEC;
+  std::vector<double> m_pixPullEtaEC;
+  std::vector<double> m_sctResEC;
+  std::vector<double> m_sctPullEC;
+
+  // Monitoring member functions 
+
+  static inline double monPt(const TrigInDetTrack *t){return t->param()->pT(); }
+  static inline double monA0(const TrigInDetTrack *t){return t->param()->a0(); }
+  static inline double monZ0(const TrigInDetTrack *t){return t->param()->z0(); }
+  static inline double monPhi0(const TrigInDetTrack *t){return t->param()->phi0();}
+  static inline double monEta(const TrigInDetTrack *t){return t->param()->eta(); }
+  static inline double monNHit_Si(const TrigInDetTrack *t){return (t->NPixelSpacePoints()+t->NSCT_SpacePoints());}
+  static inline double monChi2(const TrigInDetTrack *t){return (t->chi2()>9e7)?-9.9:t->chi2();}
+
+  void fill_a0();
+
+  void calculateRecoEfficiency(const std::vector<TrigSiSpacePointBase>&,
+			       const std::map<int,int>&,
+			       const std::map<int,int>&,
+			       const std::map<int,int>&);
+
+
+  //Setup functions
+  void clearMembers();
+  void getBeamSpot();
+  void getMagField();
+  HLT::ErrorCode getRoI(const HLT::TriggerElement* inputTE, const IRoiDescriptor*& roi);
+
+  // Timers 
+
+  TrigTimer* m_SpacePointConversionTimer;
+  TrigTimer* m_PatternRecoTimer; 
+  TrigTimer* m_SpacePointSortingTimer;
+  TrigTimer* m_DoubletFindingTimer; 
+  TrigTimer* m_TripletFindingTimer; 
+  TrigTimer* m_TripletClusterTimer; 
+  TrigTimer* m_TripletMakingTimer; 
+  TrigTimer* m_CombTrackingTimer; 
+  TrigTimer* m_TrackFitterTimer; 
+
+  // Other member functions
+
+  StatusCode storeSpacePoints(const std::vector<TrigSiSpacePointBase>& convertedSpacePoints, TrigSpacePointStorage& spacePointStorage);
+
+  // Internal bookkeeping
+
+  std::string m_instanceName, m_attachedFeatureName1, m_attachedFeatureName2, 
+    m_outputCollectionSuffix;
+
+  unsigned int m_l1Id;
+  unsigned int m_countTotalEvents;
+  unsigned int m_countTotalRoI;
+  unsigned int m_countRoIwithEnoughHits;
+  unsigned int m_countRoIwithTracks;
+
+  //efficiency calculations
+  std::vector<int> m_vSignalBarCodes;
+
+  int m_nSignalPresent;
+  int m_nSignalDetected;
+  int m_nSignalTracked;
+  int m_nGoodDoublets;
+  int m_nSignalClones;
+  const PixelID* m_pixelId;
+  const SCT_ID* m_sctId;
+  const AtlasDetectorID* m_idHelper;
+  StoreGateSvc*  m_detectorStore;
+
+  //Merge clone triplet clusters?
+  bool m_doCloneMerge;
+  bool m_doCloneRemove;
+  int m_numSeedsToTry;
+  int m_minSignalSPs;
+
+	bool m_roiForIDWarning=false;
+
+  Trk::ParticleHypothesis m_particleHypothesis = Trk::pion;//particle hypothesis to attach to each track - usually pion, can be set to other values
+
+  std::map<Identifier, std::vector<long int> > m_siClusterMap;
+
+};
+
+#endif

Trigger/TrigAlgorithms/TrigFastTrackFinder/cmt/requirements

+package TrigFastTrackFinder
+
+# =====================================================================
+
+public
+
+# General
+use AtlasPolicy			AtlasPolicy-*
+use GaudiInterface		GaudiInterface-*     External
+
+# Trigger
+use TrigInterfaces              TrigInterfaces-*	 Trigger/TrigSteer
+use TrigInDetEvent		TrigInDetEvent-*	 Trigger/TrigEvent
+use TrigSteeringEvent           TrigSteeringEvent-*      Trigger/TrigEvent
+use TrigInDetPattRecoEvent	TrigInDetPattRecoEvent-*   Trigger/TrigEvent
+use TrigInDetPattRecoTools      TrigInDetPattRecoTools-*   Trigger/TrigTools
+
+use TrkEventPrimitives   TrkEventPrimitives-*   Tracking/TrkEvent
+
+# =====================================================================
+
+apply_pattern declare_python_modules files="*.py"
+
+apply_pattern component_library
+
+library TrigFastTrackFinder  *.cxx components/*.cxx 
+
+private
+
+# Tracking
+
+use     AtlasTBB                        AtlasTBB-*                      External
+use 	TrkRIO_OnTrack			TrkRIO_OnTrack-*		Tracking/TrkEvent 
+use     InDetRIO_OnTrack                InDetRIO_OnTrack-*              InnerDetector/InDetRecEvent 
+use	TrkParameters			TrkParameters-*			Tracking/TrkEvent 
+use     EventPrimitives                 EventPrimitives-*               Event
+use TrkToolInterfaces           TrkToolInterfaces-*           Tracking/TrkTools
+
+use	InDetRecToolInterfaces		InDetRecToolInterfaces-* 	InnerDetector/InDetRecTools 
+use	IRegionSelector			IRegionSelector-*		DetectorDescription
+use     InDetIdentifier                 InDetIdentifier-*               InnerDetector/InDetDetDescr
+
+use TrigInDetToolInterfaces     TrigInDetToolInterfaces-*  Trigger/TrigTools
+use TrkTrack		        TrkTrack-* 		   Tracking/TrkEvent
+use TrigTimeAlgs		TrigTimeAlgs-*		   Trigger/TrigTools
+use StoreGate			StoreGate-*		   Control
+use InDetBeamSpotService        InDetBeamSpotService-*     InnerDetector/InDetConditions
+
+use AthenaBaseComps              AthenaBaseComps-*	   Control
+use InDetPrepRawData     	 InDetPrepRawData-*	   InnerDetector/InDetRecEvent 
+use SiSpacePointsSeed            SiSpacePointsSeed-*       InnerDetector/InDetRecEvent 
+use TrigInDetRecoTools           TrigInDetRecoTools-*      Trigger/TrigTools
+use MagFieldInterfaces           MagFieldInterfaces-*      MagneticField 

Trigger/TrigAlgorithms/TrigFastTrackFinder/src/components/TrigFastTrackFinder_entries.cxx

+#include "GaudiKernel/DeclareFactoryEntries.h"
+
+#include "TrigFastTrackFinder/TrigFastTrackFinder.h"
+
+DECLARE_ALGORITHM_FACTORY( TrigFastTrackFinder)
+DECLARE_FACTORY_ENTRIES( TrigFastTrackFinder )
+{
+  DECLARE_ALGORITHM( TrigFastTrackFinder )
+}

Trigger/TrigAlgorithms/TrigFastTrackFinder/src/components/TrigFastTrackFinder_load.cxx

+#include "GaudiKernel/LoadFactoryEntries.h"
+
+LOAD_FACTORY_ENTRIES( TrigFastTrackFinder )