posDiffPotTools
index
/home/schwitrs/xplor/python/posDiffPotTools.py


 
Tools to help create/analyze potential terms based on atomic coordinate
difference

 
Functions
       
analyze(potList)
perform analysis of PosDiffPot terms and return nicely formatted
summary
create_PosDiffPot(name, selection, selection2=None, pdbFile=None, psf=None, absoluteFit=False, removeMissing=False, cmpSel=None)
Create a PosDiffPot term based on the specified atom selection.
 
The name of the PosDiffPot term is given by the name argument (string).
If selection2 is not specified, it is taken to be the same as selection.
If pdbFile is specified, the comparison is made with the coordinates in that
file (using selection2).  It is an error to omit both the selection2 and
pdbFile arguments.  For this potential term, the order of atom indices in
selection2 must match that in selection. If selection is an ordered
atomSel.AtomSel object, then selection2 will also be ordered if omitted
or specified by string. If selection2 is the string "initial", the
selection string is taken from selection, while the comparison coordinates
are those loaded when create_PosDiffPot is called. 
 
This potential term is based on the posSymmPot.PosSymmPot potential.
In this case each atom is placed in its own equivalent set (as in the
NCS potential).
 
If absoluteFit=True, selection and selection2 will be physically overlaid,
as rigid body translation and rotation will not be performed before the
fit.
 
cmpSel is a selection or pair of selections (the later is required if
selection2 is supplied) to use for an additional rmsd computation.  It is
not used in an energy calculation.
 
If removeMissing=True, remove from a selection atoms not present in
the other, for cases where the selections (selection, selection2), 
(cmpSel, cmpSel2) do not have the same number of atoms. This will
fail if more than one segid (chain) is selected.
                              
If pdbFile (a string with a PDB filename) is specified, by default the
associated coordinates are loaded using protocol.loadPDB (with
deleteUnknownAtoms=True); this assumes a standard topology.  A non-standard
topology for the input PDB file can be accomodated with the psf argument,
which can be a string with the contents of a PSF file, a string with a PSF
filename, or a sequence with one of more PSF filenames.  If psf is
specified, the coordinates are loaded with protocol.initCoords, using
deleteUnknownAtoms=True.
 
The returned potential term has four members beyond those in a
 posSymmPot.PosSymmPot term:
 
    pdbFile     - the value of the pdbFile argument.
    selection   - an atomSel.AtomSel corresponding to the selection
                  argument.
    selection2  - an atomSel.AtomSel corresponding to the selection2
                  argument.
    xsim        - a xplorSimulation.XplorSimulation generated for
                  containing atomic information if pdbFile is
                  specified.
findMatchingSubSel(sel1, sel2, name1, name2)
Return selections corresponding to the inputs which 
have the same length and matching segid, resid, atom name.