Movatterモバイル変換

radsafer: Radiation Safety

Description

Provides functions for radiation safety, also known as "radiation protection" and "radiological control". The science of radiation protection is called "health physics" and its engineering functions are called "radiological engineering". Functions in this package cover many of the computations needed by radiation safety professionals. Examples include: obtaining updated calibration and source check values for radiation monitors to account for radioactive decay in a reference source, simulating instrument readings to better understand measurement uncertainty, correcting instrument readings for geometry and ambient atmospheric conditions. Many of these functions are described in Johnson and Kirby (2011, ISBN-13: 978-1609134198). Utilities are also included for developing inputs and processing outputs with radiation transport codes, such as MCNP, a general-purpose Monte Carlo N-Particle code that can be used for neutron, photon, electron, or coupled neutron/photon/electron transport (Werner et. al. (2018)doi:10.2172/1419730).

Author(s)

Maintainer: Mark Hoguemark.hogue.chp@gmail.com

See Also

Useful links:

• https://github.com/markhogue/radsafer

• Report bugs athttps://github.com/markhogue/radsafer/issues

Pipe operator

Description

See⁠magrittr::[\%>\%][magrittr::\%>\%]⁠ for details.

Usage

lhs %>% rhs

Specific Activity

Description

Provides specific activity of a radionuclide in Bq/g.

Usage

RN_Spec_Act(RN_select, numeric = "n")

Arguments

Value

specific activity in Bq / g

See Also

Other radionuclides:RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

RN_Spec_Act("Ac-230")RN_Spec_Act(c("At-219", "Es-251"))RN_Spec_Act("Pd-96", numeric = "y")RN_Spec_Act(c("Cs-137", "Ba-137m"), numeric = "y")

Search for radioisotopes that dominate a specified energy bin

Description

Identify photon emitters that represent a target range of energies, while screeningout other selected energy ranges. This may be helpful for identifying radionuclidesin low-definition spectroscopy or in selecting representative spectra for modelingshielding.

Usage

RN_bin_screen_phot(  E_min = 0,  E_max = 10,  min_prob = 0,  min_half_life_seconds = NULL,  max_half_life_seconds = NULL,  no_E_min = 0,  no_E_max = 10,  no_min_prob = 100,  no_E_min2 = 0,  no_E_max2 = 10,  no_min_prob2 = 100)

Arguments

Value

radionuclides that match selection criteria

See Also

RN_plot_spectrum()

Other radionuclides:RN_Spec_Act(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

RN_bin_screen_phot(  E_min = 0.1, E_max = 0.3,  min_prob = 0.4, min_half_life_seconds = 30 * 24 * 3600,  max_half_life_seconds = 3.153e7, no_E_min = 0.015,  no_E_max = 0.0999, no_min_prob = 0.05, no_E_min2 = 0.301, no_E_max2 = 10, no_min_prob2 = 0.01)

Find a potential precursor of a radionuclide@description Find a potential parent radionuclide by searching the progeny fields in RadData ICRP_07.NDX

Description

Find a potential precursor of a radionuclide@description Find a potential parent radionuclide by searching the progeny fields in RadData ICRP_07.NDX

Usage

RN_find_parent(RN_select)

Arguments

Value

a subset of the data frame RadData::ICRP_07.NDX

Examples

Th_230_df <- RN_find_parent("Th-230")Tl_208_df <- RN_find_parent("Tl-208")

Screen radionuclide data to find matches to decay mode, half-life, and total emission energy

Description

Provides a set of radionuclides matching screening criteria. This is a limited screening based on average energy per transformation. ConsiderRN_search_phot_by_E,RN_search_alpha_by_E, andRN_search_beta_by_E for spectroscopic measurement matching.

Usage

RN_index_screen(  dk_mode = NULL,  min_half_life_seconds = NULL,  max_half_life_seconds = NULL,  min_E_alpha = NULL,  min_E_electron = NULL,  min_E_photon = NULL)

Arguments

Value

data frame of radionuclide data from the RadData package index data (RadData::ICRP_07.NDX), matching search criteria.

See Also

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

RN_index_screen(dk_mode = "SF")RN_index_screen(dk_mode = "IT", max_half_life_seconds = 433 * 3.15e7)

Quick table of Radionuclide Data from the RadData package

Description

Access a quick summary of radionuclide data. This is for convenience only and does not replace a more comprehensive view as is available in the Radiological Toolboxdoi:10.2172/1201298

Usage

RN_info(RN_select)

Arguments

Value

a table including half-life, decay modes, decay progeny, and branch fractions

See Also

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

Es_254m <- RN_info("Es-254m") #saves output to global environmentRN_info("Cf-252")RN_info("Cs-137")RN_info("Am-241")

Plot results of RN_search functions

Description

Plots results by radionuclide with E_MeV on x-axis and prob on y-axis.

Usage

RN_plot_df(discrete_df, title = deparse(substitute(discrete_df)), log_plot = 0)

Arguments

See Also

UseRN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),orRN_search_phot_by_E() and save the results,e.g. save_results <- RN_search_phot_by_E(0.99, 1.01, 13 * 60, 15 * 60, 1e-4)

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

search_results <- RN_search_phot_by_E(0.99, 1.01, 13 * 60, 15 * 60, 1e-4)RN_plot_df(search_results, title = "example1", log_plot = 0)search_results <- RN_save_spectrum("In-115m", photon = TRUE)RN_plot_df(search_results, title = "In-115m", log_plot = 0)

Plot results of RN_search functions

Description

Plots results by radionuclide with E_MeV on x-axis and prob on y-axis.:This function is deprecatedand will be removed in a future package revision.For now, it is still usable.The replacement,RN_plot_df plots saved data frames including those savedwith search functions.

Usage

RN_plot_search_results(  discrete_df,  title = deparse(substitute(discrete_df)),  log_plot = 0)

Arguments

See Also

UseRN_search_alpha_by_E(),RN_search_beta_by_E(),orRN_search_phot_by_E() and save the results,e.g. save_results <- RN_search_phot_by_E(0.99, 1.01, 13 * 60, 15 * 60, 1e-4)

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

search_results <- RN_search_phot_by_E(0.99, 1.01, 13 * 60, 15 * 60, 1e-4)RN_plot_search_results(search_results, title = "example1", log_plot = 0)

Plot radionuclide emission spectra directly from ICRP-107 data.

Description

Plot emission spectra based on radionuclide and desired radiation type.Plot on log axes if desired.Select cutoff value for probability optional, included at 1% by default.Plot includes energy times probability for dosimetric importance comparisons.

Usage

RN_plot_spectrum(  desired_RN,  rad_type = NULL,  photon = FALSE,  log_plot = 0,  prob_cut = 0.01)

Arguments

Value

plot of spectrum

See Also

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

RN_plot_spectrum(  desired_RN = c("Sr-90", "Y-90"), rad_type = "B-",  photon = FALSE, prob_cut = 0.01)RN_plot_spectrum(  desired_RN = c("Co-60", "Ba-137m"), rad_type = NULL,  photon = TRUE, prob_cut = 0.015)RN_plot_spectrum(  desired_RN = c("Co-60", "Ba-137m"), rad_type = NULL,  photon = TRUE, log_plot = 0)RN_plot_spectrum(desired_RN = c("Co-60", "Ba-137m"), rad_type = "G")RN_plot_spectrum(  desired_RN = c("Pu-238", "Pu-239", "Am-241"), rad_type = "A",  photon = FALSE, prob_cut = 0.01, log_plot = 0)

Save radionuclide emission spectra.

Description

Save emission spectra based on radionuclide and desired radiation type.Select cutoff value for probability optional, included at 1% by default.

Usage

RN_save_spectrum(desired_RN, rad_type = NULL, photon = FALSE, prob_cut = 0)

Arguments

Value

Dataframe with energy spectra - including probability of emission quantum, or, forbeta, the probability density.

See Also

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

Sr_Y_90_df <- RN_save_spectrum(desired_RN = c("Sr-90", "Y-90"), rad_type = "B-", photon = FALSE, prob_cut = 0.01)Co_60_Ba_137m_p_df <- RN_save_spectrum(desired_RN = c("Co-60", "Ba-137m"), rad_type = NULL, photon = TRUE, prob_cut = 0.015)Co_60_Ba_137m_g_df <- RN_save_spectrum(desired_RN = c("Co-60", "Ba-137m"), rad_type = "G")actinide_a_df <- RN_save_spectrum(desired_RN = c("Pu-238", "Pu-239", "Am-241"), rad_type = "A",photon = FALSE, prob_cut = 0.01)

Search for alpha

Description

Search for alpha emission based on energy, half-life and minimum probability.

Usage

RN_search_alpha_by_E(  E_min = 0,  E_max = 10,  min_half_life_seconds = NULL,  max_half_life_seconds = NULL,  min_prob = 0)

Arguments

Value

search results in order of half-life. Recommend assigningresults to a viewable object, such as 'search_results'

See Also

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_beta_by_E(),RN_search_phot_by_E()

Examples

# between 7 and 8 MeVsearch_results <- RN_search_alpha_by_E(7, 8)# 1-4 MeV; half-life between 1 and 4 hourssearch_results <- RN_search_alpha_by_E(1, 4, 1 * 3600, 4 * 3600)# between 7 and 10 MeV with at least 1e-3 probabilitysearch_results <- RN_search_alpha_by_E(7, 10, min_prob = 1e-3)

Search for beta

Description

Search for beta emission based on maximum energy and half-life.

Usage

RN_search_beta_by_E(  E_max,  min_half_life_seconds = NULL,  max_half_life_seconds = NULL)

Arguments

Value

search results in order of half-life. Recommend assigningresults to a viewable object, such as 'search_results'

See Also

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_phot_by_E()

Examples

# Max beta at least 2 MeVsearch_results <- RN_search_beta_by_E(2)# Max beta at least 2 MeV and half-life between 1 s and 1 hsearch_results <- RN_search_beta_by_E(2, 1, 3600)# Max beta at least 1 MeV and half-life between 1 d and 2 dsearch_results <- RN_search_beta_by_E(1, 3600 * 24, 2 * 3600 * 24)

Search for photon

Description

Search for photon emission based on energy, half-life and minimum probability.

Usage

RN_search_phot_by_E(  E_min = 0,  E_max = 10,  min_half_life_seconds = NULL,  max_half_life_seconds = NULL,  min_prob = 0)

Arguments

Value

search results in order of half-life. Recommend assigningresults to a viewable object, such as 'search_results'

See Also

RN_plot_spectrum()

Other radionuclides:RN_Spec_Act(),RN_bin_screen_phot(),RN_index_screen(),RN_info(),RN_plot_df(),RN_plot_search_results(),RN_plot_spectrum(),RN_save_spectrum(),RN_search_alpha_by_E(),RN_search_beta_by_E()

Examples

# between 1 and 1.2 MeV, between 6 and 6.2 hours half-life,# ... probability at least 1e-4search_results <- RN_search_phot_by_E(1, 1.2, 6 * 3600, 6.2 * 3600, 1e-4)# between 0.1 and 0.15 MeV, between 1 and 3 million years half-lifesearch_results <- RN_search_phot_by_E(0.1, 0.15, 1e6 * 3.153e7, 3e6 * 3.153e7)

Correct for air density - useful for vented ion chambers

Description

Obtain a correction factor for ion chamber temperature andpressure vs reference calibration values.

Usage

air_dens_cf(T.actual, P.actual, T.ref = 20, P.ref = 760)

Arguments

Value

The ratio of actual to reference air density.

See Also

Other rad measurements:disk_to_disk_solid_angle(),neutron_geom_cf(),scaler_sim(),tau_estimate()

Examples

air_dens_cf(T.actual = 20, P.actual = 760, T.ref = 20, P.ref = 760)air_dens_cf(30, 750)

Calculate fractional solid angle for disk to disk

Description

Returns fractional solid angle for a geometry frequentlyencountered in health physics analysis of air samples or disk smears. Thisis useful in correcting configurations that do not exactly matchcalibration (by ratioing the respective fractional solid angles). Whileunits of steridian are used for solid angle, this function only uses afraction of the total field of view.

Usage

disk_to_disk_solid_angle(  r.source,  gap,  r.detector,  plot.opt = "n",  runs = 10000,  off_center = 0,  beep = "off")

Arguments

Value

Fractional solid angle and plot of simulation.

References

https://karthikkaranth.me/blog/generating-random-points-in-a-sphere/https://en.wikipedia.org/wiki/Algorithms_for_calculating_variance

See Also

Other rad measurements:air_dens_cf(),neutron_geom_cf(),scaler_sim(),tau_estimate()

Examples

disk_to_disk_solid_angle(r.source = 15, gap = 20, r.detector = 10, plot.opt = "n", runs = 1e3)

Correct activity-dependent value based on radioactive decay.

Description

Decay-corrected values are provided. Either a single or multiplevalues are computed. The computation is made either based on a singleradionuclide, or based on user-provided half-life, with time unit. Thedifferential time is either computed based on dates entered or time lapsedbased on the time unit. Time units must be consistent. Decay-correct asource to today's date by assigning a referencedate1 and allowingdefaultdate2, the system date.

Usage

dk_correct(  RN_select = NULL,  half_life = NULL,  time_unit = NULL,  time_lapse = NULL,  date1 = NULL,  date2 = Sys.Date(),  A1 = 1,  num = FALSE)

Arguments

Value

Decay adjusted activity or related parameter. SeeA1.

See Also

Other decay corrections:dk_pct_to_num_half_life(),dk_time()

Examples

# RN_select and date1 (saving numerical data)my_dks <- dk_correct(  RN_select = "Sr-90",  date1 = "2009-01-01",  date2 = "2019-01-01",  num = TRUE)#   RN_select and time_lapse (random sample)dk_correct(  RN_select = base::sample(RadData::ICRP_07.NDX$RN, 1),  time_lapse = 1:10,  time_unit = base::sample(c("y", "d", "h", "m", "s"), 1))#   half_life and date1dk_correct(  half_life = 10,  time_unit = "y",  date1 = "2009-01-01",  date2 = c(    "2015-01-01",    "2016-01-01",    "2017-01-01"  ))#   half_life and time_lapsedk_correct(  half_life = 10,  time_lapse = 10,  time_unit = "y")# decay to todaydk_correct(RN_select = "Sr-90", date1 = "2009-01-01")# reverse decay - find out what readings should have been in the past given today's reading of 3000dk_correct(  RN_select = "Sr-90",  date1 = "2019-01-01",  date2 = c("2009-01-01", "1999-01-01"),  A1 = 3000)

Number of half-lives past

Description

Given a percentage reduction in activity, calculate how manyhalf-lives have passed.

Usage

dk_pct_to_num_half_life(pct_lost)

Arguments

Value

Number of half-lives passed.

See Also

Other decay corrections:dk_correct(),dk_time()

Examples

dk_pct_to_num_half_life(pct_lost = 93.75)

Time to decay to target radioactivity.

Description

Calculate time for a radionuclide to decay to a target activity.

Usage

dk_time(half_life, A0, A1)

Arguments

Value

Time, in same units as half-life, to decay to target activity.

See Also

Other decay corrections:dk_correct(),dk_pct_to_num_half_life()

Examples

# A carbonaceous artifact has a C-14 measurement of 1 dpm per g pure carbon.# The reference activity is 14 dpm per g pure carbon. How old is our sample?dk_time(half_life = 5730, A0 = 14, A1 = 1)

Calculate half-life based on two data points

Description

Estimate half-life from two data points. Half-life units areconsistent with time units of input.@family rad measurements

Usage

half_life_2pt(time1, time2, N1, N2)

Arguments

Value

The calculated half-life in units of time input.

Examples

# Between the first two data points in a series of countshalf_life_2pt(time1 = 0, time2 = 1, N1 = 45, N2 = 30)## Between the second and third in the series (same intervals)half_life_2pt(time1 = 1, time2 = 2, N1 = 30, N2 = 21)## Use on a seriescount_times <- 1:5acts <- 10000 * 2^(-count_times / 10) # activitiesacts <- rpois(5, acts) # activities with counting variability appliedhalf_life_2pt(  time1 = count_times[1:4], time2 = count_times[2:5],  N1 = acts[1:4], N2 = acts[2:5])

half-value layer and tenth-value layer computations

Description

Derive hvl and tvl from radiation values through a material thickness.

Usage

hvl(x, y)

Arguments

Value

a data frame with the inputs, followed by the computed values for attenuation coefficient (listed as "mu"), half-value layer (hvl), tenth-value layer (tvl), and the homogeneity coefficient (hc) which is the ratio of a half-value layer to the following half-value layer.

Examples

H50_ex <- data.frame("mm_Al" = 0:5, "mR_h" = c(7.428, 6.272, 5.325,4.535, 3.878, 3.317))hvl(x = H50_ex$mm_Al, y = H50_ex$mR_h)

MCNP Cone Opening Parameter

Description

MCNP cone surface requires a term, t^2, which is the tangent ofthe cone angle, in radians, squared. This function takes an input indegrees and provides the parameter needed by MCNP.

Usage

mcnp_cone_angle(d)

Arguments

Value

tangent of cone angle squared

See Also

Other mcnp tools:mcnp_est_nps(),mcnp_matrix_rotations(),mcnp_mesh_bins(),mcnp_plot_out_spec(),mcnp_scan2plot(),mcnp_scan_save(),mcnp_sdef_erg_hist(),mcnp_sdef_erg_line()

Examples

mcnp_cone_angle(45)

Copy and paste MCNP tally fluctuation charts

Description

Provides quick estimate of number of particles histories,(nps) to obtain target MCNP 'error'.Paste may include up to three tallies side by side in the default MCNPorder. For example, the headers of a three tally report includes columnnames: nps, mean, error, vov, slope, fom, mean, error, vov, slope, fom,mean, error, vov, slope, fom.The structure of the tfc has been the same for versions 4 through 6,including MCNPX.

Usage

mcnp_est_nps(err_target)

Arguments

Value

estimate of number of particle histories needed

See Also

Other mcnp tools:mcnp_cone_angle(),mcnp_matrix_rotations(),mcnp_mesh_bins(),mcnp_plot_out_spec(),mcnp_scan2plot(),mcnp_scan_save(),mcnp_sdef_erg_hist(),mcnp_sdef_erg_line()

Examples

# Since this function requires the user# to copy and paste input, this example# is set up to provide data for this purpose.# To run the example, copy and paste the following# into an input file and delete the hash tags to run.# Enter '1' for number of tallies.# mcnp_est_nps(0.01)#      32768000   4.5039E+00 0.2263 0.0969  0.0 5.0E-02#      65536000   3.9877E+00 0.1561 0.0553  0.0 5.1E-02#      98304000   3.4661E+00 0.1329 0.0413  0.0 4.7E-02#     131072000   3.5087E+00 0.1132 0.0305  0.0 5.0E-02#     163840000   3.5568E+00 0.0995 0.0228  0.0 5.2E-02#     196608000   3.8508E+00 0.0875 0.0164  0.0 5.5E-02#     229376000   3.8564E+00 0.0810 0.0135  0.0 5.5E-02#     262144000   3.9299E+00 0.0760 0.0118  0.0 5.5E-02#     294912000   4.0549E+00 0.0716 0.0100  0.0 5.6E-02#     327680000   4.0665E+00 0.0686 0.0090  0.0 5.4E-02#     360448000   4.1841E+00 0.0641 0.0079  0.0 5.7E-02

Rotation matrices for transformations in MCNP

Description

Create 3 x 3 rotation matrix in cosines of the anglesbetween the main and auxiliary coordinate systems in the form:xx' yx' zx'xy' yy' zy'xz' yz' zz'

Usage

mcnp_matrix_rotations(rot.axis, angle_degrees)

Arguments

Value

rotational matrix for copy and paste to MCNP input

See Also

Other mcnp tools:mcnp_cone_angle(),mcnp_est_nps(),mcnp_mesh_bins(),mcnp_plot_out_spec(),mcnp_scan2plot(),mcnp_scan_save(),mcnp_sdef_erg_hist(),mcnp_sdef_erg_line()

Examples

mcnp_matrix_rotations("x", 30)mcnp_matrix_rotations("y", 7)mcnp_matrix_rotations("z", 15)# For combined rotations, use matrix multiplication (%*%)# rotate 45 degrees on x-axis and 45 degrees on y-axismcnp_matrix_rotations("x", 45) %*% mcnp_matrix_rotations("y", 45)

Make mesh tally size settings for MCNP

Description

#':Find the parameters needed for a rectilinear "superimposed meshtally b" in MCNP. It can be a challenge to center mesh tallybins at a desired value of x, y, or z. This function looks at asingle dimension, – in units of cm – at a time. This is a newfunction and hasn't been tested thoroughly. The idea is toidentify a single setting in the MCNP mesh tally forimesh and iints (or jmesh and jints or kmesh and kints). Itis designed only for uniform mesh bin sizes.

Usage

mcnp_mesh_bins(  target,  width,  lowest_less,  lowest_high,  highest_high,  highest_less)

Arguments

Value

a data frame providing:

low_set, the minimum dimension. This is probably best used inthe origin parameter in the MCNP mesh tally.high_set, the maximum dimension for the bin. This can beidentified in the MCNP mesh tally setting of imesh, jmesh, orkmesh.width, this is just a return of the parameter supplied to thefunction.numblocks, the number of fine meshes. This can be used in theMCNP mesh tally setting of iints, jints, or kints.

See Also

Other mcnp tools:mcnp_cone_angle(),mcnp_est_nps(),mcnp_matrix_rotations(),mcnp_plot_out_spec(),mcnp_scan2plot(),mcnp_scan_save(),mcnp_sdef_erg_hist(),mcnp_sdef_erg_line()

Examples

mcnp_mesh_bins(target = 30, width = 10, lowest_less = 0, highest_less = 15, highest_high = 304.8, lowest_high = 250)#'

Convert histogram data to average points and plot as spectrum.

Description

Model results or input source histograms from MCNP and perhapsother sources typically provide binned tally results with columns representingmaximum energy in MeV, a column with the mean tally result or bin probabilityand an uncertainty column (not used). Once the data is scanned in, or otherwiseentered 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

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

Examples

mcnp_plot_out_spec(photons_cs137_hist, "example Cs-137 well irradiator")

Copy and paste MCNP output spectral data to directly plot

Description

Provides quick copy-and-paste to plot.Paste either a source histogram distribution or tally spectrum from MCNP outputs.Three-column output tally spectra have columns of maximum energy, bin tally, andrelative Monte Carlo uncertainty for the bin tally value.Four-column source histogram distributions have columns of entry number, maximumenergy, cumulative probability, and bin probability.In either case, only the maximum energy and bin probability or result values are used.

Usage

mcnp_scan2plot(title = "", log_plot = FALSE)

Arguments

Value

spectrum file with maximum energy and MCNP bin value

See Also

Other mcnp tools:mcnp_cone_angle(),mcnp_est_nps(),mcnp_matrix_rotations(),mcnp_mesh_bins(),mcnp_plot_out_spec(),mcnp_scan_save(),mcnp_sdef_erg_hist(),mcnp_sdef_erg_line()

Examples

# Since this function requires the user# to copy and paste input, this three column example# is set up to provide data for this purpose.# To run the example, copy and paste the following# into an input file and delete the hash tags to run.# mcnp_scan2plot(title = "example1")# 0.1000000 3.133122e-05 0.3348260# 0.4222222 6.731257e-05 0.2017546# 0.7444444 5.249198e-05 0.4524577# 1.0666667 2.046046e-04 0.4201954# 1.3888889 1.525125e-03 0.8049388# 1.7111111 2.922743e-05 0.7985399# 2.0333333 5.162954e-03 0.1974694# 2.3555556 2.048186e-05 0.5011170# 2.6777778 1.468040e-04 0.7248116# 3.0000000 1.037092e-04 0.7659850

Copy and paste MCNP output spectral data for use withmcnp_plot_out_spec()

Description

Provides quick copy-and-paste conversion to data frame.Paste either a source histogram distribution or tally spectrum from MCNP outputs.Three-column output tally spectra have columns of maximum energy, bin tally, andrelative Monte Carlo uncertainty for the bin tally value.Four-column source histogram distributions have columns of entry number, maximumenergy, cumulative probability, and bin probability.Seven-column biased histogram distributions have columns of entrynumber, maximum energy, cumulative probability, biased cumulativeprobability, probability of bin, biased probability, and weightmultiplier.In all cases, only the maximum energy and bin probability or result values are used.

Usage

mcnp_scan_save()

Value

spectrum file with maximum energy and MCNP bin value

See Also

Other mcnp tools:mcnp_cone_angle(),mcnp_est_nps(),mcnp_matrix_rotations(),mcnp_mesh_bins(),mcnp_plot_out_spec(),mcnp_scan2plot(),mcnp_sdef_erg_hist(),mcnp_sdef_erg_line()

Examples

# Since this function requires the user# to copy and paste input, this example# is set up to provide data for this purpose.# To run the example, copy and paste the following# into an input file and delete the hash tags to run.# my_hist_data <- mcnp_scan_save()# 0.1000000 3.133122e-05 0.3348260# 0.4222222 6.731257e-05 0.2017546# 0.7444444 5.249198e-05 0.4524577# 1.0666667 2.046046e-04 0.4201954# 1.3888889 1.525125e-03 0.8049388# 1.7111111 2.922743e-05 0.7985399# 2.0333333 5.162954e-03 0.1974694# 2.3555556 2.048186e-05 0.5011170# 2.6777778 1.468040e-04 0.7248116# 3.0000000 1.037092e-04 0.7659850

energy distribution histogram from pasted data

Description

energy distribution histogram from pasted data

Usage

mcnp_sdef_erg_hist(  entry_mode = "scan",  my_dir = NULL,  E_MeV = NULL,  bin_prob = NULL,  write_permit = "n",  log_plot = 0)

Arguments

Details

The output includes si# and sp#. The # should be changed to the appropriate distribution number. The data is saved in the global environment and appended to a file in the user's working directory, si_sp.txt. Two plots of the data are provided to the plot window, one with two linear axes and one with two log axes.

Value

A vector of energy bins and probabilities for an energy distribution, formatted as needed for MCNP input. It is designed for copying and pasting into an MCNP input.

See Also

mcnp_sdef_erg_line() for data fromRadData

Other mcnp tools:mcnp_cone_angle(),mcnp_est_nps(),mcnp_matrix_rotations(),mcnp_mesh_bins(),mcnp_plot_out_spec(),mcnp_scan2plot(),mcnp_scan_save(),mcnp_sdef_erg_line()

Examples

## Not run: mcnp_sdef_erg_hist()## End(Not run)

Produce MCNP source terms from ICRP 107 data

Description

Obtain emission data from the RadData package and write to a file for use with the radiation transport code, MCNP.

Usage

mcnp_sdef_erg_line(  desired_RN,  rad_type = NULL,  photon = FALSE,  cut = 0.001,  erg.dist = 1,  my_dir = NULL,  write_permit = "n")

Arguments

Value

a data frame can be saved to memory if desired (i.e. by my_file <- mcnp_sdef_erg_line(...)).For use with MCNP, a text file, 'si_sp.txt' is written to working directory.If file already exists, it is appended. The file contains allemission energies in the si 'card' and the Line indicator, L is included,e.g. si1 L 0.01 (showing a first energy of 0.01 MeV).This is followed by the emission probability of each si entry.An additional text entry is made summing up the probabilities.

See Also

mcnp_sdef_erg_hist() if radioactive emission data is available in histogram form and needs formatting for MCNP input.

Other mcnp tools:mcnp_cone_angle(),mcnp_est_nps(),mcnp_matrix_rotations(),mcnp_mesh_bins(),mcnp_plot_out_spec(),mcnp_scan2plot(),mcnp_scan_save(),mcnp_sdef_erg_hist()

Examples

## Not run: mcnp_sdef_erg_line("Co-60", photon = TRUE, cut = 0.01, erg.dist = 13)mcnp_sdef_erg_line("Sr-90", rad_type = "B-", cut = 0.01, erg.dist = 15)mcnp_sdef_erg_line("Am-241", rad_type = "A", cut = 0.01, erg.dist = 23)## End(Not run)

Solid Angle Correction for Neutron Detectors with Point Source

Description

Correction factors are needed when an Neutron Rem Detector (NRD)aka "Remball" is used in close proximity to a points source. This formulais per ISO ISO 8529-2-2000 section 6.2. Note, however, that the ISO formulapredicts the response. The formula used here takes the inverse to correctfor the over-response.

Usage

neutron_geom_cf(l, r.d, del = 0.5)

Arguments

Value

The correction factor for solid angle.

See Also

Other rad measurements:air_dens_cf(),disk_to_disk_solid_angle(),scaler_sim(),tau_estimate()

Examples

neutron_geom_cf(l = 11.1, r.d = 11)neutron_geom_cf(30, 11)neutron_geom_cf(5, 4.5)

File Description:

Description

This data file was generated in MCNP from a model of Gamma Well Irradiatorwith no attenuator in place. MCNP will include in the output a histogram oftally results when there is an E Tally Energy card. Results in the output upto MCNP version 6 have no headers, but the columns are:

Usage

photons_cs137_hist

Format

Adata.frame

E_max

Maximum Energy in MeV

bin_tally

Tally result for this bin

R

Monte Carlo uncertainty for this bin

Ratemeter Simulation

Description

Plot simulated ratemeter readings once per second for 600seconds. The meter starts with a reading of zero and builds up based on thetime constant. Resolution uncertainty is established to express theuncertainty from reading an analog scale, including the instability of itsreadings. Many standard references identify the precision or resolutionuncertainty of analog readings as half of the smallest increment. Thisshould be considered the single coverage uncertainty for a very stablereading. When a reading is not very stable, evaluation of the readingfluctuation is evaluated in terms of numbers of scale increments covered bymeter indication over a reasonable evaluation period.

Usage

rate_meter_sim(  cpm_equilibrium,  meter_scale_increments,  trials = 600,  tau = 9.5,  log_opt = "")

Arguments

Value

Plot of simulated meter reading every second..

Examples

rate_meter_sim(cpm_equilibrium = 270, meter_scale_increments = seq(100, 1000, 20))rate_meter_sim(cpm_equilibrium = 2.7e5, meter_scale_increments = seq(2e5, 1e6, 2e4))rate_meter_sim(450, seq(20, 1000, 20), trials = 1200, tau = 24.8534)

Count Room Scaler Simulation

Description

Returns a plotted distribution of results for a scaler modelbased on the Poisson distribution. Inputs and outputs in counts per minute.

Usage

scaler_sim(true_bkg, true_samp, ct_time, trials = 1e+05)

Arguments

Value

A histogram of all trial results including limits for +/- 1 standarddeviation.

See Also

Other rad measurements:air_dens_cf(),disk_to_disk_solid_angle(),neutron_geom_cf(),tau_estimate()

Examples

scaler_sim(true_bkg = 5, true_samp = 10, ct_time = 1, trials = 1e5)scaler_sim(true_bkg = 50, true_samp = 30, ct_time = 1, trials = 1e5)

Stay time for radiation work.

Description

Calculate stay time for radiation work.

Usage

stay_time(dose_rate, dose_allowed, margin = 20)

Arguments

Value

Time in minutes allowed for the work.

Examples

stay_time(dose_rate = 100, dose_allowed = 50, margin = 20)

Estimate tau parameter forratemeter_sim

Description

If the time constant is not known, but the vendor specifies thatthe ratemeter will reach some percentage of equilibrium in some number ofseconds, use this function to estimate tau.

Usage

tau_estimate(pct_eq, t_eq)

Arguments

Value

tau, the time constant, in seconds.

See Also

Other rad measurements:air_dens_cf(),disk_to_disk_solid_angle(),neutron_geom_cf(),scaler_sim()

Examples

tau_estimate(pct_eq = 90, t_eq = 22)