RawTowerDigitizer.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RawTowerDigitizer.h

#ifndef G4CALO_RAWTOWERDIGITIZER_H
#define G4CALO_RAWTOWERDIGITIZER_H

#include <fun4all/SubsysReco.h>

#include <phparameter/PHParameters.h>

#include <string>

class PHCompositeNode;
class RawTowerContainer;
class RawTowerGeomContainer;
class RawTowerDeadMap;

class RawTower;

#include <gsl/gsl_rng.h>

class RawTowerDigitizer : public SubsysReco
{
 public:
  RawTowerDigitizer(const std::string &name = "RawTowerDigitizer");
  ~RawTowerDigitizer() override;

  int InitRun(PHCompositeNode *topNode) override;
  int process_event(PHCompositeNode *topNode) override;
  void Detector(const std::string &d) { m_Detector = d; _tower_params.set_name(d);}
  void TowerType(const int type) { m_TowerType = type; }
  void set_seed(const unsigned int iseed);
  unsigned int get_seed() const { return m_Seed; }
  enum enu_digi_algorithm
  {
    kNo_digitization = 0,
    kNo_digitalization = 0,

    kSimple_photon_digitization = 1,
    kSimple_photon_digitalization = 1,

    kSiPM_photon_digitization = 2
  };

  enu_digi_algorithm
  get_digi_algorithm() const
  {
    return m_DigiAlgorithm;
  }

  void
  set_digi_algorithm(enu_digi_algorithm digiAlgorithm)
  {
    m_DigiAlgorithm = digiAlgorithm;
  }

  double
  get_pedstal_central_ADC() const
  {
    return m_PedstalCentralADC;
  }

  void
  set_pedstal_central_ADC(const double pedstalCentralAdc)
  {
    m_PedstalCentralADC = pedstalCentralAdc;
  }

  double
  get_pedstal_width_ADC() const
  {
    return m_PedstalWidthADC;
  }

  void
  set_pedstal_width_ADC(const double pedstalWidthAdc)
  {
    m_PedstalWidthADC = pedstalWidthAdc;
  }

  double
  get_photonelec_ADC() const
  {
    return m_PhotonElecADC;
  }

  void
  set_photonelec_ADC(const double photonelecAdc)
  {
    m_PhotonElecADC = photonelecAdc;
  }

  double
  get_photonelec_yield_visible_GeV() const
  {
    return m_PhotonElecYieldVisibleGeV;
  }

  void
  set_photonelec_yield_visible_GeV(const double photonelecYieldVisibleGeV)
  {
    m_PhotonElecYieldVisibleGeV = photonelecYieldVisibleGeV;
  }

  double
  get_zero_suppression_ADC() const
  {
    return m_ZeroSuppressionADC;
  }

  void
  set_zero_suppression_ADC(const double zeroSuppressionAdc)
  {
    m_ZeroSuppressionADC = zeroSuppressionAdc;
  }

  void
  set_variable_zero_suppression(const bool value)
  {
    m_ZeroSuppressionFile = value;
  }

  void
  set_variable_pedestal(const bool value)
  {
    m_pedestalFile = value;
  }

  PHParameters &
  GetParameters()
  {
    return _tower_params;
  }

  std::string
  get_raw_tower_node_prefix() const
  {
    return m_RawTowerNodePrefix;
  }

  void
  set_raw_tower_node_prefix(const std::string &rawTowerNodePrefix)
  {
    m_RawTowerNodePrefix = rawTowerNodePrefix;
  }

  std::string
  get_sim_tower_node_prefix() const
  {
    return m_SimTowerNodePrefix;
  }

  void
  set_sim_tower_node_prefix(const std::string &simTowerNodePrefix)
  {
    m_SimTowerNodePrefix = simTowerNodePrefix;
  }

  // ! SiPM effective pixel per tower, only used with kSiPM_photon_digitalization
  void set_sipm_effective_pixel(const unsigned int &d) { m_SiPMEffectivePixel = d; }

  // ! SiPM effective pixel per tower, only used with kSiPM_photon_digitalization
  unsigned int get_sipm_effective_pixel() { return m_SiPMEffectivePixel; }

 private:
  void CreateNodes(PHCompositeNode *topNode);

  RawTower *simple_photon_digitization(RawTower *sim_tower);

  RawTower *sipm_photon_digitization(RawTower *sim_tower);

  enu_digi_algorithm m_DigiAlgorithm;

  RawTowerContainer *m_SimTowers;
  RawTowerContainer *m_RawTowers;
  RawTowerGeomContainer *m_RawTowerGeom;
  RawTowerDeadMap *m_DeadMap;

  std::string m_Detector;

  std::string m_SimTowerNodePrefix;
  std::string m_RawTowerNodePrefix;

  double m_PhotonElecYieldVisibleGeV;

  double m_PhotonElecADC;

  double m_PedstalCentralADC;

  double m_PedstalWidthADC;

  bool m_pedestalFile;

  double m_ZeroSuppressionADC;

  bool m_ZeroSuppressionFile;

  int m_TowerType;

  unsigned int m_Seed;

  // ! SiPM effective pixel per tower, only used with kSiPM_photon_digitalization
  unsigned int m_SiPMEffectivePixel;

  PHParameters _tower_params;

  gsl_rng *m_RandomGenerator;
};

#endif /* G4CALO_RAWTOWERDIGITIZER_H */