#include <ZDCFitWrapper.h>

Inheritance diagram for ZDCFitExpFermiInductPreExp:

Collaboration diagram for ZDCFitExpFermiInductPreExp:

Public Member Functions
	ZDCFitExpFermiInductPreExp (const std::string &tag, float tmin, float tmax, float tau1, float tau2, float defExpTau, float fixExpTau)
	ZDCFitExpFermiInductPreExp (const std::string &tag, const std::string &overRideParams, float tmin, float tmax, float tau1, float tau2)
	~ZDCFitExpFermiInductPreExp ()
virtual void	DoInitialize (float initialAmp, float initialT0, float ampMin, float ampMax) override
virtual void	SetT0FitLimits (float tMin, float tMax) override
virtual void	SetInitialExpPulse (float amp) override
virtual void	doConstrainFit () override
virtual void	doUnconstrainFit () override
virtual float	GetAmplitude () const override
virtual float	GetAmpError () const override
virtual float	GetTau1 () const override
virtual float	GetTau2 () const override
virtual float	GetTime () const override
virtual float	GetExpAmp () const override
virtual float	GetExpTau () const override
virtual unsigned int	GetNumShapeParameters () const override
virtual float	GetShapeParameter (size_t index) const override
virtual float	GetBkgdMaxFraction () const override
virtual double	operator() (const double x, const double p) override
bool	fixExpTau () const
float	getDefaultExpTau () const
bool	HaveParameterSet () const
void	OverrideParameterSet (const std::vector< FitParameter > &parameterSet)
const FitParameter &	GetParameterDescr (const std::string &name)
void	Initialize (float initialAmp, float initialT0, float ampMin, float ampMax)
void	Initialize (float initialAmp, float initialT0, float ampMin, float ampMax, float fitTmin, float fitTmax, float fitTRef)
bool	finalized () const
virtual void	Finalize ()
void	SetAmpMinMax (float minAmp, float maxAmp)
void	SetT0Range (float t0Min, float t0Max)
void	ConstrainFit ()
void	UnconstrainFit ()
float	GetMinAmp () const
float	GetMaxAmp () const
float	GetTMin () const
float	GetTMax () const
float	GetT0Min () const
float	GetT0Max () const
float	GetTMinAdjust () const
virtual std::shared_ptr< TF1 >	GetWrapperTF1 ()
virtual const TF1 *	GetWrapperTF1 () const
virtual TF1 *	GetWrapperTF1RawPtr () const

Protected Member Functions
void	LoadDefaultParameterSet (const std::vector< FitParameter > &parameterSet)

Protected Attributes
float	m_defaultTau
bool	m_fixTau

Private Member Functions
virtual void	DoFinalize () override

Private Attributes
float	m_tau1 {0}
float	m_tau2 {0}
float	m_tFuncMax {-999}
float	m_funcMax {-999}
std::shared_ptr< TF1 >	m_expFermiInductFunc {0}
std::shared_ptr< TF1 >	m_expFermiPreFunc {0}
std::shared_ptr< TF1 >	m_wrapperTF1 {}
float	m_tmin {0}
float	m_tmax {0}
float	m_ampMin {0}
float	m_ampMax {0}
float	m_t0Min {0}
float	m_t0Max {0}
bool	m_adjTLimitsEvent {false}
float	m_tminAdjust {0}
float	m_tempTmin {0}
float	m_tempTmax {0}
bool	m_haveParameterSet {false}
std::vector< FitParameter >	m_parameterSet
bool	m_finalized {false}

Static Private Attributes
static const std::map< std::string, std::vector< ZDCFitWrapper::FitParameter > >	s_ConfigParamMap

Detailed Description

Definition at line 972 of file ZDCFitWrapper.h.

Constructor & Destructor Documentation

◆ ZDCFitExpFermiInductPreExp() [1/2]

ZDCFitExpFermiInductPreExp::ZDCFitExpFermiInductPreExp	(	const std::string &	tag,
		float	tmin,
		float	tmax,
		float	tau1,
		float	tau2,
		float	defExpTau,
		float	fixExpTau )

Definition at line 1081 of file ZDCFitWrapper.cxx.

                                                             :
  ZDCPreExpFitWrapper(std::make_shared<TF1>(("ExpFermiInductPreExp" + tag).c_str(), this, tmin, tmax, 10), defExpTau, fixExpTau),
  m_tau1(tau1), m_tau2(tau2)
{
  // Create the reference function that we use to evaluate ExpFermiFit more efficiently
  //
  std::string funcNameRefFunc = "ExpFermiInductPreExpRefFunc" + tag;
 
  m_expFermiInductFunc = std::make_shared<TF1>(funcNameRefFunc.c_str(), ZDCFermiExpFitInduct, tmin, tmax, 10);
 
  // The parameters for the FermiExpInduct
  //
  m_expFermiInductFunc->SetParameter(0, 1);
  m_expFermiInductFunc->SetParameter(1, 0);
  m_expFermiInductFunc->SetParameter(2, m_tau1); // tau1 is fixed
  m_expFermiInductFunc->SetParameter(3, m_tau2); // tau2 is fixed
  m_expFermiInductFunc->SetParameter(4, 0);      // the constant term will be handled directly not here
  m_expFermiInductFunc->SetParameter(5, 14); // the period will be fixed
  m_expFermiInductFunc->SetParameter(6, 0.25);  // A will be varied as part of the fit
  m_expFermiInductFunc->SetParameter(7, 0.5);  // B will be varied as part of the fit
  m_expFermiInductFunc->SetParameter(8, 0.18);
  m_expFermiInductFunc->SetParameter(9, 0.015);
 
  // Now set up the actual TF1
  //
  std::shared_ptr<TF1> theTF1 = ZDCFitWrapper::GetWrapperTF1();
  theTF1->SetParName(0, "Amp");
  theTF1->SetParName(1, "T0");
  theTF1->SetParName(2, "Aexp");
  theTF1->SetParName(3, "tauexp");
  theTF1->SetParName(4, "bsqrtexp");
  theTF1->SetParName(5, "C");
  theTF1->SetParName(6, "Period");
  theTF1->SetParName(7, "IndA");
  theTF1->SetParName(8, "IndB");
  theTF1->SetParName(9, "bI");
 
  theTF1->SetParLimits(1, tmin, tmax);
  theTF1->SetParLimits(2, -1, 8196); // Increase the upper range to 2 times of ADC range to deal with large exponential tail case of pre-pulse.
  theTF1->SetParLimits(3, 4, 12);
  theTF1->SetParLimits(4, -0.001, 0.001);
  theTF1->SetParLimits(5, -10, 10);
  theTF1->FixParameter(6, 13);
  theTF1->SetParLimits(7, -0.1, 2);
  theTF1->SetParLimits(8, -0.1, 2);
  theTF1->SetParLimits(9, -0.05, 0.5);
  
}

◆ ZDCFitExpFermiInductPreExp() [2/2]

ZDCFitExpFermiInductPreExp::ZDCFitExpFermiInductPreExp	(	const std::string &	tag,
		const std::string &	overRideParams,
		float	tmin,
		float	tmax,
		float	tau1,
		float	tau2 )

Definition at line 1131 of file ZDCFitWrapper.cxx.

                                                                                                                                                       :
  ZDCPreExpFitWrapper(std::make_shared<TF1>(("ExpFermiInductPreExp" + tag).c_str(), this, tmin, tmax, 10), s_ConfigParamMap.at("default")),
  m_tau1(tau1), m_tau2(tau2)
{
  auto iter = s_ConfigParamMap.find(paramSet);
  if (iter == s_ConfigParamMap.end()) throw std::runtime_error("ZDCFitExpFermiVariableTausInduct:: invalid parameter set identifier: " + paramSet);
 
  std::vector<ZDCFitWrapper::FitParameter> overRideParams = iter->second;
 
  OverrideParameterSet(overRideParams);
 
  // Create the reference function that we use to evaluate ExpFermiFit more efficiently
  //
  std::string funcNameRefFunc = "ExpFermiInductPreExpRefFunc" + tag;
 
  m_expFermiInductFunc = std::make_shared<TF1>(funcNameRefFunc.c_str(), ZDCFermiExpFitInduct, tmin, tmax, 10);
 
  // The parameters for the FermiExpInduct
  //
  m_expFermiInductFunc->SetParameter(0, 1);
  m_expFermiInductFunc->SetParameter(1, 0);
  m_expFermiInductFunc->SetParameter(2, m_tau1); // tau1 is fixed
  m_expFermiInductFunc->SetParameter(3, m_tau2); // tau2 is fixed
  m_expFermiInductFunc->SetParameter(4, 0);      // the constant term will be handled directly not here
  m_expFermiInductFunc->SetParameter(5, 14); // the period will be fixed
  m_expFermiInductFunc->SetParameter(6, 0.25);  // A will be varied as part of the fit
  m_expFermiInductFunc->SetParameter(7, 0.5);  // B will be varied as part of the fit
  m_expFermiInductFunc->SetParameter(8, 0.18);
  m_expFermiInductFunc->SetParameter(9, 0.015);
}

◆ ~ZDCFitExpFermiInductPreExp()

ZDCFitExpFermiInductPreExp::~ZDCFitExpFermiInductPreExp ( )

inline

Definition at line 1005 of file ZDCFitWrapper.h.

1005{}

Member Function Documentation

◆ ConstrainFit()

void ZDCFitWrapper::ConstrainFit ( )

inherited

Definition at line 277 of file ZDCFitWrapper.cxx.

{
  if (m_haveParameterSet) {  
    auto tf1_p = GetWrapperTF1(); 
 
    for (unsigned int ipar = 2; ipar < m_parameterSet.size(); ipar++) {
      const auto & fitparDescr = m_parameterSet[ipar];
      
      if (!fitparDescr.fixed() && fitparDescr.constrainFixed()) {
    tf1_p->FixParameter(ipar, fitparDescr.constrainValue());
      }
    }
  }
  else doConstrainFit();
}

◆ doConstrainFit()

void ZDCFitExpFermiInductPreExp::doConstrainFit ( )

overridevirtual

Implements ZDCFitWrapper.

Definition at line 1195 of file ZDCFitWrapper.cxx.

1196{

1197}

◆ DoFinalize()

virtual void ZDCFitExpFermiInductPreExp::DoFinalize ( )

inlineoverrideprivatevirtual

Reimplemented from ZDCFitWrapper.

Definition at line 985 of file ZDCFitWrapper.h.

  {
    m_expFermiInductFunc->SetNpx(20);
    double t0 = GetWrapperTF1()->GetParameter(1);
    double tsubtr = getInductFuncTimeMax(m_expFermiInductFunc.get());
    if (std::isnan(tsubtr)) {
      m_tFuncMax = -100; 
      m_funcMax = 0;
    }
    else {
      m_tFuncMax = tsubtr + t0;
      m_funcMax = m_expFermiInductFunc->Eval(tsubtr);
    }
  }

◆ DoInitialize()

void ZDCFitExpFermiInductPreExp::DoInitialize	(	float	initialAmp,
		float	initialT0,
		float	ampMin,
		float	ampMax )

overridevirtual

Implements ZDCFitWrapper.

Definition at line 1162 of file ZDCFitWrapper.cxx.

{
  GetWrapperTF1()->SetParameter(0, std::max(initialAmp, ampMin));
 
  float t0 = initialT0;
  if (t0 < GetT0Min()) t0 = GetT0Min()*1.1;
  if (t0 > GetT0Max()) t0 = GetT0Max()/1.1;
  GetWrapperTF1()->SetParameter(1, t0);
 
  GetWrapperTF1()->SetParameter(2, 1);
  GetWrapperTF1()->SetParameter(3, std::max(getDefaultExpTau(), (float) 6.01));
  GetWrapperTF1()->SetParameter(6, 14);
  GetWrapperTF1()->SetParameter(7, 0.2);
  GetWrapperTF1()->SetParameter(8, 0.2);
  GetWrapperTF1()->SetParameter(9, 0.01);
 
  GetWrapperTF1()->SetParLimits(0, ampMin, ampMax);
 
  // Set parameter errors for fit step size
  //
  double ampStep = std::min(0.05*initialAmp, std::abs(ampMax - initialAmp)/2.);
  GetWrapperTF1()->SetParError(0, ampStep);
  GetWrapperTF1()->SetParError(1, 1.0);
  GetWrapperTF1()->SetParError(2, 5.0);
  GetWrapperTF1()->SetParError(3, 0.5);
  GetWrapperTF1()->SetParError(4, 0.01);
  GetWrapperTF1()->SetParError(5, 1);
  GetWrapperTF1()->SetParError(6, 1);
  GetWrapperTF1()->SetParError(7, 0.05);
  GetWrapperTF1()->SetParError(8, 0.05);
  GetWrapperTF1()->SetParError(9, 0.01);
}

◆ doUnconstrainFit()

void ZDCFitExpFermiInductPreExp::doUnconstrainFit ( )

overridevirtual

Implements ZDCFitWrapper.

Definition at line 1198 of file ZDCFitWrapper.cxx.

1199{

1200}

◆ Finalize()

virtual void ZDCFitWrapper::Finalize ( )

inlinevirtualinherited

Definition at line 136 of file ZDCFitWrapper.h.

  {
    DoFinalize();
    m_finalized = true;
  }

◆ finalized()

bool ZDCFitWrapper::finalized ( ) const

inlineinherited

Definition at line 134 of file ZDCFitWrapper.h.

134{return m_finalized;}

◆ fixExpTau()

bool ZDCPreExpFitWrapper::fixExpTau ( ) const

inlineinherited

Definition at line 256 of file ZDCFitWrapper.h.

256{return m_fixTau;}

ZDCPreExpFitWrapper::m_fixTau

bool m_fixTau

Definition ZDCFitWrapper.h:239

◆ GetAmpError()

virtual float ZDCFitExpFermiInductPreExp::GetAmpError ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1023 of file ZDCFitWrapper.h.

1023{return GetWrapperTF1()->GetParError(0); }

◆ GetAmplitude()

virtual float ZDCFitExpFermiInductPreExp::GetAmplitude ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1018 of file ZDCFitWrapper.h.

                                              {
    if (finalized()) return m_funcMax;
    return m_expFermiInductFunc->Eval(getInductFuncTimeMax(m_expFermiInductFunc.get()));
  }

◆ GetBkgdMaxFraction()

virtual float ZDCFitExpFermiInductPreExp::GetBkgdMaxFraction ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1050 of file ZDCFitWrapper.h.

  {
    const TF1* theTF1 = ZDCFitWrapper::GetWrapperTF1();
    double amp = theTF1->GetParameter(0);
    if (amp < 1e-6) return -1;
 
    double maxTime = GetTime();
    
    double expAmp = theTF1->GetParameter(2);
    double expTau = theTF1->GetParameter(3);
    double bexpSqrt = theTF1->GetParameter(4);
 
    double tRef = GetTMinAdjust();
    double expPre = expAmp * (std::exp(-maxTime/expTau - bexpSqrt*maxTime*maxTime) - std::exp(-tRef/expTau - bexpSqrt*tRef*tRef));
 
    return expPre / (amp + expPre);
  }

◆ getDefaultExpTau()

float ZDCPreExpFitWrapper::getDefaultExpTau ( ) const

inlineinherited

Definition at line 257 of file ZDCFitWrapper.h.

257{return m_defaultTau;}

ZDCPreExpFitWrapper::m_defaultTau

float m_defaultTau

Definition ZDCFitWrapper.h:238

◆ GetExpAmp()

virtual float ZDCFitExpFermiInductPreExp::GetExpAmp ( ) const

inlineoverridevirtual

Implements ZDCPreExpFitWrapper.

Definition at line 1035 of file ZDCFitWrapper.h.

1035{return GetWrapperTF1()->GetParameter(2);}

◆ GetExpTau()

virtual float ZDCFitExpFermiInductPreExp::GetExpTau ( ) const

inlineoverridevirtual

Implements ZDCPreExpFitWrapper.

Definition at line 1036 of file ZDCFitWrapper.h.

1036{return GetWrapperTF1()->GetParameter(3);}

◆ GetMaxAmp()

float ZDCFitWrapper::GetMaxAmp ( ) const

inlineinherited

Definition at line 172 of file ZDCFitWrapper.h.

172{return m_ampMax;}

ZDCFitWrapper::m_ampMax

float m_ampMax

Definition ZDCFitWrapper.h:69

◆ GetMinAmp()

float ZDCFitWrapper::GetMinAmp ( ) const

inlineinherited

Definition at line 171 of file ZDCFitWrapper.h.

171{return m_ampMin;}

ZDCFitWrapper::m_ampMin

float m_ampMin

Definition ZDCFitWrapper.h:68

◆ GetNumShapeParameters()

virtual unsigned int ZDCFitExpFermiInductPreExp::GetNumShapeParameters ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1038 of file ZDCFitWrapper.h.

1038{return 3;}

◆ GetParameterDescr()

const ZDCFitWrapper::FitParameter & ZDCFitWrapper::GetParameterDescr ( const std::string & name )

inherited

Definition at line 118 of file ZDCFitWrapper.cxx.

{
  auto iter = std::find_if(m_parameterSet.begin(), m_parameterSet.end(),
               [name](const FitParameter& parDesc){return parDesc.name() == name;});
  if (iter == m_parameterSet.end()) throw std::logic_error("ZDCFitWrapper::GetParameterDescr: invalid parameter name: " + name );
  return *iter;
}

◆ GetShapeParameter()

virtual float ZDCFitExpFermiInductPreExp::GetShapeParameter ( size_t index ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1040 of file ZDCFitWrapper.h.

  {
    if (index < 3) {
      if (index == 0) return GetWrapperTF1()->GetParameter(4);
      if (index == 1) return GetWrapperTF1()->GetParameter(6);
      if (index == 2) return GetWrapperTF1()->GetParameter(7);
    }
    throw std::runtime_error("Fit parameter does not exist.");
  }

◆ GetT0Max()

float ZDCFitWrapper::GetT0Max ( ) const

inlineinherited

Definition at line 182 of file ZDCFitWrapper.h.

                         {
    if (m_adjTLimitsEvent) return m_tempTmax;
    else return m_t0Max;
  }

◆ GetT0Min()

float ZDCFitWrapper::GetT0Min ( ) const

inlineinherited

Definition at line 177 of file ZDCFitWrapper.h.

                         {
    if (m_adjTLimitsEvent) return m_tempTmin;
    else return m_t0Min;
  }

◆ GetTau1()

virtual float ZDCFitExpFermiInductPreExp::GetTau1 ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1025 of file ZDCFitWrapper.h.

1025{return m_tau1;}

◆ GetTau2()

virtual float ZDCFitExpFermiInductPreExp::GetTau2 ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1026 of file ZDCFitWrapper.h.

1026{return m_tau2;}

◆ GetTime()

virtual float ZDCFitExpFermiInductPreExp::GetTime ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1028 of file ZDCFitWrapper.h.

                                         {
    if (finalized()) return m_tFuncMax;
    double timeMax = getInductFuncTimeMax(m_expFermiInductFunc.get());
    if (std::isnan(timeMax)) return -100;
    else return + GetWrapperTF1()->GetParameter(1);
  }

◆ GetTMax()

float ZDCFitWrapper::GetTMax ( ) const

inlineinherited

Definition at line 175 of file ZDCFitWrapper.h.

175{return m_tmax;}

ZDCFitWrapper::m_tmax

float m_tmax

Definition ZDCFitWrapper.h:66

◆ GetTMin()

float ZDCFitWrapper::GetTMin ( ) const

inlineinherited

Definition at line 174 of file ZDCFitWrapper.h.

174{return m_tmin;}

ZDCFitWrapper::m_tmin

float m_tmin

Definition ZDCFitWrapper.h:65

◆ GetTMinAdjust()

float ZDCFitWrapper::GetTMinAdjust ( ) const

inlineinherited

Definition at line 187 of file ZDCFitWrapper.h.

187{return m_tminAdjust;}

ZDCFitWrapper::m_tminAdjust

float m_tminAdjust

Definition ZDCFitWrapper.h:75

◆ GetWrapperTF1() [1/2]

virtual std::shared_ptr< TF1 > ZDCFitWrapper::GetWrapperTF1 ( )

inlinevirtualinherited

Definition at line 196 of file ZDCFitWrapper.h.

196{return m_wrapperTF1;}

ZDCFitWrapper::m_wrapperTF1

std::shared_ptr< TF1 > m_wrapperTF1

Definition ZDCFitWrapper.h:63

◆ GetWrapperTF1() [2/2]

virtual const TF1 * ZDCFitWrapper::GetWrapperTF1 ( ) const

inlinevirtualinherited

Definition at line 197 of file ZDCFitWrapper.h.

197{return m_wrapperTF1.get();}

◆ GetWrapperTF1RawPtr()

virtual TF1 * ZDCFitWrapper::GetWrapperTF1RawPtr ( ) const

inlinevirtualinherited

Definition at line 198 of file ZDCFitWrapper.h.

198{return m_wrapperTF1.get();}

◆ HaveParameterSet()

bool ZDCFitWrapper::HaveParameterSet ( ) const

inlineinherited

Definition at line 127 of file ZDCFitWrapper.h.

127{return m_haveParameterSet;}

◆ Initialize() [1/2]

void ZDCFitWrapper::Initialize	(	float	initialAmp,
		float	initialT0,
		float	ampMin,
		float	ampMax )

inherited

Definition at line 199 of file ZDCFitWrapper.cxx.

{
  // Clear the errors on the TF1 because retaining the errors seems to affect 
  //   the convergence of fits 
  //
  auto tf1_p = GetWrapperTF1(); 
  unsigned int npar = tf1_p->GetNpar();
  std::vector<double> zeroVec(npar, 0.0);
  GetWrapperTF1()->SetParErrors(&zeroVec[0]);
 
  // If we adjusted the time range on the previous event, restore to default
  //
  if (m_adjTLimitsEvent) {
    SetT0FitLimits(m_t0Min, m_t0Max);
 
    m_adjTLimitsEvent = false;
    m_tminAdjust = 0;
  }
 
  SetAmpMinMax(ampMin, ampMax);
 
  //
  // Do we have parameter sets?
  //
  if (m_haveParameterSet) {  
    for (unsigned int ipar = 0; ipar < m_parameterSet.size(); ipar++) {
      //
      // Yes, use parameter descriptors to initialize the TF1
      //
      const auto & fitparDescr = m_parameterSet[ipar];
      
      if (!fitparDescr.fixed()) {
    tf1_p->SetParameter(ipar, fitparDescr.defaultValue());
    tf1_p->SetParError(ipar, fitparDescr.stepSize());
      }
    }
 
    tf1_p->SetParameter(0, initialAmp);    
    tf1_p->SetParameter(1, initialT0);    
  }
  else {
    //
    // No, perform the wrapper-implemented initialize
    //
    DoInitialize(initialAmp, initialT0, ampMin, ampMax);
 
    // Clear the errors on the TF1 because retaining the errors seems to affect 
    //   the convergence of fits
    //
    unsigned int npar = GetWrapperTF1()->GetNpar();
    std::vector<double> zeroVec(npar, 0.0);
    tf1_p->SetParErrors(&zeroVec[0]);
 
  }
}

◆ Initialize() [2/2]

void ZDCFitWrapper::Initialize	(	float	initialAmp,
		float	initialT0,
		float	ampMin,
		float	ampMax,
		float	fitTmin,
		float	fitTmax,
		float	fitTRef )

inherited

Definition at line 255 of file ZDCFitWrapper.cxx.

{
 
  // First do the nominal initializatioon
  //
  Initialize(initialAmp, initialT0, ampMin, ampMax);
 
  // Now account for the modified time range
  //
  m_adjTLimitsEvent = true;
 
  // Handle the specific limits on T0 
  //
  m_tminAdjust = fitTRef; // Note: this is the time corresponding to the sample used for presample subtraction 
 
  m_tempTmin = std::max(m_t0Min, fitTmin);
  m_tempTmax = std::min(m_t0Max, fitTmax);
  
  SetAmpMinMax(ampMin, ampMax);
  SetT0FitLimits(m_tempTmin, m_tempTmax);
}

◆ LoadDefaultParameterSet()

void ZDCFitWrapper::LoadDefaultParameterSet ( const std::vector< FitParameter > & parameterSet )

protectedinherited

Definition at line 126 of file ZDCFitWrapper.cxx.

{
  // First step is to check consistency with number of parameters in TF1
  //
  if (int(parameterSet.size()) != m_wrapperTF1->GetNpar()) throw std::logic_error("ZDCFitWrapper::LoadDefaultParameterSet:: Incorrect number of parameters in default parameter set");
 
  m_parameterSet = parameterSet;
  m_haveParameterSet = true;
 
  std::shared_ptr<TF1> theTF1 = ZDCFitWrapper::GetWrapperTF1();
 
  // Loop over passed descriptors
  //
  unsigned parIndex = 0;
  for (const auto & parDescr : parameterSet) {
    //
    // Load the settings for this parameter in the TF1
    //
    std::string name = parDescr.name();
    theTF1->SetParName(parIndex, name.c_str());
    
    bool fixed = parDescr.fixed();
    
    if (fixed) {
      float fixedValue = parDescr.defaultValue();
      theTF1->FixParameter(parIndex, fixedValue);
    }
    else {
      float lowerLimit = parDescr.lowerLimit();
      float upperLimit = parDescr.upperLimit();
      theTF1->SetParLimits(parIndex, lowerLimit, upperLimit);
    }
 
    parIndex++;
  }
}

◆ operator()()

virtual double ZDCFitExpFermiInductPreExp::operator()	(	const double *	x,
		const double *	p )

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1068 of file ZDCFitWrapper.h.

  {
    double t = x[0];
 
    double amp = p[0];
    double t0 = p[1];
 
    double expAmp = p[2];
    double expTau = p[3];
    double bexpSqrt = p[4];
    double C = p[5];
    double period = p[6];
    double inductA = p[7];
    double inductB = p[8];
    double bI = p[9];
 
    m_expFermiInductFunc->SetParameter(0, amp);
    m_expFermiInductFunc->SetParameter(5, period);
    m_expFermiInductFunc->SetParameter(6, inductA);
    m_expFermiInductFunc->SetParameter(7, inductB);
    m_expFermiInductFunc->SetParameter(8, bI);
 
    double deltaT = t - t0;
    double pulse = m_expFermiInductFunc->operator()(deltaT);
 
    // We subtract off the value of the exponential pulse at the minimum time (nominally 0),
    //   because it would have been included in the baseline subtraction
    //
    double tRef = GetTMinAdjust();
    double expPre = 0;
 
    expPre = expAmp * (std::exp(-t/expTau - bexpSqrt*t*t) - std::exp(-tRef/expTau - bexpSqrt*tRef*tRef));
      
    return pulse + C + expPre;
   }

◆ OverrideParameterSet()

void ZDCFitWrapper::OverrideParameterSet ( const std::vector< FitParameter > & parameterSet )

inherited

Definition at line 163 of file ZDCFitWrapper.cxx.

{
  std::shared_ptr<TF1> theTF1 = ZDCFitWrapper::GetWrapperTF1();
 
  // Loop over the provided parameter descriptors
  //
  for (const auto & parDescr : overParamSet) {
    std::string inputName = parDescr.name();
 
    // Fit the parameter name in the current list
    //
    auto iter = std::find_if(m_parameterSet.begin(), m_parameterSet.end(),
                 [inputName](const FitParameter& parDesc){return parDesc.name() == inputName;});
    if (iter == m_parameterSet.end()) throw std::logic_error("ZDCFitWrapper::OverrideParameterSet: invalid parameter name: " +inputName );
 
    // Update the descriptor for this parameter
    //
    (*iter) = parDescr;
    
    // Update the settings in the TF!
    //
    unsigned int parIndex = std::distance(m_parameterSet.begin(), iter);
    bool fixed = parDescr.fixed();
 
    if (fixed) {
      float fixedValue = parDescr.defaultValue();
      theTF1->FixParameter(parIndex, fixedValue);
    }
    else {
      float lowerLimit = parDescr.lowerLimit();
      float upperLimit = parDescr.upperLimit();
      theTF1->SetParLimits(parIndex, lowerLimit, upperLimit);
    }
  }
}

◆ SetAmpMinMax()

void ZDCFitWrapper::SetAmpMinMax	(	float	minAmp,
		float	maxAmp )

inlineinherited

Definition at line 147 of file ZDCFitWrapper.h.

  {
    m_ampMin = minAmp;
    m_ampMax = maxAmp;
  }

◆ SetInitialExpPulse()

virtual void ZDCFitExpFermiInductPreExp::SetInitialExpPulse ( float amp )

inlineoverridevirtual

Implements ZDCPreExpFitWrapper.

Definition at line 1010 of file ZDCFitWrapper.h.

  {
    GetWrapperTF1()->SetParameter(2, std::max(amp, (float) 0.5)); //0.5 here ensures that we're above lower limit (0)   
  }

◆ SetT0FitLimits()

void ZDCFitExpFermiInductPreExp::SetT0FitLimits	(	float	tMin,
		float	tMax )

overridevirtual

Implements ZDCFitWrapper.

Definition at line 1202 of file ZDCFitWrapper.cxx.

{
  std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
  theTF1->SetParLimits(1, t0Min, t0Max);
}

◆ SetT0Range()

void ZDCFitWrapper::SetT0Range	(	float	t0Min,
		float	t0Max )

inlineinherited

Definition at line 153 of file ZDCFitWrapper.h.

  {
    m_t0Min = t0Min;
    m_t0Max = t0Max;
 
    SetT0FitLimits(t0Min, t0Max);
  }

◆ UnconstrainFit()

void ZDCFitWrapper::UnconstrainFit ( )

inherited

Definition at line 293 of file ZDCFitWrapper.cxx.

{
  if (m_haveParameterSet) {  
    auto tf1_p = GetWrapperTF1();
    
    for (unsigned int ipar = 2; ipar < m_parameterSet.size(); ipar++) {
      const auto & fitparDescr = m_parameterSet[ipar];
      
      if (!fitparDescr.fixed() && fitparDescr.constrainFixed()) {
    tf1_p->ReleaseParameter(ipar);
    tf1_p->SetParameter(ipar, fitparDescr.defaultValue());
 
    float lowerLimit = fitparDescr.lowerLimit();
    float upperLimit = fitparDescr.upperLimit();
    tf1_p->SetParLimits(ipar, lowerLimit, upperLimit);
 
    tf1_p->SetParError(ipar, fitparDescr.stepSize());
      }
    }
  }
  else doUnconstrainFit();
}

Member Data Documentation

◆ m_adjTLimitsEvent

bool ZDCFitWrapper::m_adjTLimitsEvent {false}

privateinherited

Definition at line 74 of file ZDCFitWrapper.h.

74{false};

◆ m_ampMax

float ZDCFitWrapper::m_ampMax {0}

privateinherited

Definition at line 69 of file ZDCFitWrapper.h.

69{0};

◆ m_ampMin

float ZDCFitWrapper::m_ampMin {0}

privateinherited

Definition at line 68 of file ZDCFitWrapper.h.

68{0};

◆ m_defaultTau

float ZDCPreExpFitWrapper::m_defaultTau

protectedinherited

Definition at line 238 of file ZDCFitWrapper.h.

◆ m_expFermiInductFunc

std::shared_ptr<TF1> ZDCFitExpFermiInductPreExp::m_expFermiInductFunc {0}

private

Definition at line 982 of file ZDCFitWrapper.h.

982{0};

◆ m_expFermiPreFunc

std::shared_ptr<TF1> ZDCFitExpFermiInductPreExp::m_expFermiPreFunc {0}

private

Definition at line 983 of file ZDCFitWrapper.h.

983{0};

◆ m_finalized

bool ZDCFitWrapper::m_finalized {false}

privateinherited

Definition at line 82 of file ZDCFitWrapper.h.

82{false};

◆ m_fixTau

bool ZDCPreExpFitWrapper::m_fixTau

protectedinherited

Definition at line 239 of file ZDCFitWrapper.h.

◆ m_funcMax

float ZDCFitExpFermiInductPreExp::m_funcMax {-999}

private

Definition at line 980 of file ZDCFitWrapper.h.

980{-999};

◆ m_haveParameterSet

bool ZDCFitWrapper::m_haveParameterSet {false}

privateinherited

Definition at line 79 of file ZDCFitWrapper.h.

79{false};

◆ m_parameterSet

std::vector<FitParameter> ZDCFitWrapper::m_parameterSet

privateinherited

Definition at line 80 of file ZDCFitWrapper.h.

◆ m_t0Max

float ZDCFitWrapper::m_t0Max {0}

privateinherited

Definition at line 72 of file ZDCFitWrapper.h.

72{0};

◆ m_t0Min

float ZDCFitWrapper::m_t0Min {0}

privateinherited

Definition at line 71 of file ZDCFitWrapper.h.

71{0};

◆ m_tau1

float ZDCFitExpFermiInductPreExp::m_tau1 {0}

private

Definition at line 976 of file ZDCFitWrapper.h.

976{0};

◆ m_tau2

float ZDCFitExpFermiInductPreExp::m_tau2 {0}

private

Definition at line 977 of file ZDCFitWrapper.h.

977{0};

◆ m_tempTmax

float ZDCFitWrapper::m_tempTmax {0}

privateinherited

Definition at line 77 of file ZDCFitWrapper.h.

77{0};

◆ m_tempTmin

float ZDCFitWrapper::m_tempTmin {0}

privateinherited

Definition at line 76 of file ZDCFitWrapper.h.

76{0};

◆ m_tFuncMax

float ZDCFitExpFermiInductPreExp::m_tFuncMax {-999}

private

Definition at line 979 of file ZDCFitWrapper.h.

979{-999};

◆ m_tmax

float ZDCFitWrapper::m_tmax {0}

privateinherited

Definition at line 66 of file ZDCFitWrapper.h.

66{0};

◆ m_tmin

float ZDCFitWrapper::m_tmin {0}

privateinherited

Definition at line 65 of file ZDCFitWrapper.h.

65{0};

◆ m_tminAdjust

float ZDCFitWrapper::m_tminAdjust {0}

privateinherited

Definition at line 75 of file ZDCFitWrapper.h.

75{0};

◆ m_wrapperTF1

std::shared_ptr<TF1> ZDCFitWrapper::m_wrapperTF1 {}

privateinherited

Definition at line 63 of file ZDCFitWrapper.h.

63{};

◆ s_ConfigParamMap

const std::map< std::string, std::vector< ZDCFitWrapper::FitParameter > > ZDCFitExpFermiInductPreExp::s_ConfigParamMap

staticprivate

Initial value:

                                                          = 
{
  {"default",
   {{{"A", 1000, false, 2, 5000, 50, false, 0}},
    {{"T0", 40, false,  0, 100, 1, false, 0}},
    {{"Aexp", 0, false,  0.5, 3, 0.05, false, 0}},  
    {{"tauexp", 0, false,  3.5, 8, 0.25, false, 0}}, 
    {{"bsqrtexp", 0, false,  -50, 50, 1, true, 0}}, 
    {{"C", 0, false,  -50, 50, 1, true, 0}},
    {{"period",  13.5,  true,  0, 0, 0, false, 0}}, 
    {{"IndA", 0.2, false,  -2, 2, 0.1, false, 0}}, 
    {{"IndB", 0.2, false,  -2, 2, 0.1, false, 0}},
    {{"bI",  -0.02, true,  -0.5, 0.5, 0.025, true, -0.02}}
   }
  },
  {"PbPb2025",
   {{{"period",  13.5,  true,  0, 0, 0, false, 0}}, 
    {{"IndA",     0.1,  false,  -0.5, 2, 0.1, true, 0}}, 
    {{"IndB",     0.2,   false,  -0.5, 2, 0.1, false, 0}}, 
    {{"bI",     -0.02,   true,  -0.5, 0.5, 0.025, true, -0.02}}, 
   }
  },
  {"PbPb2025FreeInduct",
   {{{"period",  13,     false,  10, 30, 1, true, 13}}, 
    {{"IndA",     0.25,  false,  -1, 2, 0.1, true, 0}}, 
    {{"IndB",     0.8,   false,  -1, 2, 0.1, false, 0}}, 
    {{"bI",     -0.025,  false,  -0.5, 0.5, 0, false, 0}},
   }
  }
}

Definition at line 57 of file ZDCFitWrapper.h.

                                {return std::get<6>(m_parData);}
    float constrainValue() const {return std::get<7>(m_parData);}
  };
 
 
private:
  std::shared_ptr<TF1> m_wrapperTF1{};
  
  float m_tmin{0};
  float m_tmax{0};
 
  float m_ampMin{0};
  float m_ampMax{0};
 
  float m_t0Min{0};
  float m_t0Max{0};
 
  bool  m_adjTLimitsEvent{false};
  float m_tminAdjust{0};
  float m_tempTmin{0};
  float m_tempTmax{0};
 
  bool m_haveParameterSet{false};
  std::vector<FitParameter> m_parameterSet;
 
  bool m_finalized{false};
  
protected:
  // Actually sets the t0 parameter limits in the fit function
  //

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

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ZDCFitExpFermiInductPreExp() [1/2]

◆ ZDCFitExpFermiInductPreExp() [2/2]

◆ ~ZDCFitExpFermiInductPreExp()

Member Function Documentation

◆ ConstrainFit()

◆ doConstrainFit()

◆ DoFinalize()

◆ DoInitialize()

◆ doUnconstrainFit()

◆ Finalize()

◆ finalized()

◆ fixExpTau()

◆ GetAmpError()

◆ GetAmplitude()

◆ GetBkgdMaxFraction()

◆ getDefaultExpTau()

◆ GetExpAmp()

◆ GetExpTau()

◆ GetMaxAmp()

◆ GetMinAmp()

◆ GetNumShapeParameters()

◆ GetParameterDescr()

◆ GetShapeParameter()

◆ GetT0Max()

◆ GetT0Min()

◆ GetTau1()

◆ GetTau2()

◆ GetTime()

◆ GetTMax()

◆ GetTMin()

◆ GetTMinAdjust()

◆ GetWrapperTF1() [1/2]

◆ GetWrapperTF1() [2/2]

◆ GetWrapperTF1RawPtr()

◆ HaveParameterSet()

◆ Initialize() [1/2]

◆ Initialize() [2/2]

◆ LoadDefaultParameterSet()

◆ operator()()

◆ OverrideParameterSet()

◆ SetAmpMinMax()

◆ SetInitialExpPulse()

◆ SetT0FitLimits()

◆ SetT0Range()

◆ UnconstrainFit()

Member Data Documentation

◆ m_adjTLimitsEvent

◆ m_ampMax

◆ m_ampMin

◆ m_defaultTau

◆ m_expFermiInductFunc

◆ m_expFermiPreFunc

◆ m_finalized

◆ m_fixTau

◆ m_funcMax

◆ m_haveParameterSet

◆ m_parameterSet

◆ m_t0Max

◆ m_t0Min

◆ m_tau1

◆ m_tau2

◆ m_tempTmax

◆ m_tempTmin

◆ m_tFuncMax

◆ m_tmax

◆ m_tmin

◆ m_tminAdjust

◆ m_wrapperTF1

◆ s_ConfigParamMap