LauMinuit.hh

Go to the documentation of this file.

// Copyright University of Warwick 2013 - 2013.
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

// Authors:
// Thomas Latham
// John Back
// Paul Harrison

#ifndef LAU_MINUIT
#define LAU_MINUIT

#include <utility>
#include <vector>

#include "Rtypes.h"
#include "TMatrixD.h"

#include "LauAbsFitter.hh"

class LauParameter;
class TVirtualFitter;


class LauMinuit : public LauAbsFitter {

        public:
                virtual ~LauMinuit();


                virtual void initialise( LauFitObject* fitObj, const std::vector<LauParameter*>& parameters );

                virtual LauFitObject* getFitObject();

                virtual UInt_t nParameters() const { return nParams_; }

                virtual UInt_t nFreeParameters() const { return nFreeParams_; }

                virtual Bool_t twoStageFit() const {return twoStageFit_;}


                virtual void twoStageFit(Bool_t doTwoStageFit) {twoStageFit_ = doTwoStageFit;}

                virtual Bool_t useAsymmFitErrors() const {return useAsymmFitErrors_;}


                virtual void useAsymmFitErrors(Bool_t useAsymmErrors) {useAsymmFitErrors_ = useAsymmErrors;}


                virtual std::pair<Int_t,Double_t> minimise();

                virtual void fixSecondStageParameters();

                virtual void releaseSecondStageParameters();

                virtual void updateParameters();

                virtual const TMatrixD& covarianceMatrix() const { return covMatrix_; }


        private:
                friend class LauFitter;

                explicit LauMinuit( Int_t maxPar = 100 );

                LauMinuit( const LauMinuit& );

                LauMinuit& operator=( const LauMinuit& rhs );

                TVirtualFitter* minuit_;

                const UInt_t maxPar_;

                std::vector<LauParameter*> params_;

                UInt_t nParams_;

                UInt_t nFreeParams_;

                Bool_t twoStageFit_; 

                Bool_t useAsymmFitErrors_; 

                Int_t fitStatus_;

                Double_t NLL_;

                TMatrixD covMatrix_;

                ClassDef(LauMinuit,0);
};

#endif