Lab moderna: add p14 R script
Change-Id: If7c1358de4b3738f52cc4da462dae6a19fa348f5
diff --git a/quad12/labmoderna/p14/main.r b/quad12/labmoderna/p14/main.r
new file mode 100644
index 0000000..621a461
--- /dev/null
+++ b/quad12/labmoderna/p14/main.r
@@ -0,0 +1,32 @@
+file_path <- 'data.txt'
+
+xA <- 100
+xB <- 500
+
+c <- 299792458 # m/s
+h <- 6.62607015e-34 # m^2 kg s^{-1}
+
+data <- read.table(file_path, sep="\t", header=TRUE)
+data$TA <- data$VA / data$VSenseMostra
+data$TB <- data$VB / data$VSenseMostra
+tmaxA <- max(data$TA)
+tmaxB <- max(data$TB)
+data$AlphaASquared <- ((1/xA)*log(tmaxA/data$TA))^2
+data$AlphaBSqrt <- sqrt((1/xB)*log(tmaxB/data$TB))
+data$E <- h*c/(data$Longitud * 1e-9) * 6.242e18 # eV
+
+dataFitA = data[16:18,]
+dataFitB1 <- data[56:62,]
+dataFitB2 <- data[33:55,]
+
+fitA = lm(AlphaASquared~E, data = dataFitA)
+summary(fitA)
+
+fitB1 = lm(AlphaBSqrt~E, data = dataFitB1)
+summary(fitB1)
+
+fitB2 = lm(AlphaBSqrt~E, data = dataFitB2)
+summary(fitB2)
+
+plot(dataFitA$E, dataFitA$AlphaASquared)
+