Chi2ForConicalWaveFit.h

Go to the documentation of this file.

#ifndef _Chi2ForConicalWaveFit_h_
#define _Chi2ForConicalWaveFit_h_

#include <iostream>
#include <fstream>
#include <vector>
#include <TVector3.h>
#include <TMath.h>
// #include <utl/BasicVector.h>
#include <utl/Vector.h>
#include <utl/Point.h>
#include <utl/TimeStamp.h>
// #include <utl/AugerUnits.h>                  // already included by PhysicalConstants
// #include <utl/MathConstants.h>               // already included by PhysicalConstants
#include <utl/PhysicalConstants.h>

#include <Math/IFunction.h>

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

        // Use cartesic coordinate system
        void UseCartesicCoordinates();
        // Use cylindrical coordinate system (default)
        void UseCylindricalCoordinates();

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

        // Mandatory functions of Clone, NDim ROOT::Math::IBaseFunctionMultiDim:
        virtual ROOT::Math::IBaseFunctionMultiDim* Clone() const {
                Chi2ForConicalWaveFit* foo = new Chi2ForConicalWaveFit();
                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 Cylindrical;
};
#endif /* CHI2FORCONICALWAVEFIT_H_ */