Add hbonds

import molsysmt as msm
import numpy as np
from molsysmt import pyunitwizard as puw

molsys = msm.convert(msm.systems['Barnase-Barstar']['barnase_barstar.h5msm'])

hbonds, distance, angles = msm.hbonds.get_luzard_chandler_hbonds(molsys, selection='molecule_name=="BARNASE"',
                                              selection_2='molecule_name=="BARSTAR"')

view = msm.view(molsys)

msm.thirds.nglview.add_hbonds(view, hbonds[0])

view

view._ngl_msg_archive[-len(hbonds[0]):]

view = msm.view(molsys)
msm.thirds.nglview.add_hbonds(view, hbonds, hbond_level='group')
view

view._ngl_msg_archive[-1]