def calc_energy_position(energy_edges, position_rate=0.5):
    """Calculate energies at which position_rate corresponds to."""
    out = []
    for i in range(len(energy_edges)-1):
        width = energy_edges[i+1] - energy_edges[i]
        out.append(energy_edges[i] + position_rate*width)
    return out