MuonCalib::MTTmaxPatternRecognition Class Reference

#include <MTTmaxPatternRecognition.h>

Collaboration diagram for MuonCalib::MTTmaxPatternRecognition:

Public Member Functions
	MTTmaxPatternRecognition ()
	Default constructor. More...
bool	Initialize (TH1F *hist, double t0, const T0MTSettings *settings)
	Initialize class. More...
double	GetBackground () const
	get the background level More...
double	GetA () const
	get parameter a in exp-function representing the end of the spectrum More...
double	GetB () const
	get parameter a in exp-function representing the end of the spectrum More...
double	GetEstimatedTMax () const
	get estimated t0 More...
double	GetFitRangeMin () const
	get fit range More...
double	GetFitRangeMax () const
	get fit range More...
bool	GetError () const
	return error flag More...

Private Member Functions
bool	estimate_background (TH1F *hist)
	estimates the background level More...
bool	estimate_height (TH1F *hist)
	estimates the height of the spectrum. More...

Private Attributes
const T0MTSettingsTMax *	m_settings
	settings More...
bool	m_draw_debug_graph
double	m_background
	background level More...
double	m_a
	parameters a, b in a*exp(b*x) More...
double	m_b
double	m_fit_min
	fit range More...
double	m_fit_max
double	m_tmax_est
	t0 estimate More...
VariableBinwidthHistogram	m_vbh
	Variable binwidth histogram. More...
double	m_t0
	t0 More...
bool	m_error
	error flag More...

Detailed Description

Definition at line 36 of file MTTmaxPatternRecognition.h.

Constructor & Destructor Documentation

MTTmaxPatternRecognition()

MuonCalib::MTTmaxPatternRecognition::MTTmaxPatternRecognition ( )

inline

Default constructor.

Definition at line 40 of file MTTmaxPatternRecognition.h.

                                          :
            m_settings(NULL),
            m_draw_debug_graph(false),
            m_background(0.),
            m_a(0.),
            m_b(0.),
            m_fit_min(0.),
            m_fit_max(0.),
            m_tmax_est(0.),
            m_t0(0.),
            m_error(false) {}

Member Function Documentation

estimate_background()

bool MuonCalib::MTTmaxPatternRecognition::estimate_background ( TH1F *  hist )

private

estimates the background level

Parameters

hist	input histogram

Definition at line 39 of file MTTmaxPatternRecognition.cxx.

                                                                 {
        // smooth VBH in several steps
        for (double i = 1.0; i <= 4.0; i *= 2.0) {
            for (int j = 0; j < 2; j++)
                if (!m_vbh.Smooth(hist->GetBinWidth(1) / i)) {
                    m_error = true;
                    return false;
                }
        }
        if (m_draw_debug_graph) {
            m_vbh.DenistyGraph()->Write("tmax_density");
            m_vbh.BinWidthGraph()->Write("tmax_binwidth");
            m_vbh.DiffDensityGraph()->Write("tmax_diffdensity");
            m_vbh.DiffBinwidthGraph()->Write("tmax_diffbinwidth");
        }
        // find falling edge - search for steepenst slope in a range of 600 to 800 ns after the t0
        double peak_falling(0.0), peak_falling_slope(-9e9);
        for (unsigned int i = 0; i < m_vbh.GetNumberOfBins() - 1; i++) {
            const VariableBinwidthHistogramBin &bin1(m_vbh.GetBin(i)), &bin2(m_vbh.GetBin(i + 1));
            // select region
            if (bin1.Center() < m_t0 + m_settings->EndMin() || bin1.Center() > m_t0 + m_settings->EndMax()) continue;
            if (bin2.Width() - bin1.Width() > peak_falling_slope && bin2.Right() < hist->GetBinLowEdge(hist->GetNbinsX()) - 10) {
                peak_falling = bin1.Right();
                peak_falling_slope = bin2.Width() - bin1.Width();
            }
        }
        // check region size
        int firstbin = hist->FindBin(peak_falling + m_settings->DistBackground());
        int lastbin = hist->FindBin(m_vbh.GetBin(m_vbh.GetNumberOfBins() - 1).Right());
        if (lastbin - firstbin < m_settings->MinBackgroundBins()) lastbin = hist->GetNbinsX();
        if (lastbin - firstbin < m_settings->MinBackgroundBins()) {
            MsgStream log(Athena::getMessageSvc(), "MTTmaxPatternRecognition");
            log << MSG::WARNING << "estimate_background() - Falling edge is to glose to upper histogram range!" << endmsg;
            m_error = true;
            return false;
        }
        // calcul;ate mean
        m_background = 0.0;
        double n_bins = 0.0;
        if (lastbin > hist->GetNbinsX()) lastbin = hist->GetNbinsX();
        for (int i = firstbin; i <= lastbin; i++) {
            m_background += hist->GetBinContent(i);
            n_bins++;
        }
        // mark selected range
        if (m_draw_debug_graph) {
            (new TLine(hist->GetBinCenter(firstbin), 0, hist->GetBinCenter(firstbin), hist->GetMaximum()))->Write("tmax_back_left");
            (new TLine(hist->GetBinCenter(lastbin), 0, hist->GetBinCenter(lastbin), hist->GetMaximum()))->Write("tmax_back_right");
        }
 
        m_background /= n_bins;
        m_tmax_est = peak_falling;
        m_fit_max = hist->GetBinCenter(lastbin);
        return true;
    }

estimate_height()

bool MuonCalib::MTTmaxPatternRecognition::estimate_height ( TH1F *  hist )

private

estimates the height of the spectrum.

The height is representet by an exponetial function a*exp(b*t). It returns the upper end of the height region.

Parameters

hist	input histogram

Definition at line 99 of file MTTmaxPatternRecognition.cxx.

                                                             {
        // get start value for fit
        Double_t left = m_tmax_est - m_settings->DistAB() - m_settings->WidthAB(), right = m_tmax_est - m_settings->DistAB();
        // fit
        TF1* fun = new TF1("myexp", myexp, left, right, 4);
        fun->FixParameter(2, m_background);
        fun->FixParameter(3, m_t0);
        hist->Fit("myexp", "Q0+", "", left, right);
        // store parameters
        m_a = fun->GetParameter(0);
        m_b = fun->GetParameter(1);
        m_fit_min = m_tmax_est - m_settings->DistAB() - m_settings->WidthAB();
        // mark selected range
        if (m_draw_debug_graph) {
            (new TLine(left, 0, left, hist->GetMaximum()))->Write("tmax_height_left");
            (new TLine(right, 0, right, hist->GetMaximum()))->Write("tmax_height_right");
        }
        return true;
    }

GetA()

double MuonCalib::MTTmaxPatternRecognition::GetA ( ) const

inline

get parameter a in exp-function representing the end of the spectrum

Definition at line 63 of file MTTmaxPatternRecognition.h.

{ return m_a; }

GetB()

double MuonCalib::MTTmaxPatternRecognition::GetB ( ) const

inline

get parameter a in exp-function representing the end of the spectrum

Definition at line 66 of file MTTmaxPatternRecognition.h.

{ return m_b; }

GetBackground()

double MuonCalib::MTTmaxPatternRecognition::GetBackground ( ) const

inline

get the background level

Definition at line 60 of file MTTmaxPatternRecognition.h.

{ return m_background; }

GetError()

bool MuonCalib::MTTmaxPatternRecognition::GetError ( ) const

inline

return error flag

Definition at line 78 of file MTTmaxPatternRecognition.h.

{ return m_error; }

GetEstimatedTMax()

double MuonCalib::MTTmaxPatternRecognition::GetEstimatedTMax ( ) const

inline

get estimated t0

Definition at line 69 of file MTTmaxPatternRecognition.h.

{ return m_tmax_est; }

GetFitRangeMax()

double MuonCalib::MTTmaxPatternRecognition::GetFitRangeMax ( ) const

inline

get fit range

Definition at line 75 of file MTTmaxPatternRecognition.h.

{ return m_fit_max; }

GetFitRangeMin()

double MuonCalib::MTTmaxPatternRecognition::GetFitRangeMin ( ) const

inline

get fit range

Definition at line 72 of file MTTmaxPatternRecognition.h.

{ return m_fit_min; }

Initialize()

bool MuonCalib::MTTmaxPatternRecognition::Initialize	(	TH1F *	hist,
		double	t0,
		const T0MTSettings *	settings
	)

Initialize class.

Parameters

hist	Histogram which is to be fitted
t0	comes from the t0 fit

Definition at line 20 of file MTTmaxPatternRecognition.cxx.

                                                                                                 {
        m_settings = settings->TMaxSettings();
        m_draw_debug_graph = settings->DrawDebugGraphs();
        m_error = false;
        m_t0 = t0;
        m_error = !m_vbh.Initialize(hist, m_settings->VBHBinContent(), m_settings->MaxBinwidth(), t0 + m_settings->VBHLow());
        if (m_error || m_vbh.GetNumberOfBins() < 2) {
            MsgStream log(Athena::getMessageSvc(), "MTTmaxPatternRecognition");
            log << MSG::WARNING << "Initialize() - Initialization of VBH failed!" << endmsg;
        }
        if (!m_error) estimate_background(hist);
        if (!m_error) estimate_height(hist);
        return !m_error;
    }

Member Data Documentation

m_a

double MuonCalib::MTTmaxPatternRecognition::m_a

private

parameters a, b in a*exp(b*x)

Definition at line 90 of file MTTmaxPatternRecognition.h.

m_b

double MuonCalib::MTTmaxPatternRecognition::m_b

private

Definition at line 90 of file MTTmaxPatternRecognition.h.

m_background

double MuonCalib::MTTmaxPatternRecognition::m_background

private

background level

Definition at line 87 of file MTTmaxPatternRecognition.h.

m_draw_debug_graph

bool MuonCalib::MTTmaxPatternRecognition::m_draw_debug_graph

private

Definition at line 84 of file MTTmaxPatternRecognition.h.

m_error

bool MuonCalib::MTTmaxPatternRecognition::m_error

private

error flag

Definition at line 105 of file MTTmaxPatternRecognition.h.

m_fit_max

double MuonCalib::MTTmaxPatternRecognition::m_fit_max

private

Definition at line 93 of file MTTmaxPatternRecognition.h.

m_fit_min

double MuonCalib::MTTmaxPatternRecognition::m_fit_min

private

fit range

Definition at line 93 of file MTTmaxPatternRecognition.h.

m_settings

const T0MTSettingsTMax* MuonCalib::MTTmaxPatternRecognition::m_settings

private

settings

Definition at line 83 of file MTTmaxPatternRecognition.h.

m_t0

double MuonCalib::MTTmaxPatternRecognition::m_t0

private

t0

Definition at line 102 of file MTTmaxPatternRecognition.h.

m_tmax_est

double MuonCalib::MTTmaxPatternRecognition::m_tmax_est

private

t0 estimate

Definition at line 96 of file MTTmaxPatternRecognition.h.

m_vbh

VariableBinwidthHistogram MuonCalib::MTTmaxPatternRecognition::m_vbh

private

Variable binwidth histogram.

Definition at line 99 of file MTTmaxPatternRecognition.h.

---

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

• MTTmaxPatternRecognition.h
• MTTmaxPatternRecognition.cxx