---

CSCFindPeakTime Class Reference

This is CSCFindPeakTime. More...

#include <RecoLocalMuon/CSCRecHitD/src/CSCFindPeakTime.h>

List of all members.

Public Member Functions
	CSCFindPeakTime ()
bool	FindPeakTime (const int &tmax, const float *adc, float &t_zero, float &t_peak)
	Member functions.
void	FitCharge (const int &tmax, const float *adc, const float &t_zero, const float &t_peak, std::vector< float > &adcsFit)
	Fitting the charge for a given t_zero and t_peak.
	~CSCFindPeakTime ()

---

Detailed Description

This is CSCFindPeakTime.

Author:

Dominique Fortin

Based on RecoLocalMuon/CSCStandAlone/interface/PulseTime.h by S. Durkin

Fast fit adc values to
N*(p0**2/256/exp(-4))*(t-t0)**4*exp(-p0*(t-t0))

PulseTime has peaking time fixed to 133 nsec.
It actually fits for the t0

Fit the log of data
chi2=(log(Data)-log(N)-4*log(t-t0)+p0*(t-t0))**2

Becomes least square fit in p0 and log(N); and binary search in t0
yields N, p0, t0(nsec)

Note: tpeak=4/p0 (nsec) and adc[0] is arbitrarily defined a time of 0.0 nsec.

Finally, a fit to the charge deposition for each time bin is performed using
the fitted t0.

Definition at line 28 of file CSCFindPeakTime.h.

---

Constructor & Destructor Documentation

CSCFindPeakTime::CSCFindPeakTime

(

)

[inline]

Definition at line 32 of file CSCFindPeakTime.h.

{}; 

CSCFindPeakTime::~CSCFindPeakTime

(

)

[inline]

Definition at line 34 of file CSCFindPeakTime.h.

{}; 

---

Member Function Documentation

bool CSCFindPeakTime::FindPeakTime	(	const int &	tmax,
		const float *	adc,
		float &	t_zero,
		float &	t_peak
	)

Member functions.

Finding the peak time and zeroth time

Definition at line 22 of file CSCFindPeakTime.cc.

References funct::exp(), j, N, and x.

Referenced by CSCXonStrip_MatchGatti::findXOnStrip().

                                                                                                    {
  
  // Initialize parameters in case fit fails
  float t0       = 0.;
  float N        = adc[1];
  t_peak         = 133.;
  float p0       = 4./t_peak;

  // If outside physical range, exit
  if ( tmax < 2 || tmax > 6 ) return false;

  float tb[4];
  for ( int time=0; time<4; ++time ){
    tb[time] = (tmax + time -1) * 50.;
  }

  int n_fit  = 4;
  if ( tmax == 6 ) n_fit = 3;

  float chi_min  = 1.e10;
  float chi_last = 1.e10;
  float tt0      = 0.;
  float chi2     = 0.;
  float del_t    = 100.;

  float x[4];
  float sx2 = 0.;
  float sxy = 0.;
  float NN = 0.;

  while ( del_t > 1. ) {
    sx2 = 0.;
    sxy = 0.;
        
    for ( int j=0; j < n_fit; ++j ) {
      x[j] = (tb[j] - tt0) * (tb[j] - tt0) * (tb[j] - tt0) * (tb[j] - tt0) * exp( -p0 * (tb[j] - tt0) );
      sx2  += x[j] * x[j];
      sxy  += x[j] * adc[j];
    }
    NN = sxy / sx2;
    
    // Compute chi^2
    chi2 = 0.0;
    for (int j=0; j < n_fit; ++j) chi2 += (adc[j] - NN * x[j]) * (adc[j] - NN * x[j]);

    // Test on chi^2 to decide what to do    
    if ( chi_last > chi2 ) {
      if (chi2 < chi_min ){
        t0      = tt0;
        N       = NN;
      }
      chi_last  = chi2;
      tt0       = tt0 + del_t;
    } else {
      tt0      = tt0 - 2. * del_t;
      del_t    = del_t / 2.;
      tt0      = tt0 + del_t;
      chi_last = 1.0e10;
    }
  }

  t_peak = t_peak;
  t_zero = tt0;

  return true;
}

void CSCFindPeakTime::FitCharge	(	const int &	tmax,
		const float *	adc,
		const float &	t_zero,
		const float &	t_peak,
		std::vector< float > &	adcsFit
	)

Fitting the charge for a given t_zero and t_peak.

Definition at line 93 of file CSCFindPeakTime.cc.

References funct::exp(), i, j, N, t, x, and y.

                                                                                                                                        {

  float p0  = 4./t_peak;
  float tt0 = t_zero;
  int n_fit = 4;
  if ( tmax == 6 ) n_fit=3;
  
  float tb[4], y[4];
  for ( int t = 0; t < 4; ++t ){
    tb[t] = (tmax + t - 1) * 50.;
    y[t] = adc[t];
  }

  // Find the normalization factor for the function
  float x[4];    
  float sx2 = 0.;
  float sxy = 0.;
  for ( int j=0; j < n_fit; ++j ) {
    float t = tb[j];
    x[j] = (t-tt0)*(t-tt0)*(t-tt0)*(t-tt0) * exp( -p0 * (t-tt0) );
    sx2  = sx2 + x[j] * x[j];
    sxy  = sxy + x[j] * y[j];
  }
  float N = sxy / sx2;
    

  // Now compute charge for a given t  --> only need charges at: tpeak-50, tpeak and tpeak+50
  for ( int i = 0; i < 3; ++i ) {
    float t = t_peak + (i - 1) * 50.;
    float q_fitted = N * (t-tt0)*(t-tt0)*(t-tt0)*(t-tt0) * exp( -p0 * (t-tt0) );
    adcsFit.push_back(q_fitted);
  }
  return;
}

---

The documentation for this class was generated from the following files:

• RecoLocalMuon/CSCRecHitD/src/CSCFindPeakTime.h
• RecoLocalMuon/CSCRecHitD/src/CSCFindPeakTime.cc

---