LowPassFilterTiming Class Reference

#include <LowPassFilterTiming.h>

Public Member Functions
unsigned	nParameters () const
double	operator() (double currentIn, const double *params, unsigned nParams) const

Detailed Description

Definition at line 4 of file LowPassFilterTiming.h.

Member Function Documentation

unsigned LowPassFilterTiming::nParameters

(

)

const

Definition at line 37 of file LowPassFilterTiming.cc.

References NPARAMETERS.

{
    return NPARAMETERS;
}

double LowPassFilterTiming::operator()	(	double	currentIn,
		const double *	params,
		unsigned	nParams
	)		const

Definition at line 50 of file LowPassFilterTiming.cc.

References delta, JetChargeProducer_cfi::exp, cmsBatch::log, NPARAMETERS, mps_fire::result, and two_joined_lines().

{
    assert(nParams == NPARAMETERS);
    assert(pars);
    unsigned ipar = 0;
    const double logVbias = pars[ipar++];
    const double logTauMin = pars[ipar++];

    // The middle of the join region. Not a log actually,
    // it is simple in the log space.
    const double logV0     = pars[ipar++];

    // Log of the width of the join region
    const double logDelta  = pars[ipar++];

    // Log of the minus negative slope
    const double slopeLog  = pars[ipar++];

    // Log of the maximum delay time (cutoff)
    const double tauMax    = pars[ipar++];
    assert(ipar == NPARAMETERS);

    // What happens for large (in magnitude) negative voltage inputs?
    const double Vbias = exp(logVbias);
    const double shiftedV = v + Vbias;
    if (shiftedV <= 0.0)
        return tauMax;

    const double lg = log(shiftedV);
    const double delta = exp(logDelta);
    const double tauMin = exp(logTauMin);
    double result = two_joined_lines(lg, -exp(slopeLog), tauMin,
                                     0.0, tauMin, logV0, delta);
    if (result > tauMax)
        result = tauMax;
    return result;
}