I assume you are asking how LAMMPS can compute this. Not sure what
you mean by bond angle (angle between 2 bonds?). But the compute angle/local
command can calculate the theta for angles you have defined. Its output
can be used as input to fix histo/ave to histogram those values for all angles.