MonteCarlo.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef MONTECARLO_H
#define MONTECARLO_H
 
#include "InputParser.h"
#include "MonteCarloEngine.h"
 
#define NBINSMODELPARS 100
 
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 RunMinuitOnly; 
    std::string CalculateNormalization; 
    int NIterationNormalizationMC; 
    bool PrintAllMarginalized; 
    bool PrintCorrelationMatrix; 
    bool PrintKnowledgeUpdatePlots; 
    bool PrintParameterPlot; 
    bool PrintTrianglePlot; 
    bool WritePreRunData; 
    bool checkrun; 
    bool writechains; 
    
    void ParseMCMCConfig(std::string file);
    void ReadPreRunData(std::string file);
};
 
#endif  /* MONTECARLO_H */