LikeFCN.h

Go to the documentation of this file.

#ifndef _MdLDFFinderAG_LikeFCN_h_
#define _MdLDFFinderAG_LikeFCN_h_

#include <Minuit2/FCNBase.h>
#include <vector>
#include <utility>


namespace MdLDFFinderAG {

  typedef std::vector<std::pair<unsigned int, unsigned int> > SegmentsOnMulti;


  class LikeFCN : public ROOT::Minuit2::FCNBase {

    static const double EPSILON; 
    // set the dependent parameter to this value when pushed below 0 by the independent parameters

  public:
    LikeFCN(const unsigned int, const std::vector<unsigned int>&, const double);
    // Constructor without multiplicity

    LikeFCN(const unsigned int, const SegmentsOnMulti&, const double);
    // Constructor using multiplicity (optimisation)   

    double Up() const { return 0.5; } // error definition corresponding to a likelihood 

    double operator()(const std::vector<double>&) const;

    double operator()(const double* const mup, const double* const) const;
    // required by ROOT to implement a TF1 with a functor - used only for drawing purposes
 
    double GetMu() const { return fMuons; }

    unsigned int GetN() const { return fNumberSegments; }

    SegmentsOnMulti GetK() const { return fSegmentsOnMulti; }
  
  private:
    unsigned int fNumberSegments;   // number of segments

    SegmentsOnMulti fSegmentsOnMulti;
    // holds the number of segments on and their multiplicity  

    double fMuons; // total number of muons 

  };

}


#endif