Abstract : The VAMOS spectrometer operational at GANIL is a large acceptance variable mode spectrometer designed for nuclear reaction studies using radioactive and stable ion beams. The spectrometer coupled with ancillary detectors like EXOGAM has been successfully used in recent experiments on (in)elastic, few nucleon transfer reactions in inverse kinematics and search for nuclei far from stability using deep inelastic transfer reactions In large acceptance spectrometers, the experimental resolution is worsened by aberrations of third and higher orders. Hardware corrections are limited and cannot completely correct the non-linear effects. Thus software reconstruction of trajectories (ray tracing) is essential to obtain the required resolution and identification of the products. A numerical method has been developed for reconstruction of ion trajectories and correction of aberrations in VAMOS. We have devised a procedure to select an optimum subset of closest trajectories for each focal plane event (x, θ, y, ϕ) from the database (generated by an ion-optics calculation). A polynomial fit to the momentum vector of the reaction product in terms of (x, θ, y, ϕ) is made only to this subset. Such an approach is found to give improved resolutions compared to fitting a single polynomial over the entire phase space. Extraction of charge state and angular distributions are rather difficult due to the variation of acceptance over the large phase space. Mass identification of the reaction products and characteristics of the spectrometer acceptance with its variation for different rigidities have been obtained. Applications to 238U+48Ca and 238U+58Ni systems at 5.5 MeV/u will be presented.