compute occupancies with sirius