#   ### Site Tools

##### SELF ABSORTION CORRECTION

Sometimes to subtract the sample container from the sample+container data is not right, you may have a hint that something wrong happens if, for example, after container subtraction negative intensities are obtained. That happens because neutrons are absorbed by the sample and do not scatter any more at the backside of the container. We must then calculate the coefficients Assc and Arel of the following equation:
Sample(2theta,w)=1/Assc(Ssample+container(2theta,w)-Arel*Scontainer(2theta,w)) (confer to Bee)
The calculations of these coefficients can be done by FRIDA1, and in general depend on 2theta and w. Before doing that some hints:

• To be sure that you do everything right, first you should calculate the transmission coefficient of your sample holder. To do that, we will calculate the sample holder as if it were the sample, i.e., setting the transmission coefficients of the sample holder to 1, putting the thickness of your sample holder as the thickness of the sample, and finally putting the information of the material of the sample holder (number density of scatterers, scattering cross section and absorption cross section).
• BE CAREFUL The number density of scatterers and the scattering cross section MUST BE IN THE SAME “SCATTERER UNITS”, i.e., per molecule or per atom. If you have a molecule the simplest is to calculate the number density of molecules, and then give the scattering cross section of the whole molecule (SUM Xscatt total).
• Another problem can arise when having in to account the coherent and incoherent scattering cross section. When the program is asking, it says “(coh+incoh)” but that is not really true. If your wavelength is big, you are looking (normally) at the incoherent part, therefore you should give only the incoherent part. But that is not true for bigger wavelengths, then you have incoherent plus coherent scattering. To do it correctly you should calculate which portion of your scattering is coherent and which portion is incoherent… Yes! not easy at all, but who said that neutrons experiments were easy??????
• Be careful again When you calculate the self absortion coefficients, you are getting the points of the (2theta,w) grid from the file you are starting the _sac correction from. Let's say that we do it using the sample file… then when you multiply the Arel*Scontainer you will probably get an error, because it is quite unlikely that sample and container were at the same position when measuring – therefore the points are most probably at different places. Therefore before continuing (or even before doing _sac) you should group the channels in the same grid. To do that use the mgr option.
• If your container has a transmission equal to 1, you can incorrectly think that you do not need to perform the SAC, that is (unfortunately) not true. In this case is Assc=Arel, and therefore:

Ssample(2theta,w)=1/Assc*Ssample+container(2theta,w)-Scontainer(2theta,w))
and because 0<Ascc<1, we will substract “less container”… In the following, We suppose that you have already calculated the transmission of your sample holder.

You should get a transmission of the sample of 0.89349… Now we have to make the calculation:
Ssample(2theta,w)=1/Assc(Ssample+container(2theta,w)-Arel*Scontainer(2theta,w)). We will assume that 1 you have in file Nf your sample, in file Nas your coefficient, and in Nec your emptycan. Since we have used trnasmission=1, we will not use Arel.

 Ns We will act on the sample file oy / y2 we will act con y axis, and divide by an y of a second file Nas We will divide by the Assc coefficient * We have now in file No1 the divission, and we should now substract the empty can No1 We will act in file No1 (Ssample/Assc) oy - y2 We will substract the y of another file Nec The empty can is in file Nec Nresult That is the result… Please do compare it with a simple substraction to get a feeling of what we have done
##### transmission

To calculate the transmission of the sample container… Please do not use without reading the introduction, the values are for our case, and are absolutely not correct in general

 _sac to call the Self Absortion Correction function. You can call it form any file. The program will simply make the calculus in the (2theta,w) space of your file) 0 We will suppose that the whole sample was iluminated 3 we will calculate it form a hollow cylinder 1.2 outer radius in cm of your sample holder 0.1 thickness of the sample layer, which in this case is the thickness of your sample holder 1 1 transmission of the outer container set to one for the same reason 0 Fraction of back-scattered beam transversing 60.26 Number density of scatterers in units of 10 h 21 cm-3 , in our case was Aluminium… in fact that changes with T… that should not be so important… 1.734 Scattering cross section of Aluminium in barns 0.231 Absortion cross section of Aluminum for wavelength=1.798 in barns, Thermal neutrons…

Then read what FRIDA says, transmission = 0.999857659876, and that is the transmission of your sample holder.

##### mgr

We will now try to group the channels in the same grid for two files…

 mgr to strat grouping the files 1 we will do it in a linear grid, and we will need to give the values of the 2theta,w grid 0 from this energy 0.02 in steps of 0.02 meV 200 to 200 meV. An error message will appear if too many points are chosen 3 You must decide wether to or to ***

You may get an error message if the number of points is too big. If you want to study the inelastic part of the spectra is better to perform the grid in log scale…

… or 