| 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) |
| |