mcnp_plot_out_spec: Convert histogram data to average points and plot as...

In radsafer: Radiation Safety

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")

radsafer documentation built on July 26, 2023, 6:03 p.m.