Any
Introduction
any is a program to analyse .y-files (output
files from eval14 and
eval15) and
to create various reflection files.
Start
Start any with the command any.
Once any is started, it checks whether the file
any.init exists in the directory where
any is installed. If this file exists, it will be
executed. Next, the program checks whether the file
any.init exists in the current directory.
If this file exists, it will be executed.
Read an y-file with read.
Commands may be specified on the command line.
Note: any can handle multiple yfiles, but processes one file at a time.
Help
With the command help a webbrowser
(iceweasel, firefox, mozilla, netscape) will be started displaying this page.
(You may set an environment variable BROWSER to specify a specific browser).
Also available are:
Analysis
The analysis options of the program consist of:
- R-values
Rmerge values (for the selected pointgroup) are calculated for
various selections of reflections.
You can use reflection filters (more about this in the
Reflection Filters paragraph).
Available r-values are:
- requiv. Rmerge/equivalents
- rexp. Rmerge/Experiment
- rint. Rmerge/IntensityShell
- rlaue. Rmerge/LaueGroup
- rpg. Rmerge/PointGroup
- rset. Rmerge/set
- rshell. Rmerge/ResolutionShell
- rvar. Rmerge/Variable
- total. Global Rmerge
- Intensity Control
The command icr creates a plot
with the deviations of all reflections from their mean value
as a function of framenumber, impact position or any other
reflection variable.
This is useful to inspect scaling, decay or absorption.
- Plot
To create scatter plots of 2
reflection variables use
the command plot.
- Histograms
To create histograms of 2
reflection variables use
the commands
histo and
histcnt.
- Plot Shifts
Use plot2d to display
the (magnified) horizontal and vertical shifts as a function of the impact
position on the detector.
- Plot del/sig
plotdelsig
to display
reflections with extreme del/sig values as a function of the impact position
on the detector.
- Plot 2d histogram
Use plotgrid the plot the
average value of a
reflection variable.
The variable is selected by
colourtype var and the colours
can be controlled with varlimit.
- Plot impacts
Use plotrefl to display selected
reflections as a function of the impact position on the detector.
The colours reflect weak/nonweak reflections. To use one of the
reflection variables for
the colour scheme, use
colourtype var and
varlimit or
colourtype sign.
Reflection Filters
There are various ways to select/deselect reflections:
List reflections
- hkl. Show specific reflection
- listrefl. Show reflections
- asym. Show reciprocal lattice
- plotrefl. Plots reflections
- frame. Show reflections within one frame
Prepare Refinement
A pkfile can be created with the
pk command. This file can
be used for a refinement of experimental parameters, cell dimensions
and crystal orientation by peakref.
Reflection output
Various types of reflection files can be written. Specify the filename
with hklfilename,
write the data with
hklf4,
hklf4merge,
hklf5,
hklf5merge,
mtz,
sadabs or
twinabs
and close the file with hklclose.
- hklf4
This is the normal shelx hklf4 format. You may change the shelx
batch number with hklf4batch.
The default filetype is .hkl.
- hklf4merge
Reflections will be merged. No direction cosines are written.
This is the normal shelx hklf4 format. You may change the shelx
batch number with hklf4batch.
The default filetype is .hkl.
- hklf5
You must specify the shelx batch number for all lattices with
hklf5batch.
If an overlapping reflection has an non-overlapping equivalent the
overlapping reflection can be discarded (see
omitsum).
The default filetype is .hklf5.
See the TwinPage for an example of
handling twin datasets.
- hklf5merge
Reflections will be merged. No direction cosines are written.
You must specify the shelx batch number for all lattices with
hklf5batch.
If an overlapping reflection has an non-overlapping equivalent the
overlapping reflection can be omitted (see
omitsum).
The default filetype is .hklf5.
See the TwinPage
for an example of handling twin datasets.
- mtz
The output file can be read by CCP4 programs.
The default filetype is .mtz.
- sadabs
The output file can be read by SADABS version 2.03.
The default filetype is .sad.
- twinabs
The output file can be read by TWINABS.
The default filetype is .sam. You have to set
msa to 0.
See the TwinabsPage for an
example of handling twin datasets with twinabs.
Also available is the jana
command to create
a jana-type reflection file for incommensurate datasets.
-
Direction cosines will be written if
cosines
is on. Direction cosines are calculated relative
to the main lattice. You may change this with
coscell.
-
If more than one lattice is available, you can set the variable
twinlaw to on.
When the shelx file is created, you are prompted for the relation between
the lattices. The hkl transformation matrix will then be written to the
shelx.par file.
- eval15 and eval14-like data can be combined using
hybrid.
Scaling
- Experiment dependent scaling.
Reflection intensities and sigma's will be scaled if
speedcorrect=on.
A scale factor
is calculated for each experiment. The scale factor depends on:
scan interval(in degrees), exposure time(in seconds) and
generator current(in mA).
ScaleFactorexpnr =
1000*abs(ScanInterval)/ExposureTime/GeneratorCurrent
With these scalefactors, intensities and sigma's are brought to
a comparable scale.
Use list exp
to print the scale factors.
The Experiment dependent scalefactors are applied to the reflection
intensities and sigma's while reading the input file.
- Presentation
For presentation of intensities and sigma's on the screen another
scale factor
(see presentationscale)
is used to set the median of the sigma's to 2.0. This avoids too large
or too small numbers. If the presentationscale is switched off
(nopresentationscale),
the presentationscale will be set to 1.0/ScaleFactor of Experiment 1.
In that case the presented intensities and sigma's will not
differ from the values in the input file (unless different experiments have
different scale factors).
(Use presentationscale -1
to force the re-evaluation of the (default) presentationscale).
- Output
Scaling is also done in the various output files.
The scale factor depends on the values of the Maximum Intensity,
Minimum Intensity and Minimum Sigma.
If 0.0001<Minimum Sigma<1.0 the scale factor
is set to 1.0/MinimumSigma, and the values of MaximumIntensity
and MinimumIntensity are multiplied with this scale factor.
Then, if the MaximumIntensity>9999999.0 or the
MinIntensity<-999999.0, a final scale factor is calculated
to avoid too large numbers. It is possible to set the scale factors
with hklf4sf or
hklf5sf.
Mean and Sigma
Each reflection has an intensity (I)
and a sigma (σ). The
msa,
kfactor,
msapower and
sigmafactor
commands will increase
σ.
σ = kfactor *
√( (sigmafactor*σorig)² +
( msa*intensity)msapower )
There are various ways to calculate mean values and corresponding sigma's
of a group of n (equivalent) reflections.
The following commands allow you to control
these calculations.
- weightint
Sets the weight v for intensity calculations.
- unit weight
v = 1.0
- statistical weight
v = 1/σ²
- shelx weight
v = max(I/σ,3.0)/σ
Imean = ∑
(v·I) /
∑(v)
- weightsig
Sets the weight w for sigma calculations.
- unit weight
w = 1.0
- statistical weight
w = 1/σ²
- statistical weight
w = 1/σ²
- formulasigext
Sets the formula for external sigma calculations.
- shelx scheme
SigExt=√{∑
(w·σ²) /
[n·∑(w)] }
- crystals, sortav scheme
SigExt=√{ ∑
(w·σ²) /
∑(w) }
Note: SigExt(1) =
(1/√n)·SigExt(2)
- formulasigint
Sets the formula for internal sigma calculations.
Note: Imean = ∑
(v·I) /
∑(v)
- standard scheme.
SigInt=√{ ∑
[w·(I−Imean)²] /
[(n−1)·∑(w)] }
- shelx scheme
SigInt=∑
[abs(I−Imean)] /
[n·√(n−1)]
- crystals scheme
SigInt=√{ ∑
[w·(I−Imean)²] /
∑(w) }
- sortav scheme
SigInt=√{
n·∑
[w·(I−Imean)²] /
[(n−1)·∑ (w)] }
Notes:
SigInt(4) = SigInt(1) ·
√n
SigInt(3) = SigInt(1) ·
√(n−1)
- resultsig
Sets the final sigma
- use SigInt
- use SigExt
- use maximum of SigInt and SigExt
- use minimum of SigInt and SigExt
The defaults are:
weightint 3 weightsig 3 formulasigext 1 formulasigint 2 resultsig 3 (shelxscheme)
various schemes
statscheme |
2 |
2 |
- |
1 |
1 |
shelxscheme |
3 |
3 |
1 |
2 |
3 |
unitscheme |
1 |
1 |
- |
1 |
1 |
sortavscheme |
2 |
2 |
2 |
4 |
1 |
martinscheme |
2 |
2 |
1 |
1 |
3 |
These schemes are used by
delsig,
hkl,
hklf4merge,
hklf5merge,
hklf6merge,
icr,
list intensity,
listrefl,
mtz merge,
mtz separate,
reject,
rejectmanual
and
tooweak.
Rsym
Three merge-R's are calculated for every reflection with more than
one equivalent.
Imean is calculated
as described before:
Imean = ∑
(v·I) /
∑(v)
with v set by
weightint.
Use meantype to control
which values of Imean and
nmean will be used.
- Rsym
Rsymhkl =
∑(|Ii−Imean|)/
∑(Ii)
- Rmeas
Rmeashkl =
√(n/(n−1)) ·
∑(|Ii−Imean|)/
∑(Ii)
- Rpim
Rpimhkl =
√(1/(n−1)) ·
∑(|Ii−Imean|)/
∑(Ii)
- χ²
χi² =
nmean/(nmean−1)·(Ii−Imean)/σi )²
nmean = number of equivalent reflections
χ² = ∑χi² / n
Note: the calculations of Rsym,
Rmeas and Rpim involve a division by
∑(Ii). This sum can be zero or even
negative. This can lead to absurd high or negative values. To avoid these,
the denominator in these expressions is actually
MAX(1.0,∑(Ii)).
Manipulations
- pg. Change point group.
- msa. Increase sigma.
- badmother. Set BADMOTHER flag.
- bgqual. Set bad background flag
- reject.
Reject reflections based on mean intensity and sigma
- rejectman.
Manual reject reflections.
- sigma. Set WEAK and NEGATIVE flags.
- unreject. Remove REJECT flags.
Graphics
The graphics options of the program consist of:
- Appearance
- Draw
- Scaling
- Window
- Misc
Miscellaneous
Changes in version 1.6:
Changes in version 1.5 (until 20110331) :
- One graphics window
- Removed commands: close, plotfunc, rduration, rtest, xyfile
- Added commands:
colourbw,
colourtype,
dotsizemax,
dotsizetype,
euler.
evpylog.
fractionalsize,
gridsize,
histoautoscale.
histobar.
histobarfraction.
histofile.
histoline.
histoprint.
hklfilename.
limitfilter.
minmax.
msapower.
overwrite,
pkrestfrac,
plot2dfixscale.
plotdelsig,
plotgrid,
plotray,
plotshell,
readsca,
readfcf,
readhk6,
readicalc.
readraw.
rvar.
shcolumn,
speedcorrect.
step.
stepfull.
xyicr.
xyplot.
window.
write.
wulffvec.
- Renamed commands:
Changes in version 1.4 (until 2006062700):
- Added variables: equiv2, gonio1, gonio2, gonio3, shiftframes, tau
(see reflection variables)
- Added commands:
colourvar,
colourscheme,
drx,
drxfile,
expression,
extract,
extractsubset,
funclimit,
plotfunc
- Removed variable gravrot
- Added option inside/outside to
region
Changes in version 1.3 (until 2004070900):
- Renamed commands:
- Changed commands:
icr plots as a function of one of the
variables
- New commands: field,
plot
- Removed commands: plotshift
EVPY Suite Overview