To test the performance of DC functionals constructed using the four methods above, we parametrized hybrid DC functionals based on the HCTH functional. We used the same procedure and the same multireference database as used to parameterize the rcHCTHh functional in our previous work. The multireference database consists of 59 bond energies and 60 reaction barrier heights calculated from CASSCF wave functions using an automatic active space selection scheme. The ma-TZVP basis set is used for all calculations. The CASSCF calculations were performed in Molpro or OpenMolcas. A Python module based on PySCF and Libxc 44 was used to perform MC-DCFT calculations. All MC-DCFT calculations are performed with a grid with 99 radial shells and 590 angular points per shell. The parameters in the DC functionals are optimized using a Python program based on PyTorch, Numpy, and Pandas. The accuracy of each functional is characterized by the mean unsigned error (MUE) of the predicted bond energies and reaction barrier heights. The MUE is calculated using all 119 data points, as in our previous work. Note that we used a slightly different grid as compared to our previous work. We therefore reparametrize the rcHCTHh functional using the new grid and label it as "original".