SiSmearedDigitizationTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// SiSmearedDigitizationTool.h
//   Header file for class SiSmearedDigitizationTool
// (c) ATLAS Detector software
// Top algorithm class for Pixel fast smeared digitization
 
#ifndef FASTSIDIGITIZATION_SISMEAREDDIGITIZATIONTOOL_H
#define FASTSIDIGITIZATION_SISMEAREDDIGITIZATIONTOOL_H
 
#include "PileUpTools/PileUpToolBase.h"
 
#include "GaudiKernel/ToolHandle.h"
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/AlgTool.h"
#include "xAODEventInfo/EventInfo.h"
#include "xAODEventInfo/EventAuxInfo.h"
#include "InDetSimData/InDetSimDataCollection.h"
#include <string>
 
#include "HitManagement/TimedHitCollection.h"
#include "InDetSimEvent/SiHit.h"
#include "InDetSimEvent/SiHitCollection.h"
 
#include "InDetPrepRawData/SCT_ClusterContainer.h"  // typedef
#include "InDetPrepRawData/PixelClusterContainer.h"
 
#include "TrkTruthData/PRD_MultiTruthCollection.h"
 
#include "AthenaKernel/IAthRNGSvc.h"
#include "GaudiKernel/ITHistSvc.h"
#include "EventPrimitives/EventPrimitives.h"
 
#include "InDetReadoutGeometry/SiDetectorElementCollection.h"
#include "StoreGate/ReadCondHandleKey.h"
 
#include <tuple>
typedef std::tuple< Amg::Vector2D, InDet::SiWidth, Amg::MatrixX> ClusterInfo;
 
class PixelID;
class SCT_ID;
class TTree;
class TFile;
 
namespace InDetDD{class SiDetectorElement;}
namespace CLHEP {class HepRandomEngine;}
 
namespace InDet {
  class ClusterMakerTool;
  class PixelCluster;
  class SCT_Cluster;
  class SiCluster;
}
 
class SiSmearedDigitizationTool : virtual public PileUpToolBase
{
 
public:
 
  SiSmearedDigitizationTool(
                               const std::string& type,
                               const std::string& name,
                               const IInterface* parent
                               );
 
 
  StatusCode initialize();
  StatusCode prepareEvent(const EventContext& ctx, unsigned int);
  StatusCode processBunchXing( int bunchXing,
                               SubEventIterator bSubEvents,
                               SubEventIterator eSubEvents );
  StatusCode processAllSubEvents(const EventContext& ctx);
  StatusCode mergeEvent(const EventContext& ctx);
 
  typedef std::multimap<IdentifierHash, InDet::PixelCluster*> Pixel_detElement_RIO_map;
  typedef std::multimap<IdentifierHash, InDet::SCT_Cluster*> SCT_detElement_RIO_map;
  StatusCode mergeClusters(Pixel_detElement_RIO_map * cluster_map);
  StatusCode mergeClusters(SCT_detElement_RIO_map * cluster_map);
 
  StatusCode digitize(const EventContext& ctx,
                      TimedHitCollection<SiHit>& thpcsi);
  StatusCode createAndStoreRIOs(const EventContext& ctx);
  StatusCode retrieveTruth();
  StatusCode finalize();
 
  template<typename CLUSTER>
    double calculateDistance(CLUSTER * clusterA, CLUSTER * clusterB);
 
  template<typename CLUSTER>
    double calculateSigma(CLUSTER * clusterA, CLUSTER * clusterB);
 
  template<typename CLUSTER>
    ClusterInfo calculateNewCluster(CLUSTER * clusterA, CLUSTER * clusterB);
 
  template<typename CLUSTER>
    StatusCode FillTruthMap(PRD_MultiTruthCollection*, CLUSTER*, const TimedHitPtr<SiHit>&);
 
 private:
 
  ServiceHandle<IAthRNGSvc> m_rndmSvc{this, "RndmSvc", "AthRNGSvc", ""};  
 
  SG::ReadCondHandleKey<InDetDD::SiDetectorElementCollection> m_pixelDetEleCollKey{this, "PixelDetEleCollKey", "PixelDetectorElementCollection", "Key of SiDetectorElementCollection for Pixel"};
  SG::ReadCondHandleKey<InDetDD::SiDetectorElementCollection> m_SCTDetEleCollKey{this, "SCTDetEleCollKey", "SCT_DetectorElementCollection", "Key of SiDetectorElementCollection for SCT"};
 
  const PixelID* m_pixel_ID;                             
  const SCT_ID* m_sct_ID;                             
 
  std::string                m_randomEngineName;         
 
  float m_pitch_X;
  float m_pitch_Y;
 
  bool m_merge;
  double m_nSigma;
 
  bool m_useDiscSurface;
 
  InDet::PixelClusterContainer*  m_pixelClusterContainer;               
 
  InDet::SCT_ClusterContainer*  m_sctClusterContainer;               
 
  ServiceHandle<PileUpMergeSvc> m_mergeSvc;      
  int                       m_HardScatterSplittingMode; 
  bool                      m_HardScatterSplittingSkipper;
 
  PRD_MultiTruthCollection* m_pixelPrdTruth{};
  std::string               m_prdTruthNamePixel;
 
  PRD_MultiTruthCollection* m_SCTPrdTruth{};
  std::string               m_prdTruthNameSCT;
 
  SiHitCollection* m_simHitColl{};
  std::string      m_inputObjectName;     
 
  std::vector<SiHitCollection*> m_siHitCollList;
 
  Pixel_detElement_RIO_map* m_pixelClusterMap{};
 
  SCT_detElement_RIO_map* m_sctClusterMap{};
 
  bool m_SmearPixel;
 
  bool m_emulateAtlas;
 
  std::string                           m_pixel_SiClustersName;
  std::string                           m_Sct_SiClustersName;
 
  bool       m_checkSmear;
 
  ITHistSvc* m_thistSvc;
  TFile*              m_outputFile;  
  TTree*              m_currentTree; 
 
  double           m_x_pixel;
  double           m_y_pixel;
  double           m_x_exit_pixel;
  double           m_y_exit_pixel;
  double           m_z_exit_pixel;
  double           m_x_entry_pixel;
  double           m_y_entry_pixel;
  double           m_z_entry_pixel;
  double           m_x_pixel_global;
  double           m_y_pixel_global;
  double           m_z_pixel_global;
 
  double           m_x_SCT;
  double           m_x_exit_SCT;
  double           m_y_exit_SCT;
  double           m_z_exit_SCT;
  double           m_x_entry_SCT;
  double           m_y_entry_SCT;
  double           m_z_entry_SCT;
  double           m_x_SCT_global;
  double           m_y_SCT_global;
  double           m_z_SCT_global;
 
  double           m_x_pixel_smeared;
  double           m_y_pixel_smeared;
  double           m_x_SCT_smeared;
 
  double           m_Err_x_pixel;
  double           m_Err_y_pixel;
  double           m_Err_x_SCT;
  double           m_Err_y_SCT;
 
 
  SiSmearedDigitizationTool();
  SiSmearedDigitizationTool(const SiSmearedDigitizationTool&);
 
  SiSmearedDigitizationTool& operator=(const SiSmearedDigitizationTool&);
 
 
};
 
#endif // SISMEAREDDDIGITIZATION_SISMEAREDDIGITIZATIONTOOL_H