**
** ROUTINE: taQuasarsTargetSelect
**
**
** Set the appropriate bit in the primary target bit mask for each
** Quasars Working Group target selection criterion this object
** satisfies. This routine is responsible for setting the following
** bits:
**
** It is assumed that the above bits are unset on entry, although other
** bits may be set.
**
** This routine returns the priority of the quasar candidate, in the
** range 0-255. On error, -1 is returned. Zero means it should not
** be targetted. We would like to give the same priority (200) to all
** QSO targets except FIRST QSO targets which are selected as being
** on a line between two FIRST sources. The latter type of QSO
** receives lower priority (100). In the future the lower priority
** will be replaced by a range of priority.
**
** NOTE: Update me after all changes are made.
**
**
**
**
******************************************************************************
*******************************************************************************/
int
taQuasarsTargetSelect(TA_CALIB_OBJ *obj,
/* MOD: calibrated object to test */
const TA_QUASARS_PARAMS *params
/* IN: tunable parameters*/)
{
float mag[5],err[5];
double syserr[5];
double lowzlimit[5];
/* double object_k,object_l,object_m; */
int colorindex;
/* double limitcorrection; */
int capobject, index, tmpflag,tmpflag2, maglimit, tsflag;
int inlocus, inGalaxylocus = 0;
int ImageProb=0, badCol=0;
int ta[5], det[5];
int priority=0;
float fseed, random, deltamag;
int closestpoint; /* added for debugging purposes - not used to target */
double ellipsedist; /* added for debugging purposes - not used to target */
int targetmidZ = 0;
int inAstarbox = 0;
int inWDbox = 0;
int inWMbox = 0;
int validtarget = 1;
/**** temp output ****/
FILE *ftmp;
int i, nchild;
int intercenter, intercenterI;
/*
*************************************************************************
** FATAL Errors
**
** These next sections look for objects with fatal errors. If the
** object is bad, we return 0, which means the object is never
** given a chance to be targeted.
**
*************************************************************************
**/
/**********************************************************************
** Cut "bad" objects
**
** Don't even bother to check objects which are bad. Inlcuding
** bright objects, edge objects and those that are flagged as
** saturated.
**
*********************************************************************/
if (obj->objc_flags&(AR_DFLAG_BRIGHT|AR_DFLAG_EDGE|AR_DFLAG_SATUR)) {
return 0;
}
/**********************************************************************
** Cut blended objects
**
** Do not select things which are blended but not deblended.
**
** NOTE: This will reject small separation gravitational lenses.
**
** NOTE: Who commented out the NODEBLEND thing? Please don't embed
** comments in the code in this fashion. Rather, comment out the
** entire old line, replace the line and give a comment explaining
** the change. Do we really want to reject ALL blended objects?
**
*********************************************************************/
if ((obj->objc_flags & AR_DFLAG_BLENDED) /* && !(obj->objc_flags & AR_DFLAG_NODEBLEND) */ ) {
return 0;
}
/**********************************************************************
** Reject objc_types other than stellar or galaxy
**
** Reject things which don't have either objc_type == 3 (galaxy) or
** objc_type == 6 (stellar), because the other types are all error
** codes.
**
*********************************************************************/
if ((obj->objc_type != 3) && (obj->objc_type != 6)) {
return 0;
}
/**********************************************************************
** Disregard moving objects
**
** Objects which are deblended as moving are probably asteriods,
** but certainly not quasars. However, their colors can be
** non-solar if the deblend isn't perfect, so we have to reject
** them explicitly.
**
*********************************************************************/
if((obj->objc_flags2 & AR_DFLAG2_DEBLENDED_AS_MOVING))
return 0;
/**********************************************************************
** Disregard more moving objects
**
** Some objects appear to be moving, but the deblending failed.
** Again we must reject such objects explicitly. Throw away things
** that have large velocities (if those velocities are significant
** with respect to the error).
**
*********************************************************************/
if( !(obj->objc_flags2 & AR_DFLAG2_BAD_MOVING_FIT) &&
( (fabs(obj->rowv.val) > maxV ) && (fabs(obj->rowv.val)/obj->rowv.err > maxVSig) ||
(fabs(obj->colv.val) > maxV ) && (fabs(obj->colv.val)/obj->colv.err > maxVSig)) )
return 0;
/**********************************************************************
** Correct magnitudes and errors for galactic reddening
**
** We use the Schlegel maps to correct for galactic reddening. We
** then must add the error in the reddening in quadrature with the
** error in the magnitude.
**
** NOTE: Schlegel says that the error in the reddening correction
** is 15% of the value, which is what is coded. It may be that
** since the errors in different filters are correlated that we now
** overestimate the errors in colors.
**
*********************************************************************/
ta[0] = TA_U;
ta[1] = TA_G;
ta[2] = TA_R;
ta[3] = TA_I;
ta[4] = TA_Z;
for (colorindex=0; colorindex<5; colorindex++) {
mag[colorindex] = obj->detection[ta[colorindex]].psfCounts.val-obj->reddening[ta[colorindex]];
err[colorindex] = sqrt(obj->detection[ta[colorindex]].psfCounts.err\
*obj->detection[ta[colorindex]].psfCounts.err\
+0.15*0.15*obj->reddening[ta[colorindex]]\
*obj->reddening[ta[colorindex]]);
}
/**********************************************************************
** Identify objects with image problems
**
** NOTE: Xiaohui or Michael, please comment (in further detail).
** Also, please explain why you set ImageProb=1 for some of these
** rather than simply returning 0.
**
*********************************************************************/
if((obj->objc_flags & AR_DFLAG_CHILD)!=0) {
nchild=1;
} else {
nchild=0;
}
intercenter=0;
intercenterI=0;
/* if object is a child and has PEAKCENTER or NOTCHECKED flags in high S/N band,
* there must be a problem of deblending */
for (colorindex=0; colorindex<5; colorindex++) {
tmpflag=obj->detection[ta[colorindex]].flags;
tmpflag2=obj->detection[ta[colorindex]].flags2;
if (tmpflag & AR_DFLAG_BINNED1) det[colorindex]=1;
else det[colorindex]=0;
if(mag[colorindex]<23 && err[colorindex]