EcalShapeBase.h

Go to the documentation of this file.

#ifndef EcalSimAlgos_EcalShapeBase_h
#define EcalSimAlgos_EcalShapeBase_h
 
#include <vector>
//#include<stdexcept>
#include "SimCalorimetry/CaloSimAlgos/interface/CaloVShape.h"
 
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "CondFormats/EcalObjects/interface/EcalSimPulseShape.h"
#include "CondFormats/DataRecord/interface/EcalSimPulseShapeRcd.h"
 
#include <iostream>
#include <fstream>
class EcalShapeBase : public CaloVShape {
public:
  typedef std::vector<double> DVec;
 
  EcalShapeBase(bool);
 
  ~EcalShapeBase() override;
 
  double operator()(double aTime) const override;
 
  double timeOfThr() const;
  double timeOfMax() const;
  double timeToRise() const override;
 
  double threshold() const;
 
  double derivative(double time) const;  // appears to not be used anywhere
 
  void m_shape_print(const char* fileName);
  void setEventSetup(const edm::EventSetup& evtSetup);
 
protected:
  unsigned int timeIndex(double aTime) const;
 
  void buildMe(const edm::EventSetup* = nullptr);
 
  virtual void fillShape(float& time_interval,
                         double& m_thresh,
                         EcalShapeBase::DVec& aVec,
                         const edm::EventSetup* es) const = 0;
  bool m_useDBShape;
 
private:
  unsigned int m_firstIndexOverThreshold;
  double m_firstTimeOverThreshold;
  unsigned int m_indexOfMax;
  double m_timeOfMax;
  double m_thresh;
  unsigned int m_kNBinsPerNSec;
  DVec m_shape;
  DVec m_deriv;
 
  unsigned int m_arraySize;
  unsigned int m_denseArraySize;
  double m_qNSecPerBin;
};
 
#endif