/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/TrkReco/TrkReco-00-08-59-patch4-slc6tag/src/TTrackHEP.cxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// $Id: TTrackHEP.cxx,v 1.4 2010/03/31 09:58:59 liucy Exp $
//-----------------------------------------------------------------------------
// Filename : TTrackHEP.cc
// Section  : Tracking
// Owner    : Yoshi Iwasaki
// Email    : yoshihito.iwasaki@kek.jp
//-----------------------------------------------------------------------------
// Description : A class to represent a GEN_HEPEVT particle in tracking.
//               See http://bsunsrv1.kek.jp/~yiwasaki/tracking/
//-----------------------------------------------------------------------------

#include "TrkReco/TTrackHEP.h"
#include <vector>
#include "MdcTables/HepevtTables.h"

AList<TTrackHEP>
TTrackHEP::_list = AList<TTrackHEP>();

void
TTrackHEP::update(void) {
    if (_list.length()) HepAListDeleteAll(_list);

//    unsigned n = BsCouTab(GEN_HEPEVT);
//    std::vector<Gen_hepevt>* genhep = GenHepevtCol::getGenHepevtCol();
    unsigned n = GenHepevtCol::getGenHepevtCol()->size();
//    unsigned n = genhep -> size();
   
    for (unsigned i = 0; i < n; i++) {
//      Gen_hepevt * h = *genhep[i];
        Gen_hepevt* h = &(*GenHepevtCol::getGenHepevtCol())[i];
//          (struct gen_hepevt *) BsGetEnt(GEN_HEPEVT, i + 1, BBS_No_Index);
        if (! h) {
            std::cout << "TTrackHEP::update !!! can not access to GEN_HEPEVT";
            std::cout << std::endl;
            break;
        }

        _list.append(new TTrackHEP(h));
    }
}

TTrackHEP::TTrackHEP(const Gen_hepevt * h)
: _hep(h),
  _mother(0),
  _p(h->P[0], h->P[1], h->P[2], h->P[3]),
  _v(h->V[0], h->V[1], h->V[2]) {
    if (_hep->mother != 0) {
//      _mother = _list[_hep->mother-1];
        _mother = _list[_hep->mother->id];
        _mother->_children.append(this);
    }
}

TTrackHEP::~TTrackHEP() {
}

void
TTrackHEP::dump(const std::string & msg, const std::string & pre) const {
    std::cout << pre;

    std::cout << id() << ":";
    std::cout << pType() << ":";
    if (_mother) std::cout << _mother->id();
    else         std::cout << "-";
    std::cout << ":";
    std::cout << _p << ":" << _v;
    std::cout << std::endl;
}

const AList<TTrackHEP> &
TTrackHEP::list(void) {
    return _list;
}

---