**
** 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.
**
**
**
******************************************************************************
*******************************************************************************/
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, inlocus, index, tmpflag,tmpflag2;
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 */
/**** temp output ****/
FILE *ftmp;
int i,nchild;
int intercenter, intercenterI;
/**** ---- ****/
/* First, look at the flags to see if this is an object we might
like to select */
if (obj->objc_flags&(AR_DFLAG_BRIGHT|AR_DFLAG_EDGE|AR_DFLAG_SATUR)) {
return 0;
}
/* Do not select things which are blended but not deblended */
if ((obj->objc_flags & AR_DFLAG_BLENDED) /* && !(obj->objc_flags & AR_DFLAG_NODEBLEND) */ ) return 0;
/* Get magnitudes with reddening applied */
ta[0] = TA_U;
ta[1] = TA_G;
ta[2] = TA_R;
ta[3] = TA_I;
ta[4] = TA_Z;
/* 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. */
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]]);
}
/* disregard moving objects */
if((obj->objc_flags2 & AR_DFLAG2_DEBLENDED_AS_MOVING) )
return 0;
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;
if((obj->objc_flags & AR_DFLAG_CHILD)!=0) nchild=1;
else nchild=0;
intercenter=0;
intercenterI=0;
/***** identify objects with image problems *****/
/* 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]