Example #1
0
 def visualize(self):
     if self.ExperimentalSAXS=="":
         fnInt=self.workingDirPath("pseudoatoms00.int")
     else:
         fnInt=self.workingDirPath("pseudoatoms00.fit")
     
     import numpy
     x=numpy.loadtxt(fnInt,skiprows=1)
     xplotter = XmippPlotter(*[1,1],windowTitle="SAXS Curves")
     a = xplotter.createSubPlot('SAXS curves', 'Armstrongs^-1', 'log(SAXS)', yformat=False)
     a.plot(x[:,0], numpy.log(x[:,1]))
     a.plot(x[:,0], numpy.log(x[:,2]))
     if self.ExperimentalSAXS=="":
         xplotter.showLegend(['SAXS in solution','SAXS in vacuo'])
     else:
         xplotter.showLegend(['Experimental SAXS','SAXS from volume'])
     xplotter.draw()
     xplotter.show()