-
Notifications
You must be signed in to change notification settings - Fork 6
/
plot-au.py
executable file
·42 lines (28 loc) · 997 Bytes
/
plot-au.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
#!/usr/bin/python
# -*- coding: utf-8 -*-
from __future__ import unicode_literals
import numpy as np
import pylab as plt
from getData import *
import matplotlib
matplotlib.rcParams.update({'font.size': 30})
range=(1,280)
lambdas=np.linspace(range[0],range[1],100)
ag=getData("database/main/Au/Ordal.yml",lambdas);
fig=plt.figure()
plt.plot(lambdas,[ -(x*x).real for x in ag ],'r-',label=r'-Re($\varepsilon_{Au}$)',lw=4)
plt.plot(lambdas,[ (x*x).imag for x in ag] ,'b-',label=r'Im($\varepsilon_{Au}$)',lw=4)
#plt.plot(lambdasTiO2*1e3,[ (x*x).real for x in tio2],'g-',label=r'Real($\varepsilon_{TiO_{2}}$)')
c=299792458
lambdas=lambdas*10e-6
freq=c/lambdas
freq=freq/10e12
#plt.plot(freq,[ (x*x).real for x in ag ],'r-',label=r'Real($\varepsilon_{Ag}$)')
plt.legend(loc=4)
#plt.title(r"Współczynnik przenikalnośći elektrycznej ($\varepsilon$)");
plt.xlabel('wavelength [um]')
plt.xlim([1,280])
plt.tight_layout()
plt.yscale('log')
plt.savefig("../phd/images/aueps.png")
plt.show()