probcut=-1.47 ;log10 of final probability ptot for cut (higher than this is a galc)
mag4cut=14.2
f6cut=400. ;mJy at 70um
bp=0.7 ;baseprob of being a galaxy in any high prob area

;I2-I4 vs M1 color-mag diagram
xp=mag2-mag4+0.8-(10.-mag5)/1.8 ;diagonal line with slope of 1.8, similar to line in figure
x0=where(xp le 0.) ; really bright objects, just stop at xp=0
if (x0(0) ne -1) then xp(x0)=0.
; probability from this diagram
pcc=bp*(1.-exp(-1.*xp)) ;this is the fuzzy cutoff below the line

;I2-I4 vs I2 color-mag diagram
n=n_elements(ra)
d=fltarr(n) ;the cutoff mag which depends on color
d(*)=13.05 ;cutoff is 13mag for objects redder than I2-I4=1.4
low=where(abs(mag2-mag4) lt 1.4 and t1 ne 2 and t2 ne 2) ;pointlike objects with .5 < I2-I4 < 1.4
lowe=where(abs(mag2-mag4) lt 1.4 and (t1 eq 2 or t2 eq 2)) ;extended objects with .5 < I2-I4 < 1.4
d(low)=14.5 ;cutoff for pointlike objs
d(lowe)=12.5 ;cutoff for extended objs
loww=where(abs(mag2-mag4) lt .5) ; blue objects, I2-I4 < .5
d(loww)=14.5
xp=1.2+.5*(mag2-d); is I2 above the cutoff set by "d" or not (with offset of 1.2)
x0=where(xp le 0.) ;anything below zero is just set to zero
xp(x0)=0.

pcm2_4=bp*(1.-exp(-1.*xp^3)) ;fuzzy cutoff
df=where(mag4 gt 10.5 and pcm2_4 lt .1) ;fainter objects that are above the line are set to have 10% probability of being a gal
pcm2_4(df) = 0.1

;I4-M1 vs M1 color-mag diagram
pcm245=fltarr(n)
r245=((mag4-mag5-3.5)/1.7)^2 + ((mag5-9.0)/2.)^2 ; circular region around the center of the bad box in the figure
pcm245(*)= exp(-1.* r245)
t=where(pcm245 lt .1) ;anything well outside the box still has a minimum 10% prob of being a gal
pcm245(t) = 0.1


; relatively high prob automatically for fainter objects near [8]-[24]=4 
ym=where((finite(mag5) and (mag4-mag5) gt 0.7 and (abs((mag4-mag5)-4.) lt 1.6) and mag5 gt 7.4) or (finite(mag5) ne 1))
pcm245(ym)=.7*bp ; i.e. =.49

ptot=pcc*pcm2_4*pcm245 ;product of probabilities

ym2=where((finite(mag5) and mag5 gt 10.) or mag4 gt mag4cut) ;if fainter than [24]=10 or [8.0]=14.2, then high prob
ptot(ym2)=.35

ext=where(t1 eq 2 or t2 eq 2); double prob for extended objs
ptot(ext)=2.*ptot(ext)

fnt5=where(mag5 gt 10. or finite(mag5) ne 1)
ptot(fnt5) = 0.1 ;this insures anything without a 24um flux or faint 24um will be a galc

f6hi=where(f6 gt f6cut and f3 gt 0. and f4 gt 0. and f5 gt 0.) ; mult prob by 0.1 for bright 70um sources that have I3,I4 fluxes
ptot(f6hi)=.1*ptot(f6hi)