An eigenmode solver based on particle-in-cell simulations is proposed to find the hybrid beam-electromagnetic eigenmodes in a "hot" slow-wave structure (SWS), i.e., where the electromagnetic mode interacts with an electron beam. The proposed scheme is based on the determination of the transfer matrix of the unit-cell of the hot SWS that takes into account the interaction of the electromagnetic mode with the electron beam, followed by the determination of its eigenvalues and eigenvectors using Floquet theory. The proposed method is applied to find the hot modes, with complex wavenumber, in a traveling wave tube (TWT) amplifier with helix SWS. We show dispersion relations of the modal complex wavenumber of the electromagnetic wave interacting with the electron beam when varying frequency and beam voltage, with results in agreement with Pierce theory. The method is also applied to find the complex-wavenumber modes of the hot SWS of a millimeter wave TWT amplifier based on a serpentine waveguide. The technique is general and can be applied to any SWS geometry where the electromagnetic modes interact with an electron beam.