RtCalibrationIntegration.cxx File Reference

#include "MdtCalibRt/RtCalibrationIntegration.h"
#include <algorithm>
#include <fstream>
#include <iostream>
#include "AthenaKernel/getMessageSvc.h"
#include "GaudiKernel/MsgStream.h"
#include "MdtCalibData/RtFromPoints.h"
#include "MdtCalibRt/RtParabolicExtrapolation.h"
#include "MdtCalibT0/T0MTHistos.h"
#include "TF1.h"
#include "TList.h"

Go to the source code of this file.

Functions
Double_t	slope_function_C (Double_t *x, Double_t *par)
void	update_parameter_on_mttmax (TH1 *h, TF1 *f, const float &b, const float &T, const T0MTSettingsTMax &tmax_settings)

Function Documentation

slope_function_C()

Double_t slope_function_C ( Double_t *  x, Double_t *  par  )

inline

Definition at line 21 of file RtCalibrationIntegration.cxx.

                                                             {
    Double_t &t(x[0]);
    Double_t &a(par[0]), &b(par[1]), &t_0(par[2]), &back(par[3]);
    return back + std::exp(a + b * (t - t_0));
}

update_parameter_on_mttmax()

void update_parameter_on_mttmax ( TH1 *  h, TF1 *  f, const float &  b, const float &  T, const T0MTSettingsTMax &  tmax_settings  )

inline

Definition at line 27 of file RtCalibrationIntegration.cxx.

                                                                                                                              {
    Double_t rmin, rmax;
    f->GetRange(rmin, rmax);
    std::unique_ptr<TF1> slope_function =
        std::make_unique<TF1>("slope_function", slope_function_C, rmin, rmin + tmax_settings.WidthAB(), 4);
    slope_function->SetParameter(0, f->GetParameter(T0MTHistos::TMAX_PAR_NR_A));
    slope_function->FixParameter(1, b);
    slope_function->FixParameter(2, f->GetParameter(T0MTHistos::TMAX_PAR_NR_T0));
    slope_function->FixParameter(3, f->GetParameter(T0MTHistos::TMAX_PAR_NR_BACK));
    h->Fit("slope_function", "R+", "");
    f->FixParameter(T0MTHistos::TMAX_PAR_NR_A, slope_function->GetParameter(0));
    f->FixParameter(T0MTHistos::TMAX_PAR_NR_B, b);
    f->FixParameter(T0MTHistos::TMAX_PAR_NR_T, T);
}