Commit c29dd74f authored by ILBERT Olivier's avatar ILBERT Olivier
Browse files

Improve the documentation, and clean few SED directories

parent 5fd7e95f
##############################################################################################
########### CREATION OF LIBRARIES FROM SEDs List #############
########### Run : $ZPHOTDIR/source/sedtolib -t (S/Q/G) -c zphot.para #############
##############################################################################################
#
#------------------- STELLAR LIBRARY (ASCII SEDs) ---------------------------
STAR_SED STAR_SWIRE.list # STAR list (full path)
STAR_LIB LIB_STAR # Binary STAR LIBRARY (-> $ZPHOTWORK/lib_bin/*)
STAR_FSCALE 3.432E-09 # Arbitrary Flux Scale
#
#------------------- QSO LIBRARY (ASCII SEDs) ---------------------------
QSO_SED AGN_LONSDALE.list # QSO list (full path)
QSO_LIB LIB_QSO # Binary QSO LIBRARY (-> $ZPHOTWORK/lib_bin/*)
QSO_FSCALE 1. # Arbitrary Flux Scale
#
#------------------- GALAXY LIBRARY (ASCII or BINARY SEDs) ---------------------------
GAL_SED CE_MOD.list # GALAXMuzzin09_SEDY list (full path)
GAL_LIB LIB_CE # Binary GAL LIBRARY (-> $ZPHOTWORK/lib_bin/*)
GAL_FSCALE 1. # Arbitrary Flux Scale
#SEL_AGE /data/zphot_vers25_03_03/sed/GAL/HYPERZ/AGE_GISSEL_HZ.dat # List of Age for GISSEL(full path)
AGE_RANGE 0.,15.e9 # Age Min-Max in yr
#
##############################################################################################
########### FILTERS #############
########### Run : $ZPHOTDIR/source/filter -c zphot.para #############
##############################################################################################
#
FILTER_REP $LEPHAREDIR/filt # Repository in which the filters are stored
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,hsc/yHSC.pb,vista/Y.lowres,vista/J.lowres,vista/H.lowres,vista/K.lowres,cosmos/H_wircam.lowres,cosmos/K_wircam.lowres,cosmos/IB427.lowres,cosmos/IB464.lowres,cosmos/IB484.lowres,cosmos/IB505.lowres,cosmos/IB527.lowres,cosmos/IB574.lowres,cosmos/IB624.lowres,cosmos/IB679.lowres,cosmos/IB709.lowres,cosmos/IB738.lowres,cosmos/IB767.lowres,cosmos/IB827.lowres,cosmos/NB711.lowres,cosmos/NB816.lowres,cosmos/irac_ch1.lowres,cosmos/irac_ch2.lowres,cosmos/irac_ch3.lowres,cosmos/irac_ch4.lowres,galex/NUV.lowres,galex/F1600.pb,cfht/u_megaprime_sagem.res,cfht/g_megaprime_sagem.res,cfht/r_megaprime_sagem.res,cfht/i_megaprime_sagem.res,cfht/iMegaNew.fil,cfht/z_megaprime_sagem.res
TRANS_TYPE 1 # TRANSMISSION TYPE
# 0[-def]: Energy, 1: Nb of photons
FILTER_CALIB 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0 # 0[-def]: fnu=ctt
# 1 : nu.fnu=ctt
# 2 : fnu=nu
# 3 : fnu=Black Body @ T=10000K
# 4 : for MIPS (leff with nu fnu=ctt and flux with BB @ 10000K
FILTER_FILE filter_cosmos # name of file with filters (-> $ZPHOTWORK/filt/)
#
##############################################################################################
########### THEORETICAL MAGNITUDES #############
########### Run : $ZPHOTDIR/source/mag_star -c zphot.para (star only) #######
########### Run : $ZPHOTDIR/source/mag_gal -t (Q or G) -c zphot.para (Gal. & QSO) #######
##############################################################################################
#
#------------------- From STELLAR LIBRARY ------------------------------
STAR_LIB_IN LIB_STAR # Input STELLAR LIBRARY (in $ZPHOTWORK/lib_bin/*)
STAR_LIB_OUT STAR_COSMOS # Output STELLAR MAGN. (-> $ZPHOTWORK/lib_mag/*)
#
#------------------- From QSO LIBRARY ------------------------------
QSO_LIB_IN LIB_QSO # Input QSO LIBRARY (in $ZPHOTWORK/lib_bin/*)
QSO_LIB_OUT QSO_COSMOS # Output QSO MAGN. (-> $ZPHOTWORK/lib_
#
#------------------- From GALAXY LIBRARY ------------------------------
GAL_LIB_IN LIB_CE # Input GALAXY LIBRARY (in $ZPHOTWORK/lib_bin/*)
GAL_LIB_OUT CE_COSMOS # Output GALAXY LIBRARY (-> $ZPHOTWORK/lib_mag/*)
#
#------------------ MAG + Z_STEP + COSMO + EXTINCTION -----------------------------
MAGTYPE AB # Magnitude type (AB or VEGA)
ZGRID_TYPE 0 # Define the kind of redshift grid (0: linear ; 1: dz*(1+z))
Z_STEP 0.01,0.,6. # dz, zmin, zmax
COSMOLOGY 70,0.3,0.7 # H0,om0,lbd0 (if lb0>0->om0+lbd0=1)
MOD_EXTINC 0,0 # model range for extinction
EXTINC_LAW SB_calzetti.dat # ext. law (in $ZPHOTDIR/ext/*)
EB_V 0.,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.5 # E(B-V) (<50 values)
EM_LINES EMP_UV # [NO/EMP_UV/EMP_SFR/PHYS] choice between emission line prescription
EM_DISPERSION 0.5,0.75,1.,1.5,2. # Dispersion allowed in the emission line flux factor
#
#------------------ ASCII OUTPUT FILES OPTION -------------------------------
LIB_ASCII NO # Writing Output in ASCII file (-> working directory)
#
##############################################################################################
########### PHOTOMETRIC REDSHIFTS #############
########### Run : $ZPHOTDIR/source/zphot -c zphot.para #############
##############################################################################################
#
#------------------ Input Catalog Informations -------------------------------
CAT_IN bidon # Input catalog (full path)
INP_TYPE F # Input type (F:Flux or M:MAG)
CAT_MAG AB # Input Magnitude (AB or VEGA)
CAT_FMT MEME # MEME: (Mag,Err)i
# MMEE: (Mag)i,(Err)i
CAT_LINES 0,1000000000 # MIN and MAX RANGE of ROWS used in input cat [def:-99,-99]
CAT_TYPE LONG # Input Format (LONG,SHORT-def)
BD_SCALE 0 # Bands used for scaling (Sum 2^n; n=0->nbd-1, 0->all bands)
GLB_CONTEXT 0 # Overwrite Context (Sum 2^n; n=0->nbd-1, 0->all bands, -1[-def] used context per object)
FORB_CONTEXT -1 # Not consider these filters in the fit (Sum 2^n; n=0->nbd-1)
ERR_SCALE 0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.05,0.05,0.05,0.05,0.05,0.05,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.05,0.1,0.2,0.2,0.1,0.3,0.02,0.02,0.02,0.02,0.02,0.02 # Systematic errors per band
ERR_FACTOR 1.5
#
#------------------ Theoretical libraries --------------------------------
ZPHOTLIB CE_COSMOS,STAR_COSMOS,QSO_COSMOS # Library used for Chi2 (max:3)
ADD_EMLINES 0,100 # range of model in which to apply emission lines
#
#------------------ Output catalog Name ---------------------------------
CAT_OUT zphot.out # Output catalog (-> working directory)
PARA_OUT output.para # Ouput parameters (full path)
#
################# PHOTOMETRIC REDSHIFTS OPTIONS
# FIR LIBRARY
FIR_LIB NONE
FIR_LMIN 7.0 # Lambda Min (micron) for FIR analysis
FIR_CONT -1
FIR_SCALE -1
FIR_FREESCALE YES # ALLOW FOR FREE SCALING
FIR_SUBSTELLAR NO
# PHYSICAL LIBRARY with Stochastic models from BC07
PHYS_LIB NONE
PHYS_CONT -1
PHYS_SCALE -1
PHYS_NMAX 100000
#
#------------------ Priors ---------------------------------
#MASS_SCALE 6.,16. # Lg(Scaling) min,max [0,0-def]
MAG_ABS -24,-5 # Mabs_min , Mabs_max [0,0-def]
MAG_ABS_AGN -30,-10 # Mabs_min , Mabs_max for AGN [0,0-def]
MAG_REF 2 # Reference number for band used by Mag_abs
#ZFORM_MIN 1,5,5,5,5,5,3,1 # Min. Zformation per SED -> Age constraint
Z_RANGE 0.,99.99 # Z min-max used for the Galaxy library
EBV_RANGE 0,9 # E(B-V) MIN-MAX RANGE of E(B-V) used
#NZ_PRIOR 4,5 # I Band for prior on N(z), the second number is the band to be used if the first is missing.
#
#------------------ Fixed Z (need format LONG for input Cat) ----------------------
ZFIX NO # fixed z and search best model [YES,NO-def]
#
#------------------ Parabolic interpolation for Zbest ----------------------------
Z_INTERP YES # redshift interpolation
#
Z_METHOD BEST # BEST/ML absolute magnitude, .spec, scaling will be given at this redshift
#------------------ Normalized ML(exp-(0.5*Chi^2)) curve analysis -------------------
#------- Secondary peak analysis
DZ_WIN 0.3 # Window search for 2nd peaks [0->5; 0.25-def]
MIN_THRES 0.02 # Lower theshold for 2nd peaks [0->1; 0.1-def]
#
#------- Probability (in %) per redshift intervals --------------------
#PROB_INTZ 0.01,0.1,0.1,0.2,0.2,0.3,0.3,0.4,0.4,0.5,0.5,0.6,0.6,0.7,0.7,0.8,0.8,0.9,0.9,1.,1.,1.1,1.1,1.3,1.3,1.5,1.5,2.,2.,3.,3.,4.,4.,5 # (even number)
######### ABSOLUTE MAGNITUDES COMPUTATION ###########
#
MABS_METHOD 1 # 0[-def] : obs->Ref
# 1 : best obs->Ref
# 2 : fixed obs->Ref
# 3 : mag from best SED
# 4 : Zbin
MABS_CONTEXT 2281702809
MABS_REF 11 # 0[-def]: filter obs chosen for Mabs :
# ONLY USED IF MABS_METHOD=2
MABS_FILT 1,2,3,4 # Chosen filters per redshift bin (MABS_ZBIN)
# ONLY USED IF MABS_METHOD=4
MABS_ZBIN 0,0.5,1,1.5,2,3,3.5,4 # Redshift bins (even number)
# ONLY USED IF MABS_METHOD=4
RF_COLORS 32,4,4,13
ADDITIONAL_MAG none # name of file with filters (-> $ZPHOTWORK/filt/). Predicted mag and abs mag in additional filters.
######### Z MAX #####################################
LIMITS_ZBIN 0,99 # Redshifts used to split in N bins, separated by a coma. Need N+1 values (start with the minimum redshift).
LIMITS_MAPP_REF 1 # Band in which the absolute magnitude is computed
LIMITS_MAPP_SEL 1 # Give the selection band in each redshift bin. Need 1 or N values.
LIMITS_MAPP_CUT 90 # Magnitude cut used in each redshift bin. Need 1 or N values.
######### OUTPUT SPECTRA ###########
#
SPEC_OUT NO # spectrum for each object? [YES,NO-def]
CHI2_OUT NO # output file with all values : z,mod,chi2,E(B-V),...
# BE CAREFUL can take a lot of space !!
######### OUTPUT PDZ ANALYSIS
PDZ_OUT NONE # pdz output file name [def-NONE]
# add automatically PDZ_OUT[.pdz/.mabsx/.mod/.zph]
PDZ_MABS_FILT 2,10,14 # MABS for REF FILTERS to be extracted
#
######### FAST MODE : color-space reduction #####
#
FAST_MODE NO # Fast computation [NO-def]
COL_NUM 3 # Number of colors used [3-def]
COL_SIGMA 3 # Enlarge of the obs. color-errors[3-def]
COL_SEL AND # Combination between used colors [AND/OR-def]
#
######### MAGNITUDE SHIFTS applied to libraries ######
#
# APPLY_SYSSHIFT 0. # Apply systematic shifts in each band
# used only if number of shifts matches
# with number of filters in the library
#
######### ADAPTIVE METHOD using Z spectro sample ###
#
AUTO_ADAPT NO # Adapting method with spectro [NO-def]
ADAPT_BAND 5 # Reference band, band1, band2 for color
ADAPT_LIM 1.5,20.5 # Mag limits for spectro in Ref band [18,21.5-def]
ADAPT_CONTEXT -1 # Context for bands used for training
# -1[-def] used context per object
ADAPT_ZBIN 0.01,6 # Redshift's interval used for training
# [0.001,6-Def]
ADAPT_MODBIN 1,1000 # Model's interval used for training
# [1,1000-Def]
ERROR_ADAPT NO # [YES,NO-def]
# Add error in quadrature according to
# the difference between observed
# and predicted apparent magnitudes
#
# BASIC PHOTOZ INFORMATION
IDENT
Z_BEST
#Z_BEST68_LOW
#Z_BEST68_HIGH
#Z_BEST90_LOW
#Z_BEST90_HIGH
#Z_BEST99_LOW
#Z_BEST99_HIGH
Z_ML
Z_ML68_LOW
Z_ML68_HIGH
#Z_ML90_LOW
#Z_ML90_HIGH
#Z_ML99_LOW
#Z_ML99_HIGH
# INFORMATIONS ON THE BEST FIT TEMPLATE
CHI_BEST
MOD_BEST
EXTLAW_BEST
EBV_BEST
PDZ_BEST
SCALE_BEST
NBAND_USED
NBAND_ULIM
####### Z- 2nd solution - GALAXY ###########
Z_SEC
CHI_SEC
MOD_SEC
#AGE_SEC
#EBV_SEC
#EXTLAW_SEC
#PDZ_SEC
#SCALE_SEC
######## Best fit STAR #############
CHI_STAR
MOD_STAR
######## QSO #############
CHI_QSO
Z_QSO
MOD_QSO
#Z_ML_QSO
#Z_ML68_LOW_QSO
#Z_ML68_HIGH_QSO
######## MODEL MAGNITUDES FOR GALAXY ####
######## K-COR , ABSOLUTE MAG, FILTER CHOICE FOR GALAXY ##########
#MAG_MOD()
#MAG_PRED()
#K_COR()
#MAG_ABS()
#MABS_FILT()
#EMAG_ABS()
#ABSMAG_PRED()
# PHYSICAL PARAMETER FROM THE BEST FIT TEMPLATE (AT MINIMUM CHI2)
#EBV_BEST
#EBV_MED
#EBV_INF
#EBV_SUP
#AGE_BEST
#AGE_MED
#AGE_INF
#AGE_SUP
#LDUST_BEST
#LDUST_MED
#LDUST_INF
#LDUST_SUP
#MASS_BEST
#MASS_MED
#MASS_INF
#MASS_SUP
#SFR_BEST
#SFR_MED
#SFR_INF
#SFR_SUP
#SSFR_BEST
#SSFR_MED
#SSFR_INF
#SSFR_SUP
#COL1_MED
#COL1_INF
#COL1_SUP
#COL2_MED
#COL2_INF
#COL2_SUP
#MREF_MED
#MREF_INF
#MREF_SUP
#LUM_NUV_BEST
#LUM_R_BEST
#LUM_K_BEST
## LUM TIR FROM FIR FIT
#LUM_TIR_BEST
#LUM_TIR_MED
#LUM_TIR_INF
#LUM_TIR_SUP
######## LIMITS FROM MAG CUT ####
#LIMITS_ZMAX
#LIMITS_MFAINT
####### INPUT INFO ##############
#CONTEXT
ZSPEC
MAG_OBS()
ERR_MAG_OBS()
STRING_INPUT
// GALAXY
stout.unsetf(ios::fixed); //There is a problem with Z_BEST: for rare pathological galaxies is =-1.e9
if((outkeywords[k]).compare("Z_BEST") == 0)stout << setw(16) << setprecision(4) << zmin[0] << " ";
stout.setf(ios::fixed, ios::floatfield); stout.precision(4); //set float format and the n. of decimal digits
if((outkeywords[k]).compare("Z_BEST68_LOW") == 0)stout << setw(12) << z68[0] << " ";
if((outkeywords[k]).compare("Z_BEST68_HIGH") == 0)stout << setw(12) << z68[1] << " ";
if((outkeywords[k]).compare("Z_BEST90_LOW") == 0)stout << setw(12) << z90[0] << " ";
if((outkeywords[k]).compare("Z_BEST90_HIGH") == 0)stout << setw(12) << z90[1] << " ";
if((outkeywords[k]).compare("Z_BEST99_LOW") == 0)stout << setw(12) << z99[0] << " ";
if((outkeywords[k]).compare("Z_BEST99_HIGH") == 0)stout << setw(12) << z99[1] << " ";
if((outkeywords[k]).compare("Z_ML") == 0)stout << setw(12) << zml << " ";
if((outkeywords[k]).compare("Z_ML68_LOW") == 0)stout << setw(12) << zml68[0] << " ";
if((outkeywords[k]).compare("Z_ML68_HIGH") == 0)stout << setw(12) << zml68[1] << " ";
if((outkeywords[k]).compare("Z_ML90_LOW") == 0)stout << setw(12) << zml90[0] << " ";
if((outkeywords[k]).compare("Z_ML90_HIGH") == 0)stout << setw(12) << zml90[1] << " ";
if((outkeywords[k]).compare("Z_ML99_LOW") == 0)stout << setw(12) << zml99[0] << " ";
if((outkeywords[k]).compare("Z_ML99_HIGH") == 0)stout << setw(12) << zml99[1] << " ";
// Best fit values
double massBest = -99.9;
double sfrBest = -99.9;
double ssfrBest = -99.9;
double ldustBest= -99.9;
double lirBest = -99.9;
double ageBest = -99.9;
double ebvBest = -99.9;
double extlBest = -99.9;
double luvBest = -99.9;
double loptBest = -99.9;
double lnirBest = -99.9;
if(indmin[0]>=0){
massBest =log10((fulllib[indmin[0]])->mass * dmmin[0]);
sfrBest =log10((fulllib[indmin[0]])->sfr * dmmin[0]);
ssfrBest =log10((fulllib[indmin[0]])->sfr / (fulllib[indmin[0]])->mass);
ldustBest=(fulllib[indmin[0]])->ltir + log10(dmmin[0]);
ageBest =(fulllib[indmin[0]])->age ;
ebvBest =(fulllib[indmin[0]])->ebv;
extlBest =(fulllib[indmin[0]])->extlawId;
luvBest =(fulllib[indmin[0]])->luv + log10(3.e18*400/pow(2300,2)) -log10(Lsol)+log10(dmmin[0]);
loptBest =(fulllib[indmin[0]])->lopt+ log10(3.e18*1000/pow(6000,2)) -log10(Lsol)+log10(dmmin[0]);
lnirBest =(fulllib[indmin[0]])->lnir+ log10(3.e18*2000/pow(22000,2))-log10(Lsol)+log10(dmmin[0]);
}
if(indminIR>=0)lirBest =(fulllibIR[indminIR])->ltir + log10(dmminIR);
//Integer
stout.precision(0); // in fixed format, number of decimal digits in the output
if((outkeywords[k]).compare("MOD_BEST") == 0) stout << setw(6) << std::fixed << imasmin[0] << " " ;
if((outkeywords[k]).compare("EXTLAW_BEST") == 0) stout << setw(6) << std::fixed << extlBest << " ";
//Float
stout.setf(ios::fixed, ios::floatfield);
if((outkeywords[k]).compare("EBV_BEST") == 0) stout << setw(12) << setprecision(3) << ebvBest << " ";
if((outkeywords[k]).compare("MASS_BEST") == 0) stout << setw(12) << setprecision(5) << massBest << " ";
if((outkeywords[k]).compare("LDUST_BEST") == 0) stout << setw(12) << setprecision(3) << ldustBest << " ";
if((outkeywords[k]).compare("LUM_TIR_BEST")== 0) stout << setw(12) << setprecision(3) << lirBest << " ";
if((outkeywords[k]).compare("LUM_NUV_BEST") == 0)stout << setw(12) << setprecision(3) << luvBest << " ";
if((outkeywords[k]).compare("LUM_R_BEST") == 0) stout << setw(12) << setprecision(3) << loptBest << " ";
if((outkeywords[k]).compare("LUM_K_BEST") == 0) stout << setw(12) << setprecision(3) << lnirBest << " ";
//Scientific
stout.setf(ios::scientific, ios::floatfield);
//note: PDZ_BEST has to be fixed
if((outkeywords[k]).compare("PDZ_BEST") == 0) stout << setw(14) << setprecision(5) << 99. << " ";
if((outkeywords[k]).compare("CHI_BEST") == 0) stout << setw(14) << setprecision(5) << chimin[0]<< " ";
if((outkeywords[k]).compare("SCALE_BEST") == 0) stout << setw(14) << setprecision(5) << dmmin[0] << " ";
if((outkeywords[k]).compare("AGE_BEST") == 0) stout << setw(12) << setprecision(3) << ageBest << " " ;
if((outkeywords[k]).compare("SFR_BEST") == 0) stout << setw(12) << setprecision(3) << sfrBest << " ";
if((outkeywords[k]).compare("SSFR_BEST") == 0) stout << setw(12) << setprecision(3) << ssfrBest << " ";
/*
* The following parameters are not implemented yet:
* PDZ()
*/
// GENERAL
if((outkeywords[k]).compare("STRING_INPUT") == 0)stout << " " << str_inp << " " ; //string-like
stout.setf(ios::fixed); stout.precision(0); //Integer:
if((outkeywords[k]).compare("CONTEXT") == 0)stout << setw(16) << cont << " ";
if((outkeywords[k]).compare("NBAND_USED") == 0)stout << setw(5) << nbused << " ";
if((outkeywords[k]).compare("NBAND_ULIM") == 0)stout << setw(3) << nbul << " ";
stout.setf(ios::fixed, ios::floatfield); //Float:
if((outkeywords[k]).compare("ZSPEC") == 0)stout << setw(14) << setprecision(5) << zs << " " ;
if((outkeywords[k]).compare("MAG_OBS()") == 0){for (size_t l=0; l<mab.size(); l++)stout << setw(10) << setprecision(3) << mab[l] << " ";}
if((outkeywords[k]).compare("ERR_MAG_OBS()") == 0){for (size_t l=0; l<msab.size(); l++)stout << setw(10) << setprecision(3) << msab[l] << " ";}
if((outkeywords[k]).compare("MAG_MOD()") == 0){for (size_t l=0; l<magm.size(); l++)stout << setw(10) << setprecision(3) << magm[l] << " ";}
// QSO and STARS
stout.setf(ios::fixed); stout.precision(0); //Integer:
if((outkeywords[k]).compare("MOD_QSO") == 0)stout << setw(12) << imasmin[1] << " ";
if((outkeywords[k]).compare("MOD_STAR") == 0)stout << setw(12) << imasmin[2] << " ";
stout.setf(ios::fixed, ios::floatfield); stout.precision(4); //Float:
if((outkeywords[k]).compare("Z_QSO") == 0)stout << setw(12) << zmin[1] << " " ;
if((outkeywords[k]).compare("Z_ML_QSO") == 0)stout << setw(12) << zqml << " ";
if((outkeywords[k]).compare("Z_ML68_LOW_QSO") == 0)stout << setw(12) << zqml68[0] << " ";
if((outkeywords[k]).compare("Z_ML68_HIGH_QSO") == 0)stout << setw(12) << zqml68[1] << " ";
stout.setf(ios::scientific, ios::floatfield); stout.precision(5); //Scientific:
if((outkeywords[k]).compare("CHI_QSO") == 0)stout << setw(14) << chimin[1] << " ";
if((outkeywords[k]).compare("CHI_STAR") == 0)stout << setw(14) << chimin[2] << " ";
// SECONDARY GALAXY PEAK
stout.setf(ios::fixed); stout.precision(0); //Integer:
if((outkeywords[k]).compare("MOD_SEC") == 0)stout << setw(12) << zsecMod << " ";
if((outkeywords[k]).compare("EXTLAW_SEC")== 0)stout << setw(12) << zsecExtlaw<< " ";
stout.setf(ios::fixed, ios::floatfield); stout.precision(4); //Float:
if((outkeywords[k]).compare("Z_SEC") == 0)stout << setw(12) << zsec << " ";
if((outkeywords[k]).compare("EBV_SEC") == 0)stout << setw(12) << zsecEbv<< " ";
stout.setf(ios::scientific, ios::floatfield); stout.precision(5); //Scientific:
if((outkeywords[k]).compare("PDZ_SEC") == 0)stout << setw(14) << zsecProb << " " ;
if((outkeywords[k]).compare("CHI_SEC") == 0)stout << setw(14) << zsecChi2 << " " ;
if((outkeywords[k]).compare("SCALE_SEC") == 0)stout << setw(14) << zsecScale << " " ;
if((outkeywords[k]).compare("AGE_SEC") == 0)stout << setw(14) << zsecAge << " " ;
// PHYSICAL PARAMETERS
stout.setf(ios::fixed, ios::floatfield); stout.precision(5); //Float:
if((outkeywords[k]).compare("AGE_MED") == 0)stout << setw(12) << ageml << " " ;
if((outkeywords[k]).compare("AGE_INF") == 0)stout << setw(12) << ageml68[0] << " ";
if((outkeywords[k]).compare("AGE_SUP") == 0)stout << setw(12) << ageml68[1] << " ";
if((outkeywords[k]).compare("LDUST_MED") == 0)stout << setw(12) << Ldustml << " " ;
if((outkeywords[k]).compare("LDUST_INF") == 0)stout << setw(12) << Ldustml68[0] << " ";
if((outkeywords[k]).compare("LDUST_SUP") == 0)stout << setw(12) << Ldustml68[1] << " ";
if((outkeywords[k]).compare("LUM_TIR_MED") == 0)stout << setw(12) << LIRml << " " ;
if((outkeywords[k]).compare("LUM_TIR_INF") == 0)stout << setw(12) << LIRml68[0] << " ";
if((outkeywords[k]).compare("LUM_TIR_SUP") == 0)stout << setw(12) << LIRml68[1] << " ";
if((outkeywords[k]).compare("MASS_MED") == 0)stout << setw(12) << massml << " ";
if((outkeywords[k]).compare("MASS_INF") == 0)stout << setw(12) << massml68[0] << " ";
if((outkeywords[k]).compare("MASS_SUP") == 0)stout << setw(12) << massml68[1] << " ";
if((outkeywords[k]).compare("SFR_MED") == 0)stout << setw(12) << SFRml << " ";
if((outkeywords[k]).compare("SFR_INF") == 0)stout << setw(12) << SFRml68[0] << " ";
if((outkeywords[k]).compare("SFR_SUP") == 0)stout << setw(12) << SFRml68[1] << " ";
if((outkeywords[k]).compare("SSFR_MED") == 0)stout << setw(12) << sSFRml << " ";
if((outkeywords[k]).compare("SSFR_INF") == 0)stout << setw(12) << sSFRml68[0] << " ";
if((outkeywords[k]).compare("SSFR_SUP") == 0)stout << setw(12) << sSFRml68[1] << " ";
if((outkeywords[k]).compare("COL1_MED") == 0)stout << setw(12) << col1ml << " ";
if((outkeywords[k]).compare("COL1_INF") == 0)stout << setw(12) << col1ml68[0] << " ";
if((outkeywords[k]).compare("COL1_SUP") == 0)stout << setw(12) << col1ml68[1] << " ";
if((outkeywords[k]).compare("COL2_MED") == 0)stout << setw(12) << col2ml << " ";
if((outkeywords[k]).compare("COL2_INF") == 0)stout << setw(12) << col2ml68[0] << " ";
if((outkeywords[k]).compare("COL2_SUP") == 0)stout << setw(12) << col2ml68[1] << " ";
if((outkeywords[k]).compare("EBV_MED") == 0)stout << setw(12) << ebvml << " ";
if((outkeywords[k]).compare("EBV_INF") == 0)stout << setw(12) << ebvml68[0] << " ";
if((outkeywords[k]).compare("EBV_SUP") == 0)stout << setw(12) << ebvml68[1] << " ";
if((outkeywords[k]).compare("MREF_MED") == 0)stout << setw(12) << Mrefml << " ";
if((outkeywords[k]).compare("MREF_INF") == 0)stout << setw(12) << Mrefml68[0] << " ";
if((outkeywords[k]).compare("MREF_SUP") == 0)stout << setw(12) << Mrefml68[1] << " ";
// ABSOLUTE MAGNITUDES
stout.setf(ios::fixed); stout.precision(0); //Integer:
if((outkeywords[k]).compare("MABS_FILT()") == 0){for (size_t l=0; l<absfilt.size(); l++)stout << setw(4) << absfilt[l] << " " ;}
stout.setf(ios::fixed, ios::floatfield); stout.precision(3); //Float:
if((outkeywords[k]).compare("K_COR()") == 0){for (size_t l=0; l<kap.size(); l++)stout << setw(13) << setprecision(5) << kap[l] << " " ;}
if((outkeywords[k]).compare("MAG_ABS()") == 0){for (size_t l=0; l<mabs.size(); l++)stout << setw(12) << setprecision(5) << mabs[l] << " ";}
if((outkeywords[k]).compare("EMAG_ABS()") == 0){for (size_t l=0; l<emabs.size(); l++)stout << setw(13) << setprecision(5) << emabs[l] << " ";}
// PREDICTED MAGNITUDES IN ADDITIONAL FILTERS
stout.setf(ios::scientific, ios::floatfield); stout.precision(5);
if((outkeywords[k]).compare("MAG_PRED()") == 0){for (size_t l=0; l<magPred.size(); l++)stout << setw(15) << magPred[l] << " " ;}
if((outkeywords[k]).compare("ABSMAG_PRED()") == 0){for (size_t l=0; l<absmagPred.size(); l++)stout << setw(15) << absmagPred[l] << " " ;}
// EMISSION LINES
stout.setf(ios::scientific, ios::floatfield); stout.precision(5);
if((outkeywords[k]).compare("EM_FLUX_LYA") == 0) stout << setw(15) << setprecision(5) << flux_em[0] << " " ;
if((outkeywords[k]).compare("EM_FLUX_OII") == 0) stout << setw(15) << setprecision(5) << flux_em[1] << " ";
if((outkeywords[k]).compare("EM_FLUX_HB") == 0) stout << setw(15) << setprecision(5) << flux_em[2] << " ";
if((outkeywords[k]).compare("EM_FLUX_OIIIA") == 0) stout << setw(15) << setprecision(5) << flux_em[3] << " " ;
if((outkeywords[k]).compare("EM_FLUX_OIIIB") == 0) stout << setw(15) << setprecision(5) << flux_em[4] << " " ;
if((outkeywords[k]).compare("EM_FLUX_HA") == 0) stout << setw(15) << setprecision(5) << flux_em[5] << " " ;
if((outkeywords[k]).compare("EM_FLUX_SIIIA") == 0) stout << setw(15) << setprecision(5) << flux_em[6] << " " ;
if((outkeywords[k]).compare("EM_FLUX_SIIIB") == 0) stout << setw(15) << setprecision(5) << flux_em[7] << " " ;
if((outkeywords[k]).compare("EM_EW_LYA") == 0) stout << setw(15) << setprecision(5) << EW_em[0] << " " ;
if((outkeywords[k]).compare("EM_EW_OII") == 0) stout << setw(15) << setprecision(5) << EW_em[1] << " " ;
if((outkeywords[k]).compare("EM_EW_HB") == 0) stout << setw(15) << setprecision(5) << EW_em[2] << " " ;
if((outkeywords[k]).compare("EM_EW_OIIIA") == 0) stout << setw(15) << setprecision(5) << EW_em[3] << " " ;
if((outkeywords[k]).compare("EM_EW_OIIIB") == 0) stout << setw(15) << setprecision(5) << EW_em[4] << " " ;
if((outkeywords[k]).compare("EM_EW_HA") == 0) stout << setw(15) << setprecision(5) << EW_em[5] << " " ;
if((outkeywords[k]).compare("EM_EW_SIIIA") == 0) stout << setw(15) << setprecision(5) << EW_em[6] << " " ;
if((outkeywords[k]).compare("EM_EW_SIIIB") == 0) stout << setw(15) << setprecision(5) << EW_em[7] << " " ;
// LIMITS FOR VMAX
stout.setf(ios::fixed, ios::floatfield); stout.precision(4); //Float:
if((outkeywords[k]).compare("LIMITS_ZMAX") == 0)stout << setw(12) << limits_zmax << " " ;
if((outkeywords[k]).compare("LIMITS_MFAINT") == 0)stout << setw(12) << limits_Mfaint << " " ;
\ No newline at end of file
macro read plot_sed.sm
plot 0 2.9 4 30 35
macro read plot_sed.sm
plot 0 2.9 4 30 35
macro read plot_sed.sm
plot 0 2.9 4 30 35
macro read plot_sed.sm
plot 0 2.9 4 0 35
macro read plot_sed.sm
plot 0 2.9 4 0 35
plot 0 2.9 4 0 10
plot 0 2.9 4 -5 10
macro read plot_sed.sm
plot 0 2.9 4 -1 3
plot 0 2.9 5 -1 3
plot 0 2.9 4.3 -1 3
plot 0.8 2.9 4.3 -1 3
macro read plot_sed.sm
plot 0.8 2.9 4.3 -1 3
plot 0. 2.9 4.3 -1 3
plot 0. 2.9 40.3 -1 3
plot 0. 2.9 400.3 -1 3
macro read plot_sed.sm
plot 0. 2.9 400.3 -1 3
plot 0. 2.9 4.3 -1 10
macro read plot_sed.sm
plot 0.8 2.9 4.3 -1 10
plot 0.8 3.2 4.3 5 10
macro read plot_sed.sm
plot 0.8 3.2 4.3 5 10
plot 0.9 3.2 4.3 5 10
plot 0.9 3.6 4.1 6 9
macro read plot_sed.sm
plot 0.9 3.6 4.1 6 9
plot 0.8 3.6 4.1 6 9
plot 0.6 3.6 4.1 6 9
macro read plot_sed.sm
plot 0.6 3.6 4.1 6 9
!pwd
plot 0.6 3.6 4.1 6 9
CFHTLS_SED/Ell_cww.sed
CFHTLS_SED/Ell_cww_0.sed
CFHTLS_SED/Ell_cww_1.sed
CFHTLS_SED/Ell_cww_2.sed
CFHTLS_SED/Ell_cww_3.sed
CFHTLS_SED/Ell_cww_4.sed
CFHTLS_SED/Ell_cww_5.sed
CFHTLS_SED/Ell_cww_6.sed
CFHTLS_SED/Ell_cww_7.sed
CFHTLS_SED/Ell_cww_8.sed
CFHTLS_SED/Ell_cww_9.sed
CFHTLS_SED/Ell_cww_10.sed
CFHTLS_SED/Ell_cww_11.sed
CFHTLS_SED/Ell_cww_12.sed
CFHTLS_SED/Ell_cww_13.sed
CFHTLS_SED/Ell_cww_14.sed
CFHTLS_SED/Ell_cww_15.sed
CFHTLS_SED/Ell_cww_16.sed
CFHTLS_SED/Ell_cww_17.sed
CFHTLS_SED/Ell_cww_18.sed
CFHTLS_SED/Ell_cww_19.sed
CFHTLS_SED/Ell_cww_20.sed
CFHTLS_SED/Sbc_cww.sed
CFHTLS_SED/Sbc_cww_0.sed
CFHTLS_SED/Sbc_cww_1.sed
CFHTLS_SED/Sbc_cww_2.sed
CFHTLS_SED/Sbc_cww_3.sed
CFHTLS_SED/Sbc_cww_4.sed
CFHTLS_SED/Sbc_cww_5.sed
CFHTLS_SED/Sbc_cww_6.sed
CFHTLS_SED/Sbc_cww_7.sed
CFHTLS_SED/Sbc_cww_8.sed
CFHTLS_SED/Sbc_cww_9.sed
CFHTLS_SED/Sbc_cww_10.sed
CFHTLS_SED/Sbc_cww_11.sed
CFHTLS_SED/Sbc_cww_12.sed
CFHTLS_SED/Sbc_cww_13.sed
CFHTLS_SED/Sbc_cww_14.sed
CFHTLS_SED/Sbc_cww_15.sed
CFHTLS_SED/Scd_cww.sed
CFHTLS_SED/Scd_cww_0.sed
CFHTLS_SED/Scd_cww_1.sed
CFHTLS_SED/Scd_cww_2.sed
CFHTLS_SED/Scd_cww_3.sed
CFHTLS_SED/Scd_cww_4.sed
CFHTLS_SED/Scd_cww_5.sed
CFHTLS_SED/Scd_cww_6.sed
CFHTLS_SED/Scd_cww_7.sed
CFHTLS_SED/Scd_cww_8.sed
CFHTLS_SED/Scd_cww_9.sed
CFHTLS_SED/Scd_cww_10.sed
CFHTLS_SED/Im_cww_0.sed
CFHTLS_SED/Im_cww_1.sed
CFHTLS_SED/Im_cww_2.sed
CFHTLS_SED/Im_cww_3.sed
CFHTLS_SED/Im_cww_4.sed
CFHTLS_SED/Im_cww_5.sed
CFHTLS_SED/Im_cww_6.sed
CFHTLS_SED/Im_cww_7.sed
CFHTLS_SED/Im_cww_8.sed
CFHTLS_SED/Im_cww_9.sed
CFHTLS_SED/Im_cww.sed
CFHTLS_SED/SB1.sed
CFHTLS_SED/SB1_1.sed
CFHTLS_SED/SB1_2.sed
CFHTLS_SED/SB2.sed
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.