BCMTFAnalysisFacility.h

#ifndef __BCMTFANALYSISFACILITY__H
#define __BCMTFANALYSISFACILITY__H

/*
 * Copyright (C) 2007-2014, 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; };

        int PerformSingleChannelAnalyses(const char * dirname, const char * options = "");

        int PerformSingleSystematicAnalyses(const char * dirname, const char * options = "");

        int PerformCalibrationAnalysis(const char * 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, std::string options="");

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

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

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

        TTree * PerformEnsembleTest(TTree * tree, int nensembles, int start = 0, 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