BCMTFAnalysisFacility.h

Go to the documentation of this file.

#ifndef __BCMTFANALYSISFACILITY__H
#define __BCMTFANALYSISFACILITY__H

/*
 * Copyright (C) 2008-2012, Daniel Kollar and Kevin Kroeninger.
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

// ---------------------------------------------------------
#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 SetFlagMCMC(bool flag)
         { fFlagMCMC = 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);

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

      TTree * BuildEnsembles(TTree * tree, int nensembles);

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

      TTree * PerformEnsembleTest(TTree * tree, int nensembles, int start = 0);

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

 private:

      BCMTF * fMTF;

      TRandom3 * fRandom;

      bool fFlagMCMC;

      BCLog::LogLevel fLogLevel;

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

#endif