  The orientational correlation scheme introduced earlier for tetrahedral
molecules is extended for being able to classify orientational correlations
between pairs of high symmetry molecules. It is based on the number of corner
atoms of each molecule found within the orientation cones and provides the
probability of 6 orientation groups versus the centre-centre distance. The apex
of the orientation cone coincides with the molecular centre and the axis points
towards the other centre, while the apex angle is chosen, that the average
coordination number of the corner atoms inside would be 2. The method is
applicable to molecules whose corner atoms are spaced evenly and only 1, 2, or
3 corner atoms must be found inside the cone. To demonstrate the applicability
of the method, the plastic crystalline phase of the truncated icosahedral
shaped C$_{60}$ molecules and the octahedral-shaped SF$_6$ molecules in its
gaseous, supercritical fluid, liquid and plastic crystalline phases have been
modeled using classical molecular dynamics. For SF$_6$, after the best
performing forcefield, in terms of the density and the total scattering
structure factor, was found, orientational correlations have been analyzed. In
the non-crystalline phases, close contact orientation and medium-range order
regions could be identified. While the former is invariant to the changes of
the density, the latter showed longer-ranged correlations as density raised. In
the plastic crystalline state, fluorine atoms are oriented along the lattice
directions with higher probability. In C$_{60}$, local order in the first and
second coordination shells is determined quantitatively for three force fields.
