# current * live time for each file
it1 = 1000 * 959
it2 = 1000 * 960
it3 = 1000 * 959
# all lists above indexed by file names
peaks = {}
errors = {}
irradiation_parameters = {}
# parser files
keys = csvs.keys() # or txts.keys()
for key in keys:
irradiation_parameters[key] = parseCsv(csvs[key])
txt_content = parseTxt(txts[key])
peaks[key] = txt_content['K']['peaks']
errors[key] = txt_content['K']['errors']
class blankCorrectionTest(unittest.TestCase):
def test_13(self):
# teste peak
testcase = blankCorrection(irradiation_parameters, peaks,
errors)['peaks_correction'][13]
calculed = (0 / it1 + 226 / it2 + 212 / it3) / 3
self.assertAlmostEqual(testcase, calculed)
# test error
testcase = blankCorrection(irradiation_parameters, peaks,
errors)['errors_correction'][13]
import math
sys.path.append('./lib/')
from responseFactor import responseFactor
from shimadzu import parseCsv
from winqxas import parseTxt
import numpy as np
# test data
micromatter = pathlib.Path(
'data/calibration/micromatter-table-iag.csv').read_text()
csv_34671 = pathlib.Path('data/calibration/2010-10/csv/34671.csv').read_text()
txt_34671 = pathlib.Path('data/calibration/2010-10/txt/34671.txt').read_text()
irradiation_parameters = parseCsv(csv_34671)
txt_content = parseTxt(txt_34671)
peaks = txt_content['K']['peaks']
errors = txt_content['K']['errors']
i = irradiation_parameters['current']
t = irradiation_parameters['livetime']
N = peaks[22]
sigma_N = errors[22]
class calculateResponseFactorTest(unittest.TestCase):
def test_Ti(self):
R, sigma_R = responseFactor(N, 49.4, i, t, sigma_N)
R_calculated = 454712 / (268 * 179 * 49.4)
sigma_calculated = R_calculated * np.sqrt((676 / 454712)**2 +
(0.05)**2 + (0.05)**2 +
d = df['density1'][i]
# ler contagens do arquivos txt
import sys
import pathlib
sys.path.append('elemental_analysis_tools')
from winqxas import parseTxt
#transformar nome do documento em str
serial = str(df['serial'][i])
doc_txt = 'data/calibration/2010-10/txt/'
for j in range(len(serial)):
doc_txt += serial[j]
for j in range(len(pto_txt)):
doc_txt += pto_txt[j]
file_content = pathlib.Path(doc_txt).read_text()
txt = parseTxt(file_content)
try:
N = txt['K']['peaks'][df['element1'][i]]
sigma_N = txt['K']['errors'][df['element1'][i]]
except:
pass
#print(txt['K']['errors'][11])
# ler corrente e tempo dos arquivos csv
import shimadzu as sd
doc_csv = 'data/calibration/2010-10/csv/'
for j in range(len(serial)):
doc_csv += serial[j]
for j in range(len(pto_csv)):
doc_csv += pto_csv[j]
file_content = pathlib.Path(doc_csv).read_text()
def test_56(self):
self.assertEqual(parseTxt(file_content)['L']['peaks'][56], 1182)
self.assertEqual(parseTxt(file_content)['L']['errors'][56], 214)
def test_42(self):
self.assertEqual(parseTxt(file2_content)['K']['peaks'][42], 372924)
self.assertEqual(parseTxt(file2_content)['L']['peaks'][42], 30490)
self.assertEqual(parseTxt(file2_content)['K']['errors'][42], 625)
self.assertEqual(parseTxt(file2_content)['L']['errors'][42], 293)
def test_22(self):
self.assertEqual(parseTxt(file_content)['K']['peaks'][22], 19355)
self.assertEqual(parseTxt(file_content)['K']['errors'][22], 222)
def test_14(self):
self.assertEqual(parseTxt(file_content)['K']['peaks'][14], 38268)
self.assertEqual(parseTxt(file_content)['K']['errors'][14], 282)
def test_13(self):
self.assertEqual(parseTxt(file_content)['K']['peaks'][13], 8871)
self.assertEqual(parseTxt(file_content)['K']['errors'][13], 133)
def test_11(self):
self.assertEqual(parseTxt(file_content)['K']['peaks'][11], 389)
self.assertEqual(parseTxt(file_content)['K']['errors'][11], 71)
import pathlib
import sys
import math
sys.path.append('./eas')
from elementarDensity import elementarDensity
from shimadzu import parseCsv
from winqxas import parseTxt
# test data
csv = pathlib.Path('data/ghana/EDX720/AFR390.20110406223141.csv').read_text()
txt1 = pathlib.Path(
'data/ghana/winqxas/txt1/[1]AFR39020110406223141.txt').read_text()
irradiation_parameters = parseCsv(csv)
txt1_content = parseTxt(txt1)
peaks1 = txt1_content['K']['peaks']
errors1 = txt1_content['K']['errors']
i = irradiation_parameters['current']
t = irradiation_parameters['livetime']
N = float(peaks1[11])
class calculateResponseFactorTest(unittest.TestCase):
def test_11(self):
testcase = elementarDensity(N, 0.00034021562, i, t)
calculated = N / (0.00034021562 * i * t)
self.assertAlmostEqual(testcase, calculated)