We present a classical treatment of the k-r-k' vector correlation for atom-diatom reactions, involving the directions of the initial and final relative velocities and the internuclear axis of the diatom. The formalism is based on the expansion of the joint probability distributions in multipolar moments, and it is analogous to that developed for the vector correlations implying the initial or final rotational angular momentum. Within the framework of classical mechanics, the present treatment allows the determination of the differential cross section for any preparation of the internuclear axis and for any initial rotational angular momentum of the reactants. This methodology has been applied to the study of the steric effect in the Li+HF(v = 1,j = 1,m = 0)-->LiF+H reaction and the theoretical results have been compared via simulation with the recent experimental determination by Loesch and co-workers of laboratory angular distributions (LAB ADs) for several distributions of the internuclear axis of HF. Very good agreement has been found between experimental and simulated LAB AD, allowing an interpretation of the experimental results. Although the integral steric effect is very small, the differential and/or state resolved steric effects are more pronounced. Moreover, it has been found that by varying the preparation of the internuclear axis of the HF molecule, the population of final states of the LiF changes noticeably, which represents a clear case of control of the outcome of a chemical reaction.