Reading TTrees for DC2

[brunollago]

Reading trees as suggested on the script DC2Alpha/PlotEventTraces.py by Matias.

Simulation directory: GP300_Xi_Sib_Proton_1.25_53.9_60.52_13020/

*** Issues:

1 - The site height does not seem to be available on the root trees. It would be nice to have it in meters somewhere.

2 - gt.tshower.xmax_pos is missing. gt.tshower.xmax_pos_shc is available though.

3 - gt.tshower.shower_core_pos seems to be wrong
print(gt.tshower.shower_core_pos)
[-3541.518 5465.649 0. ]

3a - gt.tshower.direction is missing missing. I was able to evaluate it using
my_core = [0, 0, 1264]
ev = my_core - gt.tshower.xmax_pos_shc
print(gt.tshower.zenith, np.arccos(-ev[2]/np.linalg.norm(ev))*180/np.pi)
53.89 53.88996881224265

4 - du_xyz is supposed to be in site's referential
print(gt.trun.du_xyz[0])
[-788.5800170898438, -1965.6500244140625, 1336.5]
Is the output compatible with site's referential? Shouldn't the site height be subtracted here?

[mjtueros]

1- The site height is available on the geoid
2- This is expected behaviour, we still dont have all the variables filled up.
3- i dont see anythng wrong with the core position...its in meters, in the site reference frame.
3a (see 2)
4- du_xyz might be in the shower referential. This needs checking

I leave the issue open for points 2,3a and 4