TPNFit Class Reference

#include <TPNFit.h>

Inheritance diagram for TPNFit:

Public Member Functions
double	doFit (int, double *)
double	getAmpl ()
double	getTimax ()
void	init (int, int, int)
	TPNFit ()
	~TPNFit () override

Static Public Attributes
static constexpr int	NMAXSAMP2 = 50

Private Attributes
double	ampl
int	firstsample
int	fNsamples
int	fNum_samp_after_max
int	fNum_samp_bef_max
double	fv1 [NMAXSAMP2]
double	fv2 [NMAXSAMP2]
double	fv3 [NMAXSAMP2]
int	lastsample
double	t [NMAXSAMP2]
double	timeatmax
double	val [NMAXSAMP2]

Detailed Description

Definition at line 6 of file TPNFit.h.

Constructor & Destructor Documentation

TPNFit()

TPNFit::TPNFit

(

)

Definition at line 15 of file TPNFit.cc.

References fNsamples, fNum_samp_after_max, and fNum_samp_bef_max.

               {
  fNsamples = 0;
  fNum_samp_bef_max = 0;
  fNum_samp_after_max = 0;
}

~TPNFit()

TPNFit::~TPNFit ( )

override

Definition at line 22 of file TPNFit.cc.

{}

Member Function Documentation

doFit()

double TPNFit::doFit	(	int	maxsample,
		double *	adc
	)

Definition at line 39 of file TPNFit.cc.

References a, gpuClustering::adc, ampl, b, DummyCfis::c, isoTrack_cff::chi2, validate-o2o-wbm::f1, validate-o2o-wbm::f2, firstsample, fNsamples, fNum_samp_after_max, fNum_samp_bef_max, fv1, fv2, fv3, mps_fire::i, lastsample, timeatmax, val, and y.

Referenced by EcalPerEvtLaserAnalyzer::analyze(), EcalTestPulseAnalyzer::analyze(), EcalLaserAnalyzer::analyze(), and EcalLaserAnalyzer2::analyze().

                                               {
  double chi2 = 0.;
  ampl = 0.;
  timeatmax = 0.;

  //printf("maxsample=%d\n",maxsample);
  firstsample = maxsample - fNum_samp_bef_max;
  lastsample = maxsample + fNum_samp_after_max;

  if (firstsample <= 0)
    return 101;
  if (lastsample >= fNsamples)
    lastsample = fNsamples - 1;
  if (lastsample - firstsample < 1)
    return 102;
  int nval = lastsample - firstsample + 1;
  //printf("firstsample=%d lastsample=%d nval=%d\n",
  //                        firstsample,lastsample,nval);
  int testneg = 0;
  for (int kn = firstsample; kn <= lastsample; kn++) {
    //printf("adc[%d]=%f\n",kn,adc[kn]);
    if (adc[kn] < 0.)
      testneg = 1;
  }
  if (testneg == 1)
    return 103;

  for (int i = firstsample; i <= lastsample; i++) {
    val[i - firstsample] = adc[i];
    fv1[i - firstsample] = 1.;
    fv2[i - firstsample] = (double)(i);
    fv3[i - firstsample] = ((double)(i)) * ((double)(i));
  }

  TVectorD y(nval, val);
  //y.Print();
  TVectorD f1(nval, fv1);
  //f1.Print();
  TVectorD f2(nval, fv2);
  //f2.Print();
  TVectorD f3(nval, fv3);
  //f3.Print();

  double bj[3];
  bj[0] = f1 * y;
  bj[1] = f2 * y;
  bj[2] = f3 * y;
  TVectorD b(3, bj);
  //b.Print();

  double aij[9];
  aij[0] = f1 * f1;
  aij[1] = f1 * f2;
  aij[2] = f1 * f3;
  aij[3] = f2 * f1;
  aij[4] = f2 * f2;
  aij[5] = f2 * f3;
  aij[6] = f3 * f1;
  aij[7] = f3 * f2;
  aij[8] = f3 * f3;
  TMatrixD a(3, 3, aij);
  //a.Print();

  double det1;
  a.InvertFast(&det1);
  //a.Print();

  TVectorD c = a * b;
  //c.Print();

  double *par = c.GetMatrixArray();
  //printf("par0=%f par1=%f par2=%f\n",par[0],par[1],par[2]);
  if (par[2] != 0.) {
    timeatmax = -par[1] / (2. * par[2]);
    ampl = par[0] - par[2] * (timeatmax * timeatmax);
  }
  //printf("amp=%f time=%f\n",amp_max,timeatmax);

  if (ampl <= 0.) {
    ampl = adc[maxsample];
    return 1.;
  }

  if ((int)timeatmax > lastsample)
    return 103;
  if ((int)timeatmax < firstsample)
    return 103;

  return chi2;
}

getAmpl()

double TPNFit::getAmpl ( )

inline

Definition at line 32 of file TPNFit.h.

References ampl.

Referenced by EcalPerEvtLaserAnalyzer::analyze(), EcalTestPulseAnalyzer::analyze(), EcalLaserAnalyzer::analyze(), and EcalLaserAnalyzer2::analyze().

{ return ampl; }

getTimax()

double TPNFit::getTimax ( )

inline

Definition at line 33 of file TPNFit.h.

References timeatmax.

{ return timeatmax; }

init()

void TPNFit::init	(	int	nsamples,
		int	firstsample,
		int	lastsample
	)

Definition at line 24 of file TPNFit.cc.

References firstsample, fNsamples, fNum_samp_after_max, fNum_samp_bef_max, isotrackApplyRegressor::k, lastsample, NMAXSAMP2, and t.

Referenced by EcalPerEvtLaserAnalyzer::analyze(), EcalTestPulseAnalyzer::analyze(), EcalLaserAnalyzer::analyze(), and EcalLaserAnalyzer2::analyze().

                                                               {
  fNsamples = nsamples;
  fNum_samp_bef_max = firstsample;
  fNum_samp_after_max = lastsample;
  //printf("nsamples=%d firstsample=%d lastsample=%d\n",nsamples,firstsample,lastsample);

  if (fNsamples > NMAXSAMP2)
    printf("number of PN samples exceed maximum\n");

  for (int k = 0; k < NMAXSAMP2; k++)
    t[k] = (double)k;

  return;
}

Member Data Documentation

ampl

double TPNFit::ampl

private

Definition at line 18 of file TPNFit.h.

Referenced by doFit(), and getAmpl().

firstsample

int TPNFit::firstsample

private

Definition at line 15 of file TPNFit.h.

Referenced by doFit(), and init().

fNsamples

int TPNFit::fNsamples

private

Definition at line 11 of file TPNFit.h.

Referenced by doFit(), init(), and TPNFit().

fNum_samp_after_max

int TPNFit::fNum_samp_after_max

private

Definition at line 13 of file TPNFit.h.

Referenced by doFit(), init(), and TPNFit().

fNum_samp_bef_max

int TPNFit::fNum_samp_bef_max

private

Definition at line 12 of file TPNFit.h.

Referenced by doFit(), init(), and TPNFit().

fv1

double TPNFit::fv1[NMAXSAMP2]

private

Definition at line 17 of file TPNFit.h.

Referenced by doFit().

fv2

double TPNFit::fv2[NMAXSAMP2]

private

Definition at line 17 of file TPNFit.h.

Referenced by doFit().

fv3

double TPNFit::fv3[NMAXSAMP2]

private

Definition at line 17 of file TPNFit.h.

Referenced by doFit().

lastsample

int TPNFit::lastsample

private

Definition at line 15 of file TPNFit.h.

Referenced by doFit(), and init().

NMAXSAMP2

constexpr int TPNFit::NMAXSAMP2 = 50

static

Definition at line 8 of file TPNFit.h.

Referenced by init().

t

double TPNFit::t[NMAXSAMP2]

private

Definition at line 16 of file TPNFit.h.

Referenced by init().

timeatmax

double TPNFit::timeatmax

private

Definition at line 19 of file TPNFit.h.

Referenced by doFit(), and getTimax().

val

double TPNFit::val[NMAXSAMP2]

private

Definition at line 16 of file TPNFit.h.

Referenced by value.Value::__eq__(), value.Value::__iadd__(), value.Value::__idiv__(), value.Value::__isub__(), value.Value::__str__(), doFit(), and value.Value::relerr().