From 24e30d0fca5ee516b9388bb350ca967202fda0f3 Mon Sep 17 00:00:00 2001 From: michaelJwilson Date: Tue, 18 Jan 2022 19:23:12 +0000 Subject: [PATCH] Fix ups. --- gen_gold_lf.py | 2 +- logs/gama_lumfn.log | 2 + logs/gama_lumfn_d8_G12.log | 242 +++++++++++++++++++++++++++++++++++++ logs/gama_lumfn_d8_G15.log | 242 +++++++++++++++++++++++++++++++++++++ logs/gama_lumfn_d8_G9.log | 242 +++++++++++++++++++++++++++++++++++++ lumfn.py | 10 +- 6 files changed, 735 insertions(+), 5 deletions(-) diff --git a/gen_gold_lf.py b/gen_gold_lf.py index cb540377..3a61ec9b 100644 --- a/gen_gold_lf.py +++ b/gen_gold_lf.py @@ -53,7 +53,7 @@ def process_cat(fpath, vmax_opath, field=None): ## Luminosity fn. opath = opath.replace('vmax', 'lumfn') - result = lumfn(gama_vmax, VV) + result = lumfn(gama_vmax) print('Writing {}.'.format(opath)) diff --git a/logs/gama_lumfn.log b/logs/gama_lumfn.log index 006f64cf..b2e0c543 100644 --- a/logs/gama_lumfn.log +++ b/logs/gama_lumfn.log @@ -6,3 +6,5 @@ Found redshift limits: 0.039 < z < 0.263 Retrieved area 180.0000 [sq. deg.] WARNING: Found 0.075% with zmax < 0.0 Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_lumfn.fits. +Done. diff --git a/logs/gama_lumfn_d8_G12.log b/logs/gama_lumfn_d8_G12.log index e8a5f597..aadb00ea 100644 --- a/logs/gama_lumfn_d8_G12.log +++ b/logs/gama_lumfn_d8_G12.log @@ -8,3 +8,245 @@ Found redshift limits: 0.040 < z < 0.248 Retrieved area 60.0000 [sq. deg.] WARNING: Found 0.000% with zmax < 0.0 Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_0_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_0_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... ---- -------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.740951649692917 1.3959422305590999e-05 ... 3.0 0.01711404067596624 +-16.659826078008454 3.257198537971233e-05 ... 7.0 0.015012879240296317 +-16.561016966744237 4.653140768530333e-05 ... 10.0 0.012429700350083138 + -16.47329403935433 1.3959422305590999e-05 ... 3.0 0.010470633869063644 +-16.349817938275983 2.7918844611181998e-05 ... 6.0 0.008138253432840999 +-16.231860848664375 9.306281537060667e-06 ... 2.0 0.006280074504121661 +-16.139512804674098 9.306281537060667e-06 ... 2.0 0.004865738155288612 + -16.06060559061654 4.653140768530333e-06 ... 1.0 0.003986100337996181 +-15.959100249556663 4.653140768530333e-06 ... 1.0 0.002837409274042242 + -- 0.0 ... 0.0 -- +-15.759701100596667 4.653140768530333e-06 ... 1.0 0.001052332197901758 +Length = 99 rows +Found vol renormalisation scale of 0.301 +Found d8 renormalisation scale of -0.915 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... ---- -------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.740951649692917 4.6428299804523936e-05 ... 3.0 0.01711404067596624 +-16.659826078008454 0.0001083326995438892 ... 7.0 0.015012879240296317 +-16.561016966744237 0.00015476099934841312 ... 10.0 0.012429700350083138 + -16.47329403935433 4.6428299804523936e-05 ... 3.0 0.010470633869063644 +-16.349817938275983 9.285659960904787e-05 ... 6.0 0.008138253432840999 +-16.231860848664375 3.0952199869682624e-05 ... 2.0 0.006280074504121661 +-16.139512804674098 3.0952199869682624e-05 ... 2.0 0.004865738155288612 + -16.06060559061654 1.5476099934841312e-05 ... 1.0 0.003986100337996181 +-15.959100249556663 1.5476099934841312e-05 ... 1.0 0.002837409274042242 + -- 0.0 ... 0.0 -- +-15.759701100596667 1.5476099934841312e-05 ... 1.0 0.001052332197901758 +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_0_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_1.fits +Found fields: [b'G12'] +Found redshift limits: 0.039 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.058% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_1_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_1_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- --------------------- ... ---- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.749798383390832 5.581790706076068e-05 ... 12.0 0.017526295044526786 +-16.636688775758195 2.790895353038034e-05 ... 6.0 0.014333657175260698 + -16.57595900591837 2.790895353038034e-05 ... 6.0 0.012779287125229681 + -16.45834523178861 1.860596902025356e-05 ... 4.0 0.01022516017150832 +-16.384969299837056 1.395447676519017e-05 ... 3.0 0.008684354016209565 + -16.26855564991519 1.395447676519017e-05 ... 3.0 0.006797678586306946 +-16.146461558177364 1.860596902025356e-05 ... 4.0 0.0051454448356892295 +-16.059952645794596 1.860596902025356e-05 ... 4.0 0.004104126782888122 +-15.946487415059092 4.65149225506339e-06 ... 1.0 0.0028093178511210694 +-15.883966293879894 9.30298451012678e-06 ... 2.0 0.0022595563171305165 +-15.724478023364936 4.65149225506339e-06 ... 1.0 0.0009319416757073396 +Length = 99 rows +Found vol renormalisation scale of 0.155 +Found d8 renormalisation scale of -0.224 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... ---- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.749798383390832 0.00036076544385949303 ... 12.0 0.017526295044526786 +-16.636688775758195 0.00018038272192974652 ... 6.0 0.014333657175260698 + -16.57595900591837 0.00018038272192974652 ... 6.0 0.012779287125229681 + -16.45834523178861 0.00012025514795316433 ... 4.0 0.01022516017150832 +-16.384969299837056 9.019136096487326e-05 ... 3.0 0.008684354016209565 + -16.26855564991519 9.019136096487326e-05 ... 3.0 0.006797678586306946 +-16.146461558177364 0.00012025514795316433 ... 4.0 0.0051454448356892295 +-16.059952645794596 0.00012025514795316433 ... 4.0 0.004104126782888122 +-15.946487415059092 3.0063786988291084e-05 ... 1.0 0.0028093178511210694 +-15.883966293879894 6.012757397658217e-05 ... 2.0 0.0022595563171305165 +-15.724478023364936 3.0063786988291084e-05 ... 1.0 0.0009319416757073396 +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_1_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_2.fits +Found fields: [b'G12'] +Found redshift limits: 0.042 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.054% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_2_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_2_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.72716766679641 2.327808668983176e-05 ... 5.0 0.015947328819278253 + -16.65606635089379 1.8622469351865406e-05 ... 4.0 0.013912916004290763 + -16.54096769997666 2.327808668983176e-05 ... 5.0 0.01114373937443421 + -16.44783744420218 2.327808668983176e-05 ... 5.0 0.009034487013940686 +-16.346886683136404 2.7933704027798108e-05 ... 6.0 0.007350956265230744 +-16.268333154532286 9.311234675932703e-06 ... 2.0 0.006036069915045775 + -16.1189057603977 9.311234675932703e-06 ... 2.0 0.0038615933242749545 +-16.041177884254317 1.8622469351865406e-05 ... 4.0 0.002999595653967763 + -- 0.0 ... 0.0 -- +-15.878800659290333 9.311234675932703e-06 ... 2.0 0.0013060776318467098 + -- 0.0 ... 0.0 -- +Length = 99 rows +Found vol renormalisation scale of 0.093 +Found d8 renormalisation scale of 1.066 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.72716766679641 0.00024976310426716615 ... 5.0 0.015947328819278253 + -16.65606635089379 0.0001998104834137329 ... 4.0 0.013912916004290763 + -16.54096769997666 0.00024976310426716615 ... 5.0 0.01114373937443421 + -16.44783744420218 0.00024976310426716615 ... 5.0 0.009034487013940686 +-16.346886683136404 0.00029971572512059935 ... 6.0 0.007350956265230744 +-16.268333154532286 9.990524170686645e-05 ... 2.0 0.006036069915045775 + -16.1189057603977 9.990524170686645e-05 ... 2.0 0.0038615933242749545 +-16.041177884254317 0.0001998104834137329 ... 4.0 0.002999595653967763 + -- 0.0 ... 0.0 -- +-15.878800659290333 9.990524170686645e-05 ... 2.0 0.0013060776318467098 + -- 0.0 ... 0.0 -- +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_2_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_3.fits +Found fields: [b'G12'] +Found redshift limits: 0.045 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.093% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_3_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_3_lumfn.fits. +PROCESS CAT FINISHED. +MEDIAN_M PHI_N PHI_N_ERROR PHI_IVMAX PHI_IVMAX_ERROR N V_ON_VMAX +-------- ----- ----------- --------- --------------- --- --------- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + ... ... ... ... ... ... ... + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- +Length = 99 rows +Found vol renormalisation scale of 0.013 +Found d8 renormalisation scale of 4.887 +MEDIAN_M PHI_N PHI_N_ERROR PHI_IVMAX PHI_IVMAX_ERROR N V_ON_VMAX +-------- ----- ----------- --------- --------------- --- --------- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + ... ... ... ... ... ... ... + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- + -- 0.0 0.0 0.0 0.0 0.0 -- +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G12_ddp_n8_d0_3_lumfn.fits +Done. diff --git a/logs/gama_lumfn_d8_G15.log b/logs/gama_lumfn_d8_G15.log index 8a484858..49bef476 100644 --- a/logs/gama_lumfn_d8_G15.log +++ b/logs/gama_lumfn_d8_G15.log @@ -8,3 +8,245 @@ Found redshift limits: 0.039 < z < 0.248 Retrieved area 60.0000 [sq. deg.] WARNING: Found 0.000% with zmax < 0.0 Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_0_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_0_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.730496175545333 1.860565005959082e-05 ... 4.0 0.017001733460170056 + -16.64773627009072 3.2559887604283934e-05 ... 7.0 0.014738668792876348 +-16.573224386892278 1.3954237544693114e-05 ... 3.0 0.012865201125392604 + -16.4307771508108 2.3257062574488523e-05 ... 5.0 0.009827052614482277 +-16.378993603481526 4.651412514897705e-06 ... 1.0 0.008537782481565134 +-16.257965319344848 1.3954237544693114e-05 ... 3.0 0.006825007736907174 +-16.104241714494194 4.651412514897705e-06 ... 1.0 0.004669571909317657 + -16.03331084244899 1.860565005959082e-05 ... 4.0 0.0037822705547863465 +-15.949469150250081 4.651412514897705e-06 ... 1.0 0.0029067378011918246 + -- 0.0 ... 0.0 -- + -15.752584708355 9.30282502979541e-06 ... 2.0 0.001134031516874485 +Length = 99 rows +Found vol renormalisation scale of 0.301 +Found d8 renormalisation scale of -0.915 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.730496175545333 6.188140741889888e-05 ... 4.0 0.017001733460170056 + -16.64773627009072 0.00010829246298307302 ... 7.0 0.014738668792876348 +-16.573224386892278 4.641105556417415e-05 ... 3.0 0.012865201125392604 + -16.4307771508108 7.735175927362359e-05 ... 5.0 0.009827052614482277 +-16.378993603481526 1.547035185472472e-05 ... 1.0 0.008537782481565134 +-16.257965319344848 4.641105556417415e-05 ... 3.0 0.006825007736907174 +-16.104241714494194 1.547035185472472e-05 ... 1.0 0.004669571909317657 + -16.03331084244899 6.188140741889888e-05 ... 4.0 0.0037822705547863465 +-15.949469150250081 1.547035185472472e-05 ... 1.0 0.0029067378011918246 + -- 0.0 ... 0.0 -- + -15.752584708355 3.094070370944944e-05 ... 2.0 0.001134031516874485 +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_0_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_1.fits +Found fields: [b'G15'] +Found redshift limits: 0.039 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.000% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_1_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_1_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.78240138655025 9.301249704840612e-06 ... 2.0 0.016989309435522872 + -16.64367414015333 1.8602499409681225e-05 ... 4.0 0.014507997379673636 +-16.568445025015194 1.8602499409681225e-05 ... 4.0 0.012505038420551472 + -16.45148413490669 2.325312426210153e-05 ... 5.0 0.009530535894777009 + -16.38273104137842 1.8602499409681225e-05 ... 4.0 0.008909475135399181 + -16.23965885231288 1.3951874557260918e-05 ... 3.0 0.006376720420415395 + -16.16852648703361 9.301249704840612e-06 ... 2.0 0.0055442810448458404 +-16.005495811168814 4.650624852420306e-06 ... 1.0 0.0035264513345085223 +-15.980592987817111 9.301249704840612e-06 ... 2.0 0.003229953430606198 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Found vol renormalisation scale of 0.155 +Found d8 renormalisation scale of -0.224 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.78240138655025 6.0116361485250436e-05 ... 2.0 0.016989309435522872 + -16.64367414015333 0.00012023272297050087 ... 4.0 0.014507997379673636 +-16.568445025015194 0.00012023272297050087 ... 4.0 0.012505038420551472 + -16.45148413490669 0.0001502909037131261 ... 5.0 0.009530535894777009 + -16.38273104137842 0.00012023272297050087 ... 4.0 0.008909475135399181 + -16.23965885231288 9.017454222787564e-05 ... 3.0 0.006376720420415395 + -16.16852648703361 6.0116361485250436e-05 ... 2.0 0.0055442810448458404 +-16.005495811168814 3.0058180742625218e-05 ... 1.0 0.0035264513345085223 +-15.980592987817111 6.0116361485250436e-05 ... 2.0 0.003229953430606198 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_1_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_2.fits +Found fields: [b'G15'] +Found redshift limits: 0.040 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.048% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_2_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_2_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- -------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.72682914593935 1.3961344661634786e-05 ... 3.0 0.01659949976911146 + -16.62295928131076 1.3961344661634786e-05 ... 3.0 0.013911390614858188 +-16.528178967831813 2.3268907769391308e-05 ... 5.0 0.011474482075181583 +-16.447969913530688 2.3268907769391308e-05 ... 5.0 0.009894790628060122 + -16.31300356989683 9.307563107756524e-06 ... 2.0 0.007258436644083143 + -16.27461579635706 1.8615126215513048e-05 ... 4.0 0.006783196087707653 +-16.174963513721586 4.653781553878262e-06 ... 1.0 0.005170895205740605 + -16.03945267955754 9.307563107756524e-06 ... 2.0 0.003541071120284103 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Found vol renormalisation scale of 0.093 +Found d8 renormalisation scale of 1.066 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- -------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.72682914593935 0.0001497987712176071 ... 3.0 0.01659949976911146 + -16.62295928131076 0.0001497987712176071 ... 3.0 0.013911390614858188 +-16.528178967831813 0.0002496646186960118 ... 5.0 0.011474482075181583 +-16.447969913530688 0.0002496646186960118 ... 5.0 0.009894790628060122 + -16.31300356989683 9.986584747840473e-05 ... 2.0 0.007258436644083143 + -16.27461579635706 0.00019973169495680947 ... 4.0 0.006783196087707653 +-16.174963513721586 4.993292373920237e-05 ... 1.0 0.005170895205740605 + -16.03945267955754 9.986584747840473e-05 ... 2.0 0.003541071120284103 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_2_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_3.fits +Found fields: [b'G15'] +Found redshift limits: 0.041 < z < 0.246 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.044% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_3_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_3_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... ---- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.765493312056343 0.000214594086542654 ... 45.0 0.01751206888405877 + -16.65713566913222 0.00018121278419157448 ... 38.0 0.014285373227107377 +-16.550078632156833 0.00024320663141500786 ... 51.0 0.011672886256267185 + -16.4548449496194 0.00023366911645755659 ... 49.0 0.009111652498541286 +-16.368201667144547 0.00022890035897883093 ... 48.0 0.007953749328284573 + -16.27101051348568 6.199384722343338e-05 ... 13.0 0.006332521785004211 + -16.17242433919671 1.4306272436176933e-05 ... 3.0 0.004941810115297417 + -16.03645377362567 1.4306272436176933e-05 ... 3.0 0.003149130892600848 +-15.947351778938028 3.338130235107951e-05 ... 7.0 0.0022843398819084792 + -15.81915174837022 1.4306272436176933e-05 ... 3.0 0.0010191612882905394 +-15.774172226965995 4.768757478725644e-06 ... 1.0 0.0006302898890691272 +Length = 99 rows +Found vol renormalisation scale of 0.013 +Found d8 renormalisation scale of 4.887 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... ---- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.765493312056343 0.016178281439687435 ... 45.0 0.01751206888405877 + -16.65713566913222 0.01366165988240272 ... 38.0 0.014285373227107377 +-16.550078632156833 0.018335385631645756 ... 51.0 0.011672886256267185 + -16.4548449496194 0.017616350900992984 ... 49.0 0.009111652498541286 +-16.368201667144547 0.017256833535666594 ... 48.0 0.007953749328284573 + -16.27101051348568 0.004673725749243036 ... 13.0 0.006332521785004211 + -16.17242433919671 0.001078552095979162 ... 3.0 0.004941810115297417 + -16.03645377362567 0.001078552095979162 ... 3.0 0.003149130892600848 +-15.947351778938028 0.0025166215572847116 ... 7.0 0.0022843398819084792 + -15.81915174837022 0.001078552095979162 ... 3.0 0.0010191612882905394 +-15.774172226965995 0.00035951736532638737 ... 1.0 0.0006302898890691272 +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G15_ddp_n8_d0_3_lumfn.fits +Done. diff --git a/logs/gama_lumfn_d8_G9.log b/logs/gama_lumfn_d8_G9.log index 2609e614..a6db0ee5 100644 --- a/logs/gama_lumfn_d8_G9.log +++ b/logs/gama_lumfn_d8_G9.log @@ -8,3 +8,245 @@ Found redshift limits: 0.043 < z < 0.248 Retrieved area 60.0000 [sq. deg.] WARNING: Found 0.000% with zmax < 0.0 Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_0_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_0_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.757530035162436 2.3336539326931244e-05 ... 5.0 0.01631836015638575 +-16.631379546591454 1.4001923596158748e-05 ... 3.0 0.012822800568102252 + -16.52807614224665 4.667307865386249e-06 ... 1.0 0.010210696092201984 +-16.417128789933166 1.4001923596158748e-05 ... 3.0 0.007607699242488683 +-16.300399356494676 4.667307865386249e-06 ... 1.0 0.005902233950808325 +-16.252453820268638 4.667307865386249e-06 ... 1.0 0.005178948636583736 + -- 0.0 ... 0.0 -- + -16.05538488090431 4.667307865386249e-06 ... 1.0 0.0017669671728029594 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Found vol renormalisation scale of 0.301 +Found d8 renormalisation scale of -0.915 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- --------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.757530035162436 7.761609474604659e-05 ... 5.0 0.01631836015638575 +-16.631379546591454 4.656965684762796e-05 ... 3.0 0.012822800568102252 + -16.52807614224665 1.552321894920932e-05 ... 1.0 0.010210696092201984 +-16.417128789933166 4.656965684762796e-05 ... 3.0 0.007607699242488683 +-16.300399356494676 1.552321894920932e-05 ... 1.0 0.005902233950808325 +-16.252453820268638 1.552321894920932e-05 ... 1.0 0.005178948636583736 + -- 0.0 ... 0.0 -- + -16.05538488090431 1.552321894920932e-05 ... 1.0 0.0017669671728029594 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_0_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_1.fits +Found fields: [b'G9'] +Found redshift limits: 0.039 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.160% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_1_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_1_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.769238354145056 4.653279074321804e-06 ... 1.0 0.01817428089201988 + -16.69397717036171 1.3959837222965412e-05 ... 3.0 0.01593443343570443 +-16.562882007537468 4.653279074321804e-06 ... 1.0 0.012725980562386057 +-16.458402460416913 9.306558148643608e-06 ... 2.0 0.010462018462594154 + -16.39486591203184 4.653279074321804e-06 ... 1.0 0.009192218863322043 + -16.21355277376218 4.653279074321804e-06 ... 1.0 0.006066828227062017 +-16.157693863414032 4.653279074321804e-06 ... 1.0 0.005327476934110926 +-16.067487161907557 4.653279074321804e-06 ... 1.0 0.0042649396729809354 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Found vol renormalisation scale of 0.155 +Found d8 renormalisation scale of -0.224 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- --------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... +-16.769238354145056 3.007533565926072e-05 ... 1.0 0.01817428089201988 + -16.69397717036171 9.022600697778216e-05 ... 3.0 0.01593443343570443 +-16.562882007537468 3.007533565926072e-05 ... 1.0 0.012725980562386057 +-16.458402460416913 6.015067131852144e-05 ... 2.0 0.010462018462594154 + -16.39486591203184 3.007533565926072e-05 ... 1.0 0.009192218863322043 + -16.21355277376218 3.007533565926072e-05 ... 1.0 0.006066828227062017 +-16.157693863414032 3.007533565926072e-05 ... 1.0 0.005327476934110926 +-16.067487161907557 3.007533565926072e-05 ... 1.0 0.0042649396729809354 + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_1_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_2.fits +Found fields: [b'G9'] +Found redshift limits: 0.039 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.056% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_2_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_2_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.74440852995106 1.3951874557260918e-05 ... 3.0 0.01757321265286468 + -16.65933113704299 1.8602499409681225e-05 ... 4.0 0.015058953483199607 + -16.52443036950733 1.8602499409681225e-05 ... 4.0 0.011576659630882672 + -16.47012166970818 2.325312426210153e-05 ... 5.0 0.01062266391979281 + -16.35732955076115 1.3951874557260918e-05 ... 3.0 0.008465191812360285 +-16.261091627218327 9.301249704840612e-06 ... 2.0 0.006279696963337101 +-16.135172685292325 1.8602499409681225e-05 ... 4.0 0.005076216295062213 +-16.036449159684274 4.650624852420306e-06 ... 1.0 0.0038713716255760824 +-15.948274125651977 1.8602499409681225e-05 ... 4.0 0.0027615766564690436 + -- 0.0 ... 0.0 -- +-15.714032676721347 4.650624852420306e-06 ... 1.0 0.0009354002067586114 +Length = 99 rows +Found vol renormalisation scale of 0.093 +Found d8 renormalisation scale of 1.066 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... --- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.74440852995106 0.00014969716137751724 ... 3.0 0.01757321265286468 + -16.65933113704299 0.000199596215170023 ... 4.0 0.015058953483199607 + -16.52443036950733 0.000199596215170023 ... 4.0 0.011576659630882672 + -16.47012166970818 0.00024949526896252876 ... 5.0 0.01062266391979281 + -16.35732955076115 0.00014969716137751724 ... 3.0 0.008465191812360285 +-16.261091627218327 9.97981075850115e-05 ... 2.0 0.006279696963337101 +-16.135172685292325 0.000199596215170023 ... 4.0 0.005076216295062213 +-16.036449159684274 4.989905379250575e-05 ... 1.0 0.0038713716255760824 +-15.948274125651977 0.000199596215170023 ... 4.0 0.0027615766564690436 + -- 0.0 ... 0.0 -- +-15.714032676721347 4.989905379250575e-05 ... 1.0 0.0009354002067586114 +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_2_lumfn.fits + +Reading: /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_3.fits +Found fields: [b'G9'] +Found redshift limits: 0.040 < z < 0.248 +Retrieved area 60.0000 [sq. deg.] +WARNING: Found 0.117% with zmax < 0.0 +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_3_vmax.fits. +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_3_lumfn.fits. +PROCESS CAT FINISHED. + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... ---- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.73814624237567 0.00012561827536222467 ... 27.0 0.016440356987521974 + -16.64898133689124 0.00015353344766494127 ... 33.0 0.014233359651233603 +-16.557272790775766 0.00015353344766494127 ... 33.0 0.011748184024985437 +-16.457533287555453 0.000172143562533419 ... 37.0 0.009856827436396022 + -16.35542173060181 0.00015353344766494127 ... 33.0 0.00787997781426617 +-16.258053943921098 0.00012096574664510523 ... 26.0 0.006359474861188578 + -16.1736565311367 3.256770101983602e-05 ... 7.0 0.005164046188948062 + -16.04352617198227 1.3957586151358296e-05 ... 3.0 0.0035466618054241653 + -- 0.0 ... 0.0 -- +-15.894236113006684 4.6525287171194314e-06 ... 1.0 0.002015023201173478 +-15.771159111866156 1.3957586151358296e-05 ... 3.0 0.0008969211010103589 +Length = 99 rows +Found vol renormalisation scale of 0.013 +Found d8 renormalisation scale of 4.887 + MEDIAN_M PHI_N ... N V_ON_VMAX +------------------- ---------------------- ... ---- --------------------- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + -- 0.0 ... 0.0 -- + ... ... ... ... ... + -16.73814624237567 0.009470381246382923 ... 27.0 0.016440356987521974 + -16.64898133689124 0.011574910412245796 ... 33.0 0.014233359651233603 +-16.557272790775766 0.011574910412245796 ... 33.0 0.011748184024985437 +-16.457533287555453 0.012977929856154377 ... 37.0 0.009856827436396022 + -16.35542173060181 0.011574910412245796 ... 33.0 0.00787997781426617 +-16.258053943921098 0.009119626385405777 ... 26.0 0.006359474861188578 + -16.1736565311367 0.0024552840268400166 ... 7.0 0.005164046188948062 + -16.04352617198227 0.001052264582931436 ... 3.0 0.0035466618054241653 + -- 0.0 ... 0.0 -- +-15.894236113006684 0.00035075486097714526 ... 1.0 0.002015023201173478 +-15.771159111866156 0.001052264582931436 ... 3.0 0.0008969211010103589 +Length = 99 rows +Writing /cosma5/data/durham/dc-wils7/GAMA4//gama_gold_G9_ddp_n8_d0_3_lumfn.fits +Done. diff --git a/lumfn.py b/lumfn.py index e02ad0dd..36ca2fa6 100644 --- a/lumfn.py +++ b/lumfn.py @@ -4,15 +4,17 @@ from cosmo import volcom -def lumfn(dat, vol, Ms=np.arange(-25.5, -15.5, 0.1), Mcol='MCOLOR_0P0'): +def lumfn(dat, Ms=np.arange(-25.5, -15.5, 0.1), Mcol='MCOLOR_0P0'): idxs = np.digitize(dat[Mcol], bins=Ms) result = [] - ds = np.diff(Ms) - dM = ds[0] + ds = np.diff(Ms) + dM = ds[0] assert np.all(ds == dM) + vol = dat.meta['VOLUME'] + for idx in np.arange(len(Ms) - 1): sample = dat[idxs == idx] nsample = len(sample) @@ -35,7 +37,7 @@ def lumfn(dat, vol, Ms=np.arange(-25.5, -15.5, 0.1), Mcol='MCOLOR_0P0'): result = Table(np.array(result), names=names) result.meta.update(dat.meta) - result.meta['MS'] = Ms.tolist() + result.meta['MS'] = str(Ms.tolist()) result.meta['VOLUME'] = vol return result