import TensorScatteringClass as ten
import numpy as np
#cifpath = '/home/spc93/spc_cifs/'
cifpath = ''
ciffile=cifpath+'MnF2.cif'
en = MnK = 6.539
t=ten.TensorScatteringClass(CIFfile=ciffile, Site='Mn1');
#non-resonant magnetic scattering and resonant scattering
hkln = hkln=np.array([1,0,0])
for ref in [[0,0,3], [0,3,0], [0,2,3], [0,3,2]]:
ref = np.array(ref)
t.PlotIntensityInPolarizationChannels('NonResMag', lam=12.4/en, hkl = ref, hkln = hkln, lk=np.array([0,0,0]), sk=np.array([0,0,1]), sigmapi='sigma'); print('=== tth: %.2f deg' % (t.theta * 2 *180./np.pi))
t.PlotIntensityInPolarizationChannels('E1E1', lam=12.4/en, hkl = ref, hkln = hkln, K=2, Time=1, Parity=1, mk=None, lk=None, sk=None, sigmapi='sigma'); print('=== tth: %.2f deg' % (t.theta * 2 *180./np.pi))
t.print_tensors()