def photon_generation(wave_len, energy):
    """Determine number of photons in a bundle
    
    Parameters
    ----------
    wave_len : float
        Wavelength associated with a particular bundle
    energy : float
        Energy in a bundle

    Returns
    -------
    nphoton : float
        Number of photons in a bundle as a function of wavelength 

    Notes
    -----
    """
    
    nphoton = 0
    h = 6.626e-34  # [joule*s] planck's constant
    c = 2.998e8    # [m/s] speed of light
    ephoton = (h*c)/(wave_len*1e-9)  # find energy in a photon
    if wave_len <= 1107:
        nphoton = energy/ephoton
    
    return nphoton