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    : [email protected]
//-----------------------------------------------------------------------------
// 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;
}