ModelConfig* b_model = new ModelConfig("B_model");  
  b_model->SetWorkspace(*wspace);
  b_model->SetPdf(*wspace->pdf("model"));
  b_model->SetObservables(*wspace->set("obs"));
  b_model->SetParametersOfInterest(*wspace->set("parOfInterest")); 
  wspace->var("BR")->setVal(0.0);  
  b_model->SetSnapshot(*wspace->set("parOfInterest"));
  b_model->SetNuisanceParameters(*wspace->set("nuisPar"));

  ModelConfig* sb_model = new ModelConfig("SB_model");  
  sb_model->SetWorkspace(*wspace);
  sb_model->SetPdf(*wspace->pdf("model"));
  sb_model->SetObservables(*wspace->set("obs"));
  sb_model->SetParametersOfInterest(*wspace->set("parOfInterest")); 
  wspace->var("BR")->setVal(0.0032); 
  sb_model->SetSnapshot(*wspace->set("parOfInterest"));
  sb_model->SetNuisanceParameters(*wspace->set("nuisPar"));


  //Right order: b_model, sb_model
  HybridCalculator* hypoCalcFlip = new HybridCalculator(*data, *b_model, *sb_model);

  hypoCalcFlip->SetToys(8000, 8000);
  hypoCalcFlip->ForcePriorNuisanceAlt(*wspace->pdf("nuisPDF"));
  hypoCalcFlip->ForcePriorNuisanceNull(*wspace->pdf("nuisPDF"));

  ToyMCSampler* samplerFlip = (ToyMCSampler*)hypoCalcFlip->GetTestStatSampler();
  samplerFlip->SetNEventsPerToy(1);

  HypoTestInverter* inverter = new HypoTestInverter(*hypoCalcFlip, wspace->var("BR"));
  inverter->SetConfidenceLevel(0.95);
  inverter->RunOnePoint(PointValue);