ThirdOrderDelayODE.h

Go to the documentation of this file.

#ifndef CalibCalorimetry_HcalAlgos_ThirdOrderDelayODE_h_
#define CalibCalorimetry_HcalAlgos_ThirdOrderDelayODE_h_
 
//
// Equation a/6*tau^3*V_out''' + b/2*tau^2*V_out'' + c*tau*V_out' + V_out = V_in,
// with parameters "a", "b", and "c". a = 1, b = 1, c = 1 corresponds to the
// Pade table delay equation with row = 0 and column = 3.
//
class ThirdOrderDelayODE {
public:
  inline ThirdOrderDelayODE(unsigned /* r */, unsigned /* c */) : a_(1.0) {}
 
  void calculate(double tau,
                 double inputCurrent,
                 double dIdt,
                 double d2Id2t,
                 const double* x,
                 unsigned lenX,
                 unsigned firstNode,
                 double* derivative) const;
 
  inline unsigned getPadeRow() const { return 0U; }
  inline unsigned getPadeColumn() const { return 3U; }
  inline unsigned nParameters() const { return 3U; }
 
  // The parameters should be set to the logs of their actual values
  void setParameters(const double* pars, unsigned nPars);
 
private:
  double a_;
  double b_;
  double c_;
};
 
#endif  // CalibCalorimetry_HcalAlgos_ThirdOrderDelayODE_h_