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.

Steve