55def linear_fit_in_zx_plane(muon_tagger_track, drifttubes_track, debug):
56
57
58 z_coords = []
59 x_coords = []
60
61 for ahit in muon_tagger_track['hits_x']:
62 z_coords.append([ahit['z']])
63 x_coords.append(ahit['xtop'])
64
65 for ahit in drifttubes_track['34']:
66 z_coords.append([ahit['z']])
67 x_coords.append(ahit['xtop'])
68
69 if debug:
70 print('z_coords:', z_coords)
71 print('x_coords:', x_coords)
72
73
74 reg = LinearRegression()
75 reg.fit(z_coords, x_coords)
76
77
78 x_coords_pred = reg.predict(z_coords)
79
80
81 rmse = np.sqrt(mean_squared_error(x_coords, x_coords_pred))
82
83
84 return rmse
85
86
87
88
89