static void
redraw( /* mark portion of a tree for redraw */
register RTREE *tp,
int x0,
int y0,
int x1,
int y1,
int l[2][2]
)
{
int quads = CH_ANY;
int mx, my;
register int i;
/* compute midpoint */
mx = (x0 + x1) >> 1;
my = (y0 + y1) >> 1;
/* see what to do */
if (l[0][0] > mx)
quads &= ~(CHF(UL)|CHF(DL));
else if (l[0][1] < mx)
quads &= ~(CHF(UR)|CHF(DR));
if (l[1][0] > my)
quads &= ~(CHF(DR)|CHF(DL));
else if (l[1][1] < my)
quads &= ~(CHF(UR)|CHF(UL));
tp->flgs |= quads; /* mark quadrants for update */
/* climb the branches */
for (i = 0; i < 4; i++)
if (tp->flgs & BRF(i) && quads & CHF(i))
redraw(tp->k[i].b, i&01 ? mx : x0, i&02 ? my : y0,
i&01 ? x1 : mx, i&02 ? y1 : my, l);
}
bool Jet::isGood() const {
// bool passesPt = pt() > 30; // Bristol original value
bool passesPt = pt() > 35; // 19.07.11 Chris's value
bool passesEta = fabs(eta()) < 2.4;
bool jetID = false;
//if (usedAlgorithm == JetAlgorithm::ParticleFlow || usedAlgorithm == JetAlgorithm::PF2PAT) {
if (usedAlgorithm == JetAlgorithm::CA08PF || usedAlgorithm == JetAlgorithm::PF2PAT) {
bool passNOD = NOD() > 1;
bool passCEF = CEF() < 0.99;
bool passNHF = NHF() < 0.99;
bool passNEF = NEF() < 0.99;
bool passCHF = true;
bool passNCH = true;
if (fabs(eta()) < 2.4) {
passCHF = CHF() > 0;
passNCH = NCH() > 0;
}
jetID = passNOD && passCEF && passNHF && passNEF && passCHF && passNCH;
}
else{
bool passesEMF = emf() > 0.01;
bool passesN90Hits = n90Hits() > 1;
bool passesFHPD = fHPD() < 0.98;
jetID = passesEMF && passesN90Hits && passesFHPD;
}
return passesPt && passesEta && jetID;
}
static void
update( /* update tree display as needed */
uby8 ca[3], /* returned average color */
register RTREE *tp,
int x0,
int y0,
int x1,
int y1
)
{
int csm[3], nc;
register uby8 *cp;
uby8 rgb[3];
double dpth2[4], d2;
int gaps = 0;
int mx, my;
register int i;
/* compute leaf depths */
d2 = FHUGE*FHUGE;
for (i = 0; i < 4; i++)
if (tp->flgs & LFF(i)) {
FVECT dv;
register float *wp = qtL.wp[tp->k[i].li];
dv[0] = wp[0] - odev.v.vp[0];
dv[1] = wp[1] - odev.v.vp[1];
dv[2] = wp[2] - odev.v.vp[2];
dpth2[i] = DOT(dv,dv);
if (dpth2[i] < d2)
d2 = dpth2[i];
}
d2 *= (1.+qtDepthEps)*(1.+qtDepthEps);
/* compute midpoint */
mx = (x0 + x1) >> 1;
my = (y0 + y1) >> 1;
/* draw leaves */
csm[0] = csm[1] = csm[2] = nc = 0;
for (i = 0; i < 4; i++) {
if (tp->flgs & LFF(i) && dpth2[i] <= d2) {
cp = qtL.rgb[tp->k[i].li];
csm[0] += cp[0]; csm[1] += cp[1]; csm[2] += cp[2];
nc++;
if (tp->flgs & CHF(i))
dev_paintr(cp, i&01 ? mx : x0, i&02 ? my : y0,
i&01 ? x1 : mx, i&02 ? y1 : my);
} else if ((tp->flgs & CHBRF(i)) == CHF(i))
gaps |= 1<<i; /* empty stem */
}
/* do branches */
for (i = 0; i < 4; i++)
if ((tp->flgs & CHBRF(i)) == CHBRF(i)) {
update(rgb, tp->k[i].b, i&01 ? mx : x0, i&02 ? my : y0,
i&01 ? x1 : mx, i&02 ? y1 : my);
csm[0] += rgb[0]; csm[1] += rgb[1]; csm[2] += rgb[2];
nc++;
}
if (nc > 1) {
ca[0] = csm[0]/nc; ca[1] = csm[1]/nc; ca[2] = csm[2]/nc;
} else {
ca[0] = csm[0]; ca[1] = csm[1]; ca[2] = csm[2];
}
/* fill in gaps with average */
for (i = 0; gaps && i < 4; gaps >>= 1, i++)
if (gaps & 01)
dev_paintr(ca, i&01 ? mx : x0, i&02 ? my : y0,
i&01 ? x1 : mx, i&02 ? y1 : my);
tp->flgs &= ~CH_ANY; /* all done */
}
