# to run type: python ratios.py
import ROOT

File1 = ROOT.TFile("hist_S1.root","READ")
myhisto = File1.Higgs_smeared
myhisto.SetDirectory(0)
File1.Close()


File2 = ROOT.TFile("hist_S2.root","READ")
myhisto2 = File2.Higgs_smeared
myhisto2.SetDirectory(0)
#This puts the ratio into myhisto:
myhisto.Divide(myhisto2)

#This adds a fit to a constant. 
myhisto.Fit("pol0")
#

myhisto.GetYaxis().SetTitle("Ratio")
myhisto.Draw()
File2.Close()


OutputFile = ROOT.TFile("hist_RatioS1toS2.root","RECREATE")
OutputFile.cd()
myhisto.Write("Ratio_S1S2")
OutputFile.Close()