plot.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python


	from numpy import *
	import sys
	import coolinglib


	##############################################################################
	# init object
	##############################################################################

	ct = coolinglib.CTable("HM01/IonRates_Hydrogen.HM01.hdf5")


	ion = 0


	# recombinaison rate

	Ts = 10**arange(1,9,0.1)
	As = zeros(len(Ts))
	Bs = zeros(len(Ts))
	Cools = zeros(len(Ts))

	for i,T in enumerate(Ts):
	As[i] = ct.RecombinasionRateCoefficient(ion,T)

	for i,T in enumerate(Ts):
	Bs[i] = ct.CollisionalIonisationRateCoefficient(ion,T)

	for i,T in enumerate(Ts):
	Cools[i] = ct.CoolingEfficiency(ion,T)


	# photo heating efficiency

	zs = arange(0,20,0.01)
	Hs = zeros(len(zs))
	Cs = zeros(len(zs))

	for i,z in enumerate(zs):
	Hs[i] = ct.PhotoHeatingEfficiency(ion,z)


	for i,z in enumerate(zs):
	Cs[i] = ct.PhotoIonisationRate(ion,z)



	##############################################################################
	# some plots
	##############################################################################
	import Ptools as pt



	#pt.scatter(ct.Temperatures,ct.Alpha[ion])
	#pt.plot(Ts,As)

	#pt.scatter(ct.Temperatures,ct.BetaColl[ion])
	#pt.plot(Ts,Bs)

	#pt.scatter(ct.Temperatures,ct.Cool[ion])
	#pt.plot(Ts,Cools)

	#pt.scatter(ct.Redshifts,ct.Heat[ion])
	#pt.plot(zs,Hs)

	pt.scatter(ct.Redshifts,ct.Heat[ion])
	pt.plot(zs,Hs)

	pt.show()