*Generate a Numpy array of random host extinctions*

    **Key Arguments:**
        - ``log`` -- logger
        - ``sampleNumber`` -- the sample number, i.e. array size
        - ``extinctionType`` -- constant or random?
        - ``extinctionConstant`` -- the constant value (when extinctionType is constant)
        - ``hostExtinctionDistributions`` -- the host extinction distribution (when extinctionType is random)
        - ``plot`` -- generate plot?

    **Return:**
        - ``hostExtinctionArray``

def random_host_extinction(
        log,
        sampleNumber,
        extinctionType,
        extinctionConstant,
        hostExtinctionDistributions,
        plot=False):
    """
    *Generate a Numpy array of random host extinctions*

    **Key Arguments:**
        - ``log`` -- logger
        - ``sampleNumber`` -- the sample number, i.e. array size
        - ``extinctionType`` -- constant or random?
        - ``extinctionConstant`` -- the constant value (when extinctionType is constant)
        - ``hostExtinctionDistributions`` -- the host extinction distribution (when extinctionType is random)
        - ``plot`` -- generate plot?

    **Return:**
        - ``hostExtinctionArray``
    """
    ################ > IMPORTS ################
    ## STANDARD LIB ##
    ## THIRD PARTY ##
    import numpy as np
    ## LOCAL APPLICATION ##

    ################ > VARIABLE SETTINGS #######

    ################ >ACTION(S) ################
    # xxx come back here and add in the random extinctions generation -- will
    # need to account for what type of SN we have ##
    if extinctionType == "constant":
        hostExtinctionArray = np.zeros(sampleNumber) * extinctionConstant
    else:
        log.error('host extiction distributions not included yet')

    return hostExtinctionArray