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() ;

      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;
      const edm::EventSetup *m_es;
};
  


#endif