TMDCWire Class Reference

A class to represent a wire in MDC. More...

#include <TMDCWire.h>

List of all members.

Public Member Functions
	TMDCWire (const MdcGeoWire *, TMDCLayer *)
	Constructor.
virtual	~TMDCWire ()
	Destructor.
unsigned	id (void) const
	returns id.
unsigned	localId (void) const
	returns local id in a wire layer.
unsigned	layerId (void) const
	returns layer id.
unsigned	localLayerId (void) const
	returns local layer id in a super layer.
unsigned	superLayerId (void) const
	returns super layer id.
unsigned	axialStereoLayerId (void) const
	returns id of axial or stereo id.
const MdcGeoWire *	geocdc (void) const
	returns a pointer to GEOMDC_WIR.
const TMDCLayer *const	layer (void) const
	returns a pointer to a layer.
const TMDCWireHit *const	hit (void) const
	returns a pointer to a TMDCWireHit.
const TMDCTsf *const	tsf (void) const
	returns a pointer to a TMDCTsf.
unsigned	state (void) const
	returns state.
bool	axial (void) const
	returns true if this wire is in an axial layer.
bool	stereo (void) const
	returns true if this wire is in a stereo layer.
bool	innerPart (void) const
	returns true if this wire is in the inner part.
bool	mainPart (void) const
	returns true if this wire is in the main part.
const TMDCWire *const	neighbor (unsigned) const
	returns a pointer to a neighbor wire.
bool	consective (const TMDCWire &) const
	returns true if a given wire is consective in a layer.
bool	adjacent (const TMDCWire &) const
	returns true if a given wire is adjacent.
int	localIdForPlus (void) const
	returns localId but if maxLocalId, return -1.
int	localIdForMinus (void) const
	returns localId but if 0, return maxLocalId + 1.
const HepPoint3D &	forwardPosition (void) const
	returns position in forward endplate.
const HepPoint3D &	backwardPosition (void) const
	returns position in backward endplate.
double *	backwardPosition (double p[3]) const
const double &	tension (void) const
const HepPoint3D &	xyPosition (void) const
	returns middle position of a wire. z componet is 0.
double *	xyPosition (double p[3]) const
const HepPoint3D	xyPosition (double z) const
const HepVector3D &	direction (void) const
	returns direction vector of the wire.
void	wirePosition (float zPosition, HepPoint3D &xyPosition, HepPoint3D &backwardPosition, HepVector3D &direction) const
	calculates position and direction vector with sag correction.
float	cellSize (void) const
	returns cell size in phi.
int	localIdDifference (const TMDCWire &) const
	returns local id difference.
std::string	name (void) const
	returns name.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
const TMDCWireHit *const	hit (const TMDCWireHit *const)
	sets a pointer to TMDCWireHit.
const TMDCWireHitMC *const	hit (TMDCWireHitMC *const)
	appends a pointer to TMDCWireHitMC.
TMDCTsf *	tsf (TMDCTsf *)
	sets a pointer of TMDCTsf.
unsigned	state (unsigned newState)
	sets state.
void	clear (void)
	clears information.
Private Attributes
unsigned	_state
const TMDCWireHit *	_hit
AList< TMDCWireHitMC >	_mcHits
unsigned	_id
unsigned	_localId
unsigned	_layerId
const TMDCLayer *	_layer
TMDCTsf *	_tsf
HepPoint3D	_xyPosition
HepPoint3D	_forwardPosition
HepPoint3D	_backwardPosition
double	_tension
HepVector3D	_direction
const MdcGeoWire *	_geo
Static Private Attributes
static TMDC *	_cdc = 0

---

Detailed Description

A class to represent a wire in MDC.

Definition at line 55 of file TMDCWire.h.

---

Constructor & Destructor Documentation

TMDCWire::TMDCWire	(	const MdcGeoWire *	,
		TMDCLayer *
	)

Constructor.

Definition at line 32 of file TMDCWire.cxx.

References _backwardPosition, _direction, _forwardPosition, _geo, _hit, _id, _layer, _layerId, _localId, _state, _tsf, _xyPosition, and w.

: _forwardPosition(w->Backward()*0.1), 
  _backwardPosition(w->Forward()*0.1), 
  _tension(w->Tension()),
  _direction(w->Backward() -w->Forward()) { 
//: _forwardPosition(w->Forward()*0.1),
//  _backwardPosition(w->Backward()*0.1),
//  _direction(w->Forward() - w->Backward()) {
    _state = 0;
    _geo = w;
    _id = w->Id();
    _localId = w->Cell();
    _layerId = w->Layer();
    _layer = l;
    _tsf = 0;
    _hit = 0;
    _xyPosition = 0.5 * (_forwardPosition + _backwardPosition);
    _xyPosition.setZ(0.);
    _direction = _direction.unit();
}

TMDCWire::~TMDCWire

(

)

[virtual]

Destructor.

Definition at line 53 of file TMDCWire.cxx.

                    {
}

---

Member Function Documentation

bool TMDCWire::adjacent

(

const TMDCWire &

)

const [inline]

returns true if a given wire is adjacent.

Definition at line 403 of file TMDCWire.h.

References genRecEmupikp::i, neighbor(), and w.

                                           {
    for (unsigned i = 0; i < 6; i++)
        if (neighbor(i) == & w)
            return true;
    return false;
}

bool TMDCWire::axial

(

void

)

const [inline]

returns true if this wire is in an axial layer.

Definition at line 348 of file TMDCWire.h.

References _layer, and TMDCLayer::axial().

Referenced by name().

                          {
    return _layer->axial();
}

unsigned TMDCWire::axialStereoLayerId

(

void

)

const [inline]

returns id of axial or stereo id.

Definition at line 360 of file TMDCWire.h.

References _layer, and TMDCLayer::axialStereoLayerId().

                                       {
    return _layer->axialStereoLayerId();
}

double * TMDCWire::backwardPosition

(

double

p[3]

)

const [inline]

Definition at line 318 of file TMDCWire.h.

References _backwardPosition.

                                            {
    p[0] = _backwardPosition.x();
    p[1] = _backwardPosition.y();
    p[2] = _backwardPosition.z();
    return p;
}

const HepPoint3D & TMDCWire::backwardPosition

(

void

)

const [inline]

returns position in backward endplate.

Definition at line 306 of file TMDCWire.h.

References _backwardPosition.

Referenced by TCurlFinder::distance(), THelixFitter::drift(), TTrack::stereoHitForCurl(), and TTrack::szPosition().

                                     {
    return _backwardPosition;
}

float TMDCWire::cellSize

(

void

)

const [inline]

returns cell size in phi.

Definition at line 380 of file TMDCWire.h.

References _layer, and TMDCLayer::cellSize().

                             {
    if (_layer) {
        return _layer->cellSize();
    }
    else {
#ifdef TRKRECO_DEBUG
        std::cout << "TMDCWire::cellSize !!! can't return cell size";
        std::cout << " because no pointer to a layer" << std::endl;
#endif
        return 0.;
    }
}

void TMDCWire::clear

(

void

)

[inline]

clears information.

Definition at line 273 of file TMDCWire.h.

References _hit, _mcHits, and _state.

                    {
    _state = 0;
    _hit = 0;
    _mcHits.removeAll();
}

bool TMDCWire::consective

(

const TMDCWire &

)

const [inline]

returns true if a given wire is consective in a layer.

Definition at line 395 of file TMDCWire.h.

References neighbor(), and w.

Referenced by TBuilder0::consectiveHits().

                                             {
    if (neighbor(2) == & w) return true;
    else if (neighbor(3) == & w) return true;
    return false;
}

const HepVector3D & TMDCWire::direction

(

void

)

const [inline]

returns direction vector of the wire.

Definition at line 342 of file TMDCWire.h.

References _direction.

Referenced by TCurlFinder::distance(), THelixFitter::main(), TTrack::stereoHitForCurl(), and TTrack::szPosition().

                              {
    return _direction;
}

void TMDCWire::dump	(	const std::string &	message = std::string(""),
		const std::string &	prefix = std::string("")
	)			const

dumps debug information.

Definition at line 57 of file TMDCWire.cxx.

References _id, _layerId, _localId, _tsf, genRecEmupikp::i, localLayerId(), neighbor(), RealDBUtil::npos, and superLayerId().

                                                                 {
    std::cout << pre;
    std::cout << "w " << _id;
    std::cout << ",local " << _localId;
    std::cout << ",layer " << _layerId;
    std::cout << ",super layer " << superLayerId();
    std::cout << ",local layer " << localLayerId();
    std::cout << std::endl;
    if (msg.find("neighbor") != std::string::npos || msg.find("detail") != std::string::npos) {
        for (unsigned i = 0; i < 6; i++)
            if (neighbor(i))
                neighbor(i)->dump("", pre + itostring(i) + "   ");
    }    
    if (msg.find("tsf") != std::string::npos || msg.find("detail") != std::string::npos) {
        if (_tsf) _tsf->dump("", pre + "    ");
    }
}

const HepPoint3D & TMDCWire::forwardPosition

(

void

)

const [inline]

returns position in forward endplate.

Definition at line 300 of file TMDCWire.h.

References _forwardPosition.

Referenced by THelixFitter::drift(), TTrack::stereoHitForCurl(), and TTrack::szPosition().

                                    {
    return _forwardPosition;
}

const MdcGeoWire * TMDCWire::geocdc

(

void

)

const [inline]

returns a pointer to GEOMDC_WIR.

Definition at line 237 of file TMDCWire.h.

References _geo.

                           {
    return _geo;
}

const TMDCWireHitMC *const TMDCWire::hit

(

TMDCWireHitMC *

const

)

[inline]

appends a pointer to TMDCWireHitMC.

Definition at line 293 of file TMDCWire.h.

References _mcHits.

                                     {
    _mcHits.append(a);
    return a;
}

const TMDCWireHit *const TMDCWire::hit

(

const TMDCWireHit *

const

)

[inline]

sets a pointer to TMDCWireHit.

Definition at line 281 of file TMDCWire.h.

References _hit.

                                         {
    return _hit = h;
}

const TMDCWireHit *const TMDCWire::hit

(

void

)

const [inline]

returns a pointer to a TMDCWireHit.

Definition at line 287 of file TMDCWire.h.

References _hit.

Referenced by TMDC::classification().

                        {
    return _hit;
}

unsigned TMDCWire::id

(

void

)

const [inline]

returns id.

Definition at line 207 of file TMDCWire.h.

References _id.

Referenced by TSegment::appendByLine(), TMDC::driftDistance(), TRungeFitter::fit(), TCurlFinder::searchAxialCand(), TCurlFinder::searchHits(), TMDC::updateMC(), and wirePosition().

                       {
    return _id;
}

bool TMDCWire::innerPart

(

void

)

const [inline]

returns true if this wire is in the inner part.

Definition at line 366 of file TMDCWire.h.

References _layerId.

                              {
    if (_layerId < 14) return true;
    return false;
}

const TMDCLayer *const TMDCWire::layer

(

void

)

const [inline]

returns a pointer to a layer.

Definition at line 243 of file TMDCWire.h.

References _layer.

Referenced by TBuilder0::consectiveHits(), and localIdDifference().

                          {
    return _layer;
}

unsigned TMDCWire::layerId

(

void

)

const [inline]

returns layer id.

Definition at line 219 of file TMDCWire.h.

References _layerId.

Referenced by TBuilder0::check2CnHits(), TBuilder0::consectiveHits(), TMDCTsf::createTsf(), THelixFitter::drift(), T3DLineFitter::drift(), TMDCWireHit::dump(), findIsolatedCloseHits(), TRungeFitter::fit(), TMLine::fit2(), TLine0::fit2(), TMLine::fit2p(), TLine0::fit2p(), TTrackManager::maskNormal(), TSegment::maxdDistance(), TrkReco::maxdDistance(), TTrackManager::merge(), neighbor(), offsetBorder(), SameLayer(), TCurlFinder::searchAxialCand(), and TCircle::weight().

                            {
    return _layerId;
}

unsigned TMDCWire::localId

(

void

)

const [inline]

returns local id in a wire layer.

Definition at line 213 of file TMDCWire.h.

References _localId.

Referenced by TBuilder0::buildStereo(), TBuilder0::check2CnHits(), TMDCTsf::createTsf(), THelixFitter::drift(), T3DLineFitter::drift(), TMDCWireHit::dump(), findIsolatedCloseHits(), TRungeFitter::fit(), localIdDifference(), neighbor(), TCurlFinder::searchAxialCand(), TSegment0::splitComplicated(), TSegment::splitComplicated(), and TCircle::weight().

                            {
    return _localId;
}

int TMDCWire::localIdDifference

(

const TMDCWire &

)

const

returns local id difference.

Definition at line 577 of file TMDCWire.cxx.

References layer(), localId(), TMDCLayer::nWires(), and superLayerId().

Referenced by TSegment0::splitParallel(), and TSegment::splitParallel().

                                                    {

#ifdef TRKRECO_DEBUG_DETAIL
    if (superLayerId() != a.superLayerId()) {
        std::cout << "TMDCWire::localIdDifference !!!";
        std::cout << "super layer assumption violation" << std::endl;
    }
#endif

    int diff = int(a.localId()) - int(localId());
    unsigned nWires = layer()->nWires();
    if (diff > 0) {
        int difR = nWires - diff;
        if (diff < difR) return diff;
        else return - difR;
    }
    else {
        int difR = nWires + diff;
        if (- diff < difR) return diff;
        else return difR;
    }
}

int TMDCWire::localIdForMinus

(

void

)

const

returns localId but if 0, return maxLocalId + 1.

Definition at line 355 of file TMDCWire.cxx.

References _layerId, and _localId.

                                    {
/*
  if(_layerId >= 0 && _layerId <= 5){
    if(_localId == 0)
      return 64;
    else return _localId;
  }

  if(_layerId >= 6 && _layerId <= 8){
    if(_localId == 0)
      return 80;
    else return _localId;
  }

  if(_layerId >= 9 && _layerId <= 14){
    if(_localId == 0)
      return 96;
    else return _localId;
  }
  
  if(_layerId >= 15 && _layerId <= 17){
    if(_localId == 0)
      return 128;
    else return _localId;
  }

  if(_layerId >= 18 && _layerId <= 22){
    if(_localId == 0)
      return 144;
    else return _localId;
  }

  if(_layerId >= 23 && _layerId <= 26){
    if(_localId == 0)
      return 160;
    else return _localId;
  }
  
  if(_layerId >= 27 && _layerId <= 31){
    if(_localId == 0)
      return 192;
    else return _localId;
  }
  
  if(_layerId >= 32 && _layerId <= 35){
    if(_localId == 0)
      return 208;
    else return _localId;
  }

  if(_layerId >= 36 && _layerId <= 40){
    if(_localId == 0)
      return 240;
    else return _localId;
  }

  if(_layerId >= 41 && _layerId <= 44){
    if(_localId == 0)
      return 256;
    else return _localId;
  }

  if(_layerId >= 45){
    if(_localId == 0)
      return 288;
    else return _localId;
  }
*/
//Liuqg 060915
  if(_layerId == 0){
    if(_localId == 0)
      return 40;
    else return _localId;
  }

  if(_layerId == 1){
    if(_localId == 0)
      return 44;
    else return _localId;
  }

  if(_layerId == 2){
    if(_localId == 0)
      return 48;
    else return _localId;
  }

  if(_layerId == 3){
    if(_localId == 0)
      return 56;
    else return _localId;
  }

  if(_layerId == 4){
    if(_localId == 0)
      return 64;
    else return _localId;
  }

  if(_layerId == 5){
    if(_localId == 0)
      return 72;
    else return _localId;
  }

  if(_layerId == 6 || _layerId == 7){
    if(_localId == 0)
      return 80;
    else return _localId;
  }

  if(_layerId == 8 || _layerId == 9){
    if(_localId == 0)
      return 76;
    else return _localId;
  }

  if(_layerId == 10 || _layerId == 11){
    if(_localId == 0)
      return 88;
    else return _localId;
  }

  if(_layerId == 12 || _layerId == 13){
    if(_localId == 0)
      return 100;
    else return _localId;
  }

  if(_layerId == 14 || _layerId == 15){
    if(_localId == 0)
      return 112;
    else return _localId;
  }

  if(_layerId == 16 || _layerId == 17){
    if(_localId == 0)
      return 128;
    else return _localId;
  }

  if(_layerId == 18 || _layerId == 19){
    if(_localId == 0)
      return 140;
    else return _localId;
  }

  if(_layerId >= 20 && _layerId <= 23){
    if(_localId == 0)
      return 160;
    else return _localId;
  }

  if(_layerId >= 24 && _layerId <= 27){
    if(_localId == 0)
      return 176;
    else return _localId;
  }

  if(_layerId >= 28 && _layerId <= 31){
    if(_localId == 0)
      return 208;
    else return _localId;
  }

  if(_layerId >= 32 && _layerId <= 35){
    if(_localId == 0)
      return 240;
    else return _localId;
  }

  if(_layerId >= 36 && _layerId <= 39){
    if(_localId == 0)
      return 256;
    else return _localId;
  }

  if(_layerId >= 40 && _layerId <= 42){
    if(_localId == 0)
      return 288;
    else return _localId;
  }
  return -1;
}

int TMDCWire::localIdForPlus

(

void

)

const

returns localId but if maxLocalId, return -1.

Definition at line 166 of file TMDCWire.cxx.

References _layerId, and _localId.

                                   {
/*
  if(_layerId >= 0 && _layerId <= 5){
    if(_localId == 63)
      return -1;
    else return _localId;
  }

  if(_layerId >= 6 && _layerId <= 8){
    if(_localId == 79)
      return -1;
    else return _localId;
  }

  if(_layerId >= 9 && _layerId <= 14){
    if(_localId == 95)
      return -1;
    else return _localId;
  }
  
  if(_layerId >= 15 && _layerId <= 17){
    if(_localId == 127)
      return -1;
    else return _localId;
  }

  if(_layerId >= 18 && _layerId <= 22){
    if(_localId == 143)
      return -1;
    else return _localId;
  }

  if(_layerId >= 23 && _layerId <= 26){
    if(_localId == 159)
      return -1;
    else return _localId;
  }
  
  if(_layerId >= 27 && _layerId <= 31){
    if(_localId == 175)
      return -1;
    else return _localId;
  }
  
  if(_layerId >= 32 && _layerId <= 35){
    if(_localId == 207)
      return -1;
    else return _localId;
  }

  if(_layerId >= 36 && _layerId <= 40){
    if(_localId == 239)
      return -1;
    else return _localId;
  }

  if(_layerId >= 41 && _layerId <= 44){
    if(_localId == 255)
      return -1;
    else return _localId;
  }

  if(_layerId >= 45){
    if(_localId == 287)
      return -1;
    else return _localId;
  }

  return -1;
*/
//Liuqg 060915
  if(_layerId == 0){
    if(_localId == 39)
      return -1;
    else return _localId;
  }

  if(_layerId == 1){
    if(_localId == 43)
      return -1;
    else return _localId;
  }

  if(_layerId == 2){
    if(_localId == 47)
      return -1;
    else return _localId;
  }

  if(_layerId == 3){
    if(_localId == 55)
      return -1;
    else return _localId;
  }

  if(_layerId == 4){
    if(_localId == 63)
      return -1;
    else return _localId;
  }

  if(_layerId == 5){
    if(_localId == 71)
      return -1;
    else return _localId;
  }

  if(_layerId == 6 || _layerId == 7){
    if(_localId == 79)
      return -1;
    else return _localId;
  }

  if(_layerId == 8 || _layerId == 9){
    if(_localId == 75)
      return -1;
    else return _localId;
  }

  if(_layerId == 10 || _layerId == 11){
    if(_localId == 87)
      return -1;
    else return _localId;
  }

  if(_layerId == 12 || _layerId == 13){
    if(_localId == 99)
      return -1;
    else return _localId;
  }

  if(_layerId == 14 || _layerId == 15){
    if(_localId == 111)
      return -1;
    else return _localId;
  }

  if(_layerId == 16 || _layerId == 17){
    if(_localId == 127)
      return -1;
    else return _localId;
  }

  if(_layerId == 18 || _layerId == 19){
    if(_localId == 139)
      return -1;
    else return _localId;
  }

  if(_layerId >= 20 && _layerId <= 23){
    if(_localId == 159)
      return -1;
    else return _localId;
  }

  if(_layerId >= 24 && _layerId <= 27){
    if(_localId == 175)
      return -1;
    else return _localId;
  }
  
  if(_layerId >= 28 && _layerId <= 31){
    if(_localId == 207)
      return -1;
    else return _localId;
  }

  if(_layerId >= 32 && _layerId <= 35){
    if(_localId == 239)
      return -1;
    else return _localId;
  }

  if(_layerId >= 36 && _layerId <= 39){
    if(_localId == 255)
      return -1;
    else return _localId;
  }
  
  if(_layerId >= 40 && _layerId <= 42){
    if(_localId == 287)
      return -1;
    else return _localId;
  }

  return -1;
}

unsigned TMDCWire::localLayerId

(

void

)

const [inline]

returns local layer id in a super layer.

Definition at line 231 of file TMDCWire.h.

References _layer, and TMDCLayer::localLayerId().

Referenced by TSegmentCurl::append(), dump(), TSegmentCurl::remove(), and TSegment::splitTsf().

                                 {
    return _layer->localLayerId();
}

bool TMDCWire::mainPart

(

void

)

const [inline]

returns true if this wire is in the main part.

Definition at line 373 of file TMDCWire.h.

References _layerId.

                             {
    if (_layerId > 13) return true;
    return false;
}

std::string TMDCWire::name

(

void

)

const [inline]

returns name.

Definition at line 412 of file TMDCWire.h.

References _layerId, _localId, axial(), and deljobs::string.

Referenced by TTrackBase::append(), TTrack::assign(), TMDCWireHit::dump(), TSegment0::splitComplicated(), and TSegment::splitComplicated().

                         {
    if (axial())
        return itostring(_layerId) + std::string("-") + itostring(_localId);
    return itostring(_layerId) + std::string("=") + itostring(_localId);
}

const TMDCWire *const TMDCWire::neighbor

(

unsigned

)

const

returns a pointer to a neighbor wire.

Definition at line 95 of file TMDCWire.cxx.

References _cdc, TMDC::getTMDC(), TMDC::layer(), layerId(), localId(), TMDC::nLayers(), TMDCLayer::nWires(), TMDCLayer::offset(), phi0, pi, s, TMDCLayer::superLayerId(), superLayerId(), TMDC::wire(), WireInnerLeft, WireInnerRight, WireLeft, WireOuterLeft, WireOuterRight, and WireRight.

Referenced by adjacent(), TMDC::classification(), consective(), and dump().

                                   {
    if (_cdc == 0) _cdc = TMDC::getTMDC();
    static const unsigned nLayers = _cdc->nLayers();

    //...Expensive...
    unsigned l = layerId();
    int local = (int) localId();
    float phi0 = _cdc->layer(l)->offset();
    int nWir = (int) _cdc->layer(l)->nWires();
    float phi = phi0 + local*2*pi/nWir;
    unsigned s = superLayerId();    
//    cout<<"i:"<<i<<"   layerId:"<<l<<"  local:"<<local<<"  nWire:"<<nWir<<"   superLyrId:"<<s<<endl;

    if (i == WireInnerLeft || i == WireInnerRight) {
        if (l == 0) return 0;
        --l;
        if (_cdc->layer(l)->superLayerId() != s) return 0;
        phi0 = _cdc->layer(l)->offset();
        nWir = _cdc->layer(l)->nWires();
        local = (int)((phi-phi0)/(2*pi/nWir)); 

//      cout<<"inner local"<<local<<", "<<local+1<<endl;
        
//      cout<<(phi-phi0)/(2*pi/nWir)<<"  local:"<<local<<endl;
//      cout<<"InnerLeft phi:"<<phi0 + local*2*pi/nWir
//              <<"   InnerRight phi:"<<phi0 + (local+1)*2*pi/nWir 
//              <<"   phi:"<<phi<<endl;
//      if(local<0 || local+1>nWir-1) 
//              cout<<"TMDCWire::neighbor()---Pay attation the cellNo !!!"<<endl
//                  <<l+1<<"th layer: InnerLeftRight,   local:"<<local<<"    local+1:"<<local+1<<endl;  
// LR changed to bes3
        if (i == WireInnerLeft) return _cdc->wire(l, local);
        else                    return _cdc->wire(l, local + 1);
//      if (i == WireInnerRight) return _cdc->wire(l, local);
//      else                     return _cdc->wire(l, local + 1);
        
    }
    else if (i == WireLeft || i == WireRight) {
//        cout<<"Left phi:"<<phi0 + (local-1)*2*pi/nWir  
//            <<"   Right phi:"<<phi0 + (local+1)*2*pi/nWir 
//            <<"   phi:"<<phi<<endl;
//      if(local-1<0 || local+1>nWir-1) 
//              cout<<"TMDCWire::neighbor()---Pay attation the Cell No !!!"<<endl
//                <<l<<"th layer: LeftRight,   local-1:"<<local-1<<"    local+1:"<<local+1<<endl;       
        if (i == WireLeft) return _cdc->wire(l, local - 1);
        else                return _cdc->wire(l, local + 1);
    }
    else if (i == WireOuterLeft || i == WireOuterRight) {
        if (l == nLayers - 1) return 0;
        ++l;
        if (_cdc->layer(l)->superLayerId() != s) return 0;
        phi0 = _cdc->layer(l)->offset();
        nWir = _cdc->layer(l)->nWires();
        local = (int)((phi-phi0)/(2*pi/nWir));
 
//      cout<<"outer local"<<local<<", "<<local+1<<endl;
                
//      cout<<(phi-phi0)/(2*pi/nWir)<<"  local:"<<local<<endl;          
//        cout<<"OuterLeft phi:"<<phi0 + local*2*pi/nWir 
//                <<"   OuterRight phi:"<<phi0 + (local+1)*2*pi/nWir 
//              <<"   phi:"<<phi<<endl;
//        if(local<0 || local+1>nWir-1) 
//              cout<<"TMDCWire::neighbor()---Pay attation the Cell No !!!"<<endl
//               <<l-1<<"th layer:  OuterLeftRight,   local:"<<local<<"    local+1:"<<local+1<<endl;    
        if (i == WireOuterLeft) return _cdc->wire(l, local);
        else                     return _cdc->wire(l, local + 1);
    }
    return 0;
}

unsigned TMDCWire::state

(

unsigned

newState

)

[inline]

sets state.

Definition at line 267 of file TMDCWire.h.

References _state.

                          {
    return _state = a;
}

unsigned TMDCWire::state

(

void

)

const [inline]

returns state.

Definition at line 261 of file TMDCWire.h.

References _state.

                          {
    return _state;
}

bool TMDCWire::stereo

(

void

)

const [inline]

returns true if this wire is in a stereo layer.

Definition at line 354 of file TMDCWire.h.

References _layer, and TMDCLayer::stereo().

                           {
    return _layer->stereo();
}

unsigned TMDCWire::superLayerId

(

void

)

const [inline]

returns super layer id.

Definition at line 225 of file TMDCWire.h.

References _layer, and TMDCLayer::superLayerId().

Referenced by TSegment::appendByLine(), TBuilder0::check2CnHits(), TCurlFinder::distance(), dump(), localIdDifference(), TSegment::maxdDistance(), TrkReco::maxdDistance(), neighbor(), and SameSuperLayer().

                                 {
    return _layer->superLayerId();
}

const double & TMDCWire::tension

(

void

)

const [inline]

Definition at line 312 of file TMDCWire.h.

References _tension.

                            {
        return _tension;
}

TMDCTsf * TMDCWire::tsf

(

TMDCTsf *

)

[inline]

sets a pointer of TMDCTsf.

Definition at line 255 of file TMDCWire.h.

References _tsf, and t().

                         {
    return _tsf = t;
}

const TMDCTsf *const TMDCWire::tsf

(

void

)

const [inline]

returns a pointer to a TMDCTsf.

Definition at line 249 of file TMDCWire.h.

References _tsf.

                        {
    return _tsf;
}

void TMDCWire::wirePosition	(	float	zPosition,
		HepPoint3D &	xyPosition,
		HepPoint3D &	backwardPosition,
		HepVector3D &	direction
	)			const

calculates position and direction vector with sag correction.

Definition at line 541 of file TMDCWire.cxx.

References _backwardPosition, _direction, _forwardPosition, _xyPosition, and id().

Referenced by TRunge::SetFlightLength().

                                                {
    back = _backwardPosition;

    //...Check z position...
    if (! (z > _backwardPosition.z() && z < _forwardPosition.z())) {
        xy = _xyPosition;
        dir = _direction;
    }

    //...Setup...
    int wireID = id();
    float wirePosition[3] = {0., 0., 0.};
    float dydz = 0;
    float ybSag = 0;
    float yfSag = 0;

    //...Get corrections...
//zsl    calcdc_sag3_(& wireID, & z, wirePosition, & dydz, & ybSag, & yfSag);   
        
    //...Wire position...
    xy.setX((double) wirePosition[0]);
    xy.setY((double) wirePosition[1]);
    xy.setZ((double) wirePosition[2]);
    back.setY((double) ybSag);
    HepVector3D v_aux(_forwardPosition.x() - _backwardPosition.x(),
                   yfSag - ybSag,
                   _forwardPosition.z() - _backwardPosition.z());
    dir = v_aux.unit();

    return;
}

const HepPoint3D TMDCWire::xyPosition

(

double

z

)

const

Definition at line 76 of file TMDCWire.cxx.

References _backwardPosition, _direction, _forwardPosition, _tension, EvtCyclic3::A, and x.

                                   {
    double A=47.35E-6/_tension;
    HepPoint3D x0 =(z - _backwardPosition.z()) /_direction.z() * _direction + _backwardPosition;

    double zinit = x0.z();
    HepVector3D wire;
    wire.setX(_forwardPosition.x()-_backwardPosition.x());
    wire.setY(_forwardPosition.y()-_backwardPosition.y());
    wire.setZ(_forwardPosition.z()-_backwardPosition.z());
    double lxz=sqrt(wire.z()*wire.z()+wire.x()*wire.x());
    double Zp = (zinit - _backwardPosition.z())*lxz/wire.z();
    double x=_direction.x()*(zinit-_backwardPosition.z())/_direction.z() + _backwardPosition.x();
    double y=(A*(Zp-lxz)+wire.y()/lxz)*Zp+_backwardPosition.y();
    x0.setX(x);
    x0.setY(y);
    return x0;
}

double * TMDCWire::xyPosition

(

double

p[3]

)

const [inline]

Definition at line 333 of file TMDCWire.h.

References _xyPosition.

                                      {
    a[0] = _xyPosition.x();
    a[1] = _xyPosition.y();
    a[2] = 0.;
    return a;
}

const HepPoint3D & TMDCWire::xyPosition

(

void

)

const [inline]

returns middle position of a wire. z componet is 0.

Definition at line 327 of file TMDCWire.h.

References _xyPosition.

Referenced by calVirtualCircle(), TPerfectFinder::doit(), TTrack::fit2D(), TCurlFinder::searchHits(), and TMLink::xyPosition().

                               {
    return _xyPosition;
}

---

Member Data Documentation

HepPoint3D TMDCWire::_backwardPosition [private]

Definition at line 186 of file TMDCWire.h.

Referenced by backwardPosition(), TMDCWire(), wirePosition(), and xyPosition().

TMDC * TMDCWire::_cdc = 0 [static, private]

Definition at line 191 of file TMDCWire.h.

Referenced by neighbor().

HepVector3D TMDCWire::_direction [private]

Definition at line 188 of file TMDCWire.h.

Referenced by direction(), TMDCWire(), wirePosition(), and xyPosition().

HepPoint3D TMDCWire::_forwardPosition [private]

Definition at line 185 of file TMDCWire.h.

Referenced by forwardPosition(), TMDCWire(), wirePosition(), and xyPosition().

const MdcGeoWire* TMDCWire::_geo [private]

Definition at line 189 of file TMDCWire.h.

Referenced by geocdc(), and TMDCWire().

const TMDCWireHit* TMDCWire::_hit [private]

Definition at line 176 of file TMDCWire.h.

Referenced by clear(), hit(), and TMDCWire().

unsigned TMDCWire::_id [private]

Definition at line 179 of file TMDCWire.h.

Referenced by dump(), id(), and TMDCWire().

const TMDCLayer* TMDCWire::_layer [private]

Definition at line 182 of file TMDCWire.h.

Referenced by axial(), axialStereoLayerId(), cellSize(), layer(), localLayerId(), stereo(), superLayerId(), and TMDCWire().

unsigned TMDCWire::_layerId [private]

Definition at line 181 of file TMDCWire.h.

Referenced by dump(), innerPart(), layerId(), localIdForMinus(), localIdForPlus(), mainPart(), name(), and TMDCWire().

unsigned TMDCWire::_localId [private]

Definition at line 180 of file TMDCWire.h.

Referenced by dump(), localId(), localIdForMinus(), localIdForPlus(), name(), and TMDCWire().

AList<TMDCWireHitMC> TMDCWire::_mcHits [private]

Definition at line 177 of file TMDCWire.h.

Referenced by clear(), and hit().

unsigned TMDCWire::_state [private]

Definition at line 175 of file TMDCWire.h.

Referenced by clear(), state(), and TMDCWire().

double TMDCWire::_tension [private]

Definition at line 187 of file TMDCWire.h.

Referenced by tension(), and xyPosition().

TMDCTsf* TMDCWire::_tsf [private]

Definition at line 183 of file TMDCWire.h.

Referenced by dump(), TMDCWire(), and tsf().

HepPoint3D TMDCWire::_xyPosition [private]

Definition at line 184 of file TMDCWire.h.

Referenced by TMDCWire(), wirePosition(), and xyPosition().

---