BCMTFAnalysisFacility.h

#ifndef __BCMTFANALYSISFACILITY__H
#define __BCMTFANALYSISFACILITY__H

/*
 * Copyright (C) 2007-2018, the BAT core developer team
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 * For documentation see http://mpp.mpg.de/bat
 */

// ---------------------------------------------------------
#include <string>
#include <vector>

#include <TH1D.h>

#include "../../BAT/BCLog.h"

class BCMTF;
class TTree;
class TRandom3;

// ---------------------------------------------------------
class BCMTFAnalysisFacility
{

public:

    BCMTFAnalysisFacility(BCMTF* mtf);

    ~BCMTFAnalysisFacility();

    BCMTF* GetBCMTF()
    { return fMTF; };

    void SetBCMTF(BCMTF* mtf)
    { fMTF = mtf; };

    void SetFlagMarginalize(bool flag)
    { fFlagMarginalize = flag; };

    // misc

    BCLog::LogLevel GetLogLevel()
    { return fLogLevel; };

    void SetLogLevel(BCLog::LogLevel level)
    { fLogLevel = level; };

    void PerformSingleChannelAnalyses(const std::string& dirname, const std::string& options = "");

    void PerformSingleSystematicAnalyses(const std::string& dirname, const std::string& options = "");

    void PerformCalibrationAnalysis(const std::string& dirname, const std::vector<double>& default_parameters, int index, const std::vector<double>& parametervalues, int nensembles = 1000);

    std::vector<TH1D> BuildEnsemble(const std::vector<double>& parameters, const std::string& options = "");

    TTree* BuildEnsembles(const std::vector<double>& parameters, int nensembles, const std::string& options = "");

    TTree* BuildEnsembles(TTree* tree, int nensembles, const std::string& options = "");

    TTree* PerformEnsembleTest(const std::vector<double>& parameters, int nensembles, const std::string& options = "");

    TTree* PerformEnsembleTest(TTree* tree, int nensembles, int start = 0, const std::string& options = "");

    std::vector<TH1D> MatrixToHistograms(const std::vector< std::vector<double> >& matrix);

private:

    BCMTF* fMTF;

    TRandom3* fRandom;

    bool fFlagMarginalize;

    BCLog::LogLevel fLogLevel;

};
// ---------------------------------------------------------

#endif