FCS_StepInfoSD.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_FASTCALOSIM_FCS_STEPINFOSD_H
#define ISF_FASTCALOSIM_FCS_STEPINFOSD_H
 
#include <map>
#include <vector>
 
#include "CLHEP/Units/SystemOfUnits.h"
#include "CaloDetDescr/CaloDetDescrManager.h"
#include "CaloIdentifier/CaloCell_ID.h"  // For CaloCell_ID::CaloSample
#include "CxxUtils/CachedPointer.h"
#include "G4VSensitiveDetector.hh"
#include "ISF_FastCaloSimEvent/FCS_StepInfoCollection.h"
#include "LArG4Code/ILArCalculatorSvc.h"
#include "LArG4Code/LArG4Identifier.h"
#include "LArSimEvent/LArHit.h"
#include "TileG4Interfaces/ITileCalculator.h"
 
// Forward declarations
class G4Step;
class G4TouchableHistory;
 
class LArEM_ID;
class LArFCAL_ID;
class LArHEC_ID;
class TileID;
 
class LArHitContainer;
 
class StoreGateSvc;
 
namespace FCS_Param {
 
struct Config {
  int verboseLevel = 0;
  bool shift_lar_subhit = true;
  bool shorten_lar_step = false;
  double substpsize = 0.2 * CLHEP::mm;  // size of splitting into substeps
                                        // before calling the calculators.
 
  // Merging properties
  std::vector<double> m_maxRadiusLateral{
      24, 5.};  

  std::vector<double> m_maxRadiusLongitudinal{
      24, 5.};  

 
  double m_maxTime = 25.;
  double m_maxTimeLAr = 25.;
  double m_maxTimeHEC = 25.;
  double m_maxTimeFCAL = 25.;
  double m_maxTimeTile = 25.;
 
  ILArCalculatorSvc* m_LArCalculator = nullptr;
  ITileCalculator* m_TileCalculator = nullptr;
};
 
}  // namespace FCS_Param

class FCS_StepInfoSD : public G4VSensitiveDetector {
 public:
  FCS_StepInfoSD(G4String a_name, const FCS_Param::Config& config);

  virtual G4bool ProcessHits(G4Step* a_step, G4TouchableHistory*) override;

  void EndOfAthenaEvent(
      ISF_FCS_Parametrization::FCS_StepInfoCollection* hitContnainer);

  void setupHelpers(const LArEM_ID* EM, const LArFCAL_ID* FCAL,
                    const LArHEC_ID* HEC, const TileID* tile) {
    m_larEmID = EM;
    m_larFcalID = FCAL;
    m_larHecID = HEC;
    m_tileID = tile;
  }
 
 protected:
  void getCaloDDManager();
  void update_map(const CLHEP::Hep3Vector& l_vec,
                  const Identifier& l_identifier, double l_energy,
                  double l_time, bool l_valid, int l_detector);
  FCS_Param::Config m_config;
  const LArEM_ID* m_larEmID{nullptr};
  const LArFCAL_ID* m_larFcalID{nullptr};
  const LArHEC_ID* m_larHecID{nullptr};
  const TileID* m_tileID{nullptr};
  CxxUtils::CachedPointer<const CaloDetDescrManager> m_calo_dd_man;
  std::map<Identifier, std::vector<ISF_FCS_Parametrization::FCS_StepInfo*>*>
      m_hit_map;
 
 private:
  double getMaxTime(const CaloCell_ID::CaloSample& layer) const;
};
 
#endif  // ISF_FASTCALOSIM_FCS_STEPINFOSD_H