code: def angle(self, center1_x, center1_y, center2_x, center2_y):
          
          phi_G = np.arctan2(center2_y - center1_y, center2_x - center1_x)
          return phi_G
docs: compute the rotation angle of the dipole
          :return: