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 ()
virtual	~TPNFit ()

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

Detailed Description

Definition at line 8 of file TPNFit.h.

Constructor & Destructor Documentation

TPNFit::TPNFit

(

)

Definition at line 16 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 ( )

virtual

Definition at line 25 of file TPNFit.cc.

{ 
}

Member Function Documentation

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

Definition at line 45 of file TPNFit.cc.

References a, ampl, b, EnergyCorrector::c, python.connectstrParser::f1, python.connectstrParser::f2, python.connectstrParser::f3, firstsample, fNsamples, fNum_samp_after_max, fNum_samp_bef_max, fv1, fv2, fv3, i, lastsample, timeatmax, val, and y.

{
  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;
}

double TPNFit::getAmpl ( )

inline

Definition at line 35 of file TPNFit.h.

References ampl.

{return ampl;}

double TPNFit::getTimax ( )

inline

Definition at line 36 of file TPNFit.h.

References timeatmax.

{return timeatmax;}

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

Definition at line 29 of file TPNFit.cc.

References firstsample, fNsamples, fNum_samp_after_max, fNum_samp_bef_max, relval_steps::k, lastsample, NMAXSAMP2, mathSSE::return(), and t.

{
  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

double TPNFit::ampl

private

Definition at line 20 of file TPNFit.h.

Referenced by doFit(), and getAmpl().

int TPNFit::firstsample

private

Definition at line 17 of file TPNFit.h.

Referenced by doFit(), and init().

int TPNFit::fNsamples

private

Definition at line 13 of file TPNFit.h.

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

int TPNFit::fNum_samp_after_max

private

Definition at line 15 of file TPNFit.h.

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

int TPNFit::fNum_samp_bef_max

private

Definition at line 14 of file TPNFit.h.

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

double TPNFit::fv1[50]

private

Definition at line 19 of file TPNFit.h.

Referenced by doFit().

double TPNFit::fv2[50]

private

Definition at line 19 of file TPNFit.h.

Referenced by doFit().

double TPNFit::fv3[50]

private

Definition at line 19 of file TPNFit.h.

Referenced by doFit().

int TPNFit::lastsample

private

Definition at line 17 of file TPNFit.h.

Referenced by doFit(), and init().

double TPNFit::t[50]

private

Definition at line 18 of file TPNFit.h.

Referenced by init().

double TPNFit::timeatmax

private

Definition at line 21 of file TPNFit.h.

Referenced by doFit(), and getTimax().

double TPNFit::val[50]

private

Definition at line 18 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().