HcalSiPM.h

Go to the documentation of this file.

// -*- C++ -*-
#ifndef HcalSimAlgos_HcalSiPM_h
#define HcalSimAlgos_HcalSiPM_h
 
#include <vector>
#include <algorithm>
#include <unordered_map>
 
#include "CalibCalorimetry/HcalAlgos/interface/HcalSiPMnonlinearity.h"
 
namespace CLHEP {
  class HepRandomEngine;
}
 
class HcalSiPM {
public:
  HcalSiPM(int nCells = 1, double tau = 15.);
 
  virtual ~HcalSiPM();
 
  void resetSiPM() { std::fill(theSiPM.begin(), theSiPM.end(), -999.); }
  virtual double hitCells(CLHEP::HepRandomEngine*, unsigned int pes, double tempDiff = 0., double photonTime = 0.);
 
  virtual double totalCharge() const { return totalCharge(theLastHitTime); }
  virtual double totalCharge(double time) const;
 
  int getNCells() const { return theCellCount; }
  double getTau() const { return theTau; }
  double getCrossTalk() const { return theCrossTalk; }
  double getTempDep() const { return theTempDep; }
 
  void setNCells(int nCells);
  void setTau(double tau);
  void setCrossTalk(double xtalk);  //  Borel-Tanner "lambda"
  void setTemperatureDependence(double tempDep);
  void setSaturationPars(const std::vector<float>& pars);
 
protected:
  typedef std::pair<unsigned int, std::vector<double> > cdfpair;
  typedef std::unordered_map<unsigned int, cdfpair> cdfmap;
 
  // void expRecover(double dt);
 
  double cellCharge(double deltaTime) const;
  unsigned int addCrossTalkCells(CLHEP::HepRandomEngine* engine, unsigned int in_pes);
 
  //numerical random generation from Borel-Tanner distribution
  double Borel(unsigned int n, double lambda, unsigned int k);
  const cdfpair& BorelCDF(unsigned int k);
 
  unsigned int theCellCount;
  std::vector<double> theSiPM;
  double theTau;
  double theTauInv;
  double theCrossTalk;
  double theTempDep;
  double theLastHitTime;
 
  HcalSiPMnonlinearity* nonlin;
 
  cdfmap borelcdfs;
};
 
#endif  //HcalSimAlgos_HcalSiPM_h