electrons_in_range = ma.masked_where(
        np.absolute(energy_CB - Data.CBM) > extrema_search_depth, energy_CB)