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 "MdcTables/HepevtTables.h"
#include <vector>
 
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; }