for(vector<const RecTrack*>::iterator irec = recTracks.begin();
irec != recTracks.end(); irec++){
std::vector<RecHit> hits = (*irec)->recHits();
int nhit = hits.size();
for (int i=0; i<SimpleTree::NHITMAX && i<nhit; i++) {
status[i] = hits[i].isMatched() ? 0 : 1 ;
}
}