trigscint::Bimoid Class Reference

Pulse made out of difference of two sigmoids. More...

#include <QIEInputPulse.h>

Public Member Functions
	Bimoid (float start, float qq)
	Constructor.
virtual	~Bimoid ()=default
	Default Destructor.
float	evalSingle (float T, int id) override
	Evaluate the pulse at time T.
float	iInt (float T, int id)
	Indefinite integral at time T.
float	integrate (float T1, float T2) override
	Integrate the pulse from T1 to T2.
float	max (int id) override
	maximum of the pulse
float	derivative (float T, int id) override
	Differentiate pulse at time T.
Public Member Functions inherited from trigscint::QIEInputPulse
virtual	~QIEInputPulse ()=default
	Destructor.
float	eval (float T)
	Evaluate the pulse train at time T.
void	addPulse (float toff, float ampl)
	To add a pulse to the collection.
int	getNPulses ()
	Get the number of pulses in the collection.

Private Attributes
float	rt_
	rise time
float	ft_
	fall time

Additional Inherited Members
Protected Attributes inherited from trigscint::QIEInputPulse
std::vector< float >	toff_
	collection of pulse time offsets
std::vector< float >	ampl_
	collection of pulse amplitudes

Detailed Description

Pulse made out of difference of two sigmoids.

Note

The pulse maximum is found numerically

Definition at line 74 of file QIEInputPulse.h.

Constructor & Destructor Documentation

Bimoid()

trigscint::Bimoid::Bimoid	(	float	start,
		float	qq )

Constructor.

Definition at line 32 of file QIEInputPulse.cxx.

                                     {
  rt_ = rise;
  ft_ = fall;
}

References ft_, and rt_.

Member Function Documentation

derivative()

float trigscint::Bimoid::derivative ( float T, int id )

overridevirtual

Differentiate pulse at time T.

Implements trigscint::QIEInputPulse.

Definition at line 86 of file QIEInputPulse.cxx.

                                        {
  // Normalization constant
  float nc = (ft_ - rt_) * log(2) / ampl_[id];
 
  float t = T - toff_[id];
  float e1 = exp(-t / rt_);
  float e2 = exp(-t / ft_);
 
  float v1 = e1 / (rt_ * pow(1 + e1, 2));
  float v2 = e2 / (ft_ * pow(1 + e2, 2));
 
  return ((v1 - v2) / nc);  // Actual derivative
}

References trigscint::QIEInputPulse::ampl_, ft_, rt_, and trigscint::QIEInputPulse::toff_.

Referenced by max().

evalSingle()

float trigscint::Bimoid::evalSingle ( float T, int id )

overridevirtual

Evaluate the pulse at time T.

Implements trigscint::QIEInputPulse.

Definition at line 37 of file QIEInputPulse.cxx.

                                        {
  if (T < toff_[id]) return 0;
  // Normalization constant
  float nc = (ft_ - rt_) * log(2) / ampl_[id];
 
  float y1 = 1 / (1 + exp((toff_[id] - T) / rt_));
  float y2 = 1 / (1 + exp((toff_[id] - T) / ft_));
  return ((y1 - y2) / nc);
}

References trigscint::QIEInputPulse::ampl_, ft_, rt_, and trigscint::QIEInputPulse::toff_.

iInt()

float trigscint::Bimoid::iInt	(	float	T,
		int	id )

Indefinite integral at time T.

Definition at line 57 of file QIEInputPulse.cxx.

                                  {
  if (T <= toff_[id]) return 0;
  // Normalization constant
  float nc = (ft_ - rt_) * log(2) / ampl_[id];
 
  float t = T - toff_[id];  // time relative to offset
 
  float ii =  // Integral
      rt_ * log(1 + exp((t - toff_[id]) / rt_)) -
      ft_ * log(1 + exp((t - toff_[id]) / ft_));
 
  return ii / nc;
}

References trigscint::QIEInputPulse::ampl_, ft_, rt_, and trigscint::QIEInputPulse::toff_.

Referenced by integrate().

integrate()

float trigscint::Bimoid::integrate ( float T1, float T2 )

overridevirtual

Integrate the pulse from T1 to T2.

Implements trigscint::QIEInputPulse.

Definition at line 47 of file QIEInputPulse.cxx.

                                          {
  float val = 0;
  for (int id = 0; id < ampl_.size(); id++) {
    if (ampl_[id] > 0 && T2 > toff_[id]) {
      val += iInt(T2, id) - iInt(T1, id);
    }
  }
  return val;
}

References trigscint::QIEInputPulse::ampl_, iInt(), and trigscint::QIEInputPulse::toff_.

max()

float trigscint::Bimoid::max ( int id )

overridevirtual

maximum of the pulse

Implements trigscint::QIEInputPulse.

Definition at line 71 of file QIEInputPulse.cxx.

                        {
  float a = 0;
  float b = 50;
  float mx = (a + b) / 2;  // maximum
 
  while (std::abs(derivative(mx, id)) >= 1e-5) {
    if (derivative(a, id) * derivative(mx, id) > 0) {
      a = mx;
    } else
      b = mx;
    mx = (a + b) / 2;
  }
  return (mx);
}

References derivative().

Member Data Documentation

ft_

float trigscint::Bimoid::ft_

private

fall time

Definition at line 111 of file QIEInputPulse.h.

Referenced by Bimoid(), derivative(), evalSingle(), and iInt().

rt_

float trigscint::Bimoid::rt_

private

rise time

Definition at line 109 of file QIEInputPulse.h.

Referenced by Bimoid(), derivative(), evalSingle(), and iInt().

---

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

• TrigScint/include/TrigScint/QIEInputPulse.h
• TrigScint/src/TrigScint/QIEInputPulse.cxx