Commit 6d89f2ee authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

continue cleaning

parent 730c7da7
#ifndef __MCMC_sampler_class_guard__
#define __MCMC_sampler_class_guard__
#include <math.h>
#include <stdio.h>
#include <vector>
#include <string>
#include "AbsRand.hh"
#include "AbsFlat.hh"
#include "AbsGauss.hh"
#include "Drand48.hh"
#include "planck_rng.hh"
#include "GaussW.hh"
#include "pdf.hh"
#include "Parser.hh"
#include "CLHEP/Matrix/Vector.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "Parser.hh"
#include <cmath>
#include <cstdio>
#include <vector>
#include <string>
class MCMC_sampler{
protected:
bool _verbose;
public:
int start_step;
//protected:
int _dim; //number of parameters
int _length; //number of steps performed by the algorithm
int _burn_in; //to cut the initial phase of the chain
......@@ -46,22 +36,18 @@ class MCMC_sampler{
std::vector<int> _recur_chain;
std::vector<double> _lkh_values; //no set/get methods for the moment
//no set/get methods for the moment
//min and max values for the parameters
HepVector _min_par;
HepVector _max_par;
//maybe helpfull pars. Valid only after calling set_par(const char* str);
std::vector<std::string> _par_name;
// std::vector<double> _par_err;
HepVector _par_err;
HepSymMatrix _extracted_cov;
//std::vector<HepVector> _GR_param; //this is for R(l)
HepVector _GR_param;
std::vector<HepVector > _autocorr;
HepVector _autocor_time;
//constructors
MCMC_sampler();
......
......@@ -46,7 +46,7 @@ int main(int argc, char *argv[])
const string target_cor_path=parser.params.find<string>("target_cov","");
const string proposal_cor=parser.params.find<string>("proposal_cor","");
const string proposal_cov=parser.params.find<string>("proposal_cov","");
const int t0=parser.params.find<int>("t0",2000);
const int t0=parser.params.find<int>("t0",1000);
const int ts=parser.params.find<int>("ts",10000);
double scale=parser.params.find<double>("scale",0.001);
const bool do_move=parser.params.find<bool>("do_move",false);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment