TMDCWireHit.cxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// $Id: TMDCWireHit.cxx,v 1.7 2010/03/31 09:58:59 liucy Exp $
//-----------------------------------------------------------------------------
// Filename : TMDCWireHit.cc
// Section  : Tracking MDC
// Owner    : Yoshi Iwasaki
// Email    : [email protected]
//-----------------------------------------------------------------------------
// Description : A class to represent a wire hit in MDC.
//               See http://bsunsrv1.kek.jp/~yiwasaki/tracking/
//-----------------------------------------------------------------------------
 
#define TMDCWireHit_INLINE_DEFINE_HERE
 
#include "TrkReco/TMDCUtil.h"
#include "TrkReco/TMDCWire.h"
#include "TrkReco/TMDCWireHit.h"
#include "TrkReco/TMDCWireHitMC.h"
#include "TrkReco/TTrackHEP.h"
 
//#include "panther/panther.h"
//#include MDC_H
#include "MdcTables/MdcTables.h"
 
/*
// inline to here --
const TTrack * const
TMDCWireHit::track(void) const {
    return _track;
}
 
const TTrack * const
TMDCWireHit::track(const TTrack * a) {
    return _track = a;
}
 
const TTrack * const
TMDCWireHit::track(const TTrack * a) const {
    return _track = a;
}
// end -- zsl
*/
 
TMDCWireHit::TMDCWireHit(TMDCWire * w, MdcRec_wirhit * r, float fudgeFactor)
: _wire(w),
  _xyPosition(w->xyPosition()),
  _r(r),
  _state(r->stat),
  _track(0),
  _mc(0) {
    w->hit(this);
    _drift[0] = r->ddl;
    _drift[1] = r->ddr;
    _driftError[0] = r->erddl * fudgeFactor;
    _driftError[1] = r->erddr * fudgeFactor;
    if (w->axial()) _state |= WireHitAxial;
    else            _state |= WireHitStereo;
}
 
TMDCWireHit::~TMDCWireHit() {
}
 
void
TMDCWireHit::dump(const std::string & msg, const std::string & pre) const {
    std::cout << pre;
    std::cout << _wire->name();
    if (msg.find("state") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << ",state ";
    bitDisplay(_state);
    if (track()) std::cout << ",trk ";
    }
    if (msg.find("drift") != std::string::npos || msg.find("detail") != std::string::npos) {
    if (_state & WireHitLeftMask) std::cout << ", L";
    if (_state & WireHitRightMask) std::cout << ", R";
    std::cout<<" layer "<<_wire->layerId()<<" cell "<<_wire->localId();
    std::cout << ",dl " << _drift[0] << "+-" << _driftError[0];
    std::cout << ",dr " << _drift[1] << "+-" << _driftError[1];
    }
    if (msg.find("mc") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << ",hep ";
    if (mc()) {
        if (mc()->hep()) std::cout << mc()->hep()->id();
        else std::cout << "0";
    }
    else {
        std::cout << "0";
    }
    }
    std::cout << std::endl;
}
 
HepPoint3D
TMDCWireHit::position(unsigned lr) const {
    const HepPoint3D myHepZHat(0.0, 0.0, 1.0);
    //...Left...
    if (lr == WireHitLeft) {
    return _xyPosition
        - _drift[WireHitLeft] *  myHepZHat.cross(_xyPosition.unit());
    }
 
    //...Right case...
    else {
    return _xyPosition
        + _drift[WireHitRight] * myHepZHat.cross(_xyPosition.unit());
    }
}
 
 
/*
double
TMDCWireHit::DriftTime(double tof,double z) const {
    double tprop = 0.;
    double _vprop = (_layer<8) ? Constants::vpropInner : Constants::vpropOuter;
    if (0 == _layer%2){
      tprop = (0.5*_zlen + z)/_vprop; //odd
    }else{
      tprop = (0.5*_zlen - z)/_vprop; //even
    }
    double driftT;
    driftT = fabs(_rawTime - _T0Walk -1.e9*tof - tprop);
 
    //std::cout<< "lay "<<_layer<<" cell "<<_wire<<" zhit "<<z<<" tprop "<<tprop << std::endl;
    return driftT;
}
*/