#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 (const edm::EventSetup *=nullptr)

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

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

◆ DVec

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

Definition at line 26 of file EcalShapeBase.h.

Constructor & Destructor Documentation

◆ EcalShapeBase()

EcalShapeBase::EcalShapeBase ( bool useDBShape )

Definition at line 11 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) {}

◆ ~EcalShapeBase()

EcalShapeBase::~EcalShapeBase ( )

override

Definition at line 9 of file EcalShapeBase.cc.

9 {}

Member Function Documentation

◆ buildMe()

void EcalShapeBase::buildMe ( const edm::EventSetup * evtSetup = nullptr )

protected

Definition at line 29 of file EcalShapeBase.cc.

References a, b, c, dumpMFGeometry_cfg::delta, MillePedeFileConverter_cfg::e, DeadROC_duringRun::f1, DeadROC_duringRun::f2, fillShape(), mps_fire::i, createfilelist::int, m_arraySize, m_denseArraySize, m_deriv, m_firstIndexOverThreshold, m_firstTimeOverThreshold, m_indexOfMax, m_kNBinsPerNSec, m_qNSecPerBin, m_shape, m_thresh, m_timeOfMax, DMR_cfg::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,
             evtSetup);              // 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.;
 
   m_deriv.resize(m_denseArraySize);
   m_shape.resize(m_denseArraySize);
 
   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;
 }

◆ derivative()

double EcalShapeBase::derivative ( double time ) const

Definition at line 117 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]);
 }

◆ fillShape()

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

◆ m_shape_print()

void EcalShapeBase::m_shape_print ( const char * fileName )

Definition at line 122 of file EcalShapeBase.cc.

References MillePedeFileConverter_cfg::fileName, compareTotals::fs, 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();
 }

◆ operator()()

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

overridevirtual

Implements CaloVShape.

Definition at line 110 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]);
 }

◆ setEventSetup()

void EcalShapeBase::setEventSetup ( const edm::EventSetup & evtSetup )

Definition at line 19 of file EcalShapeBase.cc.

References buildMe().

Referenced by EcalDigiProducer_Ph2::accumulate(), EcalDigiProducer::accumulate(), EcalTPGParamBuilder::analyze(), EcalUncalibRecHitWorkerWeights::set(), and EcalUncalibRecHitWorkerGlobal::set().

19 { buildMe(&evtSetup); }

EcalShapeBase::buildMe

void buildMe(const edm::EventSetup *=nullptr)

Definition: EcalShapeBase.cc:29

◆ threshold()

double EcalShapeBase::threshold ( ) const

Definition at line 27 of file EcalShapeBase.cc.

References m_thresh.

Referenced by buildMe().

27 { return m_thresh; }

EcalShapeBase::m_thresh

double m_thresh

Definition: EcalShapeBase.h:61

◆ timeIndex()

unsigned int EcalShapeBase::timeIndex ( double aTime ) const

protected

Definition at line 99 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);
 }