CbmStsParAsic.h

Go to the documentation of this file.

/* Copyright (C) 2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese [committer] */
 
#ifndef CBMSTSPARASIC_H
#define CBMSTSPARASIC_H 1
 
#include <Rtypes.h>  // for THashConsistencyHolder, ClassDefNV
 
#include <algorithm>
#include <array>
#include <set>
#include <string>  // for string
 
class TF1;
 
class CbmStsParAsic {
 
 public:
  CbmStsParAsic() {};
 
 
  CbmStsParAsic(uint16_t nChannels, uint16_t nAdc, double dynRange, double threshold, double timeResol, double deadTime,
                double noise, double znr);
 
 
  CbmStsParAsic(const CbmStsParAsic&);
 
 
  CbmStsParAsic& operator=(const CbmStsParAsic& other);
 
 
  ~CbmStsParAsic();
 
 
  double AdcToCharge(uint16_t adc) const { return fThreshold + fDynRange / double(fNofAdc) * (double(adc) - 0.5); }
 
 
  uint16_t DeactivateRandomChannels(double fraction);
 
 
  uint16_t ChargeToAdc(double charge) const
  {
    return (uint16_t) std::clamp(1 + (charge - fThreshold) * fNofAdc / fDynRange, 0.0, 31.0);
  };
 
 
  double GetDeadTime() const { return fDeadTime; }
 
 
  double GetDynRange() const { return fDynRange; }
 
 
  uint16_t GetNofAdc() const { return fNofAdc; }
 
 
  uint16_t GetNofChannels() const { return fNofChannels; }
 
 
  double GetNoise() const { return fNoise; }
 
 
  double GetNoiseRate() const;
 
 
  double GetRandomNoiseCharge() const;
 
 
  double GetThreshold() const { return fThreshold; }
 
 
  double GetTimeResol() const { return fTimeResolution; }
 
 
  double GetZeroNoiseRate() const { return fZeroNoiseRate; }
 
 
  void Init();
 
 
  Bool_t IsChannelActive(uint16_t channel) const { return fDeadChannels.find(channel) == fDeadChannels.end(); }
 
 
  void Set(uint16_t nChannels, uint16_t nAdc, double dynRange, double threshold, double timeResol, double deadTime,
           double noise, double zeroNoiseRate, std::set<uint16_t> deadChannels = {});
 
 
  void SetTimeOffset(double offset) { fTimeOffset = offset; }
 
 
  void SetWalkCoef(std::array<double, 31> par) { fWalkCoef = par; }
 
  double GetWalkCoef(uint32_t uIdx) const
  {
    if (uIdx < 31) return fWalkCoef[uIdx];
    return 0.;
  }
 
  std::string ToString() const;
 
 
 private:
  uint16_t fNofChannels            = 0;   
  uint16_t fNofAdc                 = 0;   
  double fDynRange                 = 0.;  
  double fThreshold                = 0.;  
  double fTimeResolution           = 0.;  
  double fDeadTime                 = 0.;  
  double fNoise                    = 0.;  
  double fZeroNoiseRate            = 0.;  
  double fTimeOffset               = 0.;  
  std::array<double, 31> fWalkCoef = {};  
  std::set<uint16_t> fDeadChannels{};     
 
  bool fIsInit = kFALSE;  
 
  TF1* fNoiseCharge = nullptr;  
 
  ClassDefNV(CbmStsParAsic, 4);
};
 
#endif /* CBMSTSPARASIC_H */