The exclusion limits were calculated by employing the multi-bin limit setting technique in a program based on
RooStats package [1]} with the modified frequentist approach, using the profile likelihood
as a test statistic [2]}, [3]}, [4]}.
The 90% C.L. excluded region in the two-dimensional plot \(m_a - \epsilon \)  is shown in Fig. REF .
The regions excluded by the \((g-2)_e\)  measurements are also shown, the most stringent is LKB [5]}.
The central value of this measurement has the sign opposite to possible contribution from a pseudoscalar particle \(a\)  coupled to electrons.
We used a frequentist approach to calculate the 90% C.L. limit from it.
We note that the limits from the \((g-2)_e\)  measurements are model-dependent and can be significantly less strict in some scenarios
[6]}, [7]}.
<FIGURE>