Merge pull request #13 from lephare-photoz/update-registry

[Auto] Update registry file 30-04-2024

data_registry.txt

@@ -1,16 +1,17 @@
 alloutputkeys.txt bfb021b6750479edec7d258e213d3a9959c96a4ee4800687b43df93697d8a3e6
 examples/AGE_BC03COMB.dat 84bad53c8c00baeeafdf62e557962e6cc6112f54187aa8753d08fdfcacf83113
 examples/AGE_MASATO_BC03.dat f4fe3011a1be9a5d76016699ad07ccae3bd4bb511bd6ffc50ed261b7ab8ce209
+examples/AGN_MOD.list f4acb3dcc1205261ac9ceabac18de758a588c16a122d59e7972272e2fc224f2b
 examples/BC03COMB_MOD.list 33a856eeb9090b252da8cc36d819781ca025b57b20dafc487d85dbe6738ae211
 examples/COSMOS.in 197096c8eed2a5448c742ed872d092d47ceba2b4ae44f75f99852f344544f688
 examples/COSMOS.in.header 45b3ab33340aaef61fcf9e4a17398794df6b060815220d8f16d5896ee129624b
-examples/COSMOS.para f195bf8730f816d104d3e46088aa9ffd37bbccdb29032d253371396388698cbe
+examples/COSMOS.para 4328d7f33e3d0b68d553480cd46f76525c70e235a631375bcc31a6d217f22db6
 examples/COSMOS_MOD.list f78aa98c20356943efbba7e685775e0708894d3b01b539df76948f258ceadf2a
 examples/DALE.list 6d92ea8b28d8b99f0685e46d417cf3fe67de7e10cf4870b123a82d63546ef03d
 examples/HERSCHEL.para 0e05f04ad8efda225f8dc6237780e0d4f86d6af89b9af8401d7275be2afdd438
 examples/MIPS.in 0a6be5e20ce9e40fa86854661707e0341621aa5785dd2f6ea65f641e15e7186f
 examples/MIPS.in.head 933617a5896275cb4257c0f4902ab4392c16281119529649f76446b99af41a2e
-examples/README_full a55c33bf1b19c725d918e712e1c82ba2745a97f8f1fd11e7a81c6650abe34be3
+examples/README_full 87aa3a118117154ef192a34315bf85ae81613cb09abaa1ea84bfed5b8bb6a403
 examples/README_short 7534e4d6c6f0aa4ffa500fe7fb5eb423e40ade54171917dd5db8f681fa434fce
 examples/STAR_MOD_ALL.list 13d1a023e7c31638c58c6228ce4606c12962b128294fd1e37b9b4c5cb67c4c64
 examples/UVISTA_MOD.list a0070b6f7440f03a723c0a26ff6cf998eb062f6fc6d119db4c57fbaad3d78fdd
@@ -32,7 +33,7 @@ ext/LMC_Fitzpatrick.dat 1d3f7729426c3e9228e242d0f97b5bb89d64c6d3ee48908e86e538c2
 ext/LMC_Fitzpatrick_noBump.dat 93b1c42d34a41e6107c23769b715343fe0f12793dfc4f92a312183e455835ea8
 ext/MW_Allen.dat c153ffb512231a03ff6fd6866ffc4bf544eee6ab49e426b6bae9830733db78ca
 ext/MW_seaton.dat 04233459f6a1b465098378338b86e78c165f00c9d965a57620e4a2733d329f15
-ext/README.extinc d0004adcc79d8bed5a26ceb1fd56c78c58bee5aa536357cd42147b1ff4d3edf2
+ext/README.extinc f1a1d43f4ab9559feb8d10e89f7b9c404174e2fa98282efa79e90176266e704e
 ext/SB_calzetti.dat 43ea0474f2ed66e20a44482782e04e91fcbc72cf5acd6045c22413fcf96bb890
 ext/SB_calzetti_bump1.dat 20c281833c6487116de5fbd6e49d087f03badc738c028c43791c52a44cf82120
 ext/SB_calzetti_bump2.dat f26874d5acbe8b53ee8b177965b052b807e67c373f70900fe59553860f8220fa
@@ -48,6 +49,7 @@ ext/extinc_ctio.dat 56da1af053d2594f7d85df53cf61351fc47c170ab911720c9386dde45613
 ext/extinc_etc.dat cb2f4a6e26083ad0110f5369fdf30d136ec1eb7da56eb35149d2dad4839fd8f6
 ext/extinction_laws.ps 52cd64d51be561740bd14b63025351e43176281adf4a84c4ecdbe4637fedda49
 ext/extlaw_0.9.dat b3414a2875ca9f584fe4c98c125f76ebb99598b4a42cfb2400a4f726a912d1e1
+ext/salim18.dat 18130b316208f10acf611ea3bca67259797a9378c78bcb926cbeda6f6643705b
 filt/2mass/H.pb 650367d4341e0406580fd8f7fc05fa3e1ddd403baea22513b51238449d6348ff
 filt/2mass/J.pb 9a53ed73c8df7e65fb5172b1809e62c674b8b5d2e6f178d303e4f7c976d681e5
 filt/2mass/Ks.pb 4360bd9bb7bd0c7c3ce0ed5100f737616ad9715b753a2b98884d95f211edf313

data_registry_hash.sha256

@@ -1 +1 @@
-8aa0fd98fbd133dad56fb069a7487fe4a1e2576796933b68fab93c9bbe6d3950
+8c08b496a0969a393911a6337f2c308fdba2d388420b3a669064374e8c24350c

examples/AGN_MOD.list

@@ -0,0 +1,14 @@
+SALVATO09/BQSO1_template_norm.sed A
+SALVATO09/I22491_template_norm.sed A
+SALVATO09/M82_template_norm.sed A
+SALVATO09/Mrk231_template_norm.sed A
+SALVATO09/QSO1_template_norm.sed A
+SALVATO09/QSO2_template_norm.sed A
+SALVATO09/Sey18_template_norm.sed A
+SALVATO09/Sey2_template_norm.sed A
+SALVATO09/Spi4_template_norm.sed A
+SALVATO09/TQSO1_template_norm.sed A
+SALVATO09/hybrid_gal20_agn80.sed A
+SALVATO09/hybrid_gal50_agn50.sed A
+SALVATO09/hybrid_gal80_agn20.sed A
+SALVATO09/QSO_High_lum_DR2-032_t0_ext.spec A

examples/COSMOS.para

@@ -102,7 +102,7 @@ ZPHOTLIB	CE_COSMOS,STAR_COSMOS,QSO_COSMOS  # Library used for Chi2 (max:3)
 #------------------     Output catalog Name           ---------------------------------
 CAT_OUT		zphot.out 	# Output catalog  (-> working directory) 
 PARA_OUT        $LEPHAREDIR/examples/output.para     # Ouput parameters (full path) 
-VERBOSE         NO             # Display the templates/objects which are computed
+VERBOSE         YES             # Display the templates/objects which are computed
 PDZ_TYPE        BAY_ZG           # PDZ in output [def-BAY]. BAY_ZG sum all probabilities at a given z. MIN_ZG takes exp(-chi2_min/2) at a given z, with chi2_min the minimum chi2 at this z.
                                  # You can add and combine BAY_ZG,BAY_ZQ,MIN_ZG,MIN_ZQ,MASS,SFR,SSFR,AGE
 PDZ_OUT         test             # root of the pdz output files [def-NONE]

examples/README_full

@@ -1,54 +1,43 @@
 #! /bin/csh
 
-
-# Location of the code
-export LEPHAREDIR='/home/oilbert/LEPHARE/'
-# Where to store the libraries (if not created with the make file, you must have a $LEPHAREWORK/filt, $LEPHAREWORK/lib_bin, $LEPHAREWORK/lib_mag
+# Indicate the directory containing all the LePHARE data, as the templates/filters/attenutation
+export LEPHAREDIR='/Users/ilbert/LEPHARE-data/'
+# Where to store the libraries created by filter/sedtolib/mag_gal
 export LEPHAREWORK='WORK'
-# Shift to apply to the photometry (systematic offsets) for zphota
-export SHIFT="0.051121,-0.016466,-0.084605,-0.058664,-0.039816,-0.053934,-0.072388,-0.023517,-0.001748,0.013681,-0.036044,0.041802,0.017099,0.006699,0.006539,0.021207,0.063196,0.015433,-0.163286,0.032249,0.039840,0.050121,0.025781,0.012588,0.052423,0.034833,0.018244,-0.024611,0.000000,0.000000"
-# Maximum number of threads
-export OMP_NUM_THREADS='25' 
+# Maximum number of threads. Don't put more than available. It could result in a strange behaviour of the code.
+export OMP_NUM_THREADS='10' 
 
 
 #############
 ##  LIBRARY
 
 # Read the filters and compile them within one file stored in $WORK/filt
-$LEPHAREDIR/bin/filter -c COSMOS.para
-# read a second set of filter which could be used to compute predicted magnitudes in some bands not present in the dataset (useful for simulation)
-$LEPHAREDIR/bin/filter -c COSMOS.para  -FILTER_LIST euclid/EUC_riz.lowres,euclid/EUC_Y.lowres,euclid/EUC_J.lowres,euclid/EUC_H.lowres -FILTER_CALIB 0 -FILTER_FILE filter_add
+filter -c COSMOS.para
+# Read a second set of filters which could be used to compute predicted magnitudes in some bands not present in the dataset (useful to simulate data)
+filter -c COSMOS.para  --FILTER_LIST euclid/EUC_riz.lowres,euclid/EUC_Y.lowres,euclid/EUC_J.lowres,euclid/EUC_H.lowres --FILTER_CALIB 0 --FILTER_FILE filter_add
 
-# Case of the galaxy templates used in Ilbert 2013
+# Generate the library of predicted magnitudes for galaxy templates
+# Case of the COSMOS templates used in Ilbert 2013
 # Read the galaxy templates and store them in $WORK/lib_bin
-$LEPHAREDIR/bin/sedtolib -c COSMOS.para -t G -GAL_SED COSMOS_MOD.list  -GAL_LIB LIB_VISTA 
+sedtolib -c COSMOS.para -t G --GAL_SED COSMOS_MOD.list  --GAL_LIB LIB_VISTA 
 # Use the galaxy templates + filters to derive a library of predicted magnitudes and store in $WORK/lib_mag
-# Case emission lines correlated to UV light + free factor in scaling the lines 
-$LEPHAREDIR/bin/mag_gal  -c COSMOS.para -t G -GAL_LIB_IN LIB_VISTA -GAL_LIB_OUT VISTA_COSMOS_FREE -MOD_EXTINC 18,26,26,33,26,33,26,33  -EXTINC_LAW SMC_prevot.dat,SB_calzetti.dat,SB_calzetti_bump1.dat,SB_calzetti_bump2.dat  -EM_LINES EMP_UV  -EM_DISPERSION 0.5,0.75,1.,1.5,2.
-# Case emission lines correlated to UV light + no free factor in scaling the lines 
-$LEPHAREDIR/bin/mag_gal  -c COSMOS.para -t G -GAL_LIB_IN LIB_VISTA -GAL_LIB_OUT VISTA_COSMOS -MOD_EXTINC 18,26,26,33,26,33,26,33  -EXTINC_LAW SMC_prevot.dat,SB_calzetti.dat,SB_calzetti_bump1.dat,SB_calzetti_bump2.dat -EM_LINES EMP_UV  -EM_DISPERSION 1.
+# Generate emission lines correlated to UV light + free factor in scaling the lines by a factor 2 
+mag_gal  -c COSMOS.para -t G --GAL_LIB_IN LIB_VISTA --GAL_LIB_OUT VISTA_COSMOS_FREE --MOD_EXTINC 18,26,26,33,26,33,26,33  --EXTINC_LAW SMC_prevot.dat,SB_calzetti.dat,SB_calzetti_bump1.dat,SB_calzetti_bump2.dat  --EM_LINES EMP_UV  --EM_DISPERSION 0.5,0.75,1.,1.5,2.
+# Generate emission lines correlated to UV light + no free factor in scaling the lines 
+mag_gal  -c COSMOS.para -t G --GAL_LIB_IN LIB_VISTA --GAL_LIB_OUT VISTA_COSMOS --MOD_EXTINC 18,26,26,33,26,33,26,33  --EXTINC_LAW SMC_prevot.dat,SB_calzetti.dat,SB_calzetti_bump1.dat,SB_calzetti_bump2.dat --EM_LINES EMP_UV  --EM_DISPERSION 1.
 
-# STAR
+# Generate the library of predicted magnitudes for stars
 # Read the stellar templates and store them in $WORK/lib_bin
-$LEPHAREDIR/bin/sedtolib -c COSMOS.para -t S -STAR_SED STAR_MOD_ALL.list
+sedtolib -c COSMOS.para -t S --STAR_SED STAR_MOD_ALL.list
 # Use the stellar templates + filters to derive a library of predicted magnitudes and store in $WORK/lib_mag
-$LEPHAREDIR/bin/mag_gal -c COSMOS.para -t S -LIB_ASCII YES -STAR_LIB_OUT ALLSTAR_COSMOS 
+mag_gal -c COSMOS.para -t S --LIB_ASCII YES --STAR_LIB_OUT ALLSTAR_COSMOS 
 
+# Generate the library of predicted magnitudes for AGN templates
 # AGN from Salvato 2009
 # Read the AGN templates and store them in $WORK/lib_bin
-$LEPHAREDIR/bin/sedtolib -c COSMOS.para -t Q -QSO_SED  $LEPHAREDIR/sed/QSO/SALVATO09/AGN_MOD.list
+sedtolib -c COSMOS.para -t Q --QSO_SED  AGN_MOD.list
 # Use the AGN templates + filters to derive a library of predicted magnitudes and store in $WORK/lib_mag
-$LEPHAREDIR/bin/mag_gal -c COSMOS.para -t Q -MOD_EXTINC 0,1000  -EB_V 0.,0.1,0.2,0.3 -EXTINC_LAW SB_calzetti.dat -LIB_ASCII NO
-
-
-# BC03 library as Ilbert+2015
-# Library to derive physical parameters
-# Read the galaxy BC03 templates and store them in $WORK/lib_bin
-$LEPHAREDIR/bin/sedtolib -c COSMOS.para -t G -GAL_SED BC03COMB_MOD.list -GAL_LIB LIB_COMB -SEL_AGE AGE_BC03COMB.dat
-# Use the galaxy BC03 templates + filters to derive a library of predicted magnitudes and store in $WORK/lib_mag
-# Use the "physical" recipes to generate the lines
-# Add the dust component with one template of Bethermin+12
-$LEPHAREDIR/bin/mag_gal  -c COSMOS.para -t G -GAL_LIB_IN LIB_COMB  -GAL_LIB_OUT BC03_COSMOS  -MOD_EXTINC 3,12,3,12 -EXTINC_LAW SB_calzetti.dat,extlaw_0.9.dat -EB_V 0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7 -Z_STEP 0.03,0,6 -EM_LINES PHYS -EM_DISPERSION 0.25,0.5,1,2,4 -ADD_DUSTEM YES
+mag_gal -c COSMOS.para -t Q --MOD_EXTINC 0,1000  --EB_V 0.,0.1,0.2,0.3 --EXTINC_LAW SB_calzetti.dat --LIB_ASCII NO
 
 
 
@@ -58,61 +47,80 @@ $LEPHAREDIR/bin/mag_gal  -c COSMOS.para -t G -GAL_LIB_IN LIB_COMB  -GAL_LIB_OUT
 
 
 # Compute the photo-z and derive the adapation of the zero-points
-$LEPHAREDIR/bin/zphota -c COSMOS.para -CAT_IN COSMOS.in -CAT_OUT zphot_vista_adapt.out -ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS  -ADD_EMLINES 0,100 -AUTO_ADAPT YES  
-# Compute the photo-z without the adapation of the zero-points and applying the shifts given in keyword
-$LEPHAREDIR/bin/zphota -c COSMOS.para -CAT_IN COSMOS.in  -CAT_OUT zphot_vista_shift.out -ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS   -ADD_EMLINES 0,100 -AUTO_ADAPT NO -APPLY_SYSSHIFT $SHIFT -PDZ_OUT zphot_vista_shift -PDZ_TYPE BAY_ZG
-# No emission lines
-$LEPHAREDIR/bin/zphota -c COSMOS.para -CAT_IN COSMOS.in  -CAT_OUT zphot_vista_noline.out -ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS  -ADD_EMLINES 0,0 -AUTO_ADAPT NO -APPLY_SYSSHIFT $SHIFT 
-# No shifts
-$LEPHAREDIR/bin/zphota -c COSMOS.para -CAT_IN COSMOS.in  -CAT_OUT zphot_vista_noadapt.out -ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS   -ADD_EMLINES 0,100 -AUTO_ADAPT NO 
+zphota -c COSMOS.para --CAT_IN COSMOS.in --CAT_OUT zphot_vista_adapt.out --ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS  --ADD_EMLINES 0,100 --AUTO_ADAPT YES  
 
-# Extract Id*.spec and Id*.chi for few sources
-$LEPHAREDIR/bin/zphota -c COSMOS.para -CAT_IN COSMOS.in  -CAT_OUT zphot_vista_select.out -ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS   -ADD_EMLINES 0,100 -AUTO_ADAPT NO -APPLY_SYSSHIFT $SHIFT -PDZ_OUT NONE -PDZ_TYPE BAY_ZG -CAT_LINES 12000,12010 -SPEC_OUT YES 
 
-# Derive the stellar masses
-$LEPHAREDIR/bin/zphota -c COSMOS.para -CAT_IN COSMOS.in  -CAT_OUT masses.out -ZPHOTLIB BC03_COSMOS  -ADD_EMLINES 0,100 -AUTO_ADAPT NO -APPLY_SYSSHIFT $SHIFT  -ZFIX YES -PARA_OUT output_mass.para
+# Compute the photo-z without the adapation of the zero-points and applying the shifts given in keyword
+# Assume that shifts are derived previously. Shifts will be substracted to apparent magnitudes and to be used later using  --APPLY_SYSSHIFT $SHIFT (to avoid rerunning auto-adapt)
+export SHIFT="0.051121,-0.016466,-0.084605,-0.058664,-0.039816,-0.053934,-0.072388,-0.023517,-0.001748,0.013681,-0.036044,0.041802,0.017099,0.006699,0.006539,0.021207,0.063196,0.015433,-0.163286,0.032249,0.039840,0.050121,0.025781,0.012588,0.052423,0.034833,0.018244,-0.024611,0.000000,0.000000"
+zphota -c COSMOS.para --CAT_IN COSMOS.in  --CAT_OUT zphot_vista_shift.out --ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS   --ADD_EMLINES 0,100 --AUTO_ADAPT NO --APPLY_SYSSHIFT $SHIFT --PDZ_OUT zphot_vista_shift --PDZ_TYPE BAY_ZG
 
+# Test with no emission lines
+zphota -c COSMOS.para --CAT_IN COSMOS.in  --CAT_OUT zphot_vista_noline.out --ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS  --ADD_EMLINES 0,0 --AUTO_ADAPT NO --APPLY_SYSSHIFT $SHIFT
 
-#######################################################
-#  FIGURES
-#######################################################
+# Test with no offsets from auto-adapt
+zphota -c COSMOS.para --CAT_IN COSMOS.in  --CAT_OUT zphot_vista_noadapt.out --ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS   --ADD_EMLINES 0,100 --AUTO_ADAPT NO 
+
+# Extract Id*.spec and Id*.chi for few sources
+zphota -c COSMOS.para --CAT_IN COSMOS.in  --CAT_OUT zphot_vista_select.out --ZPHOTLIB VISTA_COSMOS_FREE,ALLSTAR_COSMOS,QSO_COSMOS   --ADD_EMLINES 0,100 --AUTO_ADAPT NO --APPLY_SYSSHIFT $SHIFT --PDZ_OUT NONE --PDZ_TYPE BAY_ZG --CAT_LINES 12000,12010 --SPEC_OUT YES 
 
-# Diadgnostics on the photo-z quality
+# FIGURES
+# Diagnostics on the photo-z quality
 # results in figuresLPZ.pdf
 python figuresLPZ.py zphot_vista_adapt.out
 
-# Diadgnostics on the physical parameter quality
-# results in figuresLPP.pdf
-python figuresLPP.py masses.out
-
 # Plot the best-fit templates
 # Results in MULTISPEC.pdf
 python spec.py Id*spec -d pdf
 
 
 
-#################
+###################################################
+# PHYSICAL PARAMETERS
+###################################################
+
+
+# Generate the library of predicted magnitudes for galaxy templates to be used for physical parameters
+# BC03 library 
+# Library to derive physical parameters 
+# Read the galaxy BC03 templates and store them in $WORK/lib_bin
+sedtolib -c COSMOS.para -t G --GAL_SED BC03COMB_MOD.list --GAL_LIB LIB_COMB --SEL_AGE AGE_BC03COMB.dat
+# Use the galaxy BC03 templates + filters to derive a library of predicted magnitudes and store in $WORK/lib_mag
+# Example with a larger dz=0.03 step to save speed
+# Use the "physical" recipes to generate the lines
+# Add the dust component with one template of Bethermin+12
+mag_gal  -c COSMOS.para -t G --GAL_LIB_IN LIB_COMB  --GAL_LIB_OUT BC03_COSMOS  --MOD_EXTINC 0,12,0,12,0,12 --EXTINC_LAW SB_calzetti.dat,extlaw_0.9.dat,salim18.dat --EB_V 0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7 --Z_STEP 0.03,0,6 --EM_LINES PHYS --EM_DISPERSION 1 --ADD_DUSTEM YES
+
+# Derive the stellar masses. Can be done with BC03 templates.
+zphota -c COSMOS.para --CAT_IN COSMOS.in  --CAT_OUT masses.out --ZPHOTLIB BC03_COSMOS  --ADD_EMLINES 0,100 --AUTO_ADAPT NO --APPLY_SYSSHIFT $SHIFT  --ZFIX YES --PARA_OUT output_mass.para
+
+# FIGURES
+# Diagnostics on the physical parameter quality
+# results in figuresLPP.pdf
+python figuresLPP.py masses.out
+
+
 
 ######################################
-#  LIR
+#  Far IR
 #######################################################
 
 # Filters including the FIR
-$LEPHAREDIR/bin/filter -c HERSCHEL.para
+filter -c HERSCHEL.para
 
 #IR library from Dale
-$LEPHAREDIR/bin/sedtolib -c HERSCHEL.para -t G -GAL_SED DALE.list  -GAL_LIB LIB_DALE 
-$LEPHAREDIR/bin/mag_gal  -c HERSCHEL.para -t G -GAL_LIB_IN LIB_DALE  -GAL_LIB_OUT DALE_COSMOS -MOD_EXTINC 0,0
+sedtolib -c HERSCHEL.para -t G --GAL_SED DALE.list  --GAL_LIB LIB_DALE 
+mag_gal  -c HERSCHEL.para -t G --GAL_LIB_IN LIB_DALE  --GAL_LIB_OUT DALE_COSMOS --MOD_EXTINC 0,0
 
 # BC03 SED
-$LEPHAREDIR/bin/sedtolib -c HERSCHEL.para -t G -GAL_SED BC03COMB_MOD.list -GAL_LIB LIB_BC03CAL -SEL_AGE AGE_BC03COMB.dat
-$LEPHAREDIR/bin/mag_gal  -c HERSCHEL.para -t G -GAL_LIB_IN LIB_BC03CAL -GAL_LIB_OUT BC_COSMOS  -MOD_EXTINC 1,14,1,14 -EXTINC_LAW SB_calzetti.dat,extlaw_0.9.dat -EB_V 0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7 
+sedtolib -c HERSCHEL.para -t G --GAL_SED BC03COMB_MOD.list --GAL_LIB LIB_BC03CAL --SEL_AGE AGE_BC03COMB.dat
+mag_gal  -c HERSCHEL.para -t G --GAL_LIB_IN LIB_BC03CAL --GAL_LIB_OUT BC_COSMOS  --MOD_EXTINC 1,14,1,14 --EXTINC_LAW SB_calzetti.dat,extlaw_0.9.dat --EB_V 0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7 
 
 # Fit in FIR
-$LEPHAREDIR/bin/zphota -c HERSCHEL.para -CAT_IN MIPS.in -CAT_OUT lumFIR.out -ZPHOTLIB BC_COSMOS  -ADD_EMLINES 0,100 -AUTO_ADAPT NO  -ZFIX YES -FIR_LIB DALE_COSMOS  -FIR_SUBSTELLAR YES  -FIR_FREESCALE YES  -FIR_CONT 8257536
+zphota -c HERSCHEL.para --CAT_IN MIPS.in --CAT_OUT lumFIR.out --ZPHOTLIB BC_COSMOS  --ADD_EMLINES 0,100 --AUTO_ADAPT NO  --ZFIX YES --FIR_LIB DALE_COSMOS  --FIR_SUBSTELLAR YES  --FIR_FREESCALE YES  --FIR_CONT 8257536
 
 # Fit in FIR with SED in output
-$LEPHAREDIR/bin/zphota -c HERSCHEL.para -CAT_IN MIPS.in -CAT_OUT testFIR.out -ZPHOTLIB BC_COSMOS  -ADD_EMLINES 0,100 -AUTO_ADAPT NO  -ZFIX YES -FIR_LIB DALE_COSMOS  -FIR_SUBSTELLAR YES  -FIR_FREESCALE YES  -FIR_CONT 8257536 -CAT_LINES 11980,12010 -SPEC_OUT YES
+zphota -c HERSCHEL.para --CAT_IN MIPS.in --CAT_OUT testFIR.out --ZPHOTLIB BC_COSMOS  --ADD_EMLINES 0,100 --AUTO_ADAPT NO  --ZFIX YES --FIR_LIB DALE_COSMOS  --FIR_SUBSTELLAR YES  --FIR_FREESCALE YES  --FIR_CONT 8257536 --CAT_LINES 11980,12010 --SPEC_OUT YES
 
 # Plot the best-fit templates
 # Results in MULTISPEC.pdf
@@ -125,7 +133,7 @@ python spec.py Id*spec -d pdf
 #########################################################
 
 
-$LEPHAREDIR/bin/filter -c COSMOS.para  -FILTER_LIST cosmos/u_cfht.lowres,cosmos/B_subaru.lowres,cosmos/V_subaru.lowres,cosmos/r_subaru.lowres,cosmos/i_subaru.lowres,cosmos/suprime_FDCCD_z.res,vista/Y.lowres,vista/J.lowres,vista/H.lowres,vista/K.lowres -FILTER_CALIB 0 -FILTER_FILE filter_extinc
+filter -c COSMOS.para  --FILTER_LIST cosmos/u_cfht.lowres,cosmos/B_subaru.lowres,cosmos/V_subaru.lowres,cosmos/r_subaru.lowres,cosmos/i_subaru.lowres,cosmos/suprime_FDCCD_z.res,vista/Y.lowres,vista/J.lowres,vista/H.lowres,vista/K.lowres --FILTER_CALIB 0 --FILTER_FILE filter_extinc
 # Read the filters and compute extinctions
-$LEPHAREDIR/bin/filter_extinc -c COSMOS.para -FILTER_FILE filter_extinc.dat
+filter_extinc -c COSMOS.para --FILTER_FILE filter_extinc.dat
 

ext/README.extinc

@@ -25,12 +25,15 @@
    SB_calzetti_bump3.dat
 
 #
-# Added dust law
+# Added dust laws
 
 extlaw_0.9.dat
 was created in Arnouts+2013 
 
 
+salim18.dat
+added the low mass high-z analog from Salim+2018
+
 #
 # The atmospheric extinction curves for other purposes 
   There are also atmospheric extinctions from La Silla (ESO) and CTIO