MTDShapeBase Class Reference

#include <MTDShapeBase.h>

Inheritance diagram for MTDShapeBase:

Public Types
typedef std::vector< double >	DVec

Public Member Functions
unsigned int	indexOfMax () const
	MTDShapeBase ()
double	operator() (double aTime) const override
std::array< float, 3 >	timeAtThr (const float scale, const float threshold1, const float threshold2) const
double	timeOfMax () const
double	timeToRise () const override
	~MTDShapeBase () override
Public Member Functions inherited from CaloVShape
	CaloVShape ()
virtual	~CaloVShape ()

Static Public Attributes
static constexpr unsigned int	k1NSecBinsTotal = kReadoutTimeInterval * kNBinsPerNSec
static constexpr unsigned int	kNBinsPerNSec = 100
static constexpr unsigned int	kReadoutTimeInterval = 28

Protected Member Functions
void	buildMe ()
virtual void	fillShape (DVec &aVec) const =0
unsigned int	timeIndex (double aTime) const

Private Member Functions
double	linear_interpolation (const double &y, const double &x1, const double &x2, const double &y1, const double &y2) const

Private Attributes
unsigned int	indexOfMax_
const double	qNSecPerBin_
DVec	shape_
double	timeOfMax_

Detailed Description

Definition at line 9 of file MTDShapeBase.h.

Member Typedef Documentation

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

Definition at line 11 of file MTDShapeBase.h.

Constructor & Destructor Documentation

MTDShapeBase::MTDShapeBase

(

)

Definition at line 6 of file MTDShapeBase.cc.

    : qNSecPerBin_(1. / kNBinsPerNSec), indexOfMax_(0), timeOfMax_(0.), shape_(DVec(k1NSecBinsTotal, 0.0)) {}

MTDShapeBase::~MTDShapeBase ( )

override

Definition at line 4 of file MTDShapeBase.cc.

{}

Member Function Documentation

void MTDShapeBase::buildMe ( )

protected

Definition at line 73 of file MTDShapeBase.cc.

References fillShape(), mps_fire::i, indexOfMax_, qNSecPerBin_, shape_, and timeOfMax_.

Referenced by BTLPulseShape::BTLPulseShape().

                           {
  // --- Fill the vector with the pulse shape
  fillShape(shape_);

  // --- Find the index of maximum
  for (unsigned int i = 0; i < shape_.size(); ++i) {
    if (shape_[indexOfMax_] < shape_[i])
      indexOfMax_ = i;
  }

  if (indexOfMax_ != 0)
    timeOfMax_ = indexOfMax_ * qNSecPerBin_;
}

virtual void MTDShapeBase::fillShape ( DVec &  aVec ) const

protectedpure virtual

Implemented in BTLPulseShape.

Referenced by buildMe().

unsigned int MTDShapeBase::indexOfMax

(

)

const

Definition at line 69 of file MTDShapeBase.cc.

References indexOfMax_.

{ return indexOfMax_; }

double MTDShapeBase::linear_interpolation ( const double &  y, const double &  x1, const double &  x2, const double &  y1, const double &  y2  ) const

private

Definition at line 104 of file MTDShapeBase.cc.

References a, b, and Exception.

Referenced by timeAtThr().

                                                                                                   {
  if (x1 == x2)
    throw cms::Exception("BadValue") << " MTDShapeBase: Trying to interpolate two points with the same x coordinate!";

  double a = (y2 - y1) / (x2 - x1);
  double b = y1 - a * x1;

  return (y - b) / a;
}

double MTDShapeBase::operator() ( double  aTime ) const

overridevirtual

Implements CaloVShape.

Definition at line 98 of file MTDShapeBase.cc.

References k1NSecBinsTotal, shape_, and timeIndex().

                                                  {
  // return pulse amplitude for request time in ns
  const unsigned int index(timeIndex(aTime));
  return (k1NSecBinsTotal == index ? 0 : shape_[index]);
}

std::array< float, 3 > MTDShapeBase::timeAtThr	(	const float	scale,
		const float	threshold1,
		const float	threshold2
	)		const

Definition at line 9 of file MTDShapeBase.cc.

References mps_fire::i, indexOfMax_, linear_interpolation(), qNSecPerBin_, pileupReCalc_HLTpaths::scale, and shape_.

Referenced by BTLElectronicsSim::run().

                                                                                                                  {
  std::array<float, 3> times_tmp = {{0., 0., 0.}};

  // --- Check if the pulse amplitude is greater than threshold 2
  if (shape_[indexOfMax_] * scale < threshold2)
    return times_tmp;

  // --- Find the times corresponding to thresholds 1 and 2 on the pulse leading edge
  //     NB: To speed up the search we loop only on the rising edge
  unsigned int index_LT1 = 0;
  unsigned int index_LT2 = 0;

  for (unsigned int i = 0; i < indexOfMax_; ++i) {
    float amplitude = shape_[i] * scale;

    if (amplitude > threshold1 && index_LT1 == 0)
      index_LT1 = i;

    if (amplitude > threshold2 && index_LT2 == 0) {
      index_LT2 = i;
      break;
    }
  }

  // --- Find the time corresponding to thresholds 1 on the pulse falling edge
  unsigned int index_FT1 = 0;

  for (unsigned int i = shape_.size() - 1; i > indexOfMax_; i--) {
    float amplitude = shape_[i] * scale;

    if (amplitude > threshold1 && index_FT1 == 0) {
      index_FT1 = i + 1;
      break;
    }
  }

  if (index_LT1 != 0)
    times_tmp[0] = linear_interpolation(threshold1,
                                        (index_LT1 - 1) * qNSecPerBin_,
                                        index_LT1 * qNSecPerBin_,
                                        shape_[index_LT1 - 1] * scale,
                                        shape_[index_LT1] * scale);

  if (index_LT2 != 0)
    times_tmp[1] = linear_interpolation(threshold2,
                                        (index_LT2 - 1) * qNSecPerBin_,
                                        index_LT2 * qNSecPerBin_,
                                        shape_[index_LT2 - 1] * scale,
                                        shape_[index_LT2] * scale);

  if (index_FT1 != 0)
    times_tmp[2] = linear_interpolation(threshold1,
                                        (index_FT1 - 1) * qNSecPerBin_,
                                        index_FT1 * qNSecPerBin_,
                                        shape_[index_FT1 - 1] * scale,
                                        shape_[index_FT1] * scale);

  return times_tmp;
}

unsigned int MTDShapeBase::timeIndex ( double  aTime ) const

protected

Definition at line 87 of file MTDShapeBase.cc.

References k1NSecBinsTotal, kNBinsPerNSec, and LogDebug.

Referenced by operator()().

                                                       {
  const unsigned int index = aTime * kNBinsPerNSec;

  const bool bad = (k1NSecBinsTotal <= index);

  if (bad)
    LogDebug("MTDShapeBase") << " MTD pulse shape requested for out of range time " << aTime;

  return (bad ? k1NSecBinsTotal : index);
}

double MTDShapeBase::timeOfMax

(

)

const

Definition at line 71 of file MTDShapeBase.cc.

References timeOfMax_.

{ return timeOfMax_; }

double MTDShapeBase::timeToRise ( ) const

inlineoverridevirtual

Implements CaloVShape.

Definition at line 21 of file MTDShapeBase.h.

{ return 0.; }

Member Data Documentation

unsigned int MTDShapeBase::indexOfMax_

private

Definition at line 41 of file MTDShapeBase.h.

Referenced by buildMe(), indexOfMax(), and timeAtThr().

constexpr unsigned int MTDShapeBase::k1NSecBinsTotal = kReadoutTimeInterval * kNBinsPerNSec

static

Definition at line 27 of file MTDShapeBase.h.

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

constexpr unsigned int MTDShapeBase::kNBinsPerNSec = 100

static

Definition at line 26 of file MTDShapeBase.h.

Referenced by timeIndex().

constexpr unsigned int MTDShapeBase::kReadoutTimeInterval = 28

static

Definition at line 25 of file MTDShapeBase.h.

const double MTDShapeBase::qNSecPerBin_

private

Definition at line 40 of file MTDShapeBase.h.

Referenced by buildMe(), and timeAtThr().

DVec MTDShapeBase::shape_

private

Definition at line 43 of file MTDShapeBase.h.

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

double MTDShapeBase::timeOfMax_

private

Definition at line 42 of file MTDShapeBase.h.

Referenced by buildMe(), and timeOfMax().