COMPONENT sample = PowderN(radius=0.01,h=0.01, d_phi0=0.1, pack=0.5, DW=0.9, frac=0.5, focus_r=0.03, reflections="mylist.dat", Vc=3.86*3.86*11.82, sigma_abs=0, sigma_inc=2, barns=1) AT (0, 0, 0) RELATIVE a4You should also create a list of reflections, e.g. mylist.dat with the following contents:
# column_j 3 multiplicity 'j' # column_q 1 Scattering vector modulus [Angs^-1] # column_F2 2 Scattering factor |F^2| in [barns] 1 1000 8 1.3 1000 4
(Adjust placement relative to your component naming etc.). Also insert a banana shaped detector,
COMPONENT ND = Monitor_nD(xwidth=1, options="banana, theta limits=[-45 45],auto y") AT (0,0,0) RELATIVE a5and take away the previously inserted Beamstop. Test that neutrons reach the detector, then run a long simulation (many neutron rays) and have coffee.