TMDCWire Class Reference

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

#include <TMDCWire.h>

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.

Detailed Description

A class to represent a wire in MDC.

Definition at line 55 of file TMDCWire.h.

Constructor & Destructor Documentation

TMDCWire()

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

Constructor.

Definition at line 32 of file TMDCWire.cxx.

: _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::~TMDCWire ( )

virtual

Destructor.

Definition at line 53 of file TMDCWire.cxx.

                    {
}

Member Function Documentation

adjacent()

bool TMDCWire::adjacent ( const TMDCWire &  w ) const

inline

returns true if a given wire is adjacent.

Definition at line 403 of file TMDCWire.h.

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

axial()

bool TMDCWire::axial ( void  ) const

inline

returns true if this wire is in an axial layer.

Definition at line 348 of file TMDCWire.h.

                          {
    return _layer->axial();
}

Referenced by name(), TMDCWireHit::TMDCWireHit(), and TMDC::update().

axialStereoLayerId()

unsigned TMDCWire::axialStereoLayerId ( void  ) const

inline

returns id of axial or stereo id.

Definition at line 360 of file TMDCWire.h.

                                       {
    return _layer->axialStereoLayerId();
}

Referenced by TConformalFinder0::findSegments().

backwardPosition() [1/2]

double * TMDCWire::backwardPosition ( double  p[3] ) const

inline

Definition at line 318 of file TMDCWire.h.

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

backwardPosition() [2/2]

const HepPoint3D & TMDCWire::backwardPosition ( void  ) const

inline

returns position in backward endplate.

Definition at line 306 of file TMDCWire.h.

                                     {
    return _backwardPosition;
}

Referenced by TTrack::approach(), T3DLine::approach(), TRunge::approach(), TBuilder::initialLine(), TBuilder::initialLineOld(), TTrack::stereoHitForCurl(), and TTrack::szPosition().

cellSize()

float TMDCWire::cellSize ( void  ) const

inline

returns cell size in phi.

Definition at line 380 of file TMDCWire.h.

                             {
    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.;
    }
}

clear()

void TMDCWire::clear ( void  )

inline

clears information.

Definition at line 273 of file TMDCWire.h.

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

consective()

bool TMDCWire::consective ( const TMDCWire &  w ) const

inline

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

Definition at line 395 of file TMDCWire.h.

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

direction()

const HepVector3D & TMDCWire::direction ( void  ) const

inline

returns direction vector of the wire.

Definition at line 342 of file TMDCWire.h.

                              {
    return _direction;
}

Referenced by TTrack::approach(), T3DLine::approach(), TRunge::approach(), TCosmicFitter::fit(), TTrack::stereoHitForCurl(), and TTrack::szPosition().

dump()

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.

                                                                 {
    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 + "    ");
    }
}

Referenced by dump().

forwardPosition()

const HepPoint3D & TMDCWire::forwardPosition ( void  ) const

inline

returns position in forward endplate.

Definition at line 300 of file TMDCWire.h.

                                    {
    return _forwardPosition;
}

Referenced by TRunge::approach(), TBuilder::initialLine(), TBuilder::initialLineOld(), TTrack::stereoHitForCurl(), and TTrack::szPosition().

geocdc()

const MdcGeoWire * TMDCWire::geocdc ( void  ) const

inline

returns a pointer to GEOMDC_WIR.

Definition at line 237 of file TMDCWire.h.

                           {
    return _geo;
}

hit() [1/3]

const TMDCWireHit *const TMDCWire::hit ( const TMDCWireHit * const  h )

inline

sets a pointer to TMDCWireHit.

Definition at line 281 of file TMDCWire.h.

                                         {
    return _hit = h;
}

hit() [2/3]

const TMDCWireHitMC *const TMDCWire::hit ( TMDCWireHitMC * const  a )

inline

appends a pointer to TMDCWireHitMC.

Definition at line 293 of file TMDCWire.h.

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

hit() [3/3]

const TMDCWireHit *const TMDCWire::hit ( void  ) const

inline

returns a pointer to a TMDCWireHit.

Definition at line 287 of file TMDCWire.h.

                        {
    return _hit;
}

Referenced by TMDC::classification(), TMDCWireHit::TMDCWireHit(), TMDC::update(), and TMDC::updateMC().

id()

unsigned TMDCWire::id ( void  ) const

inline

returns id.

Definition at line 207 of file TMDCWire.h.

                       {
    return _id;
}

Referenced by TBuilderCurl::buildStereo(), TMDC::driftDistance(), TRungeFitter::fit(), TMDC::updateMC(), and wirePosition().

innerPart()

bool TMDCWire::innerPart ( void  ) const

inline

returns true if this wire is in the inner part.

Definition at line 366 of file TMDCWire.h.

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

layer()

const TMDCLayer *const TMDCWire::layer ( void  ) const

inline

returns a pointer to a layer.

Definition at line 243 of file TMDCWire.h.

                          {
    return _layer;
}

Referenced by TBuilder0::buildStereo(), Edges(), localIdDifference(), TTrackManager::removeHitsAcrossOverIp(), and Width().

layerId()

unsigned TMDCWire::layerId ( void  ) const

inline

returns layer id.

Definition at line 219 of file TMDCWire.h.

                            {
    return _layerId;
}

Referenced by TLine0::appendByszdistance(), TMLine::appendByszdistance(), TBuilderCosmic::buildStereo(), TMDCTsf::createTsf(), TPerfectFinder::doit(), TMDCTsf::dump(), TMDCWireHit::dump(), findIsolatedCloseHits(), TCosmicFitter::fit(), TRungeFitter::fit(), TLine0::fit2(), TMLine::fit2(), TLine0::fit2p(), TMLine::fit2p(), TLine0::fit2s(), TMLine::fit2s(), TLine0::fit2sp(), TMLine::fit2sp(), TTrackManager::makeTds(), TTrackManager::maskNormal(), TTrackManager::merge(), TTrack::movePivot(), neighbor(), offsetBorder(), TLine0::refine(), TMLine::refine(), TLine0::removeChits(), TMLine::removeChits(), SameLayer(), and TCircle::weight().

localId()

unsigned TMDCWire::localId ( void  ) const

inline

returns local id in a wire layer.

Definition at line 213 of file TMDCWire.h.

                            {
    return _localId;
}

Referenced by TLine0::appendByszdistance(), TMLine::appendByszdistance(), TBuilder0::buildStereo(), TBuilderCosmic::buildStereo(), TMDCTsf::createTsf(), TMDCTsf::dump(), TMDCWireHit::dump(), Edges(), findIsolatedCloseHits(), TCosmicFitter::fit(), TRungeFitter::fit(), TLine0::fit2(), TMLine::fit2(), TLine0::fit2p(), TMLine::fit2p(), localIdDifference(), TTrackManager::makeTds(), neighbor(), TLine0::refine(), TMLine::refine(), TTrackManager::removeHitsAcrossOverIp(), TCircle::weight(), and Width().

localIdDifference()

int TMDCWire::localIdDifference

(

const TMDCWire &

a

)

const

returns local id difference.

Definition at line 577 of file TMDCWire.cxx.

                                                    {
 
#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;
    }
}

localIdForMinus()

int TMDCWire::localIdForMinus

(

void

)

const

returns localId but if 0, return maxLocalId + 1.

Definition at line 355 of file TMDCWire.cxx.

                                    {
/*
  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;
}

localIdForPlus()

int TMDCWire::localIdForPlus

(

void

)

const

returns localId but if maxLocalId, return -1.

Definition at line 166 of file TMDCWire.cxx.

                                   {
/*
  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;
}

localLayerId()

unsigned TMDCWire::localLayerId ( void  ) const

inline

returns local layer id in a super layer.

Definition at line 231 of file TMDCWire.h.

                                 {
    return _layer->localLayerId();
}

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

mainPart()

bool TMDCWire::mainPart ( void  ) const

inline

returns true if this wire is in the main part.

Definition at line 373 of file TMDCWire.h.

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

name()

std::string TMDCWire::name ( void  ) const

inline

returns name.

Definition at line 412 of file TMDCWire.h.

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

Referenced by TTrackBase::append(), TTrack::approach(), TTrack::assign(), TBuilder0::buildRphi(), TBuilder0::buildStereo(), TTrackManager::divide(), Dump(), TMDCWireHit::dump(), TTrackManager::salvage(), and TTrackManager::salvageAssociateHits().

neighbor()

const TMDCWire *const TMDCWire::neighbor

(

unsigned

i

)

const

returns a pointer to a neighbor wire.

Definition at line 95 of file TMDCWire.cxx.

                                   {
    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;
}

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

state() [1/2]

unsigned TMDCWire::state ( unsigned  newState )

inline

sets state.

Definition at line 267 of file TMDCWire.h.

                          {
    return _state = a;
}

state() [2/2]

unsigned TMDCWire::state ( void  ) const

inline

returns state.

Definition at line 261 of file TMDCWire.h.

                          {
    return _state;
}

stereo()

bool TMDCWire::stereo ( void  ) const

inline

returns true if this wire is in a stereo layer.

Definition at line 354 of file TMDCWire.h.

                           {
    return _layer->stereo();
}

Referenced by TCosmicFitter::fit(), and TTrackManager::makeTds().

superLayerId()

unsigned TMDCWire::superLayerId ( void  ) const

inline

returns super layer id.

Definition at line 225 of file TMDCWire.h.

                                 {
    return _layer->superLayerId();
}

Referenced by TLine0::appendByszdistance(), TMLine::appendByszdistance(), TMDC::classification(), dump(), TLine0::fit2(), TMLine::fit2(), TLine0::fit2p(), TMLine::fit2p(), localIdDifference(), neighbor(), SameSuperLayer(), and Width().

tension()

const double & TMDCWire::tension ( void  ) const

inline

Definition at line 312 of file TMDCWire.h.

                            {
        return _tension;
}

tsf() [1/2]

TMDCTsf * TMDCWire::tsf ( TMDCTsf *  t )

inline

sets a pointer of TMDCTsf.

Definition at line 255 of file TMDCWire.h.

                         {
    return _tsf = t;
}

tsf() [2/2]

const TMDCTsf *const TMDCWire::tsf ( void  ) const

inline

returns a pointer to a TMDCTsf.

Definition at line 249 of file TMDCWire.h.

                        {
    return _tsf;
}

wirePosition()

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.

                                        {
    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;
}

Referenced by TTrack::approach(), T3DLine::approach(), TRunge::approach(), TRunge::SetFlightLength(), and wirePosition().

xyPosition() [1/3]

double * TMDCWire::xyPosition ( double  p[3] ) const

inline

Definition at line 333 of file TMDCWire.h.

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

xyPosition() [2/3]

const HepPoint3D TMDCWire::xyPosition

(

double

z

)

const

Definition at line 76 of file TMDCWire.cxx.

                                   {
    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;
}

xyPosition() [3/3]

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.

                               {
    return _xyPosition;
}

Referenced by TTrack::approach(), T3DLine::approach(), TRunge::approach(), TTrack::approach2D(), calVirtualCircle(), TRunge::DisToWire(), TPerfectFinder::doit(), TTrack::fit2D(), TTrackManager::salvageAssociateHits(), and TMLink::xyPosition().

---

The documentation for this class was generated from the following files:

• TMDCWire.h
• TMDCWire.cxx