Changes between Version 5 and Version 6 of BillsSAS39ColorColorPlots
- Timestamp:
- Aug 14, 2015, 1:44:12 PM (11 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
BillsSAS39ColorColorPlots
v5 v6 2 2 == Color Color diagrams from SAS 39 data == 3 3 4 Ken suggested that I look at the full force mean psf magnitude measurments by making color color diagrams. I thought it might be useful to compare these to the stack photometry so I did plots from those data as well. 4 Ken suggested that I look at the full force mean psf magnitude measurments by making color color diagrams. I thought it might be useful to compare these to the stack photometry so I did plots from those data as well. 5 5 6 In the following notes I describe what I did and show the resulting diagrams. 6 7 7 8 == Forced Photometry == 8 9 9 IPP psps forced photometry takes the positions of the detections made in the stack photometry and does forced photometry on the single frame warp images. 10 IPP psps forced photometry takes the positions of the detections made in the stack photometry and does forced photometry on the single frame warp images. These results are injested into DVO and 11 an averaging si done there. 10 12 11 I decided to use the "view" !ForcedMeanObjectView which is a join of the tables !ObjectThin and !ForcedMeanObject. I decided to cut on the qualityFlag requiring that the "good stack" bit be set. In order to keep they query time 12 managable I queried a subset of the survey using the projectionID column. This restricted the results to objects in projection cell 1406 (334 < ra < 338 0 < dec < 4) 13 14 The query was 13 I decided to use the "view" !ForcedMeanObjectView which is a join of the tables !ObjectThin and !ForcedMeanObject. 15 14 {{{ 16 15 select ForcedMeanObjectView.*, (ForcedMeanObjectView.qualityFlag & 1) as extended … … 18 17 where projectionID = 1406 and (ForcedMeanObjectView.qualityFlag) & 16 = 16 19 18 }}} 19 I decided to cut on the qualityFlag requiring that the "good stack" bit be set (see table ObjectQualityFlags for the bit detinitions) indicating that something (DVO ?) decided that we had good 20 stack measurements. In order to keep they query time and size of the results table 21 managable I queried a subset of the survey using the projectionID column. This restricted the results to objects in projection cell 1406 (334 < ra < 338 0 < dec < 4) 22 20 23 This yielded 562,311 objects. 21 24 22 The !ForcedMeanObjectView only has fluxes. Since this was my first set of queries in awhile I decided not toconvert the fluxes to magnitudes in the query. I simply grabbed the objects and ran the results through25 The !ForcedMeanObjectView only has fluxes. Since this was my first set of queries with the new schema, I decided to not convert the fluxes to magnitudes in the query. I simply grabbed the objects and ran the results through 23 26 various iterations of a [http://www.star.bristol.ac.uk/~mbt/stilts/ STILTS] program to convert to magnitudes.at I used this [wiki:forcedFluxToMagsExaample perl script] to facitiate the conversions. Adding 24 27 entries like … … 26 29 8.9-2.5*log10(gFPSFFlux) AS gFPSFMag 27 30 }}} 28 to the list of columns in the query would allow us to avoid the stilts step. However the database is going to be busy and keeping the math done there to a minimum seems like a reasonable thing to do.31 to the list of columns in the query would allow us to avoid the STILTS step. However the database is going to be busy and keeping the math done there to a minimum seems like a polite thing to do. 29 32 (8.9 is used as the zero point since fluxes are in Janskys). 30 33 34 35 === Star / Galaxy Separation === 36 37 I looked at the star galaxy separation discussion and consulted with experts and it was suggested that I use something simple. KronMag - PSFMag > 0 ==> star otherwise galaxy. Here is a [image:forced_dimag.jpg histogram] distribution from the 38 forced measurments. 39 40 41 31 42
