#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 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.

147 : m_dt(0), m_flag(0)

148 {}

◆ 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.

50 { return m_dt; }

◆ getFlag()

unsigned LArWave::getFlag ( ) const

inline

flag: ...

Definition at line 178 of file LArWave.h.

179 { return m_flag ; }

◆ getIndex()

unsigned LArWave::getIndex ( double aTime ) const

index for a time value

Definition at line 146 of file LArWave.cxx.

147 { 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.

53 {return m_amplitudes[i];}

◆ getSize()

size_t LArWave::getSize ( ) const

inline

number of time samples

Definition at line 62 of file LArWave.h.

62 { return m_amplitudes.size() ; }

◆ getTime()

double LArWave::getTime ( const unsigned i ) const

inline

time

Definition at line 172 of file LArWave.h.

173 { return i * m_dt ; }

◆ getWave()

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

inline

Wave parameters.

Definition at line 167 of file LArWave.h.

168 { return m_amplitudes ; }

◆ isEmpty()

bool LArWave::isEmpty ( ) const

inline

is LArWave uninitialized?

Definition at line 183 of file LArWave.h.

184 { 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.

78 {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.

56 {return m_amplitudes[i];}

◆ setDt()

void LArWave::setDt ( const double dt )

inline

set the delta time

Definition at line 193 of file LArWave.h.

194 { m_dt = dt ; }

◆ setFlag()

void LArWave::setFlag ( const unsigned flag )

inline

set flag

Definition at line 199 of file LArWave.h.

200 { 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.

87 { 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.

189 { 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:

Public Member Functions

accumulation alegebra operations

Detailed Description

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ LArWave() [1/4]

◆ LArWave() [2/4]

◆ LArWave() [3/4]

◆ ~LArWave()

◆ LArWave() [4/4]

Member Function Documentation

◆ getDt()

◆ getFlag()

◆ getIndex()

◆ getSample()

◆ getSize()

◆ getTime()

◆ getWave()

◆ isEmpty()

◆ operator%()

◆ operator*() [1/2]

◆ operator*() [2/2]

◆ operator*=()

◆ operator+() [1/2]

◆ operator+() [2/2]

◆ operator+=()

◆ operator-()

◆ operator-=()

◆ operator/()

◆ operator[]() [1/2]

◆ operator[]() [2/2]

◆ setDt()

◆ setFlag()

◆ setSample()

◆ setSize()

Member Data Documentation

◆ m_amplitudes

◆ m_dt

◆ m_flag