Scaling with SCALA (or better: aimless): Difference between revisions

Jump to navigation Jump to search

Scaling with SCALA (or better: aimless) (view source)

Revision as of 22:54, 14 December 2015

269 bytes removed ,  14 December 2015
→‎Using truncate
(7 intermediate revisions by the same user not shown)
Line 1: Line 1:
== Using combat ==
== How to switch off scaling in CORRECT ==


If one wants to completely switch scaling off in CORRECT, use
MINIMUM_I/SIGMA=50
CORRECTIONS=
NBATCH=1
in [[XDS.INP]]. MINIMUM_I/SIGMA=50 prevents the sigmas from being adjusted. This also prevents outlier rejection and geometry refinement. Since the latter may be desired it would be most appropriate to use MINIMUM_I/SIGMA=50 only after [[Optimization]] (re-cycling of GXPARM.XDS to XPARM.XDS). Prevention of outlier rejection would be also obtained by using (e.g.) WFAC1=2. NBATCH=1 removes any basic scaling done by CORRECT (as referred to by [http://journals.iucr.org/d/issues/2012/05/00/tz5002/index.html Krojer and von Delft]).


An example script (contributed by Clemens Vonrhein) to convert XDS_ASCII.HKL into a [http://www.ccp4.ac.uk CCP4] multi-record MTZ file is:
A good option might be to make CORRECT only do MODULATION correction since this is not easily available in SCALA:
CORRECTIONS= MODULATION
This can also be combined with NBATCH=1.


#!/bin/csh -f
== Using pointless and aimless ==
#
# Scale data from XDS
#
combat hklin XDS_ASCII.HKL hklout junk.mtz << end-combat
MONITOR 10000
INPUT XDSASCII
ADDBATCH 200            # optional: modify batch numbers
TITLE  Test combat
NAME PROJECT myproject CRYSTAL mycrystal DATASET native
END
end-combat
sortmtz hklout junk_sort.mtz  << end-sort
H K L M/ISYM BATCH I SIGI
junk.mtz
end-sort
scala hklin junk_sort.mtz hklout junk_scaled.mtz << end-scala
# scale using default parameters
run 1 batch 1 to 1000
scales rotation spacing 5 bfactor off
# if xds wouldn't have already scaled the data, e.g.
# scales rotation spacing 5  secondary 4  bfactor on  brotation  spacing 20
# would be more appropriate (this is what "quickscale" in iMosflm does)
anomalous on
end-scala
truncate  hklin  junk_scaled.mtz \
          hklout junk_truncated.mtz <<end-trunc
anomalous yes
nresidue  1049
labout  F=FP SIGF=SIGFP DANO=DANO_sulf SIGDANO=SIGDANO_sulf
end-trunc


In this way one gets the well-known outputs of [http://www.ccp4.ac.uk/html/scala.html SCALA] and [http://www.ccp4.ac.uk/html/truncate.html TRUNCATE], and the MTZ file junk_truncated.mtz .
[[pointless]] produces a multi-record MTZ file (suitable for SCALA/aimless and TRUNCATE) from XDS_ASCII.HKL, or the XSCALE reflection output file (see also "[ftp://ftp.ccp4.ac.uk/ccp4/6.0.2/prerelease/runningpointlessandscala.pdf Running Pointless and Scala]"). Note that this may use [[XDS]] for some or most of the scaling.
(see above for how to switch off scaling in [[XDS.INP]], partly or fully).


== Using pointless ==
pointless -copy xdsin XDS_ASCII.HKL hklout XDS_ASCII.mtz
aimless hklin XDS_ASCII.mtz hklout junk_scaled.mtz | tee aimless.log


An alternative way of going from a XDS formatted reflection file to a multi-record MTZ file (suitable for SCALA and TRUNCATE, but not directly for structure solution which requires single-record MTZ files) is via the [[pointless]] program
Obviously, the aimless run can be modified using suitable options.
(see also "[ftp://ftp.ccp4.ac.uk/ccp4/6.0.2/prerelease/runningpointlessandscala.pdf Running Pointless and Scala]"). Note that this way also means that [[XDS]] is used to do most of the scaling (because XDS_ASCII.HKL is written by [[CORRECT]]) but this could be switched off in [[XDS.INP]] (partly or fully).
It's very easy to do:


pointless -copy xdsin XDS_ASCII.HKL hklout XDS_ASCII.mtz
See below for how to run TRUNCATE.


The advantage of this route is, that the batch headers of the resulting MTZ file are correct - [http://www.ccp4.ac.uk/dist/html/combat.html COMBAT] doesn't seem to write complete and accurate headers. If one then needs to edit the batch information in this file, the CCP4 program [http://www.ccp4.ac.uk/dist/html/rebatch.html REBATCH] can be used (optional):
If one then needs to edit the batch information in this file, the CCP4 program [http://www.ccp4.ac.uk/dist/html/rebatch.html REBATCH] can be used (optional):


  rebatch hklin XDS_ASCII.mtz hklout XDS_ASCII_rebatch.mtz <<end_ip
  rebatch hklin XDS_ASCII.mtz hklout XDS_ASCII_rebatch.mtz <<end_ip
Line 62: Line 37:


Note: another program to work with those multi-record MTZ files is [http://www.ccp4.ac.uk/dist/html/reindex.html REINDEX].
Note: another program to work with those multi-record MTZ files is [http://www.ccp4.ac.uk/dist/html/reindex.html REINDEX].
== Using truncate ==
#!/bin/csh -f
truncate  hklin  junk_scaled.mtz \
          hklout junk_truncated.mtz <<end-trunc
anomalous yes
nresidue  1049
labout  F=FP SIGF=SIGFP DANO=DANO_sulf SIGDANO=SIGDANO_sulf
end-trunc
In this way one gets the well-known output of [http://www.ccp4.ac.uk/html/truncate.html TRUNCATE], and the MTZ file junk_truncated.mtz .


== See also ==
== See also ==


[[XDSCONV]] to get the MTZ file directly (i.e. circumventing SCALA)
[[XDSCONV]] to get the MTZ file directly (i.e. the usual way)
Kay
Bureaucrats
2,651

edits

Retrieved from "https://wiki.uni-konstanz.de/xds/index.php/Special:MobileDiff/2195...3271"
Cookies help us deliver our services. By using our services, you agree to our use of cookies.
More information

Navigation menu