PulseFitWithFunction Class Reference

#include <PulseFitWithFunction.h>

Inheritance diagram for PulseFitWithFunction:

Public Member Functions
virtual double	doFit (double *)
double	getAmpl ()
double	getAmpl_parab ()
double	getMax_parab ()
int	getSampMax_parab ()
double	getTime ()
double	getTime_parab ()
virtual void	init (int, int, int, int, double, double)
	PulseFitWithFunction ()
	~PulseFitWithFunction () override

Public Attributes
double	fAmp_fitted_max
double	fFunc_max
int	fNumber_samp_max
double	fSigma_ped
double	fTim_fitted_max
double	fTim_max
double	fValue_tim_max

Private Member Functions
double	Electronic_shape (double)
double	Fit_electronic (int, double *, double)
void	Fit_parab (double *, int, int, double *)

Private Attributes
double	amp_max
double	amp_parab
double	fAlpha
double	fAlpha_beam
double	fAlpha_laser
double	fBeta
double	fBeta_beam
double	fBeta_laser
int	fNb_iter
int	fNsamples
int	fNum_samp_after_max
int	fNum_samp_bef_max
int	imax
double	tim_parab

Detailed Description

Definition at line 26 of file PulseFitWithFunction.h.

Constructor & Destructor Documentation

PulseFitWithFunction()

PulseFitWithFunction::PulseFitWithFunction

(

)

Definition at line 26 of file PulseFitWithFunction.cc.

References fNsamples, fNum_samp_after_max, and fNum_samp_bef_max.

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

~PulseFitWithFunction()

PulseFitWithFunction::~PulseFitWithFunction ( )

override

Definition at line 33 of file PulseFitWithFunction.cc.

{}

Member Function Documentation

doFit()

double PulseFitWithFunction::doFit ( double *  adc )

virtual

Definition at line 60 of file PulseFitWithFunction.cc.

References gpuClustering::adc, amp_max, amp_parab, nano_mu_local_reco_cff::chi2, fFunc_max, Fit_electronic(), Fit_parab(), fNsamples, fNumber_samp_max, fTim_max, fValue_tim_max, imax, createfilelist::int, and tim_parab.

Referenced by EcalPerEvtLaserAnalyzer::endJob(), and EcalLaserAnalyzer::endJob().

                                              {
  double parout[4];  // amp_max ;
                     //double amp_parab , tim_parab ;
  double chi2;
  //int imax ;
  //
  // first  one has to get starting point first with parabolic fun// //
  Fit_parab(&adc[0], 3, fNsamples, parout);
  amp_parab = parout[0];
  tim_parab = parout[1];
  imax = (int)parout[2];
  amp_max = parout[3];
  //printf("amp_parab= %f tim_parab=%f amp_max=%f imax=%d\n",amp_parab,tim_parab,amp_max,imax);
  fNumber_samp_max = imax;
  //

  if (amp_parab < 1.) {
    tim_parab = (double)imax;
    amp_parab = amp_max;
  }
  //
  fValue_tim_max = tim_parab;
  fFunc_max = amp_parab;
  fTim_max = tim_parab;

  //  here to fit maximum amplitude and time of arrival ...

  chi2 = Fit_electronic(0, &adc[0], 8.);
  // adc is an array to be filled with samples
  // 0 is for Laser (1 for electron)
  // 8 is for sigma of pedestals
  // which (can be computed)

  //  double amplitude = fAmp_fitted_max ; // amplitude fitted
  // double time = fTim_fitted_max ; // time fitted

  return chi2;  // return amplitude fitted
}

Electronic_shape()

double PulseFitWithFunction::Electronic_shape ( double  tim )

private

Definition at line 210 of file PulseFitWithFunction.cc.

References dt, fAlpha, fBeta, fFunc_max, fTim_max, and HPSPFTaus_cff::variable.

Referenced by Fit_electronic().

                                                        {
  // electronic function (from simulation) to fit ECAL pulse shape
  double func_electronic, dtsbeta, variable, puiss;
  double albet = fAlpha * fBeta;
  if (albet <= 0)
    return ((Double_t)0.);
  double dt = tim - fTim_max;
  if (dt > -albet) {
    dtsbeta = dt / fBeta;
    variable = 1. + dt / albet;
    puiss = TMath::Power(variable, fAlpha);
    func_electronic = fFunc_max * puiss * TMath::Exp(-dtsbeta);
  } else
    func_electronic = 0.;
  //
  return func_electronic;
}

Fit_electronic()

double PulseFitWithFunction::Fit_electronic ( int  data, double *  adc_to_fit, double  sigmas_sample  )

private

Definition at line 103 of file PulseFitWithFunction.cc.

References nano_mu_local_reco_cff::chi2, data, dumpMFGeometry_cfg::delta, makePileupJSON::denom, dt, Electronic_shape(), fAlpha, fAlpha_beam, fAlpha_laser, fAmp_fitted_max, fBeta, fBeta_beam, fBeta_laser, fFunc_max, fNb_iter, fNum_samp_after_max, fNum_samp_bef_max, fNumber_samp_max, fTim_fitted_max, fTim_max, EcalMonitorTask_cff::func, fValue_tim_max, mps_fire::i, createfilelist::int, HPSPFTaus_cff::variable, and testProducerWithPsetDescEmpty_cfi::z2.

Referenced by doFit().

                                                                                              {
  // fit electronic function from simulation
  // parameters fAlpha and fBeta are fixed and fit is providing
  // the maximum amplitude ( fAmp_fitted_max ) and the time of
  // the maximum amplitude ( fTim_fitted_max)
  // initialization of parameters
  double chi2 = 0;
  double d_alpha, d_beta;
  // first initialize parameters fAlpha and fBeta ( depending of beam or laser)
  fAlpha = fAlpha_laser;
  fBeta = fBeta_laser;
  if (data == 1) {
    fAlpha = fAlpha_beam;
    fBeta = fBeta_beam;
  }
  //
  fAmp_fitted_max = 0.;
  fTim_fitted_max = 0.;
  double un_sur_sigma = 1.;
  double variation_func_max = 0.;
  double variation_tim_max = 0.;
  //
  if (fValue_tim_max > 20. || fValue_tim_max < 3.) {
    fValue_tim_max = fNumber_samp_max;
  }
  int num_fit_min = (int)(fValue_tim_max - fNum_samp_bef_max);
  int num_fit_max = (int)(fValue_tim_max + fNum_samp_after_max);
  //

  if (sigmas_sample > 0.)
    un_sur_sigma = 1. / sigmas_sample;

  double func, delta;
  //          Loop on iterations
  for (int iter = 0; iter < fNb_iter; iter++) {
    //          initialization inside iteration loop !
    chi2 = 0.;
    double d11 = 0.;
    double d12 = 0.;
    double d22 = 0.;
    double z1 = 0.;
    double z2 = 0.;
    fFunc_max += variation_func_max;
    fTim_max += variation_tim_max;
    int nsamp_used = 0;
    //

    // Then we loop on samples to be fitted

    for (int i = num_fit_min; i < num_fit_max + 1; i++) {
      // calculate function to be fitted

      func = Electronic_shape((double)i);

      // then calculate derivatives of function to be fitted
      double dt = (double)i - fTim_max;
      double alpha_beta = fAlpha * fBeta;

      if (dt > -alpha_beta) {
        double dt_sur_beta = dt / fBeta;

        double variable = (double)1. + dt / alpha_beta;
        double expo = TMath::Exp(-dt_sur_beta);

        double puissance = TMath::Power(variable, fAlpha);
        d_alpha = un_sur_sigma * puissance * expo;
        d_beta = fFunc_max * d_alpha * dt_sur_beta / (alpha_beta * variable);
      } else {
        continue;
      }

      nsamp_used++;  // number of samples used inside the fit
      // compute matrix elements D (symetric --> d12 = d21 )
      d11 += d_alpha * d_alpha;
      d12 += d_alpha * d_beta;
      d22 += d_beta * d_beta;
      // compute delta
      delta = (adc_to_fit[i] - func) * un_sur_sigma;
      // compute vector elements Z
      z1 += delta * d_alpha;
      z2 += delta * d_beta;
      chi2 += delta * delta;
    }  //             end of loop on samples
    double denom = d11 * d22 - d12 * d12;
    if (denom == 0.) {
      //printf( "attention denom = 0 signal pas fitte \n") ;
      return 101;
    }
    if (nsamp_used < 3) {
      //printf( "Attention nsamp = %d ---> no function fit provided \n",nsamp_used) ;
      return 102;
    }
    // compute variations of parameters fAmp_max and fTim_max
    variation_func_max = (z1 * d22 - z2 * d12) / denom;
    variation_tim_max = (-z1 * d12 + z2 * d11) / denom;
    chi2 = chi2 / ((double)nsamp_used - 2.);
  }  //      end of loop on iterations
  // results of the fit are calculated
  fAmp_fitted_max = fFunc_max + variation_func_max;
  fTim_fitted_max = fTim_max + variation_tim_max;
  //
  return chi2;
}

Fit_parab()

void PulseFitWithFunction::Fit_parab ( double *  , int  , int  , double *    )

private

Definition at line 227 of file PulseFitWithFunction.cc.

References makePileupJSON::denom, dt, imax, dqmdumpme::k, and TrackingDataMCValidation_Standalone_cff::nmin.

Referenced by doFit().

                                                                                             {
  /* Now we calculate the parabolic adjustement in order to get        */
  /*    maximum and time max                                           */

  double denom, dt, amp1, amp2, amp3;
  double ampmax = 0.;
  int imax = 0;
  int k;
  /*
                                                                   */
  for (k = nmin; k < nmax; k++) {
    //printf("ampl[%d]=%f\n",k,ampl[k]);
    if (ampl[k] > ampmax) {
      ampmax = ampl[k];
      imax = k;
    }
  }
  amp1 = ampl[imax - 1];
  amp2 = ampl[imax];
  amp3 = ampl[imax + 1];
  denom = 2. * amp2 - amp1 - amp3;
  /*                                         */
  if (denom > 0.) {
    dt = 0.5 * (amp3 - amp1) / denom;
  } else {
    //printf("denom =%f\n",denom)  ;
    dt = 0.5;
  }
  /*                                             */
  /* ampmax correspond au maximum d'amplitude parabolique et dt        */
  /* decalage en temps par rapport au sample maximum soit k + dt       */

  parout[0] = amp2 + (amp3 - amp1) * dt * 0.25;
  parout[1] = (double)imax + dt;
  parout[2] = (double)imax;
  parout[3] = ampmax;
  return;
}

getAmpl()

double PulseFitWithFunction::getAmpl ( )

inline

Definition at line 52 of file PulseFitWithFunction.h.

References fAmp_fitted_max.

Referenced by EcalPerEvtLaserAnalyzer::endJob(), and EcalLaserAnalyzer::endJob().

{ return fAmp_fitted_max; }

getAmpl_parab()

double PulseFitWithFunction::getAmpl_parab ( )

inline

Definition at line 49 of file PulseFitWithFunction.h.

References amp_parab.

{ return amp_parab; }

getMax_parab()

double PulseFitWithFunction::getMax_parab ( )

inline

Definition at line 55 of file PulseFitWithFunction.h.

References amp_max.

{ return amp_max; }

getSampMax_parab()

int PulseFitWithFunction::getSampMax_parab ( )

inline

Definition at line 56 of file PulseFitWithFunction.h.

References imax.

{ return imax; }

getTime()

double PulseFitWithFunction::getTime ( )

inline

Definition at line 53 of file PulseFitWithFunction.h.

References fTim_fitted_max.

Referenced by EcalPerEvtLaserAnalyzer::endJob(), and EcalLaserAnalyzer::endJob().

{ return fTim_fitted_max; }

getTime_parab()

double PulseFitWithFunction::getTime_parab ( )

inline

Definition at line 50 of file PulseFitWithFunction.h.

References tim_parab.

{ return tim_parab; }

init()

void PulseFitWithFunction::init ( int  n_samples, int  samplb, int  sampla, int  niter, double  alfa, double  beta  )

virtual

Definition at line 37 of file PulseFitWithFunction.cc.

References HLT_2024v12_cff::beta, fAlpha_beam, fAlpha_laser, fBeta_beam, fBeta_laser, fNb_iter, fNsamples, fNum_samp_after_max, and fNum_samp_bef_max.

Referenced by EcalPerEvtLaserAnalyzer::endJob(), and EcalLaserAnalyzer::endJob().

                                                                                                          {
  //printf("\n");
  //printf(" =========================================================\n");
  //printf(" ==     Initialising the function fit method            ==\n");
  //printf(" ==          PulseFitWithFunction::init                 ==\n");
  //printf(" ==                                                     ==\n");

  fNsamples = n_samples;
  fAlpha_laser = alfa;
  fBeta_laser = beta;
  fAlpha_beam = 0.98;
  fBeta_beam = 2.04;
  fNb_iter = niter;
  fNum_samp_bef_max = samplb;
  fNum_samp_after_max = sampla;
  //printf(" ==  # samples used = %3d                               ==\n",fNsamples);
  //printf(" ==  #sample before max= %1d and #sample after maximum= %1d ==\n",fNum_samp_bef_max,fNum_samp_after_max);
  //printf(" ==           alpha= %5.4f       beta= %5.4f          ==\n",fAlpha_laser,fBeta_laser);
  //printf(" =========================================================\n\n");
  return;
}

Member Data Documentation

amp_max

double PulseFitWithFunction::amp_max

private

Definition at line 59 of file PulseFitWithFunction.h.

Referenced by doFit(), and getMax_parab().

amp_parab

double PulseFitWithFunction::amp_parab

private

Definition at line 59 of file PulseFitWithFunction.h.

Referenced by doFit(), and getAmpl_parab().

fAlpha

double PulseFitWithFunction::fAlpha

private

Definition at line 68 of file PulseFitWithFunction.h.

Referenced by Electronic_shape(), and Fit_electronic().

fAlpha_beam

double PulseFitWithFunction::fAlpha_beam

private

Definition at line 66 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), and init().

fAlpha_laser

double PulseFitWithFunction::fAlpha_laser

private

Definition at line 64 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), and init().

fAmp_fitted_max

double PulseFitWithFunction::fAmp_fitted_max

Definition at line 43 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), and getAmpl().

fBeta

double PulseFitWithFunction::fBeta

private

Definition at line 69 of file PulseFitWithFunction.h.

Referenced by Electronic_shape(), and Fit_electronic().

fBeta_beam

double PulseFitWithFunction::fBeta_beam

private

Definition at line 67 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), and init().

fBeta_laser

double PulseFitWithFunction::fBeta_laser

private

Definition at line 65 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), and init().

fFunc_max

double PulseFitWithFunction::fFunc_max

Definition at line 41 of file PulseFitWithFunction.h.

Referenced by doFit(), Electronic_shape(), and Fit_electronic().

fNb_iter

int PulseFitWithFunction::fNb_iter

private

Definition at line 70 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), and init().

fNsamples

int PulseFitWithFunction::fNsamples

private

Definition at line 62 of file PulseFitWithFunction.h.

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

fNum_samp_after_max

int PulseFitWithFunction::fNum_samp_after_max

private

Definition at line 72 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), init(), and PulseFitWithFunction().

fNum_samp_bef_max

int PulseFitWithFunction::fNum_samp_bef_max

private

Definition at line 71 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), init(), and PulseFitWithFunction().

fNumber_samp_max

int PulseFitWithFunction::fNumber_samp_max

Definition at line 46 of file PulseFitWithFunction.h.

Referenced by doFit(), and Fit_electronic().

fSigma_ped

double PulseFitWithFunction::fSigma_ped

Definition at line 47 of file PulseFitWithFunction.h.

fTim_fitted_max

double PulseFitWithFunction::fTim_fitted_max

Definition at line 44 of file PulseFitWithFunction.h.

Referenced by Fit_electronic(), and getTime().

fTim_max

double PulseFitWithFunction::fTim_max

Definition at line 42 of file PulseFitWithFunction.h.

Referenced by doFit(), Electronic_shape(), and Fit_electronic().

fValue_tim_max

double PulseFitWithFunction::fValue_tim_max

Definition at line 45 of file PulseFitWithFunction.h.

Referenced by doFit(), and Fit_electronic().

imax

int PulseFitWithFunction::imax

private

Definition at line 60 of file PulseFitWithFunction.h.

Referenced by doFit(), Fit_parab(), and getSampMax_parab().

tim_parab

double PulseFitWithFunction::tim_parab

private

Definition at line 59 of file PulseFitWithFunction.h.

Referenced by doFit(), and getTime_parab().