MTDShapeBase Class Reference

#include <MTDShapeBase.h>

Inheritance diagram for MTDShapeBase:

Public Types
typedef std::vector< double >	DVec

Public Member Functions
float	fallTime () const
unsigned int	indexOfMax () const
float	maximum () 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
float	fallTime_
unsigned int	indexOfMax_
const double	qNSecPerBin_
DVec	shape_
double	timeOfMax_

Detailed Description

Definition at line 9 of file MTDShapeBase.h.

Member Typedef Documentation

DVec

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

Definition at line 11 of file MTDShapeBase.h.

Constructor & Destructor Documentation

MTDShapeBase()

MTDShapeBase::MTDShapeBase

(

)

Definition at line 6 of file MTDShapeBase.cc.

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

~MTDShapeBase()

MTDShapeBase::~MTDShapeBase ( )

override

Definition at line 4 of file MTDShapeBase.cc.

{}

Member Function Documentation

buildMe()

void MTDShapeBase::buildMe ( )

protected

Definition at line 81 of file MTDShapeBase.cc.

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

Referenced by BTLPulseShape::BTLPulseShape(), and ETLPulseShape::ETLPulseShape().

                           {
  // --- 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;
  }

  for (unsigned int i = indexOfMax_; i < shape_.size(); ++i) {
    if (shape_[indexOfMax_] * 0.01 > shape_[i]) {
      fallTime_ = i * qNSecPerBin_;
      break;
    }
  }

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

fallTime()

float MTDShapeBase::fallTime

(

)

const

Definition at line 79 of file MTDShapeBase.cc.

References fallTime_.

Referenced by ETLElectronicsSim::run().

{ return fallTime_; }

fillShape()

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

protectedpure virtual

Implemented in BTLPulseShape, and ETLPulseShape.

Referenced by buildMe().

indexOfMax()

unsigned int MTDShapeBase::indexOfMax

(

)

const

Definition at line 73 of file MTDShapeBase.cc.

References indexOfMax_.

{ return indexOfMax_; }

linear_interpolation()

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

private

Definition at line 119 of file MTDShapeBase.cc.

References a, b, Exception, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, y, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

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

maximum()

float MTDShapeBase::maximum

(

)

const

Definition at line 77 of file MTDShapeBase.cc.

References indexOfMax_, and shape_.

Referenced by ETLElectronicsSim::run().

{ return shape_[indexOfMax_]; }

operator()()

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

overridevirtual

Implements CaloVShape.

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

timeAtThr()

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

Definition at line 13 of file MTDShapeBase.cc.

References CustomPhysics_cfi::amplitude, mps_fire::i, indexOfMax_, linear_interpolation(), qNSecPerBin_, isotrackNtupler::scale, and shape_.

Referenced by BTLElectronicsSim::run(), and ETLElectronicsSim::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;
}

timeIndex()

unsigned int MTDShapeBase::timeIndex ( double  aTime ) const

protected

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

timeOfMax()

double MTDShapeBase::timeOfMax

(

)

const

Definition at line 75 of file MTDShapeBase.cc.

References timeOfMax_.

Referenced by ETLElectronicsSim::run().

{ return timeOfMax_; }

timeToRise()

double MTDShapeBase::timeToRise ( ) const

inlineoverridevirtual

Implements CaloVShape.

Definition at line 21 of file MTDShapeBase.h.

{ return 0.; }

Member Data Documentation

fallTime_

float MTDShapeBase::fallTime_

private

Definition at line 45 of file MTDShapeBase.h.

Referenced by buildMe(), and fallTime().

indexOfMax_

unsigned int MTDShapeBase::indexOfMax_

private

Definition at line 43 of file MTDShapeBase.h.

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

k1NSecBinsTotal

constexpr unsigned int MTDShapeBase::k1NSecBinsTotal = kReadoutTimeInterval * kNBinsPerNSec

static

Definition at line 29 of file MTDShapeBase.h.

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

kNBinsPerNSec

constexpr unsigned int MTDShapeBase::kNBinsPerNSec = 100

static

Definition at line 28 of file MTDShapeBase.h.

Referenced by timeIndex().

kReadoutTimeInterval

constexpr unsigned int MTDShapeBase::kReadoutTimeInterval = 28

static

Definition at line 27 of file MTDShapeBase.h.

qNSecPerBin_

const double MTDShapeBase::qNSecPerBin_

private

Definition at line 42 of file MTDShapeBase.h.

Referenced by buildMe(), and timeAtThr().

shape_

DVec MTDShapeBase::shape_

private

Definition at line 46 of file MTDShapeBase.h.

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

timeOfMax_

double MTDShapeBase::timeOfMax_

private

Definition at line 44 of file MTDShapeBase.h.

Referenced by buildMe(), and timeOfMax().