LArCaliWave Class Reference

#include <LArCaliWave.h>

Inheritance diagram for LArCaliWave:

Collaboration diagram for LArCaliWave:

Public Member Functions
structors
	LArCaliWave ()
	LArCaliWave (const std::vector< double > &theVector, double dt, int DAC, int isPulsed, unsigned flag)
	LArCaliWave (unsigned nSamples, double dt, int DAC, int isPulsed, unsigned flag)
	LArCaliWave (const std::vector< double > &vAmpl, const std::vector< double > &vErr, const std::vector< int > &vTrig, double dt, int DAC, int isPulsed, unsigned flag)
virtual	~LArCaliWave ()=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 accessors
int	m_DAC
int	m_isPulsed
int	getDAC () const
	DAC value.
int	getIsPulsedInt () const
	isPulsed value

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 44 of file LArCaliWave.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

LArCaliWave() [1/4]

LArCaliWave::LArCaliWave ( )

inline

Definition at line 108 of file LArCaliWave.h.

  : LArWaveCumul(), m_DAC(0), m_isPulsed(false)
{}

LArCaliWave() [2/4]

LArCaliWave::LArCaliWave ( const std::vector< double > & theVector, double dt, int DAC, int isPulsed, unsigned flag )

inline

Definition at line 113 of file LArCaliWave.h.

    :  
    LArWaveCumul(theVector,dt,flag),
    m_DAC(DAC),m_isPulsed(isPulsed) 
{}

LArCaliWave() [3/4]

LArCaliWave::LArCaliWave ( unsigned nSamples, double dt, int DAC, int isPulsed, unsigned flag )

inline

Definition at line 124 of file LArCaliWave.h.

    :  
    LArWaveCumul(nSamples,dt,flag), 
    m_DAC(DAC),m_isPulsed(isPulsed)
{}

LArCaliWave() [4/4]

LArCaliWave::LArCaliWave ( const std::vector< double > & vAmpl, const std::vector< double > & vErr, const std::vector< int > & vTrig, double dt, int DAC, int isPulsed, unsigned flag )

inline

Definition at line 135 of file LArCaliWave.h.

    :  
    LArWaveCumul(vAmpl,vErr,vTrig,dt,flag), 
    m_DAC(DAC),m_isPulsed(isPulsed)
{}

~LArCaliWave()

virtual LArCaliWave::~LArCaliWave ( )

virtualdefault

Member Function Documentation

addAccumulatedEvent()

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

inherited

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 )

inherited

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

getDAC()

int LArCaliWave::getDAC ( ) const

inline

DAC value.

Definition at line 156 of file LArCaliWave.h.

{
    return m_DAC;
}

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

inlineinherited

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

inlineinherited

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

getIsPulsedInt()

int LArCaliWave::getIsPulsedInt ( ) const

inline

isPulsed value

Definition at line 163 of file LArCaliWave.h.

{ 
  return m_isPulsed;
}

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

inlineinherited

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

inlineinherited

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_DAC

int LArCaliWave::m_DAC

protected

Definition at line 86 of file LArCaliWave.h.

m_dt

double LArWave::m_dt

protectedinherited

Definition at line 111 of file LArWave.h.

m_errors

std::vector<double> LArWaveCumul::m_errors

protectedinherited

Definition at line 89 of file LArWaveCumul.h.

m_flag

unsigned LArWave::m_flag

protectedinherited

Definition at line 113 of file LArWave.h.

m_isPulsed

int LArCaliWave::m_isPulsed

protected

Definition at line 87 of file LArCaliWave.h.

m_triggers

std::vector<int> LArWaveCumul::m_triggers

protectedinherited

Definition at line 90 of file LArWaveCumul.h.

---

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

• LArCaliWave.h