# pppin script for object finding for RCS1 Spitzer sample IRAC data
set iecho 0
set inoise 0
set rwd 11
set conthres .3333
set smsig 3
set imerad 1
# iobjsm = 0: box smooth, 1: tapered box
set iobjsm 1


# read in channel 1
rdis 1
231953irac1
#read in channel 2
rdis 2
231953irac2

# create the zero pixel mask images in frames 3 and 4
cmas 1 3
cmas 2 4

# subtract off sky values using carm, 
#  sky is subtracted only in the unmasked part
carm 1 -61
3

carm 2 -74
4


#sum chanel 1 and 2 and store in 1
sum 1 2 1

#generate the net mask for the zero pixles, this requires the intersection
# of mask frames 3 and 4. This is accomplished by doing a trick with
# integer arithmatics

# first add the 2 masks, so that common areas in the mask now have values=2
#   pixels with masking only in one channel have values=1
sum 3 4 3

# this makes all the 2's into 1, and 1's into 0, now the mask is the
# intersection between the masks of chanels 1 and 2
car 3 *.4

# add a constant to the ch1+ch2 image,  using the masking info above
carm 1 +200
3

# pad the edges of the mask

padm 3
1
padm 3
1
padm 3
1 
  
 
# object finding for saturated stars

# set parameters for regular object finding
set thresky 1.70
set iir 9 
set ior 17
set adusat 600000

# object finding, must have two blank lines at end
afim 1
RCS231953
RCS231953, sum ch1 and 2
3
2
3


# for cosmic ray rejection, read in unsmoothed image again

#rejc 1


quit
# save final pos file
# first apply an additional merging of close objects (merp)
posf
merp 1


save 1 
XXXXXX_YY
XXXXXX_YY, wt:R+z, sm 3, cr rej, ths=2.6