Any Commands
You may also visit the explanation of
internal commands: aliases, flowcontrol etc.
List of commands
A
B
C
D
E
F
G
H
I
J
K
L
M
N
O
P
Q
R
S
T
U
V
W
X
Y
Z
abort
Syntax: abort state
Default: abort off
If state=off,
warnings are not fatal.
If state=on a warning
will abort the program.
activate
Syntax: activate n
Activates a subplot area.
The subplots are numbered (from top to bottom) 1, 2 and 3.
In each subplot, several settings (plottype, horizontal axis, vertical axis,
colourtype, plotlimits)
are stored. The pictures in each subplot can be redrawn with
r.
Related command:
deactivate
allexpressions
Re-evaluates all expressions. The dependency
of expressions to other expressions is handled. But an expression may also
be dependent on one of the other variables. In general, these variables
do not change, with a few exceptions:
equiv2, imean or
smean. Also
the overwrite command can change a variable.
Example:
- an expression depends on the variable imean
- you change a reflection filter
- you use delsig to re-evaluate
imean.
Related command:
autoexpression
allow
Set reflection filters. See the flags page.
alwaysforbid
Set reflection filters.
Flags set by alwaysforbid are always forbidden.
The filter settings for these flags are not influenced by
the filter commands
allow,
forbid and
require.
Alwaysforbid Example
bad uniformity reflections are always forbidden. After the
command allow all they remain forbidden.
Use a minus sign to remove flags from the allwaysforbid mask.
See the flags page.
asym
Shows the reciprocal lattice in layers. You enter the name of the
reciprocal axis (h, k or l) along which the layers are printed.
Then you enter the layer number. Use quit to
return to the command prompt.
The following symbols are used:
# |
reflection not available |
. |
reflection not measured, but an equivalent one is available |
a |
reflection one time measured |
b |
reflection two times measured |
c | reflection three times measured |
etc |
reflection n times measured |
autoexpression
Syntax: autoexpression state
If on, all espressions will be re-evaluated at
some stages in the program (after reading a new rmatrix
using rmat, after changing an
existing expression, after reading another
dataset). The expressions can also be re-evaluated
with allexpressions.
autoincidence
Syntax: autoincidence
Calculates a incidence coefficient.
1-cosd(incid).
- store intensities
- execute expression 9 1-cosd(incid)
- execute histoline on
- set incidencecoefficient to 0.0
- divide intensities by 1+incidencecoefficient*(expr9)
- execute icr expr9
- determine (histoline) coefficients ( A*expr9+B )
- add A/B to the incidencecoefficient
- repeat this loop until A is less than 0.001
- restore intensities
automsa
Syntax: automsa
Calculate the msa factor using rint.
The final factor will show a constant χ² for all intensity
intervals.
Related commands:
kfactor,
msa,
msapower,
sigmafactor.
autoshell
Syntax: autoshell state
Default: autoshell off
Thetamax is the high value of
limit theta.
- on
If Thetamax is less than 10.0° calculate
nshell boundaries from 0.0 to 10.0
else if Thetamax is less than 15.0° calculate
nshell boundaries from 0.0 to 15.0
else if Thetamax is less than 27.5° calculate
nshell boundaries from 0.0 to 27.5
else calculate nshell boundaries from 0.0 to
Thetamax
If the boundaries are calculated with the values 10.0, 15.0 or 27.5
the resulting number of shells can be less than the number specified
with nshells (because there are no reflections
between Thetamax and the top boundary).
- off
Calculate nshell boundaries from 0.0 to
Thetamax
badmother
Syntax: badmother
Set the BADMOTHER flag
for equivalent reflections if
at least one of the equivalent reflections is forbidden.
Use the filter commands
allow,
forbid and
require to enable/disable the BADMOTHER flag once
it is set.
bfactor
Syntax: bfactor
- if bfactors have been applied, undo
- otherwise apply the bfactors
-
Every experiment should have a bfactor
(see bfactorread and bfactorexp)
- f =
exp(bfactor*sin(theta)²/lambda²)
- Multiply
intensity,
sigmaorig,
eval14i and
eval14s
by f
- Recalculate the sigma's of all reflections.
See the Mean and Sigma section
on the main page.
bfactorexp
Syntax: bfactorexp expname bfactor
Syntax: bfactorexp all bfactor
- If the bfactor has been applied, undo the bfactor.
- set the bfactor for experiment expname.
If expname=all the bfactor for each experiment
will be set to the same bfactor.
- If the bfactor was applied, apply again.
bfactorread
Syntax: bfactorread filename
- If the bfactor has been applied, undo the bfactor.
- read filename, line by line
- lines starting with ! are ignored
- otherwise, lines should contain experimentname bfactor
- If the bfactor was applied, apply again.
Related command:
write bfactor
bgqual
Syntax: bgqual f
Sets BADBG flag on reflections with
background quality > f.
The background quality is determined by eval14 and stored in
variable fombg.
The default for f is extracted from the datafile.
bravais
Change the bravais type. The bravais type is extracted from the
input file, but you may change it.
Syntax: bravais code
Code is one of ABCFIPR.
pointgroup overview
Related commands:
laue,
ps.
capillaryaxis
Syntax: capillaryaxis f1 f2 f3
Default: capillaryaxis 0 0 1
Defines the capillaryaxis.
f1 f2 f3 are coordinates in the
laboratory system.
Note, this definition does not use any goniostat
information. So think of the crystal as it is mounted on the goniometerhead,
and put the goniometer head with the crystal in front of you. Then
capillaryaxis 0 0 1 will point to the zenith.
- capillaryin
contains the angle between the primary beam (inray) and the
capillaryaxis.
- capillaryout
contains the angle between the diffracted beam (outray) and the
capillaryaxis.
cchalf
Calculate CC1/2 with the current two
groups of split intensities. These two groups are build when
cchalfsplit
Split and merge all intensities into 2 groups. These groups are used
to calculate CC1/2. See glossary on CC1/2
for more information.
cell
Syntax: cell a b c alpha beta gamma orx ory ora
Extracts cell parameters (a b c alpha beta gamma)
and orientation (orx ory ora) form the rmatrix.
You may then change these parameters. A new rmat will be build.
colour
Syntax: colour item colourname
Sets a colour. Used in the
histcount,
histo,
hkl,
icr,
plot,
plot2d,
plotgrid,
plotray,
plotrefl commands.
main |
darkgoldenrod |
background main window |
fg |
white |
hkl, go |
bg |
black |
all |
inray |
green |
plotray |
histobar1 |
grey40 |
histo, histcount |
histobar2 |
white |
histo, histcount |
histoline1 |
darkgreen |
lines through histogram points of histo
and histcount |
histoline2 |
green |
calculated line of histo and
histcount |
label |
yellow |
numbers along axes |
text |
yellow |
labels along axes, title, legenda |
ticks |
grey30 |
icr,
plot,
plotc,
plotd,
plotray
|
strong |
red |
colourtype weak |
weak |
blue |
colourtype weak |
negative |
blue |
colourtype weak |
renninger |
cyan |
icr, plot
(if renninger on)
|
var1 through var11 |
11 colours from red to blue |
colourtype var |
Use colourvar or
colourscheme if you want to change all
the colours of var1 - var11.
See also colourbw and
colourinvert
colourbw
Syntax: colourbw state
Default: colourbw off
If set to on, all colours will be set to white,
except colour bg, which is set to black.
See also colourinvert
colourinvert
Syntax: colourinvert state
Default: colourinvert off
If set to on, all colours will be inverted.
You may set colourinvert on before creating
snapshots with ps.
See also colourbw
colourscheme
Syntax: colourscheme bw/rainbow/bry/blueyellow
Default: colourscheme rainbow
Sets the colours of var1 to var11. These colours are used by the
colourtype var command.
- bw. From black to white
- rainbow. From darkred to darkblue
- bry. From darkblue via darkred to yellow
- blueyellow. From blue to yellow
colourtype
Syntax: colourtype weak/fixed/sign/var/flag
Default: colourtype weak
- weak
three colours: weak reflections are drawn in
colour weak, negative reflections
in colour negative and strong reflections
in colour strong.
The weak and negative flags are set with sigma.
- fixed
reflections are drawn in colour fg.
- sign varname
The value of varname is used to set
two colours: reflections with negative value are drawn in
colour weak and
reflections with a positive value
in colour strong.
- var varname
The value of varname is used to set
the colour. The 11 (or less if you used
ncolourvar)
colours (use colourvar or
colourscheme to
change the colours) are linearly mapped between the minimum and
maximum value of varname, but you can change
these limits with varlimit.
Use dotsizetype var
to change the dotsize as function of varname.
- flag flagnamelist
flagnamelist is a set of reflectionflags, combined
with + signs.
If one of the reflection flags is set, use
colour strong otherwise use
colour weak.
See the flags page for all possible flags.
Use flagfilter to obtain a (shorter) list of
allowed flags.
colourvar
Syntax: colourvar 11 colournames
Sets the colours of var1 to var11 (or less if you changed the number
of colours with ncolourvar). These colours are
used by the colourtype var command.
You are prompted for each colour. The answers can be:
- [return] Accept the proposed name
- . (a dot) Accept the proposed name
- - (a minus) Use the previous colour name.
Colourvar Examples
- to select a 2-colour scheme, you may use:
colourvar red - - - - - blue - - - -
- to change colour 5, you may use:
colourvar . . . . red . . . . . .
(but colour var5 red is much easier.
Related commands: colour,
colourscheme,
colourtype var,
ncolourvar.
common
Syntax: common
For every experiment print
- number of reflections
- number of unique reflections
- number of unique reflections common with the other experiments
compareexperiments
Syntax: compareexperiments e1 e2
Compares two experiments e1 e2 and prints the
differences (rotation axis, rotaxnumber, goniostat angles, goniostat
offsets, swing, detector distance, rotation increment, first and last
framenumbers, detzero, detrot, intensity correction, rotation start,
rotation end).
conevec
Syntax: conevec hkl/uvw/xyz f1 f2 f3
Default: conevec xyz 0 0 1
Defines the cone-vector. The
variable
cone
contains the angle between the diffracting c-vector and this
cone-vector.
- conevec hkl. Specify 3 reflection indices
- conevec uvw. Specify 3 indices in terms
of the direct lattice.
- conevec xyz. Specify coordinates in the
laboratory system.
Note, conevec is defined with respect to the crystal and does
not take into account the goniostat.
So think of the crystal as it is mounted on the goniometerhead,
and put the goniometer head with the crystal in front of you. Then
conevec xyz 0 0 1 will point to the zenith
and
conevec hkl 1 0 0 will coincide with a*.
copy
Syntax: copy var1 var2
Copies the values from var1 to var2.
Runtime variables are not allowed for
var2 (The src column in
variables
can be anything but c.
copy can be dangerous. It can completely ruin
the database.
Another approach is:
- expression 7
var1
- overwrite
var2
7
The following variables only contain data. They are not used by
any, unless you define expressions depending on them.
You can use these variables as storage for something else.
bga bgb bgc corbg corbox corpeak edgefrac
fombg fombox fompeak sigmabad volume
Related commands:
expression,
overwrite
correct
Syntax: correct none/normal/full/debug
Default: correct normal
Allows command and option corrections. Commands and options
are only corrected if input originates from the keyboard, and output
is sent to the screen.
- none
No correction will be applied
- normal
Commands and options will be corrected.
- full
As before plus possible corrections are printed
- debug
As before plus all available commands/options are printed
cos
You enter h,k,l and any calculates the
direction cosines
for this reflection and its qivalents. Use to debug the program.
coscell
Syntax: coscell n
Default: coscell 1
Direction cosines
describe the angles between
incoming-xray/reciprocal-axes and outgoing-xray/reciprocal-axes.
The coscell command defines which reciprocal axes will be used.
n is the rmat number.
You may set it to an interfering
cell, and the direction cosines will be calculated relative to this
other cell.
If you set n to zero, the direction cosines for
each hkl are calculated with respect to the corresponding cell
(for example: main reflections use rmat 1, interfering reflections use
rmat2).
See also the Reflection output section on the
main page.
cosid
Syntax: cosid n
Calculate and print the
direction cosines for one reflection with
reflection number n.
cosines
Syntax: cosines state
Default: cosines on
If set to on,
direction cosines
will be written to the various reflection files.
See also the Reflection output section on the
main page.
deactivate
Syntax: deactivate n
Deactivates a subplot area.
The subplots are numbered (from top to bottom) 1, 2 and 3.
Subplot 1 can not be deactivated.
Related command:
activate
delsig
Calculates and stores the variables
imean,
smean,
smean,
delsig,
iplus,
splus,
imin and
smin
for all reflections using the
current filter settings.
Imean = ∑
(v·I) /
∑(v)
with v selected according the
weightint
delta/sigma = (I−Imean)/σ
With colourtype sign delsig the colour of
reflections is dependent on the sign of delta/sigma.
You may change the filter settings (for example
limit delsig -3 3 no), but the stored delta/sigma's
will not change (unless you call delsig again).
deoverlap
Syntax: deoverlap f
This command modifies intensities of overlapping reflections where the
main reflection is equivalent to the overlapping one. The new intensity
will be multiplied by f.
If ia and ib are the average
intensities of the two twin lattices, then
f=ia/(ia+ib).
If g is the ratio of the two lattices
(i.e. ia/ib), then
then f=g/(1+g).
This command will only run if
allow OVERLAPSUM is set.
Related command:
listrefl overequiv.
detectorsize
Syntax: detectorsize left right bottom top
Sets the dimensions of the detector in the detectorwindow. The defaults
are extracted from the input file. The detector window is used for the
commands plotgrid, plotrefl
and
plot2d.
dict
Syntax: dict state
Default: dict off
If state=on the output of
flagall,
flagfilter,
list exp,
list intensity,
list var,
minmax and
rshell
will contain a list of
keyword/value pairs. The keywords start with dictname.
This output can be used to parse the output of any
into another program.
Also, if state=on, the
command read will output the EVAL15PROCESS
commands from the y-file as !!eval15process-xxxx.
dictname
Syntax: dictname label
Default: dictname !!any-
Sets the leading part of keywords of dictionary output
(if dict=on).
dotsize
Syntax: dotsize n
Default: dotsize 1
Sets the dotsize (in pixels). This value is used for the drawing
of reflections for
dotsizetype fixed.
dotsizemax
Syntax: dotsizemax n
Default: dotsize 10
Sets the maximum dotsize for
dotsizetype var.
dotsizetype
Syntax: dotsizetype fixed/var
Default: dotsizetype fixed
- fixed
reflections are drawn with size specified by
dotsize.
- var varname
The value of varname is used to set
the size of the dots. The sizes range from to 1 to
dotsizemax
and are linearly mapped between the minimum and
maximum value of varname, but you can change
these limits with varlimit.
Use colourtype var
to change the colour as function of varname.
drx
Syntax: drx n
Creates a drxfile. From all
reflections passing the reflection filters,
a random subset of
n reflections will be selected.
You can specify the filename
with the drxfile command.
The drxfile can be used by DirAx to re-investigate the current lattice.
Note: Dirax will only use 1000 reflections.
drxfile
Syntax: drxfile name
Sets the name of the
drxfile
which can be created with the drx command.
dump
Syntax: dump
Default: dump off
Toggles the dump flag. If set, any listing of reflections
will print all the information of each reflection. This will generate
a lot of output. This command is added for debugging purposes.
dumpid
Syntax: dumpid n
Print all
variables
of reflection n in the format
dump=on.
The reflection constants will also be printed if
output=full.
end
Exit program.
euler
Syntax: euler state
Default: euler off
Only useable with a kappa goniostat. If
state=on, the kappa angles
in the list exp command will be converted to
eulerian angles.
evpylog
Syntax: evpylog text
Write text to evpy.log if that
file exists in the current directory or up to three levels up in the
directory tree.
exit
Exit program.
experiment
Enable/Disable experiments. Experiments are complete scans. Experiments
are divided into sets.
Syntax: experiment all/none/next/prev/curr (+/-)n ....
Use a plus sign to allow experiments and a minus sign to forbid them.
You may enter more selections in one line of input.
Use nop or ; to
signal the end of options if you want to put more commands
on a single line.
all |
enable all experiments |
none |
disable all experiments |
invert |
toggle all experiments |
curr |
evaluates to the last selected experiment |
next |
evaluates to the next experiment |
prev |
evaluates to the previous experiment |
last |
evaluates to the last experiment |
n |
sets/unsets experiment n |
Some examples are given in the set section.
(experiment is one of the
reflection filter options).
Use list exp to
print the names and status of all experiments.
You could also use limit expnr, but the
experiment command is moch more flexible.
Related command: set.
expression
Syntax: expression n formula.
Defines the formula for expression n
- There are 9 variables called expr1 to
expr9. The names of these variables can be
changed using expressionalias.
- In the formula the names of
variables
can be used.
- The syntax of expressions and available functions are explained
elsewhere.
- Use list expression to print the definitions.
- Various constants are available in expressions:
Expression Examples
-
the total shift of reflections can be defined with:
expression 1 sqrt(shifthor^2+shiftver^2)
To print the shifts, use
field add expr1 and
listrefl refl. Or to
see how the shifts change in time you may use
plot expr1 nr.
The variable shiftmm
already exists, so you can
verify this example with:
plot expr1 shiftmm.
- The distance between the impact and the centre of the detector can be
calculated:
expression 1 sqrt(hor^2+ver^2)
- The incidence angle is:
expression 2 atand(expr1/dist)
The variable incid
already exists, so you can
verify this example with:
plot expr2 incid.
- Use only reflections with odd h-indices
expression 3 odd(h)
limit expr3 1 1 yes
- A pixel (at coordinates 12.6 8.7) in experiment 5 is suspect.
How to disable all reflections within 2 mm of this pixel.
Start with calculating the distance from the suspect pixel:
expression 1 sqrt((12.6-hor)^2+(8.7-ver)^2)
limit expr1 0 2 no
will remove the suspect
reflections, but also reflections from the other experiments. So we
add a large number to expression 1 if the experiment is not 5:
expression 2 1000.0*abs(expnr-5)
expression 3 expr1+expr2
limit expr3 0 2.0 no
Another approach: the lt(f,g) function returns 1.0 if f<g and
otherwise 0.0. So
expression 2 lt(expr1,2.0)
will mark the suspect pixels with a value of 1.0
The function eq(f,g,v,w) returns v if f=g, otherwise it returns w.
We can use this function to limit ourselves to experiment 5.
expression 3 eq(expnr,5,expr2,0.0)
Expression 3 evaluates to expression 2 if the experiment number=5, otherwise
it evaluates to zero. And expression 2 evaluates to 1.0 only if the pixel is
suspect.
Finally, use the limit command to remove the suspect pixels:
limit expr3 0 0.1 no
Expression 3 can be written in one line by subsequent substitution
of expr2 and expr1.
expression 3 eq(expnr,5,lt(sqrt((12.6-hor)^2+(8.7-ver)^2),2.0),0.0)
This corresponds to the following Fortran code:
if(expnr.eq.5) then
if(sqrt( (12.6-hor)**2+(8.7-ver)**2).lt.2.0) then
result = 1.0
else
result = 0.0
endif
else
result = 0.0
endif
- Calculate abs(I+-I-)/√(σ+²+σ-²)
but only for reflections where either I+ or I-
unequals zero.
For all reflections (without the test) this is:
expression 1 abs(iplus-imin)/sqrt(splus^2+smin^2)
If both I+ and I- should be nonzero we could multiply
both variables. The result is zero if either one is zero:
expression 2 ne(iplus*imin,0.0)
This will evaluate to 1.0 if both variables are nonzero.
To test whether one of them is nonzero, calculate the absolute value of both
variables:
expression 3 max(abs(iplus),abs(imin))
and check whether the maximum of these two numbers is nonzero:
expression 4 ne(expr3,0.0)
To combine expression 1 and 4, we will use the ne
function with four arguments:
expression 4 ne(expr3,0.0,expr1,0.0)
and in one line:
expression 4 ne(max(abs(iplus),abs(imin)),0.0,abs(iplus-imin)/sqrt(splus^2+smin^2),0.0)
Related command: allexpressions.
expressionalias
Syntax: expressionalias n s.
Changes the name of expression
n into s.
The default names are exprn.
Syntax: extract filename
Write reflection information to a file. Use the program
extractshoe to
extract reflectionboxes for the reflections in this file from the whole set
of reflectionboxfiles. As long as the extractfile is open, reflections
printed to the screen (using various commands as hkl,
listrefl) will be written to the file.
Close the file with noextract.
Syntax: extractsubset n filename
Write reflection information for a (random) subset of n
reflections passing the reflection filters
to a file. You can use
the program extractshoe to create a new
reflectionboxfile from the whole set of reflectionboxfiles.
fcffilter
Syntax: fcffilter args
Default: fcffilter *
Set a filter on the status of reflections read by
readfcf and readicalc.
The reflection status is read from the fcf file in
the _refln_observed_status field.
The filter consists of one or more characters. If
the reflection status is listed in the filter, the reflection will
be accepted. Use * to accept all reflections.
field
Syntax: field add/list/remove/set args
Change the layout of the reflections printed to the screen.
add |
variable(s) |
adds variable(s) to the set |
list |
|
show the current set |
remove |
variable(s) |
remove variable(s) from the set |
set |
fieldnumber
variable |
puts a variable at a specific position |
For floating point data, the number of decimals for each variable can be
set with precision.
The final column contains a description of the reflection
flags. This column
may be removed with fieldflags off.
Field Examples
field list
1 set
2 refl
3 h
4 k
5 l
6 intensity
7 sigma
8 ioversig
9 delsig
hkl 0 -2 0
multiplicity=2
hkl= 0 -2 0 hb=9002 im=1 Theta=6.796 Reso=3.003 MaxInt=757.873 Mean=734.858
(mean=734.858 sigma=8.705 ioversig=84.422 <ioversig>=48.734) Rsym=0.015 Rmeas=0.018 Rpim=0.011 Chi2=0.951 n=3
nr boxfil set refl h k l intensity sigma ioversig delsig
1128 s01f002: 2 125 0 -2 0 745.87 15.17 49.16 0.73 +
1303 s01f002: 2 300 0 2 0 757.87 15.42 49.15 1.49 -EDGEROT
1476 s02f001: 3 168 0 2 0 740.59 15.27 48.51 0.38 +
1578 s02f001: 3 270 0 -2 0 718.42 14.80 48.53 -1.11 +
3 average: 734.96 15.08 48.73 -0.00
field add hor ver rot
hkl 0 -2 0
multiplicity=2
hkl= 0 -2 0 hb=9002 im=1 Theta=6.796 Reso=3.003 MaxInt=757.873 Mean=734.858
(mean=734.858 sigma=8.705 ioversig=84.422 <ioversig>=48.734) Rsym=0.015 Rmeas=0.018 Rpim=0.011 Chi2=0.951 n=3
nr boxfil set refl h k l intensity sigma ioversig delsig hor ver rot
1128 s01f002: 2 125 0 -2 0 745.87 15.17 49.16 0.73 -16.67 9.37 82.63 +
1303 s01f002: 2 300 0 2 0 757.87 15.42 49.15 1.49 -4.55 -7.01 106.22 -EDGEROT
1476 s02f001: 3 168 0 2 0 740.59 15.27 48.51 0.38 0.74 0.85 -178.23 +
1578 s02f001: 3 270 0 -2 0 718.42 14.80 48.53 -1.11 21.28 0.85 -164.63 +
3 average: 734.96 15.08 48.73 -0.00 1.78 3.69 -86.74
fieldflags
Syntax: fieldflags state
Default: fieldflags on
Change the layout of the reflections printed to the screen.
(hkl,
hklm,
listrefl)
The contents of the table is controlled by the
field command, except for the final column containing
the reflection flags.
Use fieldflags off to remove this
column.
fillgap
Syntax: fillgap state
Default: fillgap on
While reading reflections from the inputfile, the imagename for each
reflection is stored in a list. The variable
frame
is the indexpointer of the reflection in this list.
If a gap in the images occurs, fillgap controls
whether the intermediate images will be inserted into the list.
If reflections occur on all images, the setting of
fillgap is of no importance. If some images do
not contain reflections (small
unit cells or very fine sliced images) and
fillgap=off, the frame number of the reflections
will be wrong. On the other hand, if the first image of a set does
not start with number 1, an articifial gap is introduced in the frame
numbers.
You can use list images to
inspect the frame numbers, image names and number of reflections per image.
Notes:
- If you change the value of fillgap after reading
a data file, you have to read the data again.
- If data originate from still images,
fillgap will be set off.
flagall
Show a summary of the reflection flags for all
reflections (ignoring the reflection
filters)
Flags for all reflections
Good 1365 Not Good 364 (Weak 282 Not Weak 82) Total: 1729
Total Percent NonWeak Percent
1365 78.947 1365 94.333 GOOD
282 16.310 WEAK
13 0.752 13 0.898 EDGEVER
11 0.636 10 0.691 EDGEHOR
28 1.619 25 1.728 EDGEROT
16 1.106 BADBG
11 0.636 11 0.760 NEGATIVE
13 0.752 BADUNIF
1 0.058 1 0.069 MAXSHIFT
More output will be given after output full:
Flags for all reflections
Good 1365 Not Good 364 (Weak 282 Not Weak 82) Total: 1729
Nrefl Code Flags
1365 0 GOOD (hkl 1 1 5 Set s01f001 Refl 5)
278 1 WEAK (hkl -2 -7 0 Set s01f001 Refl 6)
8 2 EDGEVER (hkl -1 -7 -1 Set s01f001 Refl 107)
8 4 EDGEHOR (hkl 1 -2 8 Set s01f001 Refl 105)
1 5 WEAK EDGEHOR (hkl -2 -3 0 Set s02f001 Refl 50)
25 8 EDGEROT (hkl -4 0 -5 Set s01f001 Refl 2)
3 9 WEAK EDGEROT (hkl -4 -5 -4 Set s01f001 Refl 1)
15 1024 BADBG (hkl -2 -7 -3 Set s01f001 Refl 96)
1 1028 EDGEHOR BADBG (hkl -3 2 3 Set s01f001 Refl 447)
11 8192 NEGATIVE (hkl 0 -6 -5 Set s01f001 Refl 271)
7 16384 BADUNIF (hkl -1 -6 -7 Set s01f001 Refl 243)
5 16386 EDGEVER BADUNIF (hkl 0 3 -12 Set s01f001 Refl 699)
1 16388 EDGEHOR BADUNIF (hkl -3 -3 -5 Set s02f001 Refl 173)
1 65536 MAXSHIFT (hkl 1 -6 -4 Set s01f001 Refl 373)
Total Percent NonWeak Percent
1365 78.947 1365 94.333 GOOD
282 16.310 WEAK
13 0.752 13 0.898 EDGEVER
11 0.636 10 0.691 EDGEHOR
28 1.619 25 1.728 EDGEROT
16 1.106 BADBG
11 0.636 11 0.760 NEGATIVE
13 0.752 BADUNIF
1 0.058 1 0.069 MAXSHIFT
flagfilter
Show a summary of the reflection flags for all
reflections using the reflection filters.
Flags for selected reflections
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
Good 1365 Not Good 289 (Weak 278 Not Weak 11) Total: 1654
Total Percent NonWeak Percent
1365 82.527 1365 99.201 GOOD
278 16.808 WEAK
11 0.665 11 0.799 NEGATIVE
More output will be given after output full:
Flags for selected reflections
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
Good 1365 Not Good 289 (Weak 278 Not Weak 11) Total: 1654
Nrefl Code Flags
1365 0 GOOD (hkl 1 1 5 Set s01f001 Refl 5)
278 1 WEAK (hkl -2 -7 0 Set s01f001 Refl 6)
11 8192 NEGATIVE (hkl 0 -6 -5 Set s01f001 Refl 271)
Total Percent NonWeak Percent
1365 82.527 1365 99.201 GOOD
278 16.808 WEAK
11 0.665 11 0.799 NEGATIVE
Related command:
limitcount.
focus
Change focus parameters. The focus parameters
focusdist,
focushor and
focusver are used to calculate
the coordinates (in the lab system) of the primary beam.
The primary beam is used to calculate
capillaryin,
capillaryout,
direction cosines,
eta,
and renninger.
forbid
Set reflection filters. See the flags page.
forbidsets
Syntax: forbidsets variable f1 f2 include
Disable complete sets if at least one of the reflections
(passing the reflection filters)
fulfills the expression variable f1 f2 include.
The expression is similar to limit.
Example:
Disable all sets where at least one reflection has a bgc larger than 25:
set all
minmax bgc
forbidsets bgc 25 highbgc yes
set all
minmax bgc
forbidsets lowbgc 25 no
- In an indirect way: first allow only reflections with bgc
larger than 25. The variable setsize will be nonzero if at least
one reflection passes this filter. Then disable all sets with a
nonzero setsize. And finally remove the limit on bgc.
set all
minmax bgc
limit bgc 25 highbgc yes
forbidsets setsize 0 0 no
unlimit bgc
Related commands:
limit,
set.
Syntax: formulasigext n
Default: formulasigext 1
Sets the formula used for calculation of the external sigma.
1 |
shelx |
2 |
crystals, sortav |
More information in the section
Mean and Sigma
on the main page.
Related command: resultsig.
Syntax: formulasigint n
Default: formulasigint 2
Sets the formula used for calculation of the internal sigma.
1 |
standard |
2 |
shelx |
3 |
crystals |
4 |
sortav |
More information in the section
Mean and Sigma
on the main page.
Related command: resultsig.
frame
Syntax: frame n
Lists reflections within frame n.
This could also be done using:
- limit frame
n n+1 yes
- nofull
- listrefl refl
frameoffset
Syntax: frameoffset n
Default: frameoffset 0
The command sadabs writes reflection information
to a file. The framenumber is part of this information. The value
of n will be added to this framenumber.
full
Use all reflections in reflection listings.
Use nofull to show only selected reflections.
fulldoc
Create a document file. The file will have filetype .doc.
In this file you can find the results of
flagall,
tables for GOOD,WEAK and OVERLAPSPLIT,
tables for GOOD and WEAK and tables for GOOD.
The tables are requiv,
rshell and
rexperiment
Related command: redirect.
gaussexp
Syntax: gaussexp name f1 f2 f3 f4
Sets four coefficients for the gauss function for a specific experiment, identified
with the experimentname (not the experiment number).
The four coefficients are available as reflection variables
gauss1,
gauss2,
gauss3,
gauss4.
Related command:
gausspar
gausspar
Syntax: gausspar f1 f2 f3 f4
Sets four coefficients for the gaussion function:
f4+f1*exp(-0.5*((x-f2)/f3)**2)
These four coefficients are
automatically calculated for each plot or
histo if histoline is on.
The four coefficients can be used by the function gauss() in
expressions.
The four coefficients are available as reflection variables
gauss1,
gauss2,
gauss3,
gauss4.
Related command:
gaussexp
gif
Alias for ps.
go
Use the mouse ( left or
middle) to pick a reflection in one of the
graphics subwindows. The selected reflection will be redrawn in
colour label, and
the equivalent reflections in
colour fg.
The reflections will also be printed.
Click-able windows are created with:
icr,
plot,
plot2d,
plotgrid,
plotray,
plotrefl.
Go back to the commandlevel by clicking the right mouse
button.
good
Allow good reflections in the filters. Use nogood
to disable these reflections.
See the flags page.
gridsize
Syntax: gridsize n
Default: gridsize 5
With the plotgrid command, the average value
of a reflection variable
can be plotted. The area over which the average is calculated is set
by gridsize. gridsize 1 corresponds to the area occupied by one
pixel in the 2d window (note: this is NOT the pixelsize of the
original image). By increasing gridsize, the number of areas decreases
while their sizes increase. The number of areas will also change if
the number of pixels in the 2d window changes (for example after
window fraction).
help
Start a browser with information about any.
histcount
Syntax: histcount variable [scale]
Create a histogram of a variable
along the horizontal axis.
- The histograms are build using
int(variable*scale).
- The vertical value is the
number of reflections in each histogram. Use
histo to plot the average or summed value of a variable.
- Normal reflection filters are applied.
- The plot area is surrounded by a box in
colour ticks.
- The values along the axes are drawn in
colour label.
- The legenda is drawn in
colour text.
- The histogram bars (drawn if
histobar=on)
are drawn in colour histobar1
and colour histobar2.
- The width of the histogram bars can be changed with
histobarfraction.
- The observed histogram lines (drawn if
histoline=on
or histogauss is nonzero)
are drawn in
colour histoline1.
- The calculated histogram line (drawn if
histoline=on
or histogauss is nonzero)
is drawn in
colour histoline2.
- The histogram is printed to the screen if
histoprint=on.
- If histofile is set, the values of the histogram
are written to the histofile.
histcountmax
Syntax: histcountmax n
Default: histcountmax -1.
Controls the scaling of the vertical axis in a
histcount plot. If n is
negative, the plot will be scaled between zero and the highest histogram
value. If n is positive, the scale will extend from
0 to n.
Related command:
histoautoscalever
histo
Syntax: histo variable1 [variable2 [multiply] ]
Create a histogram of two
variables.
variable1 along the vertical and
int(variable2*multiply) along the
horizontal axis. multiply defaults to 1.0
- If you omit variable2 it remains the same
(the initial value is shell).
- The horizontal plot boundaries are automatically set to the extreme values
of variable2, but you may use
histoautoscalehor on to
set the plot boundaries to the actual lowest and highest values.
You may set the horizontal plot boundaries also with
plotlimit.
- The horizontal axis will be divided into bins. The maximum number of bins
is limited, and the program will change multiply to
adapt to this maximum.
- This calculated number of bins can be modified with
histobinfactor
- The value of the histogram is the
average value of variable1 for every bin, unless
histonorm off is set. In that case the
histogram contains the summed value.
- If histosmooth is set to a positive number,
the histogram will be smoothed by a moving mean procedure.
- The values will be assigned to the
histo variable of
all reflections (with corresponding bin).
- Use histcount to plot
the number of reflections in each bin.
- The vertical plot boundaries are automatically set to the extreme values or
variable1, but you may use
histoautoscalever off and/or
plotlimit.
- Normal reflection filters are applied.
- The plot area is surrounded by a box in
colour ticks.
- The values along the axes are drawn in
colour label.
- The legenda is drawn in
colour text.
- The histogram bars (drawn if
histobar=on)
are drawn in
colour histobar1
and colour histobar2.
- The width of the histogram bars can be changed with
histobarfraction.
- The observed histogram lines (drawn if
histoline=on
or histogauss is nonzero)
are drawn in
colour histoline1.
- The calculated histogram line (drawn if
histoline=on
or histogauss is nonzero)
is drawn in
colour histoline2.
- The histogram is printed to the screen if
histoprint=on.
- If histofile is set, the values of the
histogram are written to the histofile.
histoautoscalehor
Syntax: histoautoscalehor state
Default: histoautoscalehor off.
If off, the horizontal scaling of
histo and histcount plots is
set to the extreme values of the horizontal variable.
You may use plotlimit to specify less extreme
values.
If on the horizontal plot boundaries are set to
the actual lowest and highest values (of reflections passing the reflection
filters).
histoautoscalever
Syntax: histoautoscalever state
Default: histoautoscalvere on.
If on, the vertical scaling of
histo and histcount plots
is set to the actual highest and lowest values.
If off, the vertical range is set to the extreme
values of the vertical variable. You may set different limits with
plotlimit.
Related command:
histcountmax
histobar
Syntax: histobar state
Default: histobar on.
If set, the histogram bars of the
histcount and histo commands
wil be drawn in colours
colour histobar1
and colour histobar2.
The width of the bars is controlled by
histobarfraction.
histobarfraction
Syntax: histobarfraction f
Default: histobarfraction 0.95.
The widths of the histogram bars is set to f
(0.0 to 1.0).
histobinfactor
Syntax: histobinfactor f
Default: histobinfactor 1.0.
Divide the number of bins by f.
Related commands:
histo and
histcount
histoempty
Syntax: histoempty f
Default: histoempty 0.0.
The histo command will assign the value of
every bin to the
histo variable of
corresponding reflections. If a bin is empty, the
value of f will be used.
histofile
Syntax: histofile filename
If set, all histogram values of
histcount and/or
histo commands will be written
to filename.
The filename can be closed with
nohistofile.
Related command:
histoprint,
histogauss
Syntax: histogauss n
Default: histogauss 0.
Fits a gauss curve through a histogram.
y = baseline + scale*exp(-0.5*(x-mean)/sigma)²)
- 0. No fit
- 1. Fit a gauss curve.
baseline is fixed to zero.
- 2. Fit a gauss curve.
baseline is fixed to to the lowest
value.
- 3. Fit a gauss curve.
baseline is refined.
If set to nonzero, the curve of the
histcount and histo commands
wil be drawn in colour
colour histoline1
(the observed lines between the maxima of the histogram) and
colour histoline2
(the fitted curve).
Related commands:
histoline,
histolinezero.
histoline
Syntax: histoline state
Default: histoline off.
Fit a line through the histogram
y = a + b*x
If set, the line of the
histcount and histo commands
wil be drawn in colour
colour histoline1
(the observed lines between the maxima of the histogram) and
colour histoline2
(the calculated line).
Whether the lines pass through the origin depends on
histolinezero
Related commands:
histogauss,
histolinezero.
histolinezero
Syntax: histolinezero state
Default: histoline off.
- If set, lines created by histoline will
pass the origin (i.e. y = 0 + b*x).
- If unset, general lines will be created by
histoline (i.e. y = a + b*x).
Related commands:
histogauss,
histoline.
histonorm
Syntax: histonorm state
Default: histonorm on
The hist command calculates histograms
by summing the values of the requested
variable per bin. With histonorm on, these sums are
divided by the number of contributing reflections and the histogram
will thus show the average value per bin.
With histonorm off, these sums are not divided,
and the histogram will show the summed value per bin.
histoprint
Syntax: histoprint state
Default: histoprint off.
If set, histograms, created by then
histcount and histo commands,
will be printed to the screen.
Related command:
histofile,
histosmooth
Syntax: histosmooth n
Default: histosmooth -1.
If set to a positive value, the histogram, created by
histo, will be smoothed by a moving mean
with size 2n+1 (from -n
to +n)
hkl
Syntax: hkl h k l
You enter the reflection indices h k
and l. The program shows
all reflections equivalent to this one.
Any> hkl 1 2 3
hkl 1 2 3 ==> hkl -1 -2 -3 multiplicity=8
hkl= -1 -2 -3 hb=407 im=1 Theta=10.038 Reso=2.039 MaxInt=197.755 Mean=186.509
(mean=186.509 sigma=1.415 ioversig=131.853 <ioversig>=44.237) Rsym=0.023 Rmeas=0.024 Rpim=0.008 Chi2=1.84 n=10
nr boxfil set refl h k l intensity sigma ioversig delsig
4 s01f001: 1 4 1 -2 3 171.75 3.59 47.86 -4.11 -EDGEROT
42 s01f001: 1 42 -1 -2 -3 184.98 4.35 42.56 -0.35 +
232 s01f001: 1 232 -1 2 3 187.02 4.26 43.90 0.12 +
417 s01f001: 1 417 1 -2 -3 182.90 4.20 43.56 -0.86 +
564 s01f001: 1 564 -1 2 -3 190.27 3.96 48.01 0.95 +
665 s01f001: 1 665 -1 -2 3 176.32 4.04 43.60 -2.52 +
843 s01f001: 1 843 1 2 -3 184.17 4.23 43.57 -0.55 +
1244 s01f002: 2 241 -1 -2 -3 182.78 4.27 42.83 -0.87 +
1333 s02f001: 3 25 1 2 3 189.10 4.23 44.66 0.61 +
1510 s02f001: 3 202 -1 2 -3 197.75 4.41 44.87 2.55 +
1549 s02f001: 3 241 -1 -2 -3 189.44 4.23 44.82 0.69 +
10 average: 186.47 4.22 44.24 -0.02
Explanation:
- First line
- The indices are converted to a standard setting according to
point group (see pg)
- Second line
- hkl = standard setting
- theta = resolution in degrees
- Reso = resolution in Å
- Maxint = maximum intensity of the list of equivalent reflections
- Mean = mean intensity calculated according to the selected
weighting scheme (see Mean and Sigma
on the main page).
This number is evaluated after reading the data, and at the
invocation of the delsig command.
- Third line (only displayed if at least two equivalent reflections contribute to the mean)
- mean = mean intensity calculated according to the selected
weighting scheme (see Mean and Sigma
on the main page)
and current filter settings.
- sigma = sigma of the mean intensity calculated according to the selected
weighting scheme (see Mean and Sigma
on the main page)
and current filter settings.
- ioversig = mean/sigma
- <ioversig> = mean of contributing individual intensity/sigma
- various Rsym values.
- n is the number of reflections used in the calculations of Rsym.
- Remaining lines: a table with the equivalent reflections.
- columns in the table can be added or removed with
field and fieldflags.
- the number of decimals for floating point numbers can be changed
with precision.
In all (present) graphics subwindows
the reflections will be highlighted with
colour fg and
size specified by dotsize.
For incommensurate datasets you may use the following commands:
- hkl
This will list all reflections h k l for all
available qvector combinations.
- qvp and hkl
First you set a filter on qvector combinations. The
hkl command (any many others) will
be limited to those reflections only.
- hklm
You specify h k l and the indices for
the available qvectors (m n p).
Related command:
limithkl.
hklclose
Closes the reflection file. (created with
hklf4,
hklf4merge,
hklf5,
hklf5merge,
hklf6merge,
sadabs or
twinabs).
If hkltrailer=on, an
XML-like trailer is appended to the reflection file if the file
was created with
hklf4,
hklf4merge,
hklf5 or
hklf5merge.
hklf4
Writes hklf4 formatted reflections to the reflection file.
- If the file is not yet opened, it will be created
(the filename may be set using hklfilename).
- The batchnumber in the reflection file is controlled by
hklf4batch.
- Direction cosinesare calculated if cosines
is on.
- Cell parameters are written to a file filename.par.
- Subsequent hklf4 commands will append reflections
to the reflection file.
- Use hklclose to close the reflection file.
- The intensities of the reflections will be scaled
(if necessary) to accomodate the shelx F8.0 format. You may
use the hklf4sf command to fix the scale factor.
See also the Reflection output section on the
main page.
Related command: hklf4merge.
hklf4batch
Syntax: hklf4batch variable
Default: hklf4batch expnr
Sets the reflection file batch number. You may choose any available
variable, but the resulting batchnumber is limited to values between
1 and 9999.
Examples:
- hklf4batch expnr
batch number is experiment number (the default)
- hklf4batch set
batch number is boxfile number.
- expression 7 1
hklf4batch exp7
A fixed value of 1 for all reflections.
hklf4merge
Writes hklf4 formatted reflections to the reflection file.
Equivalent reflections are merged.
See the section
Mean and Sigma
on the main page for more
information on the calculations of average intensity and
corresponding sigma.
- If the file is not yet opened, it will be created
(the filename may be set using hklfilename).
- The batchnumber in the reflection file is controlled by
hklf4batch.
- Direction cosines are calculated if cosines
is on.
- Cell parameters are written to the file filename.par.
- Subsequent hklf4merge commands will
append reflections to the reflection file.
- Use hklclose to close the reflection file.
- The intensities of the reflections will be scaled
(if necessary) to accomodate the shelx F8.0 format. You may
use the hklf4sf command to fix the scale factor.
See also the Reflection output section on the
main page.
Related command: hklf4.
hklf4sf
Syntax: hklf4sf f
Default: hklf4sf -1
Sets the intensity scale factor to acccomodate the shelx format.
If f is negative (the default), a scalefactor
for writing the reflection file will be calculated to accomodate the
shelx format. For every hklf4,
jana or sadabs
command, a
new scale factor will be used. If this is not desirable you may
set the scalefactor with hklf4sf.
If you select a too large value for f you may
run into formaterrors while writing the reflection file.
hklf5
Writes hklf5 formatted reflections to the reflection file.
- If the file is not yet opened, it will be created
(the filename may be set using hklfilename).
- Reflections passing the reflection filters
are written with the corresponding batch number
(see hklf5batch).
- If a (selected) reflection has the
OVERLAPSUM flag set, the
overlapping reflections will be written to the reflection file (with their
corresponding batchnumber (hklf5batch).
(See omitsum to control this behaviour).
- Direction cosines are calculated if cosines
is on.
- Cell parameters are written to a file filename.par.
- Subsequent hklf5 commands will append reflections
to the reflection file.
- Use hklclose to close the reflection file.
- The intensities of the reflections will be scaled
(if necessary) to accomodate the shelx F8.0 format. You may
use the hklf5sf command to fix the scale factor.
See also hklf5merge.
See the TwinPage
for an example of handling twin datasets.
See also the Reflection output section on the
main page.
hklf5batch
Syntax: hklf5batch n n ...
Sets the batch number for the hklf5,
hklf5merge and
twinabs output
files. Specify a batch number for each of the available matrices.
hklf5merge
Writes hklf5 formatted reflections to the reflection file.
Equivalent reflections are merged.
See the section
Mean and Sigma
on the main page
for more information on the calculations of average intensity and
corresponding sigma.
- If the file is not yet opened, it will be created
(the filename may be set using hklfilename).
- Reflections passing the reflection filters
are written with their corresponding batch number
(see hklf5batch).
- If a (selected) reflection has the
OVERLAPSUM flag set the
overlapping reflections will also be written to the output file
(with corresponding batchnumber (see hklf5batch).
(See omitsum to control this behaviour).
- Direction cosines are not written.
- Cell parameters are written to the file filename.par.
- Subsequent hklf5merge commands will append
reflections to the reflection file.
- Use hklclose to close the reflection file.
- The intensities of the reflections will be scaled
(if necessary) to accomodate the shelx F8.0 format. You may
use the hklf5sf command to fix the scale factor.
See also hklf5.
See also the Reflection output section on the
main page.
hklf5sf
Syntax: hklf5sf f
Sets the intensity scale factor to acccomodate the shelx F8.0 format.
The default value of f is set during the first
evocation of one one the hklf5,
hklf5merge or twinabs
commands.
Use this command to fix the scale factor.
Related commands:
hklf5,
hklf5merge and
twinabs.
hklf6merge
Writes hklf6 formatted reflections to the reflection file.
Equivalent reflections are merged.
See the section
Mean and Sigma
on the main page for more
information on the calculations of average intensity and
corresponding sigma.
- If the file is not yet opened, it will be created
(the filename may be set using hklfilename).
- The batchnumber in the reflection file is controlled by
hklf4batch.
- There are no direction cosines.
- Cell parameters are written to the file filename.par.
- Subsequent hklf6merge commands will
append reflections to the reflection file.
- Use hklclose to close the reflection file.
- The intensities of the reflections will be scaled
(if necessary) to accomodate the shelx F8.0 format. You may
use the hklf4sf command to fix the scale factor.
See also the Reflection output section on the
main page.
hklfilename
Syntax: hklfilename filename
Default: hklfilename shelx
Sets the filename for various reflection output files. The filetype depends
on the command with which the reflections are written.
The commands
hklf4,
hklf4merge,
hklf5,
hklf5merge and
hklf6merge
do not close the reflection file. Subsequent invocations will
append reflections to the file unless the file is closed with
hklclose.
hklm
Syntax: hklm h k l m n p
You enter the reflection indices h k l
and (depending on the number of q-vectors) m n
and p.
The output is similar to hkl
hkltrailer
Syntax: hkltrailer state
Default: hkltrailer on
If on, a XML-like trailer will be written
at the end of a reflections file created by
hklf4 or hklf5.
htmlfilter
Syntax: htmlfilter state
Default: htmlfilter off
If on, (almost) all screen output will be
filtered. The following changes are made:
- & → &
- > → >
- < → <
This allows inclusion of the output into html documents.
hybrid
Syntax: hybrid f15 f14
Default: hybrid 5 10
Combines eval14 and eval15 data. The value of variable
ioversig and
the arguments
f15 and f14 control the
merging.
- ioversig ≤f15
Use eval15 data. (f = 1.0)
- f15<ioversig <f14
Use a linear interpolation of eval15 and eval14 data
( f =
(f14-ioversig)/(f14-f15) )
- ioversig ≥f14
Use eval14 data. (f = 0.0)
The new intensity and sigma are calculated as follows:
The new intensity is f *
I15 +
(1 - f) * I14
The new sigma is √( f² *
σ15² +
(1 - f)² *
σ14² )
Note:
hybridscale
Syntax: hybridscale f
Default: hybridscale 1.0
The hybrid commands combines eval15 and eval14 results.
The values of
eval14i and
eval14s
will be multiplied by f
If f is negative, the scale factor will be calculated
from Σ(Intensity)/Σ(eval14i).
icr
Syntax: icr variable
Creates a plot with the deviations of reflections with respect to their
mean value.
Only reflections with at least icrmin
equivalents are plotted. The points are plotted as a function of
variable.
Icr Examples
- icr frame
Plot as a function the frame number. This is useful to inspect scaling, decay
or absorption.
- icr hor
Plot as a function the horizontal impact position. This is usefull to
inspect systematic errors.
- icr incid
Plot as a function the incidence angle.
The size of the dots depends on dotsizetype
and the colour depends on colourtype.
If the
deviations are too big (see icrscale),
the point will be plotted at the edge of the area
(these extreme reflections are also printed to the screen if
output=full).
Use icrtype ioversig if you want to examine
the differences of Intensity/Sigma instead of Intensity.
Use icrscale to change the vertical plotrange.
The coordinates of the icr plot can be written to a file with
the xyicr command.
Once the icr subwindow is active, reflections
selected with hkl will be highlighted.
icrmin
Syntax: icrmin n
Default: icrmin 2
Sets the minimum number of equivalents for the
icr command.
icrscale
Syntax: icrscale f
Default: icrscale 2.0
Sets the vertical scale for the
icr command.
The default value of 2 will set the display reflections from 1/2 to 2
times the average intensity.
icrtype
Syntax: icrtype
intensity/ioversig
Default: icrtype intensity
Applicable in the various icr plots.
On intensity, the deviations of all reflections
from their mean intensity is plotted.
On ioversig, the value of Intensity/Sigma
is plotted relative to the mean of Intensity/Sigma.
imagescale
Syntax: imagescale all/nr f
Change the scale factor of one (nr) or
all images to f.
Related commands:
imagescaleread,
imagescalefine,
write imagescale.
imagescalecalc
Syntax: imagescalecalc n
Calculate image scale factors.
- For each frame, two sums are calculated:
s1 = Σ Ii
s2 = Σ Imean
where Imean is the mean intensity
of the set of equivalent
reflections corresponding to Ii
(as calculated by delsig).
- For each frame, two sums are calculated
t1 = Σ s1
t2 = Σ s2
where Σ runs over +/-n neigbouring frames.
(the neighbouring frames will not cross experiment boundaries)
- The image scale factor for each frame will be set to
t2/t1
(Ii/Imean corresponds to the icr-value of a reflection).
The imagescalefactors are accessible via reflection variable
imscale
imagescaleread
Syntax: imagescaleread file
Read image scale factors from file.
Related commands:
imagescale,
imagescalerefine,
write imagescale.
imagescalerefine
Syntax: imagescalerefine f first last
Refine image scale factors starting with an initial offset of
f. Only image between first and
last are considered.
- start an image loop from first to last.
- collect all (mother) reflections with at least one daughter in this image.
- Obtain the image scale factor for this image which minimizes Rsym for the selected (mother) reflections.
Related commands:
imagescale,
imagescaleread,
write imagescale.
imagescaletop
Syntax: imagescaletop n
Change the image scale factors.
- For each image, the intensities of
reflections are collected.
- If more than n reflections are found,
the imagescale for that image
is divided by the sum of the n strongest intensities.
- If less than n contributing relfections
are found the imagescale for that image wqill be set to zero.
Related commands:
imagescaleread,
imagescalefine,
write imagescale.
jana
Syntax: jana
Creates a file jana.hkl.
This reflection file can be used by the jana program
(by Vaclav Petricek and Michal Dusek,
www-xray.fzu.cz/jana/jana.html).
jpeg
Alias for ps.
kfactor
Syntax: kfactor kf
Default: kfactor 1.0
Recalculate the sigma's of all reflections.
See the Mean and Sigma section
on the main page.
kf is available in
expressions as kfactor.
Related commands:
automsa,
msa,
msapower,
sigmafactor.
label
Syntax: label n
Default: label 5
Sets the number of requested horizontal labels in the
plot, histo and
histcount plots to n.
A value of zero disables both horizontal and vertical labels.
laue
Change the point group symmetry. The point group symmetry is extracted
from the input file, but you may change it.
Syntax: laue code
laue group overview
Related commands:
bravais,
pg.
limit
Syntax: limit variable f1 f2 include
Sets limits on one of the reflection
variables.
If include=yes, the value of
variable should lie between
f1 and
f2.
If include=no, the value of
variable should lie outside
f1 and
f2.
Unset the limits by invoking limit with (too) extreme
values, or use the unlimit command.
Note: you can use limit set and
limit expnr, but
the commands set and
experiment give
more flexibility.
Related commands:
forbidsets,
limitcount,
limithkl,
plotlimit,
varlimit.
limitcount
Show the reflection count for all limit conditions.
This command will set the constant nallowed, which
can be used with expression.
Related command:
flagfilter.
limithkl
Syntax: limithkl h k l
Limits reflection indices to h k l
and their equivalents.
From the reflection indices h k l
calculate the pointers
hb and
im of the mother reflection
and set a limit on these pointers.
Related commands:
hkl,
limitrecip,
unlimithkl.
limitrecip
Syntax: limitrecip h k l
Each of the arguments can be
- a number
limits the reflection index to this number
- a single character (corresponding to the index name)
no limit on this index.
- a plus or minus sign with a single character (corresponding to the index name)
a positive or negative limit on this index.
Examples:
- limitrecip 0 k 0
This limits to the b* reflections
- limitrecip 0 +k 0
This limits to the positive b* reflections
- limitrecip h k 0
This limits to reflections in the a*-b* plane.
- limitrecip h k 2
This limits to reflection with l=2
- limitrecip 4 1 3
This limits to reflection with with indices 4 1 3. You may
want to use limithkl to set limits
on reflection indices and their equivalents.
Related command:
unlimitrecip.
limitrmat
Syntax: limitrmat n
allow/forbid/list
If more than one matrix is set (by the input file), you may use this
command to control the influence of interfering matrices.
n is the matrix number of the interfering matrix
(you may specify limitrmat all).
- allow
Reflections of the main lattice overlapping with reflections of
lattice n are allowed.
- forbid
Reflections of the main lattice overlapping with reflections of
lattice n will be forbidden
- list
List all reflections of the main lattice overlapping with reflections of
lattice n. This
command will give only output if
allow OVERLAPSUM
is set.
Related command:
listrefl overlap.
The limitrmat command is only available of multiple cell matrices
are available.
list
Syntax: list exp/expactive/expression/images/intensity/possible/region/scans/set/var
- exp
Prints information of all experiments.
For each experiment: distance, swing, scalefactor, goniostat, rotation axis
and sets. (You may want to use euler on to convert
kappa angles to eulerian angles). See also list scans.
- expactive
Prints information of active experiments (i.e. experiments selected with the
experiment filter).
For each experiment: distance, swing, scalefactor, goniostat, rotation axis
and sets. (You may want to use euler on to convert
kappa angles to eulerian angles). See also list scans.
- expression
Prints the definitions of expressions (as defined by
expression).
- images
Print frame numbers, image scale factor, number of reflections per frame and image names.
You may specify a lower and upper limit. For example
list images 100 110
A reflection count of zero indicates a gap in the image name sequence.
You may want to set fillgap to
off
before reading an input file to avoid these gaps).
- intensity
Print intensity statistics by Resolution Shells for both unmerged and
merged data.
Note, the average values in the I s I/s columns are simple
averages (sum divided by number). In the right part of the table
(the merged values), the intensities and sigma's of equivalent
reflections are first merged (using weightint
and weightsig).
The resulting I s values are subsequently averaged.
- possible
Prints all possible
variables.
Some variables do not vary. For example the variable
dist
acts as a constant if only one experiment is available. Use
list variables to limit the output to
varying variables.
- region
Prints all regions. See region.
- rmatexp
Prints cell parameters per experiment (as defined by the y-file, or
rmatexp). If
output=full a second table
will be printed with the rmat coefficients per experiment.
Related commands:
rmatexp,
useintegratedcell,
write rmatexp.
- scans
Prints information of experiments.
For each experiment the values of distance, theta, goniostat and
rotation axis, as read from the image headers.
- set
Prints information of sets.
set number and name, experiment number
- var [variable]
Prints all available
variables.
Some variables do not vary. For example the variable
dist
acts as a constant if only one experiment is available. Use
list possible to list non-varying
variables also.
Sample output (output
normal):
name low high limit1 limit2 ins plot1 plot2 var1 var2
chi 0.00 73.20
delsig -42.65 5.92
.........
Legenda of list var
column |
what |
description |
name |
variable |
name of the variable |
low |
lowest value |
|
high |
highest value |
|
limit1 |
Lower Limit |
if limit is set |
limit2 |
Upper Limit |
if limit is set |
ins |
inside (yes or no) |
if limit is set
|
plot1 |
Lower Limit |
if plotlimit is set |
plot2 |
Upper Limit |
if plotlimit is set |
var1 |
Lower Limit |
if varlimit is set |
var2 |
Upper Limit |
if varlimit is set |
Sample output (output
full):
name pos avl ave T P FS scale ES Lim low high limit1 limit2 ins plot1 plot2 var1 var2
bga no no yes 1 22 2 0.00 1
bgb no no yes 1 23 2 0.00 1
bgc no no yes 1 24 2 0.00 1
chi yes yes yes 3 4 2 57.30 4 no 0.00 73.20 0.00 73.20 yes 0.00 73.20 0.00 73.20
delsig yes yes yes 1 16 2 1.00 4 no -42.65 5.92 -42.65 5.92 yes -42.65 5.92 -42.65 5.92
........
Legenda of list var
column |
what |
description |
related command |
name |
variable |
name of the variable |
|
pos |
possible |
yes or no |
|
avl |
available
| yes = possible, and variation in value
no = not possible or possible but no variation in value |
|
ave |
average
| yes = print average value per mother
no = do not print average value per mother |
hkl,
listrefl bad/mother/reject
|
T |
Type |
1=float 2=integer 3=on-the-fly float 4=on-the-fly integer |
|
P |
Pointer |
internal pointer |
|
DP |
decimal places |
used for display |
precision |
scale |
scale factor |
used for display (57.30=radians→degrees) |
presentationscale |
ES |
Extreme Status |
1=ignore 2=postpone 3=ToDo 4=Done |
|
Lim |
Limit Set |
yes/no |
limit |
low |
lowest value |
|
|
high |
highest value |
|
|
limit1 |
Lower Limit |
if limit is on |
limit |
limit2 |
Upper Limit |
if limit is on |
limit |
ins |
inside |
yes=inside limits
no = outside limits
| limit
|
plot1 |
Lower Limit |
for graphics |
plotlimit |
plot2 |
Upper Limit |
for graphics |
plotlimit |
var1 |
Lower Limit |
colour scheme for colourtype var |
varlimit |
var2 |
Upper Limit |
colour scheme for colourtype var |
varlimit |
listrefl
Syntax: listrefl
refl/mother/equiv/bad/reject/renninger/lost/mis/overlap/overequiv/double/2over
Various commands to show reflections.
See the section
Mean and Sigma
on the main page
for more information on calculating average intensities and
corresponding sigma's.
For most of the options, a table of reflections is printed.
With field variables can be added or removed in this
table. The number of decimals for floating
point variables can be changed with precision.
- refl
Shows reflections.
You may specify a lower and upper limit (reflection number). For example
listrefl refl 100 110
Related command:
dumpid
- mother
Shows reflection in the order of the mother list.
(For each mother, the corresponding daughters will be printed)
- equiv
Shows reflections with same number of equivalents.
You enter the number of equivalents (or q).
If the number of daughters passing the reflection filter
(allow,
forbid,
require,
good and
nogood)
equals the specified number, the reflection is printed.
- full = on.
Show all reflections.
- full = off.
Show one line for each reflection
- bad
Show bad reflections. Bad is based on the χ² value. Only the
worst 30 reflections are shown.
Add or remove variables in the table with field.
- reject
Show rejected reflections.
If at least one of the daughters is rejected all daughters are printed.
- renninger
Prints all renninger reflections.
The table contains the fields:
nr, h, k, l, intensity, nhit, nmis, renninger, and the indices of
the 2 corresponding reflections.
Related commands:
primarybeam,
renninger and
testrenninger.
- lost
Show reflections available in the dataset, but none of the
equivalents passes the reflection filter
(allow,
forbid,
require,
good and
nogood)
- mis
show reflections not available at all in the dataset.
You may use thminzero to set the minimum theta
value to zero or to the lowest encountered theta.
The table prints h, k, l, theta, resolution
- overlap
Shows all interfering reflections.
No reflection filters are used.
The table contains
h, k, l of the interfering reflections with pointers to
sisters (more interfering) and ids (the reflection in the main lattice)
with the matrixnumber and indices of the main lattice reflection.
Related command:
limitrmat 2 list
- overequiv
Shows the number of all overlapping reflections where the main
reflection is equivalent to the overlapping one.
If output=full, these
reflections are printed to the screen.
This command will only run if
allow OVERLAPSUM is set.
Related command: deoverlap.
- double
Show reflections which have been measured more than once with
different overlap status. It is possible that a reflection has the
OVERLAPSUM flag set, but at some other
point in the dataset the same reflection is not overlapping.
This command will only give output if the
allow OVERLAPSUM is set.
On writing
a hklf5
hklf5merge and
twinabs output
file, the command omitsum
controls the writing of double reflections.
- 2over
Show reflections which have been measured more than once with
different overlapping neighbours.
This command will only give output if
allow OVERLAPSUM is set.
meantype
Syntax: meantype runtime/imean/icalc
Default: meantype runtime
In the various calculations of Rsym and χ² (see
Rsym), the average intensity
Imean and the number of
contributing equivalent reflections nmean
are used.
- runtime
The average intensity is calculated
on the fly, using only the equivalent reflections passing the
current reflection filter.
- imean
Use the
(mother database)
values of
mimean and
minmean.
- icalc
Use the value of icalc. See readicalc.
minmax
Syntax: minmax variable
Prints the actual minimum and maximum value of variable.
These values are also available to expressions
as lowvariable and
highvariable.
minmax example
Select the reflections with a maximum number of equivalents. Firstly determine
the maximum (and minimum), and secondly set a limit.
minmax equiv2
limit equiv2 highequiv2 highequiv2 yes
modprint
Syntax: modprint n
Default: modprint 10000
If the number of reflections exceeds 10*n, a progress
meter will be displayed for various file-generating commands
(
hklf4,
hklf4merge,
hklf5,
hklf5merge,
hklf6merge,
mtz,
pk,
sadabs,
twinabs
). Every n'th reflection will display a
character on the screen.
more
Syntax: more state
Default: more on
Some commands display only one page of output. The user has to
enter a enter to get the next screen. By setting more to
off
all output is given without interrupt ( and you can limit the amount of
output with nnomore).
msa
Syntax: msa factor
Default: msa 0.02
Recalculate the sigma's of all reflections.
See the Mean and Sigma section
on the main page.
sadabs and twinabs
determine msa themselve. Both programs assume the
sigma's are as they are, thus with no msa correction. There is a but:
reflections with an I/σ exceeding 1000.0 are (erroneously)
flagged as too strong by Sadabs. To avoid unnecessary elimination of these
reflections, an msa of 0.001 is advised. The sadabs
command will increase msa to 0.001 if necessary.
factor is available in
expressions as msafactor.
Related commands:
automsa,
kfactor,
msapower,
sigmafactor.
msapower
Syntax: msapower f
Default: msapower 2
Recalculate the sigma's of all reflections.
See the Mean and Sigma section
on the main page.
f is available in
expressions as msapower.
Related commands:
automsa,
kfactor,
msa,
sigmafactor.
mtz
Syntax: mtz boxprofile/merge/separate/unmerged
Creates a ccp4 mtz file
(the filename may be set using hklfilename).
- If spgr has not been set, the spacegroup
number has to be given.
- A reflection file
hklfilename.hklin is created.
- boxprofile
The reflections will not be merged. The reflection file
will contain h, k, l, batchnr (=framenumber),
Iprof, σprof
(the eval15 integration results),
Ibox, σbox
(the result of the eval14 like box integration)
fractioncalc (=1.0)
and the impact positions hor, ver and rot
H K L BATCH I SIGI IPR SIGIPR FRACTIONCALC XDET YDET ROT
- merge
Reflections will be merged. The output consists of h, k, l, I, σ
H K L I SIGI
- separate
Reflections will be merged, but Friedels are kept separate.
Output consists of h, k, l, I, σ, I+,
σ+, I-, σ-
H K L I SIGI I(+) SIGI(+) I(-) SIGI(-)
- unmerged
Reflections will be not be merged. The output consists of
h, k, l, m/isym (=0), batchnr (=framenumber), I, σ
fractioncalc (=1.0), hor, ver and rot
H K L M/ISYM BATCH I SIGI FRACTIONCALC XDET YDET ROT
- a script tomtz is created. This script will
invoke a sequence of CCP4 programs:
- f2mtz
Set title, symmetry, cell parameters. Read reflections from a formatted
file and convert to MTZ format
- sftools
Set wavelength
- sortmtz
Sort reflections on h, k and l
- ctruncate
Convert intensities to structure factors.
Sometimes ctruncate aborts.
If mtzctruncate is on
the script will fall back to truncate.
If mtzctruncate is off
the script will only runtruncate.
- unique
set resolution limit and expand dataset
- freerflag
sets free Rflag for 5% of the reflections
- cad
combine unique and freerflag
- freerflag
complete free
- If the CCP4 programs are available, the script tomtz
will be executed.
Note: mtz files are not suitable for incommensurate data, nor for
overlapping reflections.
mtzctruncate
Syntax: mtzctruncate state
Default: mtzctruncate on
Controls whether mtz will run
ctruncate and/or truncate.
- on.
run ctruncate and on failure fall back to
truncate.
- off.
run truncate.
ncolourvar
Syntax: ncolourvar n
Default: ncolourvar 11
Sets the number of colours used in the
colourtype var
command. The maximum value is 11.
Related commands:
colourscheme,
colourtype var,
colourvar,
varlimit.
nintensity
Syntax: nintensity n
Sets the number of intensity intervals used by the command
rint.
nnomore
Syntax: nnomore n
Default: nnomore 500
Sets the maximum number of output lines for list commands if
more is off.
Syntax: noextract
Closes the file created with extract.
nofull
Show only selected reflections in the listings. Use full
to show all reflections.
nogood
Disables good reflections in the filters. Use good
to enable these reflections.
See the flags page.
nohistofile
Syntax: nohistofile
Closes the file, created by histofile.
noplotlimit
Syntax: noplotlimit variable(s)
Reset the limits (set by a previous plotlimit
command) to the initial values.
nopresentationscale
Disable the presentationscale.
novarlimit
Syntax: novarlimit variable(s)
Reset the limits (set by a previous varlimit
command) to the initial values.
nshell
Syntax: nshell n
Default: nshell 10
Sets the number of shells for the rshell command.
The shell boundaries are calculated to let each shell have approximately the
same number of reflections. See autoshell.
omitsum
Syntax: omitsum state
Default (eval14 data): omitsum on
Default (eval15 data): omitsum off
This command affects the writing of OVERLAPSUM reflections using the
hklf5 ,
hklf5merge and
twinabs commands.
- omitsum on
If a specific reflection has the
OVERLAPSUM flag set,
the list of sisters
of this reflection is scanned to see if a nonoverlapsum sister exists.
Then, if at least one non-overlapping reflection is found,
the overlapping reflection will be ignored.
- omitsum off
All overlapping reflections will be written to the reflection file
Use listrefl double to
show these reflections.
Note: omitsum has no effect if overlapping reflections
are not allowed (forbid overlapsum).
output
Syntax: output none/normal/full/debug
Default: output normal
Controls the amount of output.
- none
- normal
- full
- print extreme reflection indices for
icr plots
- print full output in list var
- print second table in list rmatexp
- print full output in show
- debug
overwrite
Syntax: overwrite variable expressionnr
Overwrites one of the
variables. Use this
with extreme caution. The contents
of variable will be overwritten by the value of
expression expressionnr.
Related command:
copy
pg
Syntax: pg code
Change the point group symmetry. The point group symmetry is extracted from the
input file, but you may change it.
This will change the
mother database
and consequently the variables
delsig,
equiv,
half1,
half2,
imean,
imin,
iplus.
smean,
smin,
splus.
pointgroup overview
Related commands:
bravais,
laue.
pk
Syntax: pk
Creates a pk file. This file
can be used by peakref to do a refinement of
detector positions, cell dimensions and orientation.
The format of the reflections in the
pk file depends on
- For Eval15 data:
- for Eval14 data:
- shiftframes > 1.0:
HKL h k l Intensity hor ver rot
- shiftframes < 1.0:
HKL2 format h k l Intensity hor ver rot rotmin rotmax
In peakref, the rotation residue depends on data format:
- HKL2 format
- rotcalc < rotmin
resrot = abs(rotcalc − rotmin).
This residue is used in rotoutside (and rotall)
- rotmin ≤ rotcalc ≤ rotmax
resrot = abs(rotcalc − rotobs).
This residue is used in rotinside (and rotall)
- rotcalc > rotmax
resrot = abs(rotcalc − rotmax).
This residue is used in rotoutside (and rotall)
- HKL format
- still images:
resrot = abs(theta − eta)
theta is calculated from the cell parameters and reflection indices h,k, and l.
90-eta is the angle between the reflection normal and the primary beam.
resrot is added in rotpartial (and rotall).
Each image in a series of stills will be considered as a different experiment ( implicit
pkbatch set)
- rotation images:
resrot = abs(rotcalc − rotobs)
This residue is added in rotpartial (and rotall)
Only reflections with
nshift>0 are written to
the pkfile.
Nonpositive values of nshift have the following meaning:
- -2:
the reflection was shifted, but something went wrong. For example
the reflection indices of the definitive position differ from
the original indices.
- -1:
The reflection is not shifted because the initial intensity/sigma
was less than isigrefine.
- 0:
An attempt was made to shift the reflection, but the initial
foms were well below the refinement tolerance.
If you want to write non-shifted reflections also, you may want
to use
copy nshift expr7
expression 8 1
copy expr8 nshift
pk
copy expr7 nshift
Related commands:
pkbatch,
pkexpnr,
pkfile,
pkimean,
pkimpact,
pkset.
pkbatch
Syntax: pkbatch expnr/set
Default: pkbatch expnr
Sets the pk file experiment name (similar to
hklf4batch).
- set
Every set will be an experiment in the pkfile. This allows you to
use the peakref command byexperiment per set instead of per
experiment.
- expnr
Every experiment will be an experiment in the pkfile.
pkexpnr
Syntax: pkexpnr
Creates separate pk files for all
allowed experiments. With experiment or
set the allowed sets can be defined.
The pkfiles
can be used by peakref to do a refinement of
detector positions, cell dimensions and orientation. Only reflections
with shifts are written to the
pk file.
The format of the reflections in the
pk file depends on
shiftrestfrac and
pkrestfrac.
The pkfiles will be named
expnr.pk.
Related commands:
pk,
pkset.
pkfile
Syntax: pkfile name
Sets the name of the
pk file
which can be created with the pk,
pkexp, pkimpact and
pkset commands.
pkimean
Syntax: pkimean state
Default: pkimean off
If on, the value of
imean
will also be written by pk to the pkfile. This information can be useful
in peakref to minimize delicr=abs(Ii-Imean)/Imean.
pkimpact
Syntax: pk
Creates a pk file. This file
can be used by peakref to refine for example
detector position, cell dimensions and orientation.
- The reflections are written without reflection indices, so
peakref will interpret them as peaksearch
data.
- The format of the reflections is:
Impact Hor Ver Intensity cx cy cz |c|
- Only reflections with
nshift>0 are written to
the pkfile.
Nonpositive values of nshift have the following meaning:
- -2:
the reflection was shifted, but something went wrong. For example
the reflection indices of the definitive position differ from
the original indices.
- -1:
The reflection is not shifted because the initial intensity/sigma
was less than isigrefine.
- 0:
An attempt was made to shift the reflection, but the initial
foms were well below the refinement tolerance.
- If you want to write non-shifted reflections also, you may want
to use
copy nshift expr7
expression 8 1
copy expr8 nshift
pkimpact
copy expr7 nshift
|
save the values of nshift
all expr8 values are 1
overwrite nshift with 1
create the pk file
restore nshift
|
Related commands:
pk,
pkbatch,
pkexpnr,
pkfile,
pkimean,
pkset.
pkkenau
Syntax: pkkenau
Create a pkfile (similar to pk), but for every
group of equivalent reflections, only the strongest one will be
written.
pkrestfrac
Syntax: pkrestfrac f
Default: pkrestfrac 0.0
Determines the format of reflections while writing a pkfile with
pk, pkexpnr and
pkset.
Reflections with
- shiftrestfrac>f
will be written with format HKL. These reflections span
more than one frame in the simulation in EVAL15, and their observed rotation value
will be used in the calculation of resrot in peakref (rotpartial).
- shiftrestfrac<f
will be written with format HKL2. These reflections lie
within one frame in the simulation in EVAL15, and the frame bounderies will be
used in the calculation of resrot in
peakref (rotinside/rotoutside).
shiftrestfrac
is only available with EVAL15 data.
pkset
Syntax: pkset
Creates separate pk files for all
allowed sets (i.e.boxfiles). With experiment or
set the allowed sets can be defined.
These files can be used by peakref to do a refinement of
detector positions, cell dimensions and orientation. Only reflections
with shifts are written to the
pk file.
The format of the reflections in the
pk file depends on
shiftrestfrac and
pkrestfrac.
The files will be named
set.pk.
You may want to use
pkbatch set
to create a single pkfile, where each set will be treated by peakref as
a separate experiment.
Related commands:
pk,
pkexpnr.
plot
Syntax: plot variable1 [variable2]
Create a scatterplot of two
variables.
variable1 along the vertical and
variable2 along the horizontal axis.
If you omit variable2 it remains the same
(the initial value is frame).
Normal reflection filters are applied.
You can set the plot boundaries with plotlimit.
The plot area is surrounded by a box in
colour ticks,
the values along the axes are drawn in
colour label and
the legenda is drawn in
colour text.
The colour of the dots is set by colourtype,
and the size is set by dotsizetype.
The coordinates of the plot can be written to a file with
the xyplot command.
Once the plotrefl subwindow is active, reflections
selected with hkl will be highlighted.
plot2d
Syntax: plot2d
Plots the horizontal and vertical shifts on the detector as a function
of the impact position. Normal reflection filters
are applied.
Shifts are multiplied with the value
of plot2dmagnify.
The (magnified) shift is drawn according to
colourtype.
The final impact is drawn as a dot in
colour fg with size
dotsizetype.
The dimensions of the detector can be set by
detectorsize.
plot2dfixscale
Syntax: plot2dfixscale state
Default: plot2dfixscale on
- If on, the scaling of the
plot2d
graphics will be fixed. This allows comparison of plot2d's with various
filter settings.
- If off, every plot2d command will try to fill the
detector window. The scaling of the shifts on the window only depends
on the actual shifts of the (selected) reflections.
plot2dmagnify
Syntax: plot2dmagnify n
Default: plot2dmagnify 10
The magnification of the shifts in the
plot2d command.
plotc
Syntax: plotc
Plots the c-vectors in a position with all goniostat angles zero.
You will see stereographic projections
of two hemispheres. The crystal is in the center of the spheres.
Use wulffvec to define the view direction.
Related command:
plotd,
plotray.
plotd
Syntax: plotd
Plots the c-vectors in diffacting position. You will see stereographic
projections of two hemispheres. The crystal is in the center of the spheres.
Use wulffvec to define the view direction.
Related commands:
plotc,
plotray.
plotdelsig
Syntax: plotdelsig f
colourtype and
limit delsig will be restored to their original values.
plotgrid
Syntax: plotgrid
Plots a 2-dimensional histogram. The average value of the variable,
selected with colourtype var, will be
plotted. The value is averaged over an area on the detector. The size of
this area is controlled by gridsize.
Reflection filters are used.
The colour of the histogram is controlled by
colourtype.
The dimensions of the detector can be set by
detectorsize.
plotlimit
Syntax: plotlimit variable f1 f2
Sets plot limits on one of the reflection
variables.
Plot limits are used only to change the
appearance of the plots (see plot,
histcount and
histo combined with
histoautoscalehor off or
histoautoscalever off).
(Use limit if you want to filter
out reflections).
noplotlimit restores the initial
values.
Related commands:
limit,
noplotlimit,
plotlimitactual,
varlimit.
plotlimitactual
Syntax: plotlimitactual variable
Set plot limits on one of the reflection
variables
to the actual minimum and maximum value.
Related command:
plotlimit
plotray
Syntax: plotray
Plots the InRays and OutRays. You will see stereographic projections
of two hemispheres. The crystal is in the center of the spheres. The
InRays are drawn in colour inray, the outrays
follow the strong/weak colour scheme.
Use wulffvec to define the view direction.
Related command2:
plotc,
plotd.
plotrefl
Syntax: plotrefl
Plots the 2-dimensional impacts of reflections.
Reflection filters are used.
The colour of the reflections is controlled by
colourtype.
The dotsize can be set with dotsize.
The dimensions of the detector can be set by
detectorsize.
plotshell
Syntax: plotshell
Runs the rshell command, and draws a histogram
of one of the output columns (selectable by shcolumn).
precision
Syntax: precision variable n
Sets the number of decimals places of a (floating point)
reflection variable
to n.
The default value is 2 with the following exceptions:
- 1: exptim, frame,
inc, shifttau
tau and swing
- 3: shiftang, shifthorang and
shiftverang
The actual values are displayed in output full
mode with list var.
presentationscale
Syntax: presentationscale f
Intensities are presented on a arbitrary scale (the median of
the sigma's will be scaled to 2.0). You may change this scalefactor
with the presentationscale command.
Specify a negative value if you want to reset to the default. If you
want to disable the scalefactor use the
nopresentationscale command.
The presentationscale is only used for printing intensities and sigma's
to the screen
(
intensity,
sigma,
sigmaorig,
renninger,
eval14i,
eval14s,
imean,
smean,
imin,
smin,
iplus,
splus,
icalc,
iobs,
half1,
half2
)
Other scale factors control the scaling in output
files. See also the Scaling section
on the main page.
primarybeam
Syntax: primarybeam f
Sets the intensity of the primary beam. It is used to calculate
the renningerscore to
guess the intensity change.
The default value is the intensity of the strongest reflection.
f is scaled using
presentationscale (similar to intensities).
Note: increasing the value of f lowers the number of
possible renninger reflections, unless you lower the value of
renningerthreshold too.
Related commands:
renninger,
renningerthreshold,
testrenninger
ps
Syntax: ps window filename.
A hardcopy of one of the windows will be written to
filename using the ImageMagick command
import.
The format of the hardcopy file is determined by filename.
Examples:
ps 1 hard.ps
ps 2 hard.jpg
ps 2 hard.gif
psitype
Syntax: psitype n
Default: psitype 1
Set the algorithm for the calculation of reflection variable psi.
- 1
- c = cvector in diffracting position
- w = laboratory z_axis in diffracting position
- View along c
- v = laboratory x_axis in view position
- psiv = atan2(v2,v1)
- psiw = atan2(w2,w1)
- psi = psiw-psiv
- 2
- c = cvector in diffracting position
- w = laboratory x_axis in diffracting position
- View along c
- psi = atan2(w2,w1)
- 3
- c = cvector in diffracting position
- g = Vecprod(x_axis,c)
- e = Vecprod(c,g)
- u = (k-l, l-h, h-k)
- if u = (0,0,0): u = (h,-h,0)
- bring u in diffracting position
- x=scalprod(u,e)
- y=scalprod(u,g)
- psi = atan2(x,y)
- 4
- c = cvector in diffracting position
- w = laboratory x_axis in diffracting position
- View along c with w horizontal
- psi=atan2(c2,c1)
quality
Prints some numbers describing the quality of the data. This output
is also appended the various hkl files.
Quality Example
<COUNTERS INPUT="1729" INPUTUNIQUE="434" LOST="2" NOTMEASURED="0" SELECTED="1643" SELECTEDUNIQUE="432" />
<QUALITY R="0.024" RMEAS="0.027" CHI2="1.899" MEANINTENSITY="6603.374" MEANSIGMA="1.415" />
qvc
Syntax: qvc
Sets qvector filter.
You specify the order for all available qvectors. The resulting
qvp is set in the qvp filter.
You may also use qvp.
(qvc is one of the
reflection filter options).
qvp
Syntax: qvp n
If qvectors are used in the dataset, you can select specific reflections
with this command.
Use qvp all to disable
the filter.
(qvp is one of the
reflection filter options).
Related command: qvc.
r
Syntax: r
Redraw all pictures.
read
Syntax: read filename
Read reflections from filename.
Existing reflections will be discarded (use readapend
to add reflections). The default
filetype is .y. Compressed files
.Z .gz .bz2 are allowed.
Use fillgap to control the insertion of images
with no reflections.
After reading filename, the files
filename.anc and/or
any.post will be executed if they exist.
Related commands:
readappend,
readfcalc
readfcf,
readhkl,
readhklis,
readhklv,
readhk6,
readicalc
readraw,
readsad,
xyzzy.
readappend
Syntax: readappend filename
Read reflections from filename. The reflections
are added to the existing set of reflections (Use
read to read new file after removing existing
reflections first). The default
filetype is .y. Compressed files
.Z .gz .bz2 are allowed.
If a file filename.anc exists,
it will be executed.
After reading filename, the files
filename.anc and/or
any.post will be executed if they exist.
readfcalc
Syntax: readfcalc filename
Add Fcalc and Fobs from filename to the existing
reflection database. This procedure
is the same as readicalc, but all Fcalc
and Fobs values will be squared to Iobs and Icalc.
readfcf
Syntax: readfcf filename
Read reflections from fcf file filename.
Start with entering an rmat file
and the wavelength.
(You may avoid this by first reading an yfile;
the rmat and wavelength are then assumed to be correct).
Reflections will be filtered using fcffilter
using the value of _refln_observed_status.
fcf field |
any variable |
_refln_index_h |
h |
_refln_index_k |
k |
_refln_index_l |
l |
_refln_F_squared_meas |
intensity |
_refln_F_squared_sigma |
sigma
sigmaorig |
_refln_F_squared_calc |
icalc |
_refln_observed_status |
status |
The available variables after this command are
h k l intensity sigma icalc and the derived
variables ioversig nr refl reso shell sigmaorig theta
readhk6
Syntax: readhk6 filename
Read reflections from (shelx hk6 type) filename.
Every line should at least contain: h k l m n p intensity sigma
with format 6i4,2f8.2.
Start with entering an rmat file
and the wavelength.
(You may avoid this by first reading an yfile;
the rmat and wavelength are then assumed to be correct).
The available variables after this command are
delsig equiv equiv2 h imean intensity ioversig
k l m (n,p) nr refl reso rsym set shell sigma sigmaorig smean theta
readhkl
Syntax: readhkl filename
Read reflections from (shelx hklf4 or hklf5 type) filename.
Every line should at least contain: h k l intensity sigma batch
with format 3i4,2f8.2,i4.
Start with entering an rmat file
and the wavelength.
(You may avoid this by first reading an yfile;
the rmat and wavelength are then assumed to be correct).
The available variables after this command are
delsig equiv equiv2 expnr h icr imean intensity ioversig
k l nr refl reso rmeas rpim rsym set shell sigma sigmaorig smean theta
readhklis
Syntax: readhklis filename
Read reflections from filename.
Every line should contain: h k l intensity sigma
with format 3i5,2f10.2.
Start with entering an rmat file
and the wavelength.
(You may avoid this by first reading an yfile;
the rmat and wavelength are then assumed to be correct).
The available variables after this command are
delsig equiv equiv2 expnr h icr imean intensity ioversig
k l nr refl reso rmeas rpim rsym set shell sigma sigmaorig smean theta
readhklv
Syntax: readhklv filename [batchoffset]
Read intensities and sigma's from a reflections file, processed
by Sadabs and created with the Sadabs V-option. Such a file
contains the frame number of the original reflections.
Every line should contain: h k l intensity sigma batch frame
with format 3i4,2f8.2,i4,f8.1.
The optional batchoffset will be added to batch.
readhklv replaces intensity and sigmaorig.
readicalc
Syntax: readicalc filename
Add Icalc and Iobs from filename to the reflection database.
The expected format depends on the filetype.
- .fcf
Read Icalc and Iobs from an fcf file and add these two values to
the existing reflections.
Reflections will be filtered using fcffilter.
- .hkl
Read Icalc and Iobs from an hkl file and add these two values to
the existing reflections.
The input format for reflections (h,k,l,icalc,iobs) is
(3i4,2f8.2). Reflections are read until a blanc line
in encountered, or until a line with h=k=l=0 is read.
The data are stored in the
mother-database
variables pMicalc
and pMiobs
and available for all reflections as
iobs
and icalc.
Related commands:
meantype,
readfcalc.
readraw
Syntax: readraw filename
Read reflections from (saint raw file) filename.
The following variables
will be available:
delsig equiv equiv2 expnr frame h hor intensity
imean ioversig k l nr refl reso rmat rsym set shell sigma sigmaorig
sister smean swing theta ver.
readsad
Syntax: readsad filename
Read reflections from (sadabs input file) filename.
The following variables
will be available:
delsig equiv equiv2 expnr frame h hor intensity
imean ioversig k l nr refl reso rmat rsym set shell sigma sigmaorig
sister smean swing theta ver.
readsca
Syntax: readsca filename
Read reflections from (scalepack) filename.
The file has to be written by scalepack with the option
no merge original index. The first
13 lines will be ignored.
Start with entering an rmat file
and the wavelength.
(You may avoid this by first reading an yfile;
the rmat and wavelength are then assumed to be correct).
The following variables
will be available after this command:
delsig equiv equiv2 h imean intensity ioversig
k l nr refl reso rsym shell sigma sigmaorig smean theta
redirect
Syntax: redirect screen/filename
You specify the name of an outputfile; all subsequent output will be
written to this file (appended if the file already exists) . Return to
normal status with redirect screen.
region
Syntax: region n
box x1 y1 x2 y2 state
Syntax: region n
circle x y f state
Syntax: region n
reso rmin rmax state
Syntax: region n
delete
Define a region on the detector. The
BADPOS flag
is set for reflections inside this region
(state=yes) or
outside this region (state=no).
n is a region number (1-10).
- box
Specify the lower left and upper right coordinates in mm.
- circle
specify the centre of a circle and the radius in mm.
- reso
specify lower and upper resolution limits (in Å)
- delete
disable the region
Use list region to
print all regions.
(If output=debug
the flagged reflections will be printed)
reject
Syntax: reject f
Set REJECT flag for reflections if
| Ii−Imean | > f*σi.
The number of equivalent reflections must be at least 3.
Imean is calculated for the group of equivalent reflections (passing the
reflection filters).
Then for every reflection, Δ/σ is evaluated.
If the largest Δ/σ exceeds f, that reflection
is rejected. This process is repeated as long as reflections are being rejected
and the number of remaining reflections is larger than 2.
Note, the evalations of Imean and Δ/σ is done with
temporary variables. The values of the reflection variables
imean and
delsig are
not affected by the
rejection process. You may want to invoke the delsig
command after using reject.
See the section
Mean and Sigma
on the main page
for more information on Imean and sigma.
f is available in
expressions as rejectfactor.
Related commands:
rejectmanual,
rejectoverlap,
tooweak and
unreject.
rejectmanual
Syntax: rejectmanual f
Same as reject but you will be prompted to
set the REJECT flag.
Related commands:
reject,
rejectoverlap and
unreject.
rejectoverlap
Syntax: rejectoverlap
Reject all reflections with more than one overlap.
Related commands:
reject,
rejectmanual and
unreject.
renameexperiments
Syntax: renameexperiments s
Renames all experiments to sXXX with XXX the
experimentnumber.
You can use mergey to combine
several eval15 y-files, even from various datasets. There is a chance
that different sets actuallly have the same name, f.e s01f001, but
were merged together from different directories. These differences
are handled by any.
In any you normally use only the
experiment number
(expnr) or
set number (set),
so the fact that two sets have identical names is no problem.
But if these names are used during import (for example
rmatexp, rmatexpfile)
or export (for example pk or
write imagescale)
the sets are possibly identified by other programs by their name, not
their number. renameexperiments can be used
to create unique names.
renninger
Syntax: renninger state
Default: renninger off
If state=on
the variable renninger
will be activated. The variable contains the
renninger score.
Use testrenninger to investigate a
specific reflection. Use
listrefl renninger to
print all renninger reflections.
If you want to exclude reflections whose intensities can be influenced
by the renninger effect, use
- renninger on
- limit renninger 0.0 0.001 yes
Related commands:
primarybeam and
testrenninger.
renningernormalize
Syntax: renningernormalize state
Default: renningernormalize on
If on the intensities of reflections will be divided
by the length of the reflection vector (sin(θ)/λ) while
calculating the renninger score.
See the explanation on the
renninger score.
renningersimult
Syntax: renningersimult state
Default: renningersimult 0.001
Sets the criterium for the simultane occurance of renninger reflections.
See the explanation on the
renninger score.
renningerthreshold
Syntax: renningerthreshold f
Default: renningerthreshold 0.1
Sets a minimum threshold (fraction of primarybeam)
for renninger reflections.
See the explanation on the
renninger score.
require
Set reflection filters. See the flags page.
If a specific flag is not in the allow list, it will be moved from
the forbid list to the allow list.
Unselectable flags are not added to the require list.
requiv
Show rmerge for different number of equivalent reflections.
Rmerge for different number of equivalents
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
Nequi nUni nRsym Rsym Rmeas Rpim Chi2
11 2 22 0.026 0.027 0.008 1.49
10 5 50 0.022 0.023 0.007 1.63
9 8 72 0.026 0.027 0.009 2.01
8 14 112 0.022 0.024 0.008 2.20
7 17 119 0.022 0.024 0.009 2.20
6 29 174 0.024 0.026 0.011 1.50
5 44 220 0.032 0.036 0.016 2.89
4 89 356 0.026 0.030 0.015 1.79
3 104 312 0.021 0.026 0.015 1.73
2 96 192 0.015 0.021 0.015 1.91
=======================================
Total 408 1629 0.024 0.027 0.012 1.97
1 25 25
Total 433 1654
More information in the section Rsym
on the main page.
resultsig
Syntax: resultsig n
Default: resultsig 3
Sets the formula used for calculation of the sigma of the mean.
More information in the section
Mean and Sigma
on the main page.
rexp
Show rmerge for the different experiments.
Rmerge for different experiments
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
Exp set1 set2 Uni1 Uni2+ nRsym nMean R Rmeas Rpim Chi2 <I> <Sig> <I/s>
1 s01f 1 2 58 365 1202 1260 0.024 0.029 0.015 2.01 237.72 5.72 26.59
2 s02f 3 3 96 96 222 318 0.016 0.022 0.014 1.78 220.55 5.35 24.69
3 s03f 4 4 38 19 38 76 0.015 0.021 0.015 1.35 237.34 5.74 28.48
More information in the section Rsym
on the main page.
rint
Show rmerge for different intensity intervals.
The number of intervals is set by nintensity.
This command uses
imean as data source.
Negative values are ignored.
400 <intensities> from 0.041 to 4991.264
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
<Int1> <Int2> R Rmeas Rpim Chi2 nUni nRsym <I> nMean
0.0 0.8 2.058 2.464 1.310 2.68 36 140 0.30 143
0.8 5.4 0.388 0.460 0.240 3.49 35 124 2.70 129
5.4 20.2 0.104 0.126 0.069 2.95 40 151 12.59 152
20.2 44.7 0.047 0.055 0.029 1.75 39 152 32.77 152
44.7 81.3 0.026 0.030 0.015 1.09 39 155 64.27 156
81.3 123.0 0.025 0.029 0.014 1.37 40 161 99.92 162
123.0 199.4 0.019 0.022 0.011 1.08 39 151 162.20 151
199.4 330.0 0.019 0.022 0.010 1.34 40 180 256.44 181
330.0 576.1 0.017 0.020 0.010 1.36 34 136 443.92 142
576.1 4991.3 0.024 0.027 0.012 2.53 39 184 1190.90 185
j= 10 a= 2.865 slope= -0.19 abdev= 0.544
More information in the section Rsym
on the main page.
rlaue
Show rmerge for all different laue groups. You may examine the results
and decide to change the symmetry with pg.
pg R Rmeas Rpim Chi2 Uni1 Uni2+ nRsym
-1 0.020 0.027 0.017 1.88 812 369 842 -1 Triclinic
2/m 0.023 0.028 0.016 2.12 215 517 1439 2/m Monoclinic
mmm 0.024 0.027 0.012 1.97 25 408 1629 mmm Orthorhombic
4/m 0.515 0.570 0.237 8903.52 52 358 1602 4/m Tetragonal low
4/mmm 0.518 0.555 0.191 8353.94 17 258 1637 4/mmm Tetragonal high
-3 0.713 0.814 0.382 17879.88 257 390 1397 -3 Trigonal low
-3m1 0.746 0.829 0.345 16231.69 83 348 1571 -3m1 Trigonal high
-31m 0.752 0.826 0.322 15946.24 104 312 1550 -31m Trigonal high
6/m 0.684 0.732 0.252 14071.47 52 319 1602 6/m Hexagonal low
6/mmm 0.745 0.790 0.251 14346.60 34 237 1620 6/mmm Hexagonal high
R-3 0.694 0.795 0.373 20950.19 368 350 1286 R-3 Rhombohedral low
R-3m 0.755 0.840 0.344 15704.20 115 341 1539 R-3m Rhombohedral high
m3 0.634 0.679 0.234 15166.54 58 344 1596 m3 Cubic low
m3m 0.756 0.786 0.202 13243.59 9 160 1645 m3m Cubic high
More information in the section Rsym
on the main page.
Related command: rpg.
rmat
Syntax: rmat filename
Read rmat from filename. This rmat will override
the rmat from the yfile. Some variables will change:
rmatexp
Syntax: rmatexp name (9 rmatrix parameters)
Set the orientation matrix for a specific experiment.
name is the experiment name
(not the experiment number). The matrix is filled with the 9 arguments,
row by row, so the order is (1,1) (1,2) (1,3) (2,1) (2,2) (2,3) (3,1)
(3,2) (3,3).
The global transformation matrix and laue/pointgroup settings are used.
This experiment-dependent rmat is used in
Related commands:
list rmatexp,
rmatexpfile,
useintegratedcell,
write rmatexp.
rmatexpfile
Syntax: rmatexpfile name rmatfile
Set the orientation matrix for a specific experiment.
name is the experiment name
(not the experiment number). The matrix is read from rmatfile.
The global transformation matrix and laue/pointgroup settings are used.
This experiment-dependent rmat is used in
Related commands:
list rmatexp,
rmatexp,
useintegratedcell,
write rmatexp.
rpg
Show rmerge for all different point groups. You may examine the results
and decide to change the symmetry with pg.
More information in the section Rsym
on the main page.
Related command: rlaue.
rq
Show rmerge for all different qvector combinations.
More information in the section Rsym
on the main page.
rset
Show rmerge for the different sets.
Rmerge for different sets
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
Set Uni1 Uni2+ nRsym nMean R Rmeas Rpim Chi2 <I> <Sig> <I/s>
1 s01f001 84 292 886 970 0.025 0.030 0.016 2.07 236.82 5.70 26.49
2 s01f002 111 79 179 290 0.016 0.021 0.014 2.15 240.72 5.79 26.93
3 s02f001 96 96 222 318 0.016 0.022 0.014 1.78 220.55 5.35 24.69
4 s03f001 38 19 38 76 0.015 0.021 0.015 1.35 237.34 5.74 28.48
More information in the section Rsym
on the main page.
rshell
Show completeness and rmerge for the different resolution shells.
The number of shells is controlled
by nshell and autoeshell.
Completeness and Rmerge for Shells
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
theta from 0.0 to 27.485
Sh Theta Reso Meas Equi Obs Mis Lost Total Perc Cum Uni1 Uni2+ Nrsym Redun Rsym Rmeas Rpim Chi2
1 12.37 1.659 247 19 266 0 0 266 100.0 100.0 0 55 317 5.76 0.027 0.029 0.012 2.56
2 15.66 1.317 210 68 278 0 0 278 100.0 100.0 1 46 247 5.37 0.024 0.027 0.012 2.02
3 18.00 1.150 167 91 258 0 0 258 100.0 100.0 3 41 186 4.54 0.021 0.023 0.011 2.15
4 19.88 1.045 175 119 294 0 0 294 100.0 100.0 1 44 183 4.16 0.020 0.023 0.011 1.66
5 21.49 0.970 122 114 236 0 0 236 100.0 100.0 2 36 120 3.33 0.025 0.029 0.015 1.96
6 22.91 0.913 150 148 298 0 2 300 99.3 99.9 2 44 148 3.36 0.017 0.021 0.011 1.03
7 24.19 0.867 119 125 244 0 0 244 100.0 99.9 0 35 119 3.40 0.026 0.031 0.017 1.39
8 25.37 0.829 128 184 312 0 0 312 100.0 99.9 8 40 120 3.00 0.021 0.026 0.015 2.43
9 26.46 0.798 87 129 216 0 0 216 100.0 99.9 2 30 85 2.83 0.028 0.034 0.020 1.52
10 27.48 0.770 110 174 284 0 0 284 100.0 99.9 6 37 104 2.81 0.030 0.038 0.022 2.11
======================================================================================================
27.48 0.770 1515 1171 2686 0 2 2688 99.9 99.9 25 408 1629 3.99 0.024 0.027 0.012 1.97
Resolution 10.247-0.77 (0.798-0.77)
Rsym 0.024 (0.03) Rmeas 0.027 (0.038) Rpim 0.012 (0.022) Completeness 99.926 (100.0)
To evaluate the completeness, the whole sphere (in space group P1) of
possible reflections is investigated.
- The lower limit of theta is limit theta, if set.
Otherwise the lower theta
is set to zero, unless
thminzero=off, in which case
the lowest encountered theta value of the whole dataset is used.
- If a specific reflection has been measured at least one time (and
passes the reflection filters)
increment Meas.
- If a specific reflection has not been measured (or is rejected by the
reflection filters), but a laue-equivalent
(of the actual pg)
is available, increase Equi.
- Obs is the sum of Meas and Equi
- Unavailable reflections are counted in Mis
- Available reflections not passing the reflection
filters are counted in Lost
- Total = Obs + Mis + Lost
- Perc = 100.0*Obs/Total
- Cum = 100.0*Obs/Total for the sum of the counters
starting with the first shell.
- Uni1 is the number of unique reflections which have been
measured one time
- Uni2+ is the number of unique reflections with at least
2 equivalents.
Of these groups of equivalents, the various
Rsym are calculated.
- Nrsym is the total number of reflections with at least
2 equivalents.
- Redun = Nrsym/Uni2+
A small interpretation from the last line in the table:
- 2 reflections lost (because of filter
settings).
- 25 reflections measured one time.
- 408 unique reflections have been measured 1629 times. The average
redundancy for these reflections is 4 (1629/408)
If dict=on, some key/value
pairs will be printed.
dict on
rshell
Completeness and Rmerge for Shells
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
theta from 0.0 to 27.485
Sh Theta Reso Meas Equi Obs Mis Lost Total Perc Cum Uni1 Uni2+ Nrsym Redun Rsym Rmeas Rpim Chi2
1 12.37 1.659 247 19 266 0 0 266 100.0 100.0 0 55 317 5.76 0.027 0.029 0.012 2.56
2 15.66 1.317 210 68 278 0 0 278 100.0 100.0 1 46 247 5.37 0.024 0.027 0.012 2.02
3 18.00 1.150 167 91 258 0 0 258 100.0 100.0 3 41 186 4.54 0.021 0.023 0.011 2.15
4 19.88 1.045 175 119 294 0 0 294 100.0 100.0 1 44 183 4.16 0.020 0.023 0.011 1.66
5 21.49 0.970 122 114 236 0 0 236 100.0 100.0 2 36 120 3.33 0.025 0.029 0.015 1.96
6 22.91 0.913 150 148 298 0 2 300 99.3 99.9 2 44 148 3.36 0.017 0.021 0.011 1.03
7 24.19 0.867 119 125 244 0 0 244 100.0 99.9 0 35 119 3.40 0.026 0.031 0.017 1.39
8 25.37 0.829 128 184 312 0 0 312 100.0 99.9 8 40 120 3.00 0.021 0.026 0.015 2.43
9 26.46 0.798 87 129 216 0 0 216 100.0 99.9 2 30 85 2.83 0.028 0.034 0.020 1.52
10 27.48 0.770 110 174 284 0 0 284 100.0 99.9 6 37 104 2.81 0.030 0.038 0.022 2.11
======================================================================================================
27.48 0.770 1515 1171 2686 0 2 2688 99.9 99.9 25 408 1629 3.99 0.024 0.027 0.012 1.97
Resolution 10.247-0.77 (0.798-0.77)
Rsym 0.024 (0.03) Rmeas 0.027 (0.038) Rpim 0.012 (0.022) Completeness 99.926 (100.0)
!!any-resolution 10.247-0.77
!!any-resolutionlastshell 0.798-0.77
!!any-rsym 0.024
!!any-rsymlastshell 0.03
!!any-rmeas 0.027
!!any-rmeaslastshell 0.038
!!any-rpim 0.012
!!any-rpimlastshell 0.022
!!any-completeness 99.926
!!any-completenesslastshell 100.0
!!any-reflections 1654
!!any-unique 433
!!any-uniquelastshell 43
!!any-redundancy 3.993
!!any-redundancylastshell 2.811
Related commands:
plotshell,
shcolumn,
thminzero.
rvar
Syntax: rvar variable
Show rmerge for variable shells
(See list of variables).
The number of shells is set by nshell.
More information in the section Rsym
on the main page.
sadabs
Writes sadabs input file. The format depends on the availability
of q-vectors.
Without qvectors: Sadabs 2.0.3 format, with qvectors: Saint ram format
- If the file is not yet opened, it will be created
(the filename may be set using hklfilename).
- The batchnumber in the file is controlled by
hklf4batch.
- Cell parameters are written to the file
filename.par.
- The intensities of the reflections will be scaled
(if necessary) to accomodate the shelx F8.0 format. You may
use the hklf4sf command to fix the scale factor.
- msa has to be set to 0.0 before
using the sadabs command.
Note: sadabs ignores reflections with an I/σ exceeding 1000.0.
To avoid this elimination, msa will be set to 0.001
if necessary.
- You can use frameoffset to add a constant
to the framenumber in the output file.
- The file will be closed (in contrary to the other
reflection file generate commands)
See also the Reflection output section on the
main page.
Related commands:
readhklv
save
Syntax: save filename
Writes the program status to filename.
The file will have filetype .anc.
You can later reload the saved settings with the
@filename
Related command:
write
set
Enable/Disable sets. Experiments are complete scans. Experiments
are divided into sets.
Syntax: set all/none/next/prev/curr (+/-)n ....
Use a plus sign to allow sets and a minus sign to forbid them.
You may enter more selections in one line of input.
Use nop or ; to
signal the end of options if you want to put more commands
on a single line.
all |
enable all sets |
none |
disable all sets |
invert |
toggle all sets |
curr |
evaluates to the last selected set |
next |
evaluates to the next set |
prev |
evaluates to the previous set |
last |
evaluates to the last set |
n |
sets/unsets set n |
Set Examples
- set all -3
use all sets except nr 3
- set none 3 4
use only sets 3 and 4
- set none 2 to 5
use only sets 2, 3, 4 and 5
- set all -3 to 6
use all sets except 3, 4, 5 and 6
- set -3 to 6 4
disable sets 3, 5 and 6 (first disable 3, 4, 5 and 6; then enable 4)
- set none 1; plotrefl
select only set 1. Plot selected reflections.
- set -curr +next; plotrefl
disable current set and enable the next one. Then plot
selected reflections.
- \repeat 5 set -curr +next; plotrefl \wait 5
Plot selected reflections for the next 5 sets with a 5 second
interval.
(set is one of the
reflection filter options).
Use list set to
print the names and status of all sets.
You could also use limit set, but
the set
command is moch more flexible.
Related commands:
experiment,
forbidsets,
setminrefl.
setminrefl
Syntax: setminrefl n
Disable sets with a number of reflections less than n
shcolumn
Syntax: shcolumn colname
Default: shcolumn Perc
colname is one of
Theta/Reso/Meas/Equi/Obs/Mis/Lost/Total/Perc/Cum/Uni1/Uni2+/Nrsym/Rsym/Rmeas/Rpim/Chi2.
Selects the column of the output of the
rshell command.
A histogram of the values in this column can be drawn with
plotshell.
show
Shows information of the various variables in the program.
Use save to write this information to a file.
Use status to display information on the
data.
colour main darkgoldenrod fg white bg black label yellow ticks grey30 text yellow strong red weak skyblue negative blue inray green
histobar1 grey40 histobar2 white histoline1 darkgreen histoline2 green var1 #0000AA var2 #0000FF var3 #0077FF var4 #00FFFF
var5 #00FF99 var6 #00FF00 var7 #99FF00 var8 #FFFF00 var9 #FF7700 var10 #FF0000 var11 #AA0000
colourbw off colourinvert off colourscheme rainbow colourtype weak
abort off autoshell off bgqual 1.0 correct normal coscell 1 cosines on dotsize 1 dotsizemax 10 dotsizetype fixed dump off euler off
fillgap on formulasigext 1 formulasigint 2 frameoffset 0.0 full good gridsize 10 histoautoscalehor off histoautoscalever on histobar on
histobarfraction 0.95 histoprint off hklf4batch expnr hklfilename shelx icrmin 2 icrscale 2.0 icrtype intensity kfactor 1.0 label 5
meantype runtime modprint 10000 more on msa 0.02 msapower 0.02 ncolourvar 11 nintensity 10 nnomore 500 nshell 10 omitsum on output normal
pkrestfrac 0.0 plot2dmagnify 10.0 primarybeam 5097.8667 renninger off renningerthreshold 0.1 renningernormalize on
renningersimult 0.001 resultsig 3 shcolumn Perc sigma 2.5 sigmafactor 1.0 speedcorrect on text on thminzero on trace off
weightint 3 weightsig 3
window border 0 fraction 0.7 padding 1 pos 168 391 size 1175 391
read example1.y.Z forbid EDGEVER+EDGEHOR+EDGEROT+BADBG+MAXSHIFT allow GOOD+WEAK+NEGATIVE presentationscale 0.0322754
sigma
Syntax: sigma f
Sets WEAK flag on reflections with I/sig < f.
f is available in
expressions as weakfactor.
Related command:
tooweak.
sigmafactor
Syntax: sigmafactor sf
Default: sigmafactor 1.0
Recalculate the sigma's of all reflections.
See the Mean and Sigma section
on the main page.
sf is available in
expressions as sigmafactor.
Related commands:
automsa,
kfactor,
msa,
msapower.
speedcorrect
Syntax: speedcorrect state
Default: speedcorrect on
- on
While reading reflections from the input file, intensities and sigma's
will be corrected for integration time, rotation angle and generator
settings.
- off
No corrections will be applied.
See also the Scaling section on the
main page.
Related command:
tooweak
spgr
Syntax: spgr n
Sets the spacegroupnumber to n.
If n is not specified, the program proposes
some numbers, based on the current Bravais lattice type and pointgroup.
Related command:
mtz
status
Shows information on the data and the filters
Use show to display information on various program
variables.
RMAT 1 ir
RMAT DMAT
0.1002587 0.0943176 0.0654995 2.2806106 -3.5076888 2.2895818
-0.1542028 -0.0265447 0.0642691 3.4023077 -0.9575436 -4.8559546
0.1006531 -0.1346151 0.0332188 6.8771739 6.7479835 3.4878356
Determinant: 0.3406898E-02 293.5221
Cell: 4.76941 6.00607 10.24675 90.0000 90.0000 90.0000 V= 293.52
Sigma 0.0006 0.0008 0.0014 0.011 0.009 0.012 Volume 0.06
Bravais=P pg=mmm
Lambda=0.71073 pointgroup=mmm
reflections=1729 unique=434 weak=282 negative=11 badbg=16
Unique: H from -6 to 0 K from -7 to 0 L from -13 to 0
Full Sphere: H from -6 to 6 K from -7 to 7 L from -13 to 13
Max equivalents 11 datalimit theta 1.9875 27.4845 datalimit reso 0.77 10.2467 datalimit chi 0.0002 73.1994
datalimit psi -176.6129 179.6856
datalimit intensity -6.7563 5097.8667 datalimit duration 1.0 4.95
3 Experiments 4 Sets 231 Images
presentationscale 0.0322754
Input Flags: GOOD WEAK EDGEVER EDGEHOR EDGEROT BADBG NEGATIVE BADUNIF MAXSHIFT
Always Forbidden: BADUNIF
Selectable: GOOD WEAK EDGEVER EDGEHOR EDGEROT BADBG NEGATIVE MAXSHIFT
Forbid: EDGEVER EDGEHOR EDGEROT BADBG MAXSHIFT
Allow: GOOD WEAK NEGATIVE
Require: NONE
step
Syntax: step variable command(s)
Starts a loop over variable. Each cycle will set
new limits on
this variable. The first lower limit is the minimum
value of variable. The first upper limit is the
lower limit plus stepsize. Then the
command(s) will be executed.
Before each cycle you are asked how to proceed. Possible answers are:
- no : skip
- quit : stop stepping
- all : execute further with no prompts
- delay : prompts for a delaytime. execute
further with no prompts but with delay
- any other: execute
The (sub)strings stepname and stepnumber
in command(s) will be substituted:
- variable=expnr and
stepsize=1: stepname becomes the name of the experiment.
- variable=set and
stepsize=1: stepname becomes the name of the set.
- otherwise: stepname becomes the name of the
variable with the stepnumber attached to it.
Related command: stepfull.
stepfull
Syntax: stepfull variable inc start stop command(s)
Starts a loop over variable. Each cycle will set
new limits on
this variable. The first lower limit is the
start. The first upper limit is the
lower limit plus stepsize. Then the
command(s) will be executed.
Before each cycle you are asked to continue. Possible answers are:
- no : skip
- quit : stop stepping
- all : execute further with no prompts
- any other: execute
The (sub)strings stepname and stepnumber
in command(s) will be substituted:
- variable=expnr and
stepsize=1: stepname becomes the name of the experiment.
- variable=set and
stepsize=1: stepname becomes the name of the set.
- otherwise: stepname becomes the name of the
variable with the stepnumber attached to it.
Related command: step.
testrenninger
Syntax: testrenninger n
n is a reflection number.
Investigate renninger reflections for reflection n.
Use listrefl renninger
to print all renninger reflections.
Related command: renninger.
text
Syntax: text state
Default: text on
Controls the drawing of textlabels in various plots.
thminzero
Syntax: thminzero state
Default: thminzero on
rshell counts missing and lost reflections in order
to evaluate the completeness of the dataset. Reflections are missing if
they are not available in the input file. Reflections get lost if flags
and filters rule them out. Also, listrefl mis
prints a list of unavailable reflections.
The choice of the ultimate lower resolution limit is arbitrary: the lowest
encountered value from the input file, or zero.
- off. Use the lowest theta value from the
reflection file as lower limit for theta.
- on. Use 0.0 as minimum theta value.
If a lower theta limit has been set (limit theta or
limit reso), thminzero
is not used.
title
Syntax: title text
Default: title ""
Draw text in the top of the graphics window
tooweak
Syntax: tooweak speed/exp/exptest/ioversig/ioversigtest
This command tries to identify reflections which worsen the
average I/σ of a group of equivalent reflections.
The idea is that some reflections do not significantly contribute to the final
I/σ.
- speed nmin
Try to remove reflections from the fastest scan.
The number of equivalent reflections in the non-fastest scan should be at
least nmin.
The average I/σ for a each group of equivalent reflections is calculated.
Then the reflections of the fastest scan are excluded one by one and
I/σ is re-evaluated. If I/σ improves, the excluded reflections
are flagged TOOWEAK.
Note: only useful if different exposure times have been used.
- exp nmin
At least nmin reflections should remain.
Try to remove reflections from the final experiment.
The average I/σ for a each group of equivalent reflections is calculated.
Then the reflections in the latest scan are excluded and
I/σ is re-evaluated. If I/σ improves, the excluded reflections
are flagged TOOWEAK.
- exptest h k l nmin
Similar to exp, but only for the (equivalent) reflections
with indices nmin h k l
- ioversig
Try to remove reflections based on I/σ.
At least nmin reflections should remain.
The average I/σ for a each group of equivalent reflections is calculated.
Then the reflections with the worst I/σ are excluded one by one and
I/σ is re-evaluated. If I/σ improves, the excluded reflections
are flagged TOOWEAK.
- ioversigtest h k l nmin
Similar to ioversig, but only for the (equivalent) reflections
with indices nmin h k l
total
Prints the global rmerge.
Rsym=0.024 Rmeas=0.027 Rpim=0.012 Chi2=1.968 nRsym=1629 Unique1=25 Unique2+=408
<I>=234.399 <s>=5.651 <I/s>=26.314 <I>/<s>=41.481 nMean=1654
More information in the section Rsym
on the main page.
twinabs
Syntax: twinabs
Writes twinabs formatted reflections to the reflection file.
- If the file is not yet opened, it will be created
(the filename may be set using hklfilename).
- The batchnumber in the reflection file is controlled by
hklf4batch.
- The twin component
numbers are set by hklf5batch.
- The intensities of the reflections will be scaled
(if necessary) to accomodate the shelx F8.0 format. You may
use the hklf5sf command to fix the scale factor.
- Cell parameters are written to the file
filename.par.
- msa has to be set to 0.0 before
using the twinabs command.
- Use hklclose to close the reflection file.
See the TwinabsPage
for an example of handling twin datasets.
See also the Reflection output section on the
main page.
twinlaw
Syntax: twinlaw state
Default: twinlaw off
If twinlaw=on and more
than one matrix is available, the hkl transformation matrix will be
written when a shelx outputfile is created.
unlimit
Syntax: unlimit variable
Unset the limits (set by limit) on one of the
reflection variables.
unlimithkl
Syntax: unlimithkl
Unset the limits (set by limithkl).
This is an alias for unlimit hb im.
unlimitrecip
Syntax: unlimitrecip
Unset the limits (set by limitrecip).
This is an alias for unlimit h k l.
unreject
Removes the reject flags from all reflections.
Related commands
reject,
rejectmanual and
rejectoverlap.
update
Make some constants available to the built-in calculator
(used by expression or
\evaluate)
The constant names and values will also be printed.
Related commands:
\addeval and
\listeval
useicalc
Syntax: useicalc state
Default: useicalc off
If set to on, the various rmerge calculations on
a group of equivalent reflections will use variable
icalc instead of
determining the mean intensity.
Related command:
readicalc.
useintegratedcell
Syntax: useintegratedcell state
Default: useintegratedcell on
During integration with eval15, every experiment uses its
own orientationmatrix. These matrices are kept in the y-files.
If any reads a y-file, only the the first
orientationmatrix is used. The remaining ones are ignored.
With useintegratedcell on, all orientation
matrices are read. They are not used by any,
but the commands
list rmatexp and
write rmatexp use them.
varexp
Syntax: varexp name f
Sets the value of
href="anydatastructure.html#varexp">varexp for a specific
experiment, identified with the experimentname (not
the experiment number).
Related command:
varexpalias
varexpalias
Syntax: varexpalias name
Default varexpalias varexp
Change the name of reflection variable
varexp.
Related command:
varexp.
varlimit
Syntax: varlimit variable f1 f2
Sets the varlimits on one of the reflection
variables.
Varlimits are used only for the colourscheme of the
colourtype var command
(Use limit if you want to filter
out reflections).
You may use novarlimit to restore the initial
values.
Related commands:
limit,
novarlimit,
plotlimit.
wait
Syntax: wait f
Wait f seconds. Useful in repeated constructs like
\repeat and
step.
weightint
Syntax: weightint n
Default: weightint 3
Sets the intensity weights for the calculations of mean and sigma.
1 |
unit weights |
2 |
statistical weights |
3 |
shelx weights |
More information in the section
Mean and Sigma
on the main page.
weightsig
Syntax: weightsig n
Default: weightsig 3
Sets the weight for sigma merging.
1 |
unit weights |
2 |
statistical weights |
3 |
statistical weights |
More information in the section
Mean and Sigma
on the main page.
window
Syntax: window border/fraction/info/log/padding/pos/size/title
window border
Syntax: window border n
Default: window border 0
Sets the border margin of the subwindows to n pixels.
Related commands:
window fraction,
window padding,
window size.
window fraction
Syntax: window fraction f
Default: window fraction 0.7
Sets the horizontal window dimension to fraction f of the
screen. The vertical dimension is 2/3 of the horizontal.
Related commands:
window border,
window padding,
window size.
window info
Syntax: windowinfo
Prints window information.
window log
Syntax: windowlog n
Sets the amount of graphics related output.
0 |
Emergency |
1 |
Alert |
2 |
Critical |
3 |
Error |
4 |
Warning |
5 |
Notice |
6 |
Info |
7 |
Debug |
window padding
Syntax: window padding n
Default: window padding 1
Sets the padding margin between subwindows to
n pixels.
Related commands:
window border,
window fraction,
window size.
window pos
Syntax: window pos nx ny
Sets the position (in pixels) of the main window.
The defaults depend on the screen dimensions and the size
of the window.
window size
Syntax: window size nx ny
Sets the window dimensions to
nx * ny pixels. The default
dimensions are set with window fraction.
Related commands:
window border,
window fraction,
window padding.
window title
Syntax: window title s
Sets the window title to s
write
Syntax: write bfactor/detalign/goniostat/imagescale/impact/rmatexp/rmatX/scanplot name
- bfactor
Writes bfactors for all experiments to name (bfactor.txt).
This file can be read with bfactorread.
In peakref the command is readbfactor.
- detalign
Writes detector position and orientation to name
(detalign.vic).
- goniostat
Writes goniostat information to name
(goniostat.vic).
- imagescale
Write image scale factors for all frames to name (imagescale.txt).
This file can be read with imagescaleread.
In peakref the command is readimagescale.
- impact
Creates name (impact.dat)
For each mother reflection with more than 2 daughters the mean intensity
is calculated. To the file is written for each individual reflection
impacthor impactver del/sig reflnr
where del = I−Imean.
- rmatexp
Write orientation matrices for all experiments to a single
file name.
(default is rmatexp.pec).
This file can be executed (using @rmatexp.pec)
by peakref and any.
Related commands:
list rmatexp,
rmatexp,
useintegratedcell.
- rmatX
Writes the rmatX to name. If X is not specified then X=1.
name defaults to the name of the original rmatfile.
- scanplot
Writes the start- and endpositions of the various experiments to
name (setup.scc). This file is to be used by the
scanplot program.
Related command:
save
wulffexpand
Syntax: wulffexpand state
Default: wulffexpand off
Controls which reflections will be plotted with the
plotc and plotd commands.
- off
Draw only available reflections.
- on
Expand the reflections using pointgroup (pg) symmetry.
This will generate a whole sphere a reflections. The
colourtype of the first member
of a group of equivalent reflections determines the
colourtype of all members of that group.
wulffhor
Syntax: wulffhor none
Syntax: wulffhor hkl/uvw/xyz f1 f2 f3
Default: wulffhor none
Defines the horizontal axis in the wulff projection windows.
In the wulff projection windows, the crystal is in the center, and you
will see the reciprocal lattice (plotc and
plotd) or the
In- and Outrays (plotray) with respect
to the crystal.
- wulffhor none. No specific horizontal axis
- wulffhor hkl. Specify 3 reflection indices
- wulffhor uvw. Specify 3 indices in terms
of the direct lattice.
- wulffhor xyz. Specify coordinates in the
laboratory system.
wulffvec
Syntax: wulffvec hkl/uvw/xyz f1 f2 f3
Default: wulffvec xyz 0 0 1
Defines the view direction of the wulff projection windows.
In the wulff projection windows, the crystal is in the center, and you
will see the reciprocal lattice (plotc and
plotd) or the
In- and Outrays (plotray) with respect
to the crystal.
- wulffvec hkl. Specify 3 reflection indices
- wulffvec uvw. Specify 3 indices in terms
of the direct lattice.
- wulffvec xyz. Specify coordinates in the
laboratory system.
Note, this definition does not use any goniostat
information. So think of the crystal as it is mounted on the goniometerhead,
and put the goniometer head with the crystal in front of you. Then
wulffvec xyz 0 0 1 will point to the zenith.
You may want to use wulffhor to specify the horizontal axis.
xyicr
Syntax: xyicr filename
Repeat the previous icr command, and write the
x,y coordinates of the plot to filename.
xyplot
Syntax: xyplot filename
Repeat the previous plot command, and write the
x,y coordinates of the plot to filename.
xyzzy
Syntax: xyplot read/write filename
Create or read a binary data file.
- Binary files can be read much faster than a normal yfile,
but they lack all the flexibility of the xml format.
- Binary files are useful with very, very large datasets.
- xyzzy write creates
filename.xyzzy
- xyzzy read reads
filename.xyzzy
Any
EVPY Suite Overview