BCModelOutput.h

#ifndef __BCMODELOUTPUT__H
#define __BCMODELOUTPUT__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 <vector>

// BAT classes
class BCModel;

// ROOT classes
class TTree;
class TFile;
class TObject;

// ---------------------------------------------------------

class BCModelOutput
{
   public:

      BCModelOutput();

      BCModelOutput(BCModel * model, const char * filenname);

      BCModelOutput(const BCModelOutput & modeloutput);

      virtual ~BCModelOutput();

      BCModelOutput & operator = (const BCModelOutput & modeloutput);

      TTree * GetAnalysisTree()
         { return fAnalysisTree; };

      TFile * GetFile()
         { return fOutputFile; };

      void SetModel(BCModel * model);


      void SetFile(const char * filename);

      void WriteMarkovChain(bool flag = true);

      void FillAnalysisTree();

      void WriteMarginalizedDistributions();

      void Write(TObject * o);

      void Close();

   private:

      void Init();

      void Copy(BCModelOutput & modeloutput) const;

      void InitializeAnalysisTree();

      void InitializeSATree();

      TTree * fAnalysisTree;

      std::vector<TTree *> fMarkovChainTrees;

      TTree * fTreeSA;

      char * fFileName;

      TFile * fOutputFile;

      BCModel * fModel;

      int fIndex;
      unsigned int fNParameters;
      double fProbability_apriori;
      double fProbability_aposteriori;
      std::vector<double> fMode_global;
      std::vector<double> fMode_marginalized;
      std::vector<double> fMean_marginalized;
      std::vector<double> fMedian_marginalized;
      std::vector<double> fQuantile_05;
      std::vector<double> fQuantile_10;
      std::vector<double> fQuantile_16;
      std::vector<double> fQuantile_84;
      std::vector<double> fQuantile_90;
      std::vector<double> fQuantile_95;

};

// ---------------------------------------------------------

#endif