EclConfiguration::signalsample_t Struct Reference

a struct for a signal sample More...

#include <EclConfiguration.h>

Public Member Functions
void	InitSample (const float *MP, double unitscale=-1)
void	InitSample (const double *MPd, double unitscale=-1)
double	Accumulate (const double a, const double t0, double *s) const

Public Attributes
double	m_sumscale
	energy deposit in fitting window scale factor
double	m_ft [m_nl *m_ns]
	Simulated signal shape.

Detailed Description

a struct for a signal sample

Definition at line 63 of file EclConfiguration.h.

Member Function Documentation

Accumulate()

double Accumulate	(	const double	a,
		const double	t0,
		double *	s
	)		const

Parameters

[in]	a	Signal amplitude
[in]	t0	Signal offset
[out]	s	Output array with added signal

Returns

Energy deposition in ADC units

Definition at line 79 of file EclConfiguration.cc.

{
  // input parameters
  // a -- signal amplitude
  // t -- signal offset
  // output parameter
  // s -- output array with added signal
  const double itick = getRF() / s_clock;  // reciprocal to avoid division in usec^-1 (has to be evaluated at compile time)
  const double  tlen = m_nl - 1.0 / m_ns;   // length of the sampled signal in ADC clocks units
  const double  tmax = m_tmin + m_nsmp - 1; // upper range of the fit region
 
  double t = t0 * itick; // rescale time in usec to ADC clocks
  double x0 = t, x1 = t + tlen;
 
  if (x0 > tmax) return 0; // signal starts after the upper range of output array -- do nothing
  if (x0 < m_tmin) {
    if (x1 < m_tmin) return 0; // signal ends before lower range of output array -- do nothing
    x0 = m_tmin; // adjust signal with range of output array
  }
 
  int imax = m_nsmp; // length of sampled signal is long enough so
  // the last touched element is the last element
  // of the output array
  if (x1 < tmax) { // if initial time is too early we need to adjust
    // the last touched element of output array to avoid
    // out-of-bound situation in m_ft
    imax = x1 - m_tmin; // imax is always positive so floor can be
    // replace by simple typecast
    imax += 1; // like s.end()
  }
 
  double imind = ceil(x0 - m_tmin); // store result in double to avoid int->double conversion below
  // the ceil function today at modern CPUs is surprisingly fast (before it was horribly slow)
  int imin = imind; // starting point to fill output array
  double w = ((m_tmin - t) + imind) * double(m_ns);
  int jmin = w; // starting point in sampled signal array
  w -= jmin;
 
  // use linear interpolation between samples. Since signal samples
  // are aligned with output samples only two weights are need to
  // calculate to fill output array
  const double w1 = a * w, w0 = a - w1;
  double sum = 0;
  for (int i = imin, j = jmin; i < imax; i++, j += m_ns) {
    double amp = w0 * m_ft[j] + w1 * m_ft[j + 1];
    s[i] += amp;
    sum  += amp;
  }
  return sum * m_sumscale;
}

InitSample() [1/2]

void InitSample	(	const double *	MPd,
		double	unitscale = -1
	)

Parameters

MPd	Double array of waveform parameters
unitscale	Normalization of template waveform

Definition at line 53 of file EclConfiguration.cc.

{
  const int N = m_ns * m_nl;
  vector<double> MP(MPd, MPd + 10);
  ShaperDSP_t dsp(MP, u);
  dsp.settimestride(m_step / m_ns);
  dsp.fillarray(0.0, N, m_ft);
 
  double sum = 0;
  for (int i = 0; i < N; i++) sum += m_ft[i];
  m_sumscale = m_ns / sum;
}

InitSample() [2/2]

void InitSample	(	const float *	MP,
		double	unitscale = -1
	)

Parameters

MP	Float array of waveform parameters
unitscale	Normalization of template waveform

Definition at line 67 of file EclConfiguration.cc.

{
  double MPd[10];
  for (int i = 0; i < 10; i++) MPd[i] = MP[i];
  InitSample(MPd, u);
}

Member Data Documentation

m_ft

double m_ft[m_nl *m_ns]

Simulated signal shape.

Definition at line 82 of file EclConfiguration.h.

m_sumscale

double m_sumscale

energy deposit in fitting window scale factor

Definition at line 81 of file EclConfiguration.h.

---

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

• ecl/digitization/include/EclConfiguration.h
• ecl/digitization/src/EclConfiguration.cc