def test_standard_plots(xsec_file, outdir='./'):
from pisa.utils.plotter import Plotter
xsec = genie.get_combined_xsec(xsec_file)
e_bins = MultiDimBinning(
[OneDimBinning(name='true_energy', tex=r'E_\nu', num_bins=150,
domain=(1E-1, 1E3)*ureg.GeV, is_log=True)]
)
xsec.compute_maps(e_bins)
logging.info('Making plots for genie xsec_maps')
plot_obj = Plotter(outdir=outdir, stamp='Cross-Section', fmt='png',
log=True, size=(12, 8),
label=r'Cross-Section ($m^{2}$)')
maps = xsec.return_mapset()
plot_obj.plot_xsec(maps, ylim=(1E-43, 1E-37))
def test_per_e_plot(xsec_file, outdir='./'):
from pisa.utils.plotter import Plotter
xsec = genie.get_combined_xsec(xsec_file)
e_bins = MultiDimBinning(
[OneDimBinning(name='true_energy', tex=r'E_\nu', num_bins=200,
domain=(1E-1, 1E3)*ureg.GeV, is_log=True)]
)
xsec.compute_maps(e_bins)
xsec.scale_maps(1/e_bins.true_energy.bin_widths.magnitude)
logging.info('Making plots for genie xsec_maps per energy')
plot_obj = Plotter(outdir=outdir, stamp='Cross-Section / Energy',
fmt='png', log=False, size=(12, 8),
label=r'Cross-Section / Energy ($m^{2} GeV^{-1}$)')
maps = xsec.return_mapset()
plot_obj.plot_xsec(maps, ylim=(3.5E-41, 3E-40))
