blob: 4c9857b55ea96ec4e0e979ef5614c20bc60705b2 [file] [log] [blame]
Adrià Vilanova Martínez20de4092023-03-21 17:53:14 +01001library(svglite)
2
3filename = 'KCl'
4scale <- 1.5 # Scale compared to the Na_smoothest (for the arrows and text)
5
6args = commandArgs(trailingOnly=TRUE)
7
8table <- read.table(paste0('dades/', filename, '.txt'), sep="\t", dec='.', header=TRUE)
9
10# Plot the data
11plotall <- function() {
12 plot(table$ChannelNumber, table$Impulses, type="l", col="#EEEEFF", xlab="Número de canal", ylab="Comptes")
13 lines(table$ChannelNumber, table$ImpulsesSmooth, type="l")
14
15 # Add arrows and text labels to the peaks
16 i <- 2499 + 1
17 arrows(table$ChannelNumber[i], table$ImpulsesSmooth[i] + scale*(3 + 2), table$ChannelNumber[i], table$ImpulsesSmooth[i] + scale*0.2, length = 0.05, code = 1, col = "#514E7E")
18 text(table$ChannelNumber[i], table$ImpulsesSmooth[i] + scale*(5.2 + 2), bquote(atop(gamma["1,0"] * (Ar), 1409 ~ keV)), col = "#514E7E", cex = 0.75)
19
20 i <- 423 + 1
21 arrows(table$ChannelNumber[i], table$ImpulsesSmooth[i] + scale*3, table$ChannelNumber[i], table$ImpulsesSmooth[i] + scale*0.2, length = 0.05, code = 1, col = "#799B6F")
22 text(table$ChannelNumber[i], table$ImpulsesSmooth[i] + scale*5.2, bquote(atop(BS ~ gamma["1,0"] * (Ar), 217 ~ keV)), col = "#799B6F", cex = 0.75)
23}
24
25pdf(paste0('output/', filename, '.pdf'), width = 6, height = 5)
26plotall()
27dev.off()
28
29if (length(args) < 1 || args[1] != "preview") {
30 svg(paste0('output/', filename, '.svg'), width = 6, height = 5)
31 plotall()
32 dev.off()
33}