SdMonteCarloPropagator.h

Go to the documentation of this file.

#ifndef _SdMonteCarloPropagatorKVI_SdMonteCarloPropagator_h_
#define _SdMonteCarloPropagatorKVI_SdMonteCarloPropagator_h_

#include <fwk/VModule.h>
#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/TimeStamp.h>
#include <utl/CovarianceMatrix.h>
#include <evt/Event.h>
#include <fwk/VModule.h>

#include <boost/random.hpp>
#include <boost/random/normal_distribution.hpp>

#include <cmath>


namespace SdMonteCarloPropagatorKVI {

  class SdMonteCarloPropagator : public fwk::VModule {
  public:
    SdMonteCarloPropagator() : fNormalNoise(0, 1), fRandomGenerator(fRandgen, fNormalNoise) { }

    fwk::VModule::ResultFlag Init() override;
    fwk::VModule::ResultFlag Run(evt::Event& event) override;
    fwk::VModule::ResultFlag Finish() override { return eSuccess; }

    enum InfoLevel {
      eNone = 0,
      eFinal = 1,
      eIntermediate = 2,
      eObscure = 3
    };

  private:
    int fInfoLevel = eNone;
    int fRandomSeed = 0;
    // Not sure obout the random engine.
    // utl::RandomEngine* fRandomEngine;
    // using the boost Mersenne twister for now...
    // DV: you are not allowed to use private random number sources
    boost::mt19937 fRandgen;
    boost::normal_distribution<double> fNormalNoise;
    boost::variate_generator<boost::mt19937, boost::normal_distribution<double>> fRandomGenerator;

    REGISTER_MODULE("SdMonteCarloPropagatorKVI", SdMonteCarloPropagator);

  };

}


#endif