Supernova remnants often evolve in material with high abundance of heavy elements such as carbon or oxygen. Hadronic collisions in these enriched media spawn the production of secondary particles such as gamma rays, neutrinos, and secondary electrons with spectra that cannot be scaled from those calculated for pp collisions, potentially leading to erroneous results. We used Monte-Carlo event generators to calculate the differential production rate of secondary particles such as gamma rays, neutrinos, and secondary electrons for H, He, C, and O nuclei as projectiles and as target material. The cross sections and the multiplicity spectra are separately computed for each of the 16 combinations of projectile and target. We describe characteristic effects in the form and normalization of spectra.