# Save canvas and root file with graph, const term and sampling term if calo_init.output(0): canv.SaveAs(calo_init.output(0)+".eps") canv.SaveAs(calo_init.output(0)+".png") plots = TFile(calo_init.output(0)+".root","RECREATE") if calo_init.args.preview: cPreview.SaveAs("preview_"+calo_init.output(0)+".png") else: canv.SaveAs("sampling_fraction_plots.eps") canv.SaveAs("sampling_fraction_plots.png") plots = TFile("sampling_fraction.root","RECREATE") if calo_init.args.preview: cPreview.SaveAs("preview_sampling_fraction.png") for g in all_graphs: g.Write() mean = numpy.zeros(1, dtype=float) std = numpy.zeros(1, dtype=float) t = TTree("samplingFraction", "Sampling fraction for detector layers") t.Branch("mean", mean, "mean/D"); t.Branch("std", std, "std/D"); for islice in range(0, Nslicesmerged): for ilay in range(0, calo_init.args.merge[islice]): mean[0] = gSF.GetY()[islice] std[0] = gSF.GetErrorY(islice) t.Fill() plots.Write() plots.Close() raw_input("Press ENTER to exit")