ThirdOrderDelayODE Class Reference

#include <ThirdOrderDelayODE.h>

Public Member Functions
void	calculate (double tau, double inputCurrent, double dIdt, double d2Id2t, const double *x, unsigned lenX, unsigned firstNode, double *derivative) const
unsigned	getPadeColumn () const
unsigned	getPadeRow () const
unsigned	nParameters () const
void	setParameters (const double *pars, unsigned nPars)
	ThirdOrderDelayODE (unsigned, unsigned)

Private Attributes
double	a_
double	b_
double	c_

Detailed Description

Definition at line 9 of file ThirdOrderDelayODE.h.

Constructor & Destructor Documentation

ThirdOrderDelayODE::ThirdOrderDelayODE ( unsigned  , unsigned    )

inline

Definition at line 12 of file ThirdOrderDelayODE.h.

References calculate(), funct::derivative(), metsig::tau, and x.

: a_(1.0) {}

Member Function Documentation

void ThirdOrderDelayODE::calculate	(	double	tau,
		double	inputCurrent,
		double	dIdt,
		double	d2Id2t,
		const double *	x,
		unsigned	lenX,
		unsigned	firstNode,
		double *	derivative
	)		const

Definition at line 7 of file ThirdOrderDelayODE.cc.

References a_, b_, c_, Exception, and mitigatedMETSequence_cff::U.

Referenced by ThirdOrderDelayODE().

{
    // Check input sanity
    if (lenX < firstNode + 3U) throw cms::Exception(
        "In ThirdOrderDelayODE::calculate: insufficient number of variables");
    if (tau <= 0.0) throw cms::Exception(
        "In ThirdOrderDelayODE::calculate: delay time is not positive");
    assert(x);
    assert(derivative);

    derivative[firstNode]   = x[firstNode+1];
    derivative[firstNode+1] = x[firstNode+2];
    derivative[firstNode+2] = 6.0/a_*(currentIn - x[firstNode] -
                                      c_*tau*x[firstNode+1] -
                                      b_/2.0*tau*tau*x[firstNode+2])/tau/tau/tau;
}

unsigned ThirdOrderDelayODE::getPadeColumn ( ) const

inline

Definition at line 20 of file ThirdOrderDelayODE.h.

References mitigatedMETSequence_cff::U.

{return 3U;}

unsigned ThirdOrderDelayODE::getPadeRow ( ) const

inline

Definition at line 19 of file ThirdOrderDelayODE.h.

References mitigatedMETSequence_cff::U.

{return 0U;}

unsigned ThirdOrderDelayODE::nParameters ( ) const

inline

Definition at line 21 of file ThirdOrderDelayODE.h.

References setParameters(), and mitigatedMETSequence_cff::U.

{return 3U;}

void ThirdOrderDelayODE::setParameters	(	const double *	pars,
		unsigned	nPars
	)

Definition at line 28 of file ThirdOrderDelayODE.cc.

References a_, b_, c_, JetChargeProducer_cfi::exp, and mitigatedMETSequence_cff::U.

Referenced by nParameters().

{
    assert(nPars == 3U);
    assert(pars);
    a_ = exp(pars[0]);
    b_ = exp(pars[1]);
    c_ = exp(pars[2]);
}

Member Data Documentation

double ThirdOrderDelayODE::a_

private

Definition at line 27 of file ThirdOrderDelayODE.h.

Referenced by calculate(), and setParameters().

double ThirdOrderDelayODE::b_

private

Definition at line 28 of file ThirdOrderDelayODE.h.

Referenced by calculate(), and setParameters().

double ThirdOrderDelayODE::c_

private

Definition at line 29 of file ThirdOrderDelayODE.h.

Referenced by calculate(), and setParameters().