GenerateEvent.h

Go to the documentation of this file.

/*
 * Copyright (C) 2013 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef GENERATEEVENT_H
#define GENERATEEVENT_H
 
#include "InputParser.h"
#include "Observable.h"
#include "Observable2D.h"
#include "CorrelatedGaussianObservables.h"
#include "ModelParameter.h"
#include "Model.h"
#include "boost/shared_ptr.hpp"
#include "boost/make_shared.hpp"
#include <TF1.h>
#ifdef _MPI
#include <mpi.h>
#endif
 
class GenerateEvent {
public:
    GenerateEvent(ModelFactory& ModelF, ThObsFactory& ThObsF, 
                  const std::string& ModelConf_i,
                  const std::string& OutDirName_i,
                  const std::string& JobTag_i
                  /*const bool noMC_i*/);
    
    virtual ~GenerateEvent();
    
    void generate(int unsigned nIteration_i, int seed = 0, bool weight_i = false);
    
    void addCustomObservableType(const std::string name, boost::function<Observable*() > funct);
    
private:
    void defineParameterDistributions();
    
    void generateRandomEvent(int iterationNo);
    
    void createDirectories();
    
    void initModel();
    
    void setDParsFromParameters(const std::vector<double>& parameters, std::map<std::string,double>& DPars_i);
    
    InputParser myInputParser; 
    std::map<std::string, double> DPars; 
    std::map<std::string, TF1*> DDist; 
    std::map<std::string, boost::shared_ptr<std::ofstream> > ObsOut;
    std::map<std::string, boost::shared_ptr<std::ofstream> > CGOOut;
    std::map<std::string, boost::shared_ptr<std::ofstream> > ParsOut;
    StandardModel* Mod; 
    std::vector<ModelParameter> ModPars; 
    std::vector<ModelParameter> ModParsVar; 
    boost::ptr_vector<Observable> Obs; 
    std::vector<Observable2D> Obs2D; 
    std::vector<CorrelatedGaussianObservables> CGO; 
    std::vector<CorrelatedGaussianParameters> CGP; 
    std::string ModelConf; 
    std::string OutDirName; 
    std::string OldOutDirName; 
    std::string ObsDirName; 
    std::string CGODirName; 
    std::string SMDebugDirName; 
    std::string SUSYDebugDirName; 
    std::string ParsDirName; 
    std::string JobTag; 
    //bool noMC;///< Flag to initiate noMC mode.
    bool outputTerm; 
    int rank;
    int procnum;
#ifdef _MPI
    int nameLen;
    char processorName[MPI_MAX_PROCESSOR_NAME];
#endif
    std::ofstream summary;
    int nIteration;
    int positionID;
    int buffersize;
    double *sendbuff;
    double *sendbuff_w;
    double **buff;
    double **buff_w;
    int *sendbuff_int;
    int **buff_int;
    std::string ModelName;
    bool weight;
};
 
#endif  /* GENERATEEVENT_H */