FitAllB is tailored to do centre-of-mass (COM) refinements of grain orientations, positions and strain tensors from far-field images of a polycrystalline material. The name FitAllB originates because the routine fits all Bi (1), the grain specific reciprocal space metrics, which contain information about the strain states of the individual grains.

It was written by Jette Oddershede, formerly at DTU in Denmark and now working for Xnovo Technology ApS. The source code, along with a complete manual is available online at https://github.com/FABLE-3DXRD/FitAllB.

Jette Oddershede, Søren Schmidt, Henning Friis Poulsen, Henning Osholm Sørensen, Jonathan Wright and Walter Reimers, Determining grain resolved stresses in polycrystalline materials using three-dimensional X-ray diffraction, J. Appl. Cryst. 43, 539-549 (2010) doi: 10.1107/S0021889810012963

Taken from the manual:

In DAC experiments, FitAllB is a second stage refinement. You use FitAllB after finding and indexeding your grains with another piece of software.

fitallb.py is the script that allows you to refine grain positions and orientations and to compute the strain tensor in the sample. To make it work, you need :

• the GrainSpotter output file (.log)
• the list of filtered peaks (.flt)
• the ImageD11 parameter file of the phase (.prm)
• the structure file for your phase (.cif)
• optionally, a res file, but i don't really know what that is so I don't use it;
• the input file (.inp)

It furthermore requires PolyXSim and Iminuit. Both can be downloaded with Anaconda or is found on GitHub.

The input file should have this look :

title 'title'
log_file grains.log
flt_file peaks_t50.flt
par_file Olivine.prm
structure_file Forsterite.cif
#res_file .cif
sgno 62
#
w_step 0.5
w_limit -28 28
bg 0
dety_size 2048
detz_size 2048
#beampol_factor 0
#beampol_factor 0
#
crystal_system orthorombic
stress 0
#abs_mu 0
abs_xlim 0.1
abs_ylim 0.1
#
xyz 1 # Fit cms positions on farfield
rod 1 # Fit orientations and thus Rodrigues vectors on farfield
eps 1 # Fit strain tensors on farfield
do 0# Fit cell parameters
#
fixx 0
fixy 0
fixz 0
#rej_ia 1
rej_vol 42
rej_resmean 10
rej_resmedian 5
min_refl 7
#
tol_grain 1e-2

Now, about the parameters of the refinement, these are the four lines starting with rej_ and min_refl (maybe also the tol_grain line but it didn't do much in my tests..).

rej_ia is for internal angle, meaning the difference between the position of the measured peak and the position of the corresponding theoretical peak. In my data, I had a hard time decreasing it below 2° but I read that C. Langrand happened to have a maximum of 0.2° in internal angle misorientation.

rej_vol is linked to the intensity relativly to the volume of the grain (?). This means that the peaks have a too weak intensity. By default, it is set to 10 but I needed to increase it for my data. I succeeded in increasing up to 41 before FitAllB crashed (and it still rejected a lot of peaks because of intensity).

rej_resmean and rej_resmedian are refined quantity (as I understand it, it corresponds to how important the correction FitAllB needs to do during the refinement is) and should be as small as possible. The difference between the two is the reference they used: rej_resmedian rejects peaks that contribute to the residual by a value more than ten times the median; rej_resmean rejects peaks that contribute to the residual by a value ten times higher the mean.

min_refl is the minimum number of peaks you need to validate a grain. As it is a refinement step, during each run the grains will loose the “bad” peaks. If they go below the min_refl threshold, they are rejected.

Edit this input file to your purposes (modify file names, experimental conditions, crystallographic info, …) and then type this to your command line:

fitallb.py -i inputfilename.inp

As a result, FitAllB will run several fits in the terminal and if it succeeds to refine your grains, it will finish with normal termination of fitallb.

Your current working directory should now contain a new folder named after your input file. In this new folder you can find several output files, which FitAllB created during the run. As I understand it now, the most usefull output files are the _final.gff and the _rej files.

The _final.gff file contains the final results of the refinement, with the new positions, orientations, stress tensor, and so on of the grains. This is the one you need to keep going with the data processing. The _rej files list all the peaks rejected during the refinement and tell you why they were rejected. Most of the time it is a problem of intensity or internal angle (marked with ia).

• “maximum call reached”: Happens when FitAllB needs too much refinement loops to complete the calculation and is not powerfull enough to do it completly (?). Decrease your rej_ia, rej_vol or increase min_refl to fix it.
• “singular matrix”: Mean that fitallb does not have enough information on a grain to refine stress tensor. But same as before, decrease your rej_ia, rej_vol or increase min_refl to fix it. You can also solve this by removing the 'bad' grain in your input log file.
• “ValueError: invalid literal for int() with base 10: 'x,'” (with x = any number): This happens when you use log-files, which were produced by merging several other log-files with timelessGrainspotterMerge.py. The problem can be solved when you search the log-file for two empty lines in a row and reduce them to one. The value of x gives you the grain number where to search for this in the log-file.
• “KeyError: u'No such item: _atom_type_symbol'”: Your CIF-file is required to have both _atom_type_symbol and _atom_site_type_symbol. If one of them is missing, this error will pop up. Check other CIF-files to see how to implement these two things.
• “UnboundLocalError: local variable 'cf' referenced before assignment”: Again the CIF-file: Some old CIF-files contain _chemical_formula_sum twice (not sure why), one of them empty. Delete the empty one and try again.
• “UnboundLocalError: local variable 'unit_cell' referenced before assignment”: You probably have a typo in your crystal system (in your input file). I got rid of it, after I changed “orthorombic” to “orthorhombic”.
• “KeyError: u'p4mnm'” (or similar): It means that FitAllB cannot read the space group indicated in the CIF-file (wrong abbreviation). If you want to know, how your space group is abbreviated in a way that FitAllB can read it, check here: https://github.com/FABLE-3DXRD/xfab/blob/master/xfab/sg.py By the way: It is not case sensitive and ignores space.
• “TypeError: 'NoneType' object has no attribute 'getitem'” This error sometimes appears during the actual refinement of certain grains. Sometimes it can happen also quite late in the refinement (on grain 300 or so). So far, the only solution for this is to delete the current grain from the log-file (yes, the entire info on the grain!) and change the number of how many grains were found by 1 (first line of the log-file). Then, run FitAllB again.

fitgloball.py is a script linked to FitAllB that does the same kind of things: refining grains. I used it only for refining cell parameters.

It is working the same way as FitAllB, with the same input file, just change d0 to 1 instead of 0 to ask the script to refine the cell parameters.

A good advice to use it to refine cell parameters is to index grains with only the lowest tth (where it's easy for GrainSpotter to index peaks) and then use Fitgloball on these “perfect” grains. You can then go back to ImageD11 to modify the parameters and calculate new (and hopefully better) G-vectors, then continue with GrainSpotter, FitAllB,…