https://wiki.uni-konstanz.de/xds/index.php?title=Thaumatin_ACA2014&feed=atom&action=history Thaumatin ACA2014 - Revision history 2024-03-29T01:11:03Z Revision history for this page on the wiki MediaWiki 1.39.6 https://wiki.uni-konstanz.de/xds/index.php?title=Thaumatin_ACA2014&diff=4099&oldid=prev Kay: ftp -> https 2020-03-24T16:36:55Z <p>ftp -&gt; https</p> <table style="background-color: #fff; color: #202122;" data-mw="interface"> <col class="diff-marker" /> <col class="diff-content" /> <col class="diff-marker" /> <col class="diff-content" /> <tr class="diff-title" lang="en"> <td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">← Older revision</td> <td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">Revision as of 17:36, 24 March 2020</td> </tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l1">Line 1:</td> <td colspan="2" class="diff-lineno">Line 1:</td></tr> <tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. Raw data are at <del style="font-weight: bold; text-decoration: none;">ftp</del>://{{SERVERNAME}}/pub/datasets/ACA2014_thaumatin/ . What is special about these data:</div></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. Raw data are at <ins style="font-weight: bold; text-decoration: none;">https</ins>://{{SERVERNAME}}/pub/datasets/ACA2014_thaumatin/ . What is special about these data:</div></td></tr> <tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr> <tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>*    tiny spots (some are only one pixel)</div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>*    tiny spots (some are only one pixel)</div></td></tr> </table> Kay https://wiki.uni-konstanz.de/xds/index.php?title=Thaumatin_ACA2014&diff=4032&oldid=prev Karsten at 16:24, 19 December 2019 2019-12-19T16:24:31Z <p></p> <table style="background-color: #fff; color: #202122;" data-mw="interface"> <col class="diff-marker" /> <col class="diff-content" /> <col class="diff-marker" /> <col class="diff-content" /> <tr class="diff-title" lang="en"> <td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">← Older revision</td> <td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">Revision as of 17:24, 19 December 2019</td> </tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l1">Line 1:</td> <td colspan="2" class="diff-lineno">Line 1:</td></tr> <tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. Raw data are at ftp://<del style="font-weight: bold; text-decoration: none;">turn5.biologie.uni-konstanz.de</del>/pub/datasets/ACA2014_thaumatin/ . What is special about these data:</div></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. Raw data are at ftp://<ins style="font-weight: bold; text-decoration: none;">{{SERVERNAME}}</ins>/pub/datasets/ACA2014_thaumatin/ . What is special about these data:</div></td></tr> <tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr> <tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>*    tiny spots (some are only one pixel)</div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>*    tiny spots (some are only one pixel)</div></td></tr> </table> Karsten https://wiki.uni-konstanz.de/xds/index.php?title=Thaumatin_ACA2014&diff=3371&oldid=prev Kay: add link to raw data 2016-07-28T11:31:03Z <p>add link to raw data</p> <table style="background-color: #fff; color: #202122;" data-mw="interface"> <col class="diff-marker" /> <col class="diff-content" /> <col class="diff-marker" /> <col class="diff-content" /> <tr class="diff-title" lang="en"> <td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">← Older revision</td> <td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">Revision as of 12:31, 28 July 2016</td> </tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l1">Line 1:</td> <td colspan="2" class="diff-lineno">Line 1:</td></tr> <tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. What is special about these data:</div></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source<ins style="font-weight: bold; text-decoration: none;">. Raw data are at ftp://turn5.biologie.uni-konstanz.de/pub/datasets/ACA2014_thaumatin/ </ins>. What is special about these data:</div></td></tr> <tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr> <tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>*    tiny spots (some are only one pixel)</div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>*    tiny spots (some are only one pixel)</div></td></tr> </table> Kay https://wiki.uni-konstanz.de/xds/index.php?title=Thaumatin_ACA2014&diff=2962&oldid=prev Karsten: Created page with "1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. What is special about these data: * tiny spots (some are only one pixel) * ..." 2014-05-26T11:13:25Z <p>Created page with &quot;1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. What is special about these data: * tiny spots (some are only one pixel) * ...&quot;</p> <p><b>New page</b></p><div>1800 frames of 0.2° were collected on a Pilatus 6M detector at 24-IDC, Advanced Photon Source. What is special about these data:<br /> <br /> * tiny spots (some are only one pixel)<br /> * fast: 0.2s total time for each frame. Since the readout time is only 3ms, a shutter is not used (&quot;shutterless data collection&quot;. <br /> * the header is incorrect: the direct beam position is at 1288.9 1262.6 rather than at 1265.47 1273.87 (as the header claims). This may lead to mis-indexing. Below it is demonstrated how to identify this problem.<br /> <br /> We start processing by creating an empty directory, and changing to this directory. This can be conveniently done using [[XDSGUI]] / Projects / navigate to the directory with the Thaumatin frames / create empty folder.<br /> <br /> <br /> == Running [[XDSGUI]] ==<br /> <br /> We click on the &quot;Frame tab&quot; and load one of the frames. Clicking &quot;Generate XDS.INP&quot; creates XDS.INP. This un-greys the XDS.INP tab, and the blind areas of the detector are shown to be masked (red rectangles). Furthermore, a green &quot;x&quot; appears where the header says that the beam position should be.<br /> <br /> Next, we click the &quot;XDS.INP&quot; tab. We may manually change parameters, like the resolution cutoff used for processing (INCLUDE_RESOLUTION_RANGE= ). In this case we leave it at its default (50), and don&#039;t specify a high-resolution cutoff (0).<br /> <br /> After reviewing and possibly changing parameters, we then &quot;save&quot;, &quot;run XDS&quot;. After a few seconds, the XYCORR, INIT, COLSPOT, IDXREF tabs become black, indicating that XYCORR.LP, INIT.LP, COLSPOT.LP and [[IDXREF.LP]] have been created.<br /> <br /> We may take a look at the XYCORR, INIT, and COLSPOT tabs to see what XDS has to say about these steps.<br /> <br /> == [[IDXREF.LP|IDXREF]] (indexing) ==<br /> <br /> Next, we click the IDXREF tab and find, after scrolling down:<br /> <br /> &lt;pre&gt;<br /> # COORDINATES OF REC. BASIS VECTOR LENGTH 1/LENGTH<br /> <br /> 1 0.0048784-0.0028812 0.0034946 0.0066567 150.22<br /> 2 -0.0111804-0.0117066 0.0059615 0.0172507 57.97<br /> 3 0.0035904-0.0102197-0.0133848 0.0172188 58.08<br /> <br /> CLUSTER COORDINATES AND INDICES WITH RESPECT TO REC. LATTICE BASIS VECTORS <br /> <br /> # COORDINATES OF VECTOR CLUSTER FREQUENCY CLUSTER INDICES <br /> 1 0.0048672-0.0029036 0.0034869 2326. 1.00 0.00 0.00<br /> 2 -0.0209391-0.0058398-0.0010916 2077. -2.01 0.99 -0.00<br /> 3 0.0097412-0.0057891 0.0069462 2065. 2.00 0.00 0.00<br /> 4 0.0160686 0.0087446-0.0024072 2051. 1.01 -1.00 0.00<br /> 5 -0.0258097-0.0029199-0.0045794 1995. -3.01 0.99 -0.00<br /> 6 -0.0111928-0.0116435 0.0059000 1900. -0.01 1.00 0.00<br /> 7 0.0306841 0.0000146 0.0080409 1816. 4.01 -0.99 0.00<br /> 8 -0.0146229 0.0086452-0.0104174 1764. -2.99 0.00 -0.00<br /> 9 -0.0294185 0.0072888 0.0088294 1720. -3.01 0.99 -1.00<br /> 10 0.0245344-0.0044011-0.0123085 1705. 2.01 -0.99 1.00<br /> 11 0.0063073 0.0145033-0.0093667 1673. -0.99 -1.00 -0.00<br /> 12 0.0133648-0.0160316-0.0064337 1664. 2.00 0.00 1.00<br /> 13 0.0418771 0.0116728 0.0021928 1661. 4.03 -1.99 0.00<br /> 14 -0.0342963 0.0101537 0.0053161 1652. -4.01 0.99 -1.00<br /> 15 0.0035949-0.0101595-0.0132820 1627. 0.00 -0.00 0.99<br /> 16 0.0355623-0.0028693 0.0115246 1600. 5.01 -0.99 0.01<br /> 17 -0.0182374 0.0188742 0.0029450 1591. -3.00 0.00 -1.00<br /> ...<br /> &lt;/pre&gt;<br /> <br /> This shows the cell parameters as expected, and the difference vectors are indeed close to integers, as they should be. Continuing ...<br /> <br /> &lt;pre&gt;<br /> ***** RESULTS FROM LOCAL INDEXING OF 3000 OBSERVED SPOTS *****<br /> <br /> MAXIMUM MAGNITUDE OF INDEX DIFFERENCES ALLOWED 8<br /> MAXIMUM ALLOWED DEVIATION FROM INTEGERAL INDICES 0.050<br /> MIMINUM QUALITY OF INDICES FOR EACH SPOT IN A SUBTREE 0.80<br /> QUALITY OF INDICES REQUIRED TO INCLUDE SECOND SUBTREE 0.00<br /> NUMBER OF SUBTREES 16<br /> <br /> SUBTREE POPULATION<br /> <br /> 1 3000<br /> &lt;/pre&gt;<br /> <br /> This shows that the 3000 strongest reflections can be indexed with a single lattice - good ! <br /> <br /> Next ...<br /> <br /> &lt;pre&gt;<br /> NUMBER OF ACCEPTED SPOTS FROM LARGEST SUBTREE 2963<br /> <br /> ***** SELECTION OF THE INDEX ORIGIN OF THE REFLECTIONS *****<br /> The origin of the reflection indices determined so far is <br /> 0,0,0 by default which is usually correct. In certain critical<br /> cases it may happen that this automatic choice is wrong which <br /> leads to misindexing of the reflections by a constant offset. <br /> You may replace the default by specifying INDEX_ORIGIN= h k l <br /> in the input file &quot;XDS.INP&quot; and rerun the IDXREF step. <br /> Below you find a list of possible alternatives together with a<br /> measure of their likelihood.<br /> QUALITY small values mean a high likelihood for this offset <br /> DELTA is the angle between given and refined beam direction<br /> XD,YD computed direct beam position (pixels) on detector <br /> given beam position (pixel): 1265.47 1273.87<br /> X,Y,Z computed coordinates of the direct beam wave vector<br /> DH,DK,DL mean absolute difference between observed and <br /> fitted indices<br /> <br /> INDEX_ QUALITY DELTA XD YD X Y Z DH DK DL<br /> ORIGIN<br /> <br /> 0 0 0 14.4 1.0 1289.1 1262.5 0.0131 -0.0063 0.8065 0.02 0.03 0.06<br /> 0 0 1 31.4 0.7 1280.3 1266.3 0.0082 -0.0042 0.8065 0.16 0.05 0.37<br /> 1 0 3 42.8 0.9 1242.6 1267.8 -0.0127 -0.0034 0.8065 0.10 0.22 0.15<br /> 1 0 2 61.3 0.7 1251.3 1264.0 -0.0078 -0.0055 0.8065 0.22 0.23 0.34<br /> 1 0 4 74.1 1.3 1233.8 1271.6 -0.0176 -0.0013 0.8064 0.14 0.22 0.40<br /> 0 0 2 99.9 0.3 1271.5 1270.1 0.0034 -0.0021 0.8066 0.32 0.10 0.73<br /> 1 0 1 131.4 0.6 1260.2 1260.2 -0.0030 -0.0076 0.8065 0.38 0.23 0.70<br /> ...<br /> &lt;/pre&gt;<br /> <br /> In this case, the &quot;QUALITY&quot; of the default indexing (first line, starting with 0 0 0) is surprisingly bad (1.0 would be a good value), and the angular distance between the given and refined beam direction is 1.0°, which is actually a lot. Nevertheless, the &quot;mean absolute difference between observed and fitted indices&quot; DH, DK, DL are the smallest in the list, which gives us some confidence in the indexing. <br /> <br /> Since we did not specify space group and cell parameters, XDS integrates the data in P1 (the default), and we postpone space group determination to the scaling and merging ([[CORRECT]]) step (below).<br /> <br /> === Digression: a computational investigation into the robustness of the indexing ===<br /> <br /> This indexing was based on 206252 spots found (by COLSPOT) in the first half of the DATA_RANGE, i.e. using SPOT_RANGE=1 900. If we had taken e.g. a single frame only (SPOT_RANGE=1 1; we do not have to re-run COLSPOT to try this - just use JOB=IDXREF) then we would have indexed based on 392 spot positions, and the last table would be:<br /> <br /> &lt;pre&gt;<br /> INDEX_ QUALITY DELTA XD YD X Y Z DH DK DL<br /> ORIGIN<br /> <br /> 0 0 0 13.2 0.6 1279.6 1268.0 0.0078 -0.0033 0.8065 0.15 0.07 0.26<br /> 0 0 1 13.6 1.0 1288.4 1262.8 0.0127 -0.0061 0.8065 0.03 0.03 0.07<br /> 0 -1 -1 15.3 1.0 1250.7 1251.8 -0.0082 -0.0122 0.8064 0.04 0.04 0.07<br /> 1 0 -2 20.9 1.0 1255.5 1296.5 -0.0055 0.0126 0.8065 0.16 0.07 0.26<br /> 0 -1 0 21.5 1.1 1259.6 1246.8 -0.0033 -0.0150 0.8064 0.12 0.06 0.21<br /> 1 0 -3 23.1 1.3 1246.6 1301.8 -0.0105 0.0155 0.8064 0.03 0.03 0.07<br /> 0 0 -1 28.3 0.2 1270.7 1273.2 0.0029 -0.0004 0.8066 0.31 0.14 0.53<br /> ...<br /> &lt;/pre&gt;<br /> <br /> Here, the &#039;&#039;second&#039;&#039; line (starting with 0 0 1) has low values of DH, DK, DL (deviation of indices from integer values), indicating a better ORGX ORGY (namely XD=1288.4 YD=1262.8) than the one provided by the header (1265.47 1273.87) which results in XD=1279.6 YD=1268.0 (first line, starting with 0 0 0), and is going to lead to mis-indexing. This mis-indexing actually happens if one does not correct ORGX and ORGY when indexing based on a single frame only: if not corrected, in CORRECT.LP the ISa is 1.03 (instead of &gt;40), and the R-values are around 100%. <br /> <br /> We may experiment with increasing SPOT_RANGEs to see how many frames are required to make IDXREF find the correct ORGX ORGY. The result is that with SPOT_RANGE up to 1 500 the correct indexing is not found automatically, whereas with SPOT_RANGE=1 600 and higher the correct indexing is identified. It is therefore a good idea to use a large SPOT_RANGE to give IDXREF more information about the lattice, and to increase the chances of automatic identification of the correct indexing.<br /> <br /> It should be noted that even with a single frame used for indexing, the table gives a very clear indication that the provided ORGX ORGY is incorrect, and better values are suggested - but the user has to read and understand the table, and act accordingly!<br /> <br /> == Integration results, first pass ==<br /> <br /> While we were inspecting IDXREF.LP, the program has been busy processing the frames. We change to the INTEGRATE tab.<br /> <br /> [[File:thaumatin_ACA2014_integrate1a.png|left|600px|click to see larger versions]]<br /> [[File:thaumatin_ACA2014_integrate1b.png|left|600px|click to see larger versions]]<br /> &lt;br clear=all&gt;<br /> Particularly noteworthy are the following findings:<br /> # the distance decreases almost monotonically, so using the DISTANCE parameter for the REFINE(INTEGRATE) keyword is meaningful - it probably indicates that radiation damage increases the cell constants somewhat, which is compensated by the distance refinement.<br /> # the missetting angles (bottom plot) oscillate with an amplitude of up to 0.5°, so the orientation matrix or other aspects of the geometric description are inaccurate.<br /> <br /> == Checking the agreement of actual and predicted reflections ==<br /> <br /> After INTEGRATE, we man move to the Frame tab and load FRAME.cbf:<br /> <br /> [[File:Thau_ACA2014_FRAME.png|600px|click to see larger versions]]<br /> <br /> It is interesting to see that XDSGUI (which uses only the values from XDS.INP) places the low-resolution circle (green) around the header-specified ORGX, ORGY. The white circular area, however, indicates where XDS positioned the low-resolution exclusion after geometric refinement in IDXREF and INTEGRATE. In this case of particularly wrong header values, the two circels disagree strongly. If we would update XDS.INP with corrected ORGX, ORGY values, the circles would agree.<br /> <br /> The most important finding is, however, that the actual and predicted reflections match nicely. We may zoom in and inspect individual reflections. This will reveal that at large 2theta, the integration areas are elongated - this illustrates XDS&#039; internal geometric transformations that lead to undistorted reflections profiles, which is a requirement for high-accuracy profile fitting.<br /> <br /> [[File:Thau_ACA2014_FRAME_highres.png|600px|click to see larger versions]]<br /> <br /> == [[CORRECT]]: scaling and merging results (first pass) ==<br /> <br /> We move to the CORRECT tab and inspect the left-hand text. It may be helpful to use ctrl-- (control-minus) to lower the font size!<br /> <br /> The program reports about the point lattices it tried:<br /> <br /> SPACE-GROUP UNIT CELL CONSTANTS UNIQUE Rmeas COMPARED LATTICE-<br /> NUMBER a b c alpha beta gamma CHARACTER<br /> <br /> 5 81.8 81.8 150.0 90.0 90.0 90.0 3802 7.0 22754 10 mC<br /> 75 57.8 57.8 150.0 90.0 90.0 90.0 1896 11.1 24660 11 tP<br /> 89 57.8 57.8 150.0 90.0 90.0 90.0 1125 11.1 25431 11 tP<br /> * 21 81.8 81.8 150.0 90.0 90.0 90.0 2032 7.9 24524 13 oC<br /> 5 81.8 81.8 150.0 90.0 90.0 90.0 3802 7.0 22754 14 mC<br /> 1 57.8 57.8 150.0 90.0 90.0 90.0 6919 4.4 19637 31 aP<br /> 16 57.8 57.8 150.0 90.0 90.0 90.0 2095 11.1 24461 32 oP<br /> 3 57.8 57.8 150.0 90.0 90.0 90.0 3850 8.5 22706 35 mP<br /> 3 57.8 150.0 57.8 90.0 90.0 90.0 3749 7.8 22807 34 mP<br /> 1 57.8 57.8 150.0 90.0 90.0 90.0 6919 4.4 19637 44 aP<br /> <br /> <br /> ************ SELECTED SPACE GROUP AND UNIT CELL FOR THIS DATA SET ************<br /> <br /> SPACE_GROUP_NUMBER= 21<br /> UNIT_CELL_CONSTANTS= 81.77 81.78 149.99 90.000 90.000 90.000<br /> <br /> The true space group (92) was unfortunately not identified by the simple heuristic that XDS uses. This is because the Rmeas in the resolution range 10-5Å is 11.1% for the tetragonal lattice (space group 89), and this is more than 2 times (MAX_FAC_RMEAS) higher than the Rmeas in the triclinic lattice, 4.4%. The next highest symmetry is space group 21, and since its Rmeas is only 7.9%, it is chosen. If we go to the TOOLS menu / &quot;Further analysis&quot; / &quot;determine spacegroup with pointless&quot; we see in the console window (!) that [[pointless]] suggests either space group 92 (correct!) or 96 (enantiomorph of 92; equivalent at this stage).<br /> We thus have to teach XDS to use 92, by transferring this information into XDS.INP: we change<br /> <br /> JOB= XYCORR INIT COLSPOT DEFPIX INTEGRATE CORRECT<br /> ...<br /> SPACE_GROUP_NUMBER=0<br /> UNIT_CELL_CONSTANTS=70 80 90 90 90 90<br /> to<br /> JOB= CORRECT<br /> SPACE_GROUP_NUMBER=92<br /> UNIT_CELL_CONSTANTS=57.8 57.8 150.0 90.0 90.0 90.0<br /> press &quot;Save&quot; and &quot;Run XDS&quot;.<br /> <br /> This forces the correct space group, and the plots in the CORRECT tab then look like this: <br /> [[File:thaumatin_ACA2014_correct1a.png|left|600px|click to see larger versions]]<br /> [[File:thaumatin_ACA2014_correct1b.png|left|600px|click to see larger versions]]<br /> [[File:thaumatin_ACA2014_correct1c.png|left|600px|click to see larger versions]]<br /> <br /> &lt;br clear=all&gt;<br /> So the data are quite good to the edge of the detector. In principle we could have processed into the corners of the detector; this would have required <br /> TRUSTED_REGION=0 1.42<br /> but then the completeness would drop to zero at the highest resolution.<br /> <br /> == Further statistics from XDSSTAT ==<br /> <br /> After moving to the XDSSTAT tab, we click &quot;run xdsstat&quot; and obtain <br /> [[File:Thaumatin ACA2014 xdsstat1.png|600px|click to see larger versions]] <br /> (the upper two panels are left out here). <br /> <br /> We may also look at 2-dimensional representations of indicators mapped onto the detector surface, by using the &quot;view&quot; pull-down menu. However, none of them looks &quot;interesting&quot; enough to be reproduced here!<br /> <br /> What is interesting, however, is<br /> # the oscillating behaviour of the CORR plot (upper panel)<br /> # the slight rise of R_d (radiation damage indicator), from about 4% (left side) to about 6% (right side). This corresponds to a 50% increase in R_meas and indicates that the influence of radiation damage at the end of data collection is about as high as the overall level of all other errors. The argumentation goes like this: two uncorrelated sources of error, that both contribute to R_meas, have to be added by adding their variances. If their variances are equal this means that the total variance is doubled. Concerning R_meas, this means a sqrt(2)=1.4142-fold increase. Here we observe roughly 1.5-fold increase - even slightly more than sqrt(2). To get a more accurate idea of the amount of radiation damage relative to the other sources of error, one could fit a weighted least-squares line to the R_d plot, and base the calculation on that.<br /> <br /> == Improved processing (round 2) ==<br /> <br /> The TOOLS menu of XDSGUI / &quot;Saving and comparing good results&quot; offers &quot;backup files to ./save&quot;. We do this now - it is useful to go back to the current state in case any changes mess something up, or just for comparing files - there is a button &quot;compare CORRECT.LP with previous best&quot;. Please note that by default this uses the &quot;xdiff&quot; binary. If on your system this graphical tool is not installed, maybe the &quot;xxdiff&quot; or &quot;tkdiff&quot; tool is installed, or you could have one of the two [[XDSGUI#Dependencies|installed]]? In that case, just change the line below the button accordingly.<br /> <br /> Next, we click &quot;Optimizing data quality&quot; / &quot;copy latest geometry description over previous one&quot;. As can be seen in the line below (which may be edited!), this just overwrites the previous XPARM.XDS with the GXPARM.XDS that CORRECT has written: it has the space group information, and refined crystal setting and beam and spindle directions.<br /> <br /> Next, we go back to the XDS.INP tab and change<br /> JOB= DEFPIX INTEGRATE CORRECT ! DEFPIX is not really necessary but it doesn&#039;t hurt either<br /> and then klick &quot;Save&quot; and &quot;Run XDS&quot;. We may move to the INTEGRATE tab and see the plots evolving. When INTEGRATE is done, we may move to the CORRECT tab and see the plots. In this case they are so similar to the previous ones that they are not reproduce here. However, the statistics in the highest shells are slightly better (check for yourself!). What has also improved, is the CORR plot from XDSSTAT:<br /> <br /> [[File:Thaumatin_ACA2014_corr.png|600px|click to see larger versions]]<br /> <br /> It no longer oscillates, and the level is slightly higher. This shows that it is useful to<br /> mv GXPARM.XDS XPARM.XDS<br /> as an optimization step.<br /> <br /> == Solving the structure ==<br /> <br /> We run [[hkl2map]] and load XDS_ASCII.HKL - SHELXC understands the XDS format and picks up cell and space group.<br /> <br /> [[File:Thau_ACA2014_hkl2map1a.png|600px|click to see larger versions]]<br /> <br /> [[File:Thau_ACA2014_hkl2map1b.png|600px|click to see larger versions]]<br /> <br /> Next, we run SHELXC and find that the anomalous CC1/2 extends to the high resolution edge.<br /> <br /> [[File:Thau_ACA2014_hkl2map1c.png|600px|click to see larger versions]]<br /> <br /> Nevertheless, we accept the defaults (because the data are so good ...), let the GUI calculate the solvent content from the number of residues (207) <br /> <br /> [[File:Thau_ACA2014_hkl2map1d.png|600px|click to see larger versions]]<br /> <br /> [[File:Thau_ACA2014_hkl2map1f.png|600px|click to see larger versions]]<br /> <br /> ... and run SHELXD, searching for 15 sites. This is just an estimate; based on the &quot;Site occupancy versus Peak Number&quot; plot there are really 16 strong (sulfur) sites <br /> <br /> [[File:Thau_ACA2014_hkl2map1e.png|600px|click to see larger versions]]<br /> <br /> At this point, the structure can be considered solved: SHELXE built 197 residues, and it would be really easy to finish this up.</div> Karsten