IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changes between Version 5 and Version 6 of BillsSAS39ColorColorPlots


Ignore:
Timestamp:
Aug 14, 2015, 1:44:12 PM (11 years ago)
Author:
bills
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • BillsSAS39ColorColorPlots

    v5 v6  
    22== Color Color diagrams from SAS 39 data ==
    33
    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.
     4Ken 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. 
    55
     6In the following notes I describe what I did and show the resulting diagrams.
    67
    78== Forced Photometry ==
    89
    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.
     10IPP 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
     11an averaging si done there.
    1012
    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
     13I decided to use the "view" !ForcedMeanObjectView which is a join of the tables !ObjectThin and !ForcedMeanObject.
    1514{{{
    1615select ForcedMeanObjectView.*, (ForcedMeanObjectView.qualityFlag & 1) as extended   
     
    1817where projectionID = 1406 and (ForcedMeanObjectView.qualityFlag) & 16 = 16
    1918}}}
     19I 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
     20stack measurements. In order to keep they query time and size of the results table
     21managable 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
    2023This yielded 562,311 objects.
    2124
    22 The !ForcedMeanObjectView only has fluxes.  Since this was my first set of queries in awhile I decided not to convert the fluxes to magnitudes in the query. I simply grabbed the objects and ran the results through
     25The !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
    2326various 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
    2427entries like
     
    26298.9-2.5*log10(gFPSFFlux) AS gFPSFMag
    2730}}}
    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.
     31to 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.
    2932(8.9 is used as the zero point since fluxes are in Janskys).
    3033
     34
     35=== Star / Galaxy Separation  ===
     36
     37I 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
     38forced measurments.
     39
     40
     41
    3142