ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | List of all members
ZDCFitExpFermiLinearFixedTaus Class Reference

#include <ZDCFitWrapper.h>

Inheritance diagram for ZDCFitExpFermiLinearFixedTaus:
Collaboration diagram for ZDCFitExpFermiLinearFixedTaus:

Public Member Functions

 ZDCFitExpFermiLinearFixedTaus (const std::string &tag, float tmin, float tmax, float tau1, float tau2)
 ~ZDCFitExpFermiLinearFixedTaus ()
virtual void DoInitialize (float initialAmp, float initialT0, float ampMin, float ampMax) override
virtual void SetT0FitLimits (float tMin, float tMax) override
virtual void ConstrainFit () override
virtual void UnconstrainFit () 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
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

Private Attributes

float m_tau1 {}
float m_tau2 {}
float m_norm {}
float m_timeCorr {}
std::shared_ptr< TF1 > m_expFermiFunc
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 958 of file ZDCFitWrapper.h.

Constructor & Destructor Documentation

◆ ZDCFitExpFermiLinearFixedTaus()

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

Definition at line 852 of file ZDCFitWrapper.cxx.

852 :
853 ZDCFitWrapper(std::make_shared<TF1>(("ExpFermiFixedTaus" + tag).c_str(), this, tmin, tmax, 4)),
854 m_tau1(tau1), m_tau2(tau2)
855{
856 std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
857
858 // BAC, parameter 0 limits now is set in DoInitialize
859 theTF1->SetParLimits(1, tmin, tmax);
860
861 theTF1->SetParName(0, "Amp");
862 theTF1->SetParName(1, "T0");
863 theTF1->SetParName(2, "s_{b}");
864 theTF1->SetParName(3, "c_{b}");
865
866 // Now create the reference function that we use to evaluate ExpFermiFit more efficiently
867 //
868 std::string funcNameRefFunc = "ExpFermiFixedTausRefFunc" + tag;
869
870 m_expFermiFunc = std::make_shared<TF1>(funcNameRefFunc.c_str(), ZDCFermiExpFit, -50, 100, 5);
871
872 m_expFermiFunc->SetParameter(0, 1);
873 m_expFermiFunc->SetParameter(1, 0);
874 m_expFermiFunc->SetParameter(2, m_tau1);
875 m_expFermiFunc->SetParameter(3, m_tau2);
876 m_expFermiFunc->FixParameter(4, 0);
877
878 m_norm = 1. / m_expFermiFunc->GetMaximum();
879 m_timeCorr = m_tau1 * std::log(m_tau2 / m_tau1 - 1.0);
880}
ZDCFermiExpFit
double ZDCFermiExpFit(const double *xvec, const double *pvec)
Definition ZDCFitWrapper.h:1390
ZDCFitExpFermiLinearFixedTaus::m_expFermiFunc
std::shared_ptr< TF1 > m_expFermiFunc
Definition ZDCFitWrapper.h:967
ZDCFitWrapper::GetWrapperTF1
virtual std::shared_ptr< TF1 > GetWrapperTF1()
Definition ZDCFitWrapper.h:115
ZDCFitWrapper::ZDCFitWrapper
ZDCFitWrapper(const std::shared_ptr< TF1 > &wrapperTF1)
Definition ZDCFitWrapper.h:47
CaloCondBlobAlgs_fillNoiseFromASCII.tag
str tag
Definition CaloCondBlobAlgs_fillNoiseFromASCII.py:23

◆ ~ZDCFitExpFermiLinearFixedTaus()

ZDCFitExpFermiLinearFixedTaus::~ZDCFitExpFermiLinearFixedTaus ( )
inline

Definition at line 973 of file ZDCFitWrapper.h.

973{}

Member Function Documentation

◆ ConstrainFit()

void ZDCFitExpFermiLinearFixedTaus::ConstrainFit ( )
overridevirtual

Implements ZDCFitWrapper.

Definition at line 882 of file ZDCFitWrapper.cxx.

883{
884 // We force the linear terms to zero
885 //
886 std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
887
888 theTF1->FixParameter(2, 0);
889 theTF1->FixParameter(3, 0);
890}

◆ DoInitialize()

void ZDCFitExpFermiLinearFixedTaus::DoInitialize ( float initialAmp,
float initialT0,
float ampMin,
float ampMax )
overridevirtual

Implements ZDCFitWrapper.

Definition at line 899 of file ZDCFitWrapper.cxx.

900{
901 float slope = std::abs(0.1 * initialAmp / initialT0);
902 float intercept = std::abs(0.1 * initialAmp);
903 GetWrapperTF1()->SetParLimits(2, -slope , slope );
904 GetWrapperTF1()->SetParLimits(3, -intercept, intercept);
905
906 GetWrapperTF1()->SetParameter(0, initialAmp);
907 GetWrapperTF1()->SetParameter(1, initialT0);
908 GetWrapperTF1()->SetParameter(2, 0);
909 GetWrapperTF1()->SetParameter(3, 0);
910
911 GetWrapperTF1()->SetParLimits(0, ampMin, ampMax);
912
913 double ampStep = std::min(0.05*initialAmp, std::abs(ampMax - initialAmp)/2.);
914
915 GetWrapperTF1()->SetParError(0, ampStep);
916 GetWrapperTF1()->SetParError(1, 1);
917 GetWrapperTF1()->SetParError(2, 1);
918 GetWrapperTF1()->SetParError(3, 0.5);
919}

◆ GetAmpError()

virtual float ZDCFitExpFermiLinearFixedTaus::GetAmpError ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 982 of file ZDCFitWrapper.h.

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

◆ GetAmplitude()

virtual float ZDCFitExpFermiLinearFixedTaus::GetAmplitude ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 981 of file ZDCFitWrapper.h.

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

◆ GetBkgdMaxFraction()

virtual float ZDCFitExpFermiLinearFixedTaus::GetBkgdMaxFraction ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 998 of file ZDCFitWrapper.h.

999 {
1000 const TF1* theTF1 = ZDCFitWrapper::GetWrapperTF1();
1001 double amp = theTF1->GetParameter(0);
1002 if (amp <= 0) return -1;
1003
1004 double slope = theTF1->GetParameter(2);
1005
1006 double background = slope * GetTime();
1007 return background / amp;
1008 }
ZDCFitExpFermiLinearFixedTaus::GetTime
virtual float GetTime() const override
Definition ZDCFitWrapper.h:987

◆ GetMaxAmp()

float ZDCFitWrapper::GetMaxAmp ( ) const
inlineinherited

Definition at line 91 of file ZDCFitWrapper.h.

91{return m_ampMax;}

◆ GetMinAmp()

float ZDCFitWrapper::GetMinAmp ( ) const
inlineinherited

Definition at line 90 of file ZDCFitWrapper.h.

90{return m_ampMin;}

◆ GetNumShapeParameters()

virtual unsigned int ZDCFitExpFermiLinearFixedTaus::GetNumShapeParameters ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 991 of file ZDCFitWrapper.h.

991{return 2;}

◆ GetShapeParameter()

virtual float ZDCFitExpFermiLinearFixedTaus::GetShapeParameter ( size_t index) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 992 of file ZDCFitWrapper.h.

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

◆ GetT0Max()

float ZDCFitWrapper::GetT0Max ( ) const
inlineinherited

Definition at line 101 of file ZDCFitWrapper.h.

101 {
102 if (m_adjTLimitsEvent) return m_tempTmax;
103 else return m_t0Max;
104 }

◆ GetT0Min()

float ZDCFitWrapper::GetT0Min ( ) const
inlineinherited

Definition at line 96 of file ZDCFitWrapper.h.

96 {
97 if (m_adjTLimitsEvent) return m_tempTmin;
98 else return m_t0Min;
99 }

◆ GetTau1()

virtual float ZDCFitExpFermiLinearFixedTaus::GetTau1 ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 984 of file ZDCFitWrapper.h.

984{return m_tau1;}

◆ GetTau2()

virtual float ZDCFitExpFermiLinearFixedTaus::GetTau2 ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 985 of file ZDCFitWrapper.h.

985{return m_tau2;}

◆ GetTime()

virtual float ZDCFitExpFermiLinearFixedTaus::GetTime ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 987 of file ZDCFitWrapper.h.

987 {
988 return GetWrapperTF1()->GetParameter(1) + m_timeCorr; // Correct the time to the maximum
989 }

◆ GetTMax()

float ZDCFitWrapper::GetTMax ( ) const
inlineinherited

Definition at line 94 of file ZDCFitWrapper.h.

94{return m_tmax;}

◆ GetTMin()

float ZDCFitWrapper::GetTMin ( ) const
inlineinherited

Definition at line 93 of file ZDCFitWrapper.h.

93{return m_tmin;}

◆ GetTMinAdjust()

float ZDCFitWrapper::GetTMinAdjust ( ) const
inlineinherited

Definition at line 106 of file ZDCFitWrapper.h.

106{return m_tminAdjust;}

◆ 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.

10{
11 // Clear the errors on the TF1 because retaining the errors seems to affect
12 // the convergence of fits
13 //
14 unsigned int npar = GetWrapperTF1()->GetNpar();
15 std::vector<double> zeroVec(npar, 0.0);
16 GetWrapperTF1()->SetParErrors(&zeroVec[0]);
17
18 // If we adjusted the time range on the previous event, restore to default
19 //
20 if (m_adjTLimitsEvent) {
21 SetT0FitLimits(m_t0Min, m_t0Max);
22
23 m_adjTLimitsEvent = false;
24 m_tminAdjust = 0;
25 }
26
27 SetAmpMinMax(ampMin, ampMax);
28
29 DoInitialize(initialAmp, initialT0, ampMin, ampMax);
30}
ZDCFitWrapper::SetT0FitLimits
virtual void SetT0FitLimits(float tMin, float tMax)=0
ZDCFitWrapper::SetAmpMinMax
void SetAmpMinMax(float minAmp, float maxAmp)
Definition ZDCFitWrapper.h:67
ZDCFitWrapper::DoInitialize
virtual void DoInitialize(float initialAmp, float initialT0, float ampMin, float ampMax)=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 32 of file ZDCFitWrapper.cxx.

33{
34 // Clear the errors on the TF1 because retaining the errors seems to affect
35 // the convergence of fits
36 //
37 unsigned int npar = GetWrapperTF1()->GetNpar();
38 std::vector<double> zeroVec(npar, 0.0);
39 GetWrapperTF1()->SetParErrors(&zeroVec[0]);
40
41 m_adjTLimitsEvent = true;
42
43 m_tminAdjust = fitTRef; // Note: this is the time corresponding to the sample used for presample subtraction
44
45 m_tempTmin = std::max(m_t0Min, fitTmin);
46 m_tempTmax = std::min(m_t0Max, fitTmax);
47
48 SetAmpMinMax(ampMin, ampMax);
49 SetT0FitLimits(m_tempTmin, m_tempTmax);
50
51 DoInitialize(initialAmp, initialT0, ampMin, ampMax);
52}

◆ operator()()

virtual double ZDCFitExpFermiLinearFixedTaus::operator() ( const double * x,
const double * p )
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 1010 of file ZDCFitWrapper.h.

1011 {
1012 double amp = p[0];
1013 double t0 = p[1];
1014 double deltaT = x[0] - t0;
1015
1016 double bckgd = p[2] * x[0] + p[3];
1017
1018 double expFermi = amp * m_norm * m_expFermiFunc->operator()(deltaT);
1019
1020 return expFermi + bckgd;
1021 }
t0
static Double_t t0
Definition LArPhysWaveHECTool.cxx:38
x
#define x

◆ SetAmpMinMax()

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

Definition at line 67 of file ZDCFitWrapper.h.

68 {
69 m_ampMin = minAmp;
70 m_ampMax = maxAmp;
71 }

◆ SetT0FitLimits()

void ZDCFitExpFermiLinearFixedTaus::SetT0FitLimits ( float tMin,
float tMax )
overridevirtual

Implements ZDCFitWrapper.

Definition at line 921 of file ZDCFitWrapper.cxx.

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

◆ SetT0Range()

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

Definition at line 73 of file ZDCFitWrapper.h.

74 {
75 m_t0Min = t0Min;
76 m_t0Max = t0Max;
77
78 SetT0FitLimits(t0Min, t0Max);
79 }

◆ UnconstrainFit()

void ZDCFitExpFermiLinearFixedTaus::UnconstrainFit ( )
overridevirtual

Implements ZDCFitWrapper.

Definition at line 892 of file ZDCFitWrapper.cxx.

893{
894 std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
895 theTF1->ReleaseParameter(2);
896 theTF1->ReleaseParameter(3);
897}

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_expFermiFunc

std::shared_ptr<TF1> ZDCFitExpFermiLinearFixedTaus::m_expFermiFunc
private

Definition at line 967 of file ZDCFitWrapper.h.

◆ m_norm

float ZDCFitExpFermiLinearFixedTaus::m_norm {}
private

Definition at line 964 of file ZDCFitWrapper.h.

964{};

◆ 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 ZDCFitExpFermiLinearFixedTaus::m_tau1 {}
private

Definition at line 961 of file ZDCFitWrapper.h.

961{};

◆ m_tau2

float ZDCFitExpFermiLinearFixedTaus::m_tau2 {}
private

Definition at line 962 of file ZDCFitWrapper.h.

962{};

◆ 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_timeCorr

float ZDCFitExpFermiLinearFixedTaus::m_timeCorr {}
private

Definition at line 965 of file ZDCFitWrapper.h.

965{};

◆ 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:
Generated on for ATLAS Offline Software by doxygen 1.16.1