Ultra-strong electromagnetic field in relativistic heavy-ion collisions  could be quantized as a large flux of linearly polarized quasi-real photons. These photons can interact with nuclei or nucleons, leading to vector meson production ($\gamma + A \rightarrow V + A$). Notably, azimuthal asymmetries between the pair momentum and the daughter momentum could arise from the linear polarization of incident photons correlated with spin interference effect. On the other hand, the decay daughters of these vector mesons inherit polarization information of the photons, which is related to the initial collision geometry.  Thus, the measurement of azimuthal anisotropy of these decay daughters offers a novel and direct probe into both the initial collision geometry and the polarization characteristics of the photons.
 
 In this presentation,  we will report the angular modulation measurements of the photon-induced $J/\psi$  pair production in  Ru+Ru and Zr+Zr ultra peripheral collisions at $\sqrt{s_{NN}}$ = 200 GeV. Additionally, we will also present the measurements of the azimuthal anisotropy with respect to the event plane ($v_2$) of electrons decayed from photon-induced $J/\psi$ in non-central collisions from the same collision systems. 
 Our results offer novel insights into both the collision geometry and spin interference effect.