MuonCalib::MultilayerRtDifference Class Reference

#include <MultilayerRtDifference.h>

Collaboration diagram for MuonCalib::MultilayerRtDifference:

Public Member Functions
	MultilayerRtDifference (int min_hits, TDirectory *control_histogram_dir=NULL)
	MultilayerRtDifference (const MultilayerRtDifference &MLRTD)
virtual	~MultilayerRtDifference ()
MultilayerRtDifference &	operator= (const MultilayerRtDifference &MLRTD)
void	Fill (const MdtCalibHitBase &hit, const IRtRelation &rt_relation)
bool	DoFit (IRtRelation *rt_relation, const IMdtCalibration::MuonSegVec &seg)
const TF1 *	GetFunction () const

Private Attributes
TF1 *	m_polfun
MultilayerRtDifference_Histograms *	m_histograms
int	m_min_number_of_hits

Detailed Description

Definition at line 23 of file MultilayerRtDifference.h.

Constructor & Destructor Documentation

MultilayerRtDifference() [1/2]

MuonCalib::MultilayerRtDifference::MultilayerRtDifference	(	int	min_hits,
		TDirectory *	control_histogram_dir = NULL )

Definition at line 147 of file MultilayerRtDifference.cxx.

                                                                                                  :
        m_histograms(new MultilayerRtDifference_Histograms(control_histogram_dir)), m_min_number_of_hits(min_hits) {
        m_polfun = new TF1("polfun", MultilayerRtDifference_ScaleFunction, 4.0, 15.0, 1);
    }

MultilayerRtDifference() [2/2]

MuonCalib::MultilayerRtDifference::MultilayerRtDifference

(

const MultilayerRtDifference &

MLRTD

)

Definition at line 153 of file MultilayerRtDifference.cxx.

                                                                                      : m_min_number_of_hits(MLRTD.m_min_number_of_hits) {
        m_polfun = new TF1(*MLRTD.m_polfun);
        m_histograms = new MultilayerRtDifference_Histograms(*MLRTD.m_histograms);
    }

~MultilayerRtDifference()

MuonCalib::MultilayerRtDifference::~MultilayerRtDifference ( )

virtual

Definition at line 170 of file MultilayerRtDifference.cxx.

                                                    {
        delete m_histograms;
        delete m_polfun;
    }

Member Function Documentation

DoFit()

bool MuonCalib::MultilayerRtDifference::DoFit	(	IRtRelation *	rt_relation,
		const IMdtCalibration::MuonSegVec &	seg )

Definition at line 183 of file MultilayerRtDifference.cxx.

                                                                                                     {
        TProfile *prov_diff = m_histograms->GetProfileDiff(m_min_number_of_hits);
        if (!prov_diff) {
            MsgStream log(Athena::getMessageSvc(), "MultilayerRtDifference");
            log << MSG::WARNING << "MultilayerRtDifference::DoFit: Not enough hits!" << endmsg;
            return false;
        }
        if (prov_diff->Fit("polfun", "Q", "", 4., 15.) != 0) {
            MsgStream log(Athena::getMessageSvc(), "MultilayerRtDifference");
            log << MSG::WARNING << "MultilayerRtDifference: Fit of polinomial failed! Not updating scale!" << endmsg;
            return false;
        }
        if (std::abs(m_polfun->GetParameter(0)) < m_polfun->GetParError(0)) {
            MsgStream log(Athena::getMessageSvc(), "MultilayerRtDifference");
            log << MSG::WARNING << "MultilayerRtDifference: No Scale update needed! Scale correction: " << m_polfun->GetParameter(0)
                << " +- " << m_polfun->GetParError(0) << endmsg;
            return true;
        }
        if (!rt_relation) return true;
        float scale = m_polfun->GetParameter(0);
        if (std::abs(scale) > 2) scale *= 4;
        if (rt_relation->hasTmaxDiff()) { scale += rt_relation->GetTmaxDiff(); }
        rt_relation->SetTmaxDiff(scale);
        if (seg.empty()) return true;
        // update input segments
        for (const auto &segment : seg) {
            for (const MuonCalibSegment::MdtHitPtr &hit : segment->mdtHOT()) {
                float old_corr = hit->TemperatureTime();
                float corr = RtScaleFunction(hit->driftTime(), hit->identify().mdtMultilayer() == 2, *rt_relation);
                hit->setTemperatureTime(corr);
                hit->setDriftTime(hit->driftTime() - corr + old_corr);
                hit->setDriftRadius(rt_relation->radius(hit->driftTime()), hit->sigmaDriftRadius());
            }
        }
        return true;
    }

Fill()

void MuonCalib::MultilayerRtDifference::Fill	(	const MdtCalibHitBase &	hit,
		const IRtRelation &	rt_relation )

Definition at line 175 of file MultilayerRtDifference.cxx.

                                                                                                {
        int ml = hit.identify().mdtMultilayer() - 1;
        double r_track = std::abs(hit.signedDistanceToTrack());
        double res = std::abs(hit.driftRadius()) - r_track;
        double v = rt_relation.driftVelocity(hit.driftTime());
        m_histograms->GetResHist(ml)->Fill(r_track, res / v);
    }

GetFunction()

const TF1 * MuonCalib::MultilayerRtDifference::GetFunction ( ) const

inline

Definition at line 34 of file MultilayerRtDifference.h.

{ return m_polfun; }

operator=()

MultilayerRtDifference & MuonCalib::MultilayerRtDifference::operator=

(

const MultilayerRtDifference &

MLRTD

)

Definition at line 159 of file MultilayerRtDifference.cxx.

                                                                                               {
        if (this != &MLRTD) {
            m_min_number_of_hits = MLRTD.m_min_number_of_hits;
            delete m_polfun;
            m_polfun = new TF1(*MLRTD.m_polfun);
            delete m_histograms;
            m_histograms = new MultilayerRtDifference_Histograms(*MLRTD.m_histograms);
        }
        return *this;
    }

Member Data Documentation

m_histograms

MultilayerRtDifference_Histograms* MuonCalib::MultilayerRtDifference::m_histograms

private

Definition at line 38 of file MultilayerRtDifference.h.

m_min_number_of_hits

int MuonCalib::MultilayerRtDifference::m_min_number_of_hits

private

Definition at line 39 of file MultilayerRtDifference.h.

m_polfun

TF1* MuonCalib::MultilayerRtDifference::m_polfun

private

Definition at line 37 of file MultilayerRtDifference.h.

---

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

• MultilayerRtDifference.h
• MultilayerRtDifference.cxx