====== Identify sample phases ======

Before you move on to multigrain indexing, you will need to know:
  * what phases you are looking for,
  * their crystal structures and cell parameters.

If you have a decent powder diffraction pattern to work with, you can try to do so with typical software people use in DAC experiments: dioptas, maud, etc.

There is a second option:
  * collect diffraction peaks using peak-search,
  * generate a 2-theta histogram of peaks.
This will be discussed below.

===== 2-theta histograms =====

Generating a 2-theta histogram of peaks identified with a [[processing:search_for_peaks|peak search]] has several advantages:
  * most of the background signal (diamonds, shadows, etc) has been removed,
  * a strong powder signal from the pressure medium will be partly removed,
  * the signal from the sample you are trying to index is stronger,
  * the histogram is based on the peak center of mass only, and hence removes the broadening due to grain size or internal strains.

Some drawbacks
  * this is not a true powder diffraction pattern, peak broadening, in particular, is different,
  * pay attention to peak intensities
    * older versions of ImageD11 created histograms of number of measured peaks -> intensity in the histogram should be proportional to the peak multiplicity
    * new version have an option //weight_hist_intensities//. If this option is set to ''True'' or ''1'', histograms are weighted by the peak intensities. These should be closer to a true powder diffraction intensity, but also affected by very large grains in the sample.

===== Creating and exporting a 2theta histogram =====

This can be done in [[software:imaged11|ImageD11]]. Once calibrated for tilts, beam center, detector distance, etc, you can select
  * Transformation -> Load Filtered Peaks
  * Transformation -> Plot tth histogram

ImageD11 will ask for a number of bins (10000 is good) and you will get something like the screenshot below. New versions of ImageD11 (as of Feb 2021) have an option //weight_hist_intensities//. If this option is set to ''True'' or ''1'', histograms are weighted by the peak intensities. If this option is set to ''False'' or ''0'', histograms are based on the number of measured peaks.

{{ :processing:2theta-histogram.png?direct&400 |}}

Newer versions of ImageD11 (update February 2021) can then export this histogram using //Transformation → Export tth histogram //

===== Phase identification and peak fitting =====

You can use [[https://github.com/seto77/PDIndexer|PDIndexer]] to identify your phases and fit their unit cell parameters.

To do so
  * Remove the first line in the file (the header),
  * Add a constant background to all intensity values. PDIndexer does not work if your diffraction intensities go to zero. This can be done in excel, for instance,
  * Save the histogram with the new intensity values with an '.xy' extension,
  * Open the XY file in PDIndexer.

===== Loading the 2theta histogram in MAUD =====

It can be useful to load the histogram in [[software:maud|MAUD]]. Remember, however, that this is not a true powder diffraction pattern, 
  * peak broadening, in particular, is different,
  * pay attention to peak intensities
    * older versions of ImageD11 created histograms of number of measured peaks -> intensity in the histogram should be proportional to the peak multiplicity
    * new version have an option //weight_hist_intensities//. If this option is set to ''True'' or ''1'', histograms are weighted by the peak intensities. These should be closer to a true powder diffraction intensity, but also affected by very large grains in the sample.

To convert the histogram, use ''timelessTthHistogram2Maud'' from the [[software:timelesstools|TIMEleSS tools]] (see details below). Simply provide the detector distance (in mm), the name of the histogram file, and the name of the new ESG file, and you are good to go.

To load the data in MAUD
  * Create a dataset with a simple diffraction instrument,
    * no intensity calibration,
    * flat image transmission, with the proper detector distance, no tilt, no beam center correction, nothing,
    * Synchrotron source with the proper wavelength,
    * Caglioti broadening, with asymmetry parameters set to 0,
  * Load the ESG file
It should work...

Remember, intensities are messed up, but you will get your phases and cell parameters...

{{ :processing:histogram-test-maud.png?direct&400 |}}

===== timelessTthHistogram2Maud =====

<code>
timelessTthHistogram2Maud  -h
usage: timelessTthHistogram2Maud -d distance histogram.dat data.esg

Creates a MAUD ESG file from a text file with intensity vs. 2theta (in
degrees). This is part of the TIMEleSS project http://timeless.texture.rocks

positional arguments:
  txtdata               Name of file with intensity vs. 2theta
  output                Name of esg file

optional arguments:
  -h, --help            show this help message and exit
  -d DISTANCE, --distance DISTANCE
                        Detector distance (in mm, required)
</code>