MultiDimFit.h

Go to the documentation of this file.

#ifndef HiggsAnalysis_CombinedLimit_MultiDimFit_h
#define HiggsAnalysis_CombinedLimit_MultiDimFit_h

#include "../interface/FitterAlgoBase.h"
#include <RooRealVar.h>
#include <vector>

class MultiDimFit : public FitterAlgoBase {
public:
  MultiDimFit() ;
  virtual const std::string & name() const {
    static const std::string name("MultiDimFit");
    return name;
  }
  virtual void applyOptions(const boost::program_options::variables_map &vm) ;

protected:
  virtual bool runSpecific(RooWorkspace *w, RooStats::ModelConfig *mc_s, RooStats::ModelConfig *mc_b, RooAbsData &data, double &limit, double &limitErr, const double *hint);

  enum Algo { None, Singles, Cross, Grid, RandomPoints, Contour2D };
  static Algo algo_;

  static std::vector<std::string>  poi_;
  static std::vector<RooRealVar*>  poiVars_;
  static std::vector<float>        poiVals_;
  static RooArgList                poiList_; 
  static unsigned int              nOtherFloatingPoi_; // keep a count of other POIs that we're ignoring, for proper chisquare normalization
  static float                     deltaNLL_;

  // options    
  static unsigned int points_, firstPoint_, lastPoint_;
  static bool floatOtherPOIs_;
  static bool fastScan_;

  // initialize variables
  void initOnce(RooWorkspace *w, RooStats::ModelConfig *mc_s) ;

  // variables
  void doSingles(RooFitResult &res) ;
  void doGrid(RooAbsReal &nll) ;
  void doRandomPoints(RooAbsReal &nll) ;
  void doContour2D(RooAbsReal &nll) ;

  // utilities
  void doBox(RooAbsReal &nll, double cl, const char *name="box", bool commitPoints=true) ;
};


#endif