Simulation of near perfect crystals
Hi Dr. Smid,
I have a quick question regarding a simulation using a flat, near perfect crystal (in this case ADP). When trying to assume an idealised perfect crystal by setting mosaicity to 0, I continually get divide by zero errors when using the default rockingCurveFWHM and integrated_reflectivity. By setting the rockingCurveFWHM to 1 I can get results, but they don't look particularly promising.
How would you simulate this situation, assuming the crystal can be taken as perfect (where resolution is only limited by the detector resolution)?
Attached you'll find the code I'm using.ADP_Simulation.py
Thanks and best regards, Carlos Butler