From 8999b0973a55c969bf9713134ec127e6b226b541 Mon Sep 17 00:00:00 2001 From: RLumSK Date: Wed, 16 Jul 2025 22:23:11 +0200 Subject: [PATCH 1/5] Add Heidelberg reference dataset + ad raw data + ad data conversion + ad example data + ad docu + up manuals + ad NEWS --- NEWS.md | 4 +++ R/data.R | 33 +++++++++++++++++++++++ data-raw/DATASET.R | 5 ++++ data-raw/heidelberg_2024.csv | 9 +++++++ data/datalist | 1 + data/heidelberg_2024.rda | Bin 0 -> 1227 bytes man/AIX_NaI_1.Rd | 3 ++- man/BDX_LaBr_1.Rd | 3 ++- man/clermont.Rd | 3 ++- man/clermont_2024.Rd | 3 ++- man/heidelberg_2024.Rd | 49 +++++++++++++++++++++++++++++++++++ 11 files changed, 109 insertions(+), 4 deletions(-) create mode 100644 data-raw/heidelberg_2024.csv create mode 100644 data/heidelberg_2024.rda create mode 100644 man/heidelberg_2024.Rd diff --git a/NEWS.md b/NEWS.md index 49a63704..902c2ba8 100644 --- a/NEWS.md +++ b/NEWS.md @@ -2,6 +2,10 @@ ## Additions * Add user friendly coercion methods: `as_peaks()` and `as_spectra()`. +## Datasets +* Add reference datasets for calibration sites around Heidelberg as detailed +in Kreutzer et al. (2025, Geochronology) + # gamma 1.1.0 ## Bugfixes * Fix an error in the uncertainty calculation of `dose_predict()`. The returned error was too large and did not make much sense due to an internal calculation error. Along with the fix, the manual was updated to detail the uncertainty calculation (PR #42 by @RLumSK). diff --git a/R/data.R b/R/data.R index 30f32efc..e28bb3c7 100644 --- a/R/data.R +++ b/R/data.R @@ -104,6 +104,39 @@ NULL #' @keywords datasets "clermont_2024" +# ============================================================================== +#' Heidelberg Reference Data 2024 +#' +#' Reference data for the calibration sites near Heidelberg as detailed +#' in Kreutzer at al. (2025) +#' +#' @usage data("heidelberg_2024") +#' @format A [`data.frame`] with tabulated data and conversion factor reference. +#' @source +#' Adamiec, G. & Aitken, M.J. (1998). +#' Dose-rate conversion factors: update. *Ancient TL*, 16, p. 37–50. +#' +#' Cresswell, A.J., Carter, J. & Sanderson, D.C.W. (2018). Dose rate +#' conversion parameters: Assessment of nuclear data. +#' *Radiation Measurements*, 120, p. 195–201. +#' +#' Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: +#' Update. *Ancient TL*, 29(1), p. 5-8. +#' +#' Liritzis, I., Stamoulis, K., Papachristodoulou, C. & Ioannides, K. (2013). +#' A re-evaluation of radiation dose-rate conversion factors. +#' *Mediterranean Archaeology and Archaeometry*, 12, p. 1–15. +#' +#' Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. +#' (2009). The Clermont Radiometric Reference Rocks: A Convenient Tool +#' for Dosimetric Purposes. *Ancient TL*, 27(2), p. 37-44. +#' @family datasets +#' @references Kreutzer, S., Martin, L., Miallier, D., Mercier, N., 2025. +#' Environmental gamma dose rate measurements using cadmium zinc telluride (CZT) detectors. +#' Geochronology 7, 229–246. \doi{10.5194/gchron-7-229-2025} +#' @keywords datasets +"heidelberg_2024" + #' CEREGE Calibration Curve (NaI) #' #' @usage data(AIX_NaI_1) diff --git a/data-raw/DATASET.R b/data-raw/DATASET.R index bbfbb4a2..e9b0e571 100644 --- a/data-raw/DATASET.R +++ b/data-raw/DATASET.R @@ -17,6 +17,11 @@ clermont <- utils::read.table("./data-raw/clermont.csv", rownames(clermont) <- clermont[["name"]] usethis::use_data(clermont, internal = FALSE, overwrite = FALSE) +# Milieux Heidelberg =========================================================== +heidelberg_2024 <- utils::read.table("./data-raw/heidelberg_2024.csv", + header = TRUE, sep = ",", dec = ".") +usethis::use_data(heidelberg_2024, internal = FALSE, overwrite = FALSE) + # Decay data =================================================================== .decay <- utils::read.table("./data-raw/decay.csv", header = TRUE, sep = ";", dec = ",") |> diff --git a/data-raw/heidelberg_2024.csv b/data-raw/heidelberg_2024.csv new file mode 100644 index 00000000..c204e67b --- /dev/null +++ b/data-raw/heidelberg_2024.csv @@ -0,0 +1,9 @@ +name,laboratory,nature,U,U_error,Th,Th_error,K2O,K2O_error,K,K_error,gamma_dose,gamma_error,cosmic_dose,cosmic_error,ref_conversion_factors,conv_U,conv_U_error,conv_Th,conv_Th_error,conv_K,conv_K_error +FLOSSI,"Max-Planck-Institute for Nuclear Physics, Heidelberg, DE",granite,19.02,0.74,13.25,0.50,,,4.11,0.095,3778.690,92.096,203,20,Adamiec & Aitken (1998),0.113000000,0.000202147,0.047600000,0.000197921,0.243000000,0.004669335 +WH2024,"Weisse Hohl (natural site), DE",loess,3.09,0.02,9.44,0.53,,,1.20,0.01,1090.114,26.129,141,14,Adamiec & Aitken (1998),0.113000000,0.000202147,0.047600000,0.000197921,0.243000000,0.004669335 +FLOSSI,"Max-Planck-Institute for Nuclear Physics, Heidelberg, DE",granite,19.02,0.74,13.25,0.50,,,4.11,0.095,3781.108,91.445,203,20,Guérin et al. (2011),0.111600000,0.000199642,0.047900000,0.000199168,0.249100000,0.004786549 +WH2024,"Weisse Hohl (natural site), DE",loess,3.09,0.02,9.44,0.53,,,1.20,0.01,1095.940,26.435,141,14,Guérin et al. (2011),0.111600000,0.000199642,0.047900000,0.000199168,0.249100000,0.004786549 +FLOSSI,"Max-Planck-Institute for Nuclear Physics, Heidelberg, DE",granite,19.02,0.74,13.25,0.50,,,4.11,0.095,3790.439,91.634,203,20,Liritzis et al. (2013),0.111800000,0.000200000,0.048100000,0.000200000,0.249800000,0.004800000 +WH2024,"Weisse Hohl (natural site), DE",loess,3.09,0.02,9.44,0.53,,,1.20,0.01,1099.286,26.543,141,14,Liritzis et al. (2013),0.111800000,0.000200000,0.048100000,0.000200000,0.249800000,0.004800000 +FLOSSI,"Max-Planck-Institute for Nuclear Physics, Heidelberg, DE",granite,19.02,0.74,13.25,0.50,,,4.11,0.095,3797.445,92.475,203,20,Cresswell et al. (2018),0.112000000,0.001000000,0.048900000,0.000300000,0.248000000,0.003000000 +WH2024,"Weisse Hohl (natural site), DE",loess,3.09,0.02,9.44,0.53,,,1.20,0.01,1105.296,26.826,141,14,Cresswell et al. (2018),0.112000000,0.001000000,0.048900000,0.000300000,0.248000000,0.003000000 \ No newline at end of file diff --git a/data/datalist b/data/datalist index 0ebee885..59effdb8 100644 --- a/data/datalist +++ b/data/datalist @@ -2,3 +2,4 @@ AIX_NaI_1: AIX_NaI_1 BDX_LaBr_1: BDX_LaBr_1 clermont_2024: clermont_2024 clermont: clermont +heidelberg_2024: heidelberg_2024 diff --git a/data/heidelberg_2024.rda b/data/heidelberg_2024.rda new file mode 100644 index 0000000000000000000000000000000000000000..26cab643a2cd1fb4d6d4e40c0008d1f382661169 GIT binary patch literal 1227 zcmZ>Y%CIzaj8qGb?3~27l7Z=0{r~^}-|zkU{qHXzuTf@c4 z4m458km({I#bA=bz$?Xo308S1x;M@1%25tZS^nZw)eK+%;0yg8_90*9&6?-=Dx)ss z;;*dPn#-;n%QW$QH9K%qs;$2hf7ZO^lcdsDr7g%@aN*PJz-}V zIx3ynx<2!q)PrJN{JjN@-5j*^_MKDSFy&96@9a0%w>qX}OK|XQI-+em*JOr&TZ-C5 zp$AHj7&+xHUs;~ubYw3NpV_w-7spMLW~uC6R1})D&oSsd`>T%{i?oW2zGh{g3SI3R zG)L7lY;vri>E_*2IBqF&stO8qw3tn5Hu60C=PsvH#zf{j)t;+EmRt-Je44M;Q6$t- zt?AgiH7TeYbYL6S8NV@=$QQ+2Sj%p{b$jxrRsoo?(9RQT4G+$;?Ya6JILJb$aib zsjV(6o#e4D(azpxO6;uv=98+@UH)=Esn}`px8sS6?(CzlC-iK5k$g(M?eBgA?$7ce zb225k1AgUgj;Naxb2nz`eeurAJGj%Ad{V64koxr9v#C?QJUIVt=~P?k=Wpi4e69)J zT_#tivSXv4zvjc*PJgHW!d`hQmma%Xs#eTiQOndZzpe05o{jhFnXj+BGk#T;L&Y838FFEv}vi8fJ{lO|vjn2*9 zu&tu<$PUfZZ4WGK55AZmTw1zi>XS{rZ=d}ByyCEJaMMJ|-O5w4_6WRuB2w+Jep1A_ zMXq%`78PHsRVrI~H)smqb;w{zip znjI4UO=9z*MW0^%nQbb}y=c)q&*r{(=W~_snZy!q-bxm|p>gP4)|uNEinE-f>l^a> z<~BDRwhN9ppy!m@@?ZCAP?Cs*kJ9Qt&RKpRX1GbGw7xyGc814fi5)AVJQi-wu31}~ zr`a1((zQ|MZQ|XjZ;!+|wk5@FoU-za@xzYI9NV(YS$j|V{F66a*KS&}mF-xEwpX;s zZI+d@3`4_HL)m3kmISOV)mY=%5teoG;v=mPJs~Zjh=sz+55CKArp*vY)pf~}t%%vQ zp<3_Ww(zB8i3<;Z_jLKt|LnQ#4Ugwa`ttFAO>XISUQxb(!n^d})$@&V*=^zZ;%D`D zY)ROn8(O3PlDlu`?;_b*kFQ0PM!q%+PrE1OvQXmOrLxYOAqBk?b7o}Pq_gj7RsOLz z<8qes{TKGa2kYaFS)Q34{}Uw>KI6`dxiZhwANEhWD$AagY$V0faeec-dvl$)Y`j*X up3|!yxwEdk;_{W!TitGp%>||y$N#R{7j1DSRI~nLxk_*=C&x?$sR;nT-8~cl literal 0 HcmV?d00001 diff --git a/man/AIX_NaI_1.Rd b/man/AIX_NaI_1.Rd index ebb0f6a0..e89cad37 100644 --- a/man/AIX_NaI_1.Rd +++ b/man/AIX_NaI_1.Rd @@ -28,7 +28,8 @@ plot(AIX_NaI_1) Other datasets: \code{\link{BDX_LaBr_1}}, \code{\link{clermont}}, -\code{\link{clermont_2024}} +\code{\link{clermont_2024}}, +\code{\link{heidelberg_2024}} } \concept{datasets} \keyword{datasets} diff --git a/man/BDX_LaBr_1.Rd b/man/BDX_LaBr_1.Rd index 47f132d1..c9ca8e41 100644 --- a/man/BDX_LaBr_1.Rd +++ b/man/BDX_LaBr_1.Rd @@ -28,7 +28,8 @@ plot(BDX_LaBr_1) Other datasets: \code{\link{AIX_NaI_1}}, \code{\link{clermont}}, -\code{\link{clermont_2024}} +\code{\link{clermont_2024}}, +\code{\link{heidelberg_2024}} } \concept{datasets} \keyword{datasets} diff --git a/man/clermont.Rd b/man/clermont.Rd index 386dbfa7..5fd4cc2e 100644 --- a/man/clermont.Rd +++ b/man/clermont.Rd @@ -25,7 +25,8 @@ Clermont Reference Data Other datasets: \code{\link{AIX_NaI_1}}, \code{\link{BDX_LaBr_1}}, -\code{\link{clermont_2024}} +\code{\link{clermont_2024}}, +\code{\link{heidelberg_2024}} } \concept{datasets} \keyword{datasets} diff --git a/man/clermont_2024.Rd b/man/clermont_2024.Rd index 05730315..be1f986a 100644 --- a/man/clermont_2024.Rd +++ b/man/clermont_2024.Rd @@ -41,7 +41,8 @@ errors. Other datasets: \code{\link{AIX_NaI_1}}, \code{\link{BDX_LaBr_1}}, -\code{\link{clermont}} +\code{\link{clermont}}, +\code{\link{heidelberg_2024}} } \concept{datasets} \keyword{datasets} diff --git a/man/heidelberg_2024.Rd b/man/heidelberg_2024.Rd new file mode 100644 index 00000000..2f04c118 --- /dev/null +++ b/man/heidelberg_2024.Rd @@ -0,0 +1,49 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/data.R +\docType{data} +\name{heidelberg_2024} +\alias{heidelberg_2024} +\title{Heidelberg Reference Data 2024} +\format{ +A \code{\link{data.frame}} with tabulated data and conversion factor reference. +} +\source{ +Adamiec, G. & Aitken, M.J. (1998). +Dose-rate conversion factors: update. \emph{Ancient TL}, 16, p. 37–50. + +Cresswell, A.J., Carter, J. & Sanderson, D.C.W. (2018). Dose rate +conversion parameters: Assessment of nuclear data. +\emph{Radiation Measurements}, 120, p. 195–201. + +Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: +Update. \emph{Ancient TL}, 29(1), p. 5-8. + +Liritzis, I., Stamoulis, K., Papachristodoulou, C. & Ioannides, K. (2013). +A re-evaluation of radiation dose-rate conversion factors. +\emph{Mediterranean Archaeology and Archaeometry}, 12, p. 1–15. + +Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. +(2009). The Clermont Radiometric Reference Rocks: A Convenient Tool +for Dosimetric Purposes. \emph{Ancient TL}, 27(2), p. 37-44. +} +\usage{ +data("heidelberg_2024") +} +\description{ +Reference data for the calibration sites near Heidelberg as detailed +in Kreutzer at al. (2025) +} +\references{ +Kreutzer, S., Martin, L., Miallier, D., Mercier, N., 2025. +Environmental gamma dose rate measurements using cadmium zinc telluride (CZT) detectors. +Geochronology 7, 229–246. \doi{10.5194/gchron-7-229-2025} +} +\seealso{ +Other datasets: +\code{\link{AIX_NaI_1}}, +\code{\link{BDX_LaBr_1}}, +\code{\link{clermont}}, +\code{\link{clermont_2024}} +} +\concept{datasets} +\keyword{datasets} From b8f4b1e204684e6cfa2773a950c824f0810bb1f4 Mon Sep 17 00:00:00 2001 From: RLumSK Date: Wed, 16 Jul 2025 22:29:42 +0200 Subject: [PATCH 2/5] Update documentation and update references related DOIs --- R/data.R | 21 ++++++++++++++++----- man/clermont.Rd | 2 ++ man/clermont_2024.Rd | 4 ++++ man/heidelberg_2024.Rd | 15 +++++++++------ 4 files changed, 31 insertions(+), 11 deletions(-) diff --git a/R/data.R b/R/data.R index e28bb3c7..eb189cf9 100644 --- a/R/data.R +++ b/R/data.R @@ -62,10 +62,12 @@ NULL #' @source #' Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: #' Update. *Ancient TL*, 29(1), p. 5-8. +#' \doi{10.26034/la.atl.2011.443} #' #' Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. #' (2009). The Clermont Radiometric Reference Rocks: A Convenient Tool #' for Dosimetric Purposes. *Ancient TL*, 27(2), p. 37-44. +#' \doi{10.26034/la.atl.2009.428} #' @family datasets #' @keywords datasets "clermont" @@ -85,13 +87,16 @@ NULL #' @source #' Adamiec, G. & Aitken, M.J. (1998). #' Dose-rate conversion factors: update. *Ancient TL*, 16, p. 37–50. +#' \doi{10.26034/la.atl.1998.292} #' #' Cresswell, A.J., Carter, J. & Sanderson, D.C.W. (2018). Dose rate #' conversion parameters: Assessment of nuclear data. #' *Radiation Measurements*, 120, p. 195–201. +#' \doi{10.1016/j.radmeas.2018.02.007} #' #' Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: #' Update. *Ancient TL*, 29(1), p. 5-8. +#' \doi{10.26034/la.atl.2011.443} #' #' Liritzis, I., Stamoulis, K., Papachristodoulou, C. & Ioannides, K. (2013). #' A re-evaluation of radiation dose-rate conversion factors. @@ -100,6 +105,7 @@ NULL #' Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. #' (2009). The Clermont Radiometric Reference Rocks: A Convenient Tool #' for Dosimetric Purposes. *Ancient TL*, 27(2), p. 37-44. +#' \doi{10.26034/la.atl.2009.428} #' @family datasets #' @keywords datasets "clermont_2024" @@ -115,25 +121,30 @@ NULL #' @source #' Adamiec, G. & Aitken, M.J. (1998). #' Dose-rate conversion factors: update. *Ancient TL*, 16, p. 37–50. +#' \doi{10.26034/la.atl.1998.292} #' #' Cresswell, A.J., Carter, J. & Sanderson, D.C.W. (2018). Dose rate #' conversion parameters: Assessment of nuclear data. #' *Radiation Measurements*, 120, p. 195–201. +#' \doi{10.1016/j.radmeas.2018.02.007} #' #' Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: #' Update. *Ancient TL*, 29(1), p. 5-8. +#' \doi{10.26034/la.atl.2011.443} #' #' Liritzis, I., Stamoulis, K., Papachristodoulou, C. & Ioannides, K. (2013). #' A re-evaluation of radiation dose-rate conversion factors. #' *Mediterranean Archaeology and Archaeometry*, 12, p. 1–15. #' -#' Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. -#' (2009). The Clermont Radiometric Reference Rocks: A Convenient Tool -#' for Dosimetric Purposes. *Ancient TL*, 27(2), p. 37-44. #' @family datasets -#' @references Kreutzer, S., Martin, L., Miallier, D., Mercier, N., 2025. +#' @references +#' Kreutzer, S., Martin, L., Miallier, D. & Mercier, N. (2025). #' Environmental gamma dose rate measurements using cadmium zinc telluride (CZT) detectors. -#' Geochronology 7, 229–246. \doi{10.5194/gchron-7-229-2025} +#' *Geochronology*, 7, p. 229–246. \doi{10.5194/gchron-7-229-2025} +#' +#' Kreutzer, S., Martin, L., Miallier, D., Mercier, N. & Faske, M. (2024). +#' Dataset: Environmental Gamma Dose Rate Measurements using CZT Detectors (1.0.0). +#' *Zenodo*. \doi{10.5281/zenodo.13731839} #' @keywords datasets "heidelberg_2024" diff --git a/man/clermont.Rd b/man/clermont.Rd index 5fd4cc2e..f5ecb0f3 100644 --- a/man/clermont.Rd +++ b/man/clermont.Rd @@ -10,10 +10,12 @@ A \code{\link{data.frame}} with tabulated data and conversion factor reference. \source{ Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: Update. \emph{Ancient TL}, 29(1), p. 5-8. +\doi{10.26034/la.atl.2011.443} Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. (2009). The Clermont Radiometric Reference Rocks: A Convenient Tool for Dosimetric Purposes. \emph{Ancient TL}, 27(2), p. 37-44. +\doi{10.26034/la.atl.2009.428} } \usage{ data("clermont") diff --git a/man/clermont_2024.Rd b/man/clermont_2024.Rd index be1f986a..a93b0224 100644 --- a/man/clermont_2024.Rd +++ b/man/clermont_2024.Rd @@ -10,13 +10,16 @@ A \code{\link{data.frame}} with tabulated data and conversion factor reference. \source{ Adamiec, G. & Aitken, M.J. (1998). Dose-rate conversion factors: update. \emph{Ancient TL}, 16, p. 37–50. +\doi{10.26034/la.atl.1998.292} Cresswell, A.J., Carter, J. & Sanderson, D.C.W. (2018). Dose rate conversion parameters: Assessment of nuclear data. \emph{Radiation Measurements}, 120, p. 195–201. +\doi{10.1016/j.radmeas.2018.02.007} Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: Update. \emph{Ancient TL}, 29(1), p. 5-8. +\doi{10.26034/la.atl.2011.443} Liritzis, I., Stamoulis, K., Papachristodoulou, C. & Ioannides, K. (2013). A re-evaluation of radiation dose-rate conversion factors. @@ -25,6 +28,7 @@ A re-evaluation of radiation dose-rate conversion factors. Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. (2009). The Clermont Radiometric Reference Rocks: A Convenient Tool for Dosimetric Purposes. \emph{Ancient TL}, 27(2), p. 37-44. +\doi{10.26034/la.atl.2009.428} } \usage{ data("clermont_2024") diff --git a/man/heidelberg_2024.Rd b/man/heidelberg_2024.Rd index 2f04c118..d236463d 100644 --- a/man/heidelberg_2024.Rd +++ b/man/heidelberg_2024.Rd @@ -10,21 +10,20 @@ A \code{\link{data.frame}} with tabulated data and conversion factor reference. \source{ Adamiec, G. & Aitken, M.J. (1998). Dose-rate conversion factors: update. \emph{Ancient TL}, 16, p. 37–50. +\doi{10.26034/la.atl.1998.292} Cresswell, A.J., Carter, J. & Sanderson, D.C.W. (2018). Dose rate conversion parameters: Assessment of nuclear data. \emph{Radiation Measurements}, 120, p. 195–201. +\doi{10.1016/j.radmeas.2018.02.007} Guérin, G., Mercier, N. & Adamiec, G. (2011). Dose-Rate Conversion Factors: Update. \emph{Ancient TL}, 29(1), p. 5-8. +\doi{10.26034/la.atl.2011.443} Liritzis, I., Stamoulis, K., Papachristodoulou, C. & Ioannides, K. (2013). A re-evaluation of radiation dose-rate conversion factors. \emph{Mediterranean Archaeology and Archaeometry}, 12, p. 1–15. - -Miallier, D., Guérin, G., Mercier, N., Pilleyre, T. & Sanzelle, S. -(2009). The Clermont Radiometric Reference Rocks: A Convenient Tool -for Dosimetric Purposes. \emph{Ancient TL}, 27(2), p. 37-44. } \usage{ data("heidelberg_2024") @@ -34,9 +33,13 @@ Reference data for the calibration sites near Heidelberg as detailed in Kreutzer at al. (2025) } \references{ -Kreutzer, S., Martin, L., Miallier, D., Mercier, N., 2025. +Kreutzer, S., Martin, L., Miallier, D. & Mercier, N. (2025). Environmental gamma dose rate measurements using cadmium zinc telluride (CZT) detectors. -Geochronology 7, 229–246. \doi{10.5194/gchron-7-229-2025} +\emph{Geochronology}, 7, p. 229–246. \doi{10.5194/gchron-7-229-2025} + +Kreutzer, S., Martin, L., Miallier, D., Mercier, N. & Faske, M. (2024). +Dataset: Environmental Gamma Dose Rate Measurements using CZT Detectors (1.0.0). +\emph{Zenodo}. \doi{10.5281/zenodo.13731839} } \seealso{ Other datasets: From 0004cde38b4ad3d3f5a0395d4cf0943966243187 Mon Sep 17 00:00:00 2001 From: RLumSK Date: Wed, 16 Jul 2025 22:31:17 +0200 Subject: [PATCH 3/5] Update writing of nuclides in documentation --- R/data.R | 2 +- man/decay.Rd | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/R/data.R b/R/data.R index eb189cf9..61a2f2dc 100644 --- a/R/data.R +++ b/R/data.R @@ -2,7 +2,7 @@ #' Nuclear Decay Data #' #' An internal dataset containing the decay data for the following isotopes: -#' 232-Th, 235-U, 238-U and 40-K. +#' Th-232, U-235, U-238 and K-40. #' @format #' A [`data.frame`] with the following columns (absolute errors): #' \describe{ diff --git a/man/decay.Rd b/man/decay.Rd index db4958d9..b0d4180f 100644 --- a/man/decay.Rd +++ b/man/decay.Rd @@ -26,7 +26,7 @@ Nuclides And Isotopes : Chart of the Nuclides. } \description{ An internal dataset containing the decay data for the following isotopes: -232-Th, 235-U, 238-U and 40-K. +Th-232, U-235, U-238 and K-40. } \keyword{datasets} \keyword{internal} From 41bfcf10a9388e39aa364f3e23bcac90b62694d0 Mon Sep 17 00:00:00 2001 From: RLumSK Date: Wed, 16 Jul 2025 23:02:04 +0200 Subject: [PATCH 4/5] Add calibration datasets Heidelberg + ad calibration dataset + up documentation + up NEWS --- NEWS.md | 4 +-- R/data.R | 42 +++++++++++++++++++++++++++++ data/HD_CZT_1.rda | Bin 0 -> 45448 bytes data/datalist | 1 + man/AIX_NaI_1.Rd | 1 + man/BDX_LaBr_1.Rd | 1 + man/HD_CZT_1.Rd | 59 +++++++++++++++++++++++++++++++++++++++++ man/clermont.Rd | 1 + man/clermont_2024.Rd | 1 + man/heidelberg_2024.Rd | 1 + 10 files changed, 109 insertions(+), 2 deletions(-) create mode 100644 data/HD_CZT_1.rda create mode 100644 man/HD_CZT_1.Rd diff --git a/NEWS.md b/NEWS.md index 902c2ba8..3e45cef2 100644 --- a/NEWS.md +++ b/NEWS.md @@ -3,8 +3,8 @@ * Add user friendly coercion methods: `as_peaks()` and `as_spectra()`. ## Datasets -* Add reference datasets for calibration sites around Heidelberg as detailed -in Kreutzer et al. (2025, Geochronology) +* Add reference datasets for calibration sites around Heidelberg and the CZT detector calibration +files as detailed in Kreutzer et al. (2025, Geochronology) # gamma 1.1.0 ## Bugfixes diff --git a/R/data.R b/R/data.R index 61a2f2dc..943f2298 100644 --- a/R/data.R +++ b/R/data.R @@ -185,3 +185,45 @@ NULL #' @family datasets #' @keywords datasets "BDX_LaBr_1" + +#' heiLUM Calibration Curve (CZT) +#' +#' @usage data(HD_CZT_1) +#' @format +#' Objects of class [CalibrationCurve-class]. +#' \tabular{ll}{ +#' **Laboratory** \tab heiLUM, Heidelberg University, DE \cr +#' **Instrument** \tab Kromek GR1 & Kromek RayMon10 (SN: 5261)\cr +#' **Detector** \tab GR1+ (SN: 1625) & GR1 (SN: 1663) \cr +#' **Authors** \tab Sebastian Kreutzer, Didier Maillier, Norbert Mercier, heiLUM Team +#' } +#' +#' The calibration dataset consists of the following objects: +#' +#' * CAL_GR1_AdamiecAitken1998 +#' * CAL_GR1_Cresswelletal2018 +#' * CAL_GR1_Gurinetal2011 +#' * CAL_GR1_Liritzisetal2013 +#' * CAL_RAYMON_GR1_AdamiecAitken1998 +#' * CAL_RAYMON_GR1_Cresswelletal2018 +#' * CAL_RAYMON_GR1_Gurinetal2011 +#' * CAL_RAYMON_GR1_Liritzisetal2013 +#' +#' @examples +#' ## Load the curve +#' data(HD_CZT_1, package = "gamma") +#' names(HD_CZT_1) +#' plot(HD_CZT_1$CAL_GR1_AdamiecAitken1998) +#' +#' @references +#' Kreutzer, S., Martin, L., Miallier, D. & Mercier, N. (2025). +#' Environmental gamma dose rate measurements using cadmium zinc telluride (CZT) detectors. +#' *Geochronology*, 7, p. 229–246. \doi{10.5194/gchron-7-229-2025} +#' +#' Kreutzer, S., Martin, L., Miallier, D., Mercier, N. & Faske, M. (2024). +#' Dataset: Environmental Gamma Dose Rate Measurements using CZT Detectors (1.0.0). +#' *Zenodo*. \doi{10.5281/zenodo.13731839} +#' +#' @family datasets +#' @keywords datasets +"HD_CZT_1" diff --git a/data/HD_CZT_1.rda b/data/HD_CZT_1.rda new file mode 100644 index 0000000000000000000000000000000000000000..f3f7131f6aaae78f2ac3886276eb3391432165e0 GIT binary patch literal 45448 zcmV(lK=i-;H+ooF0004LBHlIv03iV!0000G&sfas&RnrfT>vQ&2UKVgRpfklJ z)Cw1%v#nn{1OYnOYyJAsqs;liV8N~ahYx;^hz9nn^RN!H1fTu0IdM}bUuqF8y3tJwaandPoNP<39+fGdM~j$#fYRY`4!kcbIXFPcal~FEYAi@I`8^BC$1v9 z=`j$><87N*pO)s3<8}$&dHGXu%%oU{U{1yrs5?9=OyXoM8!th&(QG!$OI*vwrKf;w zsM+S@{3%Q#7|ed{dX3nQR$KVN0C%S+-E)Zid3RHtv7JZdkhN+1jd+>R?M8)#t;{MO zP#B7YI}k|vU~`;WGSE;(SZ_!O4E$ zv3NA!nWda?8O`y#P^FmZHEwmZ73Y*J2?vK6eq!t=rap2npx~LZNNznma(^6*0D7u-L?HKAANfOS$LXBV7d1ql@INtsq9pAz|23o%u zf02hCyfa2%s~xb~GSf4X==0av4i=Nzd!VbaT^AIKeRnUU<7X?vfO znL~B`&EU4kyTV+-U>aXZcrs+lBrBIRa}O-2idn8Fnf8|?;za{?}>)FHrnG_8}Iv*B)J1l!JW4M_YEPMT@Jj57xn`_mM zmD_4zFVFTl18m}q&g?q2LF{tZ^Y0C{AG6zrT@MxO52RUcfMr0yKdC^v@w#akxzOfJx6e(`DMIw2;{`oCDCYv5g2-D}#=}quO$N zsaVY^$svKkS)+XV%)WC@v53_!6XG4g#iO}#zPqyMSj z@l|0qDogLHl=w#tw1E)4Fh3sh7#`)oF!fiGr%&0Bg>1Rf17u{4*WEN{V=K={1O~DQ zyWBY>cWEI}Z)0+Nlfe zW=6XO*0U@u%uFPuBG+mZ@oQ*yX?H*F3<=2jSto9EWUW<0g>v>16#dyu5qHZ`7Rk}Y znDS7H-7kDH83qA&L2TU1M7qYNHOE&RwWNA}xvg86iJH4-3|8_dsrARP5$E3l#<`0% zw03Bz;hbBjY?tjwu8Cxba=Qje>MOYq{FdC^MyHe`xA2+O1x}CATaeb8>a*0vjVZ&C z6RVym1j~*00_7jyV5<}@*NLz|=Q%N!4b@Y&>s0upW>7abY)9y3hz6-{Pyk;=b~2Z4 zTWxuK(A$oAZIuFjWQ(l!LIb|;8}GHKGreW=wz87#DlHzU%_EDqZbAXHCrV^X6<$ED zhIu%E+Vm7NJb5!Ml8HmgKsF05;7(=U@rk-D(4S6<9ZKq`O^;LjOKQ$;T$Pz!Z8YA& z!?wSaWBch|R1Ji-KqEZ(r53uX|716{CYIW=9O;4X-j+lsvQEw()Ta9D`T1wUy&%or}T^R5MpHT%EpHZk{I`gm)? zhD@d!MS+l~V_S#P&~534k5FO4`CkJ)4T3 z`Jv6Cddges7A#W|sW$)h1Z!0q&wI6}leR2x94C}@3SSdJ78=R2rNTJ4uJW1N(f6e4 zkRjFKiisj#GJv|HWN7+i;w5l&89v7fb-z{5h;v6@Ns>W%re$8^+?pjBZ{*VEiw0QE zepsF`%LrF-anVty)gDk#IW4E{Yt4IKqvWnNkTM)AV4-zkIGDJm`}b)2|H)B-~IF%N3r4=gU8-!WpDo23d zlXP>!<6js~sO2OmMrew+8o`fX0^6cDktGN1UiaYg^nUS%3Xp7;Mlgi2D!~wmGIvT? z=1-sZz)#f#16}-ebLi3AWS@-qKyUV^p(a;{_)+{eV&wfG9Z$@Y86{9!iS0gZBq%}Y zi0GU^zU7B&Hde}?(Y@U2{BwCpbm?T85J1P!el#}v{Q0&hb0=V;h^AgxSbF5$;r-NIOZ(Ia4bRdlmrh9xb2q(WA<@Vk%s;I!!eY&3!6zWc_y@OKi{qX()L;Vv){ zyVEaYd=oXTB=@~b<}gXMeme1=d1{?{!2iMLK8k0)I1-g@L`sjK7SGB8QoOKl^TaYH ziwWZst8@-VT0%*0O1#p$K@WV+kXvJ=BL*n^8k6qDX+0#@1b*L89i?~tY&#~P2JAY2 z0K#jt=~&ikPV6l6fr7)hZcc(g5J+)_VOwC4d*rPq!%sBY2D}UepX|(gR8#=~p#0Ysw1>H6Ct)!@58F_LL-K z1JPK~Z;yr5BCN()59Wg)%0EXxWMl^eBAs)b+Px?>1cUjURKip2DlyVRQ|`rnCsY&Zb`A)LCqkZ7Lst>ztZN=e}A*zG2# zR3Waajzwh1+FgU@O2KR>U;q9UNF&gJWr;mHHiuS{7CS^XE#`0(BaMyeXS4&1p6nBJ z$K51DbK^DP=Q#7&X;;8cLJD_OMfCK3!J?|2r{BJYWWb55oO0u zG}o^i#EeiYe#4;yk=_xWi6}0C2=mv{JU6JS0=hmT2w)SngRMp!lINW5;HR&aDpDK- zKwF;v!ZX{0VqZrIp!b5JcqouFhW%H1n!vx8CkQI_ay{lhcL9 zs9J6<9T@JQ2Z%Ek5gWO2DB|cO13)Kworu}i;#;D{Ke|E8afMr-L*P6EDwRgtU&6v> z)hhq%@#^HH6oUr60g)3~HKvI3r7!#=jn~%G09{t9*_#H=D@*y$CMzWx_Z3$Xj1%Iu zJXEzNGsmBKER&}EB!ekQmx%l?!W&N)HtP+!_XRtHPvjY*EN5C)1m51C$>zkb;5?YB zh64Y>kB3I9GCwKxeb?*0e0@N1=Aw99g#jAUqz-n(MAxx~(veQ252lSfj`J8;dP3gl z;ilgFR@u0NyHo2Z=xbw}d{PDRo>Q?+{=~GPcQ)h1?O$X@&-aweJeajKK?JrmlS6vV zUAyw|BbmfPN2`6U4fzp&kBtxUAfHOolDQ?C`ab9cHInfAuH>*{iuubpwJeNXQRw0M zlm!aoTq$$luVXH>fA zV)!+n>%*On)~PG)pYu&22+5{J4w-(U5aWY++j0XL(%mE%{?N5pvlhpX&GYFmVCE4Q z96YLaj*l^jp_BMAoIy}E#m|N6d4%daZ}I|TQk{xVemJ-LHU}&Iey5UhkvWMuJe@h3 zT!V;ELpuj!H|=%_TpgQS-hPKuNx5F4b}~sL7tu**t;|Pe$&T}+!CZ>$bAEGO(sDqe zj&OOO?aEk$54<{Vg(L%>it zHBtGg)}Q8pnYCu&Ka!bc7eX;BZq+x;kkw!;7w!os|>*Y4AC%nVJ zrRD8JTGgxhLR5U`;$9AY5(wP0$Y`bkMj!1Kv*o?%j1X36s?G$b0FJu06GLzn@Cz0cgGzf@g=cDriKm$Qz{i;0{T}O;6wgSb<7Z|UF#q5IPksp%8XUU$`)XR~qE7OXKW_(!{7POm ztx1*}Ak(PV-?8YA%T5)wq6zh9&o!Tpq4YvJSur~|j>GKAbQK(t;b7L5Y!lnspPIgS z_8IB_oBB8+xlx~&y#$mLNVFu!?^(USI#vtqRz5W0XwE-%)80)zr$h_R^NhH$D9x|F zYKR(dKoRLIQ|HCf;&rp45~(c6eq=|F{}*%8oZENeEY@iBxkp?>S*}s0U06Y*zai!% z7=#Fskxf<6q$~04S# z1^2>m$r=a5UbMS%`33ZeGD@F3ukXw@3d))q9xUtj7v!8pslcVC=q(K_GmRC}vQ@$# z`vG(!#mX>Um$eEnOY#u2!?W3HWgX+z-aydjgM0A7k{OM8_7kNehuv@7bnOCQW|1$R zN}UIqHWg05OC_q)FO*c_G@bfqMmp1Z8B|w%%mVMH$vSrN-1V_M1NHHxG)HzUmv{7& zc=BDGgQpmKdxLDu_NbR#in={3DwN@Eeyv>c?dTE?X038)B_N#H;AeKL=4A(u$Uk$2 zFEg%TxMwx|tWflZCa^Heswyzga#$`tcCP+z^f-L!Aahi>@*i8tdOFtlA(>#kOXbR2-`w^ zVrgUyNlyU<8#DSU=fQx=4*w3A&Op>^REk6uoS&z5Soz{@OL&BstUlGDb+Wy}48&^g zm8*}tU*pPiZsJ{ieeqj1Ixxm8hI-V%w zK;Js?$5D2{95r%JOL^EXsS^?fb-b^dPz|{42}=3E1K@a2^GIEc=yU+fpoex@0@;6n zdxBJUxKbNk>RY0-2MAEQYKf;*f{&{PEQVIg@5Ur$D;@y1t$Q^Roh4D1X84W^Vs^mq zTvII-bwHMUmeJFXs!`BAyNiI3(olXdc^n3eabzy%{eL+7)$~LqSGNYEE_8mos+N)o zYpM#wt-F5Zi^2#)8M;c~_I>m5)CY%;9?sQ2G=)dDNaKJSI4{$9x4@8iL}XI=+I$i@ zsX0j)$bXcTDQqABS5LKWEfOk%HQLn16~a*y#F~9Hv>!;JhU;ofi4_*Um*;G)TB_Wh z;@mUzkFCn>DC+QhmA(tRTV>h9@c_gKP&(cmr2=uTIS2Lg2d2Wj=>xU|jY55i1L$qe z&yP^PO3C}V4Qq*ZL>bCu1~x?j9<<7HdW|s;Rbg5hqDjL;hlN%AzIs5)d25fSw%P>B z^OS(}^Ya73e?nk7=}WiCMwGMc)ZKzjXzuF)*2fGU%HX+(cHJACa&CvJrn$m~Idugj zpoo#pR^275=-2sFO>4Zv2B_uAK5L}effgEIspFTRgj8`}2LrFj1N8{}59FKhvTof(m4a)y1}fJ; zLn(40{QN^L>;<{ItfJrn{~G@vim4rGK;E9e8_oA^@)CD4$p3X|QCHWG~_vKk_58e_EM8W~8xE5(F2=Q{|LVk;T3UQRLP^ zSX8>v&!QyKzveUO=c;kACtD0r_(@ez{T}m7>4@9YicO7&L%$DWG_AwfeNE;tLU)vr|mC$O$k-I_Rztig^?eaVl z+^?AFTeS+`St3Fav`V%KS*yLUz?&>qB;BM9p1SfunA)({Z|ni8$m?#LkH5xj9{KSt z#ti4mA~C-jaDL;1gd$4dKEPA+WkJxD1*v1Y-gU16momHLp~$5ZKsavVfd$sSI25iT z%B6ROjnT>7vFt6+5Z9X`7!d<=@>g=O4lBcrkpI)E#>oO(@7!D4Oa)<|t>taE;Eb^O zPp^V}(XcvpoJJ>!S|`d7kDeveNtZ?n_a=SCbMLXV^fhD!vFI@sGM`w|L(Md6(j9>B zD3SuN8LDaWF}n4{PI+L12;p~UsFoe?;<_IA(a^KZ9rQ=Y8jJv=O1ErJweKlfH%$#K zVYg9I(20%{D|*x-;1^>sYFA3*6$aQ_*j(ri|5A7TZKOM*OKst3G997Y z>!7lM%CJSY{sd?H{5axnNh~%&0I#@;ruKi0)+=w7uuqOKXmaiAvS5@b#-p$mwFk^RCY` zy-RWpz5Jm3^=P2G>;qTGrd!~xCM0zD+*9OAW{ZBXJhDHwL8{rvl$F}2&6Gs9DcFKR z&dX3z!B^8d_N|78TCYN0Lbf1rpJXNjv0m0`S5>MBzU7Yjz>gS2?QQC#g)Kx(#@LM&PnB~I#B zs5QBCOF1>1ane4}iaAjlOlmCTz?*}EUKw(8sPqmB=|qvS%Z0qN6C%K%s^I%6#c;It zTUjA-ud;s9s$Ze4PYar=u#=e+ckDB$_SOBPIl6SFG`Qg``+^u7RJv^28AE=KIhBFT zUyb`tM%*3z`SAoS%!wN2o=9X!EZD|)qXS7fLEb@_B#Bm0Z7FBP`}f2$1XyZ_or&vSc$xTHCu{NhAWY1OWzcSvhx>ZE+@suF z_bEz(X@jl=^M-l0B!K+3glf9UtfdnL2w@r$w~2dMQ07?$ox@`V+jXl%c_ z%i>K_lW^t;0Cvl9nCT4Z2! z6DLCd4sm#9Q2QG>Vu5W?&e^q$axPRG^X0WS6c3}Pj~ytIr;ZmyI(-`u=at=^X$4na za~e;q+SBl#`>9yDf&Czb_=Bo>C3<+I?t`_fFp>jS~}6 zT1D4Vlm>cYhz$U+?}H{>G3$BIsX|>mhposTWO+thpL1e|9Nfy?zQ^{hlsi6MDvm3n z2R<8i;(VLanAr6y)dNWU+rO|@J?y-~YrivKr(>L0M3t8VDBNxL?`lZh=vdbvct>~= z#tzx>H!_Gp>}mn^CIwaTFz$JSt~vDGZ&rg53`SEWR`E-Sxv`R!kN-XER(&n`TS4ow>bA#SxR5wt7AO=|790~2APnu@ZUO}#2&fz4JD8(nxRpu0M2>Vk?Y6Ly zL?d@<_`>C`Yy&W17iOdEU98keROGMwnsHz`9b+eRra`;RE;;+OzOs4mr#U+hk7>U*b$vJ&r`K8&i}G`j=Zyi7VEx4o+Fn4PIR#+& zKB*f&A~=H#(S#9CAme`Nbz}I3lm~G2j0oNwdL z(&CsTN-di`gT&98aiq{cmCYhLZ$q{z967vWWz;w*@b*>FDckE5ImFEU)%vSG$pUT8 zJ;GI>mWgl$;RdG4*OQ)he=&7CI8?AUxs1bEEaX=M56CY6+H0;I9Ss|v`DXV|K&kMM zMkf@6BiEqeN*ARiGx9WBC!jcsK5UC|DCh{fz5S&V>3ecfK?QxhZsa&x9gENf8I>9m zCC`r~30umCz@&kf+bGkta4#I?g&PuGmI)i)7Es&hIOShHTnN+HdNbjRa_@KD5Otc# zK2U%St?*V)nUVsfB86f?h`5bx^TD&O)U%C3WGCY^YN(3}E7(V%3qaAqTKmRc#&s#s z&uDT4Wdcqpmcd~g6e7ILdGyX`LL_18W5B=F6vfm8w8EFRYme4C0Pehom`_S;qQkwV z5%eg;TT=z|;s-Szd8s9JkGmL$S~5KdQF| zcD)vn|1TfhzmDf*PmzQ5+U8w*G{+KIho1hqz>^sm9=SURnem&9431gA(p2U;UYb<$ zQ(BWHby%aERfQf>SGs`%sn^2vNTS2YQo>6ug40hd0@}Iu1|V`bhhLw}^{<^J8rR{6 z88zEIR@KfDM7B1w@ivPH8*>&h<-HIyGHFUiOstYJ=sA$H&G!;eGfH6pzy+M(Hij)qT?! zVp2Peomq?_f!L7=ssbWG8F@z$@i%_CdS_x4HdeA|&%f_zb>9q+j9HLqq5GYV|C@rC zkI(`~4?~Wwd1M%c0o$C6v>adTg=;XD<{<&|_VFh->9@jHAEG6a2}sYv$e=Oa0#!Vh_uT98t^rZP+Ot>^*BtqZ+7qMBqw zk*mVw??@TYuu`jJZ%IjlkD4WPb~*u@`@@pxLNt#0_o@=n?X0(jH-^aTytr<14K4c7aci;mCXnYJpkv~yRDCcnp*!4$o#CGYrq+}@1+t&W2;ddK z#5WwLDIhA4gbb!*5nbmo><7SCLJX3@SE`eyt_}7=$|20C~rhFa`FH>Q{nwX z%r#uj`Q#-N$|OG~s#AF^W7fb0GwT8}FaegoJrA|OBS=qA#|o%4sV!k*`cd`$%ney7 zqj#|m*cbP7bZtEE47ZDcl1VhPi#7Ynu;c(#edppmn8F9cZGUHGhc4NRa=)sM;k54( zd)8^ng`?DEz60LA62Wct{iA^5Io=}u9BQOr<=Or9>?P}eE}h=+9;!O9hsJ47Uw-Ep z(jU}e1Jmngy6-P#Q?E0no^>$1{txAA!&oq8cxiSC$utU8MSUHJl~X+Nd+|0bILUnx z+HFYTC_T2vR4(gz^SP2pP7RZDSPe;75xzGkEz2fC(a8hP<)?Tw8ML!s5DjV446l_h z?W#vWA5C`myx4=>>2?sd!iD$X^(U9v_0}wnURfw;7UmC5*iUI0ZE4P36AA42xY9Bl zUX+2T`08hD5!JeQG^&M(D(Bm-NOoMZqxaMOC>u<3L4hhck~f3j5Joan&dBr{Sk5`s zPVz!d*1P_fJE8EI=NxeOa5>(9PNtm%qfo5&y zn2bwezGLyK9WUq=SZnz7Mqk?h95Zx~U(LJ#xs{(AiItMH@4=9GaE~)$If6sOqsMg57`Pu)4i;Sk)SU`77U~3}`3C#buF$QY7s3(`~s>wq8&4;AZAHCtd?WO=fS-{C+ z;cT8%XXGF$p;IJIxx|0#OFPQd^ONeiTa^94l`Yy2SpCywpjf?1%t3*q9wPo)k0CS^ ztr0GEK^nj!#27 z#5Zo3D`iKW)uCI49I90Z#rRZQ{!H7Z`d>xEg*-Z>ybD?t{v-(&Vi`a)NWuk(oSmI* zWu~ymDecX2H19myKqV;OVZ1@5{rVb27Aj01?lUhmpYS};c{?$5{Af{9;U$nLb@~^0m=X-FDp{C`+O#RM)@-DCnPe;h9d>Uk=+zgi|-txls zV()3(4%RbX9PwmyB9toF48#3k>p$#`4{Xn=M%{45<(>;)4-iK)T2mY941&U5$LdV7 zU1C5XhgQ$ZKV^0ni;F_yXdNf$Pxh`kjFS+GZD_5xI|IrPO~1pY*^h;8;sGI^H-dfd zE?ovwWabjfp$Ckom;gn4u4RSVz+-yFa#X6?3_AXuR49+5obC#$jh9>S$o2sKY_+H6 zo}}W}1A4rm*VfCzqoY=u*uGjH6bhfpLl*>|C2~X-F`kj$GPpI?2uB(6d~k_EAQvH< zB(GsjB;B+6-CQlihQJ?roTP)mvY)t7M#U2ooW`3HSMAhqq|DwcB_Oh>cqEP@r;*(& z)+SMwfDZ}}tuy#9^5*pXc&4-*4o@3vkBBT+0S82a%~mgW3^iY8i*RMQuZH{+CpJ9-oF?LuqvUg z9%AqLjq;8|JfCk&%48sEW(MKzS`u;Ph-eh490{ie8=xh9=WH{Bw3mykpWz@&&Z@xX zh^47pH@QrE?Fe++$zOTmKrW3{1;Sj>m#bAE$(I@gfE>mWCdmy4?Uyepq=zJdRM}#l zjTflEdFgf`W4a+OEk^DEGqhqjB&`uASHx@3eMlEr`b!_vTaKt586Bpg*8Z17htbdI zfIYjOTxc{s+Fix`y_5=EmtCjfo0?;}L@Evdn`8e}^@fOz*7ksk0TW+6eU~znEZ&iY z@McH$2rFr9c|+|1a4P#N_G?;v=I5UaO}}Kx`I#=Vv<-($L8%Qs4dr`he-W<;!AFK^ z_*UC!f+Ti>Rjf}kUI*>`?AVFrQaT!-oLFxx6WYa(qIYejUQkY{b|Yq}$dDX3JVr`b z*95cA@tuIqfG}_=Tnlnpz3;`7Sho12|M|h7I->UCl$a(hc#2(>o}Z!)aYZzlKAKZV zx{G>55-KPhBiqweiXtR+x*n30c;>?KJH1wEfS4)^?~jcZ9%Rbe4;4W2^UR$2xzfVi zw=tBgZD?>hq~sUvtQLr4=6l=x{lMRi-0>z4R@2;m9$X!`bT!RF)^_3<*~>M(2#my` zseN%-mgPxiI}#p->7_^+ej~@PhcKx#!pV*JnX3W8yYji*FXQ#>Lsfy#dKDPOnKl;T z!ND!Ir|;}eeIKA&)@B^d4vLp_jN^67DvTs5sno*{UD!h9y2$~$Hn$#qZd+@OUnf*XV*%hS7GdXJaWJG@tW?9FyNgyB@ zrzeLPV}I{XyC)2fp6CmG>7XI|XSD)#Ejpauyza%$NEgXd^PU-Z9;WrF@1(zqPDhU? zm_`;X2-ePtvHsy6BC9rHL+G9|;bfCgC%u@=(A0t}3t3A!=+WFl7z}d&M>03T7pS3l?)S_KKT=Kt`&u4MI7SOj>z= z(S;ZZI(Nb-&E;RIh$z0oi$~le>2MGK@i!U`zezYC+Wcams;^IP>ui&}Xn$pbP!wkm zg!CaP@u~({%}V*RXaI9Qq-G5cO3!`Yr1gdGHQ8Qp$3qeI?P_Mob$VC^JE83Iq61v7 z*>xP^;P`Lj^cpH)^6rOCF7o+x6vse!K`?}t3$EiyV&hP9--|@~et;kSe#Mmz4>sZA z+~*l>`|FOx;*=~pvkuPQZra}X!r8sZ2I;x-IF0w5%x8!Fwg!X~5Wm`(Jw_F#^7n5e zg0k1d&;Qnu@tzI#ZuM(VE8i|uDo6cet}WSaXNCy8C#{)`RNt>i4D*9->4z-(h+2Ku zlZlE9jb$g3LpMci0$`w_^AK7ED-5m`tAtcfJ}uGr`$C%v*EO{MuJ1oYS-R4JZ9^2i zdmH#4L&Ea8T%lNdQyebYsYYNRITj%FD_6)hU}gbY6~0M3W+S8E7DS~xgt275Z{AHK=eOqcg`p2&ohj) z4b2-Jif|m+b|t1DzwsS3=)YanBjJwa8wd-8{u9)ARqSH$>9wZO72~~3n+X&S{4H!5 zeJ4lkhBe;8aq)o;m>BhdPkXQGg+;{EH;vK$U-Ck1f8$b)9*Mzs@(Wf?=UD(UBqTy;W57tLA<$NC>x2btu#)aQ*kCm8JrzDwI$reh z{?Zj#Qcn9mE!KU`c{S|v49!$;5~*B6yvg7|I~FisihRqF%Ru!C$PD?KTzTD&FVRXR zAVV{zD_#4~Hth4athlTy<*Rf5F%^nPKS(VyZ_>tOQk9$x=X zs1txH8&XpOT`adQit1KePSH`5ztX^+?5zHA&jJ0dpbdcYu=*rtv-TSwenY%~6AU91 zeV0~#s(Z*0JUJ)_W?pL++Nv52ypH)s6f60#vtmxv`(soQTPf`=_draFSBKm*W*)!j z*L@`amcmUNBi~3~ywR&s^R4p;430gzK7Mbf@LWj=T-* zH_+xQ13B7DJdbd|5Uxj-LNoP)MI}|zch=M$)vx@?T}T7`#T(0sg?ph0%zIrHlvCgd_nS~XcIG$)jJY0^ z_=cAQfC8jWg#!!&1e2JWC$YAKTQ>9Dc=f}02<9$u+dlLS7a|&?xO;mz0xUIQucwxs zaC(`=gMu36(g~z8^$G{roMNPPX7&v;qGyvlJb544Llju+YvuL0^OLWB4&1VBPS=GL zSO8)E(+6=$#76_YuQ$-^Q^HOgN!Gsx!(zxpdGj6ai}U%|OW5p|OSYwa_3Z7GSpPiS z9l^|IQ6WYkxGMSDw1`7o{~5aAN+3}|5OM4_DcNTIicfW<8=)Xy*7}`0)3LTZdc}Wp z3?!qfL@Oc_B3rQoQI{1X5;kU8wdIEQzIW%bO!j`h_JRbF8e`jlEwL3GRZolHxB59iQmD#L#kuq2q@}^LL9MxK zSv&VHNOH1>3GNZH!78k0*MQAH#c0_f*Jm@ORfZx48Aetl{P7>W&9Re>j8H2Rq zX{MqWZR?EC=r64&1ng}&Jrom3uz0))Pxv5Ft`YEJrZK2pt-G~FI|uj=UUVvRPC?wX%& z%VB!ksBq7;57TX52o%vYua!sgZ6^Ve4fP@B5p<2eA_W#_u@tnVh*j?;5T=F-XQLRJfV+y)nK zGU*2+)&QxhDm&DXP_WeUnA`pXpnW35J2_PjbXO(eSGh_fI*4tLP5l7ctn)40NTcjZ zP0K9fYkE+#H`B(FpgLZWF$Em1hl=%+5p!xH!)lSXOy+0##z9VWE+RHj)M;w^sn+r1 zS~^L_yr0VlI$*IRi1$vJ^2+wHqKOD%-b85(@+X{p1SRVuGxy*m2v{BTVNEwdf(G_R%{YX{n zRb@nUjFs}v)`Lm`8F8qyb{B}7%YVRnRs{xQuZ%OH`s(0OFvuFSqv~KhM?-!?cr2n0 zt-b z(rjh(fjOdMRAX^}kSC>GiHw8J1>rE~CN8QN)LcF8ZF3Zfsq3y~|^P{K;0EgC=C>TI{Slu5C+ z6lmY{R7zS5+Ed92rPm@~2dJ$T_2>B~ID8+b%fC2^NU3M_vyn_ImHxcHF4fy)Pj~5E zmxon9K0VAVTx3$T(J6Y)_kg8i>=c0z6WeQ&P}+@DTgBdn7JM_PZmEXL^N2vJX15C_ zDCM#q8mYUwij^QvVpwL{-6LjL*`P!@gv*5e)`>J-KYm{>D6W_dK#3RXT}_pO%FH4$ zzC*Py$%>#Z2Jixz^9Ekk31AKwFDG*s7k_AR3u|)x&;JI4PVG7lyBIrZdkvJQB50=x zA^zlGKnIA3?2SraVXZOZBC3NUhk#Qfl2k6Xty%UCDi6Km2m%kVzu!Q%<5i(IJ0r0J zMz0kLptXx<`q#bippqw!2+o_uN+yjcq3HQ@`_|MM&)7ewB@+G%~Tm;AiMo{cY?4Xyq1S$g# z`5S~gf#B3t3Y|8T;SM^MdHP!J&NjPD=wSZjqihI4gV@X8eLC-Z)o_E&y1{qPHa=zI zHdGivaG`N^mp(0dPgrWaVAkvUC)w*ni`E$hm}ucGEe!a{vN~Sn>0)y#=d0Q}=c&CJ4pj{J^@2x{xAWL;;AY z7Z4o%dIJhAn*YZyMSvQJ59@qVO|1XukWBTy2oPip64@iW z&$8ilw!Rkle3$%E)HN@hvaBo7W%n3m%~OvdEJ~-lNOFPrgY`ci$VaKJa*nVB&QJ{{?2YP<`PRg77yHD!JG;hR&I$aY4@s;SM_EY0#!|DG46#VF`dXy>`nHVcM8st5XS?-YXhQ`-Mm{j!3gW@meOX?kD4ghCWds$Vw zzW1xS9^}o*Q5Up5ZBxuIb535`8-DZ%gB3SV_CaX`l(jC6lHftke8FkM#T@AY88EA5 zoI7{2K~A52BZg5MyXn*`_xDB$r2^g-aID!$vzr22{pXk{{b$k^)qjOSQf-PQ*;aQe z4%E|5Z!zI%?V~HHI%-t_d2< zY{xmz#$lq4Ow~{XOjsJ3Zi#-PGPzH2UF$H-1|AIEPdE)5$-*7Gj%b|tscuD}Ap4_d z^dmAXY1v;|O(2ye*D!{qhA2V~r2IpJrE^~0_^Wev8yA|dW)^=pC7W#-nnm`yBo}6} z(5N!^ zD!1-nNiG8cUr#G3&BUn;tZd~2hs*4c>e^{Z(jz91E}%*(%J#{m?(!0G%aLAKyAQ*S zsNQ{bSo;TO`3CuRwQ@8eP+@y*Z$GfuJf8w4Bl^a`N6`7{L@I})BZRo-=$$5>ciTSj zdftoj!LSq;cW8V|FX7Y*hEVlAveBxu`Uwk}4439l6(IfN^RA_akX_X>L0~hW6o2+t z+FPWqeu}fA*>ozk$i?~yH69a;Lo)~qYK4|<`sn8vr9%+36Jz7x17#ZN(^!iu_?c+6Clu5 z)HCIi1p}RNV5b%ByAwbp!Xab}9HVyN%#H-MC-%0b9UchmP6 z@+_Kg3*Rx$psTqlrV1C{|(scpg2981n;N~#kDtJPvUg@^3`|jo)FA7J=ViO zJ5+Z|mm~NSBF`coX7D|-yLWOU5d$OFGY^@GtvJo$gaPn0#BNAYJATEu_|CI{*SgqSLcG+Fw zB2VoBsx^4w_@!Y}el+yJDjZ-4t+dAqezvHuh&vU^fbUfWA_5wBwXq8$~pVv0u1;L&k_Jq5YnfNzxrgB5gS;v;*Zx zN+i}g7g&hk@(8Xo;}=Jx0!&^HSnXn**GqtnkAZq?ke;EkOtNJ&Ei~&D=_Z1FOZPXq zp)})>Yw=zsd{WW}sgl%*4gL5W>*U8asG7{i=-!$Xs4rnl(4?#4h1NIT36#d=q2~Eb z$Y0o%u{~^2C%uAX*(oX>4EA|{AfRpx#z+Ry*Z&%SPBv;>Y<4K57ptR!n2Z_bh&wS- z`q(oFwsC?>fHM;GPd%C#^?7WL;nD&gE%GY*2kwb~`($q*2sf$1#`ie{`W!6hvVPG{<*Sox@0QNzyvA{MPj40F@}GEler9HJBU6h78Pu=rEwHH zcVXlks6e$;jEf{yY>I5B+uM@l^z|fkB-fjpzQR=*yC@T%zCmF*8|TY?VZS2Y$fD5i zaooIkrtPIh6tk`X`TK~;bJ+%<8N``hIC0|6uM083m^VJ>&Bj_pe;X`u$lvcr85@KR z(yLw_$uE?Y7!^^GG(ZEvRH_!#DTRmL1|zI=Wx?DkGq86pkS_AS(zF}#Fd;?(*Gtwe z>~qtJisV6GZTSo?{A`S}^&nf}GukT>%@s6TR}O%*(^w|OzjDOu#io$KS(g6-#fV?V ziCZgFkge$^zn%Fv(91=Z{1m$QOuAMe!U)M9Qc1oIjycWA%qou#0mUqOHe#8Ne1e}b zmIoj0yP@=Z@2>hSt46vpp>f~(yGh|^m!vG!#2wr5ru+<~gf|9@aIlCwUz&Y)E&`&o zBX$VGxWmK0G9OMY#;))c6FS;^(4G(s2DyP}i9oK4i2SxXf5%ieI3dQvapSg<4+72e zh)B2!iIOI>#Kf$sU_+*4GL11Nc%-Fq3D4gHx7~Hh(=R)}RSaFGB&0s>JSnzt$kJ+F#0nCPpD(^l2P=-^1rSATdT{~CiIr_nPBE`$L)B$i*3NbRW};j*^73xyuy--D zoE;lsi>2g@WkMOwsBnYJLlOLoSpgM#eBtU#Z7o8aZ7mDYT@_J^3LpnITP&=e6#Yz> zt-w-055=hl35u>*G3RwLLl2<%uc9%pGJ?HkQ}cGsAJt_o8J*kMeBonraT~I6HIFwg zwLVj2D(L2d?Miik;3nhKdqa1}XJy0Dh|ql-!fEH929|{LugpkZ+01d1q~q1Jr`XQF)HOD#Eu^wJuIh+Z~ynoKeZPmoLS=OUKddNsE#F z7J*AOMu;UVcTo~O+t`H4YCr(I9T{)Gu`d7=F`vK{T(NKje8;l?_7qslLnQ(EE+gQy zv#03(kl-R51KO#u-5STE3p!tC;3X^o6p=je%E+XDBe8zi!0uE+BN{lxa)-L45qVi= z*zC~A3@MQo10X>pX*V+G=1RoHt;^~bzjwyXck_lL;N&uzD=Fe3Yo<4QgeCc|9Ywb8 z#M6hGQ-R;mqopf|p`?}1KTXXVF}nb}bQ3vzzE8%^D>QC_4=|u=L!b*Ou^h8dVt7dln$ z(qRCz69o<9cRmiLs=W*qYWs>1Kn>d0Pf-M?rW9tgIuQ`cDU>?(DMPgY;w zD=XhzcFGN(Z9x|+pc^5{k9)5UWungp0QepThfLLF$AI0P_pE6wyNQfrM(A1^$b|Ep z#Qw>`Rvi$gG<^=Z)dlN0a3C|ybCpO%>Ac54q@GE>|JXO{hOzs#<9y;bo!8dXH$`t+ zceBY5!gw94E)xM3P9fAiP1b&>n6-EIxBte_yB~21x7gr&vg^R(K%&ayxmkhynO-(SzEFwE zmkC(|tN3-Jm0_uCdYLJ1)cXy|yY4c;R#zhaEuyH2$(4qXPRQmfn;+q0-+Ld@+`_`m zp56L8idM0y5)^VB#G7ra{{Yh3<<{e@Q#b0CUEB3f6Ec9JdIp^@1qbQTGEwK7N2{gl zyF_*PJ~4-vCMAzB$HH6V?e?D+*kgIrH#`UZuuT}Ts+nF7^CkJmsb)3aFq1&>?0Y4l zxG<-`Q_bac_AD`OtsL0AdDDAF;}2F%Jg_n!?k6ObWoDb$Z;qcgPp^`IrR`J!D8F8Yy%Ic*3#d1DUyZ#Zcr&5oiu4)=KL8 zGo9gk#BW+eSPzN|3PI_k6D#z@ZndCi!9xQ3rbFKAcD9nMc|FdYXzz&H@76q z%700*nN~sPuh&b>7xpy#ZhNFzNYLYLw9rb7Ow}XHowx8G9;2^2` zoMj1>R#IMVhcDNA@VU+9;lq*ck5_kG==3Ug+~9k=k2azWE8HJ_3l^F@JQoe)$evsg zjV()3YarVqTGgmW^tHx9zijp|WN!>~ai!FCi^( z%#p{KBe+5C-0P9AQwnZTD=4B?7->M#EfwJbn=qnkb;oIdzvVazFcE`2X#U9n@)p}j zcBTgHb}40RU(p0mo|?d$AD%glP+;BrY}u^Zs3s3bU|i$r71k0%4maI0>eRWJ(lnpx zu9l>_02odSxd-O;oD#lrrj!D^7T@jd13hpKBO<(m$uCI2{8-SArlY9T_6v@wlJxfL z;Tvx5(2V8aNd`=zw6Aq%8KjJW%??*Ri{V=@pvjY@BP)98+#iMemTy5Q!q4hk6rs(+#BpV-(ramI?Vnfy_~gBO}#C>)XOOqm`21fElh)3Ap?Qh;1n>@riHP5 z0FW)IqDfB%%cm*!@7#^~Y*f_nNCvsx#8971m$)tYM*nX#S&$vS$5Jgw*3>44I5-DSYtbsw}eGc!?~ zDK5%$G&k%>d@82aNi^%V`7umE*eke4EhIlUzNY?+qCUHXb|_H+Efq9m)hq8Q7S^3K z5vXJaZ>&+qvsf>Z)~#~8l@KPX5u2wGyR6dl)51jD)Je~@V5n+c^UU~Uuz z7&Q^Plt;2!fprNkL?b&&6Ejts()r|nGDxU6awSW1_j5A}uFh;hJ&eqDFn7YGC0t@1 zb1M>a+CM{&zWwJ(4X|xPq2ya#O;kVP23c9hSFEeD=lyJ5xJ{ z;Iz5!SD)A>7P) zO7_(lfD@N)yK|>veE~a#OzCX!_#X5+t=s==P5??eP|`hYzA~m%B2EAGJ!4XTtXh{Sr%DcFeFt+&j6W*K+}Q`Kpg(}L(#N&QnUz!! zAv~GYR{jW*D+i`n)5Xn~N;g!Z6#B`f7ws|r@yZ&Ft{(m!d#@9rX1tTX5e9;Lz+wjO zT!_~tYhyvBS_AWVI8h{|3O=)qgJd*@sSVxOU06t7<n``n7`A2|f0X>)A7)C)mI8oS= z59~!enbQFEI=n5rCaRrfa0+S09_l*K zeCue_rqX5b2Nznm=(>f&RtZl`x|Fwgx{5Jt#Ggqr8#8>j#&2Wpf{AV5RV`YZa>47Z zR##>$YxtbRvxDO2-gDJ}f^8Fz$zMRz7aNE6S2JW)tlFpl4NR|!DA`s2Ff%{Znzjj% zq-_$X2V<7KZwmA^;ia$-nbX&pY%i77eq|rMNbCkK*AhPuwho43_;z{6rBL1HDILbx zY6zujMuJ|z+Uv0AjP<_61>&O@JS--eB7XlEq8T2E~oshaH?5XinOCs0%;FuE~$%JCbx6;rR0vkK!`4 zT%A6ZHL#$9efTm~Ik)_bNbu z7mAFwMq8UYY9H10-ccod3tGJw*uS7ZjC6X0zvpe30y1G@RXc*ug&Rj=1*B&EycnZBY+txa7{#L7I%ufMm8 zn)4<6U=%ETt`-o&SbNtC9!h%!q{1~`a6_4-nDnxsR|&F-%1YdC@Ay&ZA-N3JhZ-2E z)&5hK!WMy@jN05d$LS69UpjuG#p6q1;XLt0pcNIs!ey1f+#I}Glk4SSP4TNkO|2}) z;^lqi*=dK47R(?sKrG9Y|z0i=*cO3tT zk}Jc-0+|~56`K*&*|S~m#Csg2DX3RMS!)f{obU3JBPufBt)A^kx6X@A%UX$ z&#TP3+IJ9|{BiYs`+na86Wh~2#EfvOS<;FJpoT;mMOa6&dKPyNqVf_sgkq6sk6x;Q z9=t}5+Y=x|R3bY(H9lGg-hMzdQ1Vd1El+n9UyX4+Fpwf%UhBZ!&+^j>9iahxwW3sN zjX?8=Fc$jbND^G;%v!8zxTV64jIts7wV)l7*K}_yej~V2%Qh1DHsP;$Q)-9rQV}CwJ^`H5nvDkf%@3REc+ViVA^k5NUcnJB)Dw0Na?G`E+ry;U&~_UYzcWZIw*U0l0J_)jO0Cia8N zd&0|m27JoTtY{#V@VO%qiAiIuw z@qQ87normlh(@?xXE()Rub68Ba3`){KsID?7cR zhnr_BQyw%aAzdtQJt8-aAnSnGPuX??!jZKGc z;{@_*|M>|W`kWE)(!2xp;qk4&uzxvF-+y^Hohy0+dN^fE!CJ)OwdBeJ7&Kpr_E$es zr>hV+UVE-~I>cTzoqHI~iac?R+L=%N>$vXhi~y85%QZgz7a!-( ze&|L_we9jA2guYH1-Cz5TGFsYqmCD(P#w><6fJI(L;%0J1u}HrS6z^ja(QpBW|fMM zjVcEBpwM$K!AuN3u-e^H4qe`N^xI_ZF{2ui%+Ny?&(8D|_9Gz#vqTse+FW)|n3nni zV)Ks$QgIZ-yskMM3qx0pc^DgduR@%f)|AOJa4@ueuQLrRxLQH74|u>zXIntyGUde??*ZZ|MUOk|sW+A_p`r#}J;s5iX_0>%ue8 zCWc1Zs^MB0%H^NJH2T2D;(P`$Uhr7rN;d>NzPk;5t&R|)5Hy(g+tPp90KVlfs2w8{JK5pa`Ut6JGH zd)jkTP4XmTz&@_hoTql;P@#Tzxd-X{`YK;uw5l)@s1a>Qv)$D200GsqGFhtoV%`l} zm0Eh{t%qR)239UI5ogNKUs5F1OC@4!FErE!McKu#7j2KU!C<_Z{_1r7=dr-IYW>BhPG%85JO%}_k#KePDJM{21hgiuF1M!;9>rtR2&9`

#j2QLNiEt z0R77IW9X?ljtYWA-5rs}+)|dBE=NvrS#d6I>4a8vS#SH5UckgeD@@)sb}o8m za;W&vc*f{|GSDWb_uXT*}73%txBA>W3Z8s}J~UAdKNm2?sK( z9UfIYokkUe<&2YM{g!<0a)`9@JzRnjWe|z`iv8m4PYm2=`kR46dVON@jqzB>_TCr) zc6)O>!ru1+ZR*}|>)KaB-*KBl}t z2+yA#Q*9)@a>M1b>(#<6YjUI{;o1}56RKCtr;KnACKi>RVzwH(KL=28TR~+nVpgbF z@UrqR7R(W={ohqQoZX?%vM9nyisY3tgeGB-&n|E3S2nH@QyiNcYgl!sB0m6OuT{{=>X<-k;HwMJ zeyROKkEe#yEvt8`M|3J0hmE-X2u*iAq>b>k zrmu!Xa(?5wI!%`zyY>h1rI!%I2-+d^Z7n7I&B8mj$*x82XEOl}ux`zaBuewf<=U)r zi&pE3Q_s&-l9A5(7+$ZFKffd8H`2ziL( z-V1*PQHp&Kv+b2ylne_evp|yYV2Xt)p!h9dDMC1O{fr)?*hoRNBxh%?6h4@pC-Om` zpp+&(@E`^v2RrhyQDy%l#M)1B=>gm8u~t#f>wG8qD!xg$*T!+MRw-}77Pnu2wwP>w-v(fpb^&BNe;6SGF^n=c4QE=fmoM`aE(W{W-eufXrRD$fkTT%BNJ_Ok! z|5c7g(#t?*$7uX;suUR@WH3QBh~qObeOVBkTH4hd#t#!p`6|S4#n5Tx(&S53GY~A5 zbtp#i;&P-gcZ?6x*&B&L#e;$rGkU7pG^vocjE&6dD60Ev)$GwL$K(dL9lI(w8>KX7 z{eN{SGY`OxH9b?1@YudQAq&>-9VU|M#-3+{+~4EtKVnIJ66JA z%TseNfwlXF=+gtQaowE6mJ3ePknc+ba*!!7^m|0=c=cFt?iiWtz2uOrn!5Vfx{ceo zPRNIJTNc~TgACZN75x+5DZodom@Bu)@Ik#(t>x6_^s3542y6eRl z-}~#5iK@Y5w9SS~Q8ErosJe%4SD%wUD0h#^WQ8T(BB_wAp2|v5A7M!MBi}mTE2agr z#!s4~1bsJOUgf=fqRs$hXBV-ldQSR8S3`hYP{(0e@?5ibgAL9+;XHp19E zLWpfOCu*(_mtkU37pV#&Tcpwhh&>8g-=jXrQ zixU4{L(hbBud%(>DD*t2QR_M0EyE<}HQYK{jnf{Z1+t{?N|IA6$mili~5Gt1}l*%>Y;)Za>`4uqg%>3zr~7%5Su(Yu_)6S`TS+^Nsv!wOGj* zVIvxrzP4UQ&EfxBJzKx)41x}a=$8L%1)~~^F8%uew+jg=B5JkBKUOq>-n%(6;6q{tdvpYq$l4;B_yM zGJwuqb*P_B9fd+6M*IPZE-)xVl!_d(v`lZ}nM_2MMZZ84M||R=Q``j=K$v_)`!sQ> z;xUvg!|=yL^PdAByyA1Alu2G#r zrvm^*&DSPfLkx)X!D#d`4>B$eTAlfWmZz1>O&}CIiRHxjrgFc_2g6qn+0TuFlF-N> zsS{zr15Yd@!M)WwaUpm<|Jm3DUEdfdeQq(v4F4UFMr|`5y>4_3u&l8vnd({R5;4BW zo|}S{J23i_!4ky4Lf7?|F*zzO%4E%|*8d2>6j$LwmAkO)VfL<{{P}zo3wGi5-^8D8 z+mdJyaO=z29vr2A&>G}c|F&)z{B4RH6u{q2o9oP9CrQ-vZAkzcG3j?w86CLfishPI zk6jZf$4E5O?7tpvtqUL;>baIuxOw2v{byztMP>kJGwi&kjfUvxxyYzsS#oxgW4--U z0h#`A5w-AyC$WB_7KY%v1bi4>rKzM7>x;nXBl9|H`$hv$drh&&<%;&}k)ag`In|uT zZOfXi#QTE#tf<%bRV4R?nxvGD98iZlH1Ilx)%oy; z#R$nC+Rj2|MS!%1ARi%Yjcccyz=ui8j=+_ys*NGL?jbc}_=GPLW-?FlImawNk-(&q zckGk_W3S{__C_Y6>&5IB}`&G!FlC*uPZi7q!QpC@kPQ8^3sx2_ILy67Hln1$sz^ou@252NG=7&u zH1)Xx0E1~509dW2;*BXOR^Q?($aeavHZFXm*xcxPzj?g~3DZs|6}n-K{zz(T#&G?N zdL+h4vji8eQD&tEOf-r=>}jpxdFzSKhFSXkJr7}%E@yCpWJzGF#^D`D8f!M=@SjFG zjZhtftpo>Ec=AGfCW7eeTKkB!DtbHjpo+QQCxfs!is(Z0Zjac^-e)T$%9M@wji&kf zx>#>Hu$0|_Z;!MK6`E?Q?|&g6VNplr znw17hY^FkPpZht%q77N%3J&UmGp0Nrv93Lee7Gm#uSwv_w9zNg>*8nn(W<#PZp@g} z&b(I>4iv1+)Y3AkrYmptD>?f|ACuY18eUSULw+vl2U?lH)1LKqi5m=8rAZIbuR@HNu zAsvsbpPge9^&VegnNo%lJt?fuXPc>0Tid~wF{D++--++JDPco(p^*y9^>UM9CoByn z9qd;uPoii8>3w;(u9+W4q8{Oe_I0WlZ6j%c<^hD7c>_ zl6dnisZ66_UP#sy5ZmOR75SL>ZUvAa|DhFnvi7ODDMCFTxoV_)#}o0)F0YO0cQ%{+ zVf~u3KqK$J!PONdBIt3A2eIMS3w>J)^>g68$&27iU~NItBS>LyR+|Ja30xZ6xejO? zYhycUS)yI5h)B;Patai1uqCQ7#pca;CtH{mRfHsspy4ESOk{S3CwD?r+f`5Z8y?j$Cq>`T~76#yz8;AEhFouFXHjTwxN8O9N_}XI`Bsc+|v|1pHl`y zj=#CHN2}=*2~bqP+=+{KwH260@qD$SV}#vSB|dW6 zn(Bmnxk4?8u`p@8Z}yUvEM8~)wMVJ5a~?K5sRGjJQj~?D6|}(sE@RmUVOqE z1*KgOiJY#fKPwjg&knB885iDA!iu#;FA(Kh#lI}ud9L4Uh{ssHnby3k%&V;oOt5LD z!Rh8xIS=m=TODuit7kBZD4hYZmFs+IH2}mnMhQ_c%l(4jILxe(mJ3@FJEH`wpo*tx z*ApX><(^c#kbm}73k$8CnhrWTJr6Xg_Kw>K;E3qVOPQfGpfyNUJsYZwD9(KnRi2wx zLB$tc85U0H{hO-^pI2{#XXl`RJ=jQ(2g7J#kwIA(TqY4&{iL%XVYky;-olq-IB^YU z_j^)QE868>)V^qo*!~lz4c|nBa~ss&Ua80B9Aa$2)Ej4W0qEY=-apSeSfApkju3S9xvXhf+ie%50vb z!pN>1Vs!IV*Om(g6zR8Mm*xw)@(@?o>_cJH?+kZ%o_14J71c`fN6BJogUxV z5DrR~-#7Y2CROqmdL;+4%hX;fxZn_i{7k`GhyZ=7vpRoW~RwH`#MQPO0h!?(XP`@}ZOen!k3TLE7lsg#^e$rGQ} zpA=g2Ij@{2`Z}ZlnkN3lIeBC{Cr>~Tj%}$sIT&cw!ma6}zq6wc{w?34#WbPX+;<$1 z8b@4>!T60Qy)nRm0xfS>X7*S;K;IC4X5?=V0N5=rvfm{aY}p&$Q`=Yk1g9r88cb3) zCcY3DTo%$<9(Wf{AiEd1^aSLVDPaxaLbmVT^?Sf)mI

PQXEa4@;QAwoi9RY+I5% z3XZEnC2oBF!ut+D3AIqlwG$M)C~Oov@xf^Esr{yuI!39-$rE_m9dX2G7I2$6I`e{n z$0#dfrslSXlxqb~y8VK$`}ilW%M&wa{_k~AML*p$XVuOrbv}(JE6S3>yKcr^ZTsGkj&H*keq0dd7gF00AX4=gOBx}l`imuzT@$G-0~M*{@AaO^TVQhP zW(uD`_t)}d^f+*YYo<%lmzEOf6Q;lO#!@k6>`thiL$5df67>HnV+j1IZSQ8joo0pq zh`#5(Xa9(QaIB~Ds7oQYWmrBgL*zr=qAz>AcAu0s5k_bW;VD+#zgg2ub znk`WeB^c-GaFg9dap?YdDGOYvk?0ETyKNbML&OMi;oH~l=JMR^BFS{>vu(@yV?Pi> ztP|Lfx{QJJLVZDw#%-sFEzQ$+;Vb{6DhD32&RUCKRr+=e*Hu7X&?Oirc3*WLfi2M; zmMRL)=5aNyaZ-VxsMC}X16BMICSUT3#qi>f!h=Eq2L{PqVq zPiO(3oCxEF-opc)uj`F}*kY#xYn%a?PBNIaLe!%k*w?iLELs!L(gHSM0X^u&++oG+ z2!G*q(iHDYuYo88Zp>NNV#)`7evWFev-rMiT}42i9$lHK-15J1{|Ec?z~&0oi>EGY z&Y51X4=qxmu}FCSMvXy8{t7`pA#r&3vatx2{7y z=Um{?96_@EI@pnWSH8t~jW)I`_|3QCKntxDAIxW+&W$83eLDo%2csuWQ}F3_%}n;G z`f~2n0)_o2cP#8EwY+ADM~h@G-JBMsgSzjwM~Iw4tMQPW)Ah ztk}8U!d2{fxRngC}Ceiq!r7`4Vp0;S<^QW;-D-ga>0}? zlMT~hPfbOWr-7JIZ1(-O1LpRW#0NZaWsMg<5Np&KliBVqQayxwgbLfX*ccr_3jd+9 z8ME4|z8m9_v>P0?zTUoSbCZs10xcJ6dd9(ou`Mo+IzsAMTjj-=^4dlk4JE@mx#qjC zvu#!Zc%Fu5l%6UPnLWb++onIU?g~tC=eGq7wj9j$T0}~&A2%5K9XCGs`RF|}E)NOeQ2;aCDU`TL@{Q}jQ-t$+G6A>Lsi zTE=r`E36E7nI>E!qdZHJdE17IFEmDazU=`jl!(UMU~_=f2H(+KW>^Caow(@RLzVHv zH%Lc@o#)fV%){dx{e(y?4=@L1#RjxdsXX@TzXF8uh{IxRQ=}X;o27Owdte3aR^(;C zsZ6K%J#IK~aqNv=5v+b8KbWJ`z=Q4SC#JBEt%M*pJ~~y;SVr6nip2U5A60$F@YZhL zB@-Ed&2MIsft{9pjaOvb-+M^-6S$#MQgZW{=~8664!)6$#D)qB%I-S_{&(`i!Q?9wOI5{u^1kUGQG2j3C_AlG`I3B&cnboWN6vVDrh9Cp37`v~#S z51H|<5&87P+jK2iqBDkopSE1}v4X$kMBJIl5@y6Z00S+AhanFOa_*C$y0v@CUW8(y zI8i{K_Z!DA@_m$dhF~AR)f?5MULo7eoT)Wdfd0;(anBV?fHwndP8EBFa`-X-xd`Eh zgLBmsY0J0qOZ7nK=HL!(BK-N}XzSgSFqu#1#1kP6#sC>02JRIMfnLDi6yb`dliolg z#Ul$t5Iuwqp8Ib&K}KNXGuSJ|<}3H=4Gn(Fkh%A#5+=NY#x+{TVW$N)Jzsb`xe*yS zK~UT6mXSJJ-6$f&8Ci*PoqKtJ>lB~7aW<%r--}8iL*ifi1}wA*M~bY zp{rH1Kmlx<9yf{05b+*_JMka=-pRUxDvymib7*%7Aht^&%av)>Ytx{d@;Ee`19ZQX z&jNsW^i<{k+oz&4;3uYOa)VR9Cp4L!8z*P;27z!89NoSZ2h@EyA>ErU40Wxu`Lp=0X&aqeUtAt&3DKQM%&#`{KW?3r8kb z`Zn0`z|EpPj69Am%J!3PBFqjl0~1zJJW*@!^@n8g>MqS|q)kCM)=O&bc)52b7sk$NTZ996UpXCuD2T~ue*iud(la)sLtZe_W`(+I& z@Ep0UB)E9f9X0Ea&P|f4Y%X7C1K?G8)P%a9nWiP5%TI<|(sfpw8DTT_c<@f*u{ce; ziEw-S?6?*fD~!;4Vfga2g}gEW%L3IBzCJRNRnEkJRZ7Nqb4^kF(G6%&#CF`_VEY_) z@YDv$j+519lRXt%pP4DU?Kg6|#|+-b2Zi530%gv>K;T&kXRhFGsKX+_FJlyp^mKf~ zAUL{_L0PHBZv6#>W(iUfXhc;m1s88zQ)dUMFg|;7?;8Qc8nj;a@;2B6-(l{4l!Oa* z^*~FRCi2^r69xnOG%0ffI5{Ysq8wv#nH)U5gL%y6=CbiT#wrlV?S~H{@rK?G{LDwT!skUG=TpJB(W@z~}v}mAKM80mEvr zo+q@WySokKK3JE+?)6@>NGlOK=H96G4psde_!?5yD3#f}r$9SBb=B{e_jPhX%VzkI}q9&od;PP6!&t zw8r$nr4HJ_p=;hMep|sob%RLNYmz!}(^TE&@1bNQL9m2<{I#btC~X>|!U&cQuhZJf1(zLYjrwp+orgqI+|G9TB6@K7EjT0!Q^X|J)cu26OXnFS(XepW}J z5={k7U;XtQ%xV^mi(&|3Te22lLiF6qX)>28*joyx+^g8C(1u`1Eo-JyENr4z`?Nt;%n!RrjdYg5vpfplBP`p%;tT9mw zzm=_P61&oY4^qA-)&IKr)Rog4mkhoe(@bbu zw@)p_%<&z|+5_G>pz;-sLWzBuwh16^+WP%1=fr^hjDpx;4s={`@H)3^!44fW-V@82 zhL;%!>ni*NW(#Gzmta$@`r2#vSO8D|?R%V7mYyk1*-y~|%D2~U(8 zMsmRh{pFx3Xm$F_*RB)Kbt!Z3E9M*jaY_jGG)oY2gY&N#01?M?QX{~kB zo0=U5QS180MX^6?(7r!2@HrQi0^LW7;7H8rrXo8>cn}6k!!11_F0A67m zVYpO2dEm`q?^Qj47QkguP$!sYj~SR}?Jxf20VGJSk7OqJzo>$WyE%?nK2vQWnpJo7 z%@+&;n$x3oMRJ3E3}e#bQUtA-&1Nw6b7lCPya|8u9D~FZhkB`+cCG+r7x9&wP+D@3FVbC<2?0oF(8Kkda0&Ky>$ zz*Nfb!zY)V4|>&AQc$z5RL$e;6fQ|6yrg=1FR%=5AgTpKzBZVVzaXVgX^)aTT)t8> z*RtVEd{Ptm2$96to?c&AyngJ0S=%hB#Mv#xk*Vb;6|DZxuC+0}4PXhu;1wM6$KxAY z6V{_3*Ff>W$kA69ZPHc`_71vBM^N?3a_#|6F|p3ZyB@ldcRyAHd$V4R%c6xu^c4XY zgPwp7gJI-rcAOu;(MA<4S~Q!|stdm0bC#Y} zrUl2MNmdw{Cv&8jOc*yftGkp8Pe3QlVt$6wa#io!*)wjrP+IMeWafy(G`OtiQ58*1 z^`twNS?Id@+=hOC7X2QuWKUttcj1pX_BH7;B*tPOg_2o` zGTXVA!gF5el?SxQTZc9cN3jul#x2`fCE49I?3UU4a!!9X_bI4buo?I3S3MH~(RJH7 z^@a0osxo+dm9S;Kwi>7A$clD?7(_X~UcFzWlgE+9RbPfbQp)?a!CEi(cyIu3Z9Q$1 zCS5A}q->~icq6&5be1SkH>piNW?-IlxrdigPLID4P0_;YnCPof`Jc^he#yVdPu2{1 zg>y<@v<`o%j3l=^nHW@yBoZa8#+M^|bOY?0#(u8iH$CnRe?-=yg(OQEntmMx{rO`B zmO^7B4&T%b?3%4CbHM%a@Dtbmh(-%orW^0W!I|LD>Fft61k<3QhYpJADup+O zxk6I8N&~(RqS4-r9QnQp_TWl0sDrn#TebO4#cR|l@nF33v}HOvu>jrN-GIFO3vgKM zlAhXz5{T6Ik_;$NIq2EelkSXjf|ItuKd+h{M4WjoD0nIt!-s)Irn>TW`(k@ zIcXrDF8KQOY+>8gFVs!)w`T>EoFY0-=S+>%af+kYTmyK)o~>CDAcU+K{#Bnis_|rd zBO2>I;VYC2kJavLSTNr+PqbLxnwSkQly2Ns7+*nStca4MI5^2u{>p@hT!ogppKhSO zjZm)bu74oqIh7@ZMn-s7vZCk#p9u6{7w>4PjaEnKOXbmMPqnS^Y*f(O91VCVQZL0z zy_@DLhp`nF3~BLPIUE*0t=OHs%|u=tt{WrE#Il?E{)z~PCFt)3nofJ0Vx4;UtVM$9 ziIs*f8Qs4)cL^iYI2%$G?Wn-AwEta*TD#vi4il*J)rLMi7Q=ZTX_qRAnOa~vuS~ko z`zlZ$!KUumXwt^we`}GXc(nu&jFG?|{CEK8tOnhStRaD3!8984VF(upo05O(^E#D9 zGLhQBh;RkF{h>VA7-=NbI>wh9J;>bxat-);>rR$c8dHAG42X`b0xtOYL6XZJh5A;} zjoJ(~+W97+zu)9Hp9{^#e4HjHuBfCx=#%lTRr^gEl!~2EL=-a1lUWdBXu{7FVL^ra zwP<68DVqfO)oJhJi|wIRN9=FHPiT>2odI;XTd~hR0)!DVgVq)?K5`pRo3&p5&m8|r^J-MT z4AegCTUwgrk0PyvxCwMza~$#rA9|Wwokm4TT5Dsy;I=`vk4=6gCax1WxAIVsKBSSe zs5*UW|J_P*+GV0Tq)(T5Qz}BO?VCh^s6d}$>X&QGiH7>DMoV*>O41d7Yt9S2n5b3^ zFc{CUlo!=>9%75r$S4)a@;{HXJ_&^~5{3`y*aCQk^||GlVHHrG2ch{WI(hr{_MgA+ za2MKpvGJ?wBVf&b2`|?8`a@dH%Wg^oS^S z>5Wch`&FmRT2;H_zKDza{XQX`_u?spV;PZtI0zi?W(fNJ zlBRX>zs@i)nppor=i|eW->+?EQ}{?FG`BM;o%2&?W%A~t_IZb*9JGg!;41pq0~lt) zYt}<8{s}gjRh4NL(DQyD&{VtF?$6%Oq@v?EDCh-!75EkQG*a+VOzy3f2?wU8IS0J! zf%EK~wy?r6x&tg|?^F*~1ricM2YP#mGGTsl;KrCW67NQ$jbEW#XxiH%Cyg~-#0W4y zsI{S0DY^)J>h{&1>OJ>eSpE(ksv=#O9(TPh14St2fa--VSdiGeuN63zhPT>1QMYy| zMpePYqSRnK;0#C&xH|u6`b-xEhJc}sS`|KK5^LC>?(!>tfAZY7qWiLzxDUW-v!YZ$ zZf^>zvpSp4O4GQHl7^A;TPd1o+dqkK{_St3GOZEuzStc5DN!VIM1nyrKw}$aWqj(f z#IWZ?K7X68{h2Xcw$U<+Uf+HGpI0#H;_(&7Plp672YA1+3ARz`6F%kILrQC0nP9!j z+NeW)fPPe&#ykDkM1jefwV46@$dT;i!Lvqs?WK6KFSgHkZsgV22iHP4$-kq{W$s8c z@p>AMSzAi&!YeVX*1sXPIIvLtGvwuftvukfw{<%Gq~_zMNjoRQ#&BF6$Yfg8AFoJX zlTn!pQEL<&=P2d9ICz2loTnG-21>>QjX=p(vh0*Ijaxv`GXt{3KLF+x%<5qpR(*X~ zx5R=dqm*>4&^A&XfDea2xbHtsHT3&*JxRrn+Xa8B;)kUyE2QqfNa(;*F|K=_ zn5Jp^8w26jY<%JZtG2XvOv)_N%UTjNw%>emIC9nxdrhlZY$aO$c3OP>4 zolor*BaEBEFNSOJQApm~;0@#rF%)NdTa~y-8kVMwiwrpJ9cSB~_)9gy*DmuF7fB#% zw=3E;&8|8hIBWY~X;tjnRXnFvMpyCA&M*z?qXk?I$ihNig`RSaey-FF@o zz8JB6_#+U^UB{2l3Tlt%-pm# zC;HeY-?M?^GH#A)f zdV^6Du4$VjSCf+!T5JKlrr38Y`~-OSacEC*R}PxkJb76-PTs=&<3Icw_nv;&hXU*V z_gN_wXMgF$$mzUZZUj~9kBw$ukj1L$Q4hjpG)}mLtdo0EuzdFuC3?Q=ZO;JH!g|KG z?)ofsc1TWWmS%YLBTl{@)&QG^kRVI0D7D6^GU~;CPD3meZ}a`@5$LCa>bEZMEDl3t zYiWtLH8blEGE%gwzVP4Ee^pK{9ulgl-2ypz{_?%r(`bL#9t-$BSYK)`@A>pi?m1|ss)m!8R2ayD z2bm)^ILDw)5;8xuM;yz#ZMV-_(UGi}(ChH6y@NZ8lLrh|)fi%2^wosP2?=fEXw`rh?WH z5FVG#<32YCA)Y#D{PwHK?Tw6!hdYU2Y%WVp2emrN4~gt3B3Sogpt=aYj7oq$Q+!zz zuht_G707du6dW%^dZ~r265crTEK=1~aN^x7f=(tBKzj{)4%5Y#ZXGkkN=r$iR{?M? zY$joitu!A=XPkD7Y{6h6ESU=SPaRU$oxwVuS0Re-DN7j+vE~@BOt^&)tg7Im6C?gLK+u!4NUdzU@NL!YuZz*MrHP%Y{8N-$ z#oU4hf#|AU$O%enH|UM!RGyhX+~P~%JxRC}nP7lLM~h25Z63FEC|GB}{~X4Jvd6d= z@qzJDh9M>vGKU5=_k!jcemO1L(Z!g%gE9fqugPKgouWDap`FM@G~}C9!sm8Nu0UzB zjmP<4CKT)Ct4xM5JR=Y&<5@4ms2hkY%^RZVL7C7A##1UjJimjpz53PY9?tgB2_b|c z>?bCJt;a6sPy|!@PM>yFz(lgRolNY1$CYVTwsZX1jE)L5AWfgX)*O2zo2wl~!ucSy z$}+4FUE^a*o60`qACBYLi*A5g1vbF6xoARn-~%%K65p9Jt-LW zD&$Uet+&yfHU14`e(#AX{L5^mw#Bf;u=H|}=L$?D<~%eGAW&P$s|;`*^iC#3tBpdn zrQsMBVQqYIM*9n#O-{BB=}Kn+CZJY_lGkav09?N{yxi|%oO4*5Fy@0wodCmr&(4dK ztQ=Zvp%&b|F*3)J@u)z}R=zX{+1#r?V`-0$dX{Qnq%ilgy4_G?Zz}Kz!nQWV24Ujm zzqg{I^r}!-Cpqr?Zl_XhF77=BNgoBtG>fSf7Qn=EwxT;8M4}PpyTGmZvpn zYvGCuRCtf&%0@qOCGC$a5ZE(i_~XmYNm+P$5cQNHN0 zxhy>w(XS=hBE!O{NZ?6UGbnbdnZgkOlE6SI6Tl@L1YfOG=Z>iaalmm>`+>$$fH@T? z<`QQCt#GW;u&{=kvC>R5r1G&#{g3pst-fTPLOBJZ5`qhAntrhS4dR{#EDh4f{0%3~ zb0*MJaG96cDXip-=bg<=Gsc_Vct4MFycp#XYY^=UW^VH~mq=vmbQqqdS}ppr7x$6Dg{EW1R{Q;=K z2Yi41qT#3${Hdc;yyO%^*SzRgdx)+G@Y%_GcQBY)6aJM11*BA&S6_35?HRwGeYHz@$GSnUkp$Mw@b zGvs$a)t3F*5lB%mnjYX<>+hcM-6oD0+$qn{cj9~xC7gGF{;RLYK}g7F))lK$izeEv`O(*ap6g^3>x5 z-j?dBfNK@q)$z_HqJC2<(iRl~kYd=Z3r(4VPf9E1cWw9K<1%1Hl=pZ!+BF^)#ASm8 zXh&^Kb-$3xN2$hV`6+A=laH!bskU_qpQ_>QWn822TK2VYI@GjBym8xb_XOln+v%f$ zE&31O*1@PlB2*07^VkDVkPVAAW_62{WK?*8k+tZXn&FkQmpY?ye4A!qmm=`BS{2Zxz^xRf{Npc;_<(o$=IK zdE#Iqd&BQB5{Xeo!^*7!GoH*yUm$J6xICU~YEYji$&)ud=eB8%n^sPX(*=Fjf`YR= zg0C$}QOA~FXMzObOC5|TeHoOEBgBjAQoKNPPs}|`OMX~q^UG&LkjBv(Q`;DWBn!F? z*5epB5HK7I9k1{~=Vl+F(7@v~$~f)=aOZ)=Ed@gklVA0)#-x5Bx9+SRcGA8Ob&(ugn%a4maqeFs^-%S(PCI1I+C1eE7;bKH9Qw$6!;X1Dx1 zAcw8Wp0xOc~@@yJ8YTL7lq~xnD=*_65#5!9XPWI;?TF9TlZur^q z2g)lSpI>v>y2>aEoC$Uynq8px^VaNxG!rIscz*t&z0R#sN_=Eom!Y>Q_u%%Fcj5YD z%V13O?qa4BO{H_8!O)~Pt`o{~vHp`AQqO8b70!bW2)BQsFMRf=IYqF|9Xi6AUv3Qc zczK8UNk5xeoM$*?^`=h4?$HIshEohE`iMf5AL)@({s;r9eDkdLr5xtJG1C9jzNHCo zBn=r9N7rDP(&=jLw;@j;5YkR!GEr`6jevV79(-olz7QuK#LDD)p7z~FEhpgt1!lPxp=nSGf4SM9Y48z z`=0kdetYHD*HOsxRR#ALJmq>b|H#popq_^~!rIJVWZV74$V(LNKAFckP+2}n!jY%A zvElgXeN*iF9&h9rYX5@~_T2?(@r+B4qjI=ql;qFU(qP1nQVH1W_SfOKjDIswd%6N? zI=gXL*0Nx4Wr(`6ViLgE`8!mXh83990z}>A3hN3FkA7n|>#nx4Fhi-`%H|6X;kB_)@-cl$2DFWWEq9-3 zUI0PChqwIq=<9=w+$0J=8MDU$rcLW7U&El6e8m1dC9}E&15Mo2EB!G^S^iy=9qt{# zfFRx{%qM0kVBQXSXCdFRhoUAMl)kLZ)>25)Thy#DEHt3$Eq&hjkXcPE2vgMgFWRJF zv3N@k5p!ug#{U>=9iN|jxEDm@D4_koFTAzqcW=Y(1@7zTT2gLO*P>>=lk%B;ivw_} z=N$iAC&l+2Hoe9-P!fc$ByYCmZq+ta>4 z8ZDq04r-3LXcRP6+*PuUghnY2LC<7GmhFskqu5jjXq$?5C@cRq$UJh-{0CifF4nV` zk*QJdwP_dq6#TQZXzmcbTJj<@f=Bs;4inH$5I4pJ$16WIG?J3-2obkqr)=fvSrlm# z_kwQ_JMt2vRy_O{&9)kMVb<|PV89YGa^JG0MkCGC}&vyq^#?Ec*UPq^f)@?i@hEG?!4_LOpv2 zf`>E@yVUXw$YNbop`h}wLiEZn)>TI4OG`ad?S9kQ25plYXW(d^2S+S)#iQju9^pId z7@Z>YN7PGQ2R&zUck8J1336vm;z`5qiMdZO`D4nGgWUaZdD{&>iH#Uy6Vd9bNcUxi zO$9&Wu4<{=Gwg)M05gNH3)&w~Ot%w+)~cDnMPM6pEVT7Cjxjq^^IpguH}}aW{w8~j z>X~()jJl;vWmER=SCp}l?hxMX=*ZA@hHhevU7JKn?Q5|A$T?qFcMDME<*1$XJ6a7J1mV#Vip|GOIR4z%PM) zk2D{_XkC{+_GrH0`!p*gw(7vs8!O>M8osohPK`N1p@_9*4a1VAQlqT%gj*hCbgpJw zkRL83b<#idKO*!1j;&V2IusOCcE6;lNry8w0Y5Pz`GQFOajE^|LD3WOUHJ%{1V92H$b6F1w z!N)p2O5%^_R+GkgCUL`5t?it2Yo6CzMiig}=uD;4IPqu36Jlw4Kq0)x~8$KQ4if;##&ZI5-wnCt-YR7&v4lSYM&whWs2{kI`*F+ zt=+kdP%N?O!D^L48{6TlED+5R9?oVYnC`X|VDI7+=lwU)m37qkx%M~0>^xc3xy&p?1H!I+OpTQ^*+%#!lo%^PW zuWkUBnSS)sGQer;N4)8*d=s{ImV5+Pm#}8qk0!y_+V3ktXe1G(gyeQs$k8R8t1?oYm}h5R;9MnCjn&E>pb(*+1b{Ac zWW^7oJv-`MwTS?-zz;1vsObE&$Gi0D@F$)SJa)%LH-z%tRpha5C)!#J^0%T?h%1%V z7_ulSI65-NUmupZ{nyV^Jh=-e5LbX1B#K9kg@v9ccGPamBy3NhXJwXNni@r{C=$Tn z%Obx}ssaYn$I;-hxJ)qCK)8BYQn1+vvN+R%qCcJqd!W=UNUNmuJRzh{QVZy=ip{f8 z52yq}kdg#Zx3j839fAEu4CEYa75ZzYMn+AG^5a{UnRVg>=Yi&mM)@GvcLajEM*kU% zP`a~Z>XNd<=#m0bEW^`)HAKX|-)1yua`O){;s?$`oEW>rL7&y}ihbNwlbSxo?!9Z= zqQ8Cxd}P5peQ2~VWgyXbrQx-8q6^^~E`E0MZZe9jeP}?0?$K+##@qPuPSOIdUMzC& zkia;^_NvuWCi?-yf2H_xV{eOK6@?@s7~GWtD-u3(-OZa~o09kk932Wrov~wQ4HYa% z$~^*pJHT07f3}js2LkEivl+K9a1AoWLA|;~zq#+{&+)-ZxsC-}#Ndo21R@2haP}D{ zy>v)v+3V1-KBSJ?*80^1w+X5Cn|<`3Sn}K)WI+W~78tnU49^O5n-> zn}0A&>BFwMU)er@j_p68ay}nWK-=f#$FrMZCE-%olwOF;StQCx7He%StlB zJiy2%I>{&Y2e}O*ebAA}_$upD5O*oKW80*=V3$C`GvIt-Bzr-?P}6 zYbSKjzmS@$>ip_$o7T|Xx0@HykQ11^vIMws~%0*$viUrJ|y|?aehh;d@HwF@5z}g!%1A;pJ2{U4+-?KlT zkR|2OwRJ0eo!T0G;SkF+PXjR8)G0Z+i-R?KDjr?!nr#fO02jNaCjPiw>lt4%wcD?m zCd&s?c3WWv<-#`fzT_?T^QT3B#KaagN3XMsTtcf>CSuiovYom`^nnAj6VFx${a>Vxcc5b<`=44uw zhYOSj6f%6(vJXjPhNAdG8(9F;L+wy6g}}R1HHkRXU5CZ;~ApQKsAz7E1JFlZpSHBC(B>cBSGaQk?`Ih^|eva$Xd=fZ) zP1FLbQ^fPEpcG9XxFTL|vHND<6e@a#lHH8n?Ou-X4<ET$Ozz=b9ZO<*e|be+&?&EcycYrM6-R38cj}^IJQwv_4;7@Kvr_jD#2ah%nZ@O* zm@7-c?JEag{PhXuB}8ld(&}@C9M#B@!r?Ht+FT^B;ph`gj=dzas#=+NuSw!bC>3bZ zD)Ed-Oe`<)A=(X1RJqP0%+bSPj0&UI&!OkvP)em$v!7`Om-eHe`Md`R6^Fj?-X-EL z;i;4(u3w(Y(`{z^O3;C=RjNC=me!DFoUcXE_gF$@V=laG3KK4Hz!aR@t|rg&u%cIS zsvRPk=CKk<=pMR;uMEE~S`NFe%!H^~YtpnF^N^`>d~P39ua+|=*jF>(@6IG3N2&oJ{pO?=8AXDJ1Q^Mt-aMM;# z%oGYq6&3F{PZ59O85!Sd;-1bSvcQw;%j5^7i!*lCh79w*K2|)KyXsmm!8_+sX*jhV`rIMwU-B;}@8dG$qm}X@A{R zdLBj!_3ZiNSlU~)ZItC9-p|;NDGB|zU>clU+AQ6Ioz2u0&p0mf5d_9Dva#RFddiPE z6!G_hl;&4np;mjhA#@}L7))XB@g*$3OYdAJG3q%CuvEOt3hUi>W_1sgMa!M@?Nsn= zM;BNLv>{2&a%`r`OYiLq)Zz09iMb0J!d|NOu98~RM4}DKDLrl~9JzHRHb;{Hr;3$` z)qOE|D_hfQSA>JMR%sn=!&vB&e!eJgJLw9J<~gbdMHp!&-Fa~dYp$kF5A}u)y{H># z4)9%Zji&#I1sWHWxD3YKp6Bs18KcGJq_dE-NzkhKuc{9y@)V_OOrtP6qckgcYOXrU zE#B?J+6(Z45#GycqGuN*rj4sjEQSk1=g!_XsY`>{-LzH4Q4Q3D%cUzCw{~l8v9Tn^ zL`sBiU?J~b%`GX3ShoEi8Q*gxwigRqd&w0&r)gbt7HcdpyllM$m^!AIm>*{n$|iTO zae7`4L0Df`GNif`qF2;AcJl+^N*+cFW3}zFXH`e}EG$R2r>98mMMxs`$WV<_OaNV> z+*$O z?54VEkPNkyJKTaT^{?exJ!oiYAhfM{m=?zkYt7X=8{IFZC~N z{E3yY`g?HCSAFXQq3o7X;01JwypO0nK(hn;z|0{2HH>}MlwL6V0xf{7xaI4>P+7uD zN=n_b3`mazW6%Fi!FBBLr@-|cE=ed7WPG>9^PXn`nTKqAl=F!OShmOXo-u=E^m2gk zC!4gFwzq^*Q$Ec{+U;I3SQ0XBVr3IgYiXlnlBt+Lb^5nJBQ%Xoq-}^73<;x$j$?n- zHj7rj{|Jn@xtnwi_IS%yO61joV5{71(BnptkqS(+KJpv-&FbawEKf2$nE;M|ONTDA zO^kn~(0;||GrEAk=f_7%8k`AD;r`yH$>G*dBdnl1_4W1WlCT@ilvk6j(@CNhTcxOnG_ zx~77#p>)Yx`I36pwTG03H0~T8fjR5hhG#;C{$;&ua0;kZNNOcUTO~ zLuZkJ**AjGH?oJ}-V(FfgvF3s28H~3eG)@2L4f_R0CHxZ{oM!TB?!-5W>9O| zVMp)iL-om)6=cmmm^V)qWvWj2?VTey4~;#aW6A_`M#2D`35GN8qnMfG?vzJA@mO@> z;YimYR9EoL^2ZrBpdAEAD;Oy-0UgSpIL&a$y_+)@)Sj6xTR{0OeDkaAuy?nFFW$$MXJCfOhz_DX+b9M9}nz7*Id%yfVR3an$9(f$pJhj#Z zC6H)ML6iRa<8cC8r1QB8B2uCrw?4Kp!QoDYl#B{&pTtyyYNA{{X_b$f9@0F{%?MB3T8 z<{Uj_jX_QiV1d19V81DnLR?@bVwQX`9YdGnN{4%s0~~@@&|6R9fYhphEgLJXL6tD0 zhoIE#1?geyv!9SJQfR`TL1a9AQ^i9RAG^O9rC^i7E1YM0$~B2e@TMdMu~X`Gy1R)( zFQu$X0C;vM0fE69AECP*A^A~?wj$_xj9+M)2)FnraBt>Op*Zj;T~zT2LxRk&oH#cy zofJNg4KO+S_#=-?#}9yJws=)alLP5uZktbDEMm-qx?de<-AD&~c+M3WhHohVj zIs|0wkAOYZIyGLrCWdr3=}^sJZ{Xi7VjZFXiW(GFHFmENgZhnUUK2v4vvsHgB2eT* zKOO78AK-fzT+9;Z88Ii zkV)1%8T6J6A3KAB)%R_)df`yp80*kOmQ+DXklR&fBh4~s+=x?6s;>|V7`UqjNPjE8 zsgWn&R-g;%)tc12!pJ~XD4Colf11;a<)nKP%_76?yhHs2u5>N*uKuATzu-|^am}lt zmra;-&&lf0hu{`imS}evP3*XHRONjC#=W9!I&T|X`nCEb>Xcv?_`c2y`8d(QTrB`X z|9oTWv0GJ&C5oxF@oLInw>pT4!J{sx$q(E+R6`LmbF%HRyx~aJX4qF(S9wSj zU}mq$h65if<4rH<+ebxXSwnc_$j_323~=5Z4e>WamW*!7H0d2&wl*Kps6I+N3qK&) zqNPxw3TZcEZs^Mc6 zFIuOxIKH)T-xz3fEtwQF!q3vP1o_QHf@(NqNetkf%FIf(a@c!zX_M$x`Bpuv}ZldVhgnv9tD%AORnxMCV#e85ADlypS#>O{mlGJNG?pTva`77D_X%pF?r!pWRy_W@Jz}NvyLZN=%CZsoZ^{>#wb>s5GSt zrS1Q|+AeWY3y#_1Z!BAT=It`xtZ5=kZKp|lBqwyAMUJrw9~m?N0N~7LWdH%>;sV{A UFHw)_MLsYMivj=u00045TH*OOzyJUM literal 0 HcmV?d00001 diff --git a/data/datalist b/data/datalist index 59effdb8..5ff0fef9 100644 --- a/data/datalist +++ b/data/datalist @@ -2,4 +2,5 @@ AIX_NaI_1: AIX_NaI_1 BDX_LaBr_1: BDX_LaBr_1 clermont_2024: clermont_2024 clermont: clermont +HD_CZT_1: HD_CZT_1 heidelberg_2024: heidelberg_2024 diff --git a/man/AIX_NaI_1.Rd b/man/AIX_NaI_1.Rd index e89cad37..b3458674 100644 --- a/man/AIX_NaI_1.Rd +++ b/man/AIX_NaI_1.Rd @@ -27,6 +27,7 @@ plot(AIX_NaI_1) \seealso{ Other datasets: \code{\link{BDX_LaBr_1}}, +\code{\link{HD_CZT_1}}, \code{\link{clermont}}, \code{\link{clermont_2024}}, \code{\link{heidelberg_2024}} diff --git a/man/BDX_LaBr_1.Rd b/man/BDX_LaBr_1.Rd index c9ca8e41..080873e5 100644 --- a/man/BDX_LaBr_1.Rd +++ b/man/BDX_LaBr_1.Rd @@ -27,6 +27,7 @@ plot(BDX_LaBr_1) \seealso{ Other datasets: \code{\link{AIX_NaI_1}}, +\code{\link{HD_CZT_1}}, \code{\link{clermont}}, \code{\link{clermont_2024}}, \code{\link{heidelberg_2024}} diff --git a/man/HD_CZT_1.Rd b/man/HD_CZT_1.Rd new file mode 100644 index 00000000..afe2983d --- /dev/null +++ b/man/HD_CZT_1.Rd @@ -0,0 +1,59 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/data.R +\docType{data} +\name{HD_CZT_1} +\alias{HD_CZT_1} +\title{heiLUM Calibration Curve (CZT)} +\format{ +Objects of class \linkS4class{CalibrationCurve}. +\tabular{ll}{ +\strong{Laboratory} \tab heiLUM, Heidelberg University, DE \cr +\strong{Instrument} \tab Kromek GR1 & Kromek RayMon10 (SN: 5261)\cr +\strong{Detector} \tab GR1+ (SN: 1625) & GR1 (SN: 1663) \cr +\strong{Authors} \tab Sebastian Kreutzer, Didier Maillier, Norbert Mercier, heiLUM Team +} + +The calibration dataset consists of the following objects: +\itemize{ +\item CAL_GR1_AdamiecAitken1998 +\item CAL_GR1_Cresswelletal2018 +\item CAL_GR1_Gurinetal2011 +\item CAL_GR1_Liritzisetal2013 +\item CAL_RAYMON_GR1_AdamiecAitken1998 +\item CAL_RAYMON_GR1_Cresswelletal2018 +\item CAL_RAYMON_GR1_Gurinetal2011 +\item CAL_RAYMON_GR1_Liritzisetal2013 +} +} +\usage{ +data(HD_CZT_1) +} +\description{ +heiLUM Calibration Curve (CZT) +} +\examples{ +## Load the curve +data(HD_CZT_1, package = "gamma") +names(HD_CZT_1) +plot(HD_CZT_1$CAL_GR1_AdamiecAitken1998) + +} +\references{ +Kreutzer, S., Martin, L., Miallier, D. & Mercier, N. (2025). +Environmental gamma dose rate measurements using cadmium zinc telluride (CZT) detectors. +\emph{Geochronology}, 7, p. 229–246. \doi{10.5194/gchron-7-229-2025} + +Kreutzer, S., Martin, L., Miallier, D., Mercier, N. & Faske, M. (2024). +Dataset: Environmental Gamma Dose Rate Measurements using CZT Detectors (1.0.0). +\emph{Zenodo}. \doi{10.5281/zenodo.13731839} +} +\seealso{ +Other datasets: +\code{\link{AIX_NaI_1}}, +\code{\link{BDX_LaBr_1}}, +\code{\link{clermont}}, +\code{\link{clermont_2024}}, +\code{\link{heidelberg_2024}} +} +\concept{datasets} +\keyword{datasets} diff --git a/man/clermont.Rd b/man/clermont.Rd index f5ecb0f3..bb06076a 100644 --- a/man/clermont.Rd +++ b/man/clermont.Rd @@ -27,6 +27,7 @@ Clermont Reference Data Other datasets: \code{\link{AIX_NaI_1}}, \code{\link{BDX_LaBr_1}}, +\code{\link{HD_CZT_1}}, \code{\link{clermont_2024}}, \code{\link{heidelberg_2024}} } diff --git a/man/clermont_2024.Rd b/man/clermont_2024.Rd index a93b0224..1351a45f 100644 --- a/man/clermont_2024.Rd +++ b/man/clermont_2024.Rd @@ -45,6 +45,7 @@ errors. Other datasets: \code{\link{AIX_NaI_1}}, \code{\link{BDX_LaBr_1}}, +\code{\link{HD_CZT_1}}, \code{\link{clermont}}, \code{\link{heidelberg_2024}} } diff --git a/man/heidelberg_2024.Rd b/man/heidelberg_2024.Rd index d236463d..45ba8af5 100644 --- a/man/heidelberg_2024.Rd +++ b/man/heidelberg_2024.Rd @@ -45,6 +45,7 @@ Dataset: Environmental Gamma Dose Rate Measurements using CZT Detectors (1.0.0). Other datasets: \code{\link{AIX_NaI_1}}, \code{\link{BDX_LaBr_1}}, +\code{\link{HD_CZT_1}}, \code{\link{clermont}}, \code{\link{clermont_2024}} } From 7444fbd04d072f59506af209dd89d5c1aaf93c82 Mon Sep 17 00:00:00 2001 From: RLumSK Date: Wed, 16 Jul 2025 23:05:47 +0200 Subject: [PATCH 5/5] Add PR ID --- NEWS.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/NEWS.md b/NEWS.md index 3e45cef2..08a25502 100644 --- a/NEWS.md +++ b/NEWS.md @@ -4,7 +4,7 @@ ## Datasets * Add reference datasets for calibration sites around Heidelberg and the CZT detector calibration -files as detailed in Kreutzer et al. (2025, Geochronology) +files as detailed in Kreutzer et al. (2025, Geochronology) (PR #54 by @RLumSK) # gamma 1.1.0 ## Bugfixes