Changeset 3718
- Timestamp:
- Nov 5, 2018 5:30:35 AM (5 years ago)
- Files:
-
- 41 added
- 30 edited
Legend:
- Unmodified
- Added
- Removed
-
Tutorials/CWInstDemo/FindProfParamCW.htm
r3569 r3718 25 25 <o:Author>Brian Toby</o:Author> 26 26 <o:LastAuthor>Brian Toby</o:LastAuthor> 27 <o:Revision> 4</o:Revision>28 <o:TotalTime>3 08</o:TotalTime>27 <o:Revision>6</o:Revision> 28 <o:TotalTime>372</o:TotalTime> 29 29 <o:Created>2018-04-21T19:50:00Z</o:Created> 30 <o:LastSaved>2018- 08-23T20:44:00Z</o:LastSaved>31 <o:Pages>1 3</o:Pages>32 <o:Words> 2546</o:Words>33 <o:Characters>1 4518</o:Characters>30 <o:LastSaved>2018-11-04T06:56:00Z</o:LastSaved> 31 <o:Pages>15</o:Pages> 32 <o:Words>3160</o:Words> 33 <o:Characters>18013</o:Characters> 34 34 <o:Company>Argonne National Laboratory</o:Company> 35 <o:Lines>1 20</o:Lines>36 <o:Paragraphs> 34</o:Paragraphs>37 <o:CharactersWithSpaces> 17030</o:CharactersWithSpaces>35 <o:Lines>150</o:Lines> 36 <o:Paragraphs>42</o:Paragraphs> 37 <o:CharactersWithSpaces>21131</o:CharactersWithSpaces> 38 38 <o:Version>14.0</o:Version> 39 39 </o:DocumentProperties> … … 382 382 mso-generic-font-family:auto; 383 383 mso-font-pitch:variable; 384 mso-font-signature: -536870145 1791491579 18 0 1312310;}384 mso-font-signature:1 134676480 16 0 131072 0;} 385 385 @font-face 386 386 {font-family:"\FF2D\FF33 \30B4\30B7\30C3\30AF"; … … 520 520 {mso-style-noshow:yes; 521 521 mso-style-priority:99; 522 mso-margin-top-alt:auto; 523 margin-right:0in; 524 mso-margin-bottom-alt:auto; 525 margin-left:0in; 526 mso-pagination:widow-orphan; 522 527 font-size:10.0pt; 523 528 font-family:Times; … … 717 722 </style> 718 723 <![endif]--><!--[if gte mso 9]><xml> 719 <o:shapedefaults v:ext="edit" spidmax="10 74"/>724 <o:shapedefaults v:ext="edit" spidmax="1026"/> 720 725 </xml><![endif]--><!--[if gte mso 9]><xml> 721 726 <o:shapelayout v:ext="edit"> … … 733 738 <p class=MsoNormal><o:p> </o:p></p> 734 739 735 <p class=MsoNormal>A video version of this tutorial is available at <span 736 style='mso-field-code:" HYPERLINK \0022https\:\/\/anl\.box\.com\/v\/FindProfParamCW\0022 \\t \0022_blank\0022 "'><span 737 class=MsoHyperlink><b>https://anl.box.com/v/FindProfParamCW</b></span></span></p> 738 739 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 740 741 <p class=MsoNormal><span style='font-size:14.0pt'>The goal of this exercise is to 742 determine approximate instrument profile parameters by a quick peak fit to a 740 <p class=MsoNormal>A video version of the first two parts of this tutorial is 741 available at <a href="https://anl.box.com/v/FindProfParamCW"><b>https://anl.box.com/v/FindProfParamCW</b></a></p> 742 743 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 744 745 <p class=MsoNormal><span style='font-size:14.0pt'>The goal of this exercise is 746 to determine approximate instrument profile parameters by a quick peak fit to a 743 747 sample with a lab diffractometer. Ideally, one should use a material or mixture 744 748 of materials that has peaks over the entire range where you collect data and … … 755 759 style='font-size:14.0pt'>Note</span></b><span style='font-size:14.0pt'> that it 756 760 is highly recommended to collect reference data to a much higher 2theta angle 757 than here. What is done here serves as an example, and is sufficient only if758 on e will never collect and use data above 70 degrees (unlikely!). <o:p></o:p></span></p>761 than here. What is done in this part serves as an example, and is sufficient 762 only if one will never collect and use data above 70 degrees (unlikely!). <o:p></o:p></span></p> 759 763 760 764 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> … … 774 778 mso-layout-grid-align:none;text-autospace:none'><span style='font-size:14.0pt'>Use 775 779 <b style='mso-bidi-font-weight:normal'>Import/Powder Data/from <span 776 class=SpellE>Bru cker</span> RAW file</b> to read </span><span style='font-size:780 class=SpellE>Bruker</span> RAW file</b> to read </span><span style='font-size: 777 781 14.0pt;font-family:Courier'>file </span><span style='font-size:14.0pt; 778 782 font-family:Courier;mso-bidi-font-family:"Menlo Regular";color:black; … … 782 786 directory (download from <a 783 787 href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/CWInstDemo/data/">https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/CWInstDemo/data/</a>). 784 After selecting this file, answer yes to "Is this the file you want?"<o:p></o:p></span></p>788 After selecting this file, answer yes to "Is this the file you want?"<o:p></o:p></span></p> 785 789 786 790 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> … … 792 796 793 797 <p class=MsoNormal><span style='font-size:14.0pt'>The next dialog to appear is 794 titled, "Choose inst. <span class=SpellE>Param</span> file for LaB6_Jan2018.raw795 Scan 1 (or Cancel for default)" since we do not have a set of parameters to 796 read, we must use a default set. <b style='mso-bidi-font-weight:normal'>Press 797 Cancel</b>.<o:p></o:p></span></p>798 titled, "Choose inst. <span class=SpellE>Param</span> file for 799 LaB6_Jan2018.raw Scan 1 (or Cancel for default)" since we do not have a 800 set of parameters to read, we must use a default set. <b style='mso-bidi-font-weight: 801 normal'>Press Cancel</b>.<o:p></o:p></span></p> 798 802 799 803 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> … … 858 862 normal'>Click on the Limits</b> data tree item, and either change the <span 859 863 class=SpellE>Tmin</span> value from 10 to 20 or in the plot <b 860 style='mso-bidi-font-weight:normal'> "drag" the green line to the right</b> to861 approximately 20 degrees. <o:p></o:p></span></p>864 style='mso-bidi-font-weight:normal'>"drag" the green line to the 865 right</b> to approximately 20 degrees. <o:p></o:p></span></p> 862 866 863 867 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; … … 877 881 878 882 <p class=MsoNormal><span style='font-size:14.0pt'>To define peaks, <b 879 style='mso-bidi-font-weight:normal'>click on the "Peak List</b>" data tree item.880 Note that as below, the peak list is initially empty. <o:p></o:p></span></p>883 style='mso-bidi-font-weight:normal'>click on the "Peak List</b>" data 884 tree item. Note that as below, the peak list is initially empty. <o:p></o:p></span></p> 881 885 882 886 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; … … 893 897 style='font-size:14.0pt'>Move</span></b><span style='font-size:14.0pt'> the 894 898 mouse to any of the data points close to the <b style='mso-bidi-font-weight: 895 normal'>top of the first peak and click </b>the left mouse button. A line will 896 be drawn through the peak and the position will be added to the peak table and 897 aline marking this is added to the plot, as below. <o:p></o:p></span></p>899 normal'>top of the first peak and click </b>the left mouse button. A line will be 900 drawn through the peak and the position will be added to the peak table and a 901 line marking this is added to the plot, as below. <o:p></o:p></span></p> 898 902 899 903 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; … … 910 914 style='font-size:14.0pt'>Repeat this for all 8 peaks</span></b><span 911 915 style='font-size:14.0pt'> in the pattern. Note that if a peak is entered in the 912 wrong place <span class=GramE>it can be moved by "dragging" it with the mouse</span>,913 or use a right-click to delete it. Be careful to make sure two peaks are not 914 entered in the same place by accident<span class=GramE>,.</span><o:p></o:p></span></p>916 wrong place <span class=GramE>it can be moved by "dragging" it with 917 the mouse</span>, or use a right-click to delete it. Be careful to make sure 918 two peaks are not entered in the same place by accident.<o:p></o:p></span></p> 915 919 916 920 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; … … 988 992 style='font-size:14.0pt'><span style="mso-spacerun:yes"> </span><o:p></o:p></span></p> 989 993 990 <p class=MsoNormal><span style='font-size:14.0pt'>Select "vary all" and press991 OK. Now all peak positions and areas will be refined.<o:p></o:p></span></p>994 <p class=MsoNormal><span style='font-size:14.0pt'>Select "vary all" 995 and press OK. Now all peak positions and areas will be refined.<o:p></o:p></span></p> 992 996 993 997 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; … … 1005 1009 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1006 1010 1007 <p class=MsoNormal><span style='font-size:14.0pt'>The fit improves 1008 significantly,as below, but further improvements are needed. <o:p></o:p></span></p>1011 <p class=MsoNormal><span style='font-size:14.0pt'>The fit improves significantly, 1012 as below, but further improvements are needed. <o:p></o:p></span></p> 1009 1013 1010 1014 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; … … 1018 1022 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1019 1023 1020 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 7:1021 Refine peak areas, heights and widths (optional)<o:p></o:p></span></h3>1024 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 1025 7: Refine peak areas, heights and widths (optional)<o:p></o:p></span></h3> 1022 1026 1023 1027 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> … … 1128 1132 <p class=MsoNormal><span style='font-size:14.0pt'>Note that these U, V & W 1129 1133 values will be used to set the Gaussian peak widths for those peaks where sigma 1130 is not being refined and likewise X, Y & Z will be used to determine the 1131 Lorentzian widths where gamma is not refined for that peak. If we had any peaks 1132 that were not consistent with the width of the others, we might choose to 1133 continue to refine their sigma & gamma values so they would not affect U, <span 1134 class=GramE>V,É,</span> but here we will refine U, V, W, X & Y against all 1135 peaks. Select the <b style='mso-bidi-font-weight:normal'>Peak List data</b> 1136 tree item and <b style='mso-bidi-font-weight:normal'>remove refinement of sigma 1134 is not being refined and likewise X, Y & Z will be used to determine the Lorentzian 1135 widths where gamma is not refined for that peak. If we had any peaks that were 1136 not consistent with the width of the others, we might choose to continue to 1137 refine their individual profile terms (sigma & gamma), if that is done the 1138 peaks varied individually would be excluded from the U, V, W, X & Y 1139 refinements. Here we will refine U, V, W, X & Y against all peaks. To do 1140 this, select the <b style='mso-bidi-font-weight:normal'>Peak List data</b> tree 1141 item and <b style='mso-bidi-font-weight:normal'>remove refinement of sigma 1137 1142 & gamma</b> for all peaks by double-clicking on the refine column headers 1138 for each and select "N – vary none" so that the table appears as below:<o:p></o:p></span></p> 1143 for each and select "N – vary none" so that the table appears 1144 as below:<o:p></o:p></span></p> 1139 1145 1140 1146 <p class=MsoNormal><span style='font-size:14.0pt'><!--[if gte vml 1]><v:shape … … 1273 1279 <p class=MsoNormal><i style='mso-bidi-font-style:normal'><span 1274 1280 style='font-size:14.0pt'>These profile terms are more than adequate for most 1275 structural fitting problems, but it should be noted that these values are not 1276 sufficientto obtain quantitative measurements of <span class=SpellE>microstrain</span>1281 structural fitting problems, but it should be noted that these values are not sufficient 1282 to obtain quantitative measurements of <span class=SpellE>microstrain</span> 1277 1283 and/or crystallite size. To obtain even better terms the instrument profile is 1278 1284 best determined by a Rietveld fit for (with a standard where <span … … 1305 1311 to the prompt to save the current project (No would not hurt.) As before, <span 1306 1312 class=GramE>use <b style='mso-bidi-font-weight:normal'>Import/Powder Data/from <span 1307 class=SpellE>Bru cker</span> RAW file</b> to read <span style='font-family:Courier'>file</span></span></span><span1313 class=SpellE>Bruker</span> RAW file</b> to read <span style='font-family:Courier'>file</span></span></span><span 1308 1314 style='font-size:14.0pt;font-family:Courier'> </span><span style='font-size: 1309 1315 14.0pt;font-family:Courier;mso-bidi-font-family:"Menlo Regular";color:black; … … 1313 1319 directory (downloaded already from <a 1314 1320 href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/CWInstDemo/data/">https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/CWInstDemo/data/</a>). 1315 After selecting this file, answer yes to "Is this the file you want?"<o:p></o:p></span></p>1321 After selecting this file, answer yes to "Is this the file you want?"<o:p></o:p></span></p> 1316 1322 1317 1323 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> … … 1323 1329 1324 1330 <p class=MsoNormal><span style='font-size:14.0pt'>The next dialog to appear is 1325 titled, "Choose inst. <span class=SpellE>Param</span> file for LaB6_Jan2018.raw1326 Scan 1 unlike before, we do have profile terms to read, from the file created 1327 in Step 10, above. Select the file written before (such as <<span 1328 class=SpellE><i style='mso-bidi-font-style:normal'>MyInstrument</i></span><i1331 titled, "Choose inst. <span class=SpellE>Param</span> file for 1332 LaB6_Jan2018.raw Scan 1 unlike before, we do have profile terms to read, from 1333 the file created in Step 10, above. Select the file written before (such as 1334 <<span class=SpellE><i style='mso-bidi-font-style:normal'>MyInstrument</i></span><i 1329 1335 style='mso-bidi-font-style:normal'>></i>20-<span class=GramE>70deg.instparm 1330 1336 )</span> and <b style='mso-bidi-font-weight:normal'>press OK</b>.<span … … 1510 1516 style='font-size:14.0pt'><o:p></o:p></span></p> 1511 1517 1512 <p class=MsoNormal><span style='font-size:14.0pt'>Nothing needs to be changed of1513 the new parameters added here:<o:p></o:p></span></p>1518 <p class=MsoNormal><span style='font-size:14.0pt'>Nothing needs to be changed 1519 of the new parameters added here:<o:p></o:p></span></p> 1514 1520 1515 1521 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; … … 1674 1680 <p class=MsoNormal><span style='font-size:14.0pt'>Use the <b style='mso-bidi-font-weight: 1675 1681 normal'>Calculate/Refine</b> menu item to start another refinement.<span 1676 style="mso-spacerun:yes"> </span>A warning message that "Steepest1677 Descents dominates "this is because a high degree of parameter correlation1682 style="mso-spacerun:yes"> </span>A warning message that "Steepest 1683 Descents dominates" this is because a high degree of parameter correlation 1678 1684 occurs as the reflection intensities change. The fit does not improve very much 1679 1685 if <b style='mso-bidi-font-weight:normal'>Calculate/Refine</b> is used a second … … 1730 1736 refinement cycle brings the <span class=SpellE>Rw</span> to circa 5.8%. The fit 1731 1737 improves dramatically and <span class=SpellE>microstrain</span> refines to a 1732 non-physical negative value. This is because our somewhat na ve approach above1738 non-physical negative value. This is because our somewhat naive approach above 1733 1739 assumed no sample broadening. <o:p></o:p></span></p> 1734 1740 … … 1738 1744 style='font-size:14.0pt'>This shows why if you are planning to measure <span 1739 1745 class=SpellE>microstrain</span> and/or crystallite size quantitatively, you 1740 should determine instrument profile using a Rietveld fit with a standard with known1741 values for the <span class=SpellE>microstrain</span> and size. Fix the size and 1742 <span class=SpellE>microstrain</span> to the known values when fitting the data 1743 from the standard. <o:p></o:p></span></i></p>1746 should determine instrument profile using a Rietveld fit with a standard with 1747 known values for the <span class=SpellE>microstrain</span> and size. Fix the 1748 size and <span class=SpellE>microstrain</span> to the known values when fitting 1749 the data from the standard. <o:p></o:p></span></i></p> 1744 1750 1745 1751 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step … … 1790 1796 <v:imagedata src="FindProfParamCW_files/image091.png" o:title=""/> 1791 1797 </v:shape><![endif]--><![if !vml]><img border=0 width=533 height=397 1792 src="FindProfParamCW_files/image092.png" v:shapes="Picture_x0020_48"><![endif]></span><span 1798 src="FindProfParamCW_files/image092.png" v:shapes="Picture_x0020_48"><![endif]><o:p></o:p></span></p> 1799 1800 <h1>Part 3: Fitting data over a wider range<o:p></o:p></h1> 1801 1802 <p class=MsoNormal><o:p> </o:p></p> 1803 1804 <p class=MsoNormal><span style='font-size:14.0pt'>In this section we fit a 1805 pattern with a wider range of data. This proceeds quite simply. <o:p></o:p></span></p> 1806 1807 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 1808 3.1: Start new project and read in data<o:p></o:p></span></h3> 1809 1810 <p class=MsoNormal><o:p> </o:p></p> 1811 1812 <p class=MsoNormal><span style='font-size:14.0pt'>Restart GSAS-II or use the 1813 File/New <span class=GramE>Project<span style="mso-spacerun:yes"> 1814 </span>to</span> create a new empty project. <o:p></o:p></span></p> 1815 1816 <p class=MsoNormal style='mso-pagination:none;tab-stops:28.0pt 56.0pt 84.0pt 112.0pt 140.0pt 168.0pt 196.0pt 224.0pt 3.5in 280.0pt 308.0pt 336.0pt; 1817 mso-layout-grid-align:none;text-autospace:none'><span style='font-size:14.0pt'>Use 1818 <b style='mso-bidi-font-weight:normal'>Import/Powder Data/from GSAS powder data 1819 file</b> to read </span><span style='font-size:14.0pt;font-family:Courier'>file 1820 NIST660CBI.gsas</span><span style='font-size:14.0pt;font-family:"Menlo Regular"; 1821 color:black;mso-fareast-language:EN-US'> </span><span style='font-size:14.0pt'>from 1822 the Tutorials <span class=SpellE>CWInstDemo</span>/data directory (download 1823 from <a 1824 href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/CWInstDemo/data/">https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/CWInstDemo/data/</a>). 1825 After selecting this file, answer yes to "Is this the file you want?<span 1826 class=GramE>".</span> As before, use Cancel to select the default <span 1827 class=SpellE>CuKa</span> lab data instrument parameters. <o:p></o:p></span></p> 1828 1829 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 1830 3.2: Change Data Range<o:p></o:p></span></h3> 1831 1832 <p class=MsoNormal><o:p> </o:p></p> 1833 1834 <p class=MsoNormal><span style='font-size:14.0pt'>Since the background is quite 1835 curved at lower angles, reducing the angular range to be fit will require fewer 1836 background terms. Select Limits in the data tree and then set <span 1837 class=SpellE>Tmin</span> to 20.<span style="mso-spacerun:yes"> 1838 </span>Then select Background in the tree and change the number of background 1839 parameters to 8. <o:p></o:p></span></p> 1840 1841 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 1842 3.3: Add Peaks<o:p></o:p></span></h3> 1843 1844 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1845 1846 <p class=MsoNormal><span style='font-size:14.0pt'>Select Peak List from the 1847 tree allows adding peaks. Zooming in on the lower portion of the pattern makes 1848 it easier to see all the peaks. Be sure to click on magnification icon a second 1849 time to turn off zoom mode. Then click on a point for each of the 20 peaks in 1850 the pattern. You should have a plot that looks like this: <o:p></o:p></span></p> 1851 1852 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; 1853 mso-no-proof:yes'><!--[if gte vml 1]><v:shape id="_x0000_i1075" type="#_x0000_t75" 1854 style='width:375pt;height:330pt;visibility:visible;mso-wrap-style:square'> 1855 <v:imagedata src="FindProfParamCW_files/image093.png" o:title=""/> 1856 </v:shape><![endif]--><![if !vml]><img border=0 width=377 height=332 1857 src="FindProfParamCW_files/image094.png" v:shapes="_x0000_i1075"><![endif]></span><span 1858 style='font-size:14.0pt'><o:p></o:p></span></p> 1859 1860 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 1861 3.3: Fit Peaks<o:p></o:p></span></h3> 1862 1863 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1864 1865 <p class=MsoNormal><span style='font-size:14.0pt'>As the peaks are initially 1866 added to the table, the refine flag is set for the peak intensity. Fit the 1867 individual peak intensities by using Peak Fitting/<span class=SpellE>Peakfit</span> 1868 menu command. Provide a file name and the <span class=SpellE>Rwp</span> drops 1869 to ~40% with a reasonable background level, but not good peak shapes. <o:p></o:p></span></p> 1870 1871 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1872 1873 <p class=MsoNormal><span style='font-size:14.0pt'>Select the Instrument 1874 Parameter tree item, select the refinement flag for U, V, W, X, Y and SH/L, 1875 then select the Peak List tree entry again and use the Peak Fitting/<span 1876 class=SpellE>Peakfit</span> menu command. The fit improves a bit, but since the 1877 peak positions have not been refined, a good fit is not to be expected. Click 1878 on the refine column heading to the immediate right of the peak positions, then 1879 select ÒY – vary allÓ and again use the Peak Fitting/<span class=SpellE>Peakfit</span> 1880 menu command. The fit improves significantly with an <span class=SpellE>Rwp</span> 1881 value of ~8%. The fit shows very small difference plot:<o:p></o:p></span></p> 1882 1883 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1884 1885 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; 1886 mso-no-proof:yes'><!--[if gte vml 1]><v:shape id="_x0000_i1074" type="#_x0000_t75" 1887 style='width:348pt;height:247pt;visibility:visible;mso-wrap-style:square'> 1888 <v:imagedata src="FindProfParamCW_files/image095.png" o:title=""/> 1889 </v:shape><![endif]--><![if !vml]><img border=0 width=350 height=249 1890 src="FindProfParamCW_files/image096.png" v:shapes="_x0000_i1074"><![endif]></span><span 1891 style='font-size:14.0pt'><span style="mso-spacerun:yes"> </span><o:p></o:p></span></p> 1892 1893 <p class=MsoNormal><span style='font-size:14.0pt'>However, pressing the w key 1894 in the plot window changes the plot appearance and it can be seen that the fit 1895 at low angle is significantly worse than at high angle. This suggests that log 1896 angle asymmetry may not be well fit. <o:p></o:p></span></p> 1897 1898 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; 1899 mso-no-proof:yes'><!--[if gte vml 1]><v:shape id="_x0000_i1073" type="#_x0000_t75" 1900 style='width:338pt;height:232pt;visibility:visible;mso-wrap-style:square'> 1901 <v:imagedata src="FindProfParamCW_files/image097.png" o:title=""/> 1902 </v:shape><![endif]--><![if !vml]><img border=0 width=340 height=234 1903 src="FindProfParamCW_files/image098.png" v:shapes="_x0000_i1073"><![endif]></span><span 1904 style='font-size:14.0pt'><o:p></o:p></span></p> 1905 1906 <p class=MsoNormal><span style='font-size:14.0pt'>To improve this, we can force 1907 the SH/L value to be larger. Instead of the local minimum that was found with 1908 SH/L of ~0.003, we can select the Instrument Parameter tree item and set SH/L 1909 to a larger value, say 0.03. Then select the Peak List tree entry again and use 1910 the Peak Fitting/<span class=SpellE>Peakfit</span> menu command. <span 1911 class=GramE>The<span style="mso-spacerun:yes"> </span><span class=SpellE>Rwp</span></span> 1912 improves somewhat, to circa 5%, but now the difference plot is greatly improved 1913 at low angle with SH/L increasing to 0.063:<o:p></o:p></span></p> 1914 1915 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; 1916 mso-no-proof:yes'><!--[if gte vml 1]><v:shape id="Picture_x0020_50" o:spid="_x0000_i1072" 1917 type="#_x0000_t75" style='width:344pt;height:241pt;visibility:visible; 1918 mso-wrap-style:square'> 1919 <v:imagedata src="FindProfParamCW_files/image099.png" o:title=""/> 1920 </v:shape><![endif]--><![if !vml]><img border=0 width=346 height=243 1921 src="FindProfParamCW_files/image100.png" v:shapes="Picture_x0020_50"><![endif]></span><span 1922 style='font-size:14.0pt'><o:p></o:p></span></p> 1923 1924 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 1925 3.4: Save Instrument Parameters<o:p></o:p></span></h3> 1926 1927 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1928 1929 <p class=MsoNormal><span style='font-size:14.0pt'>Select the Instrument 1930 Parameter tree entry and then use the Operations/Save ProfileÉ menu command to 1931 write a file. This can be used as the starting point for future refinements. 1932 The plot shown at this point compares the mostly meaningless default starting 1933 profile terms (as solid lines) to the fit here (as dashed lines). Using Operations/Load 1934 ProfileÉ resets the default values to the current refinement parameters, which 1935 makes for a simpler plot. <o:p></o:p></span></p> 1936 1937 <h3><span style='font-size:14.0pt;mso-bidi-font-family:"Times New Roman"'>Step 1938 3.5: Compare individual peaks widths to fit profile<o:p></o:p></span></h3> 1939 1940 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1941 1942 <p class=MsoNormal><span style='font-size:14.0pt'>While the Instrument 1943 Parameter tree entry is selected, turn off refinement of the U, V, W X, and Y refinement 1944 flags (SH/L can be left to continue to refine) and then select the Peak List 1945 tree entry. <span style="mso-spacerun:yes"> </span>Double-click on the 1946 refine headings for the sigma and gamma columns to refine the individual peak widths 1947 and again and use the Peak Fitting/<span class=SpellE>Peakfit</span> menu 1948 command. The fit quality does not change very significantly, and selecting the Instrument 1949 Parameter tree entry generates a plot that shows that the individual peak fits 1950 are quite close to the fitted curves. <o:p></o:p></span></p> 1951 1952 <p class=MsoNormal><span style='font-size:14.0pt'><o:p> </o:p></span></p> 1953 1954 <p class=MsoNormal><span style='font-size:14.0pt;mso-fareast-language:EN-US; 1955 mso-no-proof:yes'><!--[if gte vml 1]><v:shape id="Picture_x0020_51" o:spid="_x0000_i1071" 1956 type="#_x0000_t75" style='width:330pt;height:236pt;visibility:visible; 1957 mso-wrap-style:square'> 1958 <v:imagedata src="FindProfParamCW_files/image101.png" o:title=""/> 1959 </v:shape><![endif]--><![if !vml]><img border=0 width=332 height=238 1960 src="FindProfParamCW_files/image102.png" v:shapes="Picture_x0020_51"><![endif]></span><span 1793 1961 style='font-size:14.0pt'><o:p></o:p></span></p> 1794 1962 -
Tutorials/CWInstDemo/FindProfParamCW_files/filelist.xml
r3354 r3718 94 94 <o:File HRef="image091.png"/> 95 95 <o:File HRef="image092.png"/> 96 <o:File HRef="image093.png"/> 97 <o:File HRef="image094.png"/> 98 <o:File HRef="image095.png"/> 99 <o:File HRef="image096.png"/> 100 <o:File HRef="image097.png"/> 101 <o:File HRef="image098.png"/> 102 <o:File HRef="image099.png"/> 103 <o:File HRef="image100.png"/> 104 <o:File HRef="image101.png"/> 105 <o:File HRef="image102.png"/> 96 106 <o:File HRef="filelist.xml"/> 97 107 </xml> -
Tutorials/FitPeaks/Fit Peaks.htm
r3569 r3718 1722 1722 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1723 1723 mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Calculate/Refine</span></b> 1724 menu item; that it to beused for Rietveld/Pawley refinements. You will be1724 menu item; that should used for Rietveld/Pawley refinements. You will be 1725 1725 asked to save the current project before the fit will proceed; change directory 1726 1726 if desired and provide a suitable name (we assume <b style='mso-bidi-font-weight: -
Tutorials/Magnetic-I/Magnetic Structures-I_files/image001.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image011.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image028.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image032.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image042.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image046.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image048.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image052.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image054.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image056.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image058.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image060.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image064.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image066.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image068.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image070.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image072.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image074.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image076.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image078.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image080.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image082.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image084.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
Tutorials/Magnetic-I/Magnetic Structures-I_files/image086.png
- Property svn:mime-type changed from text/html to application/octet-stream
-
trunk/GSASIIctrlGUI.py
r3714 r3718 5149 5149 '''This shows how to determine profile parameters by fitting individual peaks 5150 5150 with data collected on a standard using a lab diffractometer.'''], 5151 5151 ['FPAfit', 'FPAfit.htm', 'Determining Profile Parameters with Fundamental Parameters', 5152 '''This shows how to determine profile parameters by fitting 5153 peaks that are computed using the NIST Fundamental Parameters Python 5154 code. 5155 Input is formulated to use FPA values similar to those in Topas.'''], 5152 5156 ['TOF Calibration', 'Calibration of a TOF powder diffractometer.htm', 'Calibration of a Neutron TOF diffractometer', 5153 5157 '''This uses the fitted positions of all visible peaks in a pattern of NIST SRM 660b La11B6 -
trunk/GSASIImapvars.py
r3713 r3718 31 31 references 32 32 to numbered objects are resolved using the appropriate random ids and the 33 variable object is expressed asa string of form ``ph:hst:VARNAM:at``.33 parameter object is converted to a string of form ``ph:hst:VARNAM:at``. 34 34 35 35 Alternate parameters (New Var) … … 55 55 degree of freedom where parameter Px and Py, etc. are varied to keep 56 56 their ratio 57 fixed according the expression. A new variable is assigned to each degree of57 fixed according the expression. A new variable parameter is assigned to each degree of 58 58 freedom when refined. An example where this can be valuable is when 59 59 two parameters, P1 and P2, have similar values and are highly correlated. It is often better to create a new variable, Ps = P1 + P2, and refine Ps. … … 81 81 a formula, but this can be expanded if needed. 82 82 83 These equations set an interdependence between variables.83 These equations set an interdependence between parameters. 84 84 Common uses of parameter constraints are to set rules that decrease the number of parameters, 85 85 such as restricting the sum of fractional occupancies for atoms that share … … 154 154 ^^^^^^^^^^^^^^^^^^^^^^^^ 155 155 156 When variables are refined where a single refinement flag determines that several variables156 When parameters are refined where a single refinement flag determines that several variables 157 157 are refined at the same time (examples are: cell parameters, atom positions, anisotropic 158 158 displacement parameters, magnetic moments,...) it can be useful to specify that a 159 159 specific parameter should not be varied. These will most commonly be generated due to symmetry, 160 but under specific conditions, there may be other good reasons to constrain a variable.160 but under specific conditions, there may be other good reasons to constrain a parameter. 161 161 162 162 A "Hold" constraint is stored as described for type "h" in the … … 361 361 362 362 Rows in the matrix corresponding to "New Var" constraints and those that 363 were generated by the Gram-Schmidt process are provided with variablenames363 were generated by the Gram-Schmidt process are provided with parameter names 364 364 (this can be specified if a "New Var" entry by using a ``"_name"`` element 365 365 in the constraint dict, but at present this is not implemented.) Names are … … 488 488 489 489 fixedVarList: 490 a list of variables that have been 'fixed'490 a list of parameters that have been 'fixed' 491 491 by defining them as equal to a constant (::var: = 0). Note that 492 the constant value is ignored at present. These variables are492 the constant value is ignored at present. These parameters are 493 493 later removed from varyList which prevents them from being refined. 494 494 Unlikely to be used externally. … … 504 504 505 505 problemVars: 506 a list containing variables that show up in constraints producing errors506 a list containing parameters that show up in constraints producing errors 507 507 508 508 … … 548 548 ''' 549 549 fixedVarList = [] 550 '''List of variables that should not be refined.550 '''List of parameters that should not be refined. 551 551 ''' 552 552 symGenList = [] … … 554 554 ''' 555 555 problemVars = [] 556 '''a list of variables causing errors556 '''a list of parameters causing errors 557 557 ''' 558 558 dependentVars = [] … … 561 561 'A list of dependent variables, taken from (:data:`indParmList`).' 562 562 genVarLookup = {} 563 'provides a list of variables that are related to each generated variable'563 'provides a list of parameters that are related to each generated parameter' 564 564 paramPrefix = "::constr" 565 565 'A prefix for generated parameter names' … … 586 586 587 587 def VarKeys(constr): 588 """Finds the keys in a constraint that represent variables588 """Finds the keys in a constraint that represent parameters 589 589 e.g. eliminates any that start with '_' 590 590 … … 687 687 # varied, create a warning message. 688 688 for group,varlist in zip(groups,parmlist): 689 if len(varlist) == 1: continue 689 if len(varlist) == 1: # process fixed (held) variables 690 var = varlist[0] 691 if var not in fixedVarList: 692 fixedVarList.append(var) 693 continue 690 694 for rel in group: 691 695 varied = 0 … … 694 698 if var.startswith('_'): continue 695 699 if not re.match('[0-9]*:[0-9\*]*:',var): 696 warnmsg += "\n Variable"+str(var)+" does not begin with a ':'"700 warnmsg += "\nParameter "+str(var)+" does not begin with a ':'" 697 701 if var in varyList: 698 702 varied += 1 … … 705 709 errmsg += _FormatConstraint(constrDict[rel],fixedList[rel])+"\n" 706 710 if varied > 0 and varied != len(VarKeys(constrDict[rel])): 707 warnmsg += "\nNot all variables refined in constraint:\n\t"711 warnmsg += "\nNot all parameters refined in constraint:\n\t" 708 712 warnmsg += _FormatConstraint(constrDict[rel],fixedList[rel]) 709 713 warnmsg += '\nNot refined: ' + notvaried + '\n' … … 714 718 # non-singular square matrix 715 719 for group,varlist in zip(groups,parmlist): 716 if len(varlist) == 1: continue # a constraint group with a single variablecan be ignored720 if len(varlist) == 1: continue # a constraint group with a single parameter can be ignored 717 721 if len(varlist) < len(group): # too many relationships -- no can do 718 722 errmsg += "\nOver-constrained input. " 719 723 errmsg += "There are more constraints " + str(len(group)) 720 errmsg += "\n\tthan variables " + str(len(varlist)) + "\n"724 errmsg += "\n\tthan parameters " + str(len(varlist)) + "\n" 721 725 for rel in group: 722 726 errmsg += _FormatConstraint(constrDict[rel],fixedList[rel]) … … 744 748 mapvar = [] 745 749 group = group[:] 746 # scan through all generated and input variables747 # Check again for inconsistent variableuse748 # for new variables -- where varied and unvaried parameters get grouped750 # scan through all generated and input parameters 751 # Check again for inconsistent parameter use 752 # for new parameters -- where varied and unvaried parameters get grouped 749 753 # together. I don't think this can happen when not flagged before, but 750 754 # it does not hurt to check again. … … 764 768 fixedval = None 765 769 if fixedval is None: 766 varname = paramPrefix + str(consNum) # assign a name to a variable770 varname = paramPrefix + str(consNum) # assign a name to a parameter 767 771 mapvar.append(varname) 768 772 consNum += 1 … … 770 774 mapvar.append(fixedval) 771 775 if varied > 0 and notvaried != '': 772 warnmsg += "\nNot all variables refined in generated constraint"776 warnmsg += "\nNot all parameters refined in generated constraint" 773 777 warnmsg += '\nPlease report this unexpected error\n' 774 778 for rel in group: … … 793 797 '''Takes a list of relationship entries comprising a group of 794 798 constraints and builds the relationship lists and their inverse 795 and stores them in global variables Also checks for internal799 and stores them in global parameters Also checks for internal 796 800 conflicts or inconsistencies in parameter/variable definitions. 797 801 … … 810 814 811 815 :param int SeqHist: number of current histogram, when used in a sequential 812 refinement. None (default) otherwise. Wildcard variablenames are816 refinement. None (default) otherwise. Wildcard parameter names are 813 817 set to the current histogram, when found if not None. 814 818 ''' … … 818 822 shortmsg = '' 819 823 changed = False 824 820 825 # Process the equivalences 821 826 # If there are conflicting parameters, move them into constraints. This … … 840 845 groups,parmlist = GroupConstraints(constrDict) 841 846 for group,varlist in zip(groups,parmlist): 842 if len(varlist) == 1: continue 847 if len(varlist) == 1: # process fixed (held) variables 848 var = varlist[0] 849 if var not in fixedVarList: 850 fixedVarList.append(var) 851 continue 843 852 for rel in group: 844 853 varied = 0 … … 868 877 #if unused > 0:# and unused != len(VarKeys(constrDict[rel])): 869 878 if unused > 0 and unused != len(VarKeys(constrDict[rel])): 870 #msg += "\nSome (but not all) variables in constraint are not defined:\n\t"879 #msg += "\nSome (but not all) parameters in constraint are not defined:\n\t" 871 880 #msg += _FormatConstraint(constrDict[rel],fixedList[rel]) 872 881 #msg += '\nNot used: ' + notused + '\n' 873 882 shortmsg += notused+" not used in constraint "+_FormatConstraint(constrDict[rel],fixedList[rel]) 874 883 elif varied > 0 and varied != len(VarKeys(constrDict[rel])): 875 #msg += "\nNot all variables refined in constraint:\n\t"884 #msg += "\nNot all parameters refined in constraint:\n\t" 876 885 #msg += _FormatConstraint(constrDict[rel],fixedList[rel]) 877 886 #msg += '\nNot refined: ' + notvaried + '\n' … … 899 908 for group,varlist in zip(groups,parmlist): 900 909 if len(varlist) == 1: continue 901 # for constraints, if all included variables are refined,910 # for constraints, if all included parameters are refined, 902 911 # set the VaryFree flag, and remaining degrees of freedom will be 903 # varied (since consistency was checked, if any one variableis912 # varied (since consistency was checked, if any one parameter is 904 913 # refined, then assume that all are) 905 varsList = [] # make a list of all the referenced variables as well914 varsList = [] # make a list of all the referenced parameters as well 906 915 VaryFree = False 907 916 for rel in group: … … 919 928 msg = 'too many relationships' 920 929 break 921 # Since we checked before, if any variables are unused, then all must be.930 # Since we checked before, if any parameters are unused, then all must be. 922 931 # If so, this set of relationships can be ignored 923 932 if unused: 924 if debug: print('Constraint ignored (all variables undefined)')933 if debug: print('Constraint ignored (all parameters undefined)') 925 934 if debug: print (' '+_FormatConstraint(constrDict[rel],fixedList[rel])) 926 935 continue … … 953 962 varname = '' 954 963 fixedval = None 955 if fixedval is None: # this is a new variable, not a constraint964 if fixedval is None: # this is a new parameter, not a constraint 956 965 if not varname: 957 966 varname = paramPrefix + str(consNum) # no assigned name, create one 958 967 consNum += 1 959 968 mapvar.append(varname) 960 genVarLookup[varname] = varlist # save list of variables related to this new var969 genVarLookup[varname] = varlist # save list of parameters related to this new var 961 970 # vary the new relationship if it is a degree of freedom in 962 971 # a set of contraint equations or if a New Var is flagged to be varied. … … 964 973 unused = False 965 974 varyList.append(varname) 966 # fix (prevent varying) of all the variables inside the constraint group975 # fix (prevent varying) of all the parameters inside the constraint group 967 976 # (dependent vars) 968 977 for var in varsList: … … 971 980 unused = False 972 981 mapvar.append(fixedval) 973 if unused: # skip over constraints that don't matter (w/o fixed value or any refined variables)974 if debug: print('Constraint ignored (all variables unrefined)')982 if unused: # skip over constraints that don't matter (w/o fixed value or any refined parameters) 983 if debug: print('Constraint ignored (all parameters unrefined)') 975 984 if debug: print (' '+_FormatConstraint(constrDict[rel],fixedList[rel])) 976 985 continue … … 1132 1141 errmsg = '' 1133 1142 problemVars = [] 1134 # process fixed variables (holds)1143 # process fixed parameters (holds) 1135 1144 fixVlist = [] # list of Fixed vars 1136 1145 constrVars = [] # list of vars in constraint expressions … … 1148 1157 multdepVarList = [] 1149 1158 dropVarList = [] 1150 translateTable = {} # lookup table for wildcard referenced variables1159 translateTable = {} # lookup table for wildcard referenced parameters 1151 1160 for varlist,mapvars,multarr,invmultarr in zip( 1152 1161 dependentParmList,indParmList,arrayList,invarrayList): … … 1168 1177 notvaried += mv 1169 1178 if parmDict is not None and mv not in parmDict: 1170 print ("Dropping equivalence for variable"+str(mv)+". Not defined in this refinement")1179 print ("Dropping equivalence for parameter "+str(mv)+". Not defined in this refinement") 1171 1180 if mv not in dropVarList: dropVarList.append(mv) 1172 1181 if mv not in indepVarList: indepVarList.append(mv) … … 1188 1197 notvaried += v 1189 1198 if v in indepVarList: 1190 errmsg += '\n Variable'+v+' is used to set values in a constraint before its value is set in another constraint\n'1199 errmsg += '\nParameter '+v+' is used to set values in a constraint before its value is set in another constraint\n' 1191 1200 if v not in problemVars: problemVars.append(v) 1192 1201 if v in depVarList: … … 1195 1204 depVarList.append(v) 1196 1205 if varied > 0 and varied != len(varlist)+1: 1197 warnmsg += "\nNot all variables refined in equivalence:\n\t"1206 warnmsg += "\nNot all parameters refined in equivalence:\n\t" 1198 1207 s = "" 1199 1208 for v in varlist: … … 1380 1389 1381 1390 def GetDependentVars(): 1382 '''Return a list of dependent variables: e.g. variables that are1383 constrained in terms of other variables1384 1385 :returns: a list of variablenames1391 '''Return a list of dependent variables: e.g. parameters that are 1392 constrained in terms of other parameters 1393 1394 :returns: a list of parameter names 1386 1395 1387 1396 ''' … … 1390 1399 1391 1400 def GetIndependentVars(): 1392 '''Return a list of independent variables: e.g. variables that are1393 created by constraints of other variables1394 1395 :returns: a list of variablenames1401 '''Return a list of independent variables: e.g. parameters that are 1402 slaved to other parameters by constraints 1403 1404 :returns: a list of parameter names 1396 1405 1397 1406 ''' … … 1400 1409 1401 1410 def PrintIndependentVars(parmDict,varyList,sigDict,PrintAll=False,pFile=None): 1402 '''Print the values and uncertainties on the independent variables'''1411 '''Print the values and uncertainties on the independent parameters''' 1403 1412 global dependentParmList,arrayList,invarrayList,indParmList,fixedDict 1404 1413 printlist = [] … … 1414 1423 s3 = '' 1415 1424 pFile.write(130*'-'+'\n') 1416 pFile.write(" Variables generated by constraints\n")1425 pFile.write("Parameters generated by constraints\n") 1417 1426 printlist.append(3*[None]) 1418 1427 for name,val,esd in printlist: … … 1493 1502 s = '' 1494 1503 if len(fixedVarList) > 0: 1495 s += 'Fixed Variables:\n'1504 s += 'Fixed Parameters:\n' 1496 1505 for v in fixedVarList: 1497 1506 s += ' ' + v + '\n' 1498 1507 if not inputOnly: 1499 s += 'User-supplied variablemapping relations:\n'1508 s += 'User-supplied parameter mapping relations:\n' 1500 1509 symout = '' 1501 1510 global dependentParmList,arrayList,invarrayList,indParmList,fixedDict,symGenList … … 1536 1545 s += 'Symmetry-generated relations:\n' + symout 1537 1546 if inputOnly: return s 1538 s += 'Inverse variablemapping relations:\n'1547 s += 'Inverse parameter mapping relations:\n' 1539 1548 for varlist,mapvars,invmultarr in zip(dependentParmList,indParmList,invarrayList): 1540 1549 for i,mv in enumerate(varlist): … … 1603 1612 for varlist,mapvars,multarr,invmultarr in zip(dependentParmList,indParmList,arrayList,invarrayList): 1604 1613 for i,name in enumerate(mapvars): 1605 # grouped variables: need to add in the derv. w/r1614 # grouped parameters: need to add in the derv. w/r 1606 1615 # dependent variables to the independent ones 1607 1616 if name not in varyList: continue # skip if independent var not varied … … 1628 1637 defined using :func:`StoreEquivalence`, 1629 1638 :func:`GroupConstraints` and :func:`GenerateConstraints` and 1630 before any variablerefinement is done. This completes the parameter1639 before any parameter refinement is done. This completes the parameter 1631 1640 dictionary by defining independent parameters and it satisfies the 1632 1641 constraint equations in the initial parameters … … 1657 1666 np.dot(multarr,np.array(valuelist))) 1658 1667 ) 1659 # now remove fixed variables from the varyList1668 # now remove fixed parameters from the varyList 1660 1669 global fixedVarList 1661 1670 for item in fixedVarList: -
trunk/help/Tutorials.html
r3637 r3718 59 59 [link: <A href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/Magnetic-IV/data">Exercise files</A>]. 60 60 <blockquote><I>Analysis of a Type IV antiferromagnet with a lattice centering change using Bilbao k-SUBGROUPSMAG from CW neutron powder data</I></blockquote> 61 <LI><A href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/Magnetic-V/Magnetic-V.htm">Magnetic Structure Analysis-V</A> 62 [link: <A href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/Magnetic-V/data">Exercise files</A>]. 63 <blockquote><I>Analysis of a complex Type IV antiferromagnet with two propagation vectorse using Bilbao k-SUBGROUPSMAG from TOF neutron powder data</I></blockquote> 61 64 </UL><h4>Parametric sequential fitting</H4><UL> 62 65 <LI><A href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/SeqRefine/SequentialTutorial.htm">Sequential refinement of multiple datasets</A> … … 114 117 <blockquote><I>This shows how to determine profile parameters by fitting individual peaks 115 118 with data collected on a standard using a lab diffractometer.</I></blockquote> 119 <LI><A href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/FPAfit/FPAfit.htm">Determining Profile Parameters with Fundamental Parameters</A> 120 [no example data or video]. 121 <blockquote><I>This shows how to determine profile parameters by fitting 122 peaks that are computed using the NIST Fundamental Parameters Python 123 code. 124 Input is formulated to use FPA values similar to those in Topas.</I></blockquote> 116 125 <LI><A href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/TOF Calibration/Calibration of a TOF powder diffractometer.htm">Calibration of a Neutron TOF diffractometer</A> 117 126 [links: <A href="https://anl.box.com/v/CalibrationofaTOFpowderdiffrac">video</A>, <A href="https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/TOF Calibration/data">Exercise files</A>].
Note: See TracChangeset
for help on using the changeset viewer.