MulticanonicalMCMCHandler.h

Go to the documentation of this file.

#pragma once
#include <ostream>
#include <string>
#include <vector>
 
#include "Parameters/ParameterHandlerBase.h"
#include "Manager/Manager.h"
 
_MaCh3_Safe_Include_Start_ //{
#include "TSpline.h"
_MaCh3_Safe_Include_End_ //}
 
 
namespace M3 {
  enum class BiasFunction {
    kGaussian,
    kVonMises,
    kGeneralisedGaussian
  };
 
  constexpr double UmbrellaGaussianNormFactor = 0.906402477055;
}
 
class MulticanonicalMCMCHandler {
public:
  MulticanonicalMCMCHandler();
  virtual ~MulticanonicalMCMCHandler();
 
  #ifdef MACH3_DEBUG
  void setDebugStream(std::ostream* os, bool enabled);
  #endif
 
  void FindOscCovParams(const std::vector<ParameterHandlerBase*>& systematics);
 
  void AdjustUmbrellaStepScale(const std::vector<ParameterHandlerBase*>& systematics);
 
  void InitializeMulticanonicalHandlerConfig(Manager* fitMan, std::vector<ParameterHandlerBase*>& systematics);
 
  void InitializeMulticanonicalParams(std::vector<ParameterHandlerBase*>& systematics);
 
  double GetMulticanonicalWeight(double deltacp, double delm23_value);
 
  double GetMulticanonicalWeightSpline(double deltacp, double delm23_value);
 
  double GetMulticanonicalWeightGaussian(double deltacp);
 
  double GetMulticanonicalWeightTripleGaussian(double deltacp);
 
  double GetMulticanonicalWeightVonMises(double deltacp);
 
  double GetMulticanonicalWeightGenGaussian(double deltacp);
 
  double circularDistance(double x, double mean);
 
  double generalisedGaussian2(double x, double mean, double width);
 
  int oscCovVar;
  int multicanonicalVar;
  int multicanonicalVar_dm23;
 
  M3::BiasFunction umbrellaBiasFunction;
 
  std::string umbrellaBiasFunctionName;
 
  bool multicanonicalSpline;
 
protected:
  double multicanonicalBeta;
 
  double delta_cp_value;
  double delm23_value;
 
  TSpline3* dcp_spline_IO;
  TSpline3* dcp_spline_NO;
 
  double umbrellaMean;
  double umbrellaWidth;
  int umbrellaNumber;
  bool umbrellaOverlapMode;
  double umbrellaSigmaOverlap;
  bool umbrellaAdjustStepScale;
  double umbrellaStepScaleFactor;
  bool flipWindow;
 
  double vonMises_kappa;
  double vonMises_I0_kappa;
 
#ifdef MACH3_DEBUG
  std::ostream* debugStream = nullptr;
  bool debugEnabled = false;
#endif
};