ShaperDSP_t Class Reference

Class include function that calculate electronic response from energy deposit. More...

#include <shaperdsp.h>

Classes
struct	shaperdspshift_t
	struct for a shift of the shaper dsp More...
struct	sv123shift_t
	struct to encapsulate the electronic response from energy deposit More...

Public Member Functions
	ShaperDSP_t ()
	class constructor
	ShaperDSP_t (const std::vector< double > &s)
	calculate derivative of the response function
	ShaperDSP_t (const std::vector< double > &s, double u)
	class constructor
void	init (const std::vector< double > &s, double u=-1)
	Initialization of the parameters response function.
double	operator() (double) const
	wrapper of the function
double	operator() (const double *, double *)
	TF1 ROOT interface.
void	settimestride (double)
	set grid step for function calculation
void	setseedoffset (double)
	set timeoffset
void	settimeseed (double)
	set initial time
void	nextseed ()
	substruct toffset to tzero
void	fillvector (std::vector< double > &) const
	fill vector with response function values and its derivative
void	fillvector (std::vector< dd_t > &) const
	fill vector with response function values and its derivative
void	fillvector (double, std::vector< double > &) const
	fill vector with response function values and its derivative
void	fillvector (double, std::vector< dd_t > &) const
	fill vector with response function values and its derivative
void	fillarray (int, double *) const
	fill array for amplitude and time calculation
void	fillarray (int, dd_t *) const
	fill array for amplitude and time calculation
void	fillarray (double, int, double *) const
	fill array for amplitude and time calculation
void	fillarray (double, int, dd_t *) const
	fill array for amplitude and time calculation

Private Member Functions
void	Sv123_init (double t01, double tb1, double t02, double tb2, double td1, double ts1)
	calculate some values for Sv123 function
void	init (const double *, double)
	calculate some values for Sv123 function
void	init (const double *)
	calculate some values for Sv123 function
double	Sv123 (const sv123shift_t &) const
	calculate Sv123 function
dd_t	ddSv123 (const sv123shift_t &) const
	calculate derivative of the Sv123 function
double	Sv123_filtered (const sv123shift_t &) const
	Numerical calculation of the time convolution.
dd_t	ddSv123_filtered (const sv123shift_t &) const
	This is derivative of the confolution.
double	ShaperDSP (const shaperdspshift_t &) const
	calculate response function
dd_t	ddShaperDSP (const shaperdspshift_t &) const
	calculate derivative of the response function

Private Attributes
double	_cs0 = realNaN
	linear coefficient before sin of the first Bessel stage
double	_cc0 = realNaN
	linear coefficient before cos of the first Bessel stage
double	_cs1 = realNaN
	linear coefficient before sin of the second Bessel stage
double	_cc1 = realNaN
	linear coefficient before cos of the second Bessel stage
double	_ces = realNaN
	linear coefficient before first part of tail section
double	_ced = realNaN
	linear coefficient before second part of tail section
double	_dw0 = realNaN
	circular frequency of the first Bessel stage
double	_dw1 = realNaN
	circular frequency of the second Bessel stage
double	_dks0 = realNaN
	decrement of the first Bessel stage
double	_dks1 = realNaN
	decrement of the second Bessel stage
double	_ds = realNaN
	inverse scintillation decay time
double	_dd = realNaN
	inverse time of the differential stage
double	_dt0 = realNaN
	coefficient for first exponent factor
double	_dt1 = realNaN
	coefficient for second exponent factor
double	_toff = realNaN
	time offset
double	_w0 = realNaN
	weight coefficient at sv123(t) = (1-a)
double	_w1 = realNaN
	weight coefficient at sv123(t+_filterdt) +sv123(t-_filterdt) = a/2
double	_ccc = realNaN
	exponent factor for tail part of the signal
sv123shift_t	_tp
	_filterdt
sv123shift_t	_tm
shaperdspshift_t	_tstride
	time step of the grid for response function calculation
shaperdspshift_t	_toffset
	time offset
shaperdspshift_t	_tzero
	initial time

Static Private Attributes
static const double	_defs [] = {0.5, 0.6483, 0.4017, 0.3741, 0.8494, 0.00144547, 4.7071, 0.8156, 0.5556, 0.2752}
	parameters of the response function that use as default
static constexpr double	_filterdt = 0.2
	time shift that include in response function for numerical calculation time convolutions.

Detailed Description

Class include function that calculate electronic response from energy deposit.

Definition at line 26 of file shaperdsp.h.

Constructor & Destructor Documentation

ShaperDSP_t() [1/3]

ShaperDSP_t ( )

inline

class constructor

Definition at line 156 of file shaperdsp.h.

{ init(_defs); }

ShaperDSP_t() [2/3]

ShaperDSP_t ( const std::vector< double > & s )

inlineexplicit

calculate derivative of the response function

Definition at line 158 of file shaperdsp.h.

{ init(s); }

ShaperDSP_t() [3/3]

ShaperDSP_t ( const std::vector< double > & s, double u )

inline

class constructor

Definition at line 160 of file shaperdsp.h.

{ init(s, u); }

~ShaperDSP_t()

~ShaperDSP_t ( )

inline

Definition at line 162 of file shaperdsp.h.

{}

Member Function Documentation

ddShaperDSP()

dd_t ddShaperDSP ( const shaperdspshift_t & t0 ) const

private

calculate derivative of the response function

Definition at line 290 of file shaperdsp.cc.

{
  const sv123shift_t& t00 = static_cast<const ShaperDSP_t::sv123shift_t&>(t0);
  dd_t f = ddSv123_filtered(t00);
  if (t0.t > 0) {
    double z = t0.t * _dt1;
    double z2 = z * z;
    double odd = 1 + z2 * (1 / 6. + z2 * (1 / 120.));
    double evn = 1 + z2 * (1 / 2. + z2 * (1 / 24.));
    double texp = evn + z * odd;
    double u = t0.et0 * (1 - t0.et1 * texp);
    f.first -= _ccc * u;
 
    double up = -_dt0 * u + ((_dt1 * t0.et0) * (t0.et1 * z2)) * ((z2 * z) * (1 / 120.));
    f.second -= _ccc * up;
  }
  return f;
}

ddSv123()

dd_t ddSv123 ( const sv123shift_t & c ) const

private

calculate derivative of the Sv123 function

Definition at line 229 of file shaperdsp.cc.

{
  double f0 = (_cs0 * c.s0 + _cc0 * c.c0) * c.e0;
  double f1 = (_cs1 * c.s1 + _cc1 * c.c1) * c.e1;
  double fs =  _ces * c.es;
  double fd =  _ced * c.ed;
 
  double f0p = (_dw0 * c.e0) * (_cs0 * c.c0 - _cc0 * c.s0) - (_dks0 * f0);
  double f1p = (_dw1 * c.e1) * (_cs1 * c.c1 - _cc1 * c.s1) - (_dks1 * f1);
  double fsp =  _ds * fs;
  double fdp =  _dd * fd;
 
  double f  = (f0  + f1) + (fd  - fs);
  double fp = (f0p + f1p) - (fdp - fsp);
  return dd_t(f, fp);
}

ddSv123_filtered()

dd_t ddSv123_filtered ( const sv123shift_t & t0 ) const

private

This is derivative of the confolution.

Definition at line 260 of file shaperdsp.cc.

{
  if (t0.validshift(_tp)) {
    dd_t ft0p = ddSv123(t0 + _tp), ft0 = dd_t(0, 0), ft0m = dd_t(0, 0);
    if (t0.t > 0) {
      ft0 = ddSv123(t0);
      if (t0.validshift(_tm))
        ft0m = ddSv123(t0 + _tm);
    }
    return _w0 * ft0 + _w1 * (ft0p + ft0m);
  }
  return dd_t(0, 0);
}

fillarray() [1/4]

void fillarray	(	double	t,
		int	n,
		dd_t *	s ) const

fill array for amplitude and time calculation

Definition at line 409 of file shaperdsp.cc.

{
  shaperdspshift_t t0 = shaperdspshift_t(t  - _toff, *this);
 
  for (int i = 0; i < n; i++, t0 += _tstride) {
    s[i] = ddShaperDSP(t0);
  }
}

fillarray() [2/4]

void fillarray	(	double	t,
		int	n,
		double *	s ) const

fill array for amplitude and time calculation

Definition at line 400 of file shaperdsp.cc.

{
  shaperdspshift_t t0 = shaperdspshift_t(t  - _toff, *this);
 
  for (int i = 0; i < n; i++, t0 += _tstride) {
    s[i] = ShaperDSP(t0);
  }
}

fillarray() [3/4]

void fillarray	(	int	n,
		dd_t *	s ) const

fill array for amplitude and time calculation

Definition at line 391 of file shaperdsp.cc.

{
  shaperdspshift_t t0 = _tzero;
 
  for (int i = 0; i < n; i++, t0 += _tstride) {
    s[i] = ddShaperDSP(t0);
  }
}

fillarray() [4/4]

void fillarray	(	int	n,
		double *	s ) const

fill array for amplitude and time calculation

Definition at line 382 of file shaperdsp.cc.

{
  shaperdspshift_t t0 = _tzero;
 
  for (int i = 0; i < n; i++, t0 += _tstride) {
    s[i] = ShaperDSP(t0);
  }
}

fillvector() [1/4]

void fillvector	(	double	t,
		std::vector< dd_t > &	s ) const

fill vector with response function values and its derivative

Definition at line 433 of file shaperdsp.cc.

{
  fillarray(t, s.size(), s.data());
}

fillvector() [2/4]

void fillvector	(	double	t,
		std::vector< double > &	s ) const

fill vector with response function values and its derivative

Definition at line 428 of file shaperdsp.cc.

{
  fillarray(t, s.size(), s.data());
}

fillvector() [3/4]

void fillvector

(

std::vector< dd_t > &

s

)

const

fill vector with response function values and its derivative

Definition at line 423 of file shaperdsp.cc.

{
  fillarray(s.size(), s.data());
}

fillvector() [4/4]

void fillvector

(

std::vector< double > &

s

)

const

fill vector with response function values and its derivative

Definition at line 418 of file shaperdsp.cc.

{
  fillarray(s.size(), s.data());
}

init() [1/3]

void init ( const double * s )

private

calculate some values for Sv123 function

Definition at line 338 of file shaperdsp.cc.

{
  init(s, -1);
}

init() [2/3]

void init ( const double * s, double unitscale )

private

calculate some values for Sv123 function

Definition at line 309 of file shaperdsp.cc.

{
  double t01 = s[2];
  double tb1 = s[3];
  double t02 = s[7];
  double tb2 = s[8];
  double td1 = s[1];
  double ts1 = s[4];
  Sv123_init(t01, tb1, t02, tb2, td1, ts1);
  _dt0 = 1 / s[6];
  _dt1 = 1 / s[2];
 
  _toff = s[0];
  _w0   = 1.0 - s[9];
  _w1   = 0.5 * s[9];
  _ccc  = s[5];
 
  _ccc *= unitscale;
  _cs0 *= unitscale;
  _cc0 *= unitscale;
  _cs1 *= unitscale;
  _cc1 *= unitscale;
  _ces *= unitscale;
  _ced *= unitscale;
 
  _tp.init(_filterdt, *this);
  _tm.init(-_filterdt, *this);
}

init() [3/3]

void init	(	const std::vector< double > &	s,
		double	u = -1 )

Initialization of the parameters response function.

Parameters

s	vector of shape parameters
u	'unitscale' of the waveform template This sets the normalization for the template amplitude to be unity. If the waveform templates are modified this needs to be adjusted as well. If unitscale is set to negative value, it is calculated dynamically based on shape parameters.

Definition at line 343 of file shaperdsp.cc.

{
  if (s.size() == 10)
    init(s.data(), unitscale);
  else
    init(_defs, unitscale);
}

nextseed()

void nextseed

(

)

substruct toffset to tzero

Definition at line 377 of file shaperdsp.cc.

{
  _tzero += _toffset;
}

operator()() [1/2]

double operator()	(	const double *	x,
		double *	)

TF1 ROOT interface.

Definition at line 357 of file shaperdsp.cc.

{
  return (*this)(x[0]);
}

operator()() [2/2]

double operator()

(

double

t

)

const

wrapper of the function

Definition at line 351 of file shaperdsp.cc.

{
  shaperdspshift_t t0 = shaperdspshift_t(t  - _toff, *this);
  return ShaperDSP(t0);
}

setseedoffset()

void setseedoffset

(

double

dt

)

set timeoffset

Definition at line 367 of file shaperdsp.cc.

{
  _toffset.init(dt, *this);
}

settimeseed()

void settimeseed

(

double

t0

)

set initial time

Definition at line 372 of file shaperdsp.cc.

{
  _tzero.init(t0 - _toff, *this);
}

settimestride()

void settimestride

(

double

dt

)

set grid step for function calculation

Definition at line 362 of file shaperdsp.cc.

{
  _tstride.init(dt, *this);
}

ShaperDSP()

double ShaperDSP ( const shaperdspshift_t & t0 ) const

private

calculate response function

Definition at line 274 of file shaperdsp.cc.

{
  const sv123shift_t& t00 = static_cast<const ShaperDSP_t::sv123shift_t&>(t0);
  double f = Sv123_filtered(t00);
  if (t0.t > 0) {
    double z = t0.t * _dt1;
    double z2 = z * z;
    double odd = 1 + z2 * (1 / 6. + z2 * (1 / 120.));
    double evn = 1 + z2 * (1 / 2. + z2 * (1 / 24.));
    double texp = evn + z * odd;
    double df = _ccc * t0.et0 * (1 - t0.et1 * texp);
    f -= df;
  }
  return f;
}

Sv123()

double Sv123 ( const sv123shift_t & c ) const

private

calculate Sv123 function

Definition at line 219 of file shaperdsp.cc.

{
  double f0 = (_cs0 * c.s0 + _cc0 * c.c0) * c.e0;
  double f1 = (_cs1 * c.s1 + _cc1 * c.c1) * c.e1;
  double fs =  _ces * c.es;
  double fd =  _ced * c.ed;
 
  return (f0 + f1) + (fd - fs);
}

Sv123_filtered()

double Sv123_filtered ( const sv123shift_t & t0 ) const

private

Numerical calculation of the time convolution.

Integration with histogram is used. Int g(s-t)f(t)dt=g(s)+a*(g(s+_filterdt)-g(s-_filterdt))/2

Definition at line 246 of file shaperdsp.cc.

{
  if (t0.validshift(_tp)) {
    double ft0p = Sv123(t0 + _tp), ft0 = 0, ft0m = 0;
    if (t0.t > 0) {
      ft0 = Sv123(t0);
      if (t0.validshift(_tm))
        ft0m = Sv123(t0 + _tm);
    }
    return _w0 * ft0 + _w1 * (ft0p + ft0m);
  }
  return 0;
}

Sv123_init()

void Sv123_init ( double t01, double tb1, double t02, double tb2, double td1, double ts1 )

private

calculate some values for Sv123 function

Definition at line 125 of file shaperdsp.cc.

{
  double  dks0, dks1, dksm,
          dw0, dw1, dwp, dwm, das1, dac1, das0, dac0, dzna, dksm2, ds, dd,
          dcs0, dsn0, dzn0, td, ts, dr,
          dcs0s, dsn0s, dcs0d, dsn0d, dcs1s, dsn1s, dcs1d, dsn1d;
 
  dr = (ts1 - td1) / td1;
  if (std::abs(dr) >= 0.0000001) {
    td = td1;
    ts = ts1;
  } else {
    td = td1;
    if (ts1 > td1) {
      ts = td1 * 1.00001;
    } else {
      ts = td1 * 0.99999;
    }
  }
 
  dr = ((t01 - t02) * (t01 - t02) + (tb1 - tb2) * (tb1 - tb2)) / (t01 * t01 + tb1 * tb1);
  dks0 = 1 / t01;
  dks1 = 1 / t02;
  if (dr < 0.0000000001) {
    if (dks0 > dks1) {
      dks0 = dks1 * 1.00001;
    } else {
      dks0 = dks1 * 0.99999;
    }
  }
 
  dksm = dks1 - dks0;
 
  ds = 1 / ts;
  dd = 1 / td;
 
  dw0 = 1 / tb1;
  dw1 = 1 / tb2;
  dwp = dw0 + dw1;
  dwm = dw1 - dw0;
 
  dksm2 = dksm * dksm;
 
  dzna = (dksm2 + dwm * dwm) * (dksm2 + dwp * dwp);
 
  das0 = dw1 * (dksm2 + dwp * dwm);
  dac0 = -2 * dksm * dw0 * dw1;
  das1 = dw0 * (dksm2 - dwp * dwm);
  dac1 = -dac0;
 
  dsn0 = (ds - dks0);
  dcs0 = -dw0;
  dzn0 = dcs0 * dcs0 + dsn0 * dsn0;
 
  dsn0s = (dsn0 * das0 - dcs0 * dac0) / dzn0;
  dcs0s = (dcs0 * das0 + dsn0 * dac0) / dzn0;
 
  dsn0 = (ds - dks1);
  dcs0 = -dw1;
  dzn0 = dcs0 * dcs0 + dsn0 * dsn0;
 
  dsn1s = (dsn0 * das1 - dcs0 * dac1) / dzn0;
  dcs1s = (dcs0 * das1 + dsn0 * dac1) / dzn0;
 
  dsn0 = (dd - dks0);
  dcs0 = -dw0;
  dzn0 = dcs0 * dcs0 + dsn0 * dsn0;
 
  dsn0d = (dsn0 * das0 - dcs0 * dac0) / dzn0;
  dcs0d = (dcs0 * das0 + dsn0 * dac0) / dzn0;
 
  dsn0 = (dd - dks1);
  dcs0 = -dw1;
  dzn0 = dcs0 * dcs0 + dsn0 * dsn0;
 
  dsn1d = (dsn0 * das1 - dcs0 * dac1) / dzn0;
  dcs1d = (dcs0 * das1 + dsn0 * dac1) / dzn0;
 
  _dw0  = dw0;
  _dw1  = dw1;
  _dks0 = dks0;
  _dks1 = dks1;
  _ds   = ds;
  _dd   = dd;
 
  double sc = 1 / (dzna * (ts - td));
  _cs0 = sc * (dsn0s - dsn0d);
  _cc0 = sc * (dcs0s - dcs0d);
  _cs1 = sc * (dsn1s - dsn1d);
  _cc1 = sc * (dcs1s - dcs1d);
  _ces = sc * (dcs0s + dcs1s);
  _ced = sc * (dcs0d + dcs1d);
}

Member Data Documentation

_cc0

double _cc0 = realNaN

private

linear coefficient before cos of the first Bessel stage

Definition at line 81 of file shaperdsp.h.

_cc1

double _cc1 = realNaN

private

linear coefficient before cos of the second Bessel stage

Definition at line 85 of file shaperdsp.h.

_ccc

double _ccc = realNaN

private

exponent factor for tail part of the signal

Definition at line 115 of file shaperdsp.h.

_ced

double _ced = realNaN

private

linear coefficient before second part of tail section

Definition at line 89 of file shaperdsp.h.

_ces

double _ces = realNaN

private

linear coefficient before first part of tail section

Definition at line 87 of file shaperdsp.h.

_cs0

double _cs0 = realNaN

private

linear coefficient before sin of the first Bessel stage

Definition at line 79 of file shaperdsp.h.

_cs1

double _cs1 = realNaN

private

linear coefficient before sin of the second Bessel stage

Definition at line 83 of file shaperdsp.h.

_dd

double _dd = realNaN

private

inverse time of the differential stage

Definition at line 101 of file shaperdsp.h.

_defs

const double _defs = {0.5, 0.6483, 0.4017, 0.3741, 0.8494, 0.00144547, 4.7071, 0.8156, 0.5556, 0.2752}

staticprivate

parameters of the response function that use as default

Definition at line 73 of file shaperdsp.h.

_dks0

double _dks0 = realNaN

private

decrement of the first Bessel stage

Definition at line 95 of file shaperdsp.h.

_dks1

double _dks1 = realNaN

private

decrement of the second Bessel stage

Definition at line 97 of file shaperdsp.h.

_ds

double _ds = realNaN

private

inverse scintillation decay time

Definition at line 99 of file shaperdsp.h.

_dt0

double _dt0 = realNaN

private

coefficient for first exponent factor

Definition at line 103 of file shaperdsp.h.

_dt1

double _dt1 = realNaN

private

coefficient for second exponent factor

Definition at line 105 of file shaperdsp.h.

_dw0

double _dw0 = realNaN

private

circular frequency of the first Bessel stage

Definition at line 91 of file shaperdsp.h.

_dw1

double _dw1 = realNaN

private

circular frequency of the second Bessel stage

Definition at line 93 of file shaperdsp.h.

_filterdt

double _filterdt = 0.2

staticconstexprprivate

time shift that include in response function for numerical calculation time convolutions.

Integration with histogram is used. Int g(s-t)f(t)dt=g(s)+a*(g(s+_filterdt)-g(s-_filterdt))/2

Definition at line 76 of file shaperdsp.h.

_tm

sv123shift_t _tm

private

• _filterdt

Definition at line 121 of file shaperdsp.h.

_toff

double _toff = realNaN

private

time offset

Definition at line 108 of file shaperdsp.h.

_toffset

shaperdspshift_t _toffset

private

time offset

Definition at line 127 of file shaperdsp.h.

_tp

sv123shift_t _tp

private

_filterdt

Definition at line 118 of file shaperdsp.h.

_tstride

shaperdspshift_t _tstride

private

time step of the grid for response function calculation

Definition at line 124 of file shaperdsp.h.

_tzero

shaperdspshift_t _tzero

private

initial time

Definition at line 129 of file shaperdsp.h.

_w0

double _w0 = realNaN

private

weight coefficient at sv123(t) = (1-a)

Definition at line 110 of file shaperdsp.h.

_w1

double _w1 = realNaN

private

weight coefficient at sv123(t+_filterdt) +sv123(t-_filterdt) = a/2

Definition at line 112 of file shaperdsp.h.

---

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

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