QIEInputPulse.h

 
#ifndef TRIGSCINT_QIEINPUTPULSE_H
#define TRIGSCINT_QIEINPUTPULSE_H
 
//---< C++ >---//
#include <vector>
 
namespace trigscint {
 
class QIEInputPulse {
 public:
  virtual ~QIEInputPulse() = default;
  float Eval(float T);
 
  virtual float EvalSingle(float T, int id) = 0;
 
  virtual float Integrate(float T1, float T2) = 0;
 
  virtual float Derivative(float T, int id) = 0;
 
  virtual float Max(int id) = 0;
 
  void AddPulse(float toff, float ampl);
 
  int GetNPulses() { return ampl_.size(); }
 
 protected:
  std::vector<float> toff_;
 
  std::vector<float> ampl_;
};
 
class Bimoid : public QIEInputPulse {
 public:
  Bimoid(float start, float qq);
 
  virtual ~Bimoid() = default;
 
  float EvalSingle(float T, int id) override;
 
  float I_Int(float T, int id);
 
  float Integrate(float T1, float T2) override;
 
  float Max(int id) override;
 
  float Derivative(float T, int id) override;
 
 private:
  float rt_;
  float ft_;
};
 
class Expo : public QIEInputPulse {
 public:
  Expo();
 
  Expo(float k, float tmax);  
 
  virtual ~Expo() = default;
  float GetRise() { return rt_; }
 
  float GetFall() { return ft_; }
 
  void SetRiseFall(float rr, float ff);
 
  float EvalSingle(float T, int id) override;
 
  float Integrate(float T1, float T2) override;
 
  float Max(int id) override;
 
  float Derivative(float T, int id) override;
 
 private:
  float k_;
  float tmax_;
  float rt_{-1};
  float ft_{-1};
 
  float I_Int(float T, int id);
};
 
}  // namespace trigscint
#endif