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 = 31

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 10 of file MTDShapeBase.h.

Member Typedef Documentation

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

Definition at line 14 of file MTDShapeBase.h.

Constructor & Destructor Documentation

MTDShapeBase::MTDShapeBase

(

)

Definition at line 8 of file MTDShapeBase.cc.

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

MTDShapeBase::~MTDShapeBase ( )

override

Definition at line 5 of file MTDShapeBase.cc.

{ }

Member Function Documentation

void MTDShapeBase::buildMe ( )

protected

Definition at line 101 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 89 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 145 of file MTDShapeBase.cc.

References a, b, Exception, and globals_cff::x1.

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 135 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 15 of file MTDShapeBase.cc.

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

Referenced by BTLElectronicsSim::run(), and timeToRise().

{

  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 121 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 95 of file MTDShapeBase.cc.

References timeOfMax_.

{
  return timeOfMax_; 
}

double MTDShapeBase::timeToRise ( ) const

inlineoverridevirtual

Implements CaloVShape.

Definition at line 24 of file MTDShapeBase.h.

References Scenarios_cff::scale, and timeAtThr().

{ return 0.; }

Member Data Documentation

unsigned int MTDShapeBase::indexOfMax_

private

Definition at line 51 of file MTDShapeBase.h.

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

constexpr unsigned int MTDShapeBase::k1NSecBinsTotal = kReadoutTimeInterval*kNBinsPerNSec

static

Definition at line 32 of file MTDShapeBase.h.

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

constexpr unsigned int MTDShapeBase::kNBinsPerNSec = 100

static

Definition at line 31 of file MTDShapeBase.h.

Referenced by timeIndex().

constexpr unsigned int MTDShapeBase::kReadoutTimeInterval = 31

static

Definition at line 30 of file MTDShapeBase.h.

const double MTDShapeBase::qNSecPerBin_

private

Definition at line 50 of file MTDShapeBase.h.

Referenced by buildMe(), and timeAtThr().

DVec MTDShapeBase::shape_

private

Definition at line 53 of file MTDShapeBase.h.

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

double MTDShapeBase::timeOfMax_

private

Definition at line 52 of file MTDShapeBase.h.

Referenced by buildMe(), and timeOfMax().