タンデトロン加速器質量分析法による炭素14測定データの解析パッケージ

Karlen, I., Olsson, I. U., Kallburg, P., and Kilici, S. (1964). Absolute determination of the activity of two 14C

dating standards. Arkiv Geofysik 4, 465–471.

Godwin H. 1962. Radiocarbon dating, fifth international conference. Nature 195, 943–945.

Nakamura, T., E. Niu, E., Oda, H., Ikeda, A., Minami, H., Takahashi, H., Adachi, M., Pals, L., Gottdang, A., N.

Suya, N. (2000) The HVEE Tandetron AMS system at Nagoya University. Nuclear Instruments and Methods

in Physics Research B 172, 52-57.

46

M. Okuizumi et al. / The Nagoya University Chronological Research 4, 39-47 (2020)

ʼn ŤśƚƏƙƖƏƛŸś

#####################################

# repmesG

# plot of logfile

# note: preload the "IRAP" package

# 2020/6/6 by M. Okuizumi

#####################################

repmesG <- function(logdata, sample_name) {

plotLog <- function(ypar, title) {

df <- logdata1$BLOCK_DATA[[sample_name]]

p <- ggplot(df, aes(x = Blk, y = .data[[ypar]])) +

geom_line() + ylab(title) + xlab("") + theme_bw() +

theme(axis.title.x=element_blank(),axis.text.x=element_blank(),

axis.ticks.x=element_blank())

return(p)

p1 <- .plotLog("C14WndCnts", "C14window")

p2 <- .plotLog("C13cur", "13C current")

p3 <- .plotLog("C12cur", "12C current")

p4 <- .plotLog("C14.C12", "C14/C12")

p5 <- .plotLog("C13.C12", "C13/C12")

gridExtra::grid.arrange(p1, p2, p3, p4, p5, layout_matrix = rbind(1,2,3,4,5))

####main program#####

library(gridExtra)

library(IRAP)

logdata <- readRawdata(filepath)

sn <- "781"

#sample name

repmesG(logdata, sn)

ʼn ťśƌƏƋƖƎƏƙŸś

#####################################

# beamdepG

# plot of logfile

# note: preload the "IRAP" package

# 2020/6/6 by M. Okuizumi

#####################################

beamdepG <- function(logdata, sample_name) {

df <- logdata1$BLOCK_DATA[[sample_name]]

p <- ggplot(df, aes(x = .data[["C12cur"]], y = .data[["C14.C12"]])) +

geom_point() +

ylab("14C/12C") + xlab("C12 current") +

theme_bw()

print(p)

####main program#####

library(gridExtra)

library(IRAP)

logdata <- readRawdata(filepath)

sn <- "781"

#sample name

beamdepG(logdata, sn)

10

47

...