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/TMDCWireHit.h"
#include "TrkReco/TMDCUtil.h"
#include "TrkReco/TMDCWire.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;
}
*/