def METHOD_NAME(
    E_p: ArrayLike,
    is_12_bits_system: bool = False,
    constants: Structure = CONSTANTS_BT2020,
) -> NDArrayFloat:
    """
    Define *Recommendation ITU-R BT.2020* inverse opto-electronic transfer
    function (OETF).
    Parameters
    ----------
    E_p
        Non-linear signal :math:`E'`.
    is_12_bits_system
        *BT.709* *alpha* and *beta* constants are used if system is not 12-bit.
    constants
        *Recommendation ITU-R BT.2020* constants.
    Returns
    -------
    :class:`numpy.ndarray`
        Resulting voltage :math:`E`.
    Notes
    -----
    +------------+-----------------------+---------------+
    | **Domain** | **Scale - Reference** | **Scale - 1** |
    +============+=======================+===============+
    | ``E_p``    | [0, 1]                | [0, 1]        |
    +------------+-----------------------+---------------+
    +------------+-----------------------+---------------+
    | **Range**  | **Scale - Reference** | **Scale - 1** |
    +============+=======================+===============+
    | ``E``      | [0, 1]                | [0, 1]        |
    +------------+-----------------------+---------------+
    References
    ----------
    :cite:`InternationalTelecommunicationUnion2015h`
    Examples
    --------
    >>> oetf_inverse_BT2020(0.705515089922121)  # doctest: +ELLIPSIS
    0.4999999...
    """
    E_p = to_domain_1(E_p)
    a = constants.alpha(is_12_bits_system)
    b = constants.beta(is_12_bits_system)
    with domain_range_scale("ignore"):
        E = np.where(
            E_p < oetf_BT2020(b),
            E_p / 4.5,
            spow((E_p + (a - 1)) / a, 1 / 0.45),
        )
    return as_float(from_range_1(E))