self.left_z, self.left_r = R_h_left * np.sin(t_edges_left[:-1])/(1.*n_h_left),\
                R_h_left * np.cos(t_edges_left[:-1])/(1.*n_h_left)