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"
  
class EcalShapeBase : public CaloVShape
{
   public:

      typedef std::vector<double> DVec ;
  
      EcalShapeBase( bool   aSaveDerivative ) ;

      ~EcalShapeBase() override ;

      double operator() ( double aTime ) const override ;

      double         timeOfThr()  const ;
      double         timeOfMax()  const ;
      double timeToRise() const override ;

      virtual double threshold()             const = 0 ;
  
      double derivative ( double time ) const ; // appears to not be used anywhere

      enum { kReadoutTimeInterval = 25 , // in nsec
         kNBinsPerNSec        = 10 , // granularity of internal array
         k1NSecBins           = kReadoutTimeInterval*kNBinsPerNSec ,
         k1NSecBinsTotal      = 2*k1NSecBins ,
         kNBinsStored         = k1NSecBinsTotal*kNBinsPerNSec
      } ;

      static const double qNSecPerBin ;

   protected:

      unsigned int timeIndex( double aTime ) const ;

      void buildMe() ;

      virtual void   fillShape( DVec& aVec ) const = 0 ;

   private:

      unsigned int m_firstIndexOverThreshold ;
      double       m_firstTimeOverThreshold  ;
      unsigned int m_indexOfMax ;
      double       m_timeOfMax  ;
      DVec  m_shape ;
      DVec* m_derivPtr ;
};
  


#endif