Need help to separate slip system by RDR value.

[Patcurran]

Hi everyone,

Happy new year!

My aim with this work is to plot on an IPF triangle, if a grain has been activated (or not), and if it has been activated which slip system has been activated, for a HCP system.

I have been given some code that plot on an IPF triangle if activated, based on slip plane (shown below), i.e. basal, prismatic, pyramidal but i have both pyramidal < a > and pyramidal < a+c >. So i need to use the RDR values to distinguish which system has been activated.

I checked the DicMap_RDR.txt file and the RDR values are not saved.
I then looked at how RDR values were saved, which i found for the experimentally measured values in DefDAP/defdap/inspector.py line 559. Although, i could not find how the theoretical RDR values are calculated from the EBSD data.

Can anyone help me with any of the following:

1. How to separate activated slip systems by RDR values?
2. Where i can find the code to calculate and save the theoretical RDR values (from the EBSD data)?
3. Explain to me how the theoretical RDR values are calculated and save.

any other helpful comments are appreciated :)

best wishes,
Patrick

code:

prorilist = []; pyorilist = []; baorilist = []; multorilist = [];  allgrains=[]; 
for idx, grain in enumerate(DicMap):
   try:
       system = grain.groupsList[0][2]
   except IndexError:
       continue
       
   allgrains.append(grain.ebsdGrain.refOri)
         
   if len(system)>1 or system==[]:
       i=1
   else: 
       if system == [0]:
           baorilist.append(grain.ebsdGrain.refOri)
       if system == [1] or system == [2] or system == [3]:
           prorilist.append(grain.ebsdGrain.refOri)
       if system == [4] or system == [5] or system == [6] or system == [7] or system == [8] or system == [9]:
           pyorilist.append(grain.ebsdGrain.refOri)

f, ax = plt.subplots(figsize=(4,3))
markersize=40
quat.Quat.plotIPF(allgrains, [1, 0, 0], 'hexagonal', c='b', marker='o', s=markersize, ax=ax,label='All grains')  
quat.Quat.plotIPF(allgrains, [1,0,0], 'hexagonal', c='0.9', marker='o', s=markersize, ax=ax,label='All grains')  
quat.Quat.plotIPF(multorilist, [1,0,0], 'hexagonal', c='0.4', marker='o', s=markersize, ax=ax,label='Multi slip')
quat.Quat.plotIPF(prorilist, [1,0,0], 'hexagonal', c='g', marker='s', s=markersize, ax=ax,label='Prismatic only')
quat.Quat.plotIPF(pyorilist, [1,0,0], 'hexagonal', c='r', marker='^', s=markersize, ax=ax,label='Pyramidal only')
quat.Quat.plotIPF(baorilist, [1,0,0], 'hexagonal', c='b', marker='h', s=markersize, ax=ax,label='Basal only')
plt.legend()
f.tight_layout()
#f.savefig(output + 'B2_zone1_step13.png', dpi=1200)