In the context of SARS-CoV-2 pandemic, mathematical modelling has played a funda-mental role for making forecasts, simulating scenarios and evaluating the impact of pre-ventive political, social and pharmaceutical measures. Optimal control theory represents a useful mathematical tool to plan the vaccination campaign aimed at eradicating the pandemic as fast as possible. The aim of this work is to explore the optimal prioritisation order for planning vaccination campaigns able to achieve specific goals, as the reduction of the amount of infected, deceased and hospitalized in a given time frame, among age classes. For this purpose, we introduce an age stratified SIR-like epidemic compartmental model settled in an abstract framework for modelling two-doses vaccination campaigns and conceived with the description of COVID19 disease. Compared to other recent works, our model incorporates all stages of the COVID-19 disease, including death or recovery, without accounting for additional specific compartments that would increase computa-tional complexity and that are not relevant for our purposes. Moreover, we introduce an optimal control framework where the model is the state problem while the vaccine doses administered are the control variables. An extensive campaign of numerical tests, featured in the Italian scenario and calibrated on available data from Dipartimento di Protezione Civile Italiana, proves that the presented framework can be a valuable tool to support the planning of vaccination campaigns. Indeed, in each considered scenario, our optimization framework guarantees noticeable improvements in terms of reducing deceased, infected or hospitalized individuals with respect to the baseline vaccination policy.& COPY; 2023 The Authors. Publishing services by Elsevier B.V. on behalf of KeAi Communications Co. Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).