Skip to content
Snippets Groups Projects
Commit fd649dd5 authored by Yngve Levinsen's avatar Yngve Levinsen
Browse files

fixes in TraceWin.plt

calc_avg and calc_sigma had some errors, apologies
parent 2d994a8d
No related branches found
No related tags found
No related merge requests found
......@@ -237,16 +237,13 @@ class plt:
self.avg=dict(x=[], xp=[], y=[], yp=[], E=[], phi=[])
vals=self._columns
vals=self._columns[:-1]
for i in self.Nelp:
data=self[i]
for v in vals:
self.avg[v].append(numpy.average(data[v]))
for v in vals:
self.avg[v]=numpy.array(self.avg[v])
def calc_minmax(self,pmin=5,pmax=95):
'''
Calculates min/max values of beam coordinates
......@@ -291,10 +288,15 @@ class plt:
if not hasattr(self,'avg'):
self.calc_avg()
vals=self._columns
vals=self._columns[:-1]
self.sigma=[]
for j in xrange(len(self.Nelp)):
i=self.Nelp[j]
data=self[i]
self.sigma.append([[numpy.mean( (data[n]-self.avg[m][j]) * (data[m] - self.avg[m][j]) ) for n in vals] for m in vals])
self.sigma=numpy.array([[numpy.mean( (self[n]-self.avg[n]) * (self[m] - self.avg[m]) ) for n in vals] for m in vals])
print self.sigma.shape
self.sigma=numpy.array(self.sigma)
def calc_twiss(self):
'''
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment