Summary: We develop a model to retrospectively evaluate age-dependent counterfactual vaccine allocation strategies against the coronavirus disease 2019 (COVID-19) pandemic. To estimate the effect of allocation on the expected severe-case incidence, we employ a simulation-assisted causal modeling approach that combines a compartmental infection-dynamics simulation, a coarse-grained causal model, and literature estimates for immunity waning. We compare Israel’s strategy, implemented in 2021, with counterfactual strategies such as no prioritization, prioritization of younger age groups, or a strict risk-ranked approach; we find that Israel’s implemented strategy was indeed highly effective. We also study the impact of increasing vaccine uptake for given age groups. Because of its modular structure, our model can easily be adapted to study future pandemics. We demonstrate this by simulating a pandemic with characteristics of the Spanish flu. Our approach helps evaluate vaccination strategies under the complex interplay of core epidemic factors, including age-dependent risk profiles, immunity waning, vaccine availability, and spreading rates. The bigger picture: Learning about a complex system and simulating alternative scenarios under changed conditions or dynamics is a challenging problem. Consider the time evolution of COVID-19 cases, which depends on a combination of contact patterns, demographics, and vaccination rates. How many severe cases could have been prevented had a different vaccine allocation strategy been implemented? To answer such counterfactual questions, we propose an approach that merges (1) coarse-grained causal modeling, (2) ordinary-differential-equation-based simulation, and (3) domain knowledge, combining the advantages of different modeling paradigms. The resulting hybrid model can be viewed as a “causal digital twin” of the underlying complex system; it captures relevant features thereof and allows reasoning about novel scenarios and interventions. We hope that our hybrid causal approach can inspire modeling for other domains where causal reasoning about a complex system is of interest.