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

#include <ZDCFitWrapper.h>

Inheritance diagram for ZDCFitExpFermiFixedTaus:
Collaboration diagram for ZDCFitExpFermiFixedTaus:

Public Member Functions

 ZDCFitExpFermiFixedTaus (const std::string &tag, float tmin, float tmax, float tau1, float tau2)
 ~ZDCFitExpFermiFixedTaus ()
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)
bool finalized () const
virtual void Finalize ()
virtual void DoFinalize ()
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 {0}
float m_tau2 {0}
float m_norm {0}
float m_timeCorr {0}
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}
bool m_finalized {false}

Detailed Description

Definition at line 396 of file ZDCFitWrapper.h.

Constructor & Destructor Documentation

◆ ZDCFitExpFermiFixedTaus()

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

Definition at line 317 of file ZDCFitWrapper.cxx.

317 :
318 ZDCFitWrapper(std::make_shared<TF1>(("ExpFermiFixedTaus" + tag).c_str(), this, tmin, tmax, 3)),
319 m_tau1(tau1), m_tau2(tau2)
320{
321 std::shared_ptr<TF1> theTF1 = ZDCFitWrapper::GetWrapperTF1();
322
323 // BAC, parameter 0 limits now is set in DoInitialize
324 theTF1->SetParLimits(1, tmin, tmax);
325
326 theTF1->SetParName(0, "Amp");
327 theTF1->SetParName(1, "T0");
328 theTF1->SetParName(2, "C");
329
330 // Now create the reference function that we use to evaluate ExpFermiFit more efficiently
331 //
332 std::string funcNameRefFunc = "ExpFermiFixedTausRefFunc" + tag;
333 m_expFermiFunc = std::make_shared<TF1>(funcNameRefFunc.c_str(), ZDCFermiExpFit, -50, 100, 5);
334
335 m_expFermiFunc->SetParameter(0, 1);
336 m_expFermiFunc->SetParameter(1, 0);
337 m_expFermiFunc->SetParameter(2, m_tau1);
338 m_expFermiFunc->SetParameter(3, m_tau2);
339 m_expFermiFunc->SetParameter(4, 0);
340
341 m_norm = 1. / m_expFermiFunc->GetMaximum();
342 m_timeCorr = m_tau1 * std::log(m_tau2 / m_tau1 - 1.0);
343}
ZDCFermiExpFit
double ZDCFermiExpFit(const double *xvec, const double *pvec)
Definition ZDCFitWrapper.h:1594
ZDCFitExpFermiFixedTaus::m_expFermiFunc
std::shared_ptr< TF1 > m_expFermiFunc
Definition ZDCFitWrapper.h:405
ZDCFitWrapper::GetWrapperTF1
virtual std::shared_ptr< TF1 > GetWrapperTF1()
Definition ZDCFitWrapper.h:128
ZDCFitWrapper::ZDCFitWrapper
ZDCFitWrapper(const std::shared_ptr< TF1 > &wrapperTF1)
Definition ZDCFitWrapper.h:50
CaloCondBlobAlgs_fillNoiseFromASCII.tag
str tag
Definition CaloCondBlobAlgs_fillNoiseFromASCII.py:23

◆ ~ZDCFitExpFermiFixedTaus()

ZDCFitExpFermiFixedTaus::~ZDCFitExpFermiFixedTaus ( )
inline

Definition at line 411 of file ZDCFitWrapper.h.

411{}

Member Function Documentation

◆ ConstrainFit()

void ZDCFitExpFermiFixedTaus::ConstrainFit ( )
overridevirtual

Implements ZDCFitWrapper.

Definition at line 352 of file ZDCFitWrapper.cxx.

353{
354 // We force the constant tersm to zero
355 //
356 std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
357
358 theTF1->FixParameter(2, 0);
359}

◆ DoFinalize()

virtual void ZDCFitWrapper::DoFinalize ( )
inlinevirtualinherited

Reimplemented in ZDCFitExpFermiInductPreExp, ZDCFitExpFermiInductPrePulse, and ZDCFitExpFermiVariableTausInduct.

Definition at line 77 of file ZDCFitWrapper.h.

77{}

◆ DoInitialize()

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

Implements ZDCFitWrapper.

Definition at line 345 of file ZDCFitWrapper.cxx.

346{
347 GetWrapperTF1()->SetParameter(0, initialAmp);
348 GetWrapperTF1()->SetParameter(1, initialT0);
349
350 GetWrapperTF1()->SetParLimits(0, ampMin, ampMax);
351}

◆ Finalize()

virtual void ZDCFitWrapper::Finalize ( )
inlinevirtualinherited

Definition at line 68 of file ZDCFitWrapper.h.

69 {
70 DoFinalize();
71 m_finalized = true;
72 }
ZDCFitWrapper::DoFinalize
virtual void DoFinalize()
Definition ZDCFitWrapper.h:77

◆ finalized()

bool ZDCFitWrapper::finalized ( ) const
inlineinherited

Definition at line 66 of file ZDCFitWrapper.h.

66{return m_finalized;}

◆ GetAmpError()

virtual float ZDCFitExpFermiFixedTaus::GetAmpError ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 420 of file ZDCFitWrapper.h.

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

◆ GetAmplitude()

virtual float ZDCFitExpFermiFixedTaus::GetAmplitude ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 419 of file ZDCFitWrapper.h.

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

◆ GetBkgdMaxFraction()

virtual float ZDCFitExpFermiFixedTaus::GetBkgdMaxFraction ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 438 of file ZDCFitWrapper.h.

439 {
440 const TF1* theTF1 = ZDCFitWrapper::GetWrapperTF1();
441 double amp = theTF1->GetParameter(0);
442 if (amp <= 0) return -1;
443
444 double C = theTF1->GetParameter(2);
445 return C / amp;
446 }
C
struct color C

◆ GetMaxAmp()

float ZDCFitWrapper::GetMaxAmp ( ) const
inlineinherited

Definition at line 104 of file ZDCFitWrapper.h.

104{return m_ampMax;}

◆ GetMinAmp()

float ZDCFitWrapper::GetMinAmp ( ) const
inlineinherited

Definition at line 103 of file ZDCFitWrapper.h.

103{return m_ampMin;}

◆ GetNumShapeParameters()

virtual unsigned int ZDCFitExpFermiFixedTaus::GetNumShapeParameters ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 429 of file ZDCFitWrapper.h.

429{return 2;}

◆ GetShapeParameter()

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

Implements ZDCFitWrapper.

Definition at line 431 of file ZDCFitWrapper.h.

432 {
433 if (index == 0) return m_tau1;
434 else if (index == 1) return m_tau2;
435 else throw std::runtime_error("Fit parameter does not exist.");
436 }

◆ GetT0Max()

float ZDCFitWrapper::GetT0Max ( ) const
inlineinherited

Definition at line 114 of file ZDCFitWrapper.h.

114 {
115 if (m_adjTLimitsEvent) return m_tempTmax;
116 else return m_t0Max;
117 }

◆ GetT0Min()

float ZDCFitWrapper::GetT0Min ( ) const
inlineinherited

Definition at line 109 of file ZDCFitWrapper.h.

109 {
110 if (m_adjTLimitsEvent) return m_tempTmin;
111 else return m_t0Min;
112 }

◆ GetTau1()

virtual float ZDCFitExpFermiFixedTaus::GetTau1 ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 422 of file ZDCFitWrapper.h.

422{return m_tau1;}

◆ GetTau2()

virtual float ZDCFitExpFermiFixedTaus::GetTau2 ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 423 of file ZDCFitWrapper.h.

423{return m_tau2;}

◆ GetTime()

virtual float ZDCFitExpFermiFixedTaus::GetTime ( ) const
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 425 of file ZDCFitWrapper.h.

425 {
426 return GetWrapperTF1()->GetParameter(1) + m_timeCorr; // Correct the time to the maximum
427 }

◆ GetTMax()

float ZDCFitWrapper::GetTMax ( ) const
inlineinherited

Definition at line 107 of file ZDCFitWrapper.h.

107{return m_tmax;}

◆ GetTMin()

float ZDCFitWrapper::GetTMin ( ) const
inlineinherited

Definition at line 106 of file ZDCFitWrapper.h.

106{return m_tmin;}

◆ GetTMinAdjust()

float ZDCFitWrapper::GetTMinAdjust ( ) const
inlineinherited

Definition at line 119 of file ZDCFitWrapper.h.

119{return m_tminAdjust;}

◆ GetWrapperTF1() [1/2]

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

Definition at line 128 of file ZDCFitWrapper.h.

128{return m_wrapperTF1;}
ZDCFitWrapper::m_wrapperTF1
std::shared_ptr< TF1 > m_wrapperTF1
Definition ZDCFitWrapper.h:26

◆ GetWrapperTF1() [2/2]

virtual const TF1 * ZDCFitWrapper::GetWrapperTF1 ( ) const
inlinevirtualinherited

Definition at line 129 of file ZDCFitWrapper.h.

129{return m_wrapperTF1.get();}

◆ GetWrapperTF1RawPtr()

virtual TF1 * ZDCFitWrapper::GetWrapperTF1RawPtr ( ) const
inlinevirtualinherited

Definition at line 130 of file ZDCFitWrapper.h.

130{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:79
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 ZDCFitExpFermiFixedTaus::operator() ( const double * x,
const double * p )
inlineoverridevirtual

Implements ZDCFitWrapper.

Definition at line 448 of file ZDCFitWrapper.h.

449 {
450 double amp = p[0];
451 double t0 = p[1];
452 double C = p[2];
453
454 double deltaT = x[0] - t0;
455
456 double expFermi = amp * m_norm * m_expFermiFunc->operator()(deltaT);
457
458 return expFermi + C; // + bckgd;
459 }
t0
static Double_t t0
Definition LArPhysWaveHECTool.cxx:38
x
#define x

◆ SetAmpMinMax()

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

Definition at line 79 of file ZDCFitWrapper.h.

80 {
81 m_ampMin = minAmp;
82 m_ampMax = maxAmp;
83 }

◆ SetT0FitLimits()

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

Implements ZDCFitWrapper.

Definition at line 366 of file ZDCFitWrapper.cxx.

367{
368 GetWrapperTF1()->SetParLimits(1, t0Min, t0Max);
369}

◆ SetT0Range()

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

Definition at line 85 of file ZDCFitWrapper.h.

86 {
87 m_t0Min = t0Min;
88 m_t0Max = t0Max;
89
90 SetT0FitLimits(t0Min, t0Max);
91 }

◆ UnconstrainFit()

void ZDCFitExpFermiFixedTaus::UnconstrainFit ( )
overridevirtual

Implements ZDCFitWrapper.

Definition at line 360 of file ZDCFitWrapper.cxx.

361{
362 std::shared_ptr<TF1> theTF1 = GetWrapperTF1();
363 theTF1->ReleaseParameter(2);
364}

Member Data Documentation

◆ m_adjTLimitsEvent

bool ZDCFitWrapper::m_adjTLimitsEvent {false}
privateinherited

Definition at line 37 of file ZDCFitWrapper.h.

37{false};

◆ m_ampMax

float ZDCFitWrapper::m_ampMax {0}
privateinherited

Definition at line 32 of file ZDCFitWrapper.h.

32{0};

◆ m_ampMin

float ZDCFitWrapper::m_ampMin {0}
privateinherited

Definition at line 31 of file ZDCFitWrapper.h.

31{0};

◆ m_expFermiFunc

std::shared_ptr<TF1> ZDCFitExpFermiFixedTaus::m_expFermiFunc {}
private

Definition at line 405 of file ZDCFitWrapper.h.

405{};

◆ m_finalized

bool ZDCFitWrapper::m_finalized {false}
privateinherited

Definition at line 42 of file ZDCFitWrapper.h.

42{false};

◆ m_norm

float ZDCFitExpFermiFixedTaus::m_norm {0}
private

Definition at line 402 of file ZDCFitWrapper.h.

402{0};

◆ m_t0Max

float ZDCFitWrapper::m_t0Max {0}
privateinherited

Definition at line 35 of file ZDCFitWrapper.h.

35{0};

◆ m_t0Min

float ZDCFitWrapper::m_t0Min {0}
privateinherited

Definition at line 34 of file ZDCFitWrapper.h.

34{0};

◆ m_tau1

float ZDCFitExpFermiFixedTaus::m_tau1 {0}
private

Definition at line 399 of file ZDCFitWrapper.h.

399{0};

◆ m_tau2

float ZDCFitExpFermiFixedTaus::m_tau2 {0}
private

Definition at line 400 of file ZDCFitWrapper.h.

400{0};

◆ m_tempTmax

float ZDCFitWrapper::m_tempTmax {0}
privateinherited

Definition at line 40 of file ZDCFitWrapper.h.

40{0};

◆ m_tempTmin

float ZDCFitWrapper::m_tempTmin {0}
privateinherited

Definition at line 39 of file ZDCFitWrapper.h.

39{0};

◆ m_timeCorr

float ZDCFitExpFermiFixedTaus::m_timeCorr {0}
private

Definition at line 403 of file ZDCFitWrapper.h.

403{0};

◆ m_tmax

float ZDCFitWrapper::m_tmax {0}
privateinherited

Definition at line 29 of file ZDCFitWrapper.h.

29{0};

◆ m_tmin

float ZDCFitWrapper::m_tmin {0}
privateinherited

Definition at line 28 of file ZDCFitWrapper.h.

28{0};

◆ m_tminAdjust

float ZDCFitWrapper::m_tminAdjust {0}
privateinherited

Definition at line 38 of file ZDCFitWrapper.h.

38{0};

◆ m_wrapperTF1

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

Definition at line 26 of file ZDCFitWrapper.h.

26{};
The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.17.0