def testH():
a=ELoss()
ke=40e3
E=ke*eV/gE; Mp=131.9*u/gM
v=math.sqrt(2*E/Mp)
Zp=1
n=4.518696e25
print n,a.calSP(Zp,v,n)*gE/gL/(eV/A)
def testH():
a = ELoss()
ke = 40e3
E = ke * eV / gE
Mp = 131.9 * u / gM
v = math.sqrt(2 * E / Mp)
Zp = 1
n = 4.518696e25
print n, a.calSP(Zp, v, n) * gE / gL / (eV / A)
def testH():
a=ELoss()
ke=100e3
E=ke*eV/gE; Mp=1*u/gM
v=math.sqrt(2*E/Mp)
Zp=1
n=1e22
while n<1e32:
print n,a.calSP(Zp,v,n)*gE/gL/(1000*eV/1e-2)
n=n*10
def testH():
a=ELoss()
ke=1000e3
E=ke*eV/gE; Mp=1*u/gM
v=math.sqrt(2*E/Mp)
Zp=1
n=1e22
while n<1e32:
print n,a.calSP(Zp,v,n)/n*gE*gL*gL/(eV*1e-24)
n=n*10
def testCTi():
a=ELoss()
a.readObf("ti.orb")
ke=1.0
while ke<=10000:
E=ke*1e3*eV/gE; Mp=12*u/gM
v=math.sqrt(2*E/Mp)
Mt=47.9*u/gM
Nt=4.5189/Mt
Zp=6
print ke*1e3,a.calLDA_SP(Zp,v,Nt)*gE/gL/eV*A
ke=ke*2
def testC4C():
a=ELoss()
a.readObf("c.orb")
E=100e3*12*eV/gE; Mp=12*u/gM
v=math.sqrt(2*E/Mp)
Zp=4
dz=0.5*A/gL
dt=dz/v
z0=30*A/gL;
d0=0.1
while d0<20.0:
print d0, a.calELoss(Zp, v, dt, z0, d0*A/gL)*gE/eV
d0=d0 + 0.1
def testHAl():
a=ELoss()
a.readObf("al.orb_fit")
ke=10000e3
E=ke*eV/gE; Mp=1.0*u/gM
v=math.sqrt(2*E/Mp)
Zp=1.0
i=0
while i<len(a.density):
d=float(a.density[i][0])
n=float(a.density[i][1])
d0=d*gL/A
nm=A*10.0
sp=a.calSP(Zp,v,n)
L=a.calLFun(Zp,v,n)
#print d0,n*4*PI*d**2,L,sp,sp*4*PI*d*d*gE*gL/eV/A
print d0,L,L*n*4*PI*d*d/gL*A
i=i+1
def testHAl():
a = ELoss()
a.readObf("al.orb_fit")
ke = 10000e3
E = ke * eV / gE
Mp = 1.0 * u / gM
v = math.sqrt(2 * E / Mp)
Zp = 1.0
i = 0
while i < len(a.density):
d = float(a.density[i][0])
n = float(a.density[i][1])
d0 = d * gL / A
nm = A * 10.0
sp = a.calSP(Zp, v, n)
L = a.calLFun(Zp, v, n)
#print d0,n*4*PI*d**2,L,sp,sp*4*PI*d*d*gE*gL/eV/A
print d0, L, L * n * 4 * PI * d * d / gL * A
i = i + 1