mcnp_plot_out_spec {radsafer} | R Documentation |
Convert histogram data to average points and plot as spectrum.
Description
Model results or input source histograms from MCNP and perhaps other sources typically provide binned tally results with columns representing maximum energy in MeV, a column with the mean tally result or bin probability and an uncertainty column (not used). Once the data is scanned in, or otherwise entered into the R global environment, they can be plotted with this function.
Usage
mcnp_plot_out_spec(spec.df, title = deparse(substitute(spec.df)), log_plot = 0)
Arguments
spec.df |
A data frame with no header. Maximum energy in MeV should be in the first column, (named E_MeV), and binned results in the second column, (named prob). |
title |
Title for chart (default = name of spec.df) |
log_plot |
0 = no log axes (default), 1 = log y-axis, 2 = log both axes. |
See Also
mcnp_scan_save()
to copy and paste output spectrum.
Other mcnp tools:
mcnp_cone_angle()
,
mcnp_est_nps()
,
mcnp_matrix_rotations()
,
mcnp_mesh_bins()
,
mcnp_scan2plot()
,
mcnp_scan_save()
,
mcnp_sdef_erg_hist()
,
mcnp_sdef_erg_line()
,
mcnp_si_sp_RD()
,
mcnp_si_sp_hist_scan()
,
mcnp_si_sp_hist()
Examples
mcnp_plot_out_spec(photons_cs137_hist, "example Cs-137 well irradiator")