Histogram and fitter class for drift time and pulsehight spectra The rising slope is fitted by a fermi-function: \(f(t)=r_{u,0} + \frac{A}{1+exp((t_0 - t)/T_0}\). More...

#include <T0MTHistos.h>

Collaboration diagram for MuonCalib::T0MTHistos:

Classes
class	Slice

Public Member Functions
	T0MTHistos ()=default
	Default Constructor.
	T0MTHistos (int id, const T0MTSettings settings, const char hname=nullptr)
	Initializing constructor.
void	Initialize (int id, const T0MTSettings settings, const char hname=nullptr)
	Initialize class.
TH1F *	GetTSpec () const
	get drift time spectrum
void	SetTSpec (int id, TH1F spec, const T0MTSettings settings, bool copy_spec=true)
	set the pointer of the drift-time spectrum to an existing spectrum.
void	FillT (double t)
	fill drift time spectrum
bool	FitT0 ()
	Perform t0-fit Returns true if fit is successfull.
bool	FitTmax ()
	Performs tmax-fit Returns true if fit is successfull.
int	Id () const
	return tube id
bool	T0Ok () const
	returns true if t0-fit was successfull
int	StatusCode () const
	returns status code - the status code applies only to the t0 fit
const TF1 *	GetT0Function () const
	returns function fitted to the riding edge of the spectrum
TF1 *	GetT0Function ()
	returns function fitted to the riding edge of the spectrum
bool	TmaxOk () const
	returns true if tmax-fir was successfull
const TF1 *	GetTMaxFunction () const
	returns function fitted to the riding edge of the spectrum
TF1 *	GetTMaxFunctionNC () const
	returns function fitted to the riding edge of the spectrum
double	T0Chi2 () const
	returns t0 chi2

Static Public Attributes
static constexpr int	N_T0_FIT_PAR = 4
	number of parameters in t0 fit
static constexpr int	T0_PAR_NR_T0 = 0
	parameter numbers in t0 fit
static constexpr int	T0_PAR_NR_T = 1
static constexpr int	T0_PAR_NR_BACK = 2
static constexpr int	T0_PAR_NR_A = 3
static constexpr int	N_TMAX_FIT_PAR = 6
	number of parameters for tmax fit
static constexpr int	TMAX_PAR_NR_TMAX = 0
	parameters numbers for tmax fit
static constexpr int	TMAX_PAR_NR_T = 1
static constexpr int	TMAX_PAR_NR_BACK = 2
static constexpr int	TMAX_PAR_NR_A = 3
static constexpr int	TMAX_PAR_NR_B = 4
static constexpr int	TMAX_PAR_NR_T0 = 5

Private Member Functions
bool	NormalFit ()
	normal t0 fit
bool	T0Scramble ()
	try to get better start values from a scrambled histogram
void	TopChi2 ()
	top chi2 calculation
void	TopSlicing ()
	top slicing metyhod

Private Attributes
std::unique_ptr< TH1F >	m_time {nullptr}
	time spectrum
int	m_id {-1}
	tube id;
std::unique_ptr< TF1 >	m_t0_fermi {nullptr}
	function fitted to the riding edghe of the spectrum
bool	m_t0_ok {false}
	is true if t0 fit was successful
int	m_status_code {99}
	status code for t0 fit (0 ok, 1 not fitted, 2 low statistics, 3 failed)
std::unique_ptr< TF1 >	m_tmax_fermi {nullptr}
	function fitted to the falling edge of the spectrum
double	m_tmax_ok {FLT_MAX}
	is true if tmax fit was successful
TDirectory *	m_dir {nullptr}
	TDirectory where debug and result histograms are stored.
const T0MTSettings *	m_settings {nullptr}
	Pointer to settings class.
double	m_chi2 {FLT_MAX}
	chi2/NDF value

Detailed Description

Histogram and fitter class for drift time and pulsehight spectra The rising slope is fitted by a fermi-function: \(f(t)=r_{u,0} + \frac{A}{1+exp((t_0 - t)/T_0}\).

The falling slope is fitted by \(g(t)=r_{u,max} + \frac{e(t,a,b)}{1+exp((t_{max} - t)/T_{max}}\) where \(e(t, a, b) = a e^{bt}\). A pattern recognition determines the fit ranges and the parameters \(r_u\), \(A\), \(a\) and \(b\).

Author: Felix.nosp@m..Rau.nosp@m.scher.nosp@m.@Phy.nosp@m.sik.U.nosp@m.ni-M.nosp@m.uench.nosp@m.en.D.nosp@m.e

Date: June 2006

}

Definition at line 40 of file T0MTHistos.h.

Constructor & Destructor Documentation

◆ T0MTHistos() [1/2]

MuonCalib::T0MTHistos::T0MTHistos ( )

default

Default Constructor.

◆ T0MTHistos() [2/2]

MuonCalib::T0MTHistos::T0MTHistos	(	int	id,
		const T0MTSettings *	settings,
		const char *	hname = nullptr )

inline

Initializing constructor.

Parameters

id	tube id
settings	t0-fit settings: settings will be asked about histogram binning

Definition at line 50 of file T0MTHistos.h.

                                                                                      : m_t0_fermi(nullptr), m_tmax_fermi(nullptr) {
            Initialize(id, settings, hname);
        }

Member Function Documentation

◆ FillT()

void MuonCalib::T0MTHistos::FillT ( double t )

fill drift time spectrum

Parameters

t	drift time

Definition at line 56 of file T0MTHistos.cxx.

56{ m_time->Fill(static_cast<Axis_t>(t)); }

MuonCalib::T0MTHistos::m_time

std::unique_ptr< TH1F > m_time

time spectrum

Definition T0MTHistos.h:119

◆ FitT0()

bool MuonCalib::T0MTHistos::FitT0 ( )

Perform t0-fit Returns true if fit is successfull.

Definition at line 114 of file T0MTHistos.cxx.

                           {
        if (m_time->GetEntries() < 1000) {
            m_status_code = 2;
            return false;
        }
        if (!NormalFit()) return false;
        if (m_time->GetEntries() < 10000) {
            m_status_code = 0;
            return true;
        }
        if (m_settings->T0Settings()->ScrambleThreshold() > 0 && m_chi2 > m_settings->T0Settings()->ScrambleThreshold()) {
            if (!T0Scramble()) {
                m_status_code = 3;
                return false;
            }
        }
        if (m_settings->T0Settings()->SlicingThreshold() > 0 && m_chi2 > m_settings->T0Settings()->SlicingThreshold()) { TopSlicing(); }
#ifndef NDEBUG
        MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << m_time->GetName() << " " << m_chi2 << endmsg;
#endif
        m_status_code = 0;
        return true;
    }  // end T0MTHistos::FitT0

◆ FitTmax()

bool MuonCalib::T0MTHistos::FitTmax ( )

Performs tmax-fit Returns true if fit is successfull.

Definition at line 142 of file T0MTHistos.cxx.

                             {
        TDirectory *cwd = gDirectory;
        if (m_dir != nullptr) m_dir->cd();
        if (!m_time) {
            MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
            log << MSG::WARNING << "T0MTHistos::FitTmax: Class is not initialized!" << endmsg;
            m_tmax_ok = false;
            cwd->cd();
            return false;
        }
        // check if t0-fit was successfull t0 is needed for tmax-pattern recognition
        if (!m_t0_ok || m_t0_fermi == nullptr) {
            MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
            log << MSG::WARNING << "T0MTHistos::FitTmax for tube " << m_id << ": No valid t0-value!" << endmsg;
            m_tmax_ok = false;
            cwd->cd();
            return false;
        }
        // Create pattern Recognition Class
        MTTmaxPatternRecognition rec;
        // perform pattern recognition
        if (!rec.Initialize(m_time.get(), m_t0_fermi->GetParameter(T0_PAR_NR_T0), m_settings)) {
            MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
            log << MSG::WARNING << "T0MTHistos::FitTmax for tube " << m_id << ": Pattern recognition failed!" << endmsg;
            m_tmax_ok = false;
            cwd->cd();
            return false;
        }
        // create function object
        char buffer[100];
        snprintf(buffer, 100, "mt_tmax_fermi");
        if (!m_tmax_fermi) {
            m_tmax_fermi = std::make_unique<TF1>(buffer, mt_tmax_fermi, rec.GetFitRangeMin(), rec.GetFitRangeMax(), N_TMAX_FIT_PAR);
            // set parameter names
            m_tmax_fermi->SetParName(TMAX_PAR_NR_TMAX, "t_{max}");
            m_tmax_fermi->SetParName(TMAX_PAR_NR_T, "T");
            m_tmax_fermi->SetParName(TMAX_PAR_NR_BACK, "r_{b}");
            m_tmax_fermi->SetParName(TMAX_PAR_NR_A, "a");
            m_tmax_fermi->SetParName(TMAX_PAR_NR_B, "b");
            // set fixed values
            m_tmax_fermi->FixParameter(TMAX_PAR_NR_BACK, rec.GetBackground());
            m_tmax_fermi->FixParameter(TMAX_PAR_NR_A, rec.GetA());
            m_tmax_fermi->FixParameter(TMAX_PAR_NR_B, rec.GetB());
            m_tmax_fermi->FixParameter(TMAX_PAR_NR_T0, m_t0_fermi->GetParameter(T0_PAR_NR_T0));
            // set start values
            m_tmax_fermi->SetParameter(TMAX_PAR_NR_TMAX, rec.GetEstimatedTMax());
            m_tmax_fermi->SetParameter(TMAX_PAR_NR_T, 3.0);
        }
        // perform fit
        if (m_dir != nullptr) {
            m_tmax_fermi->SetLineColor(3);
            m_tmax_fermi->Write();
        }
        m_tmax_fermi->SetLineColor(4);
        std::string fitopt("LR");
        if (m_settings->VerboseLevel() == 0) fitopt += "Q";
        if (m_settings->AddFitfun()) {
            fitopt += "+";
        } else {
            fitopt += "0";
        }
        m_time->Fit(m_tmax_fermi.get(), fitopt.c_str());
        m_tmax_ok = true;
        cwd->cd();
        return true;
    }  // end T0MTHistos::FitTmax

◆ GetT0Function() [1/2]

TF1 * MuonCalib::T0MTHistos::GetT0Function ( )

inline

returns function fitted to the riding edge of the spectrum

Definition at line 106 of file T0MTHistos.h.

106{ return m_t0_fermi.get(); }

◆ GetT0Function() [2/2]

const TF1 * MuonCalib::T0MTHistos::GetT0Function ( ) const

inline

returns function fitted to the riding edge of the spectrum

Definition at line 104 of file T0MTHistos.h.

104{ return m_t0_fermi.get(); }

◆ GetTMaxFunction()

const TF1 * MuonCalib::T0MTHistos::GetTMaxFunction ( ) const

inline

returns function fitted to the riding edge of the spectrum

Definition at line 110 of file T0MTHistos.h.

110{ return m_tmax_fermi.get(); }

◆ GetTMaxFunctionNC()

TF1 * MuonCalib::T0MTHistos::GetTMaxFunctionNC ( ) const

inline

returns function fitted to the riding edge of the spectrum

Definition at line 112 of file T0MTHistos.h.

112{ return m_tmax_fermi.get(); }

◆ GetTSpec()

TH1F * MuonCalib::T0MTHistos::GetTSpec ( ) const

inline

get drift time spectrum

Definition at line 72 of file T0MTHistos.h.

72{ return m_time.get(); }

◆ Id()

int MuonCalib::T0MTHistos::Id ( ) const

inline

return tube id

Definition at line 97 of file T0MTHistos.h.

97{ return m_id; }

◆ Initialize()

void MuonCalib::T0MTHistos::Initialize	(	int	id,
		const T0MTSettings *	settings,
		const char *	hname = nullptr )

Initialize class.

Parameters

id	tube id
settings	t0-fit settings: settings will be asked about histogram binning

Definition at line 61 of file T0MTHistos.cxx.

                                                                                       {
        m_settings = settings;
#ifndef NDEBUG
        MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "T0MTHistos::Initialize: called" << endmsg;
#endif
        char buf[100];
        if (hname == nullptr)
            snprintf(buf, 100, "t_spec_%d", id);
        else
            snprintf(buf, 100, "t_spec_%s", hname);
#ifndef NDEBUG
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "directory=" << gDirectory->GetName() << endmsg;
#endif
        m_time = std::make_unique<TH1F>(buf, "", settings->NBinsTime(), settings->TimeMin(), settings->TimeMax());
        m_id = id;
        if (settings->DrawDebugGraphs()) {
#ifndef NDEBUG
            if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "T0MTHistos::Initialize: debug directory created" << endmsg;
#endif
            TDirectory *cwd = gDirectory;
            snprintf(buf, 100, "t0_tmax_dir_%d", id);
            m_dir = gDirectory->mkdir(buf, buf);
            cwd->cd();
        } else {
            m_dir = nullptr;
        }
        m_t0_ok = false;
        m_tmax_ok = false;
    }  // end T0MTHistos::Initialize

◆ NormalFit()

bool MuonCalib::T0MTHistos::NormalFit ( )

private

normal t0 fit

Definition at line 212 of file T0MTHistos.cxx.

                               {
#ifndef NDEBUG
        MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "T0MTHistos::FitT0(): called" << endmsg;
#endif
        TDirectory *cwd = gDirectory;
        if (m_dir) m_dir->cd();
        // check if class is initialized
        if (!m_time) {
#ifdef NDEBUG
            MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
#endif
            log << MSG::WARNING << "T0MTHistos::FitT0: Class is not initialized!" << endmsg;
            m_t0_ok = false;
            cwd->cd();
            m_status_code = 3;
            return false;
        }
        // create pattern recognition class
        MTT0PatternRecognition rec;
        // perform pattern recognition
        if (!rec.Initialize(m_time.get(), m_settings)) {
            m_t0_ok = false;
#ifdef NDEBUG
            MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
#endif
            log << MSG::WARNING << "T0MTHistos::FitT0 for tube " << m_id << ": Pattern recognition failed!" << endmsg;
            cwd->cd();
            m_status_code = 3;
            return false;
        }
        m_t0_ok = true;
        // create function class
        char buffer[100];
        //  sprintf(buffer,"mt_t0_fermi_%d", m_id);
        snprintf(buffer, 100, "mt_t0_fermi");
        m_t0_fermi = std::make_unique<TF1>(buffer, mt_t0_fermi, rec.GetFitRangeMin(), rec.GetFitRangeMax(), N_T0_FIT_PAR);
        // set parameter names
        m_t0_fermi->SetParName(T0_PAR_NR_T0, "t_{0}");
        m_t0_fermi->SetParName(T0_PAR_NR_T, "T");
        m_t0_fermi->SetParName(T0_PAR_NR_BACK, "r_{b}");
        m_t0_fermi->SetParName(T0_PAR_NR_A, "A");
        // set fixed values
        m_t0_fermi->FixParameter(T0_PAR_NR_BACK, rec.GetBackground());
        m_t0_fermi->SetParameter(T0_PAR_NR_A, rec.GetHeight() - rec.GetBackground());
        // set estimates as start values
        m_t0_fermi->SetParameter(T0_PAR_NR_T0, rec.GetEstimatedT0());
        // set resonable start value for T
        m_t0_fermi->SetParameter(T0_PAR_NR_T, 3.0);
        // perform fit - NOTE: The return value of the Fit function is not documented!
        if (m_dir != nullptr) {
            std::unique_ptr<TLine> ln = std::make_unique<TLine>(rec.GetFitRangeMin(), 0, rec.GetFitRangeMin(), m_time->GetMaximum());
            ln->Write("t0_range_min");
            ln = std::make_unique<TLine>(rec.GetFitRangeMax(), 0, rec.GetFitRangeMax(), m_time->GetMaximum());
            ln->Write("t0_range_max");
            m_t0_fermi->SetLineColor(3);
            m_t0_fermi->Write();
        }
        m_t0_fermi->SetLineColor(2);
        std::string fitopt("BLR");
        if (m_settings->VerboseLevel() == 0) fitopt += "Q";
        if (m_settings->AddFitfun()) {
            fitopt += "+";
        } else {
            fitopt += "0";
        }
        m_time->Fit(m_t0_fermi.get(), fitopt.c_str(), "", rec.GetFitRangeMin(), rec.GetFitRangeMax());
        if (m_settings->T0Settings()->UseTopChi2())
            TopChi2();
        else
            m_chi2 = m_time->GetFunction("mt_t0_fermi")->GetChisquare() / m_time->GetFunction("mt_t0_fermi")->GetNDF();
        m_t0_ok = true;
        cwd->cd();
        m_status_code = 0;
        return true;
    }  // end T0MTHistos::NormalFit

◆ SetTSpec()

void MuonCalib::T0MTHistos::SetTSpec	(	int	id,
		TH1F *	spec,
		const T0MTSettings *	settings,
		bool	copy_spec = true )

set the pointer of the drift-time spectrum to an existing spectrum.

This is for testapps

Parameters

id	tube id
spec	Pointer to an existing spectrum

Definition at line 95 of file T0MTHistos.cxx.

                                                                                              {
        m_settings = settings;
        TDirectory *cwd = gDirectory;
        if (copy_spec)
            m_time = std::make_unique<TH1F>(*spec);
        else
            m_time.reset(spec);
        m_id = id;
        if (settings->DrawDebugGraphs()) {
            char buffer[100];
            snprintf(buffer, 100, "t0_tmax_dir_%d", id);
            m_dir = gDirectory->mkdir(buffer, buffer);
            cwd->cd();
        }
    }  // end T0MTHistos::SetTSpec

◆ StatusCode()

int MuonCalib::T0MTHistos::StatusCode ( ) const

inline

returns status code - the status code applies only to the t0 fit

Definition at line 102 of file T0MTHistos.h.

102{ return m_status_code; }

◆ T0Chi2()

double MuonCalib::T0MTHistos::T0Chi2 ( ) const

inline

returns t0 chi2

Definition at line 114 of file T0MTHistos.h.

114{ return m_chi2; }

◆ T0Ok()

bool MuonCalib::T0MTHistos::T0Ok ( ) const

inline

returns true if t0-fit was successfull

Definition at line 100 of file T0MTHistos.h.

100{ return m_t0_ok; }

◆ T0Scramble()

bool MuonCalib::T0MTHistos::T0Scramble ( )

private

try to get better start values from a scrambled histogram

Definition at line 292 of file T0MTHistos.cxx.

                                {
#ifndef NDEBUG
        MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "Scrambling for " << m_time->GetName() << endmsg;
#endif
        std::string fitopt("BLR");
        if (m_settings->VerboseLevel() == 0) fitopt += "Q";
        if (m_settings->AddFitfun()) {
            fitopt += "+";
        } else {
            fitopt += "0";
        }
        // create scrambled histogram
        char scramhistname[100];
        snprintf(scramhistname, 100, "%s_scram", m_time->GetName());
        std::unique_ptr<TH1F> scramhist = std::make_unique<TH1F>(scramhistname, "scrambled histogram", m_time->GetSize() - 2,
                                                                 m_time->GetXaxis()->GetXmin(), m_time->GetXaxis()->GetXmax());
 
        for (int binnr = 0; binnr < m_time->GetSize(); binnr++) {
            scramhist->SetBinContent(binnr, m_time->GetBinContent(binnr) + getTLSRandomGen()->Gaus(0, m_time->GetBinError(binnr)));
            scramhist->SetBinError(binnr, m_time->GetBinError(binnr) * 1.41421356);
            if (scramhist->GetBinContent(binnr) < 0) scramhist->SetBinContent(binnr, 0);
        }
        TDirectory *cwd = gDirectory;
        if (m_dir) m_dir->cd();
        MTT0PatternRecognition scramrec;
        // perform pattern recognition
        if (!scramrec.Initialize(scramhist.get(), m_settings)) {
#ifdef NDEBUG
            MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
#endif
            log << MSG::WARNING << "T0MTHistos::FitT0 for tube " << m_id << ": Scrambed pattern recognition failed!" << endmsg;
            cwd->cd();
            return false;
        }
        char scrambuffer[100];
        snprintf(scrambuffer, 100, "scrammt_t0_fermi");
        std::unique_ptr<TF1> scramm_t0_fermi = std::make_unique<TF1>();
        m_t0_fermi->Copy(*scramm_t0_fermi);
        scramm_t0_fermi->SetName(scrambuffer);
        scramm_t0_fermi->SetRange(scramrec.GetFitRangeMin(), scramrec.GetFitRangeMax());
        // set parameter names
        scramm_t0_fermi->SetParName(T0_PAR_NR_T0, "t_{0}");
        scramm_t0_fermi->SetParName(T0_PAR_NR_T, "T");
        scramm_t0_fermi->SetParName(T0_PAR_NR_BACK, "r_{b}");
        scramm_t0_fermi->SetParName(T0_PAR_NR_A, "A");
        // set fixed values
        scramm_t0_fermi->FixParameter(T0_PAR_NR_BACK, scramrec.GetBackground());
        scramm_t0_fermi->SetParameter(T0_PAR_NR_A, scramrec.GetHeight() - scramrec.GetBackground());
        // set estimates as start values
        scramm_t0_fermi->SetParameter(T0_PAR_NR_T0, scramrec.GetEstimatedT0());
        // set resonable start value for T
        scramm_t0_fermi->SetParameter(T0_PAR_NR_T, 3.0);
        // perform fit - NOTE: The return value of the Fit function is not documented!
        scramhist->Fit(scrambuffer, fitopt.c_str(), "", scramrec.GetFitRangeMin(), scramrec.GetFitRangeMax());
        // set parameter for the new fit of the original histogram
        m_time->GetListOfFunctions()->Clear();
        // set fixed values
        m_t0_fermi->FixParameter(T0_PAR_NR_BACK, scramm_t0_fermi->GetParameter(T0_PAR_NR_BACK));
        m_t0_fermi->SetParameter(T0_PAR_NR_A, scramm_t0_fermi->GetParameter(T0_PAR_NR_A));
        // set estimates as start values
        m_t0_fermi->SetParameter(T0_PAR_NR_T0, scramm_t0_fermi->GetParameter(T0_PAR_NR_T0));
        // set resonable start value for T
        m_t0_fermi->SetParameter(T0_PAR_NR_T, scramm_t0_fermi->GetParameter(T0_PAR_NR_T));
        m_time->GetListOfFunctions()->Clear();
        m_time->Fit("mt_t0_fermi", fitopt.c_str(), "", scramrec.GetFitRangeMin(), scramrec.GetFitRangeMax());
        if (m_settings->T0Settings()->UseTopChi2())
            TopChi2();
        else
            m_chi2 = m_time->GetFunction("mt_t0_fermi")->GetChisquare() / m_time->GetFunction("mt_t0_fermi")->GetNDF();
        cwd->cd();
        return true;
    }  // end T0MTHistos::T0Scramble

◆ TmaxOk()

bool MuonCalib::T0MTHistos::TmaxOk ( ) const

inline

returns true if tmax-fir was successfull

Definition at line 108 of file T0MTHistos.h.

108{ return m_tmax_ok; }

◆ TopChi2()

void MuonCalib::T0MTHistos::TopChi2 ( )

private

top chi2 calculation

Definition at line 366 of file T0MTHistos.cxx.

                             {
        // calculate topchi2
        m_chi2 = 0;
        int topndf = 0;
        TF1 *t0_fermi = m_time->GetFunction("mt_t0_fermi");
        Double_t min, max;
        t0_fermi->GetRange(min, max);
        int startbin = m_time->FindBin(min) + 1;
        int endbin = m_time->FindBin(max) - 1;
        for (int bin = startbin; bin < endbin; bin++) {
            float measval = m_time->GetBinContent(bin);
            float funcval = t0_fermi->Eval(m_time->GetBinCenter(bin));
            float errval = m_time->GetBinError(bin);
            // take only chi2 from top part or if the  bin content is min 10 or if the function >10
            if (measval < 10 && funcval < 10 &&
                m_time->GetBinCenter(bin) < t0_fermi->GetParameter(T0_PAR_NR_T0) + 2 * t0_fermi->GetParameter(T0_PAR_NR_T))
                continue;
            if (errval == 0)
                m_chi2 += (measval - funcval) * (measval - funcval);
            else
                m_chi2 += (measval - funcval) * (measval - funcval) / (errval * errval);
            topndf++;
        }
        if (topndf != 0) m_chi2 = m_chi2 / topndf;
    }  // end T0MTHistos::TopChi2

◆ TopSlicing()

void MuonCalib::T0MTHistos::TopSlicing ( )

private

top slicing metyhod

Definition at line 392 of file T0MTHistos.cxx.

                                {
#ifndef NDEBUG
        MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "Slicing for " << m_time->GetName() << endmsg;
#endif
        // vector with slice chi2
        std::list<Slice> slice_chi2;
        Slice current;
        current.chi_2 = 0.0;
        current.n_bins = 0;
        TF1 *t0_fermi = m_time->GetFunction("mt_t0_fermi");
        current.min_bin = m_time->FindBin(t0_fermi->GetParameter(T0_PAR_NR_T0) + 2 * t0_fermi->GetParameter(T0_PAR_NR_T));
        Double_t min, max;
        t0_fermi->GetRange(min, max);
#ifndef NDEBUG
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << current.min_bin << " " << max << endmsg;
#endif
        for (int bin = m_time->FindBin(t0_fermi->GetParameter(T0_PAR_NR_T0) + 2 * t0_fermi->GetParameter(T0_PAR_NR_T));
             bin < m_time->FindBin(max) - 1; bin++) {
            if (current.n_bins == 10) {
#ifndef NDEBUG
                if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << current.chi_2 / current.n_bins << endmsg;
#endif
                current.max_bin = bin;
                slice_chi2.push_back(current);
                current.chi_2 = 0.0;
                current.min_bin = bin;
                current.n_bins = 0;
            }
            double measval = m_time->GetBinContent(bin);
            double funcval = t0_fermi->Eval(m_time->GetBinCenter(bin));
            double errval = m_time->GetBinError(bin);
            if (errval == 0)
                current.chi_2 += std::pow(measval - funcval, 2.0);
            else
                current.chi_2 += std::pow(measval - funcval, 2.0) / std::pow(errval, 2);
            current.n_bins++;
        }
#ifndef NDEBUG
        if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "number of slices: " << slice_chi2.size() << endmsg;
#endif
        std::list<Slice>::iterator it = slice_chi2.end();
        do {
            --it;
            if (it == slice_chi2.begin()) {
#ifndef NDEBUG
                if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "No gain in slicing!" << endmsg;
#endif
                return;
            }
        } while (it->chi_2 / static_cast<double>(it->n_bins) > 3);
        max = m_time->GetBinCenter(it->min_bin);
        m_time->GetListOfFunctions()->Clear();
        std::string fitopt("BLR");
        if (m_settings->VerboseLevel() == 0) fitopt += "Q";
        if (m_settings->AddFitfun()) {
            fitopt += "+";
        } else {
            fitopt += "0";
        }
        m_time->Fit("mt_t0_fermi", fitopt.c_str(), "", min, max);
        if (m_settings->T0Settings()->UseTopChi2())
            TopChi2();
        else
            m_chi2 = m_time->GetFunction("mt_t0_fermi")->GetChisquare() / m_time->GetFunction("mt_t0_fermi")->GetNDF();
    }  // end T0MTHistos::TopSlicing

Member Data Documentation

◆ m_chi2

double MuonCalib::T0MTHistos::m_chi2 {FLT_MAX}

private

chi2/NDF value

Definition at line 137 of file T0MTHistos.h.

137{FLT_MAX};

◆ m_dir

TDirectory* MuonCalib::T0MTHistos::m_dir {nullptr}

private

TDirectory where debug and result histograms are stored.

Definition at line 133 of file T0MTHistos.h.

133{nullptr};

◆ m_id

int MuonCalib::T0MTHistos::m_id {-1}

private

tube id;

Definition at line 121 of file T0MTHistos.h.

121{-1};

◆ m_settings

const T0MTSettings* MuonCalib::T0MTHistos::m_settings {nullptr}

private

Pointer to settings class.

Definition at line 135 of file T0MTHistos.h.

135{nullptr};

◆ m_status_code

int MuonCalib::T0MTHistos::m_status_code {99}

private

status code for t0 fit (0 ok, 1 not fitted, 2 low statistics, 3 failed)

Definition at line 127 of file T0MTHistos.h.

127{99};

◆ m_t0_fermi

std::unique_ptr<TF1> MuonCalib::T0MTHistos::m_t0_fermi {nullptr}

private

function fitted to the riding edghe of the spectrum

Definition at line 123 of file T0MTHistos.h.

123{nullptr};

◆ m_t0_ok

bool MuonCalib::T0MTHistos::m_t0_ok {false}

private

is true if t0 fit was successful

Definition at line 125 of file T0MTHistos.h.

125{false};

◆ m_time

std::unique_ptr<TH1F> MuonCalib::T0MTHistos::m_time {nullptr}

private

time spectrum

Definition at line 119 of file T0MTHistos.h.

119{nullptr};

◆ m_tmax_fermi

std::unique_ptr<TF1> MuonCalib::T0MTHistos::m_tmax_fermi {nullptr}

private

function fitted to the falling edge of the spectrum

Definition at line 129 of file T0MTHistos.h.

129{nullptr};

◆ m_tmax_ok

double MuonCalib::T0MTHistos::m_tmax_ok {FLT_MAX}

private

is true if tmax fit was successful

Definition at line 131 of file T0MTHistos.h.

131{FLT_MAX};

◆ N_T0_FIT_PAR

int MuonCalib::T0MTHistos::N_T0_FIT_PAR = 4

staticconstexpr

number of parameters in t0 fit

Definition at line 55 of file T0MTHistos.h.

◆ N_TMAX_FIT_PAR

int MuonCalib::T0MTHistos::N_TMAX_FIT_PAR = 6

staticconstexpr

number of parameters for tmax fit

Definition at line 59 of file T0MTHistos.h.

◆ T0_PAR_NR_A

int MuonCalib::T0MTHistos::T0_PAR_NR_A = 3

static

Definition at line 57 of file T0MTHistos.h.

◆ T0_PAR_NR_BACK

int MuonCalib::T0MTHistos::T0_PAR_NR_BACK = 2

static

Definition at line 57 of file T0MTHistos.h.

◆ T0_PAR_NR_T

int MuonCalib::T0MTHistos::T0_PAR_NR_T = 1

static

Definition at line 57 of file T0MTHistos.h.

◆ T0_PAR_NR_T0

int MuonCalib::T0MTHistos::T0_PAR_NR_T0 = 0

staticconstexpr

parameter numbers in t0 fit

Definition at line 57 of file T0MTHistos.h.

◆ TMAX_PAR_NR_A

int MuonCalib::T0MTHistos::TMAX_PAR_NR_A = 3

static

Definition at line 61 of file T0MTHistos.h.

◆ TMAX_PAR_NR_B

int MuonCalib::T0MTHistos::TMAX_PAR_NR_B = 4

static

Definition at line 61 of file T0MTHistos.h.

◆ TMAX_PAR_NR_BACK

int MuonCalib::T0MTHistos::TMAX_PAR_NR_BACK = 2

static

Definition at line 61 of file T0MTHistos.h.

◆ TMAX_PAR_NR_T

int MuonCalib::T0MTHistos::TMAX_PAR_NR_T = 1

static

Definition at line 61 of file T0MTHistos.h.

◆ TMAX_PAR_NR_T0

int MuonCalib::T0MTHistos::TMAX_PAR_NR_T0 = 5

static

Definition at line 62 of file T0MTHistos.h.

◆ TMAX_PAR_NR_TMAX

int MuonCalib::T0MTHistos::TMAX_PAR_NR_TMAX = 0

staticconstexpr

parameters numbers for tmax fit

Definition at line 61 of file T0MTHistos.h.

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

Classes

Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ T0MTHistos() [1/2]

◆ T0MTHistos() [2/2]

Member Function Documentation

◆ FillT()

◆ FitT0()

◆ FitTmax()

◆ GetT0Function() [1/2]

◆ GetT0Function() [2/2]

◆ GetTMaxFunction()

◆ GetTMaxFunctionNC()

◆ GetTSpec()

◆ Id()

◆ Initialize()

◆ NormalFit()

◆ SetTSpec()

◆ StatusCode()

◆ T0Chi2()

◆ T0Ok()

◆ T0Scramble()

◆ TmaxOk()

◆ TopChi2()

◆ TopSlicing()

Member Data Documentation

◆ m_chi2

◆ m_dir

◆ m_id

◆ m_settings

◆ m_status_code

◆ m_t0_fermi

◆ m_t0_ok

◆ m_time

◆ m_tmax_fermi

◆ m_tmax_ok

◆ N_T0_FIT_PAR

◆ N_TMAX_FIT_PAR

◆ T0_PAR_NR_A

◆ T0_PAR_NR_BACK

◆ T0_PAR_NR_T

◆ T0_PAR_NR_T0

◆ TMAX_PAR_NR_A

◆ TMAX_PAR_NR_B

◆ TMAX_PAR_NR_BACK

◆ TMAX_PAR_NR_T

◆ TMAX_PAR_NR_T0

◆ TMAX_PAR_NR_TMAX