Difference between revisions of "Working with L1688"

From CoolWiki
Jump to navigationJump to search
m
m
Line 25: Line 25:
 
[[How can I make a color composite image using Spitzer and/or other data? | Go here]] and follow the tutorial (or use your own software) to explore the Spitzer mosaics as well as the mosaics at other wavelengths.  Think about the questions at the end of that page as they pertain to these data.
 
[[How can I make a color composite image using Spitzer and/or other data? | Go here]] and follow the tutorial (or use your own software) to explore the Spitzer mosaics as well as the mosaics at other wavelengths.  Think about the questions at the end of that page as they pertain to these data.
  
Look for image artifacts in the mosaics - they are certainly still there.  What are real features of the sky and what is a feature of the instrument and/or reduction?   
+
Look for image artifacts in the mosaics - they are certainly still there.  What are real features of the sky and what is a feature of the instrument and/or reduction?  Which objects are saturated, in which bands?
  
 
=Previously identified sources=
 
=Previously identified sources=
Line 32: Line 32:
  
 
[[How can I find out what scientists already know about a particular astronomy topic or object? | Go here]] and follow the tutorial to explore the literature. ([[I'm ready to go on to the "Advanced" Literature Searching section | Move on to these examples and questions when you're ready]].) Find some old and some recent papers on L1688.  For how many decades have people been studying this region?  Find some of the famous objects in your images.  Note that this is a pretty famous region, so it's doubtful that you will easily be able to construct a list of ALL the famous objects in this region in a short time, even using SIMBAD.  Keep in mind too that people work in different coordinate systems.  A paper reporting on radio data may have a slightly different coordinate system than Spitzer (which is tied to 2MASS).  You will have to use your judgment to decide if any given object is really the same as the object in the Spitzer image, or if the Spitzer image didn't really detect the object.
 
[[How can I find out what scientists already know about a particular astronomy topic or object? | Go here]] and follow the tutorial to explore the literature. ([[I'm ready to go on to the "Advanced" Literature Searching section | Move on to these examples and questions when you're ready]].) Find some old and some recent papers on L1688.  For how many decades have people been studying this region?  Find some of the famous objects in your images.  Note that this is a pretty famous region, so it's doubtful that you will easily be able to construct a list of ALL the famous objects in this region in a short time, even using SIMBAD.  Keep in mind too that people work in different coordinate systems.  A paper reporting on radio data may have a slightly different coordinate system than Spitzer (which is tied to 2MASS).  You will have to use your judgment to decide if any given object is really the same as the object in the Spitzer image, or if the Spitzer image didn't really detect the object.
 +
 +
=Doing photometry=
 +
 +
Using MOPEX to do many sources at once: [[I'm ready to advance to a highly technical and in-depth discussion on downloading Spitzer data. | Go here]] and get the doc file that is linked in partway down. Find the part on using Apex-1Frame.  Follow that for each channel to get the photometry.
 +
 +
Using APT to do one source at a time: COMING SOON.  If you do this, I recommend you start with a short list of famous objects from the above step.
 +
 +
Reduced (and bandmerged) photometry is included in the tarball above if you wish to short-cut this process, and I promise I won't think any less of you.
 +
 +
=Bandmerging the photometry=
 +
 +
Now that you have a measurement of the brightness of each object in each frame, you need to match the same object across multiple bands, called "bandmerging."  Maybe, if you used APT, you kept track of that manually. If you did the photometry automatically, you need to have the computer match the sources.  Some of the sources you extracted might not even be real.
 +
 +
There is no GUI tool right now to bandmerge the photometry; there might be one eventually.
 +
 +
The code I use is written in IDL (expensive to have a site license, even for academics), and for each band, starting with 2MASS JHK, looks for matches within 1 arcsecond in IRAC-1, then IRAC-2, then IRAC-3, then IRAC-4, and then within 2 arcseconds in MIPS-1, and then within 10 arcseconds in MIPS-2.  (Bonus question: why should this matching distance change?)
 +
 +
Bandmerged photometry (including 2MASS data!) is included in the tarball above if you wish to short-cut this process, and I promise I won't think any less of you.
 +
  
 
=Working with the data tables=
 
=Working with the data tables=
  
 
I've provided the data tables in something called "IPAC table format" in the tarball above.  There is a long version with many columns, and a short version with many fewer columns.  Spot, Leopard, and Skyview all understand IPAC table format.  You can also import the file into Excel to play with it there.  NOTE THAT there are a LOT of sources in this relatively small mosaic.  You will want the computer to do as much work for you as possible - you don't want to be counting sources manually.
 
I've provided the data tables in something called "IPAC table format" in the tarball above.  There is a long version with many columns, and a short version with many fewer columns.  Spot, Leopard, and Skyview all understand IPAC table format.  You can also import the file into Excel to play with it there.  NOTE THAT there are a LOT of sources in this relatively small mosaic.  You will want the computer to do as much work for you as possible - you don't want to be counting sources manually.
 +
 +
You can use Spot, Leopard, or Skyview to overplot the sources onto an image of your choice.  NOTE AGAIN THAT there are a LOT of sources in this relatively small mosaic.
  
 
Import the data table into Excel, as that's probably the easiest way for you to work with it.  (DO I NEED TO WRITE A TUTORIAL ON THIS?)
 
Import the data table into Excel, as that's probably the easiest way for you to work with it.  (DO I NEED TO WRITE A TUTORIAL ON THIS?)
 +
 +
=Making color-color and color-magnitude plots=

Revision as of 23:31, 13 May 2008

There are a LOT of things on this Wiki, very few of which are arranged linearly. The Wiki is a web, after all, so there are links interwoven all over the place. However, we do realize that this can be frustrating for new users, so here is a single analysis chain, start to finish, with the products already done so that you can step into and out of the process as you want or need.

Download the tarball with all my reduced images and data tables here.

Downloading the data

Go here and follow the tutorial to investigate what observations exist and download the data. Search on L1688. Find the AORs, download the IRAC and MIPS mosaics. The data I'm using for this example are all from pid 177, but you should investigate the coverage of the other programs too.

It turns out, though, that L1688 is part of the c2d Legacy program (pid 177 is part of that), so you may also wish to consult the page on already-reduced Spitzer data.

Making the mosaics

Go here and follow the tutorial to make the mosaics from the post-BCD pieces you have downloaded.

Mosaics from all 4 IRAC channels and 2 of the 3 MIPS channels are included in the tarball above if you wish to short-cut this process.

Getting data from other wavelengths

Go here and follow the tutorial to find data (images and data tables) at other wavelengths. Think about the questions at the end of that page as they pertain to these data.

A POSS image is included in the tarball above if you wish to short-cut this process. Reduced photometry from 2MASS is included in the data tables as well.

Investigating the mosaics

Go here and follow the tutorial (or use your own software) to explore the Spitzer mosaics as well as the mosaics at other wavelengths. Think about the questions at the end of that page as they pertain to these data.

Look for image artifacts in the mosaics - they are certainly still there. What are real features of the sky and what is a feature of the instrument and/or reduction? Which objects are saturated, in which bands?

Previously identified sources

How many of these sources are 'famous' and how many are new that you have just discovered?

Go here and follow the tutorial to explore the literature. ( Move on to these examples and questions when you're ready.) Find some old and some recent papers on L1688. For how many decades have people been studying this region? Find some of the famous objects in your images. Note that this is a pretty famous region, so it's doubtful that you will easily be able to construct a list of ALL the famous objects in this region in a short time, even using SIMBAD. Keep in mind too that people work in different coordinate systems. A paper reporting on radio data may have a slightly different coordinate system than Spitzer (which is tied to 2MASS). You will have to use your judgment to decide if any given object is really the same as the object in the Spitzer image, or if the Spitzer image didn't really detect the object.

Doing photometry

Using MOPEX to do many sources at once: Go here and get the doc file that is linked in partway down. Find the part on using Apex-1Frame. Follow that for each channel to get the photometry.

Using APT to do one source at a time: COMING SOON. If you do this, I recommend you start with a short list of famous objects from the above step.

Reduced (and bandmerged) photometry is included in the tarball above if you wish to short-cut this process, and I promise I won't think any less of you.

Bandmerging the photometry

Now that you have a measurement of the brightness of each object in each frame, you need to match the same object across multiple bands, called "bandmerging." Maybe, if you used APT, you kept track of that manually. If you did the photometry automatically, you need to have the computer match the sources. Some of the sources you extracted might not even be real.

There is no GUI tool right now to bandmerge the photometry; there might be one eventually.

The code I use is written in IDL (expensive to have a site license, even for academics), and for each band, starting with 2MASS JHK, looks for matches within 1 arcsecond in IRAC-1, then IRAC-2, then IRAC-3, then IRAC-4, and then within 2 arcseconds in MIPS-1, and then within 10 arcseconds in MIPS-2. (Bonus question: why should this matching distance change?)

Bandmerged photometry (including 2MASS data!) is included in the tarball above if you wish to short-cut this process, and I promise I won't think any less of you.


Working with the data tables

I've provided the data tables in something called "IPAC table format" in the tarball above. There is a long version with many columns, and a short version with many fewer columns. Spot, Leopard, and Skyview all understand IPAC table format. You can also import the file into Excel to play with it there. NOTE THAT there are a LOT of sources in this relatively small mosaic. You will want the computer to do as much work for you as possible - you don't want to be counting sources manually.

You can use Spot, Leopard, or Skyview to overplot the sources onto an image of your choice. NOTE AGAIN THAT there are a LOT of sources in this relatively small mosaic.

Import the data table into Excel, as that's probably the easiest way for you to work with it. (DO I NEED TO WRITE A TUTORIAL ON THIS?)

Making color-color and color-magnitude plots