We discuss the generation of the directed flow $$v_1(p_T,y_z)$$
v
1
(
p
T
,
y
z
)
induced by the electromagnetic field as a function of $$p_T$$
p
T
and $$y_z$$
y
z
. Despite the complex dynamics of charged particles due to strong interactions generating several anisotropies in the azimuthal angle, it is possible at $$p_T > m$$
p
T
>
m
to directly correlate the splitting in $$v_1$$
v
1
of heavy quarks with different charges to some main features of the magnetic field, and in particular its values at formation and freeze-out time. We further found that the slope of the splitting $$d\Delta v_1/dy_z|_{y_z=0}$$
d
Δ
v
1
/
d
y
z
|
y
z
=
0
of positively and negatively charged particles at high $$p_T$$
p
T
can be formulated as $$d\Delta v_1/dy_z|_{y_z=0}=-\alpha \frac{\partial \ln f}{\partial p_T}+\frac{2\alpha -\beta }{p_T}$$
d
Δ
v
1
/
d
y
z
|
y
z
=
0
=
-
α
∂
ln
f
∂
p
T
+
2
α
-
β
p
T
, where f is the $$p_T$$
p
T
spectra of the charged particles and the constants $$\alpha $$
α
and $$\beta $$
β
(order of MeV) are constrained by the y component of magnetic fields and the sign of $$\alpha $$
α
is simply determined by the difference $$\Delta [tB_y(t)]$$
Δ
[
t
B
y
(
t
)
]
in the center of colliding systems at the formation time of particles and at the time when particles leave the effective range of electromagnetic fields or freeze out. The formula is derived from general considerations and is confirmed by several related numerical simulations; it supplies a useful guide to quantify the effect of different magnetic field configurations and provides an evidence of why the measurement of $$\Delta v_1$$
Δ
v
1
of charm, bottom and leptons from $$Z^0$$
Z
0
decay and their correlations are a powerful probe of the initial e.m. fields in ultra-relativistic collisions.