getEscore function

Calculate expected scores at theta