How can I calculate hydrogen bonds in a specific direction w.r.t. a C-alpha atom in a protein?

I need to calculate H-bonds in a specific direction of a C-alpha atom of a protein. And, I need to calculate that from a PDB file.

Can anyone give me a general direction/idea regarding how to do that?

N.B. I need to write a python program.