Averaged XAS spectrum using site_weighted_spectrum in pymatgen.analysis.xas?

I have obtained several XAS/EELS spectra from the same type of elements in amorphous structures. We need site averaged spectra so that I compare it with the experimentally observed spectra.

Is it possible to get a superposition of several XAS/EELS spectra using site_weighted_spectrum in pymatgen.analysis.xas?

The baseline of the simulation spectra are adjusted and the data points of E(eV), which is the x-axis, are not aligned among the spectra. So, data manipulation such as interporation between the two points is required during the calculation for actual use.