def METHOD_NAME(job, jobplus):
    """Calculate the difference between HOMO and IP.
    *jobplus* should be the counterpart of *job* with one less electron."""
    homo = job.results.readrkf('Properties','HOMO', file='engine')
    IP = jobplus.results.get_energy() - job.results.get_energy()
    return IP + homo