T0MTHistos.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
//   T0MTHistos.h
//   Header file for class T0MTHistos
// (c) ATLAS Detector software
//
 
#ifndef T0MTHISTOS_H
#define T0MTHISTOS_H
 
// c - c++
#include <iostream>
 
// root
#include "TF1.h"
#include "TH1.h"
class TDirectory;
 
// this
#include "T0MTSettings.h"
 
namespace MuonCalib {

    class T0MTHistos {
    public:
        //---------------------------constructor----------------------------------------
        T0MTHistos() = default;
        T0MTHistos(int id, const T0MTSettings *settings, const char *hname = nullptr) : m_t0_fermi(nullptr), m_tmax_fermi(nullptr) {
            Initialize(id, settings, hname);
        }
        //---------------------------static constants-----------------------------------
        static constexpr int N_T0_FIT_PAR = 4;
        static constexpr int T0_PAR_NR_T0 = 0, T0_PAR_NR_T = 1, T0_PAR_NR_BACK = 2, T0_PAR_NR_A = 3;
        static constexpr int N_TMAX_FIT_PAR = 6;
        static constexpr int TMAX_PAR_NR_TMAX = 0, TMAX_PAR_NR_T = 1, TMAX_PAR_NR_BACK = 2, TMAX_PAR_NR_A = 3, TMAX_PAR_NR_B = 4,
                         TMAX_PAR_NR_T0 = 5;
        //---------------------------public member functions----------------------------
        void Initialize(int id, const T0MTSettings *settings, const char *hname = nullptr);

        inline TH1F *GetTSpec() const { return m_time.get(); }

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

        void FillT(double t);

        bool FitT0();

        bool FitTmax();

        inline int Id() const { return m_id; }

        inline bool T0Ok() const { return m_t0_ok; }
        inline int StatusCode() const { return m_status_code; }
        inline const TF1 *GetT0Function() const { return m_t0_fermi.get(); }
        inline TF1 *GetT0Function() { return m_t0_fermi.get(); }
        inline bool TmaxOk() const { return m_tmax_ok; }
        inline const TF1 *GetTMaxFunction() const { return m_tmax_fermi.get(); }
        inline TF1 *GetTMaxFunctionNC() const { return m_tmax_fermi.get(); }
        inline double T0Chi2() const { return m_chi2; }
 
    private:
        //---------------------------private data members-------------------------------
        std::unique_ptr<TH1F> m_time{nullptr};
        int m_id{-1};
        std::unique_ptr<TF1> m_t0_fermi{nullptr};
        bool m_t0_ok{false};
        int m_status_code{99};
        std::unique_ptr<TF1> m_tmax_fermi{nullptr};
        double m_tmax_ok{FLT_MAX};
        TDirectory *m_dir{nullptr};
        const T0MTSettings *m_settings{nullptr};
        double m_chi2{FLT_MAX};
        bool NormalFit();
        bool T0Scramble();
        void TopChi2();
        void TopSlicing();
        class Slice {
        public:
            double chi_2{FLT_MAX};
            int min_bin{0};
            int max_bin{0};
            int n_bins{0};
        };
 
    };  // end  class T0MTHistos
 
}  // end namespace MuonCalib
 
#endif