Everything up to the [spoiler ahead] in this comment is (as far as I can tell) exactly how things work in standard formulations of Bell's inequality. There's nothing weird or crackpot there.
Your numerical code is impossible for me to read without some basic idea of what you're trying to show, but I'd like to point out that numpy has functions like np.radians, and np.deg2rad to convert from degrees to radians, you don't have to make your own.
Your numerical code is impossible for me to read without some basic idea of what you're trying to show, but I'd like to point out that numpy has functions like np.radians, and np.deg2rad to convert from degrees to radians, you don't have to make your own.