From dd374aab24eb9490387d04cf7fd8b52331292cc2 Mon Sep 17 00:00:00 2001 From: Alexander Mendez Date: Fri, 17 Jan 2014 11:09:28 -0800 Subject: [PATCH 1/2] New examples directory that contains isedfit example and setup checks --- examples/isedfit/cosmos_example.fits.gz | Bin 0 -> 1881 bytes examples/isedfit/example.bashrc | 25 +++++++ examples/isedfit/isedfit_buildexample.pro | 27 +++++++ examples/isedfit/isedfit_check.pro | 26 +++++++ examples/isedfit/isedfit_quick.pro | 84 ++++++++++++++++++++++ 5 files changed, 162 insertions(+) create mode 100644 examples/isedfit/cosmos_example.fits.gz create mode 100755 examples/isedfit/example.bashrc create mode 100644 examples/isedfit/isedfit_buildexample.pro create mode 100644 examples/isedfit/isedfit_check.pro create mode 100644 examples/isedfit/isedfit_quick.pro diff --git a/examples/isedfit/cosmos_example.fits.gz b/examples/isedfit/cosmos_example.fits.gz new file mode 100644 index 0000000000000000000000000000000000000000..9211e7ded0bd49e7b1e775b06968977815b088d2 GIT binary patch literal 1881 zcmV-f2d4NRiwFn}c-c|_17mMs8@Bc)?f4-8NJf;P9#;Uw+UcdiZI%w_eb~pPqO>G``PE+@3Xtl zzC{{sK|X;XPaz}7fS@Y3&E`RJcow=?6GD3s?K6eB8bcAX3`k`o9Fnau6lfOSfhQZ{ zmGYbo@*c#(LKDj)%=hs_(kT~e0$U^UFa*^(5EFDkgmWU$>2xq$5LC{&EQre@>>ND9 zZ4nC%gidwm`4SZ_*6!pSb|+$k9wfe;vsz~-AT%rF%sN@SYuF({m?%-AM2Qk5N|dOd zi~{Xskw&k33W>|s=nTs2d?M~!VKkew2O3=$q^&Gc&T}rIy&SEL9;<&YWXs$fFSNMJ zpw-Qeo!_YaGk%XeFd{ECJzStDs1qM^3i$t$oYZq@C1r;?EDIH z^b2%`(Y_aW8?;(I!vA=?Lzp%1;p}Dv@;vl76TW|ldNdxm`|D@O)e!lrqLF#&-<~hy zu)1x+X9#B>tez2-qt|MQK!PFCZC+wxB0|_rflO@BE_`D-gfcGSSgXkyjNL6S)cJ*) z7noB{@RB2TXk+hcl3~oxU|3W@1n;-FLgm1&?>t}bi5I-z;#5RVhNX zREUE|%)u2zK4Hun4|zSSGy!`9*DpQ=PeK4l_(9k2r?#G`px^7CeL&7C#{F&crF7}q ztoeJsW?(6No@7caF4D%DVPLVCC7saB+DJ}X3i(mu?r-8wHw}5W@Dq$P&_vWF&{O{1&&0!x^j=`?`omgG< z7MOhYBZfHE#>}oZ(a$e#r2W~e>HAl&G@;!Xv*~RsMK|_Q$KKMx>C!~uzTAeMSlOeH zboD$VBOB#=&R2I23!d4q}auO~qcyypEllk%unrV6do3#aP=DJ0b=J zRpTcyJo64S&oqawICGQUzT{W*eLldPbSi#Git!hN7ki4%o7B@G_ z|GdT_dKq3XLC*0+OzGBRT{+i4-pL{`>xv3I8{2^ri*KRp4}J!;t6vLvgY8LunA$c| zO&>i$(-Tcsp)RTnZnS*~Y?7XP2ITD)%>f=N& z!|bI$iNFdnd{|50n|MT*AC$jh$AA9WcBUH<^2oOdds6p?U3sm>h z4NT3|XlBlNDXg7xjV66bFu_~Gobom>+ACF5YEh=KtiOSje7;ri5^2h~zEL5uk*6$Q zMP^()UpR&`KWNkMdm|zsVlrZJKc?!wO zO}ibhACRBB=oh^VvzM0N)nJvzsaVsiI{^R5WgvfV7pSgXg&OBg#OBUoLFJ`L@bY*W zu>N5hRa57sHyn?km5-;w*iU;P{c$;Uy113GUC=NOL~F=1tCGmJsgY#=A&1cW^u1WT z{jfr^`M6tFTPR=m$m^n)VfOOVLyuw~UTy^`Z+C(-Di8Smwi~!4!H530y$;LIn*(~` zIB check_isedfit +; - + +PRO _req, cmd, name, version + ;; a little wrapper that runs a command to ensure that the installed + ;; version matches the required version passed into the function. + value = call_function(cmd) + status = (value EQ version) ? 'ok' : 'Failed!' + print, 'iSEDfit requires '+version+' version of '+name+'. Found: '+value+'. '+status + + IF value NE version THEN BEGIN + message, ' !!! Version did not match the required version, please update' + ENDIF +END + + +PRO isedfit_check + ;; might also want to check directories ($ISEDFIT_DIR, + ;; $KCORRECT_DIR) as well. + _req, 'idlutils_version', 'IDLUtils', 'trunk' + _req, 'k_version', 'kCorrect', 'NOCVS:kcorrect' +END + diff --git a/examples/isedfit/isedfit_quick.pro b/examples/isedfit/isedfit_quick.pro new file mode 100644 index 0000000..f2a3292 --- /dev/null +++ b/examples/isedfit/isedfit_quick.pro @@ -0,0 +1,84 @@ +;+ +; NAME: +; ISEDFIT_QUICK +; +; PURPOSE: +; Quick start script that goes through the process. +; Overwrites previous run results +; This is just to test that you can run isedfit quickly. +;- + +PRO isedfit_quick, extras=extras + + ;; output directory -- depending on parameters may be large [few Gb] + prefix = 'example' + + isedfit_dir = '/tmp/isedfit_example/' + montegrids_dir = isedfit_dir+'montegrids/' + file_mkdir, montegrids_dir + + isedfit_paramfile = isedfit_dir+prefix+'_paramfile.par' + + ; -------------------------------------------------- + ;; catalog to be processed + example = im_mrdfits('cosmos_example.fits.gz') + filterlist = strtrim(example[0].filterlist, 2) + + + ; -------------------------------------------------- + ; choose your priors: write the iSEDfit parameter file + write_isedfit_paramfile, params=params, isedfit_dir=isedfit_dir, $ + prefix=prefix, filterlist=filterlist, spsmodels='fsps_v2.4_miles', $ + imf='chab', redcurve='charlot', /igm, zminmax=[0.1, 1.3], zbin=0.05, $ + nmodel=500L, age=[0.1, 13.0], tau=[0.1, 5.0], Zmetal=[0.004, 0.03], $ + AV=[0.35, 2.0], mu=[0.1, 4.0], pburst=0.2, interval_pburst=2.0, $ + tburst=[0.1, 13.0], /delayed, galchunksize=250L, /clobber + + ; -------------------------------------------------- + ; build the Monte Carlo grids [time to go read astro-ph] + isedfit_montegrids, isedfit_paramfile, isedfit_dir=isedfit_dir, $ + montegrids_dir=montegrids_dir, /clobber + + ; -------------------------------------------------- + ; calculate the model photometry + isedfit_models, isedfit_paramfile, isedfit_dir=isedfit_dir, $ + montegrids_dir=montegrids_dir, /clobber + + ; -------------------------------------------------- + ; fit the catalog + isedfit, isedfit_paramfile, example.maggies, example.ivarmaggies, $ + example.z, ra=example.ra, dec=example.dec, $ + isedfit_dir=isedfit_dir, /clobber, $ + isedfit_results=mass + + ;; simple plot + plot, mass.mstar_50, mass.sfr100_50, psym=6, $ + xtitle='Mass [M_sun]', ytitle='SFR [M_sun / yr]' + + + print, 'Huzza! iSEDfit finished.' + print, ' Check the plot, and make sure that these sfr/masses make sense' + + + ;; Here be dragons... + + IF keyword_set(extras) THEN BEGIN + ;; Get some fancy plots [beta] that show the properties of the + ;; models and the sort. + + isedfit_qaplot_models, isedfit_paramfile, example.maggies, $ + example.ivarmaggies, example.z, isedfit_dir=isedfit_dir, $ + thesefilters=filterlist, /clobber + + isedfit_kcorrect, isedfit_paramfile, isedfit_dir=isedfit_dir, $ + montegrids_dir=montegrids_dir, /clobber + + isedfit_qaplot_sed, isedfit_paramfile, nrandom=50, $ + isedfit_dir=isedfit_dir, montegrids_dir=montegrids_dir, $ + /xlog, /clobber + + ENDIF + + + +END From 415a09e2594fa979d7829889054b26934fbae44b Mon Sep 17 00:00:00 2001 From: Alexander Mendez Date: Fri, 17 Jan 2014 11:14:03 -0800 Subject: [PATCH 2/2] ssp directory did not make it into the commit --- examples/isedfit/example.bashrc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/examples/isedfit/example.bashrc b/examples/isedfit/example.bashrc index cee4e3f..8a21c20 100755 --- a/examples/isedfit/example.bashrc +++ b/examples/isedfit/example.bashrc @@ -22,4 +22,4 @@ export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$KCORRECT_DIR/lib # IMPRO / ISEDFIT export IMPRO_DIR=$HOME/idl/impro export IDL_PATH=$IDL_PATH:+$IMPRO_DIR - +export ISEDFIT_SSP_DIR=$HOME/raid/isedfit2/isedfit_ssp