We perform the multidimensional integration using the VEGAS Monte Carlo
method [1]}.
The centrality of collision is determined by the values of the impact parameter, \((b_{\rm min}, b_{\rm max})\) . We use an optical Glauber model [2]} to estimate geometric quantities. The electric field depends on the charge form factor. In our ”Wigner function approach” we use a charge form factor of the nucleus obtained from a realistic charge distribution described in [3]}. The cut off-forward
\(A_1 A_2 \rightarrow A_1 A_2\)  amplitude is shown in Fig.REF .
Details of helicity structure of the \(\gamma \gamma \rightarrow l^+ l^-\)  amplitude and a lot more details about the research methodology can be found in the Ref. [4]}, [5]}.
