LowPassFilterTiming.cc File Reference

#include <cmath>
#include <cassert>
#include "FWCore/Utilities/interface/Exception.h"
#include "CalibCalorimetry/HcalAlgos/interface/LowPassFilterTiming.h"

Go to the source code of this file.

Macros
#define	NPARAMETERS   6U

Functions
static double	two_joined_lines (const double x_in, const double aleft, const double bleft, const double aright, const double bright, const double join_location, const double h)

Macro Definition Documentation

#define NPARAMETERS   6U

Definition at line 7 of file LowPassFilterTiming.cc.

Referenced by LowPassFilterTiming::nParameters(), and LowPassFilterTiming::operator()().

Function Documentation

static double two_joined_lines ( const double  x_in, const double  aleft, const double  bleft, const double  aright, const double  bright, const double  join_location, const double  h  )

static

Definition at line 14 of file LowPassFilterTiming.cc.

References a, b, EnergyCorrector::c, edmIntegrityCheck::d, h, and genVertex_cff::x.

Referenced by LowPassFilterTiming::operator()().

{
    assert(h >= 0.0);
    const double x = x_in - join_location;
    if (x <= -h)
        return aleft*x + bleft;
    else if (x >= h)
        return aright*x + bright;
    else
    {
        const double vleft = -h*aleft + bleft;
        const double vright = aright*h + bright;
        const double b = (aright - aleft)/4.0/h;
        const double d = (vright + vleft - 2*b*h*h)/2.0;
        const double a = (d + aright*h - b*h*h - vright)/(2.0*h*h*h);
        const double c = -(3*d + aright*h + b*h*h - 3*vright)/(2.0*h);
        return ((a*x + b)*x + c)*x + d;
    }
}