Graphene Saturable Absorption (Multilayered System)

Hello,

Does anyone know how to model the graphene saturable absorption on a multilayered system that already calculates analytically the reflectance from the transfer matrix model?

Thanks in advance,