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

{}

Member Function Documentation

buildMe()

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

protected

Definition at line 29 of file EcalShapeBase.cc.

References a, b, HltBtagPostValidation_cff::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().

{ buildMe(&evtSetup); }

threshold()

double EcalShapeBase::threshold

(

)

const

Definition at line 27 of file EcalShapeBase.cc.

References m_thresh.

Referenced by buildMe().

{ return m_thresh; }

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

timeOfMax()

double EcalShapeBase::timeOfMax

(

)

const

Definition at line 23 of file EcalShapeBase.cc.

References m_timeOfMax.

Referenced by timeToRise().

{ return m_timeOfMax; }

timeOfThr()

double EcalShapeBase::timeOfThr

(

)

const

Definition at line 21 of file EcalShapeBase.cc.

References m_firstTimeOverThreshold.

Referenced by timeToRise().

{ return m_firstTimeOverThreshold; }

timeToRise()

double EcalShapeBase::timeToRise ( ) const

overridevirtual

Implements CaloVShape.

Definition at line 25 of file EcalShapeBase.cc.

References timeOfMax(), and timeOfThr().

Referenced by EcalUncalibRecHitRecChi2Algo< C >::EcalUncalibRecHitRecChi2Algo().

{ return timeOfMax() - timeOfThr(); }

Member Data Documentation

m_arraySize

unsigned int EcalShapeBase::m_arraySize

private

Definition at line 66 of file EcalShapeBase.h.

Referenced by buildMe().

m_denseArraySize

unsigned int EcalShapeBase::m_denseArraySize

private

Definition at line 67 of file EcalShapeBase.h.

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

m_deriv

DVec EcalShapeBase::m_deriv

private

Definition at line 64 of file EcalShapeBase.h.

Referenced by buildMe(), and derivative().

m_firstIndexOverThreshold

unsigned int EcalShapeBase::m_firstIndexOverThreshold

private

Definition at line 57 of file EcalShapeBase.h.

Referenced by buildMe(), and timeIndex().

m_firstTimeOverThreshold

double EcalShapeBase::m_firstTimeOverThreshold

private

Definition at line 58 of file EcalShapeBase.h.

Referenced by buildMe(), and timeOfThr().

m_indexOfMax

unsigned int EcalShapeBase::m_indexOfMax

private

Definition at line 59 of file EcalShapeBase.h.

Referenced by buildMe().

m_kNBinsPerNSec

unsigned int EcalShapeBase::m_kNBinsPerNSec

private

Definition at line 62 of file EcalShapeBase.h.

Referenced by buildMe(), and timeIndex().

m_qNSecPerBin

double EcalShapeBase::m_qNSecPerBin

private

Definition at line 68 of file EcalShapeBase.h.

Referenced by buildMe().

m_shape

DVec EcalShapeBase::m_shape

private

Definition at line 63 of file EcalShapeBase.h.

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

m_thresh

double EcalShapeBase::m_thresh

private

Definition at line 61 of file EcalShapeBase.h.

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

m_timeOfMax

double EcalShapeBase::m_timeOfMax

private

Definition at line 60 of file EcalShapeBase.h.

Referenced by buildMe(), and timeOfMax().

m_useDBShape

bool EcalShapeBase::m_useDBShape

protected

Definition at line 54 of file EcalShapeBase.h.

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