CSCAnalogSignal Class Reference

#include <CSCAnalogSignal.h>

Public Member Functions
	CSCAnalogSignal ()
	CSCAnalogSignal (int element, float binSize, std::vector< float > &binValues, float total=0., float timeOffset=0.)
float	getBinSize () const
float	getBinValue (int i) const
int	getElement () const
	constructor from time and amp shape More...
int	getSize () const
float	getTimeOffset () const
float	getTotal () const
float	getValue (float t) const
void	operator*= (float scaleFactor)
void	operator+= (float offset)
float &	operator[] (int i)
const float &	operator[] (int i) const
float	peakTime () const
	the time when the signal peaks More...
void	setElement (int element)
void	setTimeOffset (float offset)
unsigned	size () const
void	superimpose (const CSCAnalogSignal &signal2)

Private Attributes
float	invBinSize
std::vector< float >	theBinValues
int	theElement
float	theTimeOffset
float	theTotal

Friends
std::ostream &	operator<< (std::ostream &, const CSCAnalogSignal &)

Detailed Description

Simple histogram meant to represent the analog signal on a detector element.

Author

Rick Wilkinson

Last mod:
30-Jun-00 ptc Add further traps in getBinValue() and setBinValue().
06-Jul-00 ptc In fact the getBinValue trap was an important bug-fix: it trapped on > size() of stl std::vector but should have trapped >= size(). It occasionally does reach size().

Mods (performace improvements) by Vin 31/07/2000
Critical methods (getBinValue, get Value +=) inlined
bin-size stored and used as his inverse (encapulation helped in not changing interface, named changed to use compiler to catch its occurrencies)
swap input std::vector (be careful if const..)
do proper interpolation (not just /2)

Definition at line 32 of file CSCAnalogSignal.h.

Constructor & Destructor Documentation

CSCAnalogSignal::CSCAnalogSignal ( )

inline

Definition at line 34 of file CSCAnalogSignal.h.

: theElement(0), invBinSize(0.), theBinValues(0), theTotal(0), theTimeOffset(0.) {}

CSCAnalogSignal::CSCAnalogSignal ( int  element, float  binSize, std::vector< float > &  binValues, float  total = 0., float  timeOffset = 0.  )

inline

Definition at line 36 of file CSCAnalogSignal.h.

References theBinValues.

      : theElement(element), invBinSize(1. / binSize), theBinValues(), theTotal(total), theTimeOffset(timeOffset) {
    theBinValues.swap(binValues);
  }

Member Function Documentation

float CSCAnalogSignal::getBinSize ( ) const

inline

Definition at line 70 of file CSCAnalogSignal.h.

References invBinSize.

Referenced by CSCCrosstalkGenerator::getCrosstalk(), and operator<<().

{ return 1. / invBinSize; };

float CSCAnalogSignal::getBinValue ( int  i ) const

inline

Definition at line 48 of file CSCAnalogSignal.h.

References mps_fire::i, and theBinValues.

Referenced by CSCWireElectronicsSim::fillDigis(), getValue(), and CSCCrosstalkGenerator::ratio().

                                        {
    return (i >= static_cast<int>(theBinValues.size()) || i < 0) ? 0. : theBinValues[i];
  }

int CSCAnalogSignal::getElement ( ) const

inline

constructor from time and amp shape

Definition at line 46 of file CSCAnalogSignal.h.

References theElement.

Referenced by CSCBaseElectronicsSim::add(), and CSCStripConditions::noisify().

{ return theElement; };

int CSCAnalogSignal::getSize ( ) const

inline

Definition at line 69 of file CSCAnalogSignal.h.

References theBinValues.

Referenced by CSCWireElectronicsSim::fillDigis(), CSCCrosstalkGenerator::getCrosstalk(), operator*=(), operator+=(), operator[](), and CSCCrosstalkGenerator::ratio().

{ return theBinValues.size(); };

float CSCAnalogSignal::getTimeOffset ( ) const

inline

Definition at line 72 of file CSCAnalogSignal.h.

References theTimeOffset.

Referenced by CSCCrosstalkGenerator::getCrosstalk(), CSCStripConditions::noisify(), and operator<<().

{ return theTimeOffset; };

float CSCAnalogSignal::getTotal ( ) const

inline

Definition at line 71 of file CSCAnalogSignal.h.

References theTotal.

Referenced by SortSignalsByTotal().

{ return theTotal; };

float CSCAnalogSignal::getValue ( float  t ) const

inline

Definition at line 52 of file CSCAnalogSignal.h.

References f, getBinValue(), mps_fire::i, invBinSize, and theTimeOffset.

Referenced by CSCStripElectronicsSim::comparatorReading(), CSCStripElectronicsSim::createDigi(), CSCCrosstalkGenerator::getCrosstalk(), CSCStripElectronicsSim::makeNoiseSignal(), and superimpose().

                                       {
    // interpolate between bins, if necessary
    float retval = 0.;
    float f = (t - theTimeOffset) * invBinSize + 0.000000001;
    if (f >= 0.) {
      int i = static_cast<int>(f);
      f -= static_cast<float>(i);
      retval = (1. - f) * getBinValue(i) + f * getBinValue(i + 1);
    }
    return retval;
  }

void CSCAnalogSignal::operator*= ( float  scaleFactor )

inline

Definition at line 90 of file CSCAnalogSignal.h.

References getSize(), mps_fire::i, operator<<, JetEnergyShift_cfi::scaleFactor, theBinValues, and theTotal.

                                            {
    for (int i = 0; i < getSize(); ++i) {
      theBinValues[i] *= scaleFactor;
    }
    theTotal *= scaleFactor;
  }

void CSCAnalogSignal::operator+= ( float  offset )

inline

Definition at line 84 of file CSCAnalogSignal.h.

References getSize(), mps_fire::i, PFRecoTauDiscriminationByIsolation_cfi::offset, and theBinValues.

                                       {
    for (int i = 0; i < getSize(); ++i) {
      theBinValues[i] += offset;
    }
  }

float& CSCAnalogSignal::operator[] ( int  i )

inline

Definition at line 99 of file CSCAnalogSignal.h.

References getSize(), mps_fire::i, and theBinValues.

                           {
    assert(i >= 0 && i < getSize());
    return theBinValues[i];
  }

const float& CSCAnalogSignal::operator[] ( int  i ) const

inline

Definition at line 104 of file CSCAnalogSignal.h.

References getSize(), mps_fire::i, peakTime(), and theBinValues.

                                       {
    assert(i >= 0 && i < getSize());
    return theBinValues[i];
  }

float CSCAnalogSignal::peakTime

(

)

const

the time when the signal peaks

Definition at line 6 of file CSCAnalogSignal.cc.

References invBinSize, theBinValues, and theTimeOffset.

Referenced by operator[]().

                                      {
  size_t imax = std::max_element(theBinValues.begin(), theBinValues.end()) - theBinValues.begin();
  return imax / invBinSize + theTimeOffset;
}

void CSCAnalogSignal::setElement ( int  element )

inline

Definition at line 47 of file CSCAnalogSignal.h.

References theElement.

{ theElement = element; };

void CSCAnalogSignal::setTimeOffset ( float  offset )

inline

Definition at line 73 of file CSCAnalogSignal.h.

References PFRecoTauDiscriminationByIsolation_cfi::offset, and theTimeOffset.

{ theTimeOffset = offset; };

unsigned CSCAnalogSignal::size ( void  ) const

inline

Definition at line 111 of file CSCAnalogSignal.h.

References theBinValues.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

{ return theBinValues.size(); }

void CSCAnalogSignal::superimpose ( const CSCAnalogSignal &  signal2 )

inline

Definition at line 75 of file CSCAnalogSignal.h.

References getValue(), mps_fire::i, invBinSize, gen::n, lumiQTWidget::t, theBinValues, theTimeOffset, and theTotal.

Referenced by CSCBaseElectronicsSim::add(), CSCStripElectronicsSim::addCrosstalk(), and CSCStripConditions::noisify().

                                                          {
    size_t n = theBinValues.size();
    for (size_t i = 0; i < n; ++i) {
      float t = i / invBinSize + theTimeOffset;
      theBinValues[i] += signal2.getValue(t);
    }
    theTotal += signal2.theTotal;
  }

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  stream, const CSCAnalogSignal &  signal  )

friend

Definition at line 11 of file CSCAnalogSignal.cc.

Referenced by operator*=().

                                                                          {
  stream << "CSCAnalogSignal: Element " << signal.theElement << "   Total " << signal.theTotal << std::endl;
  for (int i = 0; i < int(signal.theBinValues.size()); ++i) {
    //@@ ptc 26-Feb-02 Don't both with very small amplitudes

    if (signal.theBinValues[i] > 1.E-10) {
      stream << i * signal.getBinSize() + signal.getTimeOffset() << "\t" << signal.theBinValues[i] << std::endl;
    }
  }
  return stream;
}

Member Data Documentation

float CSCAnalogSignal::invBinSize

private

Definition at line 115 of file CSCAnalogSignal.h.

Referenced by getBinSize(), getValue(), peakTime(), and superimpose().

std::vector<float> CSCAnalogSignal::theBinValues

private

Definition at line 116 of file CSCAnalogSignal.h.

Referenced by CSCAnalogSignal(), getBinValue(), getSize(), operator*=(), operator+=(), operator<<(), operator[](), peakTime(), size(), and superimpose().

int CSCAnalogSignal::theElement

private

Definition at line 114 of file CSCAnalogSignal.h.

Referenced by getElement(), operator<<(), and setElement().

float CSCAnalogSignal::theTimeOffset

private

Definition at line 118 of file CSCAnalogSignal.h.

Referenced by getTimeOffset(), getValue(), peakTime(), setTimeOffset(), and superimpose().

float CSCAnalogSignal::theTotal

private

Definition at line 117 of file CSCAnalogSignal.h.

Referenced by getTotal(), operator*=(), operator<<(), and superimpose().