%load_ext autoreload
%autoreload 2

import molsysmt as msm

Get lengths and angles from box

molsys = msm.convert(msm.systems['Met-enkephalin']['met_enkephalin.h5msm'])

molsys_cub = msm.build.solvate(molsys, box_shape='cubic', clearance='14.0 angstroms', engine='PDBFixer')
molsys_oct = msm.build.solvate(molsys, box_shape='truncated octahedral',  clearance='14.0 angstroms', engine='PDBFixer')
molsys_dod = msm.build.solvate(molsys, box_shape='rhombic dodecahedral', clearance='14.0 angstroms', engine='PDBFixer')

box = msm.get(molsys_cub, element='system', box=True)
msm.pbc.get_lengths_and_angles_from_box(box)

(<Quantity([[3.1236 3.1236 3.1236]], 'nanometer')>,
 <Quantity([[89.99998128 89.99998128 89.99998128]], 'degree')>)

box = msm.get(molsys_oct, element='system', box=True)
msm.pbc.get_lengths_and_angles_from_box(box)

(<Quantity([[3.1236 3.1236 3.1236]], 'nanometer')>,
 <Quantity([[ 70.52875545 109.4712071   70.52875545]], 'degree')>)

box = msm.get(molsys_dod, element='system', box=True)
msm.pbc.get_lengths_and_angles_from_box(box)

(<Quantity([[3.1236 3.1236 3.1236]], 'nanometer')>,
 <Quantity([[60.00002571 60.00002571 89.99998128]], 'degree')>)