In the case of DC-FSSH simulation with analytical potentials like the SBH model, energy conservation and normalization are never an issue for the following reasons. Firstly, the classical trajectories are solved with a symplectic velocity Verlet integrator, which tends to conserve total energy by construction. Secondly, after a surface hopping, we rescale the nuclear velocities to match the system's total energy before the hopping. Moreover, due to the energy-based decoherence corrections applied to the coefficients of the electronic wavefunction ansatz, we reset the wavefunction at each timestep, automatically preserving the normalization of the wavefunction.