Exemplo n.º 1
0
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)
Exemplo n.º 2
0
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)
Exemplo n.º 3
0
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
Exemplo n.º 4
0
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
Exemplo n.º 5
0
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
Exemplo n.º 6
0
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
Exemplo n.º 7
0
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
Exemplo n.º 8
0
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