jacopo-chevallard / BEAGLE-general

Basic information to get started with the galaxy spectral modelling tool Beagle
https://www.iap.fr/beagle/
9 stars 0 forks source link

BEAGLE v0.10.5 hangs up on first object #43

Closed kevinhainline closed 7 years ago

kevinhainline commented 7 years ago

Hello! I'm trying to run BEAGLE v0.10.5 on some mock catalogue objects using the HST+NIRCam filters in our deep survey configuration, and the program seems to be hanging up after MultiNest generates the points:

$ BEAGLE --fit --parameter-file /home/kevinhainline/BEAGLE/BEAGLE-general-master/params/final_survey_photo_z/BEAGLE_photo_z_test.param

 Reading Filter File: /home/kevinhainline/BEAGLE/BEAGLE-general-master/filters/FILTERBIN.RES
 294 filters defined, out of 500 maximum     ...done
 ---> fixed                         sfh_type                                   0
 ---> fitted                        mass                                       0
 ---> fitted                        redshift                                   0
 ---> dependent                     nebular_logU                               0
 ---> fixed                         nebular_xi                                 0
 ---> dependent                     nebular_Z                                  0
 ---> fixed                         attenuation_type                           0
 ---> fitted                        tauV_eff                                   0
 ---> fixed                         mu                                         0
 ---> fitted                        tau                                        1
 ---> fitted                        metallicity                                1
 n_fitted:            5
 *****************************************************
 MultiNest v3.9
 Copyright Farhan Feroz & Mike Hobson
 Release Oct 2014

 no. of live points =  300
 dimensionality =    5
 *****************************************************
 Starting MultiNest
 generating live points
           ---------------------------------------------------------------------------------
           --- LINEAR:  X0 =  2.091E+02 is outside X range --- 9.000E+02 6.010E+04      2362
           ----------- Error reported only once, but it may occur more than once. ----------

...and it sits here.

I've tested the program using the fit photometry example code, and it runs fine. I've even resampled the filters that I'm using down to 400 lines each:

# It contains the filters to be used when calculating only photmetric models

# The possible units that can be used are: Jy, milliJy, microJy, nanoJy
# Alternatively one can pass a conversion factor to convert the fluxes in the
# catalogue into [erg s^-1 cm^-2 Hz^-1] -----> units:value:<conversion factor>
units:nanoJy

object_ID:colName:ID

index:279          min_rel_err:0.02       flux:colName:F435W              fluxerr:colName:F435W_err  label:F435W
index:280          min_rel_err:0.02       flux:colName:F606W              fluxerr:colName:F606W_err  label:F606W
index:281          min_rel_err:0.02       flux:colName:F775W              fluxerr:colName:F775W_err  label:F775W
index:282          min_rel_err:0.02       flux:colName:F814W              fluxerr:colName:F814W_err  label:F814W
index:283          min_rel_err:0.02       flux:colName:F850LP             fluxerr:colName:F850LP_err label:F850LP
index:285          min_rel_err:0.02       flux:colName:F090W              fluxerr:colName:F090W_err  label:F090W
index:286          min_rel_err:0.02       flux:colName:F115W              fluxerr:colName:F115W_err  label:F115W
index:287          min_rel_err:0.02       flux:colName:F150W              fluxerr:colName:F150W_err  label:F150W
index:288          min_rel_err:0.02       flux:colName:F200W              fluxerr:colName:F200W_err  label:F200W
index:289          min_rel_err:0.02       flux:colName:F277W              fluxerr:colName:F277W_err  label:F277W
index:290          min_rel_err:0.02       flux:colName:F335M              fluxerr:colName:F335M_err  label:F335M
index:291          min_rel_err:0.02       flux:colName:F356W              fluxerr:colName:F356W_err  label:F356W
index:292          min_rel_err:0.02       flux:colName:fake               fluxerr:colName:fake_err   label:fake
index:293          min_rel_err:0.02       flux:colName:F410M              fluxerr:colName:F410M_err  label:F410M
index:294          min_rel_err:0.02       flux:colName:F444W              fluxerr:colName:F444W_err  label:F444W

I've included the data file, as well as my filterfrm.res file:

BEAGLE_github_bug_3_21_17.zip

...and here is the parameter file:

# This is a stripped down version of the BEAGLE parameter file for 
# with checking the correct installation of BEAGLE using Docker.
# This file will be used for producing artificial photometry using
# an analytic star formation history.
# Verbose mode: 0 ---> no info
#               1 ----> print only errors
#               2 ---> print errors and warnings
#               3 ----> print errors, warnings and info 
VERBOSE = 1

# Seed for the random number generator
SEED = 994189

#****************************************************************
#******************** COSMOLOGICAL PARAMETERS *******************
#****************************************************************
#
# If left blank, default parameter values will be used.

HUBBLE_0 = 

OMEGA_MATTER =

OMEGA_LAMBDA =

#****************************************************************
# ******************* TEMPLATES  *****************************
#****************************************************************

#TEMPLATES = $BEAGLE_TEMPLATES/bc03/bc03_miles_chab_spectra
TEMPLATES NEBULAR = $BEAGLE_TEMPLATES/ineb_Jan16_C100/cb2013_n2_mup100_N015_O01_deplO70_C100_Jan16
EMISSION LINES CONFIGURATION = $BEAGLE_TEMPLATES/ineb_Jan16_C100/cb2013_n2_mup100_N015_O01_deplO70_C100_Jan16_line_wavelengths_PHOTOMETRY.dat

SHRINK TEMPLATES WL RANGE = 900 60000
REBIN TEMPLATES = 10

#SPECTRAL INDICES CONFIGURATION =

#SPECTRAL INDICES CATALOGUE = 

#****************************************************************
# ******************* PHOTOMETRY *****************************
#****************************************************************

FILTERS FILE = $BEAGLE_FILTERS/filter_NIRCam_DEEP.dat

PHOTOMETRIC CATALOGUE = $BEAGLE_DATA/final_survey_photo_z/BEAGLE_test_3_21_17.fits

FIT PHOTOMETRY = T
FIT SPECTRUM = F

#****************************************************************
#******************** DUST ATTENUATION AND IGM ABSORPTION *******
#****************************************************************
#
IGM ABSORPTION = Inoue

# File containing set of radiative transfer (RT) calculations of dust attenuation of starlight
#RT CALCULATIONS = Tuffs_bulge_NEW-28-11.dat Tuffs_thick_NEW-28-11.dat Tuffs_thin_NEW-28-11.dat
#RT COMPONENTS NAMES = bulge thick_disc thin_disc

# **************************************************
# ************** SF BIN #1 **************************
# **************************************************
#SF_PARAMETER  = name:sfh_type                 type:fixed  char_value:rising
SF_PARAMETER  = name:sfh_type                 type:fixed  char_value:delayed
#options for analytic SFHs are 'ssp' (for a burst), 'constant', 'exponential', 'delayed', 'rising'

# In log(yr)
SF_PARAMETER  = name:tau                      type:fitted  order_priority:1  prior:distribution:uniform  prior:range:[7.,10.5]

# in log(Z/Z_sun)
SF_PARAMETER  = name:metallicity              type:fitted  order_priority:1  prior:distribution:uniform  prior:range:[-2.2,0.24]

# In log M/M_sun
SF_PARAMETER  = name:mass                     type:fitted  order_priority:0  prior:distribution:uniform  prior:range:[5.,12.] 

# **************************************************
# **************  OTHER SF PARAMETERS **************************
# **************************************************

# In log(M_sun/yr)
#SF_PARAMETER  = name:sfr                type:fitted  order_priority:1  prior:distribution:uniform  prior:range:[-4.,4.]

#****************************************************************
# ******************* PRINTED OUTPUT  *****************************
#****************************************************************

RESULTS DIRECTORY = $BEAGLE_RESULTS/final_survey_photo_z/

#PRIORS CATALOGUE = $BEAGLE_DATA/final_survey_photo_z/BEAGLE_test_3_21_17.fits

#MOCK CATALOGUE NAME = 

PRINT PHOTOMETRY = T
PRINT SPECTRUM = T
PRINT SPECTRAL INDICES = F
PRINT SF AND ChE HISTORIES = F

#****************************************************************
# ******************* PARAMETERS HANDLING  *****************************
#****************************************************************

#PDF SAMPLER FILE = $BEAGLE_PARAM_DIR/MCMC_example.param
PDF SAMPLER FILE = $BEAGLE_PARAM_DIR/MCMC_new.param

#PARAMETER  = name:redshift     type:fitted    prior:distribution:dirac    prior:colName:redshift
PARAMETER  = name:redshift     type:fitted    order_priority:0    prior:distribution:uniform    prior:range:[0.,15.]

PARAMETER  = name:nebular_logU type:dependent

PARAMETER  = name:nebular_xi   type:fixed    value:0.3

PARAMETER  = name:nebular_Z    type:dependent 

PARAMETER  = name:attenuation_type  type:fixed  char_value:CF00   
# values can be: CF00, Calzetti, CCWW13_universal, CCWW13_full

PARAMETER  = name:tauV_eff     type:fitted  order_priority:0      prior:distribution:uniform   prior:range:[0.,2.]

PARAMETER  = name:mu           type:fixed   value:0.4
kevinhainline commented 7 years ago

Oh, wait! I realized that I am terribly dumb, and did not put the fluxes in nJy. Everything seems to be running now. Whoops, I'm sorry.

jacopo-chevallard commented 7 years ago

No problem, good that you could solve this so easily !!