From ffd178e0bd72562a3c2cff9747b6e656edc881dc Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Fri, 27 May 2016 22:47:24 +0800 Subject: Add mass_profile tools * These tools are mainly use to calculate the total gravitational mass profile, as well as the intermediate products (e.g., surface brightness profile fitting, gas density profile, NFW fitting, etc.) * There are additional tools for calculating the luminosity and flux. * These tools mainly developed by Junhua GU, and contributed by Weitian (Aaron) LI, and Zhenghao ZHU. --- mass_profile/analyze_lx.py | 32 ++++++++++++++++++++++++++++++++ 1 file changed, 32 insertions(+) create mode 100755 mass_profile/analyze_lx.py (limited to 'mass_profile/analyze_lx.py') diff --git a/mass_profile/analyze_lx.py b/mass_profile/analyze_lx.py new file mode 100755 index 0000000..c3d4030 --- /dev/null +++ b/mass_profile/analyze_lx.py @@ -0,0 +1,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())) + -- cgit v1.2.2