23.03.2013 Views

IRAF Exercise II - Quick-look Photometry with Imexamine

IRAF Exercise II - Quick-look Photometry with Imexamine

IRAF Exercise II - Quick-look Photometry with Imexamine

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Name:________________________________________ Due Date:____________ Score: _____ / 30<br />

<strong>IRAF</strong> <strong>Exercise</strong> <strong>II</strong> - <strong>Quick</strong>-<strong>look</strong> <strong>Photometry</strong> <strong>with</strong> <strong>Imexamine</strong><br />

There is an image called "spicamfocus.fits.gz" in the directory:<br />

/astro/classes/Astro_480/exercises/exercise2/<br />

that you should copy to your working directory. It is a typical focus frame, taken <strong>with</strong> an instrument<br />

called SPIcam on the Apache Point 3.5 meter telescope. Each pixel spans 0.28 arcseconds on the sky.<br />

Thought Question: You notice that one of the brightest stars has a diameter on the chip of about 10<br />

pixels. How many arc seconds is that?<br />

The picture contains multiple images of each star each at a different focus setting. This is accomplished<br />

by doing 5 iterations of 1) open the shutter for 10 sec, 2) close the shutter, 3) shift the charge on the<br />

detector down by a few rows, 4) change the telescope focus, 5) repeat. The CCD is then completely<br />

read out. To keep track of things, there is a double-skip in the charge transfer on the last iteration. For<br />

this frame, assume the focus was set to (1100, 1200, 1300, 1400, 1500) microns <strong>with</strong> the 1500 micron<br />

setting corresponding to the images that are "double-spaced."<br />

Suggestion: Use tools <strong>with</strong>in ds9 viewer to rotate the image so that your view <strong>look</strong>s like this one. Also,<br />

there are overscan regions to the left, bottom, and right of the image that will play a role later in this<br />

exercise. Be sure to adjust parameters in <strong>IRAF</strong> display task to view image logarithmically or open the<br />

file directly from ds9. You will want to have the contrast such that you see the entire “stack” and the<br />

overscan regions.<br />

Fig. 1 It’s hard to replicate via sketches how the<br />

image on the left was created, and the actual method<br />

involved reading out just part of the CCD chip and<br />

taking another image. However, imagine that you<br />

have just 2 stars in the field. You take an image. You<br />

read out a number of rows of the CCD, and take<br />

another image. The images of any given star will be in<br />

a straight line vertically. It’s like stacking transparent<br />

sheets of paper that have the same pictures on them,<br />

but the vertical alignment is offset a bit for each sheet.<br />

Caution: Be sure to keep the focus settings straight between the image and the tables below. Also, the<br />

stars you choose cannot be saturated in any of the focus frames, so check all 5 images of the star<br />

before actually choosing it. In general, choose the faintest stars where you can still identify all 5<br />

images of them.<br />

4/17/12 iraf2_exercise 1


1. By how many rows is the image single-space shifted between focus settings? ______<br />

2. When the telescope has its best focus, the stellar images should be axisymmetric (ellipticity<br />

= 0), smallest (min FWHM), and have the maximum peak intensity. Check out the 5 images of the<br />

faint star whose 3rd image is at about (693, 470). Notice that the 3rd image seems better than the rest.<br />

This is a clue that the best focus setting will likely be around 1300 microns. We will now confirm this<br />

quantitatively. Edit the "rimexam" parameter file so that the object radius=4, background<br />

buffer width=2 and background width=1. Now choose two well-behaved (unsaturated and<br />

isolated) stars and fill in the tables below.<br />

Star 1: Double-skip centroid: x = ________ y = ________ (pixels)<br />

Focus Setting Peak Flux Ellipticity<br />

1100<br />

1200<br />

1300<br />

1400<br />

1500<br />

Focus Setting<br />

1100<br />

1200<br />

1300<br />

1400<br />

1500<br />

FWHM<br />

(pixels)<br />

Star 2: Double-skip centroid: x = ________ y = ________ (pixels)<br />

Peak<br />

Flux<br />

Ellipticity<br />

FWHM<br />

(pixels)<br />

FWHM (arcsec)<br />

FWHM (arcsec)<br />

3. What would you choose as the best focus setting for the telescope? On what basis? Which<br />

FWHM value did you use?<br />

2<br />

In your answering of the next question, we need to make sure that you<br />

understand completely how the spicamfocus image was made.<br />

•What information do the overscan regions contain? _________________<br />

•How many “layers of sky” are in the brightest part of the image? ______<br />

•In which part of the image will you get values of the real sky? ________<br />

•Since that part will also contain read noise, how will you get just the<br />

values of the sky in ADUs? ____________________________________<br />

4. Using just the ds9 pixel readout (accessed from the Analysis<br />

menu at the top), what would you estimate as the background sky value for<br />

this image, in A/D units? ______________


5. Investigate how the rimexam parameters affect the results you obtain for aperture<br />

photometry. The settings of the object radius, buffer width, and background width will affect the<br />

determination of the sky values and object flux. Edit the rimexam parameter file and set rplot=15<br />

(NOTE: this is not the aperture radius), buffer=1, width=1. For the values of object radius given, fill in<br />

the table below (see next page). As you fill in the table, notice how the sky value in the radial plot<br />

varies. Use the isolated star at (523,497) that is in focus. [Don’t forget to set up a log file to save your<br />

work. Also, make sure that the iterations parameter in rimexam is set to 1!]<br />

(radius = 1.) Object radius #change this value as you work through this step<br />

For Object at x = ___523___ y = ___497___<br />

Object radius Magnitude Flux Sky Peak FWHM<br />

1<br />

2<br />

3<br />

4<br />

8<br />

12<br />

16<br />

6. From your data in the table above, explain how the sky, peak and FWHM values depend on<br />

the object radius. There were probably unexpected results or trends. If so, list one and come up <strong>with</strong><br />

your guess as to why it occurred.<br />

7. <br />

a) Graph the FWHM and magnitude (smallest numbers at the top of the axis for both) versus<br />

aperture radius values in the graph shown on the last page (use different colors or symbols and each<br />

axis will have different values and units).<br />

b) Do they tend to asymptotic values as the object radius gets larger? Explain.<br />

8. If the "right" answer is obtained <strong>with</strong> a very large object radius, what disadvantages might<br />

dissuade you from using an object radius of, say, 25 pixels?<br />

3


9. How does the task imexamine determine the astronomical flux from a star? We are not<br />

<strong>look</strong>ing for the steps you took (such as, epar’ing, typing ‘r,’ making a log file, etc.), but the actual<br />

process going on <strong>with</strong>in the programs. What is <strong>IRAF</strong> doing to get the peaks, FWHM, sky values, and<br />

magnitudes? This means digging into phelp imexamine.<br />

10. Learning objectives are given at the top of the tutorial page for this exercise. In narrative<br />

form, using good writing style, meet the first 5 objectives. For the proficiency goal: What do you think<br />

you will easily remember as we go on to the next <strong>IRAF</strong> exercises, and why?<br />

4


Name ____________________________________<br />

5


We have intentionally put only this sentence on this page.<br />

6

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!