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
