def create_significance_files(filename, signal):
print 'starting to create significance file from', filename
work_file = TFile(filename, 'READ')
keys = work_file.GetListOfKeys()
categories_done = []
hist_sig = []
hist_other = []
masses = []
for key in keys:
key = str(key.GetName())
if signal in key: hist_sig.append(work_file.Get(key).Clone())
else: hist_other.append(work_file.Get(key).Clone())
for h in hist_sig:
name = h.GetName()
name = name.replace('10p', '')
name = name.replace('20p', '')
name = name.replace('30p', '')
name = ''.join(x for x in name if x.isdigit())
masses.append(name)
output = TFile(
filename.replace('.root', name + '_' + signal + '.root'),
'RECREATE')
output.Cd()
h.Write()
for o in hist_other:
o.Write()
output.Close()
work_file.Close()
return masses
def makeHist(name, title, Ytitle, file, vect, dir, nbins, xmin, xmax):
histFile = TFile(fileName, "UPDATE")
htest = TH1F(name, title, nbins, xmin, xmax)
for j in range(0, nbins):
htest.SetBinContent(j + 1, vect[j])
htest.SetXTitle("Luminosity Block")
htest.SetYTitle(Ytitle)
histFile.Cd(dir)
htest.Write()
histFile.Close()
