/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/TrkReco/TrkReco-00-08-59-patch4-slc6tag/src/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    : yoshihito.iwasaki@kek.jp
//-----------------------------------------------------------------------------
// 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;
}
*/

---