Chi2ForSphericalWaveFitVarC.h

Go to the documentation of this file.

#ifndef _Chi2ForSphericalWaveFitVarC_h_
#define _Chi2ForSphericalWaveFitVarC_h_

#include <iostream>
#include <fstream>
#include <vector>
#include <TVector3.h>
#include <TMath.h>

#include <utl/Vector.h>
#include <utl/Point.h>
#include <utl/TimeStamp.h>
#include <utl/PhysicalConstants.h>

#include <Math/IFunction.h>


  
class Chi2ForSphericalWaveFitVarC : public ROOT::Math::IBaseFunctionMultiDim
{
public:
  Chi2ForSphericalWaveFitVarC();
  ~Chi2ForSphericalWaveFitVarC();

  // Set Antenna Times and Positions
  void Set(const std::vector<utl::Vector> &_AntennaPositions, const std::vector<double> &_AntennaTimes, const std::vector<double> &_AntennaTimesError, const utl::CoordinateSystemPtr &_fgLocalCS);

  void SetSigmaGamma(double _sigma_gamma);

  // using penalty term in chi squared computation
  void SetPenalty(bool dP);
  
  // Mandatory functions of Clone, NDim ROOT::Math::IBaseFunctionMultiDim:
  virtual ROOT::Math::IBaseFunctionMultiDim* Clone() const {Chi2ForSphericalWaveFitVarC* foo = new Chi2ForSphericalWaveFitVarC();
  foo->Set(AntennaPositions,AntennaTimes,AntennaTimesError, fgLocalCS);
  return foo;};
  unsigned int NDim() const {int n=4;return n;};

  // Objective function
  double DoEval(const double* x) const;

private:
  std::vector<utl::Vector> AntennaPositions;
  std::vector<double> AntennaTimes;
  std::vector<double> AntennaTimesError;
  utl::CoordinateSystemPtr fgLocalCS;

  bool Spherical;
  bool doPenalty;
  bool computeSigmaGamma;
  double sigma_gamma;
};

#endif