#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 52 of file LArCalibUtils/LArCalibUtils/LArAutoCorr.h.

    {
      m_min = -1;
      m_max = -1;
      m_nped = 0;
    }

◆ ~LArAutoCorr() [1/2]

LArAutoCorr::~LArAutoCorr ( )

inline

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

60{}

◆ LArAutoCorr() [2/2]

LArAutoCorr::LArAutoCorr ( )

inline

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

    {
      m_min = -1;
      m_max = -1;
      m_nped = 0;
    }

◆ ~LArAutoCorr() [2/2]

LArAutoCorr::~LArAutoCorr ( )

inline

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

60{}

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

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

{
  return m_min;
}

◆ get_min() [2/2]

const short & LArAutoCorr::get_min ( ) const

◆ get_nentries() [1/2]

int LArAutoCorr::get_nentries ( ) const

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

{
  return m_nped;
} 

◆ get_nentries() [2/2]

int LArAutoCorr::get_nentries ( ) const

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

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

{
  m_max = max;
} 

◆ set_max() [2/2]

void LArAutoCorr::set_max ( const short max )

◆ set_min() [1/2]

void LArAutoCorr::set_min ( const short min )

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

{
  m_min = min;
}

◆ set_min() [2/2]

void LArAutoCorr::set_min ( const short min )

Member Data Documentation

◆ m_cov

std::vector< double > LArAutoCorr::m_cov

private

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

◆ m_cov_temp

std::vector< double > LArAutoCorr::m_cov_temp

private

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

◆ m_matrix

std::vector< double > LArAutoCorr::m_matrix

private

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

◆ m_max

short LArAutoCorr::m_max

private

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

◆ m_min

short LArAutoCorr::m_min

private

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

◆ m_nped

int LArAutoCorr::m_nped

private

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

◆ m_sum

std::vector< double > LArAutoCorr::m_sum

private

Definition at line 37 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