#!/usr/bin/python
import os,sys
from ROOT import *
from clean_canvas import *
gStyle = makegstyle(gStyle)

fn_data = '/home/jbueno/pass6/sum74-4-1_clk_raw.root'
fn_sim  = '/home/jbueno/pass6/sum474-1_2_21_31-1_clk_raw.root'

f_data = TFile(fn_data)
h_data = f_data.Get('MuLastPCut/MuLastP_before').Clone()

f_sim = TFile(fn_sim)
h_sim = f_sim.Get('MuLastPCut/MuLastP_before').Clone()

print 'n entries in data=', h_data.GetEntries()
#n_sim = h_sim.GetEntries()
n_data= h_data.GetMaximum()
n_sim = h_sim.GetMaximum()


c = TCanvas()
c.SetLogy()

h_data.SetMinimum(1E-1)
h_data.SetMaximum(1E9)

h_data.SetLineWidth(2)
h_data.SetLineColor(4)

h_sim.SetLineWidth(2)
h_sim.SetLineColor(2)
h_sim.SetLineStyle(7)

h_data.SetTitle(';muon last plane hit;muons')

h_data.Draw()
h_sim.Draw('same')
h_sim.Scale(n_data/n_sim)

l = TLegend(.8,.8,.9,.9)
l.SetBorderSize(0)
#l.AddEntry(h_sim,'simulation (renormalised)','l')
l.AddEntry(h_sim,'simulation','l')
l.AddEntry(h_data,'data','l')

l.Draw('same')

#c2 = TCanvas()
#diff = h_data.Clone()
#diff.Add(h_sim.Clone(),-1.)
#diff.Draw()

au = raw_input('>')