Context. Forthcoming experiments will enable us to determine tomographic shear spectra at a high precision level. Most predictions about them have until now been based on algorithms yielding the expected linear and non-linear spectrum of density fluctuations. Even when simulations have been used, so-called Halofit predictions on fairly large scales have been needed. Aims. We wish to go beyond this limitation. Methods. We perform N-body and hydrodynamical simulations within a sufficiently large cosmological volume to allow a direct connection between simulations and linear spectra. While covering large length-scales, the simulation resolution is good enough to allow us to explore the high-harmonics of the cosmic shear (up to ∼ 50 000), well into the domain where baryon physics becomes important. We then compare shear spectra in the absence and in presence of various kinds of baryon physics, such as radiative cooling, star formation, and supernova feedback in the form of galactic winds. Results. We distinguish several typical properties of matter fluctuation spectra in the different simulations and test their impact on shear spectra. Conclusions. We compare our outputs with those obtainable using approximate expressions for non-linear spectra, and identify substantial discrepancies even between our results and those of purely N-body results. Our simulations and the treatment of their outputs however enable us, for the first time, to obtain shear results that are fully independent of any approximate expression, also in the high-range, where we need to incorporate a non-linear power spectrum of density perturbations and the effects of baryon physics. This will allow us to fully exploit the cosmological information contained in future high-sensitivity cosmic shear surveys, exploring the physics of cosmic shears via weak lensing measurements.