from mymain import my_main_1, preparation
ev_dir = [
'rootfiles_noguard_06_09_2018', 'rootfiles_zcut_only_06_09_2018',
'rootfiles_guard_only_06_09_2018', 'rootfiles_guard+zcut_06_09_2018'
]
# hier in Bash Ersetzungen machen!!!
savedir = ['noguard', 'zcut_only', 'guard_only', 'guard+zcut']
eventfiles, background, scale, x_range = preparation(ev_dir[3])
background_czt = background[2]
eventfile_czt = eventfiles[2]
scale_czt = scale[2]
x_range_czt = x_range[2]
background_coating = background[1]
eventfile_coating = eventfiles[1]
scale_coating = scale[1]
x_range_coating = x_range[1]
background_plastic = background[0]
eventfile_plastic = eventfiles[0]
scale_plastic = scale[0]
x_range_plastic = x_range[0]
my_main_1(eventfile_template,
scale_template,
x_range_template,
savedir=savedir[3])
ev_dir = ['rootfiles_noguard_06_09_2018', 'rootfiles_zcut_only_06_09_2018', 'rootfiles_guard_only_06_09_2018', 'rootfiles_guard+zcut_06_09_2018' ] # hier in Bash Ersetzungen machen!!! savedir = ['noguard', 'zcut_only', 'guard_only', 'guard+zcut' ] eventfiles, background, scale, x_range = preparation(ev_dir[3]) background_czt = background[2] eventfile_czt = eventfiles[2] scale_czt = scale[2] x_range_czt = x_range[2] background_coating = background[1] eventfile_coating = eventfiles[1] scale_coating = scale[1] x_range_coating = x_range[1] background_plastic = background[0] eventfile_plastic = eventfiles[0] scale_plastic = scale[0] x_range_plastic = x_range[0] my_main_1(eventfile_coating, scale_coating, x_range_coating, savedir=savedir[3])
N_norm_coating = [] scale_coating = [] for i in range(len(iso_list)): N_norm_coating.append(sp_ac_glyptal[i] * m_detpaint_glyptal * years * n_chain[i] * 1/(n_det_glyptal* m_det)) N_norm_coating.append(sp_ac_epoxy[i] * m_detpaint_epoxy * years * n_chain[i] * 1/(n_det_epoxy* m_det)) N_simEv = 1e6 # 1 Mio simulated Events for i in range(len(N_norm_coating)): scale_coating.append(1/4 * N_norm_coating[i]/N_simEv) ########################### #### CZT ########################### datafile = './calc_solutions/calculated_events_2.txt' this_data = read_File(datafile) iso_list = this_data[:,0] # get isotope N_norm_czt = convert_str2num(this_data[:,1]) # get norming factors, number convertion necessary N_simEv = 1e6 # 1 Mio simulated Events scale_czt = [] for i in range(len(N_norm_czt)): scale_czt.append(1/4 * N_norm_czt[i] / N_simEv) # It was calculated for 9 detectors #background = ['plastic', 'coating', 'czt'] #eventfile = [eventfile_plastic, eventfile_coating, eventfile_czt] #scale = [scale_plastic, scale_coating, scale_czt] #x_range = [x_range_coating, x_range_coating, x_range_czt] my_main_1(eventfile_delrin, scale_delrin, x_range_delrin)
from mymain import my_main_1, preparation
ev_dir = [
'rootfiles_noguard_06_09_2018', 'rootfiles_zcut_only_06_09_2018',
'rootfiles_guard_only_06_09_2018', 'rootfiles_guard+zcut_06_09_2018'
]
# hier in Bash Ersetzungen machen!!!
savedir = ['noguard', 'zcut_only', 'guard_only', 'guard+zcut']
eventfiles, background, scale, x_range = preparation(ev_dir[3])
background_czt = background[2]
eventfile_czt = eventfiles[2]
scale_czt = scale[2]
x_range_czt = x_range[2]
background_coating = background[1]
eventfile_coating = eventfiles[1]
scale_coating = scale[1]
x_range_coating = x_range[1]
background_plastic = background[0]
eventfile_plastic = eventfiles[0]
scale_plastic = scale[0]
x_range_plastic = x_range[0]
my_main_1(eventfile_plastic,
scale_plastic,
x_range_plastic,
savedir=savedir[3])
ev_dir = ['rootfiles_noguard_06_09_2018', 'rootfiles_zcut_only_06_09_2018', 'rootfiles_guard_only_06_09_2018', 'rootfiles_guard+zcut_06_09_2018' ] # hier in Bash Ersetzungen machen!!! savedir = ['noguard', 'zcut_only', 'guard_only', 'guard+zcut' ] eventfiles, background, scale, x_range = preparation(ev_dir[3]) background_czt = background[2] eventfile_czt = eventfiles[2] scale_czt = scale[2] x_range_czt = x_range[2] background_coating = background[1] eventfile_coating = eventfiles[1] scale_coating = scale[1] x_range_coating = x_range[1] background_plastic = background[0] eventfile_plastic = eventfiles[0] scale_plastic = scale[0] x_range_plastic = x_range[0] my_main_1(eventfile_czt, scale_czt, x_range_czt, savedir=savedir[3])