blob: c3d4030fa02909e0884cc6bcb4afe8ffbc2ad0a9 (
plain)
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
|
#!/usr/bin/env python
import sys
import math
import numpy
lx1_array=[]
lx2_array=[]
lx3_array=[]
for i in open('summary_lx.dat'):
x1,x2,x3=i.split()
x1=float(x1)
x2=float(x2)
x3=float(x3)
lx1_array.append(x1)
lx2_array.append(x2)
lx3_array.append(x3)
lx1_array=numpy.array(lx1_array)
lx2_array=numpy.array(lx2_array)
lx3_array=numpy.array(lx3_array)
f=open('lx_result.txt','w')
f.write("Lx(bolo)= %4.2E +/- %4.2E erg/s\n"%(lx1_array[0],lx1_array.std()))
print("Lx(bolo)= %4.2E +/- %4.2E erg/s"%(lx1_array[0],lx1_array.std()))
f.write("Lx(0.7-7)= %4.2E +/- %4.2E erg/s\n"%(lx2_array[0],lx2_array.std()))
print("Lx(0.7-7)= %4.2E +/- %4.2E erg/s"%(lx2_array[0],lx2_array.std()))
f.write("Lx(0.1-2.4)= %4.2E +/- %4.2E erg/s\n"%(lx3_array[0],lx3_array.std()))
print("Lx(0.1-2.4)= %4.2E +/- %4.2E erg/s"%(lx3_array[0],lx3_array.std()))
|
