meeting 2024 06 11 n547 - JacobPilawa/TriaxSchwarzschild_wiki_5 GitHub Wiki
- Summarizing the kinematic extraction, MGE fitting, and other input choices before running the first set of models.
-
Stats:
- D = 75.5 ± 3.3, Jensen+21 SBF
- Note that this galaxy is not listed in 2MASS but is designated the BCG of A194 with MV = −22.98 in Lauer et al. (2007a).
- ______ group members
- _e = _____ km/s
- Mk = _____, log(Mstar) = _____
- eps = _____, lambda = _____, though likely overestimated due to systemic velocities (Slow rotator)
- log10(Mhalo) = _____
- Angles
- PA GMOS = _____
- PA Mitchell = _____
- PA Phot = _____
- Psi GMOS = _____ (MASSIVE XIV)
- Psi Mitchell = _____ (MASSIVE X)
- Photometry: *
-
GMOS:
- All files seem to be on BOX.
- Has both symmetric and non-symmetric binning schemes.
- Other details:
- Files on BOX say +x axis of the GMOS data is at +94 deg. E of N.
- PSF listed as: avg. FWHM = 0.50" (Gaussian)
-
Mitchell:
- All Mitchell files for all galaxies are present! There are folded and unfolded versions for everything. I'm currently sticking with our usual unfolded data.
| Plot |
|---|
 |
-
Need to do this still, just had moved straight into final set of models.
-
Irina lists the settings on BOX as:
Nominal settings for pPXF fits to NGC 547 GMOS data
---------------------------------------------------
Wavelength range (rest): 8420-8770 A
Template list: barthlist.txt
(15 stars from Barth et al. 2002 -- spectra are from MILES Calcium Triplet library)
Spectrum mask: /GMOS/galaxy_files/N547/spectra/N547_spectra_mask.fits
(bright sky lines and chip gaps as assessed from inspection of GMOS spectra)
Initial velocity shift: 5300 km/s
Multiplicative polynomial: degree 3
Additive polynomial: degree 0 (constant offset)
Bias parameter: 0
| Bin Centers | Radial Moments | GMOS Maps | Mitchell Maps | Kinematic Comparison |
|---|---|---|---|---|
And the spectra and resulting fits look beautiful!
| Plot |
|---|
- Need to do this still, had moved straight to final set of diagnostics.
- Note these are the point-symmetrized versions of the kinematics.
| Bin Centers | Sym vs. Non-Sym Comparison | Radial Moments | GMOS Maps | Mitchell Maps |
|---|---|---|---|---|
 |
 |
 |
 |
 |
- See the bullet on 06/17 for the full set of MGE diagnostics.
- Here are all the input files used in the modeling and the various choices I've made:
aperture_gmos/aperture_mitchell.dat
- The only real number that I control in these files are the
kin*gridthetavalues, which represent the difference between the MGE +x axis and the +x axis of the kinematics. - These numbers are computed as PA_MGE - PA_kinematics:
- The MGE PA comes from the fits abiove, which uses
- The kinematic PAs come from the files on BOX.
- For GMOS, these values are listed in /MASSIVE/GMOS_kinematics/N0547/N0547_PA_note.txt, which lists the +x axis of the orbit modeling files as +94 deg. E of N.
- For Mitchell, the numbers are located in /MASSIVE/Mitchell_reduced_data/NGC0545/kinematics_paperversion/unfolded/more_files/NGC0545-s2-unfolded-bininfo.txt, which lists the gal pa as 94.1 deg E of N.
- This means that the input files should have
kin1gridthetaandkin2gridthetaas 3.94 and 0.3, respectively.
bins_gmos/bins_mitchell.dat
- I don't make any choices for these files. These are generated when preparing the orbit model files and associated the apertures with specific kinematics.
kindata.dat
- This is our normal kindata.dat file, which comes from the Gh8Mh6 minimization described above.
- This file contains 12 moments total, 8 of which are measured for the GMOS data and are from my own spectral fits, and 6 moments are fit for Mitchell (which I obtained from the files on BOX).
- I've currently left in all bins, but we can decide to mask the overlapping Mitchell bins if we'd like.
- I've added dummy moments for GMOS h9 through h12 and Mitchell h7 through h12 by assigning all the values to 0.0, and uncertainties identical to the h7/h8 errors for GMOS and h5/h6 errors for Mitchell.
mge.dat
- Currently a placeholder MGE so that we can discuss the fits above.
orbitstart_loop/box.dat
- These are just auxillary files which point to other files in our models. No changes or choices.
orblib_loop/box.in
-
The only specific changes I had to make here relate to the PSFs we assumed for each dataset:
- I obtained the PSF for GMOS from the file /MASSIVE/GMOS_kinematics/N0547/N547_GMOS_PSF.txt, which says that the average weighted PSF = 0.50". The code accepts the equivalent "sigma" for this PSF, which I obtain as sigma = 0.38/(2sqrt(2ln(2))) = 0.212.
- I am using the same PSF for Mitchell as we have used for the other galaxies, which has a sigma = 0.5, which is equivalent to a 1.177" PSF.
-
Other notable parameters that I kept constant from past orbit modeling:
- innermost mass bin starts at 0.1"
- outermost mass bin starts at 100" and extends to 316"
- 200 dynamical time integration for the box library; 2000 dynamical times for the loop library
parameters_loop/box.in
- Many of these parameters are changed on the fly (mass components and angles).
- I'm using:
- nEner = 40 for both
- nI2 = 18, 9 for the loop and box library respectively
- nI3 = 9, 9 for both
- logrmin = 10^-2 arcsec (minimum IC radius in arcsec)
- logrmin = 10^2.5 arcsec (maximum IC radius in arcsec)
- interplogmin = 10^-8
- interplogmax = 10^2.6
- gnfw halo with gamma = 0
triaxmass.in
- Auxillary file which points to other files.
triaxmassbin.in
- Auxillary file which points to other files.
- Also contains the same PSF information present in orblib_loop/box.dat.
triaxnnls.in
- Mostly points to other files.
- I've ensured that this file correctly assumes we have 12 total moments.
- I've also checked to make sure that this file reads in both the loop orbit library and box orbit library.