HcalSiPMnonlinearity Class Reference

#include <HcalSiPMnonlinearity.h>

Public Member Functions
int	getPixelsFired (int inpes) const
double	getRecoCorrectionFactor (double inpixelsfired) const
	HcalSiPMnonlinearity (const std::vector< float > &pars)

Private Attributes
double	a2
double	b1
double	c0

Detailed Description

Definition at line 7 of file HcalSiPMnonlinearity.h.

Constructor & Destructor Documentation

HcalSiPMnonlinearity()

HcalSiPMnonlinearity::HcalSiPMnonlinearity ( const std::vector< float > &  pars )

inline

Definition at line 9 of file HcalSiPMnonlinearity.h.

References a2, cms::cuda::assert(), b1, and c0.

                                                     {
    assert(pars.size() == 3);
    c0 = (double)pars[0];
    b1 = (double)pars[1];
    a2 = (double)pars[2];
  }

Member Function Documentation

getPixelsFired()

int HcalSiPMnonlinearity::getPixelsFired

(

int

inpes

)

const

Definition at line 9 of file HcalSiPMnonlinearity.cc.

References a, a2, cms::cuda::assert(), b, b1, c, c0, mps_fire::i, createfilelist::int, w(), and z.

                                                        {
  gsl_complex z[3];
  double w = -inpes;
  // normalize params
  double a = a2 / w;
  double b = b1 / w;
  double c = c0 / w;
  int nroots = gsl_poly_complex_solve_cubic(a, b, c, &z[1], &z[2], &z[3]);
  assert(nroots);

  // all use cases tested over the full range of anticipated values;
  // the first root is always the right one.
  double realpix = 0;
  // find real roots
  for (int i = 0; i < 3; ++i) {
    if (z[i].dat[1] == 0) {
      realpix = z[i].dat[0];
      break;
    }
  }

  return realpix > 0 ? (int)(realpix + 0.5) : 0;
}

getRecoCorrectionFactor()

double HcalSiPMnonlinearity::getRecoCorrectionFactor ( double  inpixelsfired ) const

inline

Definition at line 17 of file HcalSiPMnonlinearity.h.

References a2, b1, c0, and x.

                                                                    {
    double x = inpixelsfired;
    return (a2 * x * x + b1 * x + c0);
  }

Member Data Documentation

a2

double HcalSiPMnonlinearity::a2

private

Definition at line 27 of file HcalSiPMnonlinearity.h.

Referenced by getPixelsFired(), getRecoCorrectionFactor(), and HcalSiPMnonlinearity().

b1

double HcalSiPMnonlinearity::b1

private

Definition at line 27 of file HcalSiPMnonlinearity.h.

Referenced by getPixelsFired(), getRecoCorrectionFactor(), and HcalSiPMnonlinearity().

c0

double HcalSiPMnonlinearity::c0

private

Definition at line 27 of file HcalSiPMnonlinearity.h.

Referenced by getPixelsFired(), getRecoCorrectionFactor(), and HcalSiPMnonlinearity().