¿Cómo simulo el lanzamiento de una moneda sesgada en Python?

En el lanzamiento de moneda imparcial H o T ocurre el 50% de las veces. Pero quiero simular una moneda que da H con probabilidad 'p' y T con probabilidad '(1-p)'. algo como esto: def flip(p): '''this function return H with probability p''' # ...