biodosetools (version 3.7.2)

fun.estimate.criticality: Calculate absorbed dose for mixed fields exposure

Description

Calculate absorbed dose for mixed fields exposure

Usage

fun.estimate.criticality(
  num_cases,
  dics,
  cells,
  coef_gamma,
  cov_gamma,
  coef_neutron,
  cov_neutron,
  ratio,
  p
)

Value

dose estimation

Arguments

num_cases

number of cases to estimate.

dics

dicentrics.

cells

total cells.

coef_gamma

Coefficients for gamma curve.

cov_gamma

Covariance matrix for gamma curve.

coef_neutron

Coefficients for neutron curve.

cov_neutron

Covariance matrix for neutron curve.

ratio

gamma-neutron ratio.

p

photon contribution in mixed curve.

Examples

#Results from the fitting module
fit_results_gamma <- system.file("extdata",
                                 "gamma_dicentrics-fitting-results.rds",
                                 package = "biodosetools") %>%
  readRDS()

fit_results_neutrons <- system.file("extdata",
                                    "neutrons-mixed-dicentrics-fitting-results.rds",
                                    package = "biodosetools") %>%
  readRDS()

#FUNCTION TO ESTIMATE DOSE IN CRITICALITY ACCIDENTS
fun.estimate.criticality(num_cases = 2,
                         dics = c(380,456),
                         cells = c(218,567),
                         coef_gamma = fit_results_gamma[["fit_coeffs"]][,1],
                         cov_gamma = fit_results_gamma[["fit_var_cov_mat"]],
                         coef_neutron = fit_results_neutrons[["fit_coeffs"]][,1],
                         cov_neutron = fit_results_neutrons[["fit_var_cov_mat"]],
                         ratio = 1.2,
                         p = 0)