LArWave Class Reference

#include <LArWave.h>

Inheritance diagram for LArWave:

Collaboration diagram for LArWave:

Public Member Functions
structors
	LArWave ()
	LArWave (const unsigned nSamples, const double dt, const unsigned flag=0)
	LArWave (const std::vector< double > &theVector, const double dt, const unsigned flag=0)
virtual	~LArWave ()=default
accumulation accessors
const std::vector< double > &	getWave () const
	Wave parameters.
const double &	getDt () const
	delta time
const double &	getSample (const unsigned int i) const
	Amplitude per time bin.
const double &	operator[] (const unsigned int i) const
	Amplitude per time bin.
double	getTime (const unsigned i) const
	time
size_t	getSize () const
	number of time samples
unsigned	getIndex (double aTime) const
	index for a time value
unsigned	getFlag () const
	flag: ...
bool	isEmpty () const
	is LArWave uninitialized?
accumulation setters
double &	operator[] (const unsigned int i)
	Amplitude per time bin - NOT const.
void	setSize (const unsigned nSize)
	resize the number of time bin samples
void	setDt (const double dt)
	set the delta time
void	setSample (const unsigned i, const double aVal)
	set the amplitude for time bin i
void	setFlag (const unsigned flag)
	set flag

accumulation alegebra operations
enum	{   meas = 01 , mwf = 02 , mwfN = 03 , dac0 = 04 ,   xtalkcorr = 05 , pred = 10 , predCali = 11 , predFitPhys = 12 ,   predMeasPar = 13 , model = 20 , patch = 30 , unknown = 0 }
double	m_dt
std::vector< double >	m_amplitudes
unsigned	m_flag
LArWave	operator+ (const LArWave &bWave) const
LArWave	operator- (const LArWave &bWave) const
LArWave	operator* (const LArWave &bWave) const
LArWave	operator/ (const LArWave &bWave) const
LArWave	operator% (const LArWave &bWave) const
LArWave	operator+ (const double aBias) const
LArWave	operator* (const double aScale) const
LArWave &	operator+= (const LArWave &bWave)
LArWave &	operator-= (const LArWave &bWave)
LArWave &	operator*= (const double aScale)
	LArWave (const LArWave &theWave1, const LArWave &theWave2)

Detailed Description

Definition at line 31 of file LArWave.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
meas
mwf
mwfN
dac0
xtalkcorr
pred
predCali
predFitPhys
predMeasPar
model
patch
unknown

Definition at line 121 of file LArWave.h.

         { 
    meas        = 01 ,  // measured (cumulated from several triggers)
    mwf         = 02 ,  // calibration master waveform
    mwfN        = 03 ,  // as mwf, but normalized to have peak=1
    dac0        = 04 ,  // dac0 wave
    xtalkcorr   = 05 ,  // cross-talk corrected wave (Strips)
 
    pred        = 10 ,  // predicted from observed calibration (generic)
    predCali    = 11 ,  // as pred, parameters calibration only (a` la Milano)
    predFitPhys = 12 ,  // as pred, parameters from fit to observed phys (a` la Annecy)
    predMeasPar = 13 ,  // as pred, parameters from independent measurements
 
    model       = 20 ,  // from some model (generic)
 
    patch       = 30 ,  // patched from another channel
 
    unknown     =  0 
    };

Constructor & Destructor Documentation

LArWave() [1/4]

LArWave::LArWave ( )

inline

Definition at line 147 of file LArWave.h.

                 : m_dt(0), m_flag(0)
{} 

LArWave() [2/4]

LArWave::LArWave ( const unsigned nSamples, const double dt, const unsigned flag = 0 )

inline

Definition at line 152 of file LArWave.h.

    : 
    m_dt(dt) , m_flag(flag)
{ m_amplitudes.resize(nSamples) ; }

LArWave() [3/4]

LArWave::LArWave ( const std::vector< double > & theVector, const double dt, const unsigned flag = 0 )

inline

Definition at line 159 of file LArWave.h.

    : 
    m_dt(dt) , m_amplitudes(theVector) , m_flag(flag)
{} 

~LArWave()

virtual LArWave::~LArWave ( )

virtualdefault

LArWave() [4/4]

LArWave::LArWave ( const LArWave & theWave1, const LArWave & theWave2 )

protected

Definition at line 9 of file LArWave.cxx.

  : m_flag(0)
{
  if ( theWave1.m_dt == theWave2.m_dt ) {
    m_dt = theWave1.m_dt ;
    m_amplitudes.resize( std::min( theWave1.getSize(),
                   theWave2.getSize() )) ; 
  } else {
    m_amplitudes.resize(0);
    m_dt = 0. ;
  }
} 

Member Function Documentation

getDt()

const double & LArWave::getDt ( ) const

inline

delta time

Definition at line 50 of file LArWave.h.

{ return m_dt; }

getFlag()

unsigned LArWave::getFlag ( ) const

inline

flag: ...

Definition at line 178 of file LArWave.h.

{ return m_flag ; } 

getIndex()

unsigned LArWave::getIndex

(

double

aTime

)

const

index for a time value

Definition at line 146 of file LArWave.cxx.

{ return (aTime>=0 && m_dt>0) ? (unsigned)(aTime/m_dt) : getSize()+1 ; } 

getSample()

const double & LArWave::getSample ( const unsigned int i ) const

inline

Amplitude per time bin.

Definition at line 53 of file LArWave.h.

{return m_amplitudes[i];}

getSize()

size_t LArWave::getSize ( ) const

inline

number of time samples

Definition at line 62 of file LArWave.h.

{ return m_amplitudes.size() ; } 

getTime()

double LArWave::getTime ( const unsigned i ) const

inline

time

Definition at line 172 of file LArWave.h.

{ return i * m_dt ; } 

getWave()

const std::vector< double > & LArWave::getWave ( ) const

inline

Wave parameters.

Definition at line 167 of file LArWave.h.

{ return m_amplitudes ; } 

isEmpty()

bool LArWave::isEmpty ( ) const

inline

is LArWave uninitialized?

Definition at line 183 of file LArWave.h.

{ return (m_amplitudes.size() == 0); } 

operator%()

LArWave LArWave::operator%

(

const LArWave &

bWave

)

const

Definition at line 91 of file LArWave.cxx.

{
  //W.L., 2-Sept-09: Speed-up: 
  //This method is called several million times
  //aggressive optimization pays off.
 
  const double* amplPtrA=&(this->m_amplitudes.front());
  const double* amplPtrB=&(bWave.m_amplitudes.front());
    
  LArWave result(*this,bWave) ;
  const size_t s=result.getSize();
  for (size_t i=0 ; i<s ; ++i ) {
    //double& resSample=result.m_amplitudes[i];
    double sum2 = 0.5 * ( (*this).m_amplitudes[0] * bWave.m_amplitudes[i] +
              (*this).m_amplitudes[i] * bWave.m_amplitudes[0] ) ;
    for (size_t k=1 ; k<i ; ++k ) {
      sum2 += amplPtrA[k] * amplPtrB[i-k] ;
    }
    result.m_amplitudes[i] = sum2*result.m_dt ;
  }
  return result ;
}

operator*() [1/2]

LArWave LArWave::operator*

(

const double

aScale

)

const

Definition at line 126 of file LArWave.cxx.

                                                   {
  LArWave result(*this) ;
  std::vector<double>::iterator it=result.m_amplitudes.begin();
  std::vector<double>::iterator it_e=result.m_amplitudes.end();
  for (;it!=it_e;++it) {
    (*it)*=aScale;
  }
  return result ;
}

operator*() [2/2]

LArWave LArWave::operator*

(

const LArWave &

bWave

)

const

Definition at line 69 of file LArWave.cxx.

{ 
  LArWave result(*this,bWave) ;
  const size_t s=result.getSize();
  for ( size_t i=0 ; i<s ; ++i ) {
    result.m_amplitudes[i] = 
      (*this).m_amplitudes[i] * bWave.m_amplitudes[i] ;
  }
  return result ;
}

operator*=()

LArWave & LArWave::operator*=

(

const double

aScale

)

Definition at line 137 of file LArWave.cxx.

                                                {
  for (double& a : m_amplitudes) {
    a*=aScale;
  }
  return *this ;
}

operator+() [1/2]

LArWave LArWave::operator+

(

const double

aBias

)

const

Definition at line 115 of file LArWave.cxx.

                                                   {
  LArWave result(*this) ;
  std::vector<double>::iterator it=result.m_amplitudes.begin();
  std::vector<double>::iterator it_e=result.m_amplitudes.end();
  for (;it!=it_e;++it) {
    (*it)+=aBias;
  }
  return result ;
}

operator+() [2/2]

LArWave LArWave::operator+

(

const LArWave &

bWave

)

const

Definition at line 24 of file LArWave.cxx.

{ 
  LArWave result(*this,bWave) ;
  const size_t s=result.getSize();
  for ( size_t i=0 ; i<s ; ++i ) {
    result.m_amplitudes[i] = 
      (*this).m_amplitudes[i] + bWave.m_amplitudes[i] ;
  }
  return result ;
}

operator+=()

LArWave & LArWave::operator+=

(

const LArWave &

bWave

)

Definition at line 35 of file LArWave.cxx.

{ 
  const size_t s=std::min(m_amplitudes.size(), bWave.getSize());
  for ( size_t i=0 ; i<s ; ++i ) {
    m_amplitudes[i] += bWave.m_amplitudes[i] ;
  }
  return *this;
}

operator-()

LArWave LArWave::operator-

(

const LArWave &

bWave

)

const

Definition at line 46 of file LArWave.cxx.

{ 
  LArWave result(*this,bWave) ;
  const size_t s=result.getSize();
  for ( size_t i=0 ; i<s ; ++i ) {
    result.m_amplitudes[i] = 
      (*this).m_amplitudes[i] - bWave.m_amplitudes[i] ;
  }
  return result ;
}

operator-=()

LArWave & LArWave::operator-=

(

const LArWave &

bWave

)

Definition at line 58 of file LArWave.cxx.

{ 
  const size_t s=std::min(m_amplitudes.size(), bWave.getSize());
  for ( size_t i=0 ; i<s ; ++i ) {
    m_amplitudes[i] -= bWave.m_amplitudes[i] ;
  }
  return *this;
}

operator/()

LArWave LArWave::operator/

(

const LArWave &

bWave

)

const

Definition at line 80 of file LArWave.cxx.

{ 
  LArWave result(*this,bWave) ;
  const size_t s=result.getSize();
  for ( size_t i=0 ; i<s ; ++i ) {
    result.m_amplitudes[i] = 
      (*this).m_amplitudes[i] / bWave.m_amplitudes[i] ;
  }
  return result ;
}

operator[]() [1/2]

double & LArWave::operator[] ( const unsigned int i )

inline

Amplitude per time bin - NOT const.

Definition at line 78 of file LArWave.h.

{return m_amplitudes[i];}

operator[]() [2/2]

const double & LArWave::operator[] ( const unsigned int i ) const

inline

Amplitude per time bin.

Definition at line 56 of file LArWave.h.

{return m_amplitudes[i];}

setDt()

void LArWave::setDt ( const double dt )

inline

set the delta time

Definition at line 193 of file LArWave.h.

{ m_dt = dt ; } 

setFlag()

void LArWave::setFlag ( const unsigned flag )

inline

set flag

Definition at line 199 of file LArWave.h.

{ m_flag = flag ; } 

setSample()

void LArWave::setSample ( const unsigned i, const double aVal )

inline

set the amplitude for time bin i

Definition at line 87 of file LArWave.h.

{ m_amplitudes[i] = aVal ; } 

setSize()

void LArWave::setSize ( const unsigned nSize )

inline

resize the number of time bin samples

Definition at line 188 of file LArWave.h.

{ m_amplitudes.resize(nSize,0) ; } 

Member Data Documentation

m_amplitudes

std::vector<double> LArWave::m_amplitudes

protected

Definition at line 112 of file LArWave.h.

m_dt

double LArWave::m_dt

protected

Definition at line 111 of file LArWave.h.

m_flag

unsigned LArWave::m_flag

protected

Definition at line 113 of file LArWave.h.

---

The documentation for this class was generated from the following files:

• LArWave.h
• LArWave.cxx