some more files I think we may need

models/galaxy_model/mbulge.dat

+9.00000 2.64842 3.57020 5.82344 6.54320 6.59600 6.61760 6.35000
+9.21429 3.27630 4.09011 6.07943 6.75774 6.77943 6.78830 6.58500
+9.42857 3.88126 4.58514 6.30572 6.95914 6.96857 6.97243 6.81000
+9.64286 4.46314 5.05933 6.51668 7.13031 7.17086 7.18744 7.03500
+9.85714 5.02174 5.51672 6.72667 7.35497 7.38986 7.40413 7.26000
+10.07143 5.55645 5.95687 6.93568 7.60431 7.61657 7.62159 7.46643
+10.28571 6.06663 6.37932 7.14368 7.83920 7.84957 7.85381 7.63571
+10.50000 6.55165 6.78361 7.35064 8.05700 8.09000 8.10350 7.80500
+10.71429 7.01088 7.16929 7.55653 8.27480 8.33043 8.35319 7.97429
+10.92857 7.44369 7.53591 7.76134 8.49260 8.57086 8.60287 8.14357
+11.14286 7.84945 7.88301 7.96503 8.71040 8.81129 8.85256 8.36000
+11.35714 8.16693 8.18434 8.22689 8.92820 9.05171 9.10224 8.60000
+11.57143 8.36672 8.42721 8.57507 9.14600 9.29214 9.35193 8.84000
+11.78571 8.56297 8.63471 8.81009 9.36380 9.53257 9.60161 9.08000
+12.00000 8.75720 8.82200 8.98040 9.58160 9.77300 9.85130 9.32000

models/galaxy_model/numdenrange0.txt

+1.00000e+09 2.72784e-03 3.12528e-03 3.60119e-03 5.88295e-03 6.42891e-03 7.06217e-03 4.78013e-03
+1.63789e+09 2.66275e-03 2.78114e-03 3.22697e-03 5.28790e-03 5.50619e-03 5.93243e-03 4.17255e-03
+2.68270e+09 2.46587e-03 2.58782e-03 2.87409e-03 4.58175e-03 4.80067e-03 4.98177e-03 3.53024e-03
+4.39397e+09 2.17857e-03 2.49682e-03 2.63025e-03 4.06094e-03 4.18183e-03 4.29177e-03 3.06332e-03
+7.19686e+09 1.89428e-03 2.18615e-03 2.37261e-03 3.45610e-03 3.65100e-03 3.83022e-03 2.69192e-03
+1.17877e+10 1.62196e-03 1.80502e-03 1.98090e-03 2.97304e-03 3.21960e-03 3.38692e-03 2.40791e-03
+1.93070e+10 1.36267e-03 1.42480e-03 1.64031e-03 2.51753e-03 2.79752e-03 2.94112e-03 2.06948e-03
+3.16228e+10 1.09572e-03 1.11537e-03 1.32336e-03 2.09491e-03 2.37264e-03 2.47930e-03 1.77874e-03
+5.17947e+10 8.07010e-04 8.74346e-04 1.03476e-03 1.69400e-03 1.93412e-03 1.99085e-03 1.37425e-03
+8.48343e+10 5.54084e-04 6.38925e-04 7.68744e-04 1.29210e-03 1.47641e-03 1.47771e-03 1.00789e-03
+1.38950e+11 3.39012e-04 4.16022e-04 4.80812e-04 8.69654e-04 9.62027e-04 1.01615e-03 6.61409e-04
+2.27585e+11 1.57343e-04 1.64916e-04 2.70334e-04 4.76540e-04 5.97957e-04 6.05548e-04 3.51921e-04
+3.72759e+11 2.92455e-05 5.75171e-05 1.04622e-04 1.92162e-04 2.70055e-04 3.16217e-04 1.41641e-04
+6.10540e+11 2.48224e-06 1.22809e-05 2.38279e-05 5.13978e-05 7.93777e-05 1.19923e-04 3.54624e-05
+1.00000e+12 1.27134e-07 1.16775e-06 2.19309e-06 9.10119e-06 1.22773e-05 2.71632e-05 3.60901e-06

models/galaxy_model/numdenrange1.txt

+1.00000e+09 4.77805e-04 1.05205e-03 1.80093e-03 3.77046e-03 4.34493e-03 5.59529e-03 2.91677e-03
+1.63789e+09 4.73215e-04 9.87925e-04 1.51395e-03 3.02984e-03 3.64515e-03 4.39292e-03 2.38134e-03
+2.68270e+09 4.56596e-04 8.91410e-04 1.26959e-03 2.65145e-03 3.05346e-03 3.44529e-03 1.87957e-03
+4.39397e+09 4.27272e-04 7.41082e-04 1.06040e-03 2.31275e-03 2.55322e-03 2.69669e-03 1.62685e-03
+7.19686e+09 3.97688e-04 5.95121e-04 8.55026e-04 2.00656e-03 2.14037e-03 2.22327e-03 1.37894e-03
+1.17877e+10 3.65141e-04 4.71354e-04 5.94477e-04 1.69619e-03 1.77239e-03 2.00177e-03 1.12597e-03
+1.93070e+10 3.26433e-04 3.55640e-04 4.30927e-04 1.40168e-03 1.49954e-03 1.76589e-03 9.08367e-04
+3.16228e+10 2.53064e-04 2.71045e-04 3.06327e-04 1.13430e-03 1.22279e-03 1.48777e-03 7.42820e-04
+5.17947e+10 1.73385e-04 1.74853e-04 2.25118e-04 8.86984e-04 9.52344e-04 1.16832e-03 5.74524e-04
+8.48343e+10 1.04288e-04 1.09690e-04 1.55932e-04 6.55706e-04 7.12802e-04 8.17659e-04 3.98385e-04
+1.38950e+11 5.57243e-05 5.98893e-05 8.57232e-05 4.42128e-04 4.83559e-04 5.06747e-04 2.43739e-04
+2.27585e+11 2.29472e-05 2.66783e-05 4.13646e-05 2.23225e-04 2.76841e-04 3.26748e-04 1.08278e-04
+3.72759e+11 5.79837e-06 6.96326e-06 1.25792e-05 9.76879e-05 1.30357e-04 1.71773e-04 4.08787e-05
+6.10540e+11 4.90687e-07 5.93903e-07 1.68313e-06 2.79966e-05 4.35177e-05 6.41435e-05 9.43911e-06
+1.00000e+12 4.97628e-09 9.47345e-09 1.03985e-07 4.30067e-06 8.21235e-06 1.36414e-05 1.39060e-06

models/galaxy_model/prior_check.py

+import numpy
+from numpy import inf
+
+def fpairtest(delta):
+	return 1./(delta+1.)*(1.-0.25**(delta+1.))
+
+def schechterf(M,M0,phi0,alpha0):
+	return phi0*numpy.power(M/M0,1.+alpha0)*numpy.exp(-M/M0)
+
+def Mbh(Mstar,alpha,beta):
+	return alpha+beta*numpy.log10(Mstar/1.e11)
+
+
+def check_p(livepoint):
+
+	data = numpy.genfromtxt('mbulge.dat')
+	data0 = numpy.genfromtxt('numdenrange0.txt')
+	data1 = numpy.genfromtxt('numdenrange1.txt')
+	m0 = numpy.logspace(9,12,15)
+	z0 = [0.4,0.6,0.8]
+	z1 = [1.5,2.,2.5]
+
+	l = 0.0
+	for i in z0:
+		phi0 = livepoint['Phi0']+livepoint['PhiI']*i
+		a0 = livepoint['alpha']+livepoint['alphaI']*i
+		model = schechterf(m0,10.**livepoint['M0'],10.**phi0,a0)
+		if (model<=data0[:,6]).all()==True and (model>=data0[:,1]).all()==True:
+			for j in z1:
+				phi0 = livepoint['Phi0']+livepoint['PhiI']*j
+				a0 = livepoint['alpha']+livepoint['alphaI']*j
+				model = schechterf(m0,10.**livepoint['M0'],10.**phi0,a0)
+				if (model<=data1[:,6]).all()==True and (model>=data1[:,1]).all()==True:
+					l += 0.0
+				else:
+					return -inf
+		else:
+			return -inf
+
+	model = Mbh(m0,livepoint['Ms'],livepoint['theta'])
+	if (model<=data[:,6]).all()==True and (model>=data[:,1]).all()==True:
+		return l
+	else:
+		return -inf
+
+