print('Isotopes:\n %s' % (iso_list))
print('')
element_list = list_elements(iso_list)
print('List of Elements: \n %s' % (element_list))
print('')
# scale to 1year
N_iso_ev = calc_events(iso_list, t_12, iso_perc, CdZnTe_fractions, time,
CdZnTe_molmass)
print('Total number of events after 1 year: \n %s' % (N_iso_ev))
print('')
# TOTAL NUMBER OF EVENTS Scaled TO detector mass [kg] for 1 year
N_iso_scaled = N_iso_ev * 1 / (mass_det_comp[0]) # 9 detectors
print('Total number of events after 1 year scaled to detector mass: \n %s' %
(N_iso_scaled))
print('')
### export data to a txt-file
# for saving data to a table, data must be transposed
new_data = np.array([iso_list, N_iso_scaled])
#write_txtfile(np.transpose(new_data), '../calc_solutions/', 'calculated_events.txt')
descriptions = [
'List of calculated variables',
'N_scaled: REAL NUMBER of events scaled per 1kg detector mass AND 1year'
]
var_names = ['isotope', 'p', 'm / [kg]', 'a / [Bq/kg]', 'N_scaled [1/(kg*yr)]']
write_detailed_txtfile(np.transpose(new_data), var_names, descriptions,
'./calc_solutions/', 'calculated_events_2.txt')
def my_main_2(eventfile, scale, x_range, background, Q_val_ret=None, Q_val_ret_2=None, savedir=None):
# 2) Values
all_contrib_at116Cd = []
all_contrib_at130Te = []
if not Q_val_ret_2 and not Q_val_ret:
for i_file in range(len(eventfile)):
thiscase = background+str(i_file)
data = read_evData(eventfile[i_file])
my_secHists , contrib_at116Cd, contrib_at116Cd_err, contrib_at130Te, contrib_at130Te_err = create_sector_hists(data, scale[i_file], k=thiscase, Q_val_returns=True)
sum = 0
contrib_at116Cd = unp.uarray(contrib_at116Cd, contrib_at116Cd_err)
#print(type(contrib_at116Cd))
#print('116Cd Beiträge \n', contrib_at116Cd)
for i in range(len(contrib_at116Cd)):
for j in range(len(contrib_at116Cd[i])):
sum = sum + contrib_at116Cd[i][j]
all_contrib_at116Cd.append(sum)
sum = 0
contrib_at130Te = unp.uarray(contrib_at130Te, contrib_at130Te_err)
#print('130Te Beiträge \n', contrib_at130Te)
for i in range(len(contrib_at130Te)):
for j in range(len(contrib_at130Te[i])):
sum = sum + contrib_at130Te[i][j]
all_contrib_at130Te.append(sum)
delete_all_sectorHists(my_secHists)
new_data = np.array([eventfile, all_contrib_at116Cd, all_contrib_at130Te])
descriptions = ['Contributions at Qvalues', 'N in 1/kg/keV/yr ']
var_names = ['File', 'N_at116Cd', 'N_at130Te']
if savedir:
write_detailed_txtfile(np.transpose(new_data), var_names, descriptions, './calc_solutions/%s/' %(savedir), 'events_at_Qvalues_%s.txt' %(background))
else:
write_detailed_txtfile(np.transpose(new_data), var_names, descriptions, './calc_solutions/', 'events_at_Qvalues_%s.txt' %(background))
elif Q_val_ret:
for i_file in range(len(eventfile)):
thiscase = background+str(i_file)
data = read_evData(eventfile[i_file])
my_secHists , contrib_at116Cd, contrib_at116Cd_err, contrib_at130Te, contrib_at130Te_err = create_sector_hists(data, scale[i_file], k=thiscase, Q_val_returns=True)
sum = 0
contrib_at116Cd = unp.uarray(contrib_at116Cd, contrib_at116Cd_err)
all_contrib_at116Cd.append(contrib_at116Cd)
contrib_at130Te = unp.uarray(contrib_at130Te, contrib_at130Te_err)
all_contrib_at130Te.append(contrib_at130Te)
delete_all_sectorHists(my_secHists)
return all_contrib_at116Cd, all_contrib_at130Te
elif Q_val_ret_2:
for i_file in range(len(eventfile)):
thiscase = background+str(i_file)
data = read_evData(eventfile[i_file])
my_secHists , contrib_at116Cd, contrib_at116Cd_err, contrib_at130Te, contrib_at130Te_err = create_sector_hists(data, scale[i_file], k=thiscase, Q_val_returns=True)
sum = 0
contrib_at116Cd = unp.uarray(contrib_at116Cd, contrib_at116Cd_err)
#print(type(contrib_at116Cd))
#print('116Cd Beiträge \n', contrib_at116Cd)
for i in range(len(contrib_at116Cd)):
for j in range(len(contrib_at116Cd[i])):
sum = sum + contrib_at116Cd[i][j]
all_contrib_at116Cd.append(sum)
sum = 0
contrib_at130Te = unp.uarray(contrib_at130Te, contrib_at130Te_err)
#print('130Te Beiträge \n', contrib_at130Te)
for i in range(len(contrib_at130Te)):
for j in range(len(contrib_at130Te[i])):
sum = sum + contrib_at130Te[i][j]
all_contrib_at130Te.append(sum)
delete_all_sectorHists(my_secHists)
return all_contrib_at116Cd, all_contrib_at130Te
sim = sim[len(sim) - 1]
sim = sim.split('.')
sim = sim[0]
dir = './calc_solutions/czt/%s/' % (sim)
transform2latex_tab_2(N_at116Cd_czt[i], dir,
'at116Cd_detsec_latextab.tex')
transform2latex_tab_2(N_at130Te_czt[i], dir,
'at130Te_detsec_latextab.tex')
descriptions = [
'Contributions at Qvalues', 'N in 1/kg/yr ',
'row: detector (1 to 9); columns: sector (1 to 4)'
]
write_detailed_txtfile(N_at116Cd_czt[i],
dir,
'at116Cd_detsec.txt',
var_names=None,
description=descriptions)
write_detailed_txtfile(N_at130Te_czt[i],
dir,
'at130Te_detsec.txt',
var_names=None,
description=descriptions)
sum_det_at116Cd = np.sum(N_at116Cd_czt[i],
axis=1) # creates sum of all 4 sectors
sum_det_at130Te = np.sum(N_at130Te_czt[i],
axis=1) # creates sum of all 4 sectors
new_data = [sum_det_at116Cd, sum_det_at130Te]
transform2latex_tab_2(np.transpose(new_data), dir,
'contrib_det_latextab.tex')