EBHitResponse.h

Go to the documentation of this file.

#ifndef EcalSimAlgos_EBHitResponse_h
#define EcalSimAlgos_EBHitResponse_h

#include "CalibFormats/CaloObjects/interface/CaloTSamples.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalHitResponse.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstantsMC.h"
#include "DataFormats/EcalDigi/interface/EcalConstants.h"

class APDSimParameters;

namespace CLHEP {
  class HepRandomEngine;
}

template <class constset>
class EBHitResponseImpl : public EcalHitResponse {
public:
  typedef CaloTSamples<float, constset::sampleSize> EBSamples;

  typedef std::vector<double> VecD;

  static constexpr size_t kNOffsets = constset::kNOffsets;

  static constexpr double kSamplePeriod = constset::Samp_Period;

  EBHitResponseImpl(const CaloVSimParameterMap* parameterMap,
                    const CaloVShape* shape,
                    bool apdOnly,
                    const APDSimParameters* apdPars = nullptr,
                    const CaloVShape* apdShape = nullptr);

  ~EBHitResponseImpl() override;

  void initialize(CLHEP::HepRandomEngine*);

  virtual bool keepBlank() const { return false; }

  void setIntercal(const EcalIntercalibConstantsMC* ical);

  void add(const PCaloHit& hit, CLHEP::HepRandomEngine*) override;

  void initializeHits() override;

  void finalizeHits() override;

  void run(MixCollection<PCaloHit>& hits, CLHEP::HepRandomEngine*) override;

  unsigned int samplesSize() const override;

  EcalSamples* operator[](unsigned int i) override;

  const EcalSamples* operator[](unsigned int i) const override;

protected:
  unsigned int samplesSizeAll() const override;

  EcalSamples* vSamAll(unsigned int i) override;

  const EcalSamples* vSamAll(unsigned int i) const override;

  EcalSamples* vSam(unsigned int i) override;

  void putAPDSignal(const DetId& detId, double npe, double time);

  void putAnalogSignal(const PCaloHit& inputHit, CLHEP::HepRandomEngine*) override;

private:
  const VecD& offsets() const { return m_timeOffVec; }

  const double nonlFunc(double enr) const {
    return (pelo > enr ? pext : (pehi > enr ? nonlFunc1(enr) : pfac * atan(log10(enr - pehi + 0.00001)) + poff));
  }

  const double nonlFunc1(double energy) const {
    const double enr(log10(energy));
    const double enr2(enr * enr);
    const double enr3(enr2 * enr);
    return (pcub * enr3 + pqua * enr2 + plin * enr + pcon);
  }

  const APDSimParameters* apdParameters() const;
  const CaloVShape* apdShape() const;

  double apdSignalAmplitude(const PCaloHit& hit, CLHEP::HepRandomEngine*) const;

  void findIntercalibConstant(const DetId& detId, double& icalconst) const;

  const bool m_apdOnly;
  const APDSimParameters* m_apdPars;
  const CaloVShape* m_apdShape;
  const EcalIntercalibConstantsMC* m_intercal;

  std::vector<double> m_timeOffVec;

  std::vector<double> m_apdNpeVec;
  std::vector<double> m_apdTimeVec;

  const double pcub, pqua, plin, pcon, pelo, pehi, pasy, pext, poff, pfac;

  std::vector<EBSamples> m_vSam;

  bool m_isInitialized;
};

typedef EBHitResponseImpl<ecalPh1> EBHitResponse;
typedef EBHitResponseImpl<ecalPh2> EBHitResponse_Ph2;
#include "EBHitResponse.icc"
#endif