Mock mode - gridded parameters not covering entire parameter range

[jacopo-chevallard]

migrating part of issue #47 here

@kevinhainline writes: (using Beagle version installed on orange, v0.10.5.)

I am having a new and related issue with trying to generate quiescent galaxies at a given stellar mass. Below, I show the parameter file that I'm using, which should produce galaxies with a mass of 10^8 solar masses, at a gridded range of redshifts between 0.25 and 4.0, in steps of 0.25.

BEAGLE is instead producing galaxies only out to z <= 1.75.

I can modify the parameter to file to go between z = 2.0 - 4.0, and get the remaining objects, but it seems strange that it doesn't go all the way between the requested limits.

In addition, I have tried to also grid in mass, and I get a similar redshift cutoff at z = 2.0, and if I try to go between z = 2.0 - 4.0, and a mass range of 10^5 - 10^8, in steps of 0.2 dex, I get this weird behavior:

I think that perhaps I should just try making a fits file and reading the parameters in from that?

# 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

#****************************************************************
# ******************* SPECTROSCOPY ***************************
#****************************************************************

#SPECTRAL INDICES CONFIGURATION = $BEAGLE_FILTERS/spec_config.dat

#SPECTRAL INDICES CATALOGUE = $BEAGLE_DATA/spec_config_17_EWs.fits

FIT SPECTRUM = F

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

FILTERS FILE = $BEAGLE_FILTERS/filter_NIRCam_DEEP.dat

FIT PHOTOMETRY = F
FIT SPECTRUM = F

# N_OBJECTS = 10

#****************************************************************
#******************** 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                  mock:type:fixed  char_value:exponential
#options for analytic SFHs are 'ssp' (for a burst), 'constant', 'exponential', 'delayed', 'rising'

# In log(yr)
SF_PARAMETER  = name:tau                       mock:type:fixed    mock:value:6.0

# in log(Z/Z_sun)
SF_PARAMETER  = name:metallicity               mock:type:fixed    mock:value:-0.5

SF_PARAMETER  = name:mass                      mock:type:fixed    mock:value:8.0

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

SF_PARAMETER  = name:max_stellar_age   mock:type:fixed    mock:value:9.5

SF_PARAMETER = name:formation_redshift mock:type:fixed    mock:value:13.0

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

RESULTS DIRECTORY = $BEAGLE_RESULTS/BEAGLE_mock_lowz_lowmass/

MOCK CATALOGUE NAME = BEAGLE_mock_lowz_lowmass.fits

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                mock:type:gridded  mock:grid:range:[0.25,4.0]  mock:grid:step:0.25

PARAMETER  = name:nebular_logU            mock:type:fixed    mock:value:-2.0

PARAMETER  = name:nebular_xi            mock:type:fixed    mock:value:0.3

PARAMETER  = name:nebular_Z      type:dependent 

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

PARAMETER  = name:tauV_eff      mock:type:fixed    mock:value:0.0

PARAMETER  = name:mu            mock:type:fixed    mock:value:0.4

[jacopo-chevallard]

@kevinhainline the problem in this case is due to the following lines in the parameter file

SF_PARAMETER  = name:max_stellar_age   mock:type:fixed    mock:value:9.5
SF_PARAMETER = name:formation_redshift mock:type:fixed    mock:value:13.0

since you're fixing the age of stars to log(t/yr)=9.5, but the time elapsed between z=13 (formation redshift) and any z >= 2.0 is < log(t/yr)=9.5. Beagle is therefore skipping all unphysical combinations of parameters, therefore repeating the allowed combinations, i.e. redshifts values from 0.25 to 1.75

If you comment out the line

SF_PARAMETER  = name:max_stellar_age   mock:type:fixed    mock:value:9.5

then you'll get the entire redshift and mass grid

[jacopo-chevallard]

Solved with above commit, see kevin's tests in #47