HyperBinningMakerPhaseBinning.h

Go to the documentation of this file.

#ifndef HYPERBINNINGMAKERPHASEBINNING_HH
#define HYPERBINNINGMAKERPHASEBINNING_HH

// HyperPlot includes
#include "Mint/MessageService.h"
#include "Mint/HyperBinningMaker.h"
#include "Mint/LoadingBar.h"
#include "Mint/CyclicPhaseBins.h"

// Root includes
#include "TMath.h"

// std includes
#include <complex>




class HyperBinningMakerPhaseBinning : public HyperBinningMaker{
  
  private:
  
  int   _numBinPairs;
  CyclicPhaseBins _binEdges;

  int    _maximumRandWalks;
  int    _numWalkers;
  double     _walkSizeFrac;
  
  int _numberOfSystematicSplits;
  int _numberOfGradientSplits;
  
  private:

  int splitByCoord(int volumeNumber, int dimension, HyperPoint& coord);
  //Functions for determining the bin boundaries / bin number from a given bin number / phase.

  int getBinNumFromFunc(HyperPoint& point);
  int getBinNumFromFuncVal(double phase);
  double getLowBinBoundary(double phase);
  double getHighBinBoundary(double phase);
  double closestBinBoundary(double val);
  //Functions to evaluate derivatives and second derivatives of the phase motion

  HyperPoint getGrad(HyperPoint& point);
  HyperPoint getGradPos(HyperPoint& point, double funcValAtPoint);
  double getSecondDerivative(HyperPoint& point, HyperPoint& vector, double funcValAtPoint, double& deriv);
  int splitDimFromGrad(int volumeNumber, HyperPoint gradient);
  virtual int gradientSplit(int binNumber   , int& dimension);
  int systematicSplit      (int volumeNumber, int dimension, double valAtCenter, HyperPoint gradient);
  HyperPointSet getSplitCorners( int volumeNumber );
  HyperPointSet getSplitFaces  ( int volumeNumber );
  HyperPointSet getSplitEdges  ( int volumeNumber );
  HyperPoint orderAndTestSplitPoints(HyperPointSet& points, HyperPoint& point, double valAtPoint, HyperPoint gradient);
  int randomWalkSplit(int volumeNumber, int dimension);
  void walkOrthogonal(HyperPoint& point, HyperCuboid& walkLimits);
  void walk(HyperPoint& point, HyperCuboid& walkLimits);
  public:
  
  HyperBinningMakerPhaseBinning(const HyperCuboid& binningRange, HyperFunction* func);
  virtual void makeBinning();
  void setNumBinPairs(int binpairs); 
  void setBinEdges(std::vector<double> binEdges);
  virtual int gradientSplitAll();


  ~HyperBinningMakerPhaseBinning();
};

#endif