velocity determination in 2dimensions, achieved by running the 1dim fit once for each access and then combining velocities
EXECUTION STEPS:
Code execution steps: (assuming mirisim has been installed and setup, assuming jwst pipeline has been installed)
clone the directory on git and cd to it
Activate mirisim env and run the simulation file 'simulation.ini'
Exit mirisim env, go back to base or env containing jwst pipeline and other libraries imported in cell 1 of velocity-determination.ipynb
Update the' velocity-determination.ipynb' notebook
correct ippath (directory of outpit from mirisim) and 'ipname' (indicating the exact output to use) in cell2 based on mirisim output you want to use
In cell 4, update maskrefpath to a relevant refmask for the mirifov in your local miri crds, our simulations use ' … /crds/references/jwst/miri/jwst_miri_mask_0023.fits' and the 'do not use' flagged pixels where bit 0 is set are all filled with nan values
Run all cells till the text cell that says 'basic fit'
After that chose to run the cells for the velocities you wish to evaluate, a new velocity will produce a new simulation and hence take some time
PLEASE NOTE :
Changes to stella's files before use for simulation -
Change_pixeldq.py :
packaged in a function to accept two arguments, 1- input file path 2- output file path
closed the open hdul at the end
function_move_source.py :
added argument refmask which is the reference mask for the miri fov stacked to be the same size as the simulated science array from mirisim
move_source_code.py :
packaged the execution part of the code into a function that accepts arguments move(inputdir, outputdir, velin, refmask, stattomoving = True):
Velin is the multiple of 0.010996 arcsec/sec to be used for the simulation
theta multiplied by 1.25 to represent asteroid velocity direction at 45 degrees to x-axis
flag_move_source passed through the stattomoving argument and flipped, True now means move a stationary source
Accepts and propagates the refmask array to function_move_source
Closed the open hdul at the end
simulation files
Point.ini : cen changed to '0 6' to ensure no artefacts exist close to our object row/column
simulation.ini updated to change the name of the scene.ini file to point.ini
Except the addition of the refmask and correcting the scene file name, other changes do not improve the code in any way, but were intended to make it easier to use for our specific implementation