ZDCFitWrapper.h File Reference

#include "CxxUtils/checker_macros.h"
#include <TF1.h>
#include <memory>
#include <cmath>
#include <stdexcept>
#include <algorithm>

Include dependency graph for ZDCFitWrapper.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	ZDCFitWrapper
class	ZDCFitWrapper::FitParameter
class	ZDCPrePulseFitWrapper
class	ZDCPreExpFitWrapper
class	ZDCFitExpFermiVariableTaus
class	ZDCFitExpFermiVariableTausLHCf
class	ZDCFitExpFermiVariableTausInduct
class	ZDCFitExpFermiVariableTausRun3
class	ZDCFitExpFermiFixedTaus
class	ZDCFitExpFermiPrePulse
class	ZDCFitExpFermiLHCfPrePulse
class	ZDCFitExpFermiPreExp
class	ZDCFitExpFermiLHCfPreExp
class	ZDCFitExpFermiInductPreExp
class	ZDCFitExpFermiInductPrePulse
class	ZDCFitExpFermiLinearFixedTaus
class	ZDCFitExpFermiLinearPrePulse
class	ZDCFitComplexPrePulse
class	ZDCFitGeneralPulse

Functions
double	ZDCFermiExpFit (const double *xvec, const double *pvecpvec)
double	ZDCFermiExpFitRefl (const double *xvec, const double *pvecpvec)
double	ZDCFermiExpFitInduct (const double *xvec, const double *pvecpvec)
double	getInductFuncTimeMax (const TF1 *theTF1)

Function Documentation

getInductFuncTimeMax()

double getInductFuncTimeMax ( const TF1 * theTF1 )

inline

Definition at line 1681 of file ZDCFitWrapper.h.

{
  float fitT0 = theTF1->GetParameter(1);
  
  float tau1 = theTF1->GetParameter(2);
  float tau2 = theTF1->GetParameter(3);
  
  // Correct the time to the maximum
  //
  double corrT0 = fitT0;
  if (tau2 > tau1) corrT0 += tau1 * std::log(tau2 / tau1 - 1.0);
  
  double deltaT = theTF1->GetParameter(5)/8;
  double tmin = corrT0 - deltaT;
  double tmax = corrT0 + deltaT;
  
  return theTF1->GetMaximumX(tmin, tmax, 1e-4, 20, false);
}

ZDCFermiExpFit()

double ZDCFermiExpFit ( const double * xvec, const double * pvec )

inline

Definition at line 1701 of file ZDCFitWrapper.h.

{
  double t = xvec[0];
 
  double amp = pvec[0];
  double t0 = pvec[1];
  double tau1 = pvec[2];
  double tau2 = pvec[3];
  double C = pvec[4];
 
  double tauRatio = tau2 / tau1;
  double tauRatioMinunsOne = tauRatio - 1;
 
  double norm = (std::pow(1. / tauRatioMinunsOne, 1. / (1 + tauRatio)) /
         ( 1 + std::pow(1. / tauRatioMinunsOne, 1. / (1 + 1 / tauRatio))));
 
  double deltaT = t - t0;
  if (deltaT < 0) deltaT = 0;
 
  double expTerm = std::exp(-deltaT / tau2);
  double fermiTerm = 1. / (1. + std::exp(-(t - t0) / tau1));
 
  return amp * expTerm * fermiTerm / norm + C; 
}

ZDCFermiExpFitInduct()

double ZDCFermiExpFitInduct ( const double * xvec, const double * pvec )

inline

Definition at line 1762 of file ZDCFitWrapper.h.

{
  double t = xvec[0];
 
  double amp = pvec[0];
  double t0 = pvec[1];
  double tau1 = pvec[2];
  double tau2 = pvec[3];
  double C = pvec[4];
  
  double period = pvec[5];
  double Acos = pvec[6];
  double Bsin = pvec[7];
  double bI = pvec[8];
  double delta = pvec[9];
 
  double twoPiOverPeriod = 2.0*M_PI/period;
 
  double deltaT = t - t0;
  double norm = 1, expTerm = 1, fermiTerm = 1, inductTerm = 1;
  
  if (deltaT > 0) {
    inductTerm = (1.0 + exp(-bI*deltaT)*(Acos*std::cos(deltaT*twoPiOverPeriod) +
                       Bsin*std::sin(deltaT*twoPiOverPeriod)))/(1+Acos);
    expTerm = delta + std::exp(-deltaT / tau2)*inductTerm;
  }
  if (deltaT/tau1 < 7) fermiTerm = 1. / (1. + std::exp(-deltaT / tau1));
 
  return amp * expTerm * fermiTerm * inductTerm/ norm + C; 
}

ZDCFermiExpFitRefl()

double ZDCFermiExpFitRefl ( const double * xvec, const double * pvec )

inline

Definition at line 1726 of file ZDCFitWrapper.h.

{
  double t = xvec[0];
 
  double amp = pvec[0];
  double t0 = pvec[1];
  double tau1 = pvec[2];
  double tau2 = pvec[3];
  double C = pvec[4];
  
  double refldelay = pvec[5];
  double reflFrac = pvec[6];
  double reflwidth = pvec[7];
  double delta = pvec[8];
  
  double tauRatio = tau2 / tau1;
  double tauRatioMinunsOne = tauRatio - 1;
 
  double norm = std::pow(1. / tauRatioMinunsOne, 1. / (1 + tauRatio)) /
    ( 1 + std::pow(1. / tauRatioMinunsOne, 1. / (1 + 1 / tauRatio))) ;
 
  double deltaT = t - t0;
  if (deltaT < 0) deltaT = 0;
 
  // Note: the small constant added here accounts for the very long tail on the pulse 
  //   which doesn't go to zero over the time range that we sample
  //
  double expTerm = delta + std::exp(-deltaT / tau2);
  double fermiTerm = 1. / (1. + std::exp(-(t - t0) / tau1));
 
  double deltaTRefl = deltaT - refldelay;
  double reflTerm =  -reflFrac*amp*std::exp(-0.5*deltaTRefl*deltaTRefl/reflwidth/reflwidth);
 
  return amp * expTerm * fermiTerm / norm + C + reflTerm; 
}