# Configuration for pyLINEAR

cpu: &CPU
  ncpu: 1                                # number of processors to use
  verbose: True
  
calib: &CALIB
  h5conf: hst_wfc3_ir_beams.h5           # HDF5 file for beam info
  h5flat: Null                           # HDF5 file for flat field
  h5siaf: hst_wfc3_ir_detector.h5        # HDF5 file for detectors
tables: &TABLES
  path: TABLES2/                          # relative path for the tables
  remake: False
  gzip: True                             # gzip the resulting data
  nsub: 5                                # number of intervals for TrapZ
  cpu: *CPU
                                         # import CPU settings
sources:
  segmap: seg.fits                       # segmentation map for sources
  obslst: obs.lst                        # file describing the bands
  detindex: 0                            # detection index in obslst
  maglim: 26.                            # magnitude limit of sources

modules:
  simulation:
    perform: True                        # flag to run the simulation
    beam: A                              # beams to simulate
    grism: G102                          # which grism to simulate
    imgtype: wcs                         # where WCS comes from: wcs or img
    gzip: True                           # gzip the outputs
    imglst: wcs.lst
    sedlst: sed.lst
    noise:
      perform: False
      exptime: 700.
      skyrate: 0.3
    cpu: *CPU
    calib: *CALIB
    tables: *TABLES

  extraction:
    perform: True
    mode: fixed
    group: False
    logdamp: -3.45161   #[-6., -1., 0.1]
    outroot: test2
    imglst: flt.lst
    beam: A
    mask: Null
    dlamb: 25.
    initialize: 'broadband'    # can be direct or broadband
    lsqr:
      atol: 1e-10          # a tol
      btol: 1e-10
      conlim: 1e10
      maxiter: 1e5
      show: False
    cpu: *CPU
    calib: *CALIB
    tables: *TABLES
    mcmcunc:
      perform: False
      nwalkers: 100
      nstep: 100
      burn: 0.2
      cpu: *CPU
    residuals:
      perform: True
      gzip: False

  cutout:
    perform: False
    beam: A
    #grism: G102
    imglst: flt.lst
    padding: [0, 2]
    cpu: *CPU
    calib: *CALIB