mlfit_result.C

double 
read_value(const char* path, const char* value, bool isBG) 
{ 
  TFile* file = new TFile(TString::Format("%s/out/mlfit.root", path));
  TTree* tree = (TTree*) file->FindObjectAny(isBG ? "tree_fit_b" : "tree_fit_sb");
  double val  = 0; tree->SetBranchAddress(value, &val); tree->GetEntry(0);				
  return val;
}

void 
mlfit_result(const char* path, const char* value="nll_min") 
{ 
  TFile* outFile = new TFile(TString::Format("%s/mlfit_result.root", path),"RECREATE");
  TTree *outTree = new TTree("limit","limit");
  char* output_value = (value =="nll_min") ? "limit" : value;
  double val_s = 0; outTree->Branch(output_value, &val_s, TString::Format("%s/D", output_value));
  double val_b = 0; outTree->Branch(TString::Format("%s_b", output_value), &val_b, TString::Format("%s_b/D", output_value));
  val_s = read_value(path, value, false);
  val_b = read_value(path, value, true );
  outFile->cd();
  outTree->Fill ();
  outTree->Write();
}