-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathtargeted_scaling_guide.html
154 lines (126 loc) · 9.01 KB
/
targeted_scaling_guide.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="X-UA-Compatible" content="IE=Edge" />
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<title>Scaling against a reference dataset with DIALS — DIALS documentation</title>
<link rel="stylesheet" href="_static/alabaster.css" type="text/css" />
<link rel="stylesheet" href="_static/pygments.css" type="text/css" />
<link rel="stylesheet" type="text/css" href="_static/dials-styles.css" />
<script type="text/javascript" id="documentation_options" data-url_root="./" src="_static/documentation_options.js"></script>
<script type="text/javascript" src="_static/jquery.js"></script>
<script type="text/javascript" src="_static/underscore.js"></script>
<script type="text/javascript" src="_static/doctools.js"></script>
<script type="text/javascript" src="_static/language_data.js"></script>
<script async="async" type="text/javascript" src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.2/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
<link rel="shortcut icon" href="_static/favicon.ico"/>
<link rel="author" title="About these documents" href="about.html" />
<link rel="index" title="Index" href="genindex.html" />
<link rel="search" title="Search" href="search.html" />
<link rel="next" title="Workshops" href="workshops/index.html" />
<link rel="prev" title="User guide for scaling data with DIALS" href="dials_scale_user_guide.html" />
<link rel="stylesheet" href="_static/custom.css" type="text/css" />
<meta name="viewport" content="width=device-width, initial-scale=0.9, maximum-scale=0.9" />
</head><body>
<div class="logoheader container">
<a href="index.html">
<img class="logoheader" alt="DIALS" src="_static/dials_header.png" />
</a>
</div>
<div class="document">
<div class="documentwrapper">
<div class="bodywrapper">
<div class="body" role="main">
<div class="section" id="scaling-against-a-reference-dataset-with-dials">
<h1>Scaling against a reference dataset with DIALS<a class="headerlink" href="#scaling-against-a-reference-dataset-with-dials" title="Permalink to this headline">¶</a></h1>
<p>This document provides a guide on how to scale a dataset against a reference
dataset, referred to as the target dataset.
The target dataset can be either a dataset scaled with dials.scale, or an mtz
file containing a scaled dataset.
Any number of datasets can be scaled against the reference dataset at once,
giving the same scaling result as if each dataset were scaled as an independent
job. The only difference is that all scaled datasets would be output in one
scaled.refl and scaled.expt, which may be more or less convenient
for further processing.</p>
<div class="section" id="scaling-against-a-dials-reference-dataset">
<h2>Scaling against a dials reference dataset<a class="headerlink" href="#scaling-against-a-dials-reference-dataset" title="Permalink to this headline">¶</a></h2>
<p>In this example, reference.refl and reference.expt are
from a dataset that has already been scaled with dials.scale. To scale another
dataset (datafiles integrated.refl, integrated.expt) against this
target/reference, one should use the following command:</p>
<div class="highlight-default notranslate"><div class="highlight"><pre><span></span><span class="n">dials</span><span class="o">.</span><span class="n">scale</span> <span class="n">only_target</span><span class="o">=</span><span class="kc">True</span> <span class="n">integrated</span><span class="o">.</span><span class="n">refl</span> <span class="n">integrated</span><span class="o">.</span><span class="n">expt</span> <span class="n">reference</span><span class="o">.</span><span class="n">refl</span> <span class="n">reference</span><span class="o">.</span><span class="n">expt</span>
</pre></div>
</div>
<p>This will scale the intensities of the dataset to agree as closely as possible
with the intensities of the reference dataset, and save the scaled dataset to
scaled.refl, scaled.expt (the reference files are unchanged).
The <code class="samp docutils literal notranslate"><span class="pre">only_target=True</span></code> command is important, else all the data will be
scaled together and output in a joint output file.</p>
</div>
<div class="section" id="scaling-against-a-reference-mtz-file">
<h2>Scaling against a reference mtz file<a class="headerlink" href="#scaling-against-a-reference-mtz-file" title="Permalink to this headline">¶</a></h2>
<p>In this case, it is assumed that the intensity and variance columns of the mtz
file have already been scaled. Targeted scaling would be run with the following
command:</p>
<div class="highlight-default notranslate"><div class="highlight"><pre><span></span><span class="n">dials</span><span class="o">.</span><span class="n">scale</span> <span class="n">integrated</span><span class="o">.</span><span class="n">refl</span> <span class="n">integrated</span><span class="o">.</span><span class="n">expt</span> <span class="n">target_mtz</span><span class="o">=</span><span class="n">scaled</span><span class="o">.</span><span class="n">mtz</span>
</pre></div>
</div>
<p>The targeted scaling algorithm is the same regardless of the target datafile type,
likewise the scaled dataset will be saved to scaled.refl and scaled.expt.</p>
</div>
<div class="section" id="general-considerations-for-suitable-options">
<h2>General considerations for suitable options<a class="headerlink" href="#general-considerations-for-suitable-options" title="Permalink to this headline">¶</a></h2>
<p>A common use case for scaling against a reference is to scale thin-wedge
datasets against a high quality full-sweep dataset. To give the best scaling, it
may be necessary to manually set the scaling model parameters: for more details
see the <em>In-depth guide to scaling options in DIALS</em>.
In the case of very thin wedge/stills datasets, or depending on the scientific question under investigation, it may be
suitable to set <code class="samp docutils literal notranslate"><span class="pre">model=KB</span></code>, to give a single global scale and relative B-factor
to each dataset. However, if significant intensity variation/decay is present in each
measurement, it may be best to use <code class="samp docutils literal notranslate"><span class="pre">model=physical</span></code>, setting <code class="samp docutils literal notranslate"><span class="pre">absorption_term=False</span></code>
and specifying values for <code class="samp docutils literal notranslate"><span class="pre">scale_interval</span></code> and <code class="samp docutils literal notranslate"><span class="pre">decay_interval</span></code>.</p>
</div>
</div>
</div>
</div>
</div>
<div class="sphinxsidebar" role="navigation" aria-label="main navigation">
<div class="sphinxsidebarwrapper"><!--<h3>Navigation</h3>-->
<ul class="current">
<li class="toctree-l1"><a class="reference internal" href="index.html">Home</a></li>
<li class="toctree-l1"><a class="reference internal" href="about.html">About</a></li>
<li class="toctree-l1"><a class="reference internal" href="installation.html">Installation</a></li>
<li class="toctree-l1"><a class="reference internal" href="documentation/index.html">Documentation</a></li>
<li class="toctree-l1"><a class="reference internal" href="documentation/tutorials/index.html">Tutorials</a></li>
<li class="toctree-l1 current"><a class="reference internal" href="howto.html">How-to</a></li>
<li class="toctree-l1"><a class="reference internal" href="workshops/index.html">Workshops</a></li>
<li class="toctree-l1"><a class="reference internal" href="publications.html">Publications</a></li>
<li class="toctree-l1"><a class="reference internal" href="links.html">Links</a></li>
<li class="toctree-l1"><a class="reference internal" href="license.html">License</a></li>
</ul>
</div>
</div>
<div class="clearer"></div>
</div>
<div class="footer container">
<a href="https://www.diamond.ac.uk/"><img class="logofooter" alt="Diamond" src="_static/diamond_logo.png" /></a>
<a href="https://www.ccp4.ac.uk/"><img class="logofooter" alt="CCP4" src="_static/CCP4-logo-plain.png" /></a>
<a href="https://www.lbl.gov/"><img class="logofooter" alt="LBL" src="_static/LBL-logo-wide.jpeg" /></a>
<a href="https://www.stfc.ac.uk/"><img class="logofooter" alt="STFC" src="_static/STFC_logo.png" /></a>
</div>
<script type="text/javascript">
$(document).ready(function() {
$(".toggle > *").hide();
$(".toggle .header").show();
$(".toggle .header").click(function() {
$(this).parent().children().not(".header").toggle(400);
$(this).parent().children(".header").toggleClass("open");
})
});
</script>
<div class="footer">
©2019, Diamond Light Source, Lawrence Berkeley National Laboratory and STFC.
</div>
</body>
</html>