5 from ROOT
import Belle2
7 sys.exit(
"No input .root file with VXDAlignment payload specified!")
9 inputroot = sys.argv[1]
10 file = ROOT.TFile(inputroot,
"OPEN")
11 vxd = file.Get(
"VXDAlignment_CORRECTIONS")
13 errfile = ROOT.TFile(str(inputroot).replace(
'VXDAlignment_CORRECTIONS',
'VXDAlignment_ERRORS'),
"OPEN")
14 errvxd = errfile.Get(
"VXDAlignment_ERRORS")
21 shells = [ying, yang, pat, mat]
22 shellnames = [
'PXD Ying',
'PXD Yang',
'SVD Pat',
'SVD Mat']
24 params = [1, 2, 3, 4, 5, 6]
25 paramnames = [
'x',
'y',
'z',
'alpha',
'beta',
'gamma']
27 for ishell, shell
in enumerate(shells):
28 print(
'Half-shell: ', shellnames[ishell])
29 for ipar, param
in enumerate(params):
30 print(paramnames[ipar],
' = ', vxd.get(shell.getID(), param),
' +/- ', errvxd.get(shell.getID(), param))