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    : [email protected]
//-----------------------------------------------------------------------------
// 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 "MdcGeomSvc/MdcGeoLayer.h"
#include "TrkReco/TMDC.h"
#include "TrkReco/TMDCLayer.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];
}