theory.plotOn(plot)
# theory2.plotOn(plot, LineColor(kRed), LineStyle(kDashed))
plot.Draw()
plots.append(plot)
## Photon Resolution
canvases.next('phoERes').SetLogy()
phoERes.SetTitle('E^{#gamma}_{reco}/E^{#gamma}_{gen} - 1')
phoERes.setRange(-0.9,3.2)
plot = phoERes.frame()
plot.SetTitle('Photon Resolution Modeling Zoom Out')
reducedData['phoERes'].plotOn(plot)
phoEResShape.plotOn(plot)
plot.Draw()
plot.GetYaxis().SetRangeUser(1e-3, 1e4)
llabels.position = (0.6, 0.83)
llabels.draw()
plots.append(plot)
## Photon Resolution Zoom In
canvases.next('phoEResZoom')
plot = phoERes.frame(Range(-0.25 + phoEScaleMC, 0.33 + phoEScaleMC))
plot.SetTitle('Photon Resolution Modeling Zoom In')
reducedData['phoERes'].plotOn(plot)
phoEResShape.plotOn(plot)
plot.Draw()
llabels.draw()
llabels.position = (0.21, 0.83)
plots.append(plot)
## ## Photon Resolution Zoom Out
