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	ZDCPrePulseFitWrapper
class	ZDCPreExpFitWrapper
class	ZDCFitExpFermiVariableTaus
class	ZDCFitExpFermiVariableTausLHCf
class	ZDCFitExpFermiVariableTausInduct
class	ZDCFitExpFermiVariableTausRun3
class	ZDCFitExpFermiFixedTaus
class	ZDCFitExpFermiPrePulse
class	ZDCFitExpFermiLHCfPrePulse
class	ZDCFitExpFermiPreExp
class	ZDCFitExpFermiLHCfPreExp
class	ZDCFitExpFermiInductPreExp
class	ZDCFitExpFermiLinearFixedTaus
class	ZDCFitExpFermiLinearPrePulse
class	ZDCFitComplexPrePulse
class	ZDCFitGeneralPulse

Functions
double	ZDCFermiExpFit (const double *xvec, const double *pvec)
double	ZDCFermiExpFitRefl (const double *xvec, const double *pvec)
double	ZDCFermiExpFitInduct (const double *xvec, const double *pvec)

Function Documentation

ZDCFermiExpFit()

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

inline

Definition at line 1390 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 1451 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 delta = pvec[8];
  
  double tauRatio = tau2 / tau1;
  double tauRatioMinunsOne = tauRatio - 1;
 
  double norm = std::pow(1. / tauRatioMinunsOne, 1. / (1.0 + tauRatio)) /
    ( 1.0 + std::pow(1. / tauRatioMinunsOne, 1. / (1.0 + 1.0 / tauRatio))) ;
 
  double deltaT = t - t0;
  double deltaTInduct = deltaT - tau1 * std::log(tauRatioMinunsOne);
  
  if (deltaT < 0) deltaT = 0;
  if (deltaTInduct < 0) deltaTInduct = 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 twoPiOverPeriod = 2.0*M_PI/period;
  double inductTerm = (1.0 + Acos*std::cos(deltaTInduct*twoPiOverPeriod) +
               Bsin*std::sin(deltaTInduct*twoPiOverPeriod))/(1+Acos);
  double expTerm = delta + std::exp(-deltaT / tau2)*inductTerm;
  double fermiTerm = 1. / (1. + std::exp(-(t - t0) / tau1));
 
  return amp * expTerm * fermiTerm / norm + C; 
}

ZDCFermiExpFitRefl()

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

inline

Definition at line 1415 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; 
}