CaloIsCMOSDigitizerAlgo.h

Go to the documentation of this file.

/*
 * CaloIsCMOSDigitizerAlgo.h
 *
 *  Created on: 11 Nov 2022
 *      Authors: Nicola Mori
 *               Zheng Quan
 */
 
#ifndef HERD_CALOISCMOSDIGITIZERALGO_H_
#define HERD_CALOISCMOSDIGITIZERALGO_H_
 
// HerdSoftware headers
#include "herd/analysis/dataobjects/CaloGeoParams.h"
#include "herd/analysis/dataobjects/CaloHits.h"
#include "herd/analysis/dataobjects/CaloIsCMOSChannelInfo.h"
 
// EventAnalysis headers
#include "algorithm/Algorithm.h"
 
// C++ header
#include <random>
 
namespace Herd {
 
// clang-format off
// clang-format on
class CaloIsCMOSDigitizerAlgo : public EA::Algorithm {
public:
  CaloIsCMOSDigitizerAlgo(const std::string &name);
 
  bool Initialize() override;
 
  bool Process() override;
 
private:
  EA::StorePtr m_evStore, m_globStore;
  std::mt19937 m_randomEngine;
  std::unique_ptr<std::normal_distribution<double>> m_distPedHG, m_distPedLG;
  CaloHits m_hgfHits, m_lgfHits;
 
  const double m_oneMIP = 0.0302; // GeV
 
  // Parameters
  unsigned int m_seed = defaultValue<decltype(m_seed)>;
  float m_stddevPedHG = 20.;
  float m_stddevPedLG = 60.;
 
  // Utility functions
  // Taken from original code provided by Z. Quan.
  double GetGrayscaleHighGain(double eDep, unsigned int volumeID, CaloIsCMOSChannelInfo &chInfo);
  double GetGrayscaleLowGain(double eDep, unsigned int volumeID, CaloIsCMOSChannelInfo &chInfo);
  double GetNonLinHighGain(double grayscale);
  double GetNonLinLowGain(double grayscale);
  double GetRMSHighGain(double grayScale);
  double GetRMSLowGain(double grayScale);
};
 
} // namespace Herd
 
#endif /* HERD_CALOISCMOSDIGITIZERALGO_H_ */