You will need to compute the nonadiabatic couplings.
There is an implementation for the MS-PDFT method in https://github.com/pyscf/pyscf-forge/blob/master/pyscf/nac/mspdft.py, which you may find useful.
But for CCSD and CISD, we don't have them at the moment.

There is also an implementation for NACs in this repository for regular SA-CASSCF. Example examples/nac/01-sacasscf_nac.py.

Also, this issue is related to #735; though I think with different outcomes desired.

Another person recently asked here about the potential for DBOCs to be calculated in PySCF. For a (very small) list of software that can do them, see here.