differences between the observed positions (x) and the model positions
epsilon = epsilonInRange(self,zmin,zmax)
zmin
zmax
epsilon