/*
 * dcresaccum.c		Jim Piper	02-02-88
 *
 * Dicentric finding system - accumulate results
 *
 * The input parameters are the "chromcands" structure for the current object,
 * the number of objects, and the monitoring file.
 *
 * MODIFICATIONS
 * -------------
 * 07-04-88	jimp	POSPRINT for information about classified positives.
 * 24-02-88	jimp	Do not permit dicentrics smaller than cell median area
 */
#include <stdio.h>
#include <dicen.h>
#include <otypes.h>

dcresaccum(chc,number,pf,r,POSPRINT)
struct chromcands *chc;
FILE *pf;
struct dcresults *r;
{
	int i;
	for (i=0; i<number; i++,chc++) {
		if (POSPRINT && (chc->btype & DICENTRIC)) {
			dcprint(chc,pf);
		}
		/* use the "UNSURE" bit to indicate small area */
		if ((chc->btype&DICENTRIC) && (chc->btype&UNSHURE)) {
			if (chc->Cbtype&DICENTRIC) {
				r->FN++;
				r->areaFN++;
			}
		}
		else if (chc->btype&DICENTRIC) {
			if (chc->Cbtype&DICENTRIC) {
				r->TP++;
				if (chc->Cbtype & UNSHURE) r->UTP++;
			}
			else if (chc->porigin & 022) {
				r->FP++;
				r->ovFP++;
			}
			else if (chc->Cbtype&CROSSARMS) {
				r->FP++;
				r->twistedFP++;
			}
			else if (chc->Cotype == 2 || chc->Cotype == 3) {
				r->FP++;
				r->compoFP++;
			}
			else
				r->FP++;
		}
		else {
			if (chc->Cbtype&DICENTRIC) {
				r->FN++;
				if (chc->Cbtype & ACROCENTRIC) r->acroFN++;
				else if (chc->porigin & 022) r->ovFN++;
				else if (chc->Cbtype&UNSHURE) r->UFN++;
				else if (chc->ncands < 2) r->onecandFN++;
			}
		}
	}
	fprintf(pf,"%2d %2d %2d [%1d %1d]",r->TP,r->FP,r->FN,r->UTP,r->UFN);
	fprintf(pf," (%1d %1d %1d) (%1d %1d %1d %1d)",r->compoFP,r->twistedFP,r->ovFP,r->acroFN,r->onecandFN,r->areaFN,r->ovFN);
	fprintf(pf," <%2d %2d>\n",r->nchrom,r->nwithcands);
}