To further obtain the switching probability of the molecular system from one polariton state to another polariton state, we follow the recently developed global flux surface hopping (GFSH) algorithm. This algorithm is shown to outperform the original fewest switches algorithm for systems with more than two electronic states. Here, we briefly describe how to apply it to polaritonic systems. First, we express the density matrix element in the polariton basis as follows