User Tools

Site Tools


examples:gs_loop_olivine

This is an old revision of the document!


GrainSpotter loop for olivine

Below is an example of a GrainSpotter loop for olivine. You can run it by typing

gsLoop.sh > gsLoop.log

The output of the loop will be saved in gsLoop.log.

Main idea

The loop can run nLoops indexings (typically 3). For each, the loop runs

  • 6 iterations of 100000 GrainSpotter orientation searches with the following threshold
  • minNPeaks = 20
  • minCompleteness = 0.5
  • Δ = 0.06°
  • Δη = 1°
  • Δω = 1°
  • maximum deviation nσ = 2
  • 5 iterations of 100000 GrainSpotter orientation searches with
  • minNPeaks = 20
  • minCompleteness = 0.5
  • Δ = 0.06°
  • Δη = 2°
  • Δω = 2°
  • maximum deviation nσ = 2
  • 5 iterations of 100000 GrainSpotter orientation searches with
  • minNPeaks = 20
  • minCompleteness = 0.3
  • Δ = 0.06°
  • Δη = 2°
  • Δω = 2°
  • maximum deviation nσ = 2
  • 5 iterations of 100000 GrainSpotter orientation searches with
  • minNPeaks = 20
  • minCompleteness = 0.3
  • Δ = 0.06°
  • Δη = 3°
  • Δω = 3°
  • maximum deviation nσ = 2
  • 5 iterations of 100000 GrainSpotter orientation searches with
  • minNPeaks = 20
  • minCompleteness = 0.3
  • Δ = 0.10°
  • Δη = 3°
  • Δω = 3°
  • maximum deviation nσ = 2

At the end of the loop, all results for the corresponding loop number are saved in a folder called Indexing-$loop, where $loop is the loop number. Merged grains are in a file called merge-grains.log and their Euler angles in merge-grains-euler.txt.

Getting some idea on indexing capabilities

A good start is to look for the number of assigned G-vectors. This command will do it

cat gsLoop.log | grep "Assigned gvectors"

It will print out something like

ssigned gvectors: 3915
Assigned gvectors: 3719
Assigned gvectors: 3719
Assigned gvectors: 3719
Assigned gvectors: 3719
Assigned gvectors: 3719
Assigned gvectors: 3719
Assigned gvectors: 3139
Assigned gvectors: 2997
Assigned gvectors: 2972
Assigned gvectors: 2950
Assigned gvectors: 2927
Assigned gvectors: 2873
Assigned gvectors: 2873
Assigned gvectors: 2873
Assigned gvectors: 2873
Assigned gvectors: 2873
Assigned gvectors: 1831
Assigned gvectors: 1831
Assigned gvectors: 1831
Assigned gvectors: 1817
Assigned gvectors: 2121
Assigned gvectors: 1454
Assigned gvectors: 1454
Assigned gvectors: 1454
Assigned gvectors: 1454
Assigned gvectors: 3915

3915 is the number of experimental G-vectors that could match our crystal structure. It decreases it time we generate a new GVE file based on a previous GrainSpotter indexing. When it increases back to 3915, it means that the loop fully started over. In the example above, we started from 3915 experimental G-Vectors, 1454 are left, we assigned approximately 60% of the experiment G-Vectors.

The ratio of assigned G-Vectors is not 100% trivial to calculate as the number of matching experimental G-vectors depends on your Δ and nσ parameters in the GrainSpotter input file.

examples/gs_loop_olivine.1575543091.txt.gz · Last modified: 2019/12/05 10:51 by smerkel