We report Monte Carlo (MC) simulations for the pair correlation functions of a hard sphere fluid in a hard sphere matrix. In particular, we evaluate the blocking function h(b)(r) which is that part of the fluid pair correlation function which accounts for the correlation between fluid particles mediated through at least one matrix particle. The results are compared with integral equation theory. A scheme for determining the bridge functions is presented which predicts excellent agreement of the pair structure with MC results.