code: def show_hbonds(self):
          
          grp = self.getPseudoBondGroup("Hydrogen Bonds-%i" % self.tid, associateWith=[self.model])
          grp.lineWidth = 3
          for i in self.plcomplex.hbonds.ldon_id:
              b = grp.newPseudoBond(self.atoms[i[0]], self.atoms[i[1]])
              b.color = self.colorbyname('blue')
              self.bs_res_ids.append(i[0])
          for i in self.plcomplex.hbonds.pdon_id:
              b = grp.newPseudoBond(self.atoms[i[0]], self.atoms[i[1]])
              b.color = self.colorbyname('blue')
              self.bs_res_ids.append(i[1])
docs: Visualizes hydrogen bonds.