Hello Cristina,
if you load each alignment to a dataframe you can use the biggest as a reference. You can iterate over the others row-by-row and test each mz/rt pair for matches against every feature in the reference. When you obtain matches you could aggregate peak areas from both dfs in a new dataframe. This is fairly straightforward in Python.
Chris