Made it through the week!
Accomplished all of the plots
Took a huge bite out of yesterday's To-Do List. Got all of the multi-plots done in both fields, including the histograms! The only persisting problem was in the smooth plots, for some reason the contour lines only showed up on the first of the 9 plots, while all of the other labels and commands, etc. worked throughout the for-loop. This is a mystery, to be looked at later if need be.
To Do:
Read up on and implement Kolmogrov-Smirnov test
Update Research Document
...COSMOS!?
Friday, July 17, 2009
Thursday, July 16, 2009
Keep to the Code
Today: Spent the majority of the day engulfed in editing and writing code.
Fixed Random Distribution
Beth pointed out to me how I had inadvertently cut off the edges of the field, so this is now fixed. For some reason I kept coming up with errors when I tried to run the code for the southern field.
Made Multi-plots!
Since the northern field worked, I proceeded with that, adding a for-loop and a multi-plot command to run numerous trials and got a 3x3 set of RA-Dec plots, and a data structure saving these trials' data.
Smooth and Significance Histograms
Next on the list was to make the smoothed RA-Dec plots, and the accompanying statistical significance histograms. When I tried to run my code on the random data though, I again ran into errors, and have isolated a problem or two, but not the solutions. I think the next thing I'll try will be to skip over the part of the code that's tripping up and see if I can get the rest of it to work or if there are further issues.
To Do:
Southern random positions, and multi-plot, data structure
Smooth RA-Dec plots for both fields, and significance histograms
-finish debugging and get smoothed RA-Dec plots
-make accompanying cumulative histograms
Kolmogrov Smirnoff Test
-to compare trends of random trials and observed data in addition to eye-balling it
Update Research Summary Document
Ditto to Alex's Nap Idea
-late night showing of Harry Potter at the iMax: movie and company were great beyond expectations, but caused slight shortage of sleep.
Fixed Random Distribution
Beth pointed out to me how I had inadvertently cut off the edges of the field, so this is now fixed. For some reason I kept coming up with errors when I tried to run the code for the southern field.
Made Multi-plots!
Since the northern field worked, I proceeded with that, adding a for-loop and a multi-plot command to run numerous trials and got a 3x3 set of RA-Dec plots, and a data structure saving these trials' data.
Smooth and Significance Histograms
Next on the list was to make the smoothed RA-Dec plots, and the accompanying statistical significance histograms. When I tried to run my code on the random data though, I again ran into errors, and have isolated a problem or two, but not the solutions. I think the next thing I'll try will be to skip over the part of the code that's tripping up and see if I can get the rest of it to work or if there are further issues.
To Do:
Southern random positions, and multi-plot, data structure
Smooth RA-Dec plots for both fields, and significance histograms
-finish debugging and get smoothed RA-Dec plots
-make accompanying cumulative histograms
Kolmogrov Smirnoff Test
-to compare trends of random trials and observed data in addition to eye-balling it
Update Research Summary Document
Ditto to Alex's Nap Idea
-late night showing of Harry Potter at the iMax: movie and company were great beyond expectations, but caused slight shortage of sleep.
Wednesday, July 15, 2009
Random-ness
Today:
Axes fixed
Contour RA-Dec plots are now beautiful.
Random position generator
Given the code to get me going. Set up a randomu and manipulated the arrays to have values within the RA's and Dec's of the GOODS-N field. Spherematched these with actual data positions to accommodate the odd shape of the footprint. Plotted 500 of the random sample of stars (and oplotted the GOODS survey), so check that they overlapped properly. Oddly, the Dec's range was fine, but the outer endges of the RA's got chopped off. Wasn't obvious to me how this happened by looking at the code, so this will be addressed tomorrow.
Tomorrow:
Random position cont...
-check out why edges got cut off
-make for southern field
-emperical comparison to observed stars' positions
Axes fixed
Contour RA-Dec plots are now beautiful.
Random position generator
Given the code to get me going. Set up a randomu and manipulated the arrays to have values within the RA's and Dec's of the GOODS-N field. Spherematched these with actual data positions to accommodate the odd shape of the footprint. Plotted 500 of the random sample of stars (and oplotted the GOODS survey), so check that they overlapped properly. Oddly, the Dec's range was fine, but the outer endges of the RA's got chopped off. Wasn't obvious to me how this happened by looking at the code, so this will be addressed tomorrow.
Tomorrow:
Random position cont...
-check out why edges got cut off
-make for southern field
-emperical comparison to observed stars' positions
Tuesday, July 14, 2009
Contoured Density Plots
Started off the day going over isochrone images with Beth
The Besancon Model still seems puzzling, showing a lot of faint red stars. The Trilegal looks better, and the isochrone fitted at 20kpc nearly follows the line of stars (great by my standards, since in my limited experience few things have lined up nearly this well, but Beth mentioned some discrepancies). We determined that the disk/halo stars cut-off should be at a v-i of 1.5.
I haven't yet started working the separated stars yet, as I got distracted today...
Distribution, smoothing, significance
Brainstorming ensues. Beth explained how the smoothing worked a little bit, which brought up some minor pixel issues. To get out the edges and contrasts we wanted, I used SEARCH2D to make a pixel map of just 1's and 0's to separate the outer area with no data. The significance ((image-mean)/stddev) was taken of the data, and plotted with the "nodata" area set to 0.0 on the grey-scale. I also made an accompanying histogram of n-sigma; it wasn't Gaussian but it had the right tailing-off-shape.
Went to talk on BLAST-thought it was cool
Contour lines
Adding to the distribution density plot are now contour lines. It took me pretty much all afternoon figuring out and working on the code from Beth to make these show up properly. In the end, with a little typo-spotting help from Gail, I got it to work. I even patted myself on the back a little for going back making the contour lines shift from black to white so that you could still distinguish them on the lighter and darker areas. Only problems are the axes... even more so when I tried to make another plot of the northern field (as I'd been working primarily with the southern one all day, and didn't want the poor guy to feel neglected).
Tomorrow:
(in no particular order)
Fixing axes
-Go through code again and see if I can't find the problem
-If not, talk to Beth about it
Implement color-cut
Random distribution generator
The Besancon Model still seems puzzling, showing a lot of faint red stars. The Trilegal looks better, and the isochrone fitted at 20kpc nearly follows the line of stars (great by my standards, since in my limited experience few things have lined up nearly this well, but Beth mentioned some discrepancies). We determined that the disk/halo stars cut-off should be at a v-i of 1.5.
I haven't yet started working the separated stars yet, as I got distracted today...
Distribution, smoothing, significance
Brainstorming ensues. Beth explained how the smoothing worked a little bit, which brought up some minor pixel issues. To get out the edges and contrasts we wanted, I used SEARCH2D to make a pixel map of just 1's and 0's to separate the outer area with no data. The significance ((image-mean)/stddev) was taken of the data, and plotted with the "nodata" area set to 0.0 on the grey-scale. I also made an accompanying histogram of n-sigma; it wasn't Gaussian but it had the right tailing-off-shape.
Went to talk on BLAST-thought it was cool
Contour lines
Adding to the distribution density plot are now contour lines. It took me pretty much all afternoon figuring out and working on the code from Beth to make these show up properly. In the end, with a little typo-spotting help from Gail, I got it to work. I even patted myself on the back a little for going back making the contour lines shift from black to white so that you could still distinguish them on the lighter and darker areas. Only problems are the axes... even more so when I tried to make another plot of the northern field (as I'd been working primarily with the southern one all day, and didn't want the poor guy to feel neglected).
Tomorrow:
(in no particular order)
Fixing axes
-Go through code again and see if I can't find the problem
-If not, talk to Beth about it
Implement color-cut
Random distribution generator
Monday, July 13, 2009
Isochrones
Finishing up from last week, I made gray-scale RA-Dec density plots, hess-style. Also considering going back to cutting out some of the faintest blue objects.
With some new direction on how to proceed in looking at the distributions of the GOODS stars, I set out with it in mind to select out the disk vs halo stars into two groups. Using our brains for this, I started with an isochrone generated online to check where on the Trilegal model it would fall. I had to convert it from Absolute to Apparent magnitude, using 20kpc as the distance. I was pleased to find that the curve fell along the stars in my hess diagram.
Next is to use magnitude and color limits to pick out the groups of stars and check out their distributions for any kind of structure or randomness.
With some new direction on how to proceed in looking at the distributions of the GOODS stars, I set out with it in mind to select out the disk vs halo stars into two groups. Using our brains for this, I started with an isochrone generated online to check where on the Trilegal model it would fall. I had to convert it from Absolute to Apparent magnitude, using 20kpc as the distance. I was pleased to find that the curve fell along the stars in my hess diagram.
Next is to use magnitude and color limits to pick out the groups of stars and check out their distributions for any kind of structure or randomness.
Friday, July 10, 2009
Enough of this stuff, It's Friday, I'm in love
Yep, it's Friday. And I like that song. ^_^
So, with my nice new star catalogs, as of this morning I have:
1) made RA-Dec plots
2) made CMDs
3) made a histogram of the distribution of apparent magnitudes
4) done the above for the Besancon Model catalog for comparison
The N and S catalogs have 1151 and 1175 stars in them, respectively. In the RA-Dec plots their distribution appears pretty homogeneous across the fields. About 35% of the stars are magnitude 26.5 or brighter, and these have overall the same kind of random distribution, but with some odd voids in areas. (In the future- may consider a random coordinates generator to see the probabilities of distributions, to see if these voids are significant.)
The magnitude distribution of the model is of a similar shape, but with more evenly increasing numbers the fainter the stars, rather than the leap in numbers from brighter to fainter in the GOODS stars.
This afternoon I did a lot:
I started using the Trilegal model as another comparison for a star catalog.
I also learned how to make a Hess diagram, and with a little help from code from Alex and some of Dylan's work, I got it to work by the end of the day.
Further project direction Monday.
So, with my nice new star catalogs, as of this morning I have:
1) made RA-Dec plots
2) made CMDs
3) made a histogram of the distribution of apparent magnitudes
4) done the above for the Besancon Model catalog for comparison
The N and S catalogs have 1151 and 1175 stars in them, respectively. In the RA-Dec plots their distribution appears pretty homogeneous across the fields. About 35% of the stars are magnitude 26.5 or brighter, and these have overall the same kind of random distribution, but with some odd voids in areas. (In the future- may consider a random coordinates generator to see the probabilities of distributions, to see if these voids are significant.)
The magnitude distribution of the model is of a similar shape, but with more evenly increasing numbers the fainter the stars, rather than the leap in numbers from brighter to fainter in the GOODS stars.
This afternoon I did a lot:
I started using the Trilegal model as another comparison for a star catalog.
I also learned how to make a Hess diagram, and with a little help from code from Alex and some of Dylan's work, I got it to work by the end of the day.
Further project direction Monday.
Wednesday, July 8, 2009
I'm Baaaaack!
I have returned from the British Isles and returned to work.
Since I had put together my research summary paper before leaving, it was easy to review and jump back in right where I left off. I spent the morning cleaning up a few things in my paper, and then spent most of the day working on polishing up and finalizing my flux-ratio cut-off. By group meeting after lunch I'd made a few plots trying to be more precise with my parameters and sort through the data, systematically eliminating fainter objects that clouded over the color-color plot and getting down to the good stars. By the end of the day I made the final call: up to 1.4 (best at magnitude 26.5 or brighter, but am including up to 28).
That being done, I now have my stars catalog and am going to look at their distributions.
Since I had put together my research summary paper before leaving, it was easy to review and jump back in right where I left off. I spent the morning cleaning up a few things in my paper, and then spent most of the day working on polishing up and finalizing my flux-ratio cut-off. By group meeting after lunch I'd made a few plots trying to be more precise with my parameters and sort through the data, systematically eliminating fainter objects that clouded over the color-color plot and getting down to the good stars. By the end of the day I made the final call: up to 1.4 (best at magnitude 26.5 or brighter, but am including up to 28).
That being done, I now have my stars catalog and am going to look at their distributions.
Subscribe to:
Posts (Atom)
