CMS 3D CMS Logo

EcalUncalibRecHitFixedAlphaBetaAlgo< C > Class Template Reference

Template used to compute amplitude, pedestal, time jitter, chi2 of a pulse using an analytical function fit, with the pulse parameters alpha and beta fixed. More...

#include <RecoLocalCalo/EcalRecAlgos/interface/EcalUncalibRecHitFixedAlphaBetaAlgo.h>

Inheritance diagram for EcalUncalibRecHitFixedAlphaBetaAlgo< C >:

EcalUncalibRecHitRecAbsAlgo< C >

List of all members.

Public Member Functions

 EcalUncalibRecHitFixedAlphaBetaAlgo (int n_iter, int n_bef_max=1, int n_aft_max=3, float sigma_ped=1.1)
 EcalUncalibRecHitFixedAlphaBetaAlgo ()
virtual EcalUncalibratedRecHit makeRecHit (const C &dataFrame, const double *pedestals, const double *gainRatios, const EcalWeightSet::EcalWeightMatrix **weights, const EcalWeightSet::EcalChi2WeightMatrix **chi2Matrix)
 Compute parameters.
void SetAlphaBeta (double alpha, double beta)
void SetDynamicPedestal (bool dyn_pede)
void SetMinAmpl (double ampl)
virtual ~EcalUncalibRecHitFixedAlphaBetaAlgo ()

Private Member Functions

void InitFitParameters (double *samples, int max_sample)
float PerformAnalyticFit (double *samples, int max_sample)
double pulseShapeFunction (double t)

Private Attributes

double alfabeta_
float alpha_table_ [36][1701]
float beta_table_ [36][1701]
HepSymMatrix DM1_
bool doFit_
bool dyn_pedestal
double fAlpha_
double fAmp_max_
double fBeta_
int fNb_iter_
int fNum_samp_after_max_
int fNum_samp_bef_max_
double fPed_max_
float fSigma_ped
double fTim_max_
double MinAmpl_
HepVector temp_
double un_sur_sigma

Detailed Description

template<class C>
class EcalUncalibRecHitFixedAlphaBetaAlgo< C >

Template used to compute amplitude, pedestal, time jitter, chi2 of a pulse using an analytical function fit, with the pulse parameters alpha and beta fixed.

It follows a fast fit algorithms devolped on test beam data by P. Jarry If the pedestal is not given, it is calculated from the first 2 pre-samples; FIXME: conversion gainID (1,2,3) with gain factor (12,6,1) is hardcoded here !

Author:
A.Ghezzi

Definition at line 31 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Constructor & Destructor Documentation

template<class C>
EcalUncalibRecHitFixedAlphaBetaAlgo< C >::EcalUncalibRecHitFixedAlphaBetaAlgo (  )  [inline]

template<class C>
EcalUncalibRecHitFixedAlphaBetaAlgo< C >::EcalUncalibRecHitFixedAlphaBetaAlgo ( int  n_iter,
int  n_bef_max = 1,
int  n_aft_max = 3,
float  sigma_ped = 1.1 
) [inline]

template<class C>
virtual EcalUncalibRecHitFixedAlphaBetaAlgo< C >::~EcalUncalibRecHitFixedAlphaBetaAlgo (  )  [inline, virtual]

Definition at line 92 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

00094                                                     { };

Member Function Documentation

template<class C>
void EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters ( double *  samples,
int  max_sample 
) [inline, private]

Definition at line 208 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

References a, b, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::doFit_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAmp_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fPed_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fTim_max_, and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::MinAmpl_.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::makeRecHit(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

00208                                                                                                                 {
00209 
00210   // in a first attempt just use the value of the maximum sample 
00211   fAmp_max_ = samples[max_sample];
00212   fTim_max_ = max_sample;
00213   fPed_max_ = 0;
00214 
00215   // amplitude too low for fit to converge
00216   // timing set correctly is assumed
00217   if(fAmp_max_ <  MinAmpl_){
00218     fAmp_max_      = samples[5];
00219     double sumA    = samples[5]+samples[4]+samples[6];
00220     if(sumA != 0) { fTim_max_ = 5+(samples[6]-samples[4])/sumA; }
00221     else{ fTim_max_ = -993; }//-999+6
00222     doFit_ = false;
00223   }
00224   // if timing very badly off, that just use max sample
00225   else if(max_sample <1 || max_sample > 7)
00226     {    doFit_ = false;}
00227   else{
00228     //y=a*(x-xM)^2+b*(x-xM)+c
00229     doFit_ = true;
00230     float a = float(samples[max_sample-1]+samples[max_sample+1]-2*samples[max_sample])/2.;
00231     if(a==0){doFit_ =false; return;}
00232     float b = float(samples[max_sample+1]-samples[max_sample-1])/2.;
00233   
00234     fTim_max_ = max_sample - b/(2*a);
00235     fAmp_max_ =  samples[max_sample] - b*b/(4*a); 
00236   }  
00237   
00238 }

template<class C>
EcalUncalibratedRecHit EcalUncalibRecHitFixedAlphaBetaAlgo< C >::makeRecHit ( const C &  dataFrame,
const double *  pedestals,
const double *  gainRatios,
const EcalWeightSet::EcalWeightMatrix **  weights,
const EcalWeightSet::EcalChi2WeightMatrix **  chi2Matrix 
) [inline, virtual]

Compute parameters.

Implements EcalUncalibRecHitRecAbsAlgo< C >.

Definition at line 107 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

References EcalUncalibRecHitFixedAlphaBetaAlgo< C >::dyn_pedestal, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAmp_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fPed_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fTim_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters(), pedestal, and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

Referenced by EcalFixedAlphaBetaFitUncalibRecHitProducer::produce().

00110                                                                                                                                                   {
00111   double chi2_(-1.);
00112   
00113   //  double Gain12Equivalent[4]={0,1,2,12};
00114   double frame[C::MAXSAMPLES];// will contain the ADC values
00115   double pedestal =0;     // carries pedestal for gain12 i.e. gainId==1
00116   
00117   int gainId0 = 1;        // expected gainId at the beginning of dataFrame 
00118   int iGainSwitch = 0;    // flags whether there's any gainId other than gainId0
00119   int GainId= 0;          // stores gainId at every sample
00120   double maxsample(-1);   // ADC value of maximal ped-subtracted sample
00121   int imax(-1);           // sample number of maximal ped-subtracted sample
00122   bool external_pede = false;
00123   bool isSaturated = 0;   // flag reporting whether gain0 has been found
00124 
00125   // Get time samples checking for Gain Switch and pedestals
00126   if(pedestals){
00127     external_pede = true;
00128     if(dyn_pedestal) { pedestal = (double(dataFrame.sample(0).adc()) + double(dataFrame.sample(1).adc()))/2.;}
00129     else{ pedestal  = pedestals[0];}
00130     for(int iSample = 0; iSample < C::MAXSAMPLES; iSample++) {
00131         //create frame in adc gain 12 equivalent
00132         GainId = dataFrame.sample(iSample).gainId();
00133 
00134         // FIX-ME: warning: the vector pedestal is supposed to have in the order G12, G6 and G1
00135         // if GainId is zero treat it as 3 temporarily to protect against undefined
00136         // frame will be set to ~max of gain1
00137         if ( GainId == 0 )
00138           { 
00139             GainId = 3;
00140             isSaturated = 1;
00141           }
00142 
00143         if (GainId != gainId0) iGainSwitch = 1;
00144 
00145         if(GainId==gainId0){frame[iSample] = double(dataFrame.sample(iSample).adc())-pedestal ;}
00146         else {frame[iSample] = (double(dataFrame.sample(iSample).adc())-pedestals[GainId-1])*gainRatios[GainId-1];}
00147 
00148         if( frame[iSample]>maxsample ) {
00149           maxsample = frame[iSample];
00150           imax = iSample;
00151         }
00152       }
00153   }
00154   else {// pedestal from pre-sample
00155     external_pede = false;
00156     pedestal = (double(dataFrame.sample(0).adc()) + double(dataFrame.sample(1).adc()))/2.;
00157 
00158     for(int iSample = 0; iSample < C::MAXSAMPLES; iSample++) {
00159       //create frame in adc gain 12 equivalent
00160       GainId = dataFrame.sample(iSample).gainId();
00161       //no gain switch forseen if there is no external pedestal
00162       if ( GainId == 0 ) 
00163         {
00164           GainId = 3;
00165           isSaturated = 1;
00166         }
00167 
00168       frame[iSample] = double(dataFrame.sample(iSample).adc())-pedestal ;
00169       // if gain has switched but no pedestals are available, no much good you can do...
00170       if (GainId > gainId0) iGainSwitch = 1;
00171       if( frame[iSample]>maxsample ) {
00172         maxsample = frame[iSample];
00173         imax = iSample;
00174       }
00175     } 
00176   }
00177 
00178   if( (iGainSwitch==1 && external_pede==false) ||  // ... thus you return dummy rechit
00179       imax ==-1 ){                                 // protect against all frames being <-1
00180     return EcalUncalibratedRecHit( dataFrame.id(), -1., -100., -1. , -1.);
00181   }
00182   
00183   InitFitParameters(frame, imax);
00184   chi2_ = PerformAnalyticFit(frame,imax);
00185   //uint32_t flags = 0;
00186   //if (isSaturated) flags = EcalUncalibratedRecHit::kSaturated;
00187 
00188   /*    std::cout << "separate fits\nA: " << fAmp_max_  << ", ResidualPed: " <<  fPed_max_
00189               <<", pedestal: "<<pedestal << ", tPeak " << fTim_max_ << std::endl;
00190   */
00191   return EcalUncalibratedRecHit( dataFrame.id(),fAmp_max_, pedestal+fPed_max_, fTim_max_ - 5 , chi2_);
00192   //return EcalUncalibratedRecHit( dataFrame.id(),fAmp_max_, pedestal+fPed_max_, fTim_max_ - 5 , chi2_, flags );
00193 }

template<class C>
float EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit ( double *  samples,
int  max_sample 
) [inline, private]

fit electronic function from simulation parameters fAlpha_ and fBeta_ are fixed and fit is providing the 3 following parameters the maximum amplitude ( fAmp_max_ ) the time of the maximum ( fTim_max_)

Loop on iterations

initialization inside iteration loop !

Then we loop on samples to be fitted

calculate function to be fitted

then calculate derivatives of function to be fitted

compute matrix elements DM1

compute delta

compute vector elements PROD

end of loop on samples

compute variations of parameters fAmp_max and fTim_max

end of loop on iterations

protection again diverging/unstable fit

results of the fit are calculated

Definition at line 240 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

References EcalUncalibRecHitFixedAlphaBetaAlgo< C >::alfabeta_, archive::db, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::DM1_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::doFit_, funct::exp(), cmsPerfPublish::fail(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAlpha_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAmp_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fBeta_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fNb_iter_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fNum_samp_after_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fNum_samp_bef_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fPed_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fTim_max_, stor::func(), i, i1, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters(), int, edm::isnan(), iter, LogDebug, funct::pow(), PROD, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::temp_, and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::un_sur_sigma.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::makeRecHit().

00240                                                                                                                  {
00241   
00242   //int fValue_tim_max = max_sample;  
00247   //| the pedestal (fPed_max_)  
00248   
00249   double chi2=-1 , db[3] ;
00250   
00251 
00252   //HepSymMatrix DM1(3) ; HepVector temp(3) ;
00253 
00254   int num_fit_min =(int)(max_sample - fNum_samp_bef_max_ ) ;
00255   int num_fit_max =(int)(max_sample + fNum_samp_after_max_) ;
00256 
00257   if (num_fit_min<0) num_fit_min = 0 ; 
00258   //if (num_fit_max>=fNsamples-1) num_fit_max = fNsamples-2 ;
00259   if (num_fit_max>= C::MAXSAMPLES ) {num_fit_max = C::MAXSAMPLES-1 ;}
00260 
00261   if(! doFit_ ) {
00262     LogDebug("EcalUncalibRecHitFixedAlphaBetaAlgo")<<"No fit performed. The amplitude of sample 5 will be used";
00263     return -1;
00264   }
00265 
00266   double func,delta ;
00267   double variation_func_max = 0. ;double variation_tim_max = 0. ; double variation_ped_max = 0. ;
00269   for (int iter=0 ; iter < fNb_iter_ ; iter ++) {
00271     chi2 = 0. ; //PROD.Zero() ;  DM1.Zero() ;
00272 
00273      for(int i1=0 ; i1<3 ; i1++) {
00274        temp_[i1]=0;
00275         for(int j1=i1 ; j1<3 ; j1++) { 
00276           DM1_.fast(j1+1,i1+1) = 0; }
00277       }
00278 
00279     fAmp_max_ += variation_func_max ;
00280     fTim_max_ += variation_tim_max ;
00281     fPed_max_ += variation_ped_max ;
00282     
00284     for( int i = num_fit_min ; i <= num_fit_max ; i++) {
00285       //if(i>fsamp_edge_fit && i<num_fit_min) continue ; // remove front edge samples
00287       func = pulseShapeFunction( (double)i  ) ;
00289       double dt =(double)i - fTim_max_ ;
00290       if(dt > -alfabeta_)  {      
00291         double dt_sur_beta = dt/fBeta_ ;
00292         double variable = (double)1. + dt/alfabeta_ ;
00293         double expo = exp(-dt_sur_beta) ;        
00294         double puissance = pow(variable,fAlpha_) ;
00295         
00296         db[0]=un_sur_sigma*puissance*expo ;
00297         db[1]=fAmp_max_*db[0]*dt_sur_beta/(alfabeta_*variable) ;
00298       }
00299       else {
00300         db[0]=0. ; db[1]=0. ; 
00301       }
00302       db[2]=un_sur_sigma ;
00304       for(int i1=0 ; i1<3 ; i1++) {
00305         for(int j1=i1 ; j1<3 ; j1++) { 
00306           //double & fast(int row, int col);
00307           DM1_.fast(j1+1,i1+1) += db[i1]*db[j1]; }
00308       }
00310       delta = (samples[i]-func)*un_sur_sigma ;
00312       for(int ii=0 ; ii<3 ;ii++) {temp_[ii] += delta*db[ii] ;}
00313       chi2 += delta * delta ;
00314     }
00315     
00316     int fail=0 ;
00317     DM1_.invert(fail) ;
00318       if(fail != 0.) {
00319       //just a guess from the value of the parameters in the previous interaction;
00320       //printf("wH4PulseFitWithFunction =====> determinant error --> No Fit Provided !\n") ;
00321       InitFitParameters(samples,max_sample);
00322       return -101 ;
00323     }
00324 /*     for(int i1=0 ; i1<3 ; i1++) { */
00325 /*       for(int j1=0 ; j1<3 ; j1++) {  */
00326 /*      //double & fast(int row, int col); */
00327 /*      std::cout<<"inverted: "<<DM1[j1][i1]<<std::endl;;} */
00328 /*     } */
00329 /*     std::cout<<"vector temp: "<< temp[0]<<" "<<temp[1]<<" "<<temp[2]<<std::endl; */
00331     HepVector PROD = DM1_*temp_ ;
00332     //    std::cout<<"vector PROD: "<< PROD[0]<<" "<<PROD[1]<<" "<<PROD[2]<<std::endl;
00333 
00334     // Probably the fastest way to protect against
00335     // +-inf value in the matrix DM1_ after inversion
00336     // (which is nevertheless flagged as successfull...)
00337     if ( isnan( PROD[0] ) ) {
00338             InitFitParameters(samples,max_sample);
00339             return -103 ;
00340     }
00341 
00342     variation_func_max = PROD[0] ;
00343     variation_tim_max = PROD[1] ;
00344     variation_ped_max = PROD[2] ;
00345     //chi2 = chi2/((double)nsamp_used - 3.) ;
00346   }
00347   
00348 
00350   if( variation_func_max > 2000. || variation_func_max < -1000. ) {
00351     InitFitParameters(samples,max_sample);
00352     return -102;
00353   }
00354   
00355 
00357   fAmp_max_ += variation_func_max ;
00358   fTim_max_ += variation_tim_max ;
00359   fPed_max_ += variation_ped_max ;
00360 
00361 
00362   // protection against unphysical results:
00363   // ampli mismatched to MaxSample, ampli largely negative, time off preselected range
00364   if( fAmp_max_>2*samples[max_sample]  ||  fAmp_max_<-100 ||  fTim_max_<0  ||  9<fTim_max_ ) {
00365     InitFitParameters(samples,max_sample);
00366     return -104;
00367   }
00368   
00369 
00370   //std::cout <<"chi2: "<<chi2<<" ampl: "<<fAmp_max_<<" time: "<<fTim_max_<<" pede: "<<fPed_max_<<std::endl;
00371   return chi2;
00372 }

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction ( double  t  )  [inline, private]

Definition at line 195 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

References EcalUncalibRecHitFixedAlphaBetaAlgo< C >::alfabeta_, funct::exp(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAlpha_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAmp_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fBeta_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fPed_max_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fTim_max_, and funct::pow().

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

00195                                                                                            {
00196   if( alfabeta_ <= 0 ) return((double)0.);
00197   double dtsbeta,variable,puiss;
00198   double dt = t-fTim_max_ ;
00199   if(dt > -alfabeta_)  {
00200     dtsbeta=dt/fBeta_ ;
00201     variable=1.+dt/alfabeta_ ;
00202     puiss=pow(variable,fAlpha_);
00203     return fAmp_max_*puiss*exp(-dtsbeta)+fPed_max_ ;
00204   }
00205   return  fPed_max_ ;
00206 }

template<class C>
void EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetAlphaBeta ( double  alpha,
double  beta 
) [inline]

Definition at line 374 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

References EcalUncalibRecHitFixedAlphaBetaAlgo< C >::alfabeta_, EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAlpha_, and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fBeta_.

Referenced by EcalFixedAlphaBetaFitUncalibRecHitProducer::produce().

00374                                                                                                      {
00375   fAlpha_ = alpha;
00376   fBeta_=  beta;
00377   alfabeta_ = fAlpha_*fBeta_;
00378 }

template<class C>
void EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetDynamicPedestal ( bool  dyn_pede  )  [inline]

Definition at line 383 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

References EcalUncalibRecHitFixedAlphaBetaAlgo< C >::dyn_pedestal.

Referenced by EcalFixedAlphaBetaFitUncalibRecHitProducer::EcalFixedAlphaBetaFitUncalibRecHitProducer().

00383                                                                                         {
00384   dyn_pedestal = p;
00385 }

template<class C>
void EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetMinAmpl ( double  ampl  )  [inline]

Definition at line 380 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

References EcalUncalibRecHitFixedAlphaBetaAlgo< C >::MinAmpl_.

Referenced by EcalFixedAlphaBetaFitUncalibRecHitProducer::EcalFixedAlphaBetaFitUncalibRecHitProducer().

00380                                                                                      {
00381   MinAmpl_ = ampl;
00382 }

Member Data Documentation

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::alfabeta_ [private]

Definition at line 43 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetAlphaBeta().

template<class C>
float EcalUncalibRecHitFixedAlphaBetaAlgo< C >::alpha_table_[36][1701] [private]

Definition at line 53 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

template<class C>
float EcalUncalibRecHitFixedAlphaBetaAlgo< C >::beta_table_[36][1701] [private]

Definition at line 54 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

template<class C>
HepSymMatrix EcalUncalibRecHitFixedAlphaBetaAlgo< C >::DM1_ [private]

Definition at line 61 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

template<class C>
bool EcalUncalibRecHitFixedAlphaBetaAlgo< C >::doFit_ [private]

Definition at line 56 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

template<class C>
bool EcalUncalibRecHitFixedAlphaBetaAlgo< C >::dyn_pedestal [private]

Definition at line 37 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::makeRecHit(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetDynamicPedestal().

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAlpha_ [private]

Definition at line 38 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetAlphaBeta().

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fAmp_max_ [private]

Definition at line 40 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::makeRecHit(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction().

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fBeta_ [private]

Definition at line 39 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetAlphaBeta().

template<class C>
int EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fNb_iter_ [private]

Definition at line 46 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

template<class C>
int EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fNum_samp_after_max_ [private]

Definition at line 48 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

template<class C>
int EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fNum_samp_bef_max_ [private]

Definition at line 47 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fPed_max_ [private]

Definition at line 42 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::makeRecHit(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction().

template<class C>
float EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fSigma_ped [private]

Definition at line 50 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::fTim_max_ [private]

Definition at line 41 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::makeRecHit(), EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::pulseShapeFunction().

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::MinAmpl_ [private]

Definition at line 36 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::InitFitParameters(), and EcalUncalibRecHitFixedAlphaBetaAlgo< C >::SetMinAmpl().

template<class C>
HepVector EcalUncalibRecHitFixedAlphaBetaAlgo< C >::temp_ [private]

Definition at line 61 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

template<class C>
double EcalUncalibRecHitFixedAlphaBetaAlgo< C >::un_sur_sigma [private]

Definition at line 51 of file EcalUncalibRecHitFixedAlphaBetaAlgo.h.

Referenced by EcalUncalibRecHitFixedAlphaBetaAlgo< C >::PerformAnalyticFit().

The documentation for this class was generated from the following file:
Generated on Tue Jun 9 18:20:03 2009 for CMSSW by  doxygen 1.5.4