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