PNe8alpha3 - adding contamination to our theoretical data

To better match the observations we have begun adding contamination to our theoretically generated data. This is the process:

Initial large catalogs for each field were generated by the preceeding methods (see PNe8alpha page). The new part now is the subsampling/contamination adding.
Add code: Virgo/FeldPNe/PNe8alpha3/subs.sm
The number of desired PNe for each field are set. For each field the correct number of PNe are sampled from the large catalog.
The contamination fraction is considered, and the (integer) value of contaminants is computed. The program steps backward from the end of that field's PNe and replaces the correct number with randomly generated coordinates that are within this field's boundaries [this requires reading in fits info/etc for each field].

Contamination is varied from 0% to 100%, in steps of 10%, for each field. Files are generated, here are the results:

NOTE: There is only ONE contaminated catalog made for each percentage... this might need to be fixed in later realisations.


code:  PNe8alpha3.sm  -  PNe8alpha3  - Main operations code, nearly identical to previous theoretical generators
           subs.sm         -       subs           -     Subsample the data and add contaminants
                                                                       Start by reading in all data files from each field
                                                                       Define number of PNe to sample from each field (now it is set based on how many have been observed, to compare)
                                                                       For each field (example: field 3) sample the proper number out of full catalog
                                                                       Read in fits header info
                                                                       Generate fake PNe (x,y) coordinates / convert to RA/DEC
                                                                       Back-substitute these into the sampled list from before
                                                                       repeat for all fields, 4, 7, Sub, etc
                                                                       write out final data file

GREEN circles are theoretical PNe
RED circles are added contaminants

Results: (mm means hand mask + sextractor mask ... it is without the surface brightness cut, basically)
                   Field:  3mm.fits    4mm.fits   7mm.fits   Submm.fits  Coremm.fits   LPCmm.fits   FCJmm.fits
Contamination          0%            0%            0%            0%                  0%                 0%                0%
        percentages:    10%          10%          10%          10%                10%               10%              10%
                                20%          20%          20%          20%                20%               20%              20%
                                30%          30%          30%          30%                30%               30%              30%
                                40%          40%          40%          40%                40%               40%              40%
                                50%          50%          50%          50%                50%               50%              50%
                                60%          60%          60%          60%                60%               60%              60%
                                70%          70%          70%          70%                70%               70%              70%
                                80%          80%          80%          80%                80%               80%              80%
                                90%          90%          90%          90%                90%               90%              90%
                              100%        100%        100%        100%               100%             100%            100%

3mm.fits
0% contaminated	10% contaminated
20% contaminated	30% contaminated
40% contaminated	50% contaminated
60% contaminated	70% contaminated
80% contaminated	90% contaminated
100% contaminated
4mm.fits
0% contaminated	10% contaminated
20% contaminated	30% contaminated
40% contaminated	50% contaminated
60% contaminated	70% contaminated
80% contaminated	90% contaminated
100% contaminated
7mm.fits
0% contaminated	10% contaminated
20% contaminated	30% contaminated
40% contaminated	50% contaminated
60% contaminated	70% contaminated
80% contaminated	90% contaminated
100% contaminated
Submm.fits
0% contaminated	10% contaminated
20% contaminated	30% contaminated
40% contaminated	50% contaminated
60% contaminated	70% contaminated
80% contaminated	90% contaminated
100% contaminated
Coremm.fits
0% contaminated	10% contaminated
20% contaminated	30% contaminated
40% contaminated	50% contaminated
60% contaminated	70% contaminated
80% contaminated	90% contaminated
100% contaminated
LPCmm.fits
0% contaminated	10% contaminated
20% contaminated	30% contaminated
40% contaminated	50% contaminated
60% contaminated	70% contaminated
80% contaminated	90% contaminated
100% contaminated
FCJmm.fits
0% contaminated	10% contaminated
20% contaminated	30% contaminated
40% contaminated	50% contaminated
60% contaminated	70% contaminated
80% contaminated	90% contaminated
100% contaminated

back to main Virgo page