In [1]:
import json
from Allineamento import *
In [2]:
# dizionario = {
#     "RunAllignment": 580223,
#     "d_12": 925, #cm
#     "d_1c": 925+28, #cm
#     "d_13": 925+1084, #cm
#     "offset_x2": 999,
#     "offset_y2": 999,
#     "offset_x3": 999,
#     "offset_y3": 999,
#     "div_12x": 999,
#     "div_12y": 999,
#     "div_13x": 999,
#     "div_13y": 999,
# }

# config_file = r"./config.json"

# with open(config_file, "w") as f:
#             json.dump(dizionario, f, indent=4)
In [3]:
config_file = r"./config.json"


with open(config_file, "r") as f:
            dizi = json.load(f)
nRun = dizi["RunAllignment"]

with h5py.File(f"/eos/project/i/insulab-como/testBeam/TB_2023_06_H8_GALORE/HDF5/run{nRun}.h5", 
               'r', libver='latest', swmr=True) as hf:
    xpos = np.array(hf['xpos'])

logico = (xpos > 0) & (xpos < 5)
logico2 = logico.all(axis = 1)

xpos = xpos[logico2]

offset12, divergenza12 = allineaCose(xpos[:,0], xpos[:,1], xpos[:,2], xpos[:,3], dizi["d_12"], True)
offset13, divergenza13 = allineaCose(xpos[:,0], xpos[:,1], xpos[:,4], xpos[:,5], dizi["d_13"], True)

dizionario = {
    "RunAllignment": nRun,
    "d_12": 925, #cm
    "d_1c": 925+28, #cm
    "d_13": 925+1084, #cm
    "offset_x2": offset12[0],
    "offset_y2": offset12[1],
    "offset_x3": offset13[0],
    "offset_y3": offset13[1],
    "div_12x": divergenza12[0],
    "div_12y": divergenza12[1],
    "div_13x": divergenza13[0],
    "div_13y": divergenza13[1],
}
with open(config_file, "w") as f:
            json.dump(dizionario, f, indent=4)
In [ ]: