Hello Santhos-- > > While I was going through example scripts, in eginput/gb3_ensemble and > comparing against eginput/gb3_rdc [single structure refinements], I found > some thing that looks weird. > > In eginput/gb3_ensemble: > dyn = ivm.IVM(esim) #manipulator for T-A dynamics > dynx = ivm.IVM(esim) #manipulator for Cartesian dynamics > > but dynx was never used as a cartesian dynamics [no call to > protocol.cartesianTopology(dynx).] > Is there a reason to do that? How to perform a cartesian minimization when > running an ensemble simulation?
These scripts were written before those helper functions existed. They are manually setup for T-A and Cartesian dynamics. > > I have tried to mention, > protocol.torsionTopology(dyn) > > protocol.cartesianTopology(dynx) > This must be done carefully, as the IVM objects are already configured... > but the resulting refinements should some bigger R-factors than ones > that did not. ...as you found out. In particular, the RDC setup has changed to use the freedom member instead of manually calling the configIVM functions. To use cartesian/torsionTopology you will need lines like rdc.oTensor.setFreedom("fixDa, fixRh") instead of the calls to configIVM*. The contents of the string passed to setFreedom can be found here: http://nmr.cit.nih.gov/xplor-nih/doc/current/python/ref/varTensorTools.html best regards-- Charles _______________________________________________ Xplor-nih mailing list Xplor-nih@cake.cit.nih.gov http://cake.cit.nih.gov/mailman/listinfo/xplor-nih