Chi2ForSphericalWaveFit.h

Go to the documentation of this file.

#ifndef _Chi2ForSphericalWaveFit_h_
#define _Chi2ForSphericalWaveFit_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 Chi2ForSphericalWaveFit : public ROOT::Math::IBaseFunctionMultiDim
{
public:
  Chi2ForSphericalWaveFit();
  ~Chi2ForSphericalWaveFit();

  // Use cartesic coordinate system
  void UseCartesicCoordinates();
  // Use spherical coordinate system (default)
  void UseSphericalCoordinates();

  // 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);
  

  // Mandatory functions of Clone, NDim ROOT::Math::IBaseFunctionMultiDim:
  virtual ROOT::Math::IBaseFunctionMultiDim* Clone() const {Chi2ForSphericalWaveFit* foo = new Chi2ForSphericalWaveFit();
  foo->Set(AntennaPositions,AntennaTimes,AntennaTimesError, fgLocalCS);
  return foo;};
  unsigned int NDim() const {int n=3;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;
};

#endif