MonteCarlo.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

#ifndef MONTECARLO_H
#define MONTECARLO_H

#include "InputParser.h"
#include "MonteCarloEngine.h"
#include <BAT/BCH1D.h>
#include <BAT/BCH2D.h>

class MonteCarlo {
public:
    MonteCarlo(ModelFactory& ModelF, ThObsFactory& ThObsF,
            const std::string& ModelConf_i,
            const std::string& MonteCarloConf_i,
            const std::string& OutFile_i,
            const std::string& JobTag_i);

    //virtual ~MonteCarlo();

    void TestRun(int rank);

    void Run(const int rank);
    
    std::map<std::string, BCH1D * > getHistograms1D() const
    {
        return MCEngine.getHistograms1D();
    }

    std::map<std::string, BCH2D * > getHistograms2D() const
    {
        return MCEngine.getHistograms2D();
    }
    
   void addCustomObservableType(const std::string name, boost::function<Observable*() > funct);
   
private:
    std::string ModelName; 
    InputParser myInputParser; 
    MonteCarloEngine MCEngine; 
    std::vector<ModelParameter> ModPars; 
    boost::ptr_vector<Observable> Obs; 
    std::vector<Observable2D> Obs2D; 
    std::vector<CorrelatedGaussianObservables> CGO; 
    std::vector<CorrelatedGaussianParameters> CGP; 
    std::string ModelConf; 
    std::string MCMCConf; 
    std::string OutFile; 
    std::string JobTag; 
    std::string ObsDirName; 
    bool FindModeWithMinuit; 
    bool CalculateNormalization; 
    bool PrintAllMarginalized; 
    bool PrintCorrelationMatrix; 
    bool PrintKnowledgeUpdatePlots; 
    bool PrintParameterPlot; 
    bool WritePreRunData; 
    bool checkrun; 
    double normalization; 
    
    void ReadPreRunData (std::string file);
};

#endif  /* MONTECARLO_H */