def analyse_command_line(output_folders, **kwargs):
# logger
gt.logging_conf(**kwargs)
# Run the analysis with the command line (click)
# the return code is 0 (fail) or 1 (success)
r = analyse_all_folders(output_folders)
print(f'Last test return is: {r}')
def analyse_all_folders(output_folders, **kwargs):
# logger
gt.logging_conf(**kwargs)
# take correct folder:
outputFolders = []
for o in output_folders:
if os.path.isdir(o):
outputFolders.append(o)
# plot
ncols = 3
nrows = 3 * len(outputFolders)
fig, ax = plt.subplots(ncols=ncols, nrows=nrows, figsize=(15, 10))
r = plot_all(outputFolders, ax)
for o in range(len(outputFolders)):
fig.delaxes(ax[3 * o + 2][1])
fig.delaxes(ax[3 * o + 2][2])
plt.tight_layout(rect=[0, 0.03, 1, 0.95])
plt.savefig('output.pdf')
plt.show()
print(outputFolders)
return r
def analyse_all_folders(output_folders, **kwargs):
# logger
gt.logging_conf(**kwargs)
#Analyze folder
previous_folder = None
r = None
for folder in output_folders:
if not os.path.isdir(folder):
continue
r = analyse_one_folder(folder, previous_folder)
previous_folder = folder
# plot
ncols=3
nrows=4
fig, ax = plt.subplots(ncols=ncols, nrows=nrows, figsize=(15, 10))
i = 0
i = plot_axes(output_folders, 'dose-photon-NPV1-Edep.mhd', ax, i)
i = plot_axes(output_folders, 'dose-photon-NPV2-Edep.mhd', ax, i)
i = plot_axes(output_folders, 'dose-photon-RV1-Edep.mhd', ax, i)
i = plot_axes(output_folders, 'dose-photon-RV2-Edep.mhd', ax, i)
plt.savefig('output.pdf')
plt.show()
return r
def analyse_all_folders(output_folders, **kwargs):
# logger
gt.logging_conf(**kwargs)
# loop over folders
for output_folder in output_folders:
fn = os.path.join(output_folder, 'output.root')
f = uproot.open(fn)
r = analysis(f)
# return only last analysis
return r
def analyse_all_folders(output_folders, **kwargs):
# logger
gt.logging_conf(**kwargs)
# take correct folder:
outputFolders = []
for o in output_folders:
if os.path.isdir(o):
outputFolders.append(o)
for o in outputFolders:
r = relative_uncertainty(['reference_data', o], [1, 10], 0.05)
return r
def analyse_all_folders(output_folders, **kwargs):
# logger
gt.logging_conf(**kwargs)
#Analyze folder
previous_folder = None
r = None
for folder in output_folders:
if not os.path.isdir(folder):
continue
r = analyse_one_folder(folder, previous_folder)
previous_folder = folder
return r
def analyse_all_folders(output_folders, **kwargs):
# logger
gt.logging_conf(**kwargs)
fig, ax = plt.subplots()
is_first = True
r = None
for folder in output_folders:
r = analyse_one_folder(folder, ax, is_first)
is_first = False
fig.tight_layout()
plt.savefig(f'output.pdf')
# plt.show()
print(f'Last result : {r}')
return r
