<p>One of the mysteries of solar energetic particle (SEP) events is the compositional variability in those events that are clearly shock-related and may be called gradual events.<span>&#160; </span>In particular, the reason for the enhancement of Fe with respect to O or C at high energies has been debated over the past two decades, and yet it is still unsettled.<span>&#160; </span>One hypothesis relates the compositional variability with whether the CME-driven shock is quasi-parallel or quad-perpendicular near the Sun, bu<span>t this may not be easily tested using</span> remote-sensing data alone. In recent years, however, CME-driven shock waves have been modelled by fitting shock-like features in EUV and white-light images with relatively simple shapes, and in combination with magnetic field models, ir is possible to compute shock parameters at the shock surface. In this presentation, we simulate a few CMEs whose associated SEP events show widely different Fe/O, using the Alfven wave Solar Model (AWSoM) that is part of the Space Weather Modeling Framework (SWMF). We constrain the input parameters of the simulations so that the observed pre-eruption corona, erup<span>tion and CME are well-reproduced. The shock surface, across which the shock parameters are highly non-uniform, is carefully traced, and the time-dependent connectivity of the shock surface with the observer at multiple spacecraft is compared with the SEP properties including composition. We discuss how much about the compositional variability of SEP events can be learned with this technique.</span></p>