#include "RooAbsPdf.h"
#include "RooWorkspace.h"
#include "TFile.h"
#include <RooStats/ModelConfig.h>
#include "RooRealVar.h"

RooAbsReal* nll = 0;
RooWorkspace* w = 0;
TFile* tf = 0;
RooAbsPdf* pdf = 0;
float doit(float input){
  if (!nll){
    tf = TFile::Open("higgsCombinenominal.MultiDimFit.mH120.root");
    w = (RooWorkspace*) tf->Get("w");
    w->loadSnapshot("MultiDimFit");
    RooAbsData &data = *w->data("data_obs");
    RooStats::ModelConfig* mc_s = (RooStats::ModelConfig*) w->obj("ModelConfig");
    pdf = mc_s->GetPdf();
    //const RooCmdArg &constrainCmdArg = RooFit::Constrain(*mc_s->GetNuisanceParameters());
    //nll = pdf.createNLL(data, constrainCmdArg, RooFit::Extended(pdf.canBeExtended()), RooFit::Offset(true));
  }
  auto var = w->var("r_ttH_2l_2tau"); 
  var->setVal(input);
  return pdf->getVal();
  

}