RVSpecFit determines radial velocities and stellar atmospheric parameters from spectra by direct pixel fitting by interpolated stellar templates. The code doesn't require spectrum normalization and can deal with non-flux calibrated spectra. RVSpecFit is able to fit multiple spectra simultaneously.
https://ui.adsabs.harvard.edu/abs/2019ascl.soft07013K and https://ui.adsabs.harvard.edu/abs/2011ApJ...736..146K