/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/TrkReco/TrkReco-00-08-59-patch4-slc6tag/src/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    : yoshihito.iwasaki@kek.jp
//-----------------------------------------------------------------------------
// 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

---