TMDCLayer.cxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// $Id: TMDCLayer.cxx,v 1.5 2010/03/31 09:58:59 liucy Exp $
//-----------------------------------------------------------------------------
// Filename : TMDCLayer.cc
// Section  : Tracking MDC
// Owner    : Yoshi Iwasaki
// Email    : yoshihito.iwasaki@kek.jp
//-----------------------------------------------------------------------------
// Description : A class to represent a wire layer.
//               See http://bsunsrv1.kek.jp/~yiwasaki/tracking/
//-----------------------------------------------------------------------------
 
#define TMDCLayer_INLINE_DEFINE_HERE
// #include "panther/panther.h"
// #include MDC_H
// #include "MdcRecGeo/MdcGeoLayer.h"
#include "TrkReco/TMDCLayer.h"
#include "MdcGeomSvc/MdcGeoLayer.h"
#include "TrkReco/TMDC.h"
#include "TrkReco/TMDCUtil.h"
 
TMDCLayer::TMDCLayer( const MdcGeoLayer* l ) {
  _geo          = l;
  _id           = l->Id();
  _superLayerId = TMDC::superLayerId( l );
  _localLayerId = TMDC::localLayerId( l );
  // liucy&jialk origin uses slant to judge type now directly uses superLayerId to judge
  if ( _id >= 36 || ( _id >= 8 && _id <= 19 ) ) _type = 0;
  else _type = 1;
  //    _offset = l->Offset();
  //    _cellSize = l->PCSiz();
  //  Unit from rad to cell for _offset; and from mm to cm for cellSize
  //    _offset = l->Offset()/(2. * M_PI/l->NCell());
  //    unit of offset is rad, which is diff from belle
  _offset   = l->Offset();
  _cellSize = l->PCSiz() / 10;
 
  _nWires             = (unsigned)l->NCell();
  _axialStereoLayerId = TMDC::axialStereoLayerId( l );
}
 
TMDCLayer::~TMDCLayer() {}
 
void TMDCLayer::dump( const std::string& msg, const std::string& pre ) const {
  std::cout << pre;
  std::cout << "layer " << _id;
  std::cout << ", super layer " << _superLayerId;
  std::cout << ", local layer " << _localLayerId;
  if ( axial() ) std::cout << ", axial ";
  else std::cout << ", stereo ";
  std::cout << _axialStereoLayerId;
  std::cout << ", " << _nWires << " wires";
  std::cout << std::endl;
}
 
const TMDCWire* const TMDCLayer::wire( int id ) const {
  if ( _nWires == 0 ) return 0;
 
  if ( id < 0 )
    while ( id < 0 ) id += _nWires;
 
  if ( id >= _nWires ) id %= _nWires;
 
  return ( *this )[id];
}