def __str__(self):
ret = "SectorMapMatrix:\n"
ret += util.prettyPrint66(self.RE)
return ret
def __str__(self):
ret = "RFCavity_Matrix: "
ret += "voltage = %10g[V], wavelength = %10g[m], phase = %10g[rad], E0 = %10g[GeV], m0 = %10g[MeV/c^2]" % (self.voltage,self.wavelength,self.phase, self.E0/1e9, self.m0/1e6)
ret += util.prettyPrint66(self.RE)
return ret
ifile.close()
#Sanity check input
if myRing==None:
print "ERROR: No ring was set up"
exit(1)
if myBeam==None:
print "No beam was set up"
exit(1)
print
print "Ring:"
print myRing
print "Total matrix:"
totalRE = myRing.getTotalMatrix()
print util.prettyPrint66(totalRE)
print "Eigenvalues and vectors:"
(w, v) = np.linalg.eig(totalRE)
for i in xrange(6):
print "#%i : abs= %16.10g, angle= %16.10g [2pi] -> "%(i,np.absolute(w[i]), np.angle(w[i])/(2*np.pi)), w[i]
#print "\t", v[i] #Complex vector, needs some better formatting...
print
print "Beam:"
print myBeam
print
#Track!
print "Tracking!"
print "# particles = ", myBeam.getNumParticles()
print "# elements = ", myRing.getNumElements()
