#include <ZDCFitWrapper.h>

Inheritance diagram for ZDCFitExpFermiVariableTausLHCf:

Collaboration diagram for ZDCFitExpFermiVariableTausLHCf:

Public Member Functions
	ZDCFitExpFermiVariableTausLHCf (const std::string &tag, float tmin, float tmax, bool fixTau1, bool fixTau2, float tau1, float tau2)
virtual void	DoInitialize (float initialAmp, float initialT0, float ampMin, float ampMax) override
virtual void	SetT0FitLimits (float tMin, float tMax) 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 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
virtual void	ConstrainFit () override
virtual void	UnconstrainFit () override
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)
void	SetAmpMinMax (float minAmp, float maxAmp)
void	SetT0Range (float t0Min, float t0Max)
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 Attributes
bool	m_fixTau1 {false}
bool	m_fixTau2 {false}
float	m_tau1 {0}
float	m_tau2 {0}

Private Attributes
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}

Detailed Description

Definition at line 236 of file ZDCFitWrapper.h.

Constructor & Destructor Documentation

◆ ZDCFitExpFermiVariableTausLHCf()

ZDCFitExpFermiVariableTausLHCf::ZDCFitExpFermiVariableTausLHCf	(	const std::string &	tag,
		float	tmin,
		float	tmax,
		bool	fixTau1,
		bool	fixTau2,
		float	tau1,
		float	tau2 )

Definition at line 125 of file ZDCFitWrapper.cxx.

                                                                                                                                                               :
  ZDCFitWrapper(std::make_shared<TF1>(("ExpFermiVariableTausLHCf" + tag).c_str(), ZDCFermiExpFitRefl, tmin, tmax, 9)),
  m_fixTau1(fixTau1), m_fixTau2(fixTau2), m_tau1(tau1), m_tau2(tau2)
{
  std::shared_ptr<TF1> theTF1 = ZDCFitWrapper::GetWrapperTF1();
 
  theTF1->SetParName(0, "Amp");
  theTF1->SetParName(1, "T0");
  theTF1->SetParName(2, "#tau_{1}");
  theTF1->SetParName(3, "#tau_{2}");
  theTF1->SetParName(4, "C");
  theTF1->SetParName(5, "refdelay");
  theTF1->SetParName(6, "reflAmpFrac");
  theTF1->SetParName(7, "reflWidth");
  theTF1->SetParName(8, "delta");
 
  theTF1->SetParLimits(1, tmin, tmax);
  theTF1->SetParLimits(4, -100, 100);
  theTF1->SetParLimits(6, -1e-4, 0.35);
 
  theTF1->FixParameter(5, 6.5);
  theTF1->FixParameter(7, 1.5);
  theTF1->FixParameter(8, 0.01);
 
  if (m_fixTau1) theTF1->FixParameter(2, m_tau1);
  else theTF1->SetParLimits(2, 1, 2);
    
  if (m_fixTau2) theTF1->FixParameter(3, m_tau2);
  else theTF1->SetParLimits(3, 3.5, 6);
}

Member Function Documentation

◆ ConstrainFit()

void ZDCFitExpFermiVariableTausLHCf::ConstrainFit ( )

overridevirtual

Implements ZDCFitWrapper.

Definition at line 199 of file ZDCFitWrapper.cxx.

{
  // We force the constant and reflection terms to zero
  //
  std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
 
  theTF1->FixParameter(4, 0);
  theTF1->FixParameter(6, 0);
}

◆ DoInitialize()

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

overridevirtual

Implements ZDCFitWrapper.

Definition at line 156 of file ZDCFitWrapper.cxx.

{
  std::shared_ptr<TF1> theTF1 = ZDCFitWrapper::GetWrapperTF1();
 
  theTF1->SetParameter(0, initialAmp);
 
  float t0 = initialT0;
  if (t0 < GetT0Min()) t0 = GetT0Min()*1.1;
  if (t0 > GetT0Max()) t0 = GetT0Max()/1.1;
  theTF1->SetParameter(1, t0);
 
  theTF1->SetParameter(1, initialT0);
  theTF1->SetParameter(4, 0);
  theTF1->SetParameter(6, 0.1);
 
  theTF1->SetParLimits(0, ampMin, ampMax);
 
  if (!m_fixTau1) theTF1->SetParameter(2, m_tau1);
  if (!m_fixTau2) theTF1->SetParameter(3, m_tau2);
 
  // Set the parameter errors which ROOT now uses <<<to set the initial step sizes>>>
  //
  double ampStep = std::min(0.05*initialAmp, std::abs(ampMax - initialAmp)/2.);
  
  theTF1->SetParError(0, ampStep);
  theTF1->SetParError(1, 1.0);
 
  if (!m_fixTau1) theTF1->SetParError(2, 0.05);
  if (!m_fixTau2) theTF1->SetParError(3, 0.25);
 
  theTF1->SetParError(4, 2);
  theTF1->SetParError(5, 0.25);
  theTF1->SetParError(6, 0.05);
}

◆ GetAmpError()

virtual float ZDCFitExpFermiVariableTausLHCf::GetAmpError ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 253 of file ZDCFitWrapper.h.

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

◆ GetAmplitude()

virtual float ZDCFitExpFermiVariableTausLHCf::GetAmplitude ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 252 of file ZDCFitWrapper.h.

252{return GetWrapperTF1()->GetParameter(0); }

◆ GetBkgdMaxFraction()

virtual float ZDCFitExpFermiVariableTausLHCf::GetBkgdMaxFraction ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 280 of file ZDCFitWrapper.h.

  {
    const TF1* theTF1 = ZDCFitWrapper::GetWrapperTF1();
    double amp = theTF1->GetParameter(0);
    double constant = theTF1->GetParameter(4);
 
    if (amp > 0) return constant / amp;
    else return 1;
  }

◆ GetMaxAmp()

float ZDCFitWrapper::GetMaxAmp ( ) const

inlineinherited

Definition at line 91 of file ZDCFitWrapper.h.

91{return m_ampMax;}

ZDCFitWrapper::m_ampMax

float m_ampMax

Definition ZDCFitWrapper.h:31

◆ GetMinAmp()

float ZDCFitWrapper::GetMinAmp ( ) const

inlineinherited

Definition at line 90 of file ZDCFitWrapper.h.

90{return m_ampMin;}

ZDCFitWrapper::m_ampMin

float m_ampMin

Definition ZDCFitWrapper.h:30

◆ GetNumShapeParameters()

virtual unsigned int ZDCFitExpFermiVariableTausLHCf::GetNumShapeParameters ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 272 of file ZDCFitWrapper.h.

272{return 5;}

◆ GetShapeParameter()

virtual float ZDCFitExpFermiVariableTausLHCf::GetShapeParameter ( size_t index ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 274 of file ZDCFitWrapper.h.

  {
    if (index < 5) return GetWrapperTF1()->GetParameter(5+index);
    else throw std::runtime_error("Fit parameter does not exist.");
  }

◆ GetT0Max()

float ZDCFitWrapper::GetT0Max ( ) const

inlineinherited

Definition at line 101 of file ZDCFitWrapper.h.

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

◆ GetT0Min()

float ZDCFitWrapper::GetT0Min ( ) const

inlineinherited

Definition at line 96 of file ZDCFitWrapper.h.

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

◆ GetTau1()

virtual float ZDCFitExpFermiVariableTausLHCf::GetTau1 ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 255 of file ZDCFitWrapper.h.

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

◆ GetTau2()

virtual float ZDCFitExpFermiVariableTausLHCf::GetTau2 ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 256 of file ZDCFitWrapper.h.

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

◆ GetTime()

virtual float ZDCFitExpFermiVariableTausLHCf::GetTime ( ) const

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 258 of file ZDCFitWrapper.h.

                                         {
    const TF1* theTF1 = GetWrapperTF1();
 
    float fitT0 =  theTF1->GetParameter(1);
 
    float tau1 = theTF1->GetParameter(2);
    float tau2 = theTF1->GetParameter(3);
 
    // Correct the time to the maximum
    //
    if (tau2 > tau1) fitT0 += tau1 * std::log(tau2 / tau1 - 1.0);
    return fitT0;
  }

◆ GetTMax()

float ZDCFitWrapper::GetTMax ( ) const

inlineinherited

Definition at line 94 of file ZDCFitWrapper.h.

94{return m_tmax;}

ZDCFitWrapper::m_tmax

float m_tmax

Definition ZDCFitWrapper.h:28

◆ GetTMin()

float ZDCFitWrapper::GetTMin ( ) const

inlineinherited

Definition at line 93 of file ZDCFitWrapper.h.

93{return m_tmin;}

ZDCFitWrapper::m_tmin

float m_tmin

Definition ZDCFitWrapper.h:27

◆ GetTMinAdjust()

float ZDCFitWrapper::GetTMinAdjust ( ) const

inlineinherited

Definition at line 106 of file ZDCFitWrapper.h.

106{return m_tminAdjust;}

ZDCFitWrapper::m_tminAdjust

float m_tminAdjust

Definition ZDCFitWrapper.h:37

◆ GetWrapperTF1() [1/2]

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

inlinevirtualinherited

Definition at line 115 of file ZDCFitWrapper.h.

115{return m_wrapperTF1;}

ZDCFitWrapper::m_wrapperTF1

std::shared_ptr< TF1 > m_wrapperTF1

Definition ZDCFitWrapper.h:25

◆ GetWrapperTF1() [2/2]

virtual const TF1 * ZDCFitWrapper::GetWrapperTF1 ( ) const

inlinevirtualinherited

Definition at line 116 of file ZDCFitWrapper.h.

116{return m_wrapperTF1.get();}

◆ GetWrapperTF1RawPtr()

virtual TF1 * ZDCFitWrapper::GetWrapperTF1RawPtr ( ) const

inlinevirtualinherited

Definition at line 117 of file ZDCFitWrapper.h.

117{return m_wrapperTF1.get();}

◆ Initialize() [1/2]

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

inherited

Definition at line 9 of file ZDCFitWrapper.cxx.

{
  // 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);
  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);
 
  DoInitialize(initialAmp, initialT0, ampMin, ampMax);
}

◆ Initialize() [2/2]

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

inherited

Definition at line 32 of file ZDCFitWrapper.cxx.

{
  // 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);
  GetWrapperTF1()->SetParErrors(&zeroVec[0]);
 
  m_adjTLimitsEvent = true;
 
  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);
 
  DoInitialize(initialAmp, initialT0, ampMin, ampMax);
}

◆ operator()()

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

inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 290 of file ZDCFitWrapper.h.

                                                                       {
    return ZDCFermiExpFitRefl(x, p);
  }

◆ SetAmpMinMax()

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

inlineinherited

Definition at line 67 of file ZDCFitWrapper.h.

  {
    m_ampMin = minAmp;
    m_ampMax = maxAmp;
  }

◆ SetT0FitLimits()

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

overridevirtual

Implements ZDCFitWrapper.

Definition at line 191 of file ZDCFitWrapper.cxx.

{
  // Set the parameter limits accordingly on the TF1
  //
  std::shared_ptr<TF1> theTF1 = ZDCFitWrapper::GetWrapperTF1();
  theTF1->SetParLimits(1, t0Min, t0Max);
}

◆ SetT0Range()

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

inlineinherited

Definition at line 73 of file ZDCFitWrapper.h.

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

◆ UnconstrainFit()

void ZDCFitExpFermiVariableTausLHCf::UnconstrainFit ( )

overridevirtual

Implements ZDCFitWrapper.

Definition at line 208 of file ZDCFitWrapper.cxx.

{
  std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
  theTF1->ReleaseParameter(4);
  theTF1->SetParameter(4, 0);
 
  theTF1->ReleaseParameter(6);
  theTF1->SetParLimits(6, -1e-4, 0.35);
  theTF1->SetParameter(6, 0.1);
}

Member Data Documentation

◆ m_adjTLimitsEvent

bool ZDCFitWrapper::m_adjTLimitsEvent {false}

privateinherited

Definition at line 36 of file ZDCFitWrapper.h.

36{false};

◆ m_ampMax

float ZDCFitWrapper::m_ampMax {0}

privateinherited

Definition at line 31 of file ZDCFitWrapper.h.

31{0};

◆ m_ampMin

float ZDCFitWrapper::m_ampMin {0}

privateinherited

Definition at line 30 of file ZDCFitWrapper.h.

30{0};

◆ m_fixTau1

bool ZDCFitExpFermiVariableTausLHCf::m_fixTau1 {false}

protected

Definition at line 239 of file ZDCFitWrapper.h.

239{false};

◆ m_fixTau2

bool ZDCFitExpFermiVariableTausLHCf::m_fixTau2 {false}

protected

Definition at line 240 of file ZDCFitWrapper.h.

240{false};

◆ m_t0Max

float ZDCFitWrapper::m_t0Max {0}

privateinherited

Definition at line 34 of file ZDCFitWrapper.h.

34{0};

◆ m_t0Min

float ZDCFitWrapper::m_t0Min {0}

privateinherited

Definition at line 33 of file ZDCFitWrapper.h.

33{0};

◆ m_tau1

float ZDCFitExpFermiVariableTausLHCf::m_tau1 {0}

protected

Definition at line 242 of file ZDCFitWrapper.h.

242{0};

◆ m_tau2

float ZDCFitExpFermiVariableTausLHCf::m_tau2 {0}

protected

Definition at line 243 of file ZDCFitWrapper.h.

243{0};

◆ m_tempTmax

float ZDCFitWrapper::m_tempTmax {0}

privateinherited

Definition at line 39 of file ZDCFitWrapper.h.

39{0};

◆ m_tempTmin

float ZDCFitWrapper::m_tempTmin {0}

privateinherited

Definition at line 38 of file ZDCFitWrapper.h.

38{0};

◆ m_tmax

float ZDCFitWrapper::m_tmax {0}

privateinherited

Definition at line 28 of file ZDCFitWrapper.h.

28{0};

◆ m_tmin

float ZDCFitWrapper::m_tmin {0}

privateinherited

Definition at line 27 of file ZDCFitWrapper.h.

27{0};

◆ m_tminAdjust

float ZDCFitWrapper::m_tminAdjust {0}

privateinherited

Definition at line 37 of file ZDCFitWrapper.h.

37{0};

◆ m_wrapperTF1

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

privateinherited

Definition at line 25 of file ZDCFitWrapper.h.

25{};

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

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ZDCFitExpFermiVariableTausLHCf()

Member Function Documentation

◆ ConstrainFit()

◆ DoInitialize()

◆ GetAmpError()

◆ GetAmplitude()

◆ GetBkgdMaxFraction()

◆ GetMaxAmp()

◆ GetMinAmp()

◆ GetNumShapeParameters()

◆ GetShapeParameter()

◆ GetT0Max()

◆ GetT0Min()

◆ GetTau1()

◆ GetTau2()

◆ GetTime()

◆ GetTMax()

◆ GetTMin()

◆ GetTMinAdjust()

◆ GetWrapperTF1() [1/2]

◆ GetWrapperTF1() [2/2]

◆ GetWrapperTF1RawPtr()

◆ Initialize() [1/2]

◆ Initialize() [2/2]

◆ operator()()

◆ SetAmpMinMax()

◆ SetT0FitLimits()

◆ SetT0Range()

◆ UnconstrainFit()

Member Data Documentation

◆ m_adjTLimitsEvent

◆ m_ampMax

◆ m_ampMin

◆ m_fixTau1

◆ m_fixTau2

◆ m_t0Max

◆ m_t0Min

◆ m_tau1

◆ m_tau2

◆ m_tempTmax

◆ m_tempTmin

◆ m_tmax

◆ m_tmin

◆ m_tminAdjust

◆ m_wrapperTF1