TWindow.cxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// $Id: TWindow.cxx,v 1.7 2010/03/31 09:58:59 liucy Exp $
//-----------------------------------------------------------------------------
// Filename : TWindow.cc
// Section  : Tracking MDC
// Owner    : Yoshi Iwasaki
// Email    : [email protected]
//-----------------------------------------------------------------------------
// Description : A class to display tracking object.
//               See http://bsunsrv1.kek.jp/~yiwasaki/tracking/
//-----------------------------------------------------------------------------
 
#ifdef TRKRECO_WINDOW
 
#include <time.h>
//#include "belle.h"
//#include BELLETDF_H
//#include "tables/bestdf.h"
#include "CLHEP/String/Strings.h"
#include "TrkReco/TMDC.h"
#include "TrkReco/TWindow.h"
#include "TrkReco/TMLink.h"
#include "TrkReco/TTrack.h"
#include "TrkReco/TMLine.h"
#include "TrkReco/TSegment.h"
#include "TrkReco/TMDCWireHitMC.h"
#include "TrkReco/TTrackHEP.h"
#include "TrkReco/TConformalFinder0.h"
 
extern struct {
    int unpack_truncate;
    int unpackAll;
} calcdc_unpack_;
 
bool
TWindow::_skipAllWindow = false;
 
TWindow::TWindow(const TWindow & a)
: _name("copy of " + a._name),
  _window(600, 600, _name.c_str()),
  _wireName(a._wireName),
  _coordinate(a._coordinate),
  _xmin(a._xmin),
  _xmax(a._xmax),
  _ymin(a._ymin),
  _display(false),
  _objects(a._objects),
  _colors(a._colors),
  _targetRun(a._targetRun),
  _targetEvent(a._targetEvent),
  _mode(a._mode),
  _axial(a._axial),
  _stereo(a._stereo),
  _skip(false),
  _segmentLink(a._segmentLink),
  _segmentProperty(a._segmentProperty),
  _text(a._text),
  _canDelete(false),
  _nEvents(a._nEvents) {
    for (unsigned i = 0; i < 4; i++) _zoomHistory[i] = a._zoomHistory[i];
    _window.set_show_coordinates(_coordinate);
    _window.buttons_per_line(7);
    _closeButton = _window.button("next");
    _wireNameButton = _window.button("wire #");
    _axialButton = _window.button("axial");
    _stereoButton = _window.button("stereo");
    _coordinateButton = _window.button("coordinate");
    _zoomInButton = _window.button("zoom in");
    _zoomOutButton = _window.button("zoom out");
    _skipButton = _window.button("skip");
    _skipAllButton = _window.button("skip all");
    _szButton = _window.button("rphi/sz");
    _confButton = _window.button("conformal");
    _segmentLinkButton = _window.button("sgmt link");
    _segmentPropertyButton = _window.button("sgmt vec");
    _copyButton = _window.button("window cp");
    _window.init(_xmin, _xmax, _ymin);
}
 
TWindow::TWindow(const std::string & name)
: _name(name),
  _window(600, 600, name.c_str()),
  _wireName(false),
  _coordinate(false),
  _xmin(-100),
  _xmax(100),
  _ymin(-100),
  _display(false),
  _targetRun(0),
  _targetEvent(0),
  _mode(0),
  _axial(true),
  _stereo(false),
  _skip(false),
  _segmentLink(false),
  _segmentProperty(false),
  _canDelete(true),
  _nEvents(0) {
    _window.set_show_coordinates(_coordinate);
    _window.buttons_per_line(7);
    _closeButton = _window.button("next");
    _wireNameButton = _window.button("wire #");
    _axialButton = _window.button("axial");
    _stereoButton = _window.button("stereo");
    _coordinateButton = _window.button("coordinate");
    _zoomInButton = _window.button("zoom in");
    _zoomOutButton = _window.button("zoom out");
    _skipButton = _window.button("skip");
    _skipAllButton = _window.button("skip all");
    _szButton = _window.button("rphi/sz");
    _confButton = _window.button("conformal");
    _segmentLinkButton = _window.button("sgmt link");
    _segmentPropertyButton = _window.button("sgmt vec");
    _copyButton = _window.button("window cp");
    _window.init(_xmin, _xmax, _ymin);
    _zoomHistory[0].push(HepVector3D(_xmin, _xmax, _ymin));
    _zoomHistory[1].push(HepVector3D(0., 7., 0.));
    _zoomHistory[2].push(HepVector3D(_xmin, _xmax, _ymin));
    _zoomHistory[3].push(HepVector3D(-0.28, 0.28, -0.28));
}
 
TWindow::~TWindow() {
    clear();
}
 
void
TWindow::wait(void) {
    if (_skip) return;
    if (_skipAllWindow) return;
    if (! target()) return;
    draw();
 
    bool loop = true;
    bool zoom = false;
 
    while (loop) {
 
    //...Read input...
    double x0, y0;
    int b = _window.read_mouse(x0, y0);
 
    //...Close...
    if (b == _closeButton) loop = false;
    else if (b == _szButton) {
        if (_mode != 0) _mode = 0;
        else if (_mode != 2) _mode = 2;
        draw();
    }
    else if (b == _confButton) {
        if (_mode != 3) _mode = 3;
        else if (_mode != 1) _mode = 1;
        draw();
    }
    else if (b == _segmentLinkButton) {
        _segmentLink = ! _segmentLink;
        draw();
    }
    else if (b == _segmentPropertyButton) {
        _segmentProperty = ! _segmentProperty;
        draw();
    }
    else if (b == _skipButton) {
        loop = false;
        _skip = true;
    }
    else if (b == _skipAllButton) {
        loop = false;
        _skipAllWindow = true;
    }
    else if (b == _coordinateButton) {
        _coordinate = ! _coordinate;
        _window.set_show_coordinates(_coordinate);
        draw();
    }
    else if (b == _wireNameButton) {
        _wireName = ! _wireName;
        draw();
    }
    else if (b == _axialButton) {
        _axial = ! _axial;
        draw();
    }
    else if (b == _stereoButton) {
        _stereo = ! _stereo;
        draw();
    }
    else if (b == _copyButton) {
        TWindow tmp = * this;
        tmp.wait();
    }
    else if (b == _zoomInButton) {
        zoom = true;
    }
    else if (b == _zoomOutButton) {
        if (zoom) continue;
        if (_zoomHistory[_mode].size() > 1) {
        _zoomHistory[_mode].pop();
        }
        else {  
        _zoomHistory[_mode].pop();
        _xmin *= 2.;
        _xmax *= 2.;
        _ymin *= 2.;
        _zoomHistory[_mode].push(HepVector3D(_xmin, _xmax, _ymin));
        }
        zoom = false;
        draw();
    }
    else if (b == MOUSE_BUTTON(1)) {
        if (zoom) {
        double x, y;
        int bb = _window.read_mouse_rect(x0, y0, x, y);
        double xx(x0 < x ? x0 : x);
        double xxx(x0 < x ? x : x0);
        double yy(y0 < y ? y0 : y);
        _zoomHistory[_mode].push(HepVector3D(xx, xxx, yy));
        draw();
        zoom = false;
        }
    }
    }
}
 
void
TWindow::draw(void) {
    if (! target()) return;
    if (! _display) open();
 
    _window.clear();
    _nYLine = 0;
    HepVector3D z(_zoomHistory[_mode].top());
    _window.init(z.x(), z.y(), z.z());
 
    double xmin = _window.xmin();
    double xmax = _window.xmax();
    double ymin = _window.ymin();
    double ymax = _window.ymax();
    double x = xmin + (xmax - xmin) * .05;
    double y = ymin + (ymax - ymin) * .03;
    _xPositionText = x;
    _yPositionText = y;
    _yPositionStep = (ymax - ymin) * .03;
 
    drawCdc();
    drawAxes();
    drawHeader();
 
    unsigned n = _objects.length();
    for (unsigned i = 0; i < n; i++) {
    const TTrackBase & track = * _objects[i];
    if (track.objectType() == Line)
        drawLine((const TMLine &) track, * _colors[i]);
    else if (track.objectType() == TrackBase)
        drawBase(track, * _colors[i]);
    else if (track.objectType() == Track)
        drawTrack((const TTrack &) track, * _colors[i]);
    else if (track.objectType() == Segment)
        drawSegment((const TSegment &) track, * _colors[i]);
    else
        std::cout << "TWindow::draw !!! can't display" << std::endl;
    }
 
    //...Text...
    _window.draw_text(_xPositionText, _yPositionText, _text.c_str());
//    _window.draw_text(_xPositionText, _yPositionText, (const char *) _text);
}
 
void
TWindow::drawBase(const TTrackBase & base, leda_color c) {
    const AList<TMLink> & links = base.links();
    unsigned n = links.length();
    for (unsigned i = 0; i < n; i++) {
    if (links[i]->wire() == NULL) continue;
    if (! _stereo)
        if (links[i]->wire()->stereo())
        continue;
    if (! _axial)
        if (links[i]->wire()->axial())
        continue;
 
    //...s-z mode...
    if (_mode == 2) {
 
        //...Points...
        HepPoint3D x = links[i]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name() + ":";
 
        if (links[i]->hit()->mc())
            n += "("
              + itostring(links[i]->hit()->mc()->hep()->id()) + ")";
            
        n += itostring(links[i]->leftRight());
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
    }
 
    //...Rphi mode...
    else if (_mode == 0) {
 
        //...Points...
        HepPoint3D x = links[i]->wire()->forwardPosition();
        double radius = links[i]->hit()->drift();
        _window.draw_circle(x.x(), x.y(), radius, c);
        if (_wireName)
        _window.draw_text(x.x(),
                  x.y(),
                  wireName(* links[i]).c_str());
    }
 
    //...Conformal mode...
    else if (_mode == 1) {
 
        //...Transformation...
        AList<TMDCWireHit> list;
        list.append((TMDCWireHit *) links[i]->hit());
        AList<TMLink> list2;
        TConformalFinder0::conformalTransformationRphi(ORIGIN, list,list2);
 
        //...Points...
        HepPoint3D x = list2[0]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name();
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
    }
 
    //...Conformal mode...
    else if (_mode == 3) {
 
        //...Points...
        HepPoint3D x = links[i]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name();
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
    }
    }
}
 
void
TWindow::drawSegment(const TSegment & base, leda_color c) {
    AList<TMLink> links = base.links();
    links.sort(SortByWireId);
    unsigned n = links.length();
    HepPoint3D x;
    HepPoint3D y;
    for (unsigned i = 0; i < n; i++) {
    if (! _stereo)
        if (links[i]->wire()->stereo())
        return;
    if (! _axial)
        if (links[i]->wire()->axial())
        return;
 
    //...s-z mode...
    if (_mode == 2) {
        x = links[i]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name() + ":"
            + itostring(links[i]->leftRight());
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
        if (i)
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);
        y = x;
    }
 
    //...Rphi mode...
    else if (_mode == 0) {
        x = links[i]->wire()->forwardPosition();
        double radius = links[i]->hit()->drift();
        _window.draw_circle(x.x(), x.y(), radius, c);
        if (_wireName)
        _window.draw_text(x.x(),
                  x.y(),
                  (const char *) wireName(* links[i]).c_str());
        if (i)
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);
        y = x;
    }
 
    //...Conformal mode...
    else if (_mode == 1) {
 
        //...Transformation...
        AList<TMDCWireHit> list;
        list.append((TMDCWireHit *) links[i]->hit());
        AList<TMLink> list2;
        TConformalFinder0::conformalTransformationRphi(ORIGIN, list,list2);
 
        //...Points...
        HepPoint3D x = list2[0]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name();
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
        if (i)
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);
        y = x;
        delete list2[0];
    }
 
    //...Conformal mode...
    else if (_mode == 3) {
 
        //...Transformation...
        AList<TMDCWireHit> list;
        list.append((TMDCWireHit *) links[i]->hit());
        AList<TMLink> list2;
        TConformalFinder0::conformalTransformation(ORIGIN, list, list2);
 
        //...Points...
        HepPoint3D x = list2[0]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name();
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
        if (i)
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);
        y = x;
        delete list2[0];
    }
    }
 
    if (_segmentProperty) {
    if (_mode == 3) {
        float in = links.first()->position().mag();
        float out = links.last()->position().mag();
        float length = out - in;
        HepPoint3D x = base.position() - 0.5 * length * base.direction();
        HepPoint3D y = base.position() + 0.5 * length * base.direction();
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);        
        _window.draw_point(base.position().x(), base.position().y(), c);
    }
    }
 
    //...Links to other segments...
    if (! _segmentLink) return;
    const AList<TSegment> & innerLinks = base.innerLinks();
    unsigned nLinks = innerLinks.length();
    for (unsigned i = 0; i < nLinks; i++) {
    if (i == 0) c = leda_blue;
    else        c = leda_pink;
 
    AList<TMLink> innerTLinks = innerLinks[i]->links();
    innerTLinks.sort(SortByWireId);
 
    //...s-z mode...
    if (_mode == 2) {
    }
 
    //...Rphi mode...
    else if (_mode == 0) {
        x = links[0]->wire()->forwardPosition();
        y = innerTLinks.last()->wire()->forwardPosition();
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);
//          if (base.state() & TSegmentCrowd)
//          _window.draw_text(x.x(), x.y(), "c");
    }
 
    //...Conformal mode...
    else if (_mode == 1) {
 
        //...Transformation...
        AList<TMDCWireHit> list;
        list.append((TMDCWireHit *) links[0]->hit());
        list.append((TMDCWireHit *) innerTLinks.last()->hit());
        AList<TMLink> list2;
        TConformalFinder0::conformalTransformation(ORIGIN, list, list2);
 
 
        x = list2[0]->position();
        y = list2[1]->position();
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);
//          if (base.state() & TSegmentCrowd)
//          _window.draw_text(x.x(), x.y(), "c");
    }
 
    //...Conformal mode...
    else if (_mode == 3) {
 
        //...Transformation...
        AList<TMDCWireHit> list;
        list.append((TMDCWireHit *) links[0]->hit());
        list.append((TMDCWireHit *) innerTLinks.last()->hit());
        AList<TMLink> list2;
        TConformalFinder0::conformalTransformation(ORIGIN, list, list2);
 
 
        x = list2[0]->position();
        y = list2[1]->position();
        _window.draw_segment(y.x(), y.y(), x.x(), x.y(), c);
//          if (base.state() & TSegmentCrowd)
//          _window.draw_text(x.x(), x.y(), "c");
    }
    }
}
 
void
TWindow::drawLine(const TMLine & line, leda_color c) {
    drawAxes();
    drawBase((const TTrackBase &) line, c);
 
    //...Line...
    double xmin = _window.xmin();
    double xmax = _window.xmax();
    double ymin = xmin * line.a() + line.b();
    double ymax = xmax * line.a() + line.b();
    _window.draw_segment(xmin, ymin, xmax, ymax, c);
}
 
bool
TWindow::target(void) const {
    struct belle_event * ev = 
    (struct belle_event *) BsGetEnt(BELLE_EVENT, 1, BBS_No_Index);
    if (_targetRun == 0 && _targetEvent == 0)
    return true;
    if (ev->m_RunNo == _targetRun && ev->m_EvtNo == _targetEvent)
    return true;
    return false;
}
 
void
TWindow::target(unsigned run, unsigned farm, unsigned event) {
    _targetRun = run;
    _targetEvent = (farm << 28) + event;
}
 
void
TWindow::drawTrack(const TTrack & t, leda_color c) {
    std::string p = TrackKinematics(t.helix()) + " " + TrackLayerUsage(t);
    _window.draw_text(_xPositionText, _yPositionText, p.c_str(), c);
    _yPositionText += _yPositionStep;
 
    const AList<TMLink> & links = t.links();
    unsigned n = links.length();
    for (unsigned i = 0; i < n; i++) {
    if (! _stereo)
        if (links[i]->wire()->stereo())
        continue;
    if (! _axial)
        if (links[i]->wire()->axial())
        continue;
 
        if (_mode == 0) {
        HepPoint3D x = links[i]->wire()->forwardPosition();
        double radius = links[i]->hit()->drift();
        _window.draw_circle(x.x(), x.y(), radius, c);
        if (_wireName)
        _window.draw_text(x.x(),
                  x.y(),
                  (const char *) wireName(* links[i]).c_str());
    }
 
    else if (_mode == 2) {
        HepPoint3D x = links[i]->positionOnTrack();
        HepPoint3D sz;
        t.szPosition(x, sz);
        _window.draw_point(sz.x(), sz.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name() + ":"
            + itostring(links[i]->leftRight());
        _window.draw_text(sz.x(), sz.y(), n.c_str());
        }
        x = links[i]->positionOnWire();
        sz;
        t.szPosition(x, sz);
        _window.draw_point(sz.x(), sz.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name() + ":"
            + itostring(links[i]->leftRight());
        _window.draw_text(sz.x(), sz.y(), n.c_str());
        }
    }
 
    //...Conformal mode...
    else if (_mode == 1) {
 
        //...Transformation...
        AList<TMDCWireHit> list;
        list.append((TMDCWireHit *) links[i]->hit());
        AList<TMLink> list2;
        TConformalFinder0::conformalTransformationRphi(ORIGIN, list,list2);
 
        //...Points...
        HepPoint3D x = list2[0]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name();
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
        delete list2[0];
    }
 
    //...Conformal mode...
    else if (_mode == 3) {
 
        //...Transformation...
        AList<TMDCWireHit> list;
        list.append((TMDCWireHit *) links[i]->hit());
        AList<TMLink> list2;
        TConformalFinder0::conformalTransformation(ORIGIN, list, list2);
 
        //...Points...
        HepPoint3D x = list2[0]->position();
        _window.draw_point(x.x(), x.y(), c);
        if (_wireName) {
        std::string n = links[i]->wire()->name();
        _window.draw_text(x.x(), x.y(), n.c_str());
        }
        delete list2[0];
    }
    }
 
    if (! n) {
    _window.draw_text(0., 0., (const char *) "can't display a track");
    return;
    }
 
    if (_mode == 0) {
    
    //...Parameters...
    Helix hIp = t.helix();
    hIp.pivot(ORIGIN);
    leda_point ip(hIp.x(0.).x(), hIp.x(0.).y());
    const HepPoint3D & pIn = InnerMost(t.cores())->positionOnTrack();
    const HepPoint3D & pOut = OuterMost(t.cores())->positionOnTrack();
    leda_point in(pIn.x(), pIn.y());
    leda_point out(pOut.x(), pOut.y());
    if (in == out) {
        _window.draw_text(0., 0., (const char *) "can't display a track");
        return;
    }
    _window.draw_arc(ip, in, out, c);
    _window.draw_text(pOut.x(), pOut.y(), t.name().c_str());
    }
 
    else if (_mode == 2) {
    Helix hIp = t.helix();
    hIp.pivot(ORIGIN);
 
    double xmin = _window.xmin();
    double xmax = _window.xmax();
    double ymin = xmin * hIp.tanl() * t.charge() + hIp.dz();
    double ymax = xmax * hIp.tanl() * t.charge() + hIp.dz();
    _window.draw_segment(xmin, ymin, xmax, ymax, c);
    }
}
 
void
TWindow::drawAxes(void) {
    double xmin = _window.xmin();
    double xmax = _window.xmax();
    _window.draw_segment(xmin, 0., xmax, 0., leda_grey2);
    double ymin = _window.ymin();
    double ymax = _window.ymax();
    _window.draw_segment(0., ymin, 0., ymax, leda_grey2);
}
 
void
TWindow::drawCdc(void) {
 
    //...Rphi mode...
    if (_mode == 0) {
    double r = 8.4;
    _window.draw_circle(0., 0., r);
    r = 88.;
    _window.draw_circle(0., 0., r);
    r = 2.;
    _window.draw_circle(0., 0., r);
    }
}
 
void
TWindow::clear(void) {
    ++_nEvents;
    // _skip = false;
    _objects.removeAll();
    if (_canDelete) {
    HepAListDeleteAll(_colors);
    HepAListDeleteAll(_selfObjects);
    HepAListDeleteAll(_selfTLinks);
    }
}
 
void
TWindow::append(const AList<TMLink> & list, leda_color c) {
    TTrackBase * t = new TTrackBase(list);
    _objects.append(t);
    _colors.append(new leda_color(c));
    _selfObjects.append(t);
}
 
std::string
TWindow::wireName(const TMLink & l) const {
    unsigned state = l.hit()->state();
    std::string flag;
    if (state & WireHitFindingValid) flag += "o";
    if (state & WireHitFittingValid) flag += "+";
    if (state & WireHitInvalidForFit) flag += "x";
    if (l.hit()->mc())
    flag += "(" + itostring(l.hit()->mc()->hep()->id()) + ")";
    return l.wire()->name() + flag;
}
 
void
TWindow::oneShot(const TTrackBase & t, leda_color c) {
    append(t, c);
    wait();
    unsigned id = _objects.length() - 1;
    _objects.remove(id);
    leda_color * b = _colors[id];
    _colors.remove(id);
    delete b;
}
 
 
void
TWindow::oneShot(const AList<TSegment> & t, leda_color c) {
    TTrackBase tmp;
    for (unsigned i = 0; i < t.length(); i++) {
    const TSegment & s = * t[i];
    for (unsigned j = 0; j < s.links().length(); j++)
        tmp.append(* s.links()[j]);
    }
 
    append(tmp, c);
    wait();
    unsigned id = _objects.length() - 1;
    _objects.remove(id);
    leda_color * b = _colors[id];
    _colors.remove(id);
    delete b;
}
 
void
TWindow::appendSz(const TTrack & t, const AList<TMLink> & list, leda_color c) {
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink * l0 = new TMLink(* list[i]);
    TMLink * l1 = new TMLink(* list[i]);
    TMLink * l2 = new TMLink(* list[i]);
    l0->leftRight(WireHitLeft);
    l1->leftRight(WireHitRight);
    l2->leftRight(2);
    int err = t.szPosition(* l0);
    err = t.szPosition(* l1);
    err = t.szPosition(* l2);
    AList<TMLink> links;
    links.append(l0);
    links.append(l1);
    links.append(l2);
    TTrackBase * b = new TTrackBase(links);
    _objects.append(b);
    _colors.append(new leda_color(c));
    _selfTLinks.append(l0);
    _selfTLinks.append(l1);
    _selfTLinks.append(l2);
    _selfObjects.append(b);
    }
}
 
void
TWindow::appendSz(const TTrack & t,
          const AList<TSegment> & list,
          leda_color c) {
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
    TSegment & s = * list[i];
    TMLink * l = new TMLink();
    t.szPosition(s, * l);
    AList<TMLink> links;
    links.append(l);
    TTrackBase * b = new TTrackBase(links);
    _objects.append(b);
    _colors.append(new leda_color(c));
    _selfTLinks.append(l);
    _selfObjects.append(b);
    }
}
 
void
TWindow::draw(const TPoint2D & p, leda_color c) {
    _window.draw_point(p.x(), p.y(), c);
}
 
void
TWindow::drawHeader(void) {
    struct belle_event * h = (struct belle_event *)
    BsGetEnt(BELLE_EVENT, 1, BBS_No_Index);
    if (! h) return;
 
    //static const Hepstd::string sp = " ";
    static const std::string sp = " ";
    std::string mc = "mc:" + itostring(h->m_ExpMC);
    std::string exp = "exp:" + itostring(h->m_ExpNo);
    std::string run = "run:" + itostring(h->m_RunNo);
    std::string frm = "frm:" + itostring((h->m_EvtNo >> 28));
    std::string evt = "evt:" + itostring((h->m_EvtNo & 268435455));
    std::string dat = "dat:" + std::string(ctime((const time_t *) & h->m_Time));
    std::string fld = "fld:" + dtostring(h->m_BField);
    std::string ler = "ler:" + dtostring(h->m_ELER);
    std::string her = "her:" + dtostring(h->m_EHER);
 
    std::string s0 = "";
    if (calcdc_unpack_.unpack_truncate)
    s0 += "truncated event";
    std::string head0 = exp + sp + run + sp + frm + sp + evt;
    std::string head1 = dat;
    std::string head2 = fld + sp + ler + sp + her;
 
    double xmin = _window.xmin();
    double xmax = _window.xmax();
    double ymin = _window.ymin();
    double ymax = _window.ymax();
    double x0 = xmin + (xmax - xmin) * .7;
    double x1 = xmin + (xmax - xmin) * .6;
    double y0 = ymin + (ymax - ymin) * .92;
    double y1 = ymin + (ymax - ymin) * .95;
    double y2 = ymin + (ymax - ymin) * .98;
    double y3 = ymin + (ymax - ymin) * -.90;
 
    _window.draw_text(x0, y2, head2.c_str());
    _window.draw_text(x0, y1, head1.c_str());
    _window.draw_text(x0, y0, head0.c_str());
 
 
    _window.draw_text(x1, y3, s0.c_str());
}
 
#endif //TRKRECO_WINDOW