Starting from Post-Newtonian predictions for a system of N infalling masses from the infinite past, we formulate and solve a scattering problem for the system of linearised gravity around Schwarzschild as introduced in [DHR19]. The scattering data are posed on a null hypersurface C emanating from a section of past null infinity I−, and on the part of I− that lies to the future of this section: Along C, we implement the Post-Newtonian theory-inspired hypothesis that the gauge-invariant components of the Weyl tensor α and α_ (a.k.a. Ψ0 and Ψ4) decay like r−3, r−4, respectively, and we exclude incoming radiation from I− by demanding the News function to vanish along I−. We also show that compactly supported gravitational perturbations along I− induce very similar data, with α, α_ decaying like r−3, r−5 along C. After constructing the unique solution to this scattering problem, we provide a complete analysis of the asymptotic behaviour of projections onto fixed spherical harmonic number ℓ near spacelike i0 and future null infinity I+. Using our results, we also give constructive corrections to popular historical notions of asymptotic flatness such as Bondi coordinates or asymptotic simplicity. In particular, confirming earlier heuristics due to Damour and Christodoulou, we find that the peeling property is violated both near I− and near I+, with e.g. α near I+ only decaying like r−4 instead of r−5. We also find that the resulting solution decays slower towards i0 than often assumed, with α decaying like r−3 towards i0. The issue of summing up the fixed angular mode estimates in ℓ is dealt with in forthcoming work.