Added cone H-bond criterion

You can do this by just rearranging `angle_triplets` -- i.e. `angle_triplets2 = angle_triplets[:, [1,0,2]]` without needing to call itertools again.