I redid the galfit simulations and found quite a big offset between input and recovered values for the Hercules field of view. I tried using the default mips prf for detection, and this worked much better. So, I now have to rerun galfit on Hercules, and this is a good opportunity to document the procedure using a python notebook.
to run on 24um images with PA and BA held fixed at SDSS values: (2015/02/03)
mkw11.run_just_get_images24(make_mask=1,review_mask=1) # if you haven't made cutouts yet
mkw11.rungalfit_first_time_24(fitBAflag=False,fitPAflag=False)
mkw11.write_galfit_sersic_parameters_24()
In [5]:
import numpy as np
from pylab import *
%matplotlib inline
import warnings
warnings.filterwarnings('ignore')
In [6]:
%run ~/Dropbox/pythonCode/LCSrungalfit.py
got here
Running on Rose's mac pro
got to the end!
Auto-detected table type: fits
Auto-detected table type: fits
Auto-detected table type: fits
Auto-detected table type: fits
Auto-detected table type: fits
conversion from F24 to L24 = 88772.9837327
Hercules : N members = 125 N field = 0
N spirals on 24um image = 36 Nspiral members = 36
N galaxies on 24um image = 125 N members on 24um image = 125
0.03639 458.5295
could not find 142589 in coma
could not find 104115 in coma
could not find 104020 in coma
could not find 104022 in coma
could not find 142662 in coma
could not find 142763 in coma
could not find 142797 in coma
could not find 162768 in coma
could not find 162797 in coma
In [7]:
herc.rungalfit_first_time_24(fitBAflag=False,fitPAflag=False)
# running from ipython window because the interactions are a bit painful
#################
just inside rungalfit
conv, fitBA, fitPA flag = True False False
0
here is what I think the start_index is: 0
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99037
snr24 = 4.66
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
y
running galfit again - hope that is what you want!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-fixedBA-1Comp-galfit-out.fits
convflag, fitBA, fitPA = True False False
###############
just inside keep_BA_fixed
fitBAflag, fitPAflag, convflag = 0 0 True
deleting outcoords if it exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/incoords /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um
Killing IRAF task `imgets'
42.0 41.0
xcenter, ycenter = 21.0 20.5
objIndex = (array([1]),) [ 8.37949521 1.31957304 12.68237967 22.62983131 23.40793447
26.0046133 6.16183414 20.50573639 13.44128126 5.37363648
10.60950786 4.30767501 13.53395567 13.76036064 12.2887608
6.09187295 8.00000025 7.87201886 19.63018991 11.24933456
13.86231124 14.92190182 8.85591492 12.14523841 15.81252681
17.10447339 1.70213777 11.3944282 15.23925044 14.41952898
25.02335879 26.23500326 10.94093031 3.45950864 14.77687524
14.44749961 7.11252866 10.09859678 10.19764698 7.88242862
8.05573808 16.5444756 7.33588638 16.55379123 5.55409219]
object number = 2
elliptical = False
spiral = True
agn = False
hit any key when ready to continue (q to quit)
%%%%%%%%%%%%%%%%%%%%%%%%%%
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24
1 42.0 1 41.0 42.0
%%%%%%%%%%%%%%%%%%%%%%%%%%
just before calling igalfit
fitBA = 0
%%%%%%%%%%%%%%%%%
inside galfit
1 42.0 1 41.0 42.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
%%%%%%%%%%%%%%%%%
conv, fitBA,fitPA = True 0 0
***%%%%%%%%%%%%%%%%%
inside galfit class
1 42.0 1 41.0 42.0
1 42.0 1 41.0 42.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
self.fitall = 0
***%%%%%%%%%%%%%%%%%
sersic n, fitsersicn = 0.794217 1
BA, fitBA = 0.72308 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-fixedBAgalfit.input.1Comp
appears like galfit did not complete
Are you happy with the fit? n=no x=quit y (or any other key) = yes
y
this is what I think you typed y
running for mips galaxies!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-1Comp-galfit-out.fits
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-1Comp-galfit-out.fits
2
3
4
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-1Comp-subcomps.fits
-----------------------------------------------------------------------------
Input image : Hercules-99037-galfit-cutout24.fits[1:41,1:41]
Init. par. file : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/galfit.input.1Comp
Restart file : galfit.01
Output image : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-1Comp-galfit-out.fits
sersic : ( 22.57, 20.91) 14.22 6.40 1.46 0.54 7.38
( 0.10, 0.12) 0.18 1.40 0.35 0.04 4.16
sky : [ 21.00, 21.00] -2.66e-02 [0.00e+00] [0.00e+00]
4.58e-03 [0.00e+00] [0.00e+00]
Chi^2 = 85063.27891, ndof = 1568
Chi^2/nu = 54.250
-----------------------------------------------------------------------------
here is what I have for diskflag[i] = [1 0 0 0 0 0]
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-1Comp-subcomps.fits[2] -> /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-diskimage-1Comp.fits
model_mag = [ 14.22]
total mag = 14.22
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-diskimage-1Comp.fits,MAG_TOTAL = 14.22
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-diskimage-1Comp.fits updated
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-diskimage-1Comp.fits,COMPONENT_FLAG = 1
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99037-24-diskimage-1Comp.fits updated
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99056
snr24 = 17.3
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
y
running galfit again - hope that is what you want!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-fixedBA-1Comp-galfit-out.fits
convflag, fitBA, fitPA = True False False
###############
just inside keep_BA_fixed
fitBAflag, fitPAflag, convflag = 0 0 True
deleting outcoords if it exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/incoords /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um
xcenter, ycenter = 34.0 34.0
objIndex = (array([1]),) [ 8.93721847 0.61641707 22.08204904 27.37036648 27.95343231
30.52390958 40.98853032 44.66779136 43.65021022 34.02718056
31.31714644 38.85751222 34.6812128 36.69162599 33.7192675
43.0436802 7.00257781 27.15807742 17.75019572 28.7057674
32.00005077 33.74888326 24.72361869 27.38469419 22.40509594
24.41158309 31.42826475 27.30798129 24.58043387 18.03419031
24.75647863 16.7270533 17.73245663 21.70587536 35.42301129
32.26756925 30.36507705 32.18019703 26.62949442 17.89673534
34.64376244 37.40800301 31.6799932 13.57687343]
object number = 2
elliptical = False
spiral = False
agn = False
hit any key when ready to continue (q to quit)
%%%%%%%%%%%%%%%%%%%%%%%%%%
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24
1 68.0 1 68.0 68.0
%%%%%%%%%%%%%%%%%%%%%%%%%%
just before calling igalfit
fitBA = 0
%%%%%%%%%%%%%%%%%
inside galfit
1 68.0 1 68.0 68.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
%%%%%%%%%%%%%%%%%
conv, fitBA,fitPA = True 0 0
***%%%%%%%%%%%%%%%%%
inside galfit class
1 68.0 1 68.0 68.0
1 68.0 1 68.0 68.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
self.fitall = 0
***%%%%%%%%%%%%%%%%%
sersic n, fitsersicn = 1.81101 1
BA, fitBA = 0.738283 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-fixedBAgalfit.input.1Comp
%%%%%%%%%%%%%%%%%%
inside display_results
self.galfit_flag = 1
1_XC : 34.89 +/- 0.04
1_YC : 34.58 +/- 0.03
1_MAG : 12.74 +/- 0.02
1_RE : 3.33 +/- 0.10
1_N : 2.04 +/- 0.14
1_AR : 0.74 +/- 0.00
1_PA : -67.66 +/- 0.00
2_SKY : -0.01 +/- 0.00
ERROR : 0.00
CHI2NU : 1.33
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-fixedBA-1Comp-subcomps.fits
frame delete 8
couldn't execute the following ds9 command : frame delete 8
Are you happy with the fit? n=no x=quit y (or any other key) = yes
y
this is what I think you typed y
running for mips galaxies!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-1Comp-galfit-out.fits
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-1Comp-galfit-out.fits
2
3
4
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-1Comp-subcomps.fits
-----------------------------------------------------------------------------
Input image : Hercules-99056-galfit-cutout24.fits[1:89,1:89]
Init. par. file : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/galfit.input.1Comp
Restart file : galfit.01
Output image : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-1Comp-galfit-out.fits
sersic : ( 45.76, 45.53) 12.81 4.38 1.18 0.87 -80.21
( 0.01, 0.01) 0.01 0.07 0.02 0.01 2.62
sky : [ 45.00, 45.00] -9.86e-03 [0.00e+00] [0.00e+00]
8.44e-04 [0.00e+00] [0.00e+00]
Chi^2 = 36782.89046, ndof = 7291
Chi^2/nu = 5.045
-----------------------------------------------------------------------------
here is what I have for diskflag[i] = [1 0 0 0 0 0]
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-1Comp-subcomps.fits[2] -> /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-diskimage-1Comp.fits
model_mag = [ 12.81]
total mag = 12.81
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-diskimage-1Comp.fits,MAG_TOTAL = 12.81
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-diskimage-1Comp.fits updated
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-diskimage-1Comp.fits,COMPONENT_FLAG = 1
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99056-24-diskimage-1Comp.fits updated
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99059
snr24 = 18.9
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
running galfit again - hope that is what you want!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-fixedBA-1Comp-galfit-out.fits
convflag, fitBA, fitPA = True False False
###############
just inside keep_BA_fixed
fitBAflag, fitPAflag, convflag = 0 0 True
deleting outcoords if it exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/incoords /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um
xcenter, ycenter = 37.5 37.5
objIndex = (array([0]),) [ 1.19364526 18.26750834 30.71463833 29.79103103 31.68623728
39.86337764 48.78567702 37.04171262 38.99322485 41.77864701
36.53737681 24.52370692 37.18401904 21.80085047 21.51677713
19.72790156 17.50748823 38.35408397 9.36855725 39.59407981
16.03188685 34.88079128 37.74898946 20.14175762 34.19422103
16.99749446 33.95068201 38.45831667 21.96999863 37.67173541
35.22233509 40.21248961 18.88754727 30.80554892 21.46500524
25.44850033 24.52848503 43.28186967 35.09801248 30.98443196
42.4192702 34.9084807 43.06571686 46.4382455 29.87393235
33.81384156 26.58252526 26.72598378 34.28402527 42.62235212]
object number = 1
elliptical = False
spiral = True
agn = True
hit any key when ready to continue (q to quit)
%%%%%%%%%%%%%%%%%%%%%%%%%%
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24
1 75.0 1 75.0 75.0
%%%%%%%%%%%%%%%%%%%%%%%%%%
just before calling igalfit
fitBA = 0
%%%%%%%%%%%%%%%%%
inside galfit
1 75.0 1 75.0 75.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
%%%%%%%%%%%%%%%%%
conv, fitBA,fitPA = True 0 0
***%%%%%%%%%%%%%%%%%
inside galfit class
1 75.0 1 75.0 75.0
1 75.0 1 75.0 75.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
self.fitall = 0
***%%%%%%%%%%%%%%%%%
sersic n, fitsersicn = 2.53406 1
BA, fitBA = 0.233371 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-fixedBAgalfit.input.1Comp
%%%%%%%%%%%%%%%%%%
inside display_results
self.galfit_flag = 1
1_XC : 36.94 +/- 0.03
1_YC : 38.92 +/- 0.02
1_MAG : 11.01 +/- 0.01
1_RE : 7.18 +/- 0.04
1_N : 0.20 +/- 0.01
1_AR : 0.23 +/- 0.00
1_PA : 75.88 +/- 0.00
2_SKY : -0.01 +/- 0.00
ERROR : 0.00
CHI2NU : 2.35
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-fixedBA-1Comp-subcomps.fits
frame delete 8
couldn't execute the following ds9 command : frame delete 8
Are you happy with the fit? n=no x=quit y (or any other key) = yes
n
this is what I think you typed n
CURRENT INPUTS:
mag = 12.73 1
Re = 7.31 1
n = 2.53 1
B/A = 0.23 0
PA = 75.88 0
fitall = 0
fitcenter = 1
What is wrong?
n = adjust sersic
r = reset Re
o = nearby object (toggle fitall)
b = B/A
p = PA
m = mag
c = recenter
f = hold values fixed
a = toggle asymmetry parameter
R = reset to original values
g = go (run galfit)
x=quit
r
radius = 10
CURRENT INPUTS:
mag = 12.73 1
Re = 10.00 1
n = 2.53 1
B/A = 0.23 0
PA = 75.88 0
fitall = 0
fitcenter = 1
What is wrong?
n = adjust sersic
r = reset Re
o = nearby object (toggle fitall)
b = B/A
p = PA
m = mag
c = recenter
f = hold values fixed
a = toggle asymmetry parameter
R = reset to original values
g = go (run galfit)
x=quit
g
sersic n, fitsersicn = 2.53406 1
BA, fitBA = 0.233371 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-fixedBAgalfit.input.1Comp
%%%%%%%%%%%%%%%%%%
inside display_results
self.galfit_flag = 1
1_XC : 36.94 +/- 0.03
1_YC : 38.92 +/- 0.02
1_MAG : 11.01 +/- 0.01
1_RE : 7.18 +/- 0.04
1_N : 0.20 +/- 0.01
1_AR : 0.23 +/- 0.00
1_PA : 75.88 +/- 0.00
2_SKY : -0.01 +/- 0.00
ERROR : 0.00
CHI2NU : 2.35
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-fixedBA-1Comp-subcomps.fits
frame delete 8
couldn't execute the following ds9 command : frame delete 8
Are you happy with the fit? n=no x=quit y (or any other key) = yes
this is what I think you typed
not sure what you entered, so setting repeat_flag = False
running for mips galaxies!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-1Comp-galfit-out.fits
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-1Comp-galfit-out.fits
2
3
4
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-1Comp-subcomps.fits
-----------------------------------------------------------------------------
Input image : Hercules-99059-galfit-cutout24.fits[1:127,1:127]
Init. par. file : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/galfit.input.1Comp
Restart file : galfit.01
Output image : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-1Comp-galfit-out.fits
sersic : ( 63.63, 64.61) 11.15 6.41 0.43 0.15 71.08
( 0.02, 0.01) 0.00 0.02 0.01 0.00 0.09
sky : [ 64.00, 64.00] -5.55e-03 [0.00e+00] [0.00e+00]
9.44e-04 [0.00e+00] [0.00e+00]
Chi^2 = 200987.09727, ndof = 15149
Chi^2/nu = 13.267
-----------------------------------------------------------------------------
here is what I have for diskflag[i] = [1 0 0 0 0 0]
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-1Comp-subcomps.fits[2] -> /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-diskimage-1Comp.fits
model_mag = [ 11.15]
total mag = 11.15
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-diskimage-1Comp.fits,MAG_TOTAL = 11.15
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-diskimage-1Comp.fits updated
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-diskimage-1Comp.fits,COMPONENT_FLAG = 1
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99059-24-diskimage-1Comp.fits updated
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99060
snr24 = 7.71
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
y
running galfit again - hope that is what you want!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-fixedBA-1Comp-galfit-out.fits
convflag, fitBA, fitPA = True False False
###############
just inside keep_BA_fixed
fitBAflag, fitPAflag, convflag = 0 0 True
deleting outcoords if it exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/incoords /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um
xcenter, ycenter = 19.0 19.0
objIndex = (array([1]),) [ 19.35889493 1.17867468 22.91042656 24.35615865 19.88218137
18.17471309 11.307555 16.86339447 23.04838606 10.22855689
15.35771601 11.01473785 4.94181283 17.51610291 18.26857685
7.93674373]
object number = 2
elliptical = False
spiral = True
agn = False
hit any key when ready to continue (q to quit)
%%%%%%%%%%%%%%%%%%%%%%%%%%
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24
1 38.0 1 38.0 38.0
%%%%%%%%%%%%%%%%%%%%%%%%%%
just before calling igalfit
fitBA = 0
%%%%%%%%%%%%%%%%%
inside galfit
1 38.0 1 38.0 38.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
%%%%%%%%%%%%%%%%%
conv, fitBA,fitPA = True 0 0
***%%%%%%%%%%%%%%%%%
inside galfit class
1 38.0 1 38.0 38.0
1 38.0 1 38.0 38.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
self.fitall = 0
***%%%%%%%%%%%%%%%%%
sersic n, fitsersicn = 1.34276 1
BA, fitBA = 0.91022 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-fixedBAgalfit.input.1Comp
appears like galfit did not complete
Are you happy with the fit? n=no x=quit y (or any other key) = yes
n
this is what I think you typed n
CURRENT INPUTS:
mag = 16.05 1
Re = 1.72 1
n = 1.34 1
B/A = 0.91 0
PA = 139.66 0
fitall = 0
fitcenter = 1
What is wrong?
n = adjust sersic
r = reset Re
o = nearby object (toggle fitall)
b = B/A
p = PA
m = mag
c = recenter
f = hold values fixed
a = toggle asymmetry parameter
R = reset to original values
g = go (run galfit)
x=quit
r
radius = 3
CURRENT INPUTS:
mag = 16.05 1
Re = 3.00 1
n = 1.34 1
B/A = 0.91 0
PA = 139.66 0
fitall = 0
fitcenter = 1
What is wrong?
n = adjust sersic
r = reset Re
o = nearby object (toggle fitall)
b = B/A
p = PA
m = mag
c = recenter
f = hold values fixed
a = toggle asymmetry parameter
R = reset to original values
g = go (run galfit)
x=quit
g
sersic n, fitsersicn = 1.34276 1
BA, fitBA = 0.91022 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-fixedBAgalfit.input.1Comp
%%%%%%%%%%%%%%%%%%
inside display_results
self.galfit_flag = 1
1_XC : 19.57 +/- 0.04
1_YC : 19.75 +/- 0.05
1_MAG : 15.42 +/- 0.03
1_RE : 1.23 +/- 0.09
1_N : 0.07 +/- 1.00
1_AR : 0.91 +/- 0.00
1_PA : -40.34 +/- 0.00
2_SKY : -0.00 +/- 0.00
ERROR : 0.00
CHI2NU : 1.17
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-fixedBA-1Comp-subcomps.fits
frame delete 8
couldn't execute the following ds9 command : frame delete 8
Are you happy with the fit? n=no x=quit y (or any other key) = yes
this is what I think you typed
not sure what you entered, so setting repeat_flag = False
running for mips galaxies!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-1Comp-galfit-out.fits
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-1Comp-galfit-out.fits
2
3
4
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-1Comp-subcomps.fits
-----------------------------------------------------------------------------
Input image : Hercules-99060-galfit-cutout24.fits[1:41,1:41]
Init. par. file : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/galfit.input.1Comp
Restart file : galfit.01
Output image : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-1Comp-galfit-out.fits
sersic : ( 21.61, 21.72) 15.61 2.12 0.69 0.88 -52.34
( 0.06, 0.04) 0.05 0.11 0.06 0.04 14.18
sky : [ 21.00, 21.00] -4.44e-04 [0.00e+00] [0.00e+00]
1.64e-03 [0.00e+00] [0.00e+00]
Chi^2 = 6421.99812, ndof = 1560
Chi^2/nu = 4.117
-----------------------------------------------------------------------------
here is what I have for diskflag[i] = [1 0 0 0 0 0]
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-1Comp-subcomps.fits[2] -> /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-diskimage-1Comp.fits
model_mag = [ 15.61]
total mag = 15.61
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-diskimage-1Comp.fits,MAG_TOTAL = 15.61
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-diskimage-1Comp.fits updated
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-diskimage-1Comp.fits,COMPONENT_FLAG = 1
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99060-24-diskimage-1Comp.fits updated
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99084
snr24 = 9.44
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
y
running galfit again - hope that is what you want!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-fixedBA-1Comp-galfit-out.fits
convflag, fitBA, fitPA = True False False
###############
just inside keep_BA_fixed
fitBAflag, fitPAflag, convflag = 0 0 True
deleting outcoords if it exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/incoords /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um
xcenter, ycenter = 23.0 22.5
objIndex = (array([1]),) [ 22.53331622 0.92917706 12.54475847 26.96150146 22.1770512
21.91247978 11.78912384 22.39392116 9.20316782 21.01940727
12.84346406 25.62917098 17.86893438 13.78489409 21.48059499
12.58684138 19.14849542 17.56023348]
object number = 2
elliptical = True
spiral = False
agn = False
hit any key when ready to continue (q to quit)
%%%%%%%%%%%%%%%%%%%%%%%%%%
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24
1 46.0 1 45.0 46.0
%%%%%%%%%%%%%%%%%%%%%%%%%%
just before calling igalfit
fitBA = 0
%%%%%%%%%%%%%%%%%
inside galfit
1 46.0 1 45.0 46.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
%%%%%%%%%%%%%%%%%
conv, fitBA,fitPA = True 0 0
***%%%%%%%%%%%%%%%%%
inside galfit class
1 46.0 1 45.0 46.0
1 46.0 1 45.0 46.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
self.fitall = 0
***%%%%%%%%%%%%%%%%%
sersic n, fitsersicn = 5.96805 1
BA, fitBA = 0.756128 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-fixedBAgalfit.input.1Comp
%%%%%%%%%%%%%%%%%%
inside display_results
self.galfit_flag = 1
1_XC : 23.50 +/- 0.06
1_YC : 23.39 +/- 0.06
1_MAG : 14.86 +/- 0.75
1_RE : 0.10 +/- 0.27
1_N : 12.94 +/- 29.21
1_AR : 0.76 +/- 0.00
1_PA : 59.53 +/- 0.00
2_SKY : -0.00 +/- 0.00
ERROR : 1.00
CHI2NU : 1.18
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-fixedBA-1Comp-subcomps.fits
frame delete 8
couldn't execute the following ds9 command : frame delete 8
Are you happy with the fit? n=no x=quit y (or any other key) = yes
this is what I think you typed
not sure what you entered, so setting repeat_flag = False
running for mips galaxies!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-1Comp-galfit-out.fits
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-1Comp-galfit-out.fits
2
3
4
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-1Comp-subcomps.fits
-----------------------------------------------------------------------------
Input image : Hercules-99084-galfit-cutout24.fits[1:59,1:59]
Init. par. file : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/galfit.input.1Comp
Restart file : galfit.01
Output image : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-1Comp-galfit-out.fits
sersic : ( 30.54, 30.45) 14.90 2.65 1.89 0.94 -88.39
( 0.02, 0.02) 0.07 0.26 0.19 0.04 20.45
sky : [ 30.00, 30.00] -5.50e-03 [0.00e+00] [0.00e+00]
1.19e-03 [0.00e+00] [0.00e+00]
Chi^2 = 15080.68030, ndof = 3220
Chi^2/nu = 4.683
-----------------------------------------------------------------------------
here is what I have for diskflag[i] = [1 0 0 0 0 0]
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-1Comp-subcomps.fits[2] -> /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-diskimage-1Comp.fits
model_mag = [ 14.9]
total mag = 14.9
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-diskimage-1Comp.fits,MAG_TOTAL = 14.9
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-diskimage-1Comp.fits updated
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-diskimage-1Comp.fits,COMPONENT_FLAG = 1
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99084-24-diskimage-1Comp.fits updated
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99478
snr24 = 5.82
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
y
running galfit again - hope that is what you want!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-fixedBA-1Comp-galfit-out.fits
convflag, fitBA, fitPA = True False False
###############
just inside keep_BA_fixed
fitBAflag, fitPAflag, convflag = 0 0 True
deleting outcoords if it exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/incoords /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um
Killing IRAF task `imgets'
38.0 38.0
xcenter, ycenter = 19.0 19.0
objIndex = (array([ 0, 15]),) [ 1.6308602 19.803814 19.07520155 13.59383128 25.37793262
17.88048872 17.64716952 13.5576845 10.47532577 5.0141404
23.3362755 12.26940296 13.22743951 21.49564628 22.15836763
1.6308602 19.803814 19.07520155 13.59383128 25.37793262
17.88048872 22.15836763 17.64716952 13.5576845 10.47532577
5.0141404 23.3362755 14.66337042 13.22743951 12.26940296
21.49564628]
object number = 1
elliptical = False
spiral = True
agn = False
hit any key when ready to continue (q to quit)
%%%%%%%%%%%%%%%%%%%%%%%%%%
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24
1 38.0 1 38.0 38.0
%%%%%%%%%%%%%%%%%%%%%%%%%%
just before calling igalfit
fitBA = 0
%%%%%%%%%%%%%%%%%
inside galfit
1 38.0 1 38.0 38.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
%%%%%%%%%%%%%%%%%
conv, fitBA,fitPA = True 0 0
***%%%%%%%%%%%%%%%%%
inside galfit class
1 38.0 1 38.0 38.0
1 38.0 1 38.0 38.0
psf_image = /Users/rfinn/research/LocalClusters/PRF/mips24_prf_mosaic_2.45_4x.fits
self.fitall = 0
***%%%%%%%%%%%%%%%%%
sersic n, fitsersicn = 1.12914 1
BA, fitBA = 0.355501 0
run the following: galfit /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-fixedBAgalfit.input.1Comp
%%%%%%%%%%%%%%%%%%
inside display_results
self.galfit_flag = 1
1_XC : 19.82 +/- 0.08
1_YC : 19.70 +/- 0.08
1_MAG : 15.60 +/- 0.06
1_RE : 2.06 +/- 4.44
1_N : 0.03 +/- 0.27
1_AR : 0.36 +/- 0.00
1_PA : -45.21 +/- 0.00
2_SKY : -0.00 +/- 0.00
ERROR : 1.00
CHI2NU : 1.27
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-fixedBA-1Comp-subcomps.fits
frame delete 8
couldn't execute the following ds9 command : frame delete 8
Are you happy with the fit? n=no x=quit y (or any other key) = yes
y
this is what I think you typed y
running for mips galaxies!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-1Comp-galfit-out.fits
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-1Comp-galfit-out.fits
2
3
4
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-1Comp-subcomps.fits
-----------------------------------------------------------------------------
Input image : Hercules-99478-galfit-cutout24.fits[1:41,1:41]
Init. par. file : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/galfit.input.1Comp
Restart file : galfit.01
Output image : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-1Comp-galfit-out.fits
sersic : ( 21.80, 21.82) 15.79 2.44 0.62 0.68 -20.69
( 0.06, 0.10) 0.05 0.15 0.08 0.03 5.95
sky : [ 21.00, 21.00] -4.50e-03 [0.00e+00] [0.00e+00]
1.67e-03 [0.00e+00] [0.00e+00]
Chi^2 = 6897.92378, ndof = 1594
Chi^2/nu = 4.327
-----------------------------------------------------------------------------
here is what I have for diskflag[i] = [1 0 0 0 0 0]
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-1Comp-subcomps.fits[2] -> /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-diskimage-1Comp.fits
model_mag = [ 15.79]
total mag = 15.79
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-diskimage-1Comp.fits,MAG_TOTAL = 15.79
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-diskimage-1Comp.fits updated
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-diskimage-1Comp.fits,COMPONENT_FLAG = 1
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99478-24-diskimage-1Comp.fits updated
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99495
snr24 = 12.2
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
y
running galfit again - hope that is what you want!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-fixedBA-1Comp-galfit-out.fits
convflag, fitBA, fitPA = True False False
###############
just inside keep_BA_fixed
fitBAflag, fitPAflag, convflag = 0 0 True
deleting outcoords if it exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/incoords /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um
Killing IRAF task `imgets'
54.0 54.0
xcenter, ycenter = 27.0 27.0
objIndex = (array([0]),) [ 0.42264051 14.30658813 27.88930003 32.37833685 18.17930717
35.08993203 25.23378935 17.63507553 17.77039381 28.0298127
27.94926949 15.66737298 26.21701526 11.31329894 24.49278261
18.16511186 35.79668187 11.32897987 19.49959225 13.5246631
22.05039909 16.25585387 22.05162271 26.8177882 21.51412341
25.70788828 10.01320109 30.17947491 24.22712942 32.72310804
8.80569889 23.55134752 33.44508634 17.69785569 13.21947715
8.69762801 11.41504661 11.5292328 10.60662585 23.95046524
6.65355431 9.583 31.27680011 30.5274222 26.30688309
17.86326387]
object number = 1
elliptical = False
spiral = False
agn = False
hit any key when ready to continue (q to quit)
q
running for mips galaxies!
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-1Comp-galfit-out.fits
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-1Comp-galfit-out.fits
2
3
4
file to display = /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-1Comp-subcomps.fits
-----------------------------------------------------------------------------
Input image : Hercules-99495-galfit-cutout24.fits[1:53,1:53]
Init. par. file : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/galfit.input.1Comp
Restart file : galfit.01
Output image : /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-1Comp-galfit-out.fits
sersic : ( 27.21, 27.80) 14.91 2.16 1.15 0.86 25.23
( 0.02, 0.03) 0.03 0.09 0.06 0.02 7.24
sky : [ 27.00, 27.00] -5.67e-03 [0.00e+00] [0.00e+00]
1.24e-03 [0.00e+00] [0.00e+00]
Chi^2 = 11262.44860, ndof = 2631
Chi^2/nu = 4.281
-----------------------------------------------------------------------------
here is what I have for diskflag[i] = [1 0 0 0 0 0]
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-1Comp-subcomps.fits[2] -> /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-diskimage-1Comp.fits
model_mag = [ 14.91]
total mag = 14.91
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-diskimage-1Comp.fits,MAG_TOTAL = 14.91
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-diskimage-1Comp.fits updated
add /Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-diskimage-1Comp.fits,COMPONENT_FLAG = 1
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99495-24-diskimage-1Comp.fits updated
*******************
*******************
*******************
running galfit analysis for MIPS image of Hercules - 99499
snr24 = 3.61
output image already exists
/Users/rfinn/research/LocalClusters/GalfitAnalysis/Hercules/24um/Hercules-99499-24-fixedBA-1Comp-galfit-out.fits
Run galfit again? y (or any key)=yes n=no x=quit
x
I think you want to quit
When finished, run
LCSmergespiralcats.py
This made a sizeable difference for Hercules. I also redid the galaxies in A2052, but this made much less difference, and so I went with using the point source from the scan rather than the mips PRF.
I will need to rerun A2052 using the correct PRF, but not today!
In [ ]:
Content source: rfinn/LCS
Similar notebooks: