/data/refman/pasoursint/CMSSW_5_3_6/src/HiggsAnalysis/CombinedLimit/interface/GoodnessOfFit.h

Go to the documentation of this file.

#ifndef HiggsAnalysis_CombinedLimit_GoodnessOfFit_h
#define HiggsAnalysis_CombinedLimit_GoodnessOfFit_h

#include "../interface/LimitAlgo.h"
#include "../interface/ProfileLikelihood.h"

class GoodnessOfFit : public LimitAlgo {
public:
  GoodnessOfFit() ;
  virtual bool run(RooWorkspace *w, RooStats::ModelConfig *mc_s, RooStats::ModelConfig *mc_b, RooAbsData &data, double &limit, double &limitErr, const double *hint);
  virtual const std::string & name() const {
    static const std::string name("GoodnessOfFit");
    return name;
  }
  virtual void applyOptions(const boost::program_options::variables_map &vm) ;

  virtual bool runSaturatedModel(RooWorkspace *w, RooStats::ModelConfig *mc_s, RooStats::ModelConfig *mc_b, RooAbsData &data, double &limit, double &limitErr, const double *hint);
 
protected:
  static std::string algo_;

  static std::string minimizerAlgo_;
  static float       minimizerTolerance_;
  static int         minimizerStrategy_;

  static float mu_;
  static bool  fixedMu_;

  // Return a pdf that matches this data perfectly.
  RooAbsPdf *makeSaturatedPdf(RooAbsData &data);
  mutable std::vector<RooAbsData*> tempData_;

};


#endif