ZdcSignalSinc Class Reference

#include <ZdcSignalSinc.h>

Collaboration diagram for ZdcSignalSinc:

Public Types
enum	Status {   e_OK =0, e_noData, e_wrongFrac, e_Overflow,   e_wrongSignal, e_noSignal, e_localMinimum }

Public Member Functions
	ZdcSignalSinc (int)
	~ZdcSignalSinc ()
	ZdcSignalSinc (const ZdcSignalSinc &)=delete
ZdcSignalSinc &	operator= (const ZdcSignalSinc &)=delete
int	process (double *, double gain=1., double ped=0., double frac=1., bool corr=true)
int	getError () const
int	getWarning () const
double	getTime () const
double	getAmp () const

Private Member Functions
double	waveform (double t)
double	findpeak (int)
double	findpeak ()
double	fraction (double, double)

Private Attributes
const int	m_n_Slices
const double	m_AmpThresh
const double	m_tClock
const double	m_Pi
double	m_Time
double	m_Amp
int	m_Err
int	m_Warn
double *	m_buf
int	m_p
int	m_np
double	m_tim [3] {}
double	m_wfm [3] {}
double	m_dt
int	m_AAAA

Detailed Description

Definition at line 14 of file ZdcSignalSinc.h.

Member Enumeration Documentation

Status

enum ZdcSignalSinc::Status

Enumerator
e_OK
e_noData
e_wrongFrac
e_Overflow
e_wrongSignal
e_noSignal
e_localMinimum

Definition at line 16 of file ZdcSignalSinc.h.

              {e_OK=0, e_noData, e_wrongFrac, e_Overflow,
           e_wrongSignal, e_noSignal, e_localMinimum};

Constructor & Destructor Documentation

ZdcSignalSinc() [1/2]

ZdcSignalSinc::ZdcSignalSinc

(

int

nSlc

)

Definition at line 11 of file ZdcSignalSinc.cxx.

  : m_n_Slices(nSlc)
    , m_AmpThresh(5.)
    , m_tClock      (25.)
    , m_Pi(4*atan(double(1.)))
{
  m_buf = new double[m_n_Slices];
  m_Time   =  0.;
  m_Amp    =  0.;
  m_Err    =  e_noData;
  m_Warn   = 0;
 
  m_p      = 0;
  m_np     = 0;
 
  m_wfm[0] = 0. ;
  m_wfm[1] = 0. ;
  m_wfm[2] = 0. ;
 
  m_tim[0] = 0. ;
  m_tim[1] = 0. ;
  m_tim[2] = 0. ;
 
  m_dt = 0.;
 
  m_AAAA = 0;
 
}

~ZdcSignalSinc()

ZdcSignalSinc::~ZdcSignalSinc

(

)

Definition at line 40 of file ZdcSignalSinc.cxx.

{ delete[] m_buf;}

ZdcSignalSinc() [2/2]

ZdcSignalSinc::ZdcSignalSinc ( const ZdcSignalSinc &  )

delete

Member Function Documentation

findpeak() [1/2]

double ZdcSignalSinc::findpeak ( )

private

Definition at line 201 of file ZdcSignalSinc.cxx.

                               {
  if (m_dt<0.001) return m_tim[1];
 
  if (m_wfm[0]<m_wfm[1]) {
    if (m_wfm[2]<=m_wfm[1]) {
      m_dt /=2.;
      double t1 = m_tim[1]-m_dt;
      double f1 = waveform(t1);
      if (f1>m_wfm[1]) {
    m_tim[2] = m_tim[1]; m_wfm[2] = m_wfm[1];
    m_tim[1] = t1;      m_wfm[1] = f1;
      } else {
    double t2 = m_tim[1]+m_dt;
    double f2 = waveform(t2);
    if (f2>m_wfm[1]) {
      m_tim[0] = m_tim[1]; m_wfm[0] = m_wfm[1];
      m_tim[1] = t2;      m_wfm[1] = f2;
    } else {
      m_tim[0] = t1;      m_wfm[0] = f1;
      m_tim[2] = t2;      m_wfm[2] = f2;
    }
      }
    } else {
      m_tim[0] = m_tim[1]; m_wfm[0] = m_wfm[1];
      m_tim[1] = m_tim[2]; m_wfm[1] = m_wfm[2];
      m_tim[2] += m_dt;     m_wfm[2] = waveform(m_tim[2]);
    }
  } else {
    if (m_wfm[2]<=m_wfm[1]) {
      m_tim[2] = m_tim[1]; m_wfm[2] = m_wfm[1];
      m_tim[1] = m_tim[0]; m_wfm[1] = m_wfm[0];
      m_tim[0] -= m_dt;    m_wfm[0] = waveform(m_tim[0]);
      return findpeak();
    } else {
      m_Err = e_localMinimum;
      return 0.;
    }
  }
  return findpeak();
}

findpeak() [2/2]

double ZdcSignalSinc::findpeak ( int  im )

private

Definition at line 192 of file ZdcSignalSinc.cxx.

                                     {
 
  m_tim[0] = im; m_tim[1]=im+0.5; m_tim[2]=im+1.;
  for (int i=0;i<3;i++) m_wfm[i]=waveform(m_tim[i]);
  m_dt = 0.5;
  double t = findpeak();
  return t;
}

fraction()

double ZdcSignalSinc::fraction ( double  frac, double  tpeak  )

private

Definition at line 166 of file ZdcSignalSinc.cxx.

                                                        {
  m_tim[0] = 0.;    m_wfm[0] = waveform(m_tim[0]);
  m_tim[1] = tpeak; m_wfm[1] = waveform(m_tim[1]);
  m_dt     = tpeak/2.;
  double thr = frac*m_wfm[1];
  while(m_wfm[0]>thr) {
    m_tim[0] -= 0.1;
    m_wfm[0] = waveform(m_tim[0]);
    if (m_tim[0]<-2.) {
      m_tim[0] = 0.;
      m_wfm[0] = waveform(m_tim[0]);
      m_Err = e_wrongSignal;
      break;
    }
  }
 
  while (m_dt>0.001) {
    double t = m_tim[0] + m_dt;
    double f = waveform(t);
    if (f>thr) {m_tim[1] = t; m_wfm[1] = f;}
    else       {m_tim[0] = t; m_wfm[0] = f;}
    m_dt /= 2.;
  }
  return (m_tim[0]+m_tim[1])/2.;
}

getAmp()

double ZdcSignalSinc::getAmp

(

)

const

Definition at line 154 of file ZdcSignalSinc.cxx.

{return m_Amp;}

getError()

int ZdcSignalSinc::getError

(

)

const

Definition at line 151 of file ZdcSignalSinc.cxx.

{return m_Err;}

getTime()

double ZdcSignalSinc::getTime

(

)

const

Definition at line 153 of file ZdcSignalSinc.cxx.

{return m_Time;}

getWarning()

int ZdcSignalSinc::getWarning

(

)

const

Definition at line 152 of file ZdcSignalSinc.cxx.

{return m_Warn;}

operator=()

ZdcSignalSinc& ZdcSignalSinc::operator= ( const ZdcSignalSinc &  )

delete

process()

int ZdcSignalSinc::process	(	double *	buf,
		double	gain = 1.,
		double	ped = 0.,
		double	frac = 1.,
		bool	corr = true
	)

Definition at line 42 of file ZdcSignalSinc.cxx.

                                       {
 
  m_Time = 0.;
  m_Amp  = 0.;
  m_Err  = e_OK;
  m_Warn = 0;
 
  if (ped==0) ped=buf[0];
  for (double *p=m_buf;p<m_buf+m_n_Slices;p++) {
    double a = *(buf++);
    //std::cout << "ZDC+++ A" << a << std::endl;
    if (a>1015) {
      m_Err = e_Overflow;
      return  m_Err;
    }
    *p = a-ped;
  }
  if (frac<=0.) {
    m_Err = e_wrongFrac;
    return m_Err;
  }
 
  int imax = -1;
  if (m_buf[0         ]-m_buf[1         ]>m_AmpThresh) {
    imax = 0;
    m_Warn += 1;
  }
 
  if (m_buf[m_n_Slices-1]-m_buf[m_n_Slices-2]>m_AmpThresh) {
    if (imax<0) {
      imax = m_n_Slices-1;
      m_Warn += 2;
    } else {
      m_Err = e_wrongSignal;
      return m_Err;
    }
  }
 
  for (int i=1;i<m_n_Slices-1;i++) {
    double da1 = m_buf[i]-m_buf[i-1];
    double da2 = m_buf[i]-m_buf[i+1];
    if (da1>0. && da2>=0.) {
      if (da1>m_AmpThresh || da2>m_AmpThresh) {
    if (imax<0) imax = i;
    else {
      m_Err = e_wrongSignal;
      return m_Err;
    }
      }
    }
  }
 
  if (imax<0) {
    m_Err = e_noSignal;
    return m_Err;
  }
 
  if (imax==1) m_Warn += 4;
 
  if (imax>=2) {
    //       printf("%d %f %f\n",imax,m_buf[imax-2],m_AmpThresh);
    if (m_buf[imax-2]>m_AmpThresh)  m_p = imax-2;
    else                            m_p = imax-1;
  } else m_p = imax;
 
  m_np = m_n_Slices-m_p;
 
  double t_peak = findpeak(imax-m_p);
  if (m_Err) return m_Err;
 
  if (frac>=1.) m_Time = t_peak;
  else          m_Time = fraction(frac,t_peak);
  m_Amp = waveform(t_peak);
 
  //std::cout << "ZDC+++ m_AMP " << m_Amp << std::endl;
 
  /*
  printf("aap: %d %d %d  ",imax, m_p, m_np);
  for (int i=0;i<m_n_Slices;i++) printf("%6.1f ",m_buf[i]);
  printf("  %6.1f %6.1f %6.1f\n",t_peak*m_tClock,m_Time*m_tClock,m_Amp);
  */
 
  if (gain==1.) {
    if (m_Amp<200.) m_Warn +=  8;
    if (m_Amp>900.) m_Warn += 16;
  } else {
    if (m_Amp< 50.) m_Warn +=  8;
    if (m_Amp>250.) m_Warn += 16;
  }
  m_Amp *= gain;
  m_Time = (m_Time+m_p)*m_tClock;
  if (corr) {
    int off = int(10.+m_Time/m_tClock)-10;
    double t = m_Time - off*m_tClock;
 
    if      (t< 5.) t =   0.056836274872111      + 1.532763686481279  *t
              +   0.229808343735056 *t*t - 0.0365636647119192 *t*t*t;
    else if (t<16.) t =  -1.28974955223497       + 2.653578699334604  *t
              -   0.1371240146140209*t*t + 0.00281211806384422*t*t*t;
    else            t = -42.18688740322650       + 8.61752055701946   *t
              -   0.4259239806065329*t*t + 0.00753849507486617*t*t*t;
 
    m_Time = off*m_tClock + t;
  }
 
  return m_Err;
}

waveform()

double ZdcSignalSinc::waveform ( double  t )

private

Definition at line 156 of file ZdcSignalSinc.cxx.

                                       {
  double f = 0.;
  for (int i=0; i<m_np; i++) {
    double x = m_Pi*(t-i);
    if (x) f += m_buf[m_p+i]*sin(x)/x;
    else   f += m_buf[m_p+i];
  }
  return f;
}

Member Data Documentation

m_AAAA

int ZdcSignalSinc::m_AAAA

private

Definition at line 54 of file ZdcSignalSinc.h.

m_Amp

double ZdcSignalSinc::m_Amp

private

Definition at line 38 of file ZdcSignalSinc.h.

m_AmpThresh

const double ZdcSignalSinc::m_AmpThresh

private

Definition at line 34 of file ZdcSignalSinc.h.

m_buf

double* ZdcSignalSinc::m_buf

private

Definition at line 44 of file ZdcSignalSinc.h.

m_dt

double ZdcSignalSinc::m_dt

private

Definition at line 53 of file ZdcSignalSinc.h.

m_Err

int ZdcSignalSinc::m_Err

private

Definition at line 39 of file ZdcSignalSinc.h.

m_n_Slices

const int ZdcSignalSinc::m_n_Slices

private

Definition at line 33 of file ZdcSignalSinc.h.

m_np

int ZdcSignalSinc::m_np

private

Definition at line 47 of file ZdcSignalSinc.h.

m_p

int ZdcSignalSinc::m_p

private

Definition at line 46 of file ZdcSignalSinc.h.

m_Pi

const double ZdcSignalSinc::m_Pi

private

Definition at line 36 of file ZdcSignalSinc.h.

m_tClock

const double ZdcSignalSinc::m_tClock

private

Definition at line 35 of file ZdcSignalSinc.h.

m_tim

double ZdcSignalSinc::m_tim[3] {}

private

Definition at line 53 of file ZdcSignalSinc.h.

m_Time

double ZdcSignalSinc::m_Time

private

Definition at line 37 of file ZdcSignalSinc.h.

m_Warn

int ZdcSignalSinc::m_Warn

private

Definition at line 40 of file ZdcSignalSinc.h.

m_wfm

double ZdcSignalSinc::m_wfm[3] {}

private

Definition at line 53 of file ZdcSignalSinc.h.

---

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

• ZdcSignalSinc.h
• ZdcSignalSinc.cxx