smf_literature.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.table import *\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from tqdm import tqdm\n",
    "from astropy.table import hstack\n",
    "import matplotlib\n",
    "import matplotlib.backends.backend_pdf"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "tags": []
   },
   "source": [
    "#### Ilbert et al. 2013"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5\n",
    "smf_il_1 = np.array([\n",
    "8.090647482014388, -1.603225806451613,\n",
    "8.350359712230215, -1.618279569892473,\n",
    "8.616906474820142, -1.7462365591397848,\n",
    "8.890287769784171, -1.8440860215053763,\n",
    "9.170503597122302, -1.972043010752688,\n",
    "9.45071942446043, -2.0774193548387094,\n",
    "9.73093525179856, -2.2053763440860212,\n",
    "10.01115107913669, -2.2806451612903222,\n",
    "10.277697841726617, -2.325806451612903,\n",
    "10.544244604316546, -2.408602150537634,\n",
    "10.831294964028777, -2.5892473118279566,\n",
    "11.077338129496402, -2.958064516129032,\n",
    "11.357553956834533, -3.4774193548387093,\n",
    "11.61726618705036, -4.245161290322581,\n",
    "12.068345323741006, -5.223655913978494])\n",
    "smf_il_1 = smf_il_1.reshape(len(smf_il_1)//2,2)\n",
    "np.save('smf_lit/smf_il_0.2_0.5',smf_il_1)\n",
    "\n",
    "# 0.5-0.8\n",
    "smf_il_2 = np.array([\n",
    "8.836756756756756, -1.9409660107334525,\n",
    "9.06954954954955, -1.9935599284436494,\n",
    "9.336576576576576, -2.098747763864043,\n",
    "9.576216216216217, -2.233989266547406,\n",
    "9.843243243243244, -2.3767441860465115,\n",
    "10.11027027027027, -2.444364937388193,\n",
    "10.349909909909911, -2.481932021466905,\n",
    "10.61009009009009, -2.572093023255814,\n",
    "10.856576576576575, -2.729874776386404,\n",
    "11.103063063063063, -2.992844364937388,\n",
    "11.315315315315315, -3.4586762075134168,\n",
    "11.568648648648647, -3.9695885509838997,\n",
    "11.904144144144144, -5.1642218246869405])\n",
    "smf_il_2 = smf_il_2.reshape(len(smf_il_2)//2,2)\n",
    "np.save('smf_lit/smf_il_0.5_0.8',smf_il_2)\n",
    "\n",
    "# 0.8-1.1\n",
    "smf_il_3 = np.array([\n",
    "9.247567567567566, -2.1363148479427547,\n",
    "9.494054054054054, -2.211449016100179,\n",
    "9.754234234234234, -2.331663685152057,\n",
    "9.993873873873873, -2.444364937388193,\n",
    "10.247207207207206, -2.5044722719141324,\n",
    "10.500540540540541, -2.534525939177102,\n",
    "10.733333333333333, -2.609660107334526,\n",
    "10.993513513513513, -2.7974955277280857,\n",
    "11.226306306306304, -3.128085867620751,\n",
    "11.47963963963964, -3.744186046511628,\n",
    "11.685045045045044, -4.698389982110912,\n",
    "12.034234234234233, -5.344543828264758])\n",
    "smf_il_3 = smf_il_3.reshape(len(smf_il_3)//2,2)\n",
    "np.save('smf_lit/smf_il_0.8_1.1',smf_il_3)\n",
    "\n",
    "# 1.1-1.5\n",
    "smf_il_4 = np.array([\n",
    "9.541981981981982, -2.286583184257603,\n",
    "9.788468468468468, -2.376744186046512,\n",
    "10.028108108108107, -2.4894454382826483,\n",
    "10.274594594594594, -2.609660107334526,\n",
    "10.541621621621621, -2.677280858676208,\n",
    "10.794954954954953, -2.7899821109123444,\n",
    "11.034594594594594, -3.037924865831843,\n",
    "11.267387387387387, -3.496243291592129,\n",
    "11.507027027027027, -4.360286225402505,\n",
    "11.732972972972973, -5.179248658318427])\n",
    "smf_il_4 = smf_il_4.reshape(len(smf_il_4)//2,2)\n",
    "np.save('smf_lit/smf_il_1.1_1.5',smf_il_4)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "from scipy.spatial.distance import cdist"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.7781842324894082\n"
     ]
    }
   ],
   "source": [
    "x_interp = np.linspace(smf_il_4[0][0],smf_il_4[-1][0],100)\n",
    "y_interp = np.interp(x_interp,smf_il_4[:,0],smf_il_4[:,1])+1\n",
    "\n",
    "smf_il_4_interp = np.transpose(np.stack([x_interp,y_interp]))\n",
    "arr = cdist(smf_il_4,smf_il_4_interp)\n",
    "print(np.mean(np.min(arr,axis=0)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "tags": []
   },
   "source": [
    "#### Muzzin et al. 2013"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "\n",
    "# 0.2-0.5\n",
    "smf_mu_1 = np.array([\n",
    "8.440740740740742, -1.5781818181818184,\n",
    "8.526748971193417, -1.6872727272727275,\n",
    "8.620576131687244, -1.7090909090909092,\n",
    "8.722222222222223, -1.741818181818182,\n",
    "8.823868312757202, -1.8181818181818183,\n",
    "8.91769547325103, -1.850909090909091,\n",
    "9.003703703703705, -1.8836363636363638,\n",
    "9.105349794238684, -1.9272727272727275,\n",
    "9.206995884773663, -1.981818181818182,\n",
    "9.316460905349794, -2.0254545454545454,\n",
    "9.425925925925927, -2.0909090909090913,\n",
    "9.519753086419753, -2.101818181818182,\n",
    "9.605761316872428, -2.1563636363636363,\n",
    "9.707407407407407, -2.221818181818182,\n",
    "9.809053497942386, -2.1890909090909094,\n",
    "9.910699588477367, -2.210909090909091,\n",
    "10.004526748971193, -2.2872727272727276,\n",
    "10.09835390946502, -2.3090909090909095,\n",
    "10.207818930041153, -2.3090909090909095,\n",
    "10.30164609053498, -2.352727272727273,\n",
    "10.40329218106996, -2.330909090909091,\n",
    "10.504938271604939, -2.3854545454545457,\n",
    "10.606584362139918, -2.3854545454545457,\n",
    "10.700411522633745, -2.4400000000000004,\n",
    "10.809876543209876, -2.5381818181818185,\n",
    "10.903703703703703, -2.6145454545454547,\n",
    "10.99753086419753, -2.690909090909091,\n",
    "11.09917695473251, -2.832727272727273,\n",
    "11.20082304526749, -3.04,\n",
    "11.310288065843622, -3.061818181818182,\n",
    "11.40411522633745, -3.2581818181818187,\n",
    "11.497942386831276, -3.5963636363636367,\n",
    "11.599588477366254, -4.207272727272727,\n",
    "11.693415637860081, -4.218181818181819,\n",
    "11.802880658436214, -4.512727272727274,\n",
    "11.89670781893004, -4.981818181818182,\n",
    "11.990534979423867, -4.698181818181819])\n",
    "smf_mu_1 = smf_mu_1.reshape(len(smf_mu_1)//2,2)\n",
    "np.save('smf_lit/smf_mu_0.2_0.5',smf_mu_1)\n",
    "\n",
    "# 0.5-1.0\n",
    "smf_mu_2 = np.array([\n",
    "8.964609053497943, -2.0897959183673462,\n",
    "9.066255144032922, -2.0680272108843534,\n",
    "9.16008230452675, -2.1006802721088427,\n",
    "9.26954732510288, -2.1442176870748293,\n",
    "9.363374485596708, -2.1442176870748293,\n",
    "9.457201646090535, -2.1442176870748293,\n",
    "9.558847736625514, -2.209523809523809,\n",
    "9.64485596707819, -2.2748299319727883,\n",
    "9.746502057613169, -2.2421768707482985,\n",
    "9.848148148148148, -2.263945578231292,\n",
    "9.934156378600823, -2.318367346938775,\n",
    "10.043621399176955, -2.3510204081632646,\n",
    "10.137448559670782, -2.4163265306122446,\n",
    "10.246913580246915, -2.372789115646258,\n",
    "10.332921810699588, -2.394557823129251,\n",
    "10.43456790123457, -2.427210884353741,\n",
    "10.528395061728396, -2.514285714285714,\n",
    "10.622222222222224, -2.503401360544217,\n",
    "10.723868312757201, -2.5469387755102035,\n",
    "10.817695473251028, -2.6340136054421763,\n",
    "10.92716049382716, -2.742857142857142,\n",
    "11.020987654320987, -2.906122448979591,\n",
    "11.130452674897118, -3.004081632653061,\n",
    "11.216460905349795, -3.2108843537414957,\n",
    "11.318106995884772, -3.428571428571428,\n",
    "11.40411522633745, -3.7006802721088428,\n",
    "11.505761316872427, -3.972789115646258,\n",
    "11.599588477366254, -4.179591836734693,\n",
    "11.709053497942387, -4.582312925170068])\n",
    "smf_mu_2 = smf_mu_2.reshape(len(smf_mu_2)//2,2)\n",
    "np.save('smf_lit/smf_mu_0.5_1.0',smf_mu_2)\n",
    "\n",
    "# 1.0-1.5\n",
    "smf_mu_3 = np.array([\n",
    "9.527572016460907, -2.394557823129251,\n",
    "9.629218106995886, -2.3836734693877544,\n",
    "9.730864197530865, -2.394557823129251,\n",
    "9.81687242798354, -2.4598639455782303,\n",
    "9.91851851851852, -2.5251700680272102,\n",
    "10.020164609053499, -2.536054421768707,\n",
    "10.129629629629632, -2.6340136054421763,\n",
    "10.223456790123457, -2.6013605442176866,\n",
    "10.325102880658438, -2.6231292517006795,\n",
    "10.411111111111111, -2.6448979591836728,\n",
    "10.512757201646092, -2.699319727891156,\n",
    "10.606584362139918, -2.7102040816326527,\n",
    "10.708230452674899, -2.797278911564625,\n",
    "10.802057613168724, -2.851700680272108,\n",
    "10.911522633744857, -2.960544217687074,\n",
    "11.005349794238683, -3.0693877551020403,\n",
    "11.106995884773664, -3.2544217687074823,\n",
    "11.208641975308645, -3.4612244897959177,\n",
    "11.30246913580247, -3.831292517006802,\n",
    "11.396296296296295, -3.9945578231292513,\n",
    "11.505761316872428, -4.429931972789115,\n",
    "11.599588477366257, -4.636734693877551,\n",
    "11.693415637860083, -5.191836734693878,\n",
    "11.795061728395064, -5.877551020408163])\n",
    "smf_mu_3 = smf_mu_3.reshape(len(smf_mu_3)//2,2)\n",
    "np.save('smf_lit/smf_mu_1.0_1.5',smf_mu_3)\n",
    "\n",
    "# 1.5-2.0\n",
    "smf_mu_4 = np.array([\n",
    "10.172689938398358, -2.7102040816326527,\n",
    "10.266324435318275, -2.8190476190476184,\n",
    "10.375564681724846, -2.8408163265306117,\n",
    "10.469199178644763, -2.8952380952380947,\n",
    "10.555030800821354, -2.9496598639455778,\n",
    "10.656468172484598, -2.9496598639455778,\n",
    "10.757905544147844, -2.9823129251700675,\n",
    "10.843737166324434, -3.1455782312925162,\n",
    "10.952977412731006, -3.1564625850340127,\n",
    "11.038809034907596, -3.297959183673469,\n",
    "11.140246406570842, -3.570068027210884,\n",
    "11.233880903490759, -3.7224489795918365,\n",
    "11.327515400410677, -4.070748299319727,\n",
    "11.42895277207392, -4.4190476190476184,\n",
    "11.52258726899384, -4.70204081632653,\n",
    "11.616221765913757, -5.474829931972788])\n",
    "smf_mu_4 = smf_mu_4.reshape(len(smf_mu_4)//2,2)\n",
    "np.save('smf_lit/smf_mu_1.5_2.0',smf_mu_4)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "tags": []
   },
   "source": [
    "#### Muzzin et al. 2013 (Q)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "smf_mu_q_1 = np.array([\n",
    "8.430878422357326, -2.4451467268623017,\n",
    "8.531637792149642, -2.607674943566591,\n",
    "8.625082284295125, -2.4884875846501124,\n",
    "8.725981987471554, -2.56433408577878,\n",
    "8.819005479464703, -2.7051918735891642,\n",
    "8.927763852151395, -2.7268623024830694,\n",
    "9.021068010912769, -2.694356659142212,\n",
    "9.121844922378099, -2.8460496613995483,\n",
    "9.222832333919596, -2.867720090293453,\n",
    "9.31596107595083, -2.9435665914221216,\n",
    "9.416895862473284, -2.9977426636568842,\n",
    "9.510287729599728, -2.9110609480812637,\n",
    "9.619063643959434, -2.921896162528216,\n",
    "9.720103680519973, -2.9110609480812637,\n",
    "9.828914678225708, -2.900225733634311,\n",
    "9.922324087025164, -2.8027088036117376,\n",
    "10.015365120691328, -2.932731376975169,\n",
    "10.116510407289951, -2.8568848758465006,\n",
    "10.217620610542546, -2.8027088036117376,\n",
    "10.318643105430072, -2.8027088036117376,\n",
    "10.411999889210488, -2.737697516930022,\n",
    "10.513039925771027, -2.7268623024830694,\n",
    "10.70738412109294, -2.683521444695259,\n",
    "10.808283824269367, -2.7593679458239273,\n",
    "10.9015353580117, -2.7593679458239273,\n",
    "10.606309001186371, -2.7160270880361166,\n",
    "11.010223564006335, -2.8243792325056427,\n",
    "11.118876686654943, -2.9110609480812637,\n",
    "11.211830011956035, -3.095259593679458,\n",
    "11.297275501207146, -3.1169300225733627,\n",
    "11.40578829047164, -3.2902934537246047,\n",
    "11.490707529532331, -3.6370203160270878,\n",
    "11.598553735222296, -4.222121896162529,\n",
    "11.699593771782835, -4.211286681715576,\n",
    "11.800142641498983, -4.503837471783296,\n",
    "11.892604799955684, -4.9914221218961625,\n",
    "11.994100920014588, -4.6988713318284425])\n",
    "smf_mu_q_1 = smf_mu_q_1.reshape(len(smf_mu_q_1)//2,2)\n",
    "np.save('smf_lit/smf_mu_q_0.2_0.5',smf_mu_q_1)\n",
    "\n",
    "# 0.5 - 1.0\n",
    "smf_mu_q_2 = np.array([\n",
    "8.967123376383082, -3.257787810383747,\n",
    "9.060460533619509, -3.322799097065463,\n",
    "9.161813085149689, -3.246952595936795,\n",
    "9.262901972393887, -3.3336343115124154,\n",
    "9.36420179106687, -3.2902934537246056,\n",
    "9.465484032120788, -3.257787810383747,\n",
    "9.550946416015986, -3.376975169300226,\n",
    "9.652263812308034, -3.322799097065463,\n",
    "9.745882211449505, -3.214446952595937,\n",
    "9.847111719646227, -3.214446952595937,\n",
    "9.948569736890796, -3.0735891647855533,\n",
    "10.042065092698813, -3.041083521444696,\n",
    "10.13550771564963, -3.041083521444696,\n",
    "10.244770195759168, -2.889390519187359,\n",
    "10.338177663471857, -2.911060948081264,\n",
    "10.431708174518004, -2.856884875846501,\n",
    "10.532867372238464, -2.9002257336343114,\n",
    "10.63414961329238, -2.867720090293454,\n",
    "10.82878658920179, -2.889390519187359,\n",
    "10.922106168819155, -2.965237020316027,\n",
    "11.023089590348963, -3.1169300225733636,\n",
    "11.116514635680717, -3.1277652370203164,\n",
    "11.209640861488362, -3.322799097065463,\n",
    "11.318358435406875, -3.506997742663657,\n",
    "11.419201235984161, -3.7453724604966143,\n",
    "11.51223957369648, -3.994582392776524,\n",
    "11.613170262369092, -4.178781038374718,\n",
    "11.705944935795433, -4.590519187358916])\n",
    "smf_mu_q_2 = smf_mu_q_2.reshape(len(smf_mu_q_2)//2,2)\n",
    "np.save('smf_lit/smf_mu_q_0.5_1.0',smf_mu_q_2)\n",
    "\n",
    "# 1.0 - 1.5\n",
    "smf_mu_q_3 = np.array([\n",
    "9.519253310883009, -3.701572293587319,\n",
    "9.62851730725642, -3.548944634357322,\n",
    "9.729853687419336, -3.4830665044569145,\n",
    "9.823385496412728, -3.4280903502157924,\n",
    "9.932420980567208, -3.41632243351636,\n",
    "10.026005523149584, -3.3288401848451077,\n",
    "10.127201280408547, -3.3496449734250477,\n",
    "10.220820978716908, -3.2404919951337083,\n",
    "10.32217493674282, -3.1637785004232577,\n",
    "10.423441005453762, -3.1412418297630236,\n",
    "10.516937658721163, -3.1079364051419884,\n",
    "10.618238883158094, -3.063729004861668,\n",
    "10.711577335658543, -3.127941863531024,\n",
    "10.812843404369486, -3.1054051928707898,\n",
    "10.921720687757016, -3.1911555594617482,\n",
    "11.007272146950935, -3.2554350289803464,\n",
    "11.108292125579954, -3.384593465660721,\n",
    "11.209206637031008, -3.578764091201356,\n",
    "11.30989263626313, -3.9137944592725553,\n",
    "11.403143199448609, -4.032184141992128,\n",
    "11.503618264324802, -4.497238887783849,\n",
    "11.59683367178429, -4.637299300123509,\n",
    "11.697062646578559, -5.2540491532558375,\n",
    "11.797344354961812, -5.838292911958037])\n",
    "smf_mu_q_3 = smf_mu_q_3.reshape(len(smf_mu_q_3)//2,2)\n",
    "np.save('smf_lit/smf_mu_q_1.0_1.5',smf_mu_q_3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "tags": []
   },
   "source": [
    "#### Muzzin et al. 2013 (SF)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5\n",
    "smf_mu_sf_1 = np.array([\n",
    "8.424370461669138, -1.6650112866817148,\n",
    "8.525252623172552, -1.7516930022573356,\n",
    "8.633993454186228, -1.7841986455981933,\n",
    "8.734980865727726, -1.8058690744920987,\n",
    "8.820286021594724, -1.9142212189616246,\n",
    "8.91352001366404, -1.9250564334085771,\n",
    "9.022243303004705, -1.9683972911963874,\n",
    "9.12321317287319, -2.0009029345372453,\n",
    "9.216412081596477, -2.0334085778781033,\n",
    "9.309558365300726, -2.0984198645598187,\n",
    "9.426017532440556, -2.163431151241534,\n",
    "9.519233982836859, -2.1851015801354396,\n",
    "9.612362724868092, -2.260948081264108,\n",
    "9.713297511390548, -2.3151241534988705,\n",
    "9.822108509096282, -2.304288939051918,\n",
    "9.915307417819571, -2.336794582392776,\n",
    "10.016242204342026, -2.390970654627539,\n",
    "10.10937094637326, -2.466817155756207,\n",
    "10.210340816241743, -2.499322799097065,\n",
    "10.311293144437213, -2.5426636568848755,\n",
    "10.412280555978711, -2.56433408577878,\n",
    "10.513180259155138, -2.640180586907449,\n",
    "10.606361626205414, -2.683521444695259,\n",
    "10.70706837097869, -2.8785553047404058,\n",
    "10.815703951954282, -2.976072234762979,\n",
    "10.908639735582362, -3.1711060948081258,\n",
    "11.009381563701663, -3.3444695259593678,\n",
    "11.109930433417812, -3.6370203160270878,\n",
    "11.202532925258625, -4.037923250564334,\n",
    "11.303695753530263, -3.951241534988713,\n",
    "11.403823623094075, -4.503837471783296,\n",
    "11.488988445576961, -4.6988713318284425])\n",
    "smf_mu_sf_1 = smf_mu_sf_1.reshape(len(smf_mu_sf_1)//2,2)\n",
    "np.save('smf_lit/smf_mu_sf_0.2_0.5',smf_mu_sf_1)\n",
    "\n",
    "# 0.5 - 1.0\n",
    "smf_mu_sf_2 = np.array([\n",
    "8.968969026384933, -2.1200902934537247,\n",
    "9.06244680457388, -2.0984198645598195,\n",
    "9.171410465159308, -2.1309255079006775,\n",
    "9.264782777633869, -2.1742663656884877,\n",
    "9.366029863449654, -2.163431151241535,\n",
    "9.459454908781408, -2.1742663656884877,\n",
    "9.552809643636902, -2.228442437923251,\n",
    "9.6538985308811, -2.3151241534988714,\n",
    "9.747393886689117, -2.282618510158014,\n",
    "9.840766199163676, -2.325959367945824,\n",
    "9.941907819265072, -2.380135440180587,\n",
    "10.035244976501499, -2.4451467268623026,\n",
    "10.144103171372535, -2.542663656884876,\n",
    "10.237580949561485, -2.520993227990971,\n",
    "10.338740147281944, -2.564334085778781,\n",
    "10.439899345002406, -2.6076749435665914,\n",
    "10.540900344151279, -2.7485327313769754,\n",
    "10.642112274728936, -2.7593679458239277,\n",
    "10.735431854346299, -2.835214446952596,\n",
    "10.8130546201384, -2.986907449209933,\n",
    "10.921807349295046, -3.1494356659142215,\n",
    "11.02272046034859, -3.3444695259593686,\n",
    "11.123510528068682, -3.615349887133183,\n",
    "11.216548865781002, -3.864559819413093,\n",
    "11.309323539207341, -4.276297968397292,\n",
    "11.402098212633684, -4.68803611738149,\n",
    "11.502220330829294, -5.370654627539504])\n",
    "smf_mu_sf_2 = smf_mu_sf_2.reshape(len(smf_mu_sf_2)//2,2)\n",
    "np.save('smf_lit/smf_mu_sf_0.5_1.0',smf_mu_sf_2)\n",
    "\n",
    "# 1.0 - 1.5\n",
    "smf_mu_sf_3 = np.array([\n",
    "9.529184803474866, -2.379591175912778,\n",
    "9.622628723153285, -2.378791845721873,\n",
    "9.723842058275244, -2.3887612694917695,\n",
    "9.832701763799777, -2.485347000892771,\n",
    "9.926057794163224, -2.538724494752084,\n",
    "10.019396246663675, -2.6029373534214395,\n",
    "10.120486536744673, -2.6887543308616397,\n",
    "10.221717449729626, -2.6878883898214934,\n",
    "10.315020746504086, -2.773771978110936,\n",
    "10.424021074932579, -2.78367479103159,\n",
    "10.517341949570035, -2.8587230145109896,\n",
    "10.618484973240015, -2.91203389752106,\n",
    "10.711735536425493, -3.0304235802406327,\n",
    "10.812773092917508, -3.1487466521109635,\n",
    "10.906023656102988, -3.2671363348305364,\n",
    "11.006920589691045, -3.4721423251812147,\n",
    "11.107782367553115, -3.6988190451519802,\n",
    "11.208573833963207, -3.9688372243629195,\n",
    "11.309084054565387, -4.412221240534553,\n",
    "11.40187759331301, -4.812330408315256,\n",
    "11.494899644279561, -5.071579833565393,\n",
    "11.602634366572456, -5.8616289128091745,\n",
    "11.696113441976859, -5.8391588529981835])\n",
    "smf_mu_sf_3 = smf_mu_sf_3.reshape(len(smf_mu_sf_3)//2,2)\n",
    "np.save('smf_lit/smf_mu_sf_1.0_1.5',smf_mu_sf_3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Davidzon et al. 2017"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5\n",
    "smf_da_all_1 = np.array([\n",
    "9.254201680672269, -1.952121871599565,\n",
    "9.445378151260503, -2.0522306855277477,\n",
    "9.634453781512605, -2.1175190424374324,\n",
    "9.823529411764705, -2.1218715995647446,\n",
    "10.016806722689076, -2.217627856365615,\n",
    "10.201680672268907, -2.335146898803047,\n",
    "10.386554621848738, -2.3743199129488577,\n",
    "10.571428571428571, -2.4047878128400435,\n",
    "10.766806722689076, -2.496191512513602,\n",
    "10.957983193277311, -2.761697497279652,\n",
    "11.140756302521007, -3.096844396082699,\n",
    "11.308823529411764, -3.405875952121872,\n",
    "11.521008403361343, -4.032644178454842,\n",
    "11.655462184873947, -4.733405875952122])\n",
    "smf_da_all_1 = smf_da_all_1.reshape(len(smf_da_all_1)//2,2)\n",
    "np.save('smf_lit/smf_da_0.2_0.5',smf_da_all_1)\n",
    "\n",
    "# 0.5 - 0.8\n",
    "smf_da_all_2 = np.array([\n",
    "9.102941176470587, -1.9651795429815018,\n",
    "9.294117647058822, -2.0609357997823725,\n",
    "9.502100840336134, -2.139281828073994,\n",
    "9.697478991596638, -2.2437431991294887,\n",
    "9.903361344537814, -2.3003264417845486,\n",
    "10.100840336134453, -2.426550598476605,\n",
    "10.304621848739496, -2.4439608269858546,\n",
    "10.497899159663865, -2.4744287268770404,\n",
    "10.695378151260503, -2.58759521218716,\n",
    "10.888655462184873, -2.7399347116430905,\n",
    "11.096638655462183, -2.8574537540805225,\n",
    "11.292016806722689, -3.28835690968444,\n",
    "11.485294117647058, -3.741022850924918,\n",
    "11.678571428571427, -4.755168661588684])\n",
    "smf_da_all_2 = smf_da_all_2.reshape(len(smf_da_all_2)//2,2)\n",
    "np.save('smf_lit/smf_da_0.5_0.8',smf_da_all_2)\n",
    "\n",
    "# 0.8 - 1.1\n",
    "smf_da_all_3 = np.array([\n",
    "9.279411764705882, -2.0783460282916217,\n",
    "9.483193277310924, -2.1828073993471167,\n",
    "9.682773109243696, -2.26550598476605,\n",
    "9.880252100840336, -2.3133841131664856,\n",
    "10.088235294117647, -2.4047878128400435,\n",
    "10.285714285714285, -2.4439608269858546,\n",
    "10.485294117647058, -2.465723612622416,\n",
    "10.686974789915965, -2.574537540805223,\n",
    "10.8781512605042, -2.6572361262241566,\n",
    "11.07563025210084, -2.935799782372144,\n",
    "11.26470588235294, -3.2274211099020675,\n",
    "11.455882352941176, -3.8585418933623505])\n",
    "smf_da_all_3 = smf_da_all_3.reshape(len(smf_da_all_3)//2,2)\n",
    "np.save('smf_lit/smf_da_0.8_1.1',smf_da_all_3)\n",
    "\n",
    "# 1.1 - 1.5\n",
    "smf_da_all_4 = np.array([\n",
    "9.283613445378151, -2.204570184983678,\n",
    "9.497899159663865, -2.309031556039173,\n",
    "9.716386554621849, -2.391730141458107,\n",
    "9.932773109243696, -2.4744287268770404,\n",
    "10.153361344537814, -2.5484221980413495,\n",
    "10.369747899159663, -2.618063112078346,\n",
    "10.581932773109243, -2.6789989118607185,\n",
    "10.798319327731091, -2.7704026115342764,\n",
    "11.008403361344538, -3.0837867247007615,\n",
    "11.214285714285714, -3.410228509249184,\n",
    "11.428571428571427, -4.058759521218716,\n",
    "11.615546218487394, -4.807399347116431,\n",
    "11.80252100840336, -5.586507072905332])\n",
    "smf_da_all_4 = smf_da_all_4.reshape(len(smf_da_all_4)//2,2)\n",
    "np.save('smf_lit/smf_da_1.1_1.5',smf_da_all_4)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Weaver et al. 2022"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5\n",
    "smf_we_all_1 = np.array([\n",
    "8.181640585254675, -1.6421725239616611,\n",
    "8.382155430292425, -1.693290734824281,\n",
    "8.630318878329705, -1.7955271565495206,\n",
    "8.888066990903726, -1.8849840255591053,\n",
    "9.136169390122301, -2.0127795527156547,\n",
    "9.38430231375023, -2.1277955271565494,\n",
    "9.622881097250769, -2.242811501597444,\n",
    "9.880598685415437, -2.3450479233226833,\n",
    "10.13834679798946, -2.4345047923322682,\n",
    "10.38663234366415, -2.4856230031948883,\n",
    "10.634765267292078, -2.6006389776357826,\n",
    "10.882745568873244, -2.7795527156549515,\n",
    "11.13042062636088, -3.0862619808306704,\n",
    "11.3870088113795, -3.6613418530351436,\n",
    "11.634012331861378, -4.2492012779552715,])\n",
    "smf_we_all_1 = smf_we_all_1.reshape(len(smf_we_all_1)//2,2)\n",
    "np.save('smf_lit/smf_we_0.2_0.5',smf_we_all_1)\n",
    "    \n",
    "# 0.5 - 0.8\n",
    "smf_we_all_2 = np.array([\n",
    "8.402896669236076, -1.7972331534747854,\n",
    "8.631252789159804, -1.8876627855085395,\n",
    "8.631252789159804, -1.8876627855085395,\n",
    "8.878656334942594, -1.9781735567365812,\n",
    "9.135553166457058, -2.0815043206620953,\n",
    "9.37340257211246, -2.1847539453933216,\n",
    "9.630299403626923, -2.288084709318836,\n",
    "9.877672522211855, -2.3913749036472063,\n",
    "10.134599780924175, -2.4819262444723917,\n",
    "10.382064181102683, -2.5468781694997764,\n",
    "10.639021867012861, -2.6246500872246337,\n",
    "10.886181995212787, -2.817396243255304,\n",
    "11.133068278631995, -3.125157207188932,\n",
    "11.379467726885471, -3.637388940727818,\n",
    "11.625380339973225, -4.354091443871962,\n",
    "11.871201671467402, -5.109132216317091])\n",
    "smf_we_all_2 = smf_we_all_2.reshape(len(smf_we_all_2)//2,2)\n",
    "np.save('smf_lit/smf_we_0.5_0.8',smf_we_all_2)\n",
    "\n",
    "# 0.8 - 1.1\n",
    "smf_we_all_3 = np.array([\n",
    "8.488270196257417, -1.92332268370607,\n",
    "8.650144530655712, -1.936102236421725,\n",
    "8.888087631218621, -2.000000000000000,\n",
    "9.135493686292408, -2.0894568690095845,\n",
    "9.39245397839647, -2.1661341853035143,\n",
    "9.63982960596379, -2.268370607028754,\n",
    "9.887235661037577, -2.3578274760383384,\n",
    "10.13467214361783, -2.4345047923322682,\n",
    "10.39163243572189, -2.511182108626198,\n",
    "10.6295755362848, -2.5750798722044728,\n",
    "10.886322835843602, -2.741214057507987,\n",
    "11.133211623307469, -3.047923322683706,\n",
    "11.389198235204622, -3.533546325878594,\n",
    "11.63514376996805, -4.236421725239616,\n",
    "11.87101780009128, -5.1693290734824275,\n",
    "12.126517571884984, -5.859424920127795])\n",
    "smf_we_all_3 = smf_we_all_3.reshape(len(smf_we_all_3)//2,2)\n",
    "np.save('smf_lit/smf_we_0.8_1.1',smf_we_all_3)\n",
    "\n",
    "# 1.1 - 1.5\n",
    "smf_we_all_4 = np.array([\n",
    "8.686708860759492, -2.0511182108626196,\n",
    "8.87658227848101, -2.0638977635782743,\n",
    "9.113924050632914, -2.1277955271565494,\n",
    "9.37974683544304, -2.230031948881789,\n",
    "9.62658227848101, -2.3322683706070286,\n",
    "9.873417721518987, -2.4472843450479234,\n",
    "10.110759493670887, -2.536741214057508,\n",
    "10.37658227848101, -2.6389776357827475,\n",
    "10.613924050632914, -2.741214057507987,\n",
    "10.860759493670887, -2.8945686900958463,\n",
    "11.117088607594937, -3.2268370607028753,\n",
    "11.373417721518987, -3.6741214057507983,\n",
    "11.620253164556964, -4.415335463258786,\n",
    "11.867088607594937, -5.6293929712460065,\n",
    "12.123417721518987, -6.07667731629393])\n",
    "smf_we_all_4 = smf_we_all_4.reshape(len(smf_we_all_4)//2,2)\n",
    "np.save('smf_lit/smf_we_1.1_1.5',smf_we_all_4)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Weaver et al. 2022 (SF)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5\n",
    "smf_we_sf_1 = np.array([\n",
    "8.137096774193548, -1.7115044247787605,\n",
    "8.387096774193548, -1.7433628318584067,\n",
    "8.629032258064516, -1.8495575221238933,\n",
    "8.887096774193548, -1.9345132743362827,\n",
    "9.129032258064516, -2.0619469026548667,\n",
    "9.379032258064516, -2.1999999999999993,\n",
    "9.629032258064516, -2.316814159292035,\n",
    "9.887096774193548, -2.4548672566371676,\n",
    "10.137096774193548, -2.5823008849557514,\n",
    "10.387096774193548, -2.6566371681415926,\n",
    "10.629032258064516, -2.8477876106194686,\n",
    "10.879032258064516, -3.060176991150442,\n",
    "11.137096774193548, -3.4849557522123886,\n",
    "11.379032258064516, -4.122123893805309])\n",
    "smf_we_sf_1 = smf_we_sf_1.reshape(len(smf_we_sf_1)//2,2)\n",
    "np.save('smf_lit/smf_we_sf_0.2_0.5',smf_we_sf_1)\n",
    "    \n",
    "# 0.5 - 0.8\n",
    "smf_we_sf_2 = np.array([\n",
    "8.371717171717172, -1.8424027509133898,\n",
    "8.622222222222222, -1.918203309692672,\n",
    "8.872727272727273, -2.015280464216635,\n",
    "9.123232323232322, -2.112357618740598,\n",
    "9.365656565656565, -2.2200300881151955,\n",
    "9.616161616161616, -2.3490221362561794,\n",
    "9.874747474747476, -2.46741886954653,\n",
    "10.117171717171717, -2.6176445304104887,\n",
    "10.375757575757575, -2.714764667956158,\n",
    "10.626262626262626, -2.865033311841823,\n",
    "10.876767676767678, -3.079131742961531,\n",
    "11.11919191919192, -3.5059531485063404,\n",
    "11.377777777777778, -4.145626477541372,\n",
    "11.62020202020202, -5.104362776703203])\n",
    "smf_we_sf_2 = smf_we_sf_2.reshape(len(smf_we_sf_2)//2,2)\n",
    "np.save('smf_lit/smf_we_sf_0.5_0.8',smf_we_sf_2)\n",
    "    \n",
    "# 0.8 - 1.1\n",
    "smf_we_sf_3 = np.array([\n",
    "8.457167369700436, -1.9231429184242512,\n",
    "8.650632730392509, -1.9540287409737616,\n",
    "8.892428683777794, -2.0059340816472435,\n",
    "9.134110245227713, -2.1003932222778294,\n",
    "9.383885036104953, -2.1841710159433756,\n",
    "9.641638664474154, -2.299821262600987,\n",
    "10.13289483091442, -2.552527346822049,\n",
    "10.632272824765856, -2.7839136340887967,\n",
    "10.881790233788518, -2.963437477657825,\n",
    "11.13085007506971, -3.3131765210552655,\n",
    "11.379280760706369, -3.8969614642167727,\n",
    "11.618788875384286, -4.799942804032315,\n",
    "11.882919854150284, -5.543218703081432])\n",
    "smf_we_sf_3 = smf_we_sf_3.reshape(len(smf_we_sf_3)//2,2)\n",
    "np.save('smf_lit/smf_we_sf_0.8_1.1',smf_we_sf_3)\n",
    "    \n",
    "# 1.1 -1.5\n",
    "smf_we_sf_4 = np.array([\n",
    "8.678787878787878, -2.053129525341916,\n",
    "8.872727272727271, -2.0748833467417542,\n",
    "9.131313131313131, -2.1499061410565856,\n",
    "9.373737373737374, -2.246124966478951,\n",
    "9.874747474747474, -2.491702869401986,\n",
    "10.125252525252524, -2.6198015553767773,\n",
    "10.383838383838384, -2.758541163850899,\n",
    "10.626262626262626, -2.9078573344060077,\n",
    "10.884848484848483, -3.1315526950925188,\n",
    "11.127272727272727, -3.461399839098955,\n",
    "11.377777777777776, -3.8549852507374647,\n",
    "11.620202020202019, -4.598991686779299,\n",
    "11.878787878787879, -5.810297666934834])\n",
    "smf_we_sf_4 = smf_we_sf_4.reshape(len(smf_we_sf_4)//2,2)\n",
    "np.save('smf_lit/smf_we_sf_1.1_1.5',smf_we_sf_4)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Weaver et al. 2022 (Q)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5\n",
    "smf_we_q_1 = np.array([\n",
    "8.129032258064516, -2.380530973451327,\n",
    "8.39516129032258, -2.635398230088495,\n",
    "8.637096774193548, -2.8265486725663713,\n",
    "8.879032258064516, -2.890265486725663,\n",
    "9.137096774193548, -2.9858407079646017,\n",
    "9.379032258064516, -3.0707964601769904,\n",
    "9.637096774193548, -3.060176991150442,\n",
    "9.879032258064516, -3.007079646017699,\n",
    "10.137096774193548, -3.0176991150442474,\n",
    "10.379032258064516, -3.007079646017699,\n",
    "10.637096774193548, -2.975221238938053,\n",
    "10.887096774193548, -3.1026548672566365,\n",
    "11.129032258064516, -3.3150442477876103,\n",
    "11.387096774193548, -3.846017699115044,\n",
    "11.629032258064516, -4.260176991150442])\n",
    "smf_we_q_1 = smf_we_q_1.reshape(len(smf_we_q_1)//2,2)\n",
    "np.save('smf_lit/smf_we_q_0.2_0.5',smf_we_q_1)\n",
    "    \n",
    "# 0.5 - 0.8\n",
    "smf_we_q_2 = np.array([\n",
    "8.460606060606061, -2.8535138620245006,\n",
    "8.614141414141415, -3.0245433053943698,\n",
    "8.872727272727273, -3.174854932301742,\n",
    "9.123232323232322, -3.2932086825703855,\n",
    "9.373737373737374, -3.3583709434773272,\n",
    "9.616161616161616, -3.3596604341285197,\n",
    "9.874747474747476, -3.180184826993339,\n",
    "10.117171717171717, -3.107006232538148,\n",
    "10.375757575757575, -3.02327530625403,\n",
    "10.618181818181819, -3.0564796905222447,\n",
    "10.876767676767678, -3.1429615301955733,\n",
    "11.11919191919192, -3.357016978293575,\n",
    "11.36969696969697, -3.794519664732431,\n",
    "11.628282828282828, -4.444831291639803,\n",
    "11.878787878787879, -5.105738233397808])\n",
    "smf_we_q_2 = smf_we_q_2.reshape(len(smf_we_q_2)//2,2)\n",
    "np.save('smf_lit/smf_we_q_0.5_0.8',smf_we_q_2)\n",
    "    \n",
    "# 0.8 - 1.1\n",
    "smf_we_q_3 = np.array([\n",
    "8.880817902337885, -3.325144777293201,\n",
    "9.130707085150497, -3.3663687710016443,\n",
    "9.388775291341958, -3.3649960677772217,\n",
    "9.630685636662616, -3.374347608493601,\n",
    "9.88111818116823, -3.213441052405805,\n",
    "10.139615357117322, -3.052491599342246,\n",
    "10.381783084292557, -2.966097090155144,\n",
    "10.639851290484021, -2.964724386930721,\n",
    "10.873468220490455, -3.05922642453707,\n",
    "11.130792879102025, -3.334453421033817,\n",
    "11.379709730464, -3.737384714377636,\n",
    "11.619989990705655, -4.353127904482733,\n",
    "11.875341388432116, -5.362407950239507,\n",
    "12.124058053907198, -5.839808393508258])\n",
    "smf_we_q_3 = smf_we_q_3.reshape(len(smf_we_q_3)//2,2)\n",
    "np.save('smf_lit/smf_we_q_0.8_1.1',smf_we_q_3)\n",
    "    \n",
    "# 1.1 -1.5\n",
    "smf_we_q_4 = np.array([\n",
    "9.123232323232322, -3.8914776079377864,\n",
    "9.373737373737374, -3.9877178868329315,\n",
    "9.624242424242425, -3.9140466613033,\n",
    "9.874747474747474, -3.6704639313488885,\n",
    "10.125252525252524, -3.4375006704210245,\n",
    "10.375757575757575, -3.2894931617055514,\n",
    "10.634343434343434, -3.2477017967283457,\n",
    "10.876767676767678, -3.343920622150711,\n",
    "11.127272727272727, -3.652550281576832,\n",
    "11.36969696969697, -4.247884151246983,\n",
    "11.628282828282828, -4.89635827299544])\n",
    "smf_we_q_4 = smf_we_q_4.reshape(len(smf_we_q_4)//2,2)\n",
    "np.save('smf_lit/smf_we_q_1.1_1.5',smf_we_q_4)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Moutard et al. 2016 (Q)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5 \n",
    "data_q1 = np.array([8.9170, -2.7438,\n",
    "9.0377, -2.8086,\n",
    "9.1584, -2.9081,\n",
    "9.2867, -2.9816,\n",
    "9.4049, -3.0508,\n",
    "9.5256, -3.1027,\n",
    "9.6513, -3.0984,\n",
    "9.7695, -3.0854,\n",
    "9.8927, -3.0422,\n",
    "10.018, -2.9773,\n",
    "10.139, -2.9341,\n",
    "10.262, -2.9168,\n",
    "10.383, -2.8692,\n",
    "10.506, -2.8692,\n",
    "10.627, -2.8605,\n",
    "10.755, -2.9038,\n",
    "10.873, -2.9946,\n",
    "10.994, -3.0897,\n",
    "11.120, -3.2800,\n",
    "11.243, -3.4703,\n",
    "11.364, -3.7946,\n",
    "11.487, -4.1665,\n",
    "11.605, -4.5081,\n",
    "11.731, -4.9276,\n",
    "11.839, -5.5330])\n",
    "data_q1 = data_q1.reshape(len(data_q1)//2,2)\n",
    "np.save('smf_lit/smf_mo_q_0.2_0.5',data_q1)\n",
    "\n",
    "# 0.5 - 0.8\n",
    "data_q2 = np.array([9.6745, -3.3924,\n",
    "9.7626, -3.2670,\n",
    "9.8557, -3.2368,\n",
    "9.9463, -3.2022,\n",
    "10.049, -3.1589,\n",
    "10.135, -3.1157,\n",
    "10.231, -3.0984,\n",
    "10.319, -3.0595,\n",
    "10.412, -3.0292,\n",
    "10.508, -3.0032,\n",
    "10.601, -3.0335,\n",
    "10.701, -3.0119,\n",
    "10.789, -3.0595,\n",
    "10.883, -3.1027,\n",
    "10.978, -3.1762,\n",
    "11.066, -3.2800,\n",
    "11.167, -3.4357,\n",
    "11.268, -3.5870,\n",
    "11.348, -3.7859,\n",
    "11.439, -4.0584,\n",
    "11.537, -4.3568,\n",
    "11.615, -4.5989,\n",
    "11.708, -5.2995,\n",
    "11.801, -5.5114,\n",
    "11.924, -6.0130])\n",
    "data_q2 = data_q2.reshape(len(data_q2)//2,2)\n",
    "np.save('smf_lit/smf_mo_q_0.5_0.8',data_q2)\n",
    "\n",
    "# 0.8 - 1.1\n",
    "data_q3 = np.array([10.116, -3.2757,\n",
    "10.290, -3.2151,\n",
    "10.476, -3.0465,\n",
    "10.657, -2.9730,\n",
    "10.836, -3.0292,\n",
    "11.020, -3.2022,\n",
    "11.199, -3.5049,\n",
    "11.383, -3.9892,\n",
    "11.567, -4.6551,\n",
    "11.746, -5.4897])\n",
    "data_q3 = data_q3.reshape(len(data_q3)//2,2)\n",
    "np.save('smf_lit/smf_mo_q_0.8_1.1',data_q3)\n",
    "\n",
    "# 1.1 - 1.5\n",
    "data_q4 = np.array([10.354, -4.0195,\n",
    "10.437, -3.8681,\n",
    "10.545, -3.7600,\n",
    "10.654, -3.7514,\n",
    "10.759, -3.5654,\n",
    "10.857, -3.5178,\n",
    "10.956, -3.4962,\n",
    "11.051, -3.5697,\n",
    "11.172, -3.6908,\n",
    "11.263, -3.8638,\n",
    "11.371, -4.2141,\n",
    "11.474, -4.5600,\n",
    "11.577, -5.0486,\n",
    "11.673, -5.4465,\n",
    "11.781, -5.9438,\n",
    "11.884, -6.0951])\n",
    "data_q4 = data_q4.reshape(len(data_q4)//2,2)\n",
    "np.save('smf_lit/smf_mo_q_1.1_1.5',data_q4)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Moutard et al. 2016 (SF)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5 \n",
    "data_sf1 = np.array([8.7888, -1.9092,\n",
    "8.9070, -1.9697,\n",
    "9.0427, -1.9784,\n",
    "9.1584, -2.0389,\n",
    "9.2816, -2.0865,\n",
    "9.3973, -2.1514,\n",
    "9.5331, -2.2162,\n",
    "9.6513, -2.2854,\n",
    "9.7921, -2.3286,\n",
    "9.9053, -2.3373,\n",
    "10.023, -2.4108,\n",
    "10.154, -2.4368,\n",
    "10.275, -2.5059,\n",
    "10.406, -2.5535,\n",
    "10.519, -2.6443,\n",
    "10.645, -2.7395,\n",
    "10.758, -2.9038,\n",
    "10.891, -3.0465,\n",
    "11.009, -3.3189,\n",
    "11.132, -3.5611,\n",
    "11.261, -3.9373,\n",
    "11.386, -4.3438,\n",
    "11.502, -4.9405,\n",
    "11.635, -5.9568])\n",
    "data_sf1 = data_sf1.reshape(len(data_sf1)//2,2)\n",
    "np.save('smf_lit/smf_mo_sf_0.2_0.5',data_sf1)  \n",
    "\n",
    "# 0.5 - 0.8\n",
    "data_sf2 = np.array([9.4857, -2.2292,\n",
    "9.6871, -2.2811,\n",
    "9.7978, -2.3373,\n",
    "9.8935, -2.3978,\n",
    "10.017, -2.4497,\n",
    "10.128, -2.4973,\n",
    "10.226, -2.5189,\n",
    "10.341, -2.5578,\n",
    "10.440, -2.5838,\n",
    "10.550, -2.6530,\n",
    "10.661, -2.7568,\n",
    "10.762, -2.8519,\n",
    "10.878, -2.9686,\n",
    "10.983, -3.1589,\n",
    "11.079, -3.3924,\n",
    "11.197, -3.6822,\n",
    "11.293, -4.0195,\n",
    "11.421, -4.3395,\n",
    "11.524, -4.8368,\n",
    "11.628, -5.3946,\n",
    "11.711, -6.0346,\n",
    "11.985, -6.2249])\n",
    "data_sf2 = data_sf2.reshape(len(data_sf2)//2,2)\n",
    "np.save('smf_lit/smf_mo_sf_0.5_0.8',data_sf2)  \n",
    "\n",
    "\n",
    "# 0.8 - 1.1\n",
    "data_sf3 = np.array([9.9622, -2.4843,\n",
    "10.038, -2.4541,\n",
    "10.108, -2.5103,\n",
    "10.194, -2.5319,\n",
    "10.262, -2.5665,\n",
    "10.325, -2.5838,\n",
    "10.406, -2.6141,\n",
    "10.489, -2.6400,\n",
    "10.549, -2.6573,\n",
    "10.625, -2.7135,\n",
    "10.713, -2.7697,\n",
    "10.786, -2.8389,\n",
    "10.864, -2.9168,\n",
    "10.945, -3.0162,\n",
    "11.005, -3.1416,\n",
    "11.081, -3.2843,\n",
    "11.146, -3.4616,\n",
    "11.222, -3.7038,\n",
    "11.297, -3.9416,\n",
    "11.380, -4.2011,\n",
    "11.456, -4.5903,\n",
    "11.597, -5.3470,\n",
    "11.748, -6.1038])\n",
    "data_sf3 = data_sf3.reshape(len(data_sf3)//2,2)\n",
    "np.save('smf_lit/smf_mo_sf_0.8_1.1',data_sf3)  \n",
    "\n",
    "# 1.1 - 1.5\n",
    "data_sf4 = np.array([10.284, -2.7827,\n",
    "10.321, -2.7611,\n",
    "10.359, -2.7308,\n",
    "10.399, -2.7178,\n",
    "10.442, -2.7827,\n",
    "10.497, -2.7957,\n",
    "10.568, -2.8476,\n",
    "10.648, -2.8822,\n",
    "10.724, -2.9211,\n",
    "10.789, -2.9643,\n",
    "10.870, -3.0378,\n",
    "10.945, -3.1200,\n",
    "11.013, -3.2151,\n",
    "11.081, -3.3319,\n",
    "11.159, -3.4876,\n",
    "11.230, -3.6562,\n",
    "11.313, -3.8941,\n",
    "11.383, -4.2097,\n",
    "11.451, -4.4995,\n",
    "11.527, -4.8368,\n",
    "11.592, -5.0400,\n",
    "11.663, -5.5459,\n",
    "11.731, -6.2076])\n",
    "data_sf4 = data_sf4.reshape(len(data_sf4)//2,2)\n",
    "np.save('smf_lit/smf_mo_sf_1.1_1.5',data_sf4)  \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 0.2 - 0.5 \n",
    "data_all1 = np.array([\n",
    "8.806797853309481, -1.8715933625316152,\n",
    "8.930232558139535, -1.9179487179487191,\n",
    "9.059033989266547, -1.9275389139061962,\n",
    "9.1824686940966, -1.9922850739210012,\n",
    "9.305903398926654, -2.0386404293381046,\n",
    "9.440071556350626, -2.1126148911233122,\n",
    "9.56350626118068, -2.1681656488392664,\n",
    "9.692307692307692, -2.205342051693296,\n",
    "9.821109123434704, -2.260909259145026,\n",
    "9.949910554561717, -2.2888902597002048,\n",
    "10.07334525939177, -2.2984640059219057,\n",
    "10.196779964221825, -2.3448193613390096,\n",
    "10.330948121645797, -2.3728168116299644,\n",
    "10.45438282647585, -2.419172167047068,\n",
    "10.583184257602863, -2.483934776797649,\n",
    "10.706618962432916, -2.5486809368124534,\n",
    "10.835420393559929, -2.6778113626549884,\n",
    "10.964221824686941, -2.852918799991776,\n",
    "11.087656529516995, -3.0280097875927874,\n",
    "11.216457960644007, -3.3042666502169316,\n",
    "11.33989266547406, -3.617288672300702,\n",
    "11.463327370304114, -4.05904632656838,\n",
    "11.592128801431127, -4.390475602985628,\n",
    "11.726296958855098, -4.951806386609915,\n",
    "11.833631484794275, -5.476273312359921])\n",
    "data_all1 = data_all1.reshape(len(data_all1)//2,2)\n",
    "np.save('smf_lit/smf_mo_0.2_0.5',data_all1)  \n",
    "\n",
    "# 0.5 - 0.8\n",
    "data_all2 = np.array([\n",
    "9.648214285714285, -2.2390804597701157,\n",
    "9.755357142857143, -2.275862068965518,\n",
    "9.8625, -2.3034482758620696,\n",
    "9.964285714285714, -2.3770114942528737,\n",
    "10.066071428571428, -2.3862068965517245,\n",
    "10.167857142857143, -2.3862068965517245,\n",
    "10.275, -2.441379310344828,\n",
    "10.382142857142856, -2.4505747126436788,\n",
    "10.48392857142857, -2.4781609195402305,\n",
    "10.591071428571428, -2.5149425287356326,\n",
    "10.692857142857143, -2.570114942528736,\n",
    "10.794642857142858, -2.6436781609195403,\n",
    "10.907142857142857, -2.7540229885057474,\n",
    "11.003571428571428, -2.8919540229885063,\n",
    "11.110714285714284, -3.1034482758620694,\n",
    "11.212499999999999, -3.3149425287356324,\n",
    "11.319642857142856, -3.563218390804598,\n",
    "11.42142857142857, -3.839080459770115,\n",
    "11.528571428571428, -4.225287356321839,\n",
    "11.625, -4.583908045977012,\n",
    "11.732142857142858, -5.227586206896552,\n",
    "11.8125, -5.586206896551724,\n",
    "11.946428571428571, -5.816091954022988])\n",
    "data_all2 = data_all2.reshape(len(data_all2)//2,2)\n",
    "np.save('smf_lit/smf_mo_0.5_0.8',data_all2)  \n",
    "\n",
    "# 0.8 - 1.1\n",
    "data_all3 = np.array([\n",
    "10.051878354203936, -2.422988505747127,\n",
    "10.223613595706619, -2.459770114942529,\n",
    "10.395348837209301, -2.4781609195402297,\n",
    "10.561717352415027, -2.496551724137931,\n",
    "10.733452593917711, -2.5793103448275856,\n",
    "10.91592128801431, -2.726436781609195,\n",
    "11.103756708407872, -3.0022988505747126,\n",
    "11.248658318425761, -3.44367816091954,\n",
    "11.436493738819319, -4.0137931034482754,\n",
    "11.613595706618963, -4.685057471264367,\n",
    "11.774597495527729, -5.678160919540227])\n",
    "data_all3 = data_all3.reshape(len(data_all3)//2,2)\n",
    "np.save('smf_lit/smf_mo_0.8_1.1',data_all3)  \n",
    "\n",
    "# 1.1 - 1.5\n",
    "data_all4 = np.array([\n",
    "10.4030225352924, -2.724108278027165,\n",
    "10.499401335577126, -2.7612707320986996,\n",
    "10.569006875383167, -2.7891343435881195,\n",
    "10.64933900041557, -2.8078143193946703,\n",
    "10.740373023481634, -2.8357437633979465,\n",
    "10.820680461321341, -2.8728568430840893,\n",
    "10.906332676379707, -2.919202932838493,\n",
    "10.986566052641326, -3.0116153241634143,\n",
    "11.066799428902943, -3.1040277154883347,\n",
    "11.15230352080513, -3.260972428520292,\n",
    "11.237795269110974, -3.4271336934920456,\n",
    "11.312474027016018, -3.6669944577252416,\n",
    "11.39778062137664, -3.971404001793937,\n",
    "11.483037841351871, -4.312679753621818,\n",
    "11.568344435712493, -4.617089297690516,\n",
    "11.648096411716542, -5.068947214667485,\n",
    "11.72786073131694, -5.511588579704661,\n",
    "11.818487415703522, -5.8436642377212085,\n",
    "11.898757822754185, -5.908426973226742])\n",
    "data_all4 = data_all4.reshape(len(data_all4)//2,2)\n",
    "np.save('smf_lit/smf_mo_1.1_1.5',data_all4)  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "####  Weaver et al. (original data)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/lejay/anaconda3/lib/python3.6/site-packages/ipykernel_launcher.py:12: RuntimeWarning: invalid value encountered in log10\n",
      "  if sys.path[0] == '':\n"
     ]
    }
   ],
   "source": [
    "zmins = [0.2,0.5,0.8,1.1]\n",
    "zmaxs = [0.5,0.8,1.1,1.5]\n",
    "\n",
    "for z_idx in range(len(zmins)):\n",
    "    for sfq in ['total','sfg','qg']:\n",
    "        # load fits file\n",
    "        smf_fits = Table.read('/Users/lejay/Desktop/SMF_Farmer_Weaver/observed_smf/SMF_Farmer_'+str(zmins[z_idx])+'z'+str(zmaxs[z_idx])+'_'+sfq+'.fits')\n",
    "        \n",
    "        # formatting\n",
    "        bin_centers = smf_fits['bin_center']\n",
    "        phi = np.log10(smf_fits['phi'])\n",
    "        phi_low = np.log10(smf_fits['phi_lo'])\n",
    "        phi_high = np.log10(smf_fits['phi_hi'])\n",
    "        smf = np.vstack([bin_centers,phi,phi-phi_low,phi_high-phi]).T\n",
    "        \n",
    "        if sfq == 'total':\n",
    "            np.save('smf_lit/smf_we_original_'+str(zmins[z_idx])+'_'+str(zmaxs[z_idx]),smf)\n",
    "        else:\n",
    "            np.save('smf_lit/smf_we_original_'+sfq+'_'+str(zmins[z_idx])+'_'+str(zmaxs[z_idx]),smf)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ 8.17798342, -1.63871806,  0.05339747,  0.04754513],\n",
       "       [ 8.375     , -1.69107886,  0.05348281,  0.04761276],\n",
       "       [ 8.625     , -1.79954317,  0.05423603,  0.04820863],\n",
       "       [ 8.875     , -1.889414  ,  0.05469125,  0.04856786],\n",
       "       [ 9.125     , -2.00528432,  0.05544331,  0.04915986],\n",
       "       [ 9.375     , -2.13554102,  0.05657339,  0.05004599],\n",
       "       [ 9.625     , -2.24290559,  0.0582318 ,  0.05133898],\n",
       "       [ 9.875     , -2.3483682 ,  0.06096254,  0.05344902],\n",
       "       [10.125     , -2.44225597,  0.06444148,  0.05610344],\n",
       "       [10.375     , -2.49265621,  0.06894686,  0.05948598],\n",
       "       [10.625     , -2.60470956,  0.07539113,  0.06421917],\n",
       "       [10.875     , -2.777747  ,  0.08456391,  0.07075111],\n",
       "       [11.125     , -3.08800363,  0.09872198,  0.08038504],\n",
       "       [11.375     , -3.66672278,  0.13144676,  0.10076852],\n",
       "       [11.625     , -4.25879855,  0.21126458,  0.14151179]])"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.load('smf_lit/smf_we_original_0.2_0.5.npy')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### cosmic density "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Weaver + 23\n",
    "\n",
    "density_we = np.array([\n",
    "0.3666933546837472, 155204842.19980335,\n",
    "0.6581265012009607, 150817004.38390854,\n",
    "0.9663730984787833, 173759232.85817063,\n",
    "1.3082465972778223, 124439762.12907402,\n",
    "1.7622097678142516, 75155569.41107331])\n",
    "\n",
    "density_we = density_we.reshape(len(density_we)//2,2)\n",
    "np.save('smf_lit/density_we',density_we) \n",
    "\n",
    "# Madau&Dickenson 2014\n",
    "\n",
    "density_md = np.array([\n",
    "0.15451122577046417, 8.758002560819463,\n",
    "0.5165466244813617, 8.742637644046095,\n",
    "0.9011983693339154, 8.723431498079385,\n",
    "1.4385502936611565, 8.68886043533931,\n",
    "1.9532279183652772, 8.642765685019207,\n",
    "2.484762761990485, 8.565941101152369,\n",
    "3.038841509905518, 8.46606914212548,\n",
    "3.530555067579192, 8.331626120358516])\n",
    "\n",
    "density_md = density_md.reshape(len(density_md)//2,2)\n",
    "np.save('smf_lit/density_md',density_md) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}