def METHOD_NAME(n_events, rnd=None):
    """
    Standard inelasticity for deep inelastic scattering used so far.
    Ported from ShelfMC
    Parameters
    ----------
    n_events: int
        Number of events to be returned
    rnd: random generator object
        if None is provided, a new default random generator object is initialized
    Returns
    -------
    inelasticies: array
        Array with the inelasticities
    """
    if(rnd is None):
        rnd = np.random.default_rng()
    R1 = 0.36787944
    R2 = 0.63212056
    inelasticities = (-np.log(R1 + rnd.uniform(0., 1., n_events) * R2)) ** 2.5
    return inelasticities