EcalShapeBase Class Reference

#include <EcalShapeBase.h>

Inheritance diagram for EcalShapeBase:

Public Types
typedef std::vector< double >	DVec

Public Member Functions
double	derivative (double time) const
	EcalShapeBase (bool)
void	m_shape_print (const char *fileName)
double	operator() (double aTime) const override
void	setEventSetup (const edm::EventSetup &evtSetup)
double	threshold () const
double	timeOfMax () const
double	timeOfThr () const
double	timeToRise () const override
	~EcalShapeBase () override
Public Member Functions inherited from CaloVShape
	CaloVShape ()
virtual	~CaloVShape ()

Protected Member Functions
void	buildMe ()
virtual void	fillShape (float &time_interval, double &m_thresh, EcalShapeBase::DVec &aVec, const edm::EventSetup *es) const =0
unsigned int	timeIndex (double aTime) const

Protected Attributes
bool	m_useDBShape

Private Attributes
unsigned int	m_arraySize
unsigned int	m_denseArraySize
DVec	m_deriv
const edm::EventSetup *	m_es
unsigned int	m_firstIndexOverThreshold
double	m_firstTimeOverThreshold
unsigned int	m_indexOfMax
unsigned int	m_kNBinsPerNSec
double	m_qNSecPerBin
DVec	m_shape
double	m_thresh
double	m_timeOfMax

Detailed Description

Definition at line 24 of file EcalShapeBase.h.

Member Typedef Documentation

typedef std::vector<double> EcalShapeBase::DVec

Definition at line 28 of file EcalShapeBase.h.

Constructor & Destructor Documentation

EcalShapeBase::EcalShapeBase

(

bool

useDBShape

)

Definition at line 13 of file EcalShapeBase.cc.

                                            :
   m_useDBShape              (useDBShape),
   m_firstIndexOverThreshold ( 0   ) ,
   m_firstTimeOverThreshold  ( 0.0 ) ,
   m_indexOfMax              ( 0   ) ,
   m_timeOfMax               ( 0.0 ) , 
   m_thresh                  ( 0.0 ) , 
   m_es                      ( nullptr )
{
}

EcalShapeBase::~EcalShapeBase ( )

override

Definition at line 9 of file EcalShapeBase.cc.

{
}

Member Function Documentation

void EcalShapeBase::buildMe ( )

protected

Definition at line 54 of file EcalShapeBase.cc.

References a, b, EnergyCorrector::c, delta, MillePedeFileConverter_cfg::e, connectstrParser::f1, connectstrParser::f2, connectstrParser::f3, fillShape(), mps_fire::i, createfilelist::int, m_arraySize, m_denseArraySize, m_deriv, m_es, m_firstIndexOverThreshold, m_firstTimeOverThreshold, m_indexOfMax, m_kNBinsPerNSec, m_qNSecPerBin, m_shape, m_thresh, m_timeOfMax, scrapingFilter_cfi::thresh, threshold(), relativeConstraints::value, and x.

Referenced by APDShape::APDShape(), EBShape::EBShape(), EEShape::EEShape(), and setEventSetup().

{
   DVec shapeArray;

   float time_interval = 0;
   fillShape(time_interval, m_thresh, shapeArray, m_es) ; // pure virtual function, implementation may vary for EB/EE/APD ...
   m_arraySize = shapeArray.size();  // original data
  
   m_denseArraySize    = 10*m_arraySize; // dense array with interpolation between data 
   m_kNBinsPerNSec     = (unsigned int) (10/time_interval); // used to be an unsigned int = 10 in  < CMSSW10X, should work for time intervals ~0.1, 0.2, 0.5, 1
   m_qNSecPerBin = time_interval/10.; 

   for(unsigned int i = 0; i < m_denseArraySize; ++i) { m_deriv.push_back(0.0); m_shape.push_back(0.0); }

   const double maxel ( *max_element( shapeArray.begin(), shapeArray.end() ) ) ;

   const double maxelt ( 1.e-5 < maxel ? maxel : 1 ) ;

   for( unsigned int i ( 0 ) ; i != shapeArray.size(); ++i )
   {
      shapeArray[i] = shapeArray[i]/maxelt ;
   } 

   const double thresh ( threshold()/maxelt ) ;

   

   const double delta ( m_qNSecPerBin/2. ) ; 

   for( unsigned int denseIndex ( 0 ) ; denseIndex != m_denseArraySize ; ++denseIndex )
   {
      const double xb ( ( denseIndex + 0.5 )*m_qNSecPerBin ) ; 

      const unsigned int ibin ( denseIndex/10 ) ; 

      double value = 0.0 ;
      double deriv = 0.0 ;

      if( 0                 ==     ibin ||
      shapeArray.size() == 1 + ibin    ) // cannot do quadratic interpolation at ends
      {
     value = shapeArray[ibin]; 
     deriv = 0.0 ;
      }
      else 
      {
     const double x  ( xb - ( ibin + 0.5 )*time_interval ) ;
     const double f1 ( shapeArray[ ibin - 1 ] ) ;
     const double f2 ( shapeArray[ ibin     ] ) ;
     const double f3 ( shapeArray[ ibin + 1 ] ) ;
     const double a  ( f2 ) ;
     const double b  ( ( f3 - f1 )/(2.*time_interval) ) ;
     const double c  ( (( f1 + f3 )/2. - f2 )/(time_interval*time_interval) );
     value = a + b*x + c*x*x;
     deriv = ( b + 2*c*x )/delta ;
      }

      m_shape[ denseIndex ] = value;
      m_deriv[ denseIndex ] = deriv;

      if( 0      <  denseIndex                         &&
      thresh <  value                     &&
      0      == m_firstIndexOverThreshold     )
      {
     m_firstIndexOverThreshold = denseIndex - 1 ;
     m_firstTimeOverThreshold  = m_firstIndexOverThreshold*m_qNSecPerBin ;
      }

      if( m_shape[ m_indexOfMax ] < value )
      {
     m_indexOfMax = denseIndex ;
      }

   }
   m_timeOfMax = m_indexOfMax*m_qNSecPerBin ;

}

double EcalShapeBase::derivative

(

double

time

)

const

Definition at line 158 of file EcalShapeBase.cc.

References m_denseArraySize, m_deriv, and timeIndex().

{
   const unsigned int index ( timeIndex( aTime ) ) ;
   return ( m_denseArraySize == index         ? 0 : m_deriv[ index ] ) ;
}

virtual void EcalShapeBase::fillShape ( float &  time_interval, double &  m_thresh, EcalShapeBase::DVec &  aVec, const edm::EventSetup *  es  ) const

protectedpure virtual

Implemented in APDShape, EBShape, and EEShape.

Referenced by buildMe().

void EcalShapeBase::m_shape_print

(

const char *

fileName

)

Definition at line 165 of file EcalShapeBase.cc.

References mps_fire::i, and m_shape.

{std::ofstream fs; fs.open(fileName); for (auto i:m_shape) fs << "vec.push_back(" << i << ");\n"; fs.close();}

double EcalShapeBase::operator() ( double  aTime ) const

overridevirtual

Implements CaloVShape.

Definition at line 149 of file EcalShapeBase.cc.

References m_denseArraySize, m_shape, and timeIndex().

{
   // return pulse amplitude for request time in ns

   const unsigned int index ( timeIndex( aTime ) ) ;
   return ( m_denseArraySize == index ? 0 : m_shape[ index ] ) ;
}

void EcalShapeBase::setEventSetup

(

const edm::EventSetup &

evtSetup

)

Definition at line 24 of file EcalShapeBase.cc.

References buildMe(), and m_es.

Referenced by EcalDigiProducer::accumulate(), EcalTPGParamBuilder::analyze(), EcalDigiProducer::beginRun(), EcalMixingModuleValidation::dqmBeginRun(), EcalUncalibRecHitWorkerWeights::set(), and EcalUncalibRecHitWorkerGlobal::set().

{ m_es = &evtSetup;  buildMe();} 

double EcalShapeBase::threshold

(

)

const

Definition at line 47 of file EcalShapeBase.cc.

References m_thresh.

Referenced by buildMe().

{
   return m_thresh; 
}

unsigned int EcalShapeBase::timeIndex ( double  aTime ) const

protected

Definition at line 133 of file EcalShapeBase.cc.

References LogDebug, m_denseArraySize, m_firstIndexOverThreshold, and m_kNBinsPerNSec.

Referenced by derivative(), and operator()().

{
   const int index ( m_firstIndexOverThreshold +
             (unsigned int) ( aTime*m_kNBinsPerNSec + 0.5 ) ) ;

   const bool bad ( (int) m_firstIndexOverThreshold >  index || 
            (int) m_denseArraySize              <= index    ) ;

   if(          (int) m_denseArraySize              <= index    )
   {
      LogDebug("EcalShapeBase") << " ECAL MGPA shape requested for out of range time " << aTime ;
   }
   return ( bad ? m_denseArraySize : (unsigned int) index ) ;
}

double EcalShapeBase::timeOfMax

(

)

const

Definition at line 34 of file EcalShapeBase.cc.

References m_timeOfMax.

Referenced by EcalTPGParamBuilder::computeWeights(), and timeToRise().

{
   return m_timeOfMax              ; 
}

double EcalShapeBase::timeOfThr

(

)

const

Definition at line 28 of file EcalShapeBase.cc.

References m_firstTimeOverThreshold.

Referenced by EcalTPGParamBuilder::computeWeights(), and timeToRise().

{
   return m_firstTimeOverThreshold ;
 }

double EcalShapeBase::timeToRise ( ) const

overridevirtual

Implements CaloVShape.

Definition at line 40 of file EcalShapeBase.cc.

References timeOfMax(), and timeOfThr().

Referenced by EcalDccWeightBuilder::computeWeights(), and EcalUncalibRecHitRecChi2Algo< C >::EcalUncalibRecHitRecChi2Algo().

{
   return timeOfMax() - timeOfThr() ;
}

Member Data Documentation

unsigned int EcalShapeBase::m_arraySize

private

Definition at line 67 of file EcalShapeBase.h.

Referenced by buildMe().

unsigned int EcalShapeBase::m_denseArraySize

private

Definition at line 68 of file EcalShapeBase.h.

Referenced by buildMe(), derivative(), operator()(), and timeIndex().

DVec EcalShapeBase::m_deriv

private

Definition at line 65 of file EcalShapeBase.h.

Referenced by buildMe(), and derivative().

const edm::EventSetup* EcalShapeBase::m_es

private

Definition at line 70 of file EcalShapeBase.h.

Referenced by buildMe(), and setEventSetup().

unsigned int EcalShapeBase::m_firstIndexOverThreshold

private

Definition at line 58 of file EcalShapeBase.h.

Referenced by buildMe(), and timeIndex().

double EcalShapeBase::m_firstTimeOverThreshold

private

Definition at line 59 of file EcalShapeBase.h.

Referenced by buildMe(), and timeOfThr().

unsigned int EcalShapeBase::m_indexOfMax

private

Definition at line 60 of file EcalShapeBase.h.

Referenced by buildMe().

unsigned int EcalShapeBase::m_kNBinsPerNSec

private

Definition at line 63 of file EcalShapeBase.h.

Referenced by buildMe(), and timeIndex().

double EcalShapeBase::m_qNSecPerBin

private

Definition at line 69 of file EcalShapeBase.h.

Referenced by buildMe().

DVec EcalShapeBase::m_shape

private

Definition at line 64 of file EcalShapeBase.h.

Referenced by buildMe(), m_shape_print(), and operator()().

double EcalShapeBase::m_thresh

private

Definition at line 62 of file EcalShapeBase.h.

Referenced by APDShape::APDShape(), buildMe(), EBShape::EBShape(), EEShape::EEShape(), and threshold().

double EcalShapeBase::m_timeOfMax

private

Definition at line 61 of file EcalShapeBase.h.

Referenced by buildMe(), and timeOfMax().

bool EcalShapeBase::m_useDBShape

protected

Definition at line 54 of file EcalShapeBase.h.

Referenced by APDShape::fillShape(), EEShape::fillShape(), and EBShape::fillShape().