LArWaveCumul Class Reference

#include <LArWaveCumul.h>

Inheritance diagram for LArWaveCumul:

Collaboration diagram for LArWaveCumul:

Public Member Functions
structors
	LArWaveCumul ()
	LArWaveCumul (const std::vector< double > &theVector, double dt, unsigned flag=0)
	LArWaveCumul (const std::vector< double > &vAmpl, const std::vector< double > &vErr, const std::vector< int > &vTrig, double dt, unsigned flag=0)
	LArWaveCumul (unsigned nSamples, double dt, unsigned flag=0)
virtual	~LArWaveCumul ()=default
accumulation accessors
const std::vector< double > &	getErrors () const
	error vector
double	getError (unsigned i) const
	error per time bin
const std::vector< int > &	getTriggers () const
	trigger vector
int	getTrigger (unsigned i) const
	trigger per time bin
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 setters
std::vector< double >	m_errors
std::vector< int >	m_triggers
void	addEvent (int delay, int step, const std::vector< double > &Samples)
	add in set of channel responses for a given step/delay
void	addAccumulatedEvent (int delay, int step, const std::vector< double > &SamplesSum, const std::vector< double > &Samples2Sum, unsigned nTriggers)
	used when building from LArAccumulatedDigits

accumulation alegebra operations
double	m_dt
std::vector< double >	m_amplitudes
unsigned	m_flag
LArWave	operator+ (const LArWave &bWave) const
LArWave	operator+ (const double aBias) const
LArWave	operator- (const LArWave &bWave) const
LArWave	operator* (const LArWave &bWave) const
LArWave	operator* (const double aScale) const
LArWave	operator/ (const LArWave &bWave) const
LArWave	operator% (const LArWave &bWave) const
LArWave &	operator+= (const LArWave &bWave)
LArWave &	operator-= (const LArWave &bWave)
LArWave &	operator*= (const double aScale)
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 }

Detailed Description

Definition at line 30 of file LArWaveCumul.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

inherited

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

LArWaveCumul() [1/4]

LArWaveCumul::LArWaveCumul ( )

inline

Definition at line 96 of file LArWaveCumul.h.

    : 
    LArWave() 
{}

LArWaveCumul() [2/4]

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

inline

Definition at line 102 of file LArWaveCumul.h.

    :
    LArWave(theVector,dt,flag)
{
    unsigned nSamples = theVector.size();
    m_errors.resize(nSamples,0); 
    m_triggers.resize(nSamples,0);
}

LArWaveCumul() [3/4]

LArWaveCumul::LArWaveCumul ( const std::vector< double > & vAmpl, const std::vector< double > & vErr, const std::vector< int > & vTrig, double dt, unsigned flag = 0 )

inline

Definition at line 114 of file LArWaveCumul.h.

    :
    LArWave(vAmpl,dt,flag), 
    m_errors(vErr),
    m_triggers(vTrig)
{} 

LArWaveCumul() [4/4]

LArWaveCumul::LArWaveCumul ( unsigned nSamples, double dt, unsigned flag = 0 )

inline

Definition at line 126 of file LArWaveCumul.h.

    :
    LArWave(nSamples,dt,flag)
{
    m_errors.resize(nSamples,0) ; 
    m_triggers.resize(nSamples,0) ;
}

~LArWaveCumul()

virtual LArWaveCumul::~LArWaveCumul ( )

virtualdefault

Member Function Documentation

addAccumulatedEvent()

void LArWaveCumul::addAccumulatedEvent	(	int	delay,
		int	step,
		const std::vector< double > &	SamplesSum,
		const std::vector< double > &	Samples2Sum,
		unsigned	nTriggers )

used when building from LArAccumulatedDigits

Definition at line 39 of file LArWaveCumul.cxx.

{
  if (!nTriggers) return;  // there should be related warnings displayed by the LArLATOMEDecoder
  const unsigned int nSwave = getSize() ;
  const unsigned int nSevt  = SamplesSum.size() ;
  
  for ( unsigned int i=0 ; i<nSevt ; i++ ) {
 
    unsigned int k=(i+1)*step-(delay+1);
    
    if (k<nSwave ) {
 
      int    N1    = m_triggers[k] ;
      double sum1  = m_amplitudes[k] * N1 ;
      double sum12 = (N1-1)*m_errors[k]*m_errors[k] + N1*m_amplitudes[k]*m_amplitudes[k] ;
 
      int    N2    = nTriggers ;
      double sum2  = SamplesSum[i] ;
      double sum22 = Samples2Sum[i] ;
      
      m_triggers[k]   = N1+N2 ;
      m_amplitudes[k] = (sum1+sum2)/m_triggers[k] ;
      m_errors[k]     = std::sqrt(std::max(sum12+sum22-m_amplitudes[k]*m_amplitudes[k]*m_triggers[k], 0.)/(m_triggers[k]-1) ) ;
 
    }
  }
}

addEvent()

void LArWaveCumul::addEvent	(	int	delay,
		int	step,
		const std::vector< double > &	Samples )

add in set of channel responses for a given step/delay

Definition at line 10 of file LArWaveCumul.cxx.

{
  const unsigned int nSwave = getSize() ;
  const unsigned int nSevt  = Samples.size() ;
  for ( unsigned int i=0 ; i<nSevt ; i++ ) {
    unsigned int k=(i+1)*step-(delay+1);
    if (k<nSwave) {
      double evt  = Samples[i] ;
      int    N    = m_triggers[k] ;
      double amp2 = m_amplitudes[k]*m_amplitudes[k];
      double sum2 = (N-1)*m_errors[k]*m_errors[k] + N*amp2;
 
      double sum  = m_amplitudes[k] * N ;
      m_amplitudes[k] = ( sum + evt ) / (N+1) ;
      double arg = N? ( ( sum2 + evt*evt ) - (N+1)*m_amplitudes[k]*m_amplitudes[k] )/N: 0.;
 
      // Can happen due to underflow if m_errors is zero and evt is
      // exactly the same as m_amplitudes.
      if (arg < 0) arg = 0;
      m_errors[k]     = sqrt(arg);
      m_triggers[k]   = N + 1 ;
 
    }
  }
}

getDt()

const double & LArWave::getDt ( ) const

inlineinherited

delta time

Definition at line 50 of file LArWave.h.

{ return m_dt; }

getError()

double LArWaveCumul::getError ( unsigned i ) const

inline

error per time bin

Definition at line 143 of file LArWaveCumul.h.

{ 
    return (i<getSize()) ? m_errors[i] : -9.9999E+99 ; 
}

getErrors()

const std::vector< double > & LArWaveCumul::getErrors ( ) const

inline

error vector

Definition at line 138 of file LArWaveCumul.h.

{ return m_errors ; }

getFlag()

unsigned LArWave::getFlag ( ) const

inlineinherited

flag: ...

Definition at line 178 of file LArWave.h.

{ return m_flag ; } 

getIndex()

unsigned LArWave::getIndex ( double aTime ) const

inherited

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

inlineinherited

Amplitude per time bin.

Definition at line 53 of file LArWave.h.

{return m_amplitudes[i];}

getSize()

size_t LArWave::getSize ( ) const

inlineinherited

number of time samples

Definition at line 62 of file LArWave.h.

{ return m_amplitudes.size() ; } 

getTime()

double LArWave::getTime ( const unsigned i ) const

inlineinherited

time

Definition at line 172 of file LArWave.h.

{ return i * m_dt ; } 

getTrigger()

int LArWaveCumul::getTrigger ( unsigned i ) const

inline

trigger per time bin

Definition at line 157 of file LArWaveCumul.h.

{ return (i<getSize()) ? m_triggers[i] : -999999999; }

getTriggers()

const std::vector< int > & LArWaveCumul::getTriggers ( ) const

inline

trigger vector

Definition at line 150 of file LArWaveCumul.h.

{
    return m_triggers;
}

getWave()

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

inlineinherited

Wave parameters.

Definition at line 167 of file LArWave.h.

{ return m_amplitudes ; } 

isEmpty()

bool LArWave::isEmpty ( ) const

inlineinherited

is LArWave uninitialized?

Definition at line 183 of file LArWave.h.

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

operator%()

LArWave LArWave::operator% ( const LArWave & bWave ) const

inherited

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

inherited

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

inherited

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 )

inherited

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

inherited

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

inherited

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 )

inherited

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

inherited

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 )

inherited

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

inherited

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 )

inlineinherited

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

inlineinherited

Amplitude per time bin.

Definition at line 56 of file LArWave.h.

{return m_amplitudes[i];}

setDt()

void LArWave::setDt ( const double dt )

inlineinherited

set the delta time

Definition at line 193 of file LArWave.h.

{ m_dt = dt ; } 

setFlag()

void LArWave::setFlag ( const unsigned flag )

inlineinherited

set flag

Definition at line 199 of file LArWave.h.

{ m_flag = flag ; } 

setSample()

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

inlineinherited

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 )

inlineinherited

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

protectedinherited

Definition at line 112 of file LArWave.h.

m_dt

double LArWave::m_dt

protectedinherited

Definition at line 111 of file LArWave.h.

m_errors

std::vector<double> LArWaveCumul::m_errors

protected

Definition at line 89 of file LArWaveCumul.h.

m_flag

unsigned LArWave::m_flag

protectedinherited

Definition at line 113 of file LArWave.h.

m_triggers

std::vector<int> LArWaveCumul::m_triggers

protected

Definition at line 90 of file LArWaveCumul.h.

---

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

• LArWaveCumul.h
• LArWaveCumul.cxx