We introduce a global analysis of relativistic heavy ion collisions using Trajectum of a significantly higher precision and including a new option to vary the normalization of the centrality estimator. We use the posterior distribution of our parameters to generate a set of high statistics samples that allows us to make precise predictions including statistical and systematic uncertainties estimated from the model parameter distribution. The results are systematically compared with experiment whereby we also include many observables not included in the global analysis. This includes in particular (extremely) ultracentral anisotropic flow and mean transverse momentum, whereby we find satisfactory agreement with experiment where data is available. Lastly, we compute spectra and anisotropic flow for oxygen-oxygen collisions performed at RHIC and to be performed at the LHC and comment on how these collisions may inform us on properties of the Quark-Gluon-Plasma. CONTENTS I. Introduction 1 II. Model 2 A. Oxygen nuclei 3 B. Continuous number of constituents 3 C. UrQMD and SMASH 3 D. Centrality determination 3 III. Precise PbPb results and OO predictions 4 A. The design and emulator validation 4 B. Posterior distribution 6 C. Postdictions with systematic uncertainties 8 D. p T -differential predictions 10 IV. Predictions for ultracentral collisions 13 V. Bayesian analysis using OO simulated data 14 VI. Discussion 16 References 18