#include <LArAutoCorr.h>

Collaboration diagram for LArAutoCorr:

Public Member Functions
	LArAutoCorr ()

	~LArAutoCorr ()

void	correl_zero ()

void	set_min (const short min)

void	set_max (const short max)

int	get_nentries () const

double	get_mean () const

double	get_rms () const

const short &	get_min () const

const short &	get_max () const

void	add (const std::vector< short > &samples, size_t maxnsamples)

const std::vector< double > &	get_cov (int m_normalize, int m_phys)

	LArAutoCorr ()

	~LArAutoCorr ()

void	correl_zero ()

void	set_min (const short min)

void	set_max (const short max)

int	get_nentries () const

double	get_mean () const

double	get_rms () const

const short &	get_min () const

const short &	get_max () const

void	add (const std::vector< short > &samples, size_t maxnsamples)

const std::vector< double > &	get_cov (int normalize, int phys)

Private Attributes
short	m_min

short	m_max

std::vector< double >	m_sum

std::vector< double >	m_matrix

std::vector< double >	m_cov

std::vector< double >	m_cov_temp

int	m_nped

Detailed Description

Definition at line 26 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

Constructor & Destructor Documentation

◆ LArAutoCorr() [1/2]

LArAutoCorr::LArAutoCorr ( )

inline

Definition at line 64 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ ~LArAutoCorr() [1/2]

LArAutoCorr::~LArAutoCorr ( )

inline

Definition at line 72 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ LArAutoCorr() [2/2]

LArAutoCorr::LArAutoCorr ( )

inline

Definition at line 64 of file LArRawEvent/LArRawEvent/LArAutoCorr.h.

◆ ~LArAutoCorr() [2/2]

LArAutoCorr::~LArAutoCorr ( )

inline

Definition at line 72 of file LArRawEvent/LArRawEvent/LArAutoCorr.h.

Member Function Documentation

◆ add() [1/2]

void LArAutoCorr::add	(	const std::vector< short > &	samples,
		size_t	maxnsamples
	)

Definition at line 222 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   unsigned int nsamples = std::min(samples.size(),maxnsamples);
   int k=0;
  
   if(m_sum.size()<nsamples)
     {
     
       m_cov.resize(nsamples - 1);
       m_cov_temp.resize((nsamples*(nsamples+1))/2);
       m_sum.resize(nsamples);
       m_matrix.resize((nsamples*(nsamples+1))/2);
     }
  
   for(unsigned int i=0; i<nsamples; i++)
     {
       for(unsigned int j=i; j<nsamples; j++,k++)
     m_matrix[k] += ((double)(samples[j]*samples[i]));
   
       m_sum[i] += ((double) samples[i]);
     }
   m_nped++;
 }

◆ add() [2/2]

void LArAutoCorr::add	(	const std::vector< short > &	samples,
		size_t	maxnsamples
	)

◆ correl_zero() [1/2]

void LArAutoCorr::correl_zero ( )

Definition at line 248 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   int j =0;
   int nsamples = m_sum.size();
   for(int l=0; l<nsamples; l++)
     {
       for(int k=l; k<nsamples; k++,j++)
     m_matrix[j] = 0;
       m_sum[l]=0;
     }
   m_nped=0;
 }

◆ correl_zero() [2/2]

void LArAutoCorr::correl_zero ( )

◆ get_cov() [1/2]

const std::vector< double > & LArAutoCorr::get_cov	(	int	m_normalize,
		int	m_phys
	)

Definition at line 149 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 { 
   int nsamples = m_sum.size();
   int k =0;
   if (m_nped==0) return m_cov;
   for(int i=0;i<nsamples;i++){
       for(int j=i;j<nsamples;j++,k++)
     m_cov_temp[k]= (double)((m_matrix[k]/(m_nped)) - ((m_sum[i]*m_sum[j])/(m_nped*m_nped)));
   }
  
   //MGV normalize covariance elements if required
   if (m_normalize == 1 && m_phys == 0 ){
     k=0;
     int si,sj;
     si=0;
     for(int i=0;i<nsamples;i++){
       sj=si;
       for(int j=i;j<nsamples;j++,k++) {
     if (m_sum[i]!=.0 && m_sum[j]!=.0) 
       m_cov_temp[k]= m_cov_temp[k]/(double)sqrt((m_matrix[si]/(m_nped)-(m_sum[i]*m_sum[i])/(m_nped*m_nped))*(m_matrix[sj]/(m_nped)-(m_sum[j]*m_sum[j])/(m_nped*m_nped)));
     sj+=nsamples-j;
       }
       si+=nsamples-i;
     }
   }
   //MGV
  
   if (m_phys==0) {
     double sum;
     int s;
     for(int diag =1;diag<nsamples;diag++)
       {
         sum =0;
         s = 0;
         for(int i=0; i < nsamples-diag;i++)
       { 
         sum += m_cov_temp[s + diag];
         s += nsamples - i;
       }
         sum = sum/(nsamples-diag);
         m_cov[diag-1]= sum;
       }
     return m_cov;
   }
   else
     return m_cov_temp;
 }

◆ get_cov() [2/2]

const std::vector<double>& LArAutoCorr::get_cov	(	int	normalize,
		int	phys
	)

◆ get_max() [1/2]

const short & LArAutoCorr::get_max ( ) const

Definition at line 126 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   return m_max;
 }

◆ get_max() [2/2]

const short& LArAutoCorr::get_max ( ) const

◆ get_mean() [1/2]

double LArAutoCorr::get_mean ( ) const

Definition at line 134 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   double mean = 0;
  
  int nsamples = m_sum.size();
   for(int i=0; i<nsamples; i++)
     mean += m_sum[i];
     mean /= ((double)(nsamples*m_nped));
  
   return mean;
 }

◆ get_mean() [2/2]

double LArAutoCorr::get_mean ( ) const

◆ get_min() [1/2]

const short& LArAutoCorr::get_min ( ) const

◆ get_min() [2/2]

const short & LArAutoCorr::get_min ( ) const

Definition at line 119 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   return m_min;
 }

◆ get_nentries() [1/2]

int LArAutoCorr::get_nentries ( ) const

◆ get_nentries() [2/2]

int LArAutoCorr::get_nentries ( ) const

Definition at line 111 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   return m_nped;
 } 

◆ get_rms() [1/2]

double LArAutoCorr::get_rms ( ) const

Definition at line 199 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   double x=0, y=0;
   int k = 0;
   int nsamples = m_sum.size();
  
   for(int i=0; i<nsamples; i++)
     {
       x += m_sum[i];
       y += m_matrix[k];
       k += nsamples - i;// Index of diagonal element
     }
   
   x/= (double) (nsamples*m_nped);
   y/=(double) (nsamples*m_nped);
   
   double noise = sqrt(y- x*x);
   
   return noise;
 }

◆ get_rms() [2/2]

double LArAutoCorr::get_rms ( ) const

◆ set_max() [1/2]

void LArAutoCorr::set_max ( const short max )

◆ set_max() [2/2]

void LArAutoCorr::set_max ( const short max )

Definition at line 104 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   m_max = max;
 } 

◆ set_min() [1/2]

void LArAutoCorr::set_min ( const short min )

◆ set_min() [2/2]

void LArAutoCorr::set_min ( const short min )

Definition at line 97 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

 {
   m_min = min;
 }

Member Data Documentation

◆ m_cov

std::vector< double > LArAutoCorr::m_cov

private

Definition at line 55 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ m_cov_temp

std::vector< double > LArAutoCorr::m_cov_temp

private

Definition at line 56 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ m_matrix

std::vector< double > LArAutoCorr::m_matrix

private

Definition at line 52 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ m_max

short LArAutoCorr::m_max

private

Definition at line 46 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ m_min

short LArAutoCorr::m_min

private

Definition at line 43 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ m_nped

int LArAutoCorr::m_nped

private

Definition at line 59 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

◆ m_sum

std::vector< double > LArAutoCorr::m_sum

private

Definition at line 49 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ LArAutoCorr() [1/2]

◆ ~LArAutoCorr() [1/2]

◆ LArAutoCorr() [2/2]

◆ ~LArAutoCorr() [2/2]

Member Function Documentation

◆ add() [1/2]

◆ add() [2/2]

◆ correl_zero() [1/2]

◆ correl_zero() [2/2]

◆ get_cov() [1/2]

◆ get_cov() [2/2]

◆ get_max() [1/2]

◆ get_max() [2/2]

◆ get_mean() [1/2]

◆ get_mean() [2/2]

◆ get_min() [1/2]

◆ get_min() [2/2]

◆ get_nentries() [1/2]

◆ get_nentries() [2/2]

◆ get_rms() [1/2]

◆ get_rms() [2/2]

◆ set_max() [1/2]

◆ set_max() [2/2]

◆ set_min() [1/2]

◆ set_min() [2/2]

Member Data Documentation

◆ m_cov

◆ m_cov_temp

◆ m_matrix

◆ m_max

◆ m_min

◆ m_nped

◆ m_sum