Exemple #1
0
def p60photsub(ofile, inlis, refimage, refstars, ot, filter):
   
    refstars.set_mag(filter)
    coofile=open('ot.coo','w')
    ot.wcs2pix(refimage)
    coofile.write('%10.3f%10.3f\n' % (ot[0].xval, ot[0].yval))
    coofile.close()

    outfile=open(ofile, 'a+')
  
    p60photutils.p60zeropt(inlis, refstars, filter)

    subimage=iraffiles('@%s' % inlis)
    for image in subimage:

        root,ext=image.split('.')
        iraf.phot('%s.sub' % root, coo='ot.coo', output='%s.grb' % root,
                  interac=no)

        stars=Starlist('%s.mag' % root)
        refstars.wcs2pix(image)
        zp,zpu=stars.zeropt(refstars,method='mean',rejout=0)

        otstars=Starlist('%s.grb' % root)    
        
        [mjd,utshut]=get_head(image,['OBSMJD','UTSHUT'])

        if (len(stars)>0):
            outfile.write('%s\t%s\t%.3f\t%.3f\t%.3f\t%.3f\n' % (mjd, utshut, otstars[0].mag+zp, sqrt(pow(zpu,2)+pow(otstars[0].magu,2)), otstars[0].magu, zpu))

    outfile.close()
Exemple #2
0
def p60sdssresults(inlis, outfile, object, filter, reflist, aperture=10.0):

    images = iraffiles(inlis)
    images.sort()

    outf = open(outfile, 'w')

    inf = open('%s.lis' % object, 'w')
    inf.write('%s-%s.fits' % (object, filter))
    inf.close()
    iraf.phot.apertures = aperture
    refstars = Starlist(reflist)
    p60photutils.p60zeropt('%s.lis' % object, refstars,
                           '%sMAG' % filter[0].upper())
    stars = Starlist('%s-%s.mag' % (object, filter))
    zp, zpu = stars.zeropt(refstars, method='mean', rejout=0)

    for image in images:

        root = image.split('.')[0]
        stars = Starlist('%s.mag' % root)
        [utshut, obsmjd, exptime] = get_head(image,
                                             ['UTSHUT', 'OBSMJD', 'EXPTIME'])
        # If well-detected, use mag/err from IRAF file
        if (len(stars) == 1) and (float(stars[0].magu) < 0.2):
            snmag = float(stars[0].mag) + zp
            snerr = sqrt(pow(float(stars[0].magu), 2) + pow(max(zpu, 0.03), 2))
        # For weak detection, provide upper limit
        elif (len(stars) == 1):
            magfile = open("%s.mag" % root)
            lines = magfile.readlines()
            sky = float(lines[len(lines) - 3].strip().split()[0])
            flux = float(lines[len(lines) - 1].strip().split()[3])
            npix = float(lines[len(lines) - 1].strip().split()[2])
            skystd = float(lines[len(lines) - 3].strip().split()[1])
            exp = float(lines[len(lines) - 2].strip().split()[0])
            if (flux > 0):
                flxerr = sqrt(flux + npix * pow(skystd, 2))
                snmag = -2.5 * log10((flux + 5.0 * flxerr) / exp) + 25.0 + zp
            else:
                flxerr = sqrt(npix) * skystd
                snmag = -2.5 * log10(5.0 * fluxerr / exp) + 25.0 + zp
            snerr = 99.00
        # For bad subtraction, just give -1.0
        else:
            snmag = -1.0
            snerr = -1.0

        outf.write('%s\t%15.7f%7.2f%10.2f%10.3f\n' %
                   (utshut, float(obsmjd), exptime, snmag, snerr))

    outf.close()
Exemple #3
0
def p60sdssresults(inlis, outfile, object, filter, reflist, aperture=10.0):

    images=iraffiles(inlis)
    images.sort()

    outf=open(outfile,'w')

    inf=open('%s.lis' % object, 'w')
    inf.write('%s-%s.fits' % (object, filter))
    inf.close()
    iraf.phot.apertures=aperture
    refstars=Starlist(reflist)
    p60photutils.p60zeropt('%s.lis' % object, refstars, 
                           '%sMAG' % filter[0].upper())
    stars=Starlist('%s-%s.mag' % (object,filter))
    zp,zpu = stars.zeropt(refstars,method='mean',rejout=0) 
    
    for image in images:

        root=image.split('.')[0]
        stars=Starlist('%s.mag' % root)
        [utshut,obsmjd,exptime] = get_head(image, ['UTSHUT','OBSMJD','EXPTIME'])
        # If well-detected, use mag/err from IRAF file
        if (len(stars)==1) and (float(stars[0].magu) < 0.2):
	    snmag = float(stars[0].mag)+zp
	    snerr = sqrt(pow(float(stars[0].magu),2)+pow(max(zpu,0.03),2))
        # For weak detection, provide upper limit
	elif (len(stars)==1):
            magfile = open("%s.mag" % root)
            lines=magfile.readlines()
            sky = float(lines[len(lines)-3].strip().split()[0])
            flux = float(lines[len(lines)-1].strip().split()[3])
            npix = float(lines[len(lines)-1].strip().split()[2])
            skystd = float(lines[len(lines)-3].strip().split()[1])
            exp = float(lines[len(lines)-2].strip().split()[0])
	    if (flux > 0):
	        flxerr = sqrt(flux + npix * pow(skystd,2))
	        snmag = -2.5 * log10((flux + 5.0 * flxerr) / exp) + 25.0 + zp
	    else:
	        flxerr = sqrt(npix) * skystd
                snmag = -2.5 * log10(5.0 * fluxerr / exp) + 25.0 + zp
            snerr = 99.00
        # For bad subtraction, just give -1.0
	else:
            snmag = -1.0; snerr = -1.0

        outf.write('%s\t%15.7f%7.2f%10.2f%10.3f\n' % 
                   (utshut, float(obsmjd), exptime, snmag, snerr))

    outf.close()
Exemple #4
0
def p60photsub(ofile, inlis, refimage, refstars, ot, filter):

    refstars.set_mag(filter)
    coofile = open('ot.coo', 'w')
    ot.wcs2pix(refimage)
    coofile.write('%10.3f%10.3f\n' % (ot[0].xval, ot[0].yval))
    coofile.close()

    outfile = open(ofile, 'a+')

    p60photutils.p60zeropt(inlis, refstars, filter)

    subimage = iraffiles('@%s' % inlis)
    for image in subimage:

        root, ext = image.split('.')
        iraf.phot('%s.sub' % root,
                  coo='ot.coo',
                  output='%s.grb' % root,
                  interac=no)

        stars = Starlist('%s.mag' % root)
        refstars.wcs2pix(image)
        zp, zpu = stars.zeropt(refstars, method='mean', rejout=0)

        otstars = Starlist('%s.grb' % root)

        [mjd, utshut] = get_head(image, ['OBSMJD', 'UTSHUT'])

        if (len(stars) > 0):
            outfile.write('%s\t%s\t%.3f\t%.3f\t%.3f\t%.3f\n' %
                          (mjd, utshut, otstars[0].mag + zp,
                           sqrt(pow(zpu, 2) + pow(otstars[0].magu, 2)),
                           otstars[0].magu, zpu))

    outfile.close()