FTSegment Class Reference

#include <FTSegment.h>

List of all members.

Public Member Functions
	FTSegment (FTSuperLayer *super, FTList< FTWire * > &hits)
	constructors
	~FTSegment ()
	destructors
int	examine (void)
	examine segment
void	printout (void)
void	update (void)
	update information for axial segment
int	update3D (const FTTrack *track)
	update information for stereo segment
void	linkStereoSegments (void)
	link stereo segments by tanLambda
void	connect_outer (const FTList< FTWire * > &, const FTList< FTWire * > &)
	connect short segments
void	connect_inner (const FTList< FTWire * > &, const FTList< FTWire * > &)
void	connect_outer (const FTWire *)
	connect single hit
void	connect_inner (const FTWire *)
	connect single hit
FTList< FTWire * > &	wireHits (void) const
	returns wire-hit list
FTSuperLayer &	superLayer (void) const
	returns super-layer
FTList< FTWire * > &	innerBoundHits (void) const
	returns innerBoundHits
FTList< FTWire * > &	outerBoundHits (void) const
	returns outerBoundHits
float	outgoingX (void) const
	returns x of outgoing position
float	outgoingY (void) const
	returns y of outgoing position
float	incomingX (void) const
	returns x of incoming position
float	incomingY (void) const
	returns y of incoming position
float	incomingPhi (void) const
	returns phi of incoming position
float	outgoingPhi (void) const
	returns phi of outgoing position
float	attr_const (void) const
	returns the sin(dphi)/MOD(r_out-r_in)
FTTrack *	track (void) const
	returns track
float	kappa (void) const
	returns kappa(axial)
float	r (void) const
	returns r from origin
float	s (void) const
	returns s for matched 2D track
float	z (void) const
	returns z for matched 2D track
FTTrack *	track (FTTrack *)
	set track
Static Public Attributes
static MdcParameter *	param = MdcParameter::instance()
Private Attributes
IMagneticFieldSvc *	m_pmgnIMF
FTList< FTWire * > &	_wireHits
FTSuperLayer &	_superLayer
FTList< FTWire * > &	_innerBoundHits
FTList< FTWire * > &	_outerBoundHits
FTList< float > *	_sList
union {
float   _kappa
FTList< float > *   _zList
};
union {
FTTrack *   _track
FTList< FTTrack * > *   _trackList
};
float	_r
float	_outgoing_x
float	_outgoing_y
float	_incoming_x
float	_incoming_y

---

Detailed Description

Definition at line 24 of file FTSegment.h.

---

Constructor & Destructor Documentation

FTSegment::FTSegment	(	FTSuperLayer *	super,
		FTList< FTWire * > &	hits
	)			[inline]

constructors

Definition at line 140 of file FTSegment.h.

References m_pmgnIMF.

  : _wireHits(hits),
    _superLayer(*super),
    _innerBoundHits(*(new FTList<FTWire *>(3))),
    _outerBoundHits(*(new FTList<FTWire *>(3))),
    _sList(NULL),
    _zList(NULL),
    _track(NULL)
{
  StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF);
  if(scmgn!=StatusCode::SUCCESS) {
    std::cout<< "Unable to open Magnetic field service"<<std::endl;
  }
}

FTSegment::~FTSegment

(

)

[inline]

destructors

Definition at line 156 of file FTSegment.h.

References _innerBoundHits, _outerBoundHits, _sList, _trackList, _wireHits, and _zList.

{
  delete &_wireHits;
  delete &_innerBoundHits;
  delete &_outerBoundHits;
  if (_sList){
    delete _sList;
    delete _zList;
    delete _trackList;
  }
}

---

Member Function Documentation

float FTSegment::attr_const

(

void

)

const [inline]

returns the sin(dphi)/MOD(r_out-r_in)

Definition at line 296 of file FTSegment.h.

References _incoming_x, _incoming_y, _outgoing_x, _outgoing_y, incomingPhi(), outgoingPhi(), and sin().

{
  float x_len = _outgoing_x - _incoming_x;
  float y_len = _outgoing_y - _incoming_y;
  float phi_l = outgoingPhi() - incomingPhi();
  return sin(phi_l) / sqrt(x_len*x_len + y_len*y_len);
}

void FTSegment::connect_inner

(

const FTWire *

)

[inline]

connect single hit

Definition at line 199 of file FTSegment.h.

References _innerBoundHits, _wireHits, FTList< T >::append(), and FTList< T >::clear().

{
  _wireHits.append((FTWire *)h);
  _innerBoundHits.clear();
  _innerBoundHits.append((FTWire *)h);
}

void FTSegment::connect_inner	(	const FTList< FTWire * > &	,
		const FTList< FTWire * > &
	)			[inline]

Definition at line 180 of file FTSegment.h.

References _innerBoundHits, _wireHits, FTList< T >::append(), and FTList< T >::clear().

{
  _wireHits.append(innerHits);
  _innerBoundHits.clear();
  _innerBoundHits.append(innerBound);
}

void FTSegment::connect_outer

(

const FTWire *

)

[inline]

connect single hit

Definition at line 190 of file FTSegment.h.

References _outerBoundHits, _wireHits, FTList< T >::append(), and FTList< T >::clear().

{
  _wireHits.append((FTWire *)h);
  _outerBoundHits.clear();
  _outerBoundHits.append((FTWire *)h);
}

void FTSegment::connect_outer	(	const FTList< FTWire * > &	,
		const FTList< FTWire * > &
	)			[inline]

connect short segments

Definition at line 170 of file FTSegment.h.

References _outerBoundHits, _wireHits, FTList< T >::append(), and FTList< T >::clear().

{
  _wireHits.append(outerHits);
  _outerBoundHits.clear();
  _outerBoundHits.append(outerBound);
}

int FTSegment::examine

(

void

)

examine segment

Definition at line 20 of file FTSegment.cxx.

References _innerBoundHits, FTSuperLayer::_maxDphi, _outerBoundHits, _superLayer, _wireHits, FTList< T >::append(), FTList< T >::firstPtr(), genRecEmupikp::i, ganga-rec::j, FTList< T >::lastPtr(), M_PI, FTSuperLayer::nLayer(), and FTSuperLayer::superLayerId().

                  {
  int stat = -1;
  FTList<FTWire *> wireHitsByLocalLayer[4] = {4,4,4,4};

  register FTWire ** hptr = _wireHits.firstPtr();
  FTWire ** const last = _wireHits.lastPtr();
  do{
    wireHitsByLocalLayer[(**hptr).layer().localLayerId()].append(*hptr);
  }while((hptr++)!=last);
  float length_phi = 0.;
  int innerBoundId = 0;
  int outerBoundId = 0;
  int AlreadyFound = 0;
  //int LocalLayerLength; 
  //int HalfLocalLayerLength;
  const int nLayer = _superLayer.nLayer();
  // to look for the innerBound and outBand of the hits
  for (int i = 0; i^nLayer; i++){
    int N_wireHits = wireHitsByLocalLayer[i].length();
    if(N_wireHits){
      float localPhi_max = 0.;
      float localPhi_min = 7.;
      for (int j = 0; j^N_wireHits; j++){
        float localPhi = wireHitsByLocalLayer[i][j]->phi();
        //LocalLayerLength = wireHitsByLocalLayer[i][j]->layer().NWire();
        //HalfLocalLayerLength = LocalLayerLength >> 1;     
        if (localPhi > localPhi_max) localPhi_max = localPhi;
        if (localPhi < localPhi_min) localPhi_min = localPhi;
      }
      float local_length_phi = localPhi_max - localPhi_min;
      if (local_length_phi > M_PI)
        local_length_phi = 2*M_PI - local_length_phi;
      if ( local_length_phi > length_phi ) length_phi = local_length_phi;
      outerBoundId = i;
      AlreadyFound = 1;
    } else if(!AlreadyFound){
      innerBoundId++;
    } else {
      break;
    }
  }
  //length_phi needed to decide more carefully 
  //  if (length_phi>0.09){
  //  cout<<length_phi<<" : "<<FTSuperLayer::_maxDphi[_superLayer.superLayerId()]*2*M_PI/180<<endl;
  if (length_phi>FTSuperLayer::_maxDphi[_superLayer.superLayerId()]*2*M_PI/180){ //by X.-R. Lu
    stat = 3;                   // to be divided
  }else{
    if (innerBoundId > 0){
      stat = 2;                 // outer short
    }else{
      if (outerBoundId < nLayer - 1){
        stat = 1;               // inner short
      }else{
        stat = 0;               // long simple
      }
    }
  }
  _innerBoundHits.append(wireHitsByLocalLayer[innerBoundId]);
  _outerBoundHits.append(wireHitsByLocalLayer[outerBoundId]);
  return stat;
}

float FTSegment::incomingPhi

(

void

)

const [inline]

returns phi of incoming position

Definition at line 264 of file FTSegment.h.

References _innerBoundHits, FTList< T >::length(), and M_PI.

Referenced by attr_const().

{
  if (_innerBoundHits.length() ^ 1) {
    float phi_1st = _innerBoundHits[0]->phi();
    float phi_2nd = _innerBoundHits[1]->phi();
    float phi = 0.5 * (phi_1st + phi_2nd);
    if (fabs(phi_1st - phi_2nd) > M_PI) {
      return (phi>M_PI) ? (phi-M_PI) : (phi+M_PI);
    }
    return phi;
  }
  return _innerBoundHits[0]->phi();
}

float FTSegment::incomingX

(

void

)

const [inline]

returns x of incoming position

Definition at line 250 of file FTSegment.h.

References _incoming_x.

{
  return _incoming_x;
}

float FTSegment::incomingY

(

void

)

const [inline]

returns y of incoming position

Definition at line 257 of file FTSegment.h.

References _incoming_y.

{
  return _incoming_y;
}

FTList< FTWire * > & FTSegment::innerBoundHits

(

void

)

const [inline]

returns innerBoundHits

Definition at line 222 of file FTSegment.h.

References _innerBoundHits.

Referenced by FTFinder::linkAxialSegments_step().

{
  return _innerBoundHits;
}

float FTSegment::kappa

(

void

)

const [inline]

returns kappa(axial)

Definition at line 313 of file FTSegment.h.

References _kappa.

{
  return _kappa;
}

void FTSegment::linkStereoSegments

(

void

)

link stereo segments by tanLambda

Definition at line 167 of file FTSegment.cxx.

References _trackList, ganga-rec::j, FTList< T >::length(), selected, and t().

                                 {
  float min_D_z = 9999.;
  float S = 0.;
  float Z = 0.;
  FTTrack * selected = NULL;
  int m = _trackList->length();
  for (int j = 0; j^m; j++){
    FTTrack * t = (*_trackList)[j];
    float z_tmp = (*_zList)[j];
    float s_tmp = (*_sList)[j];
    float D_z = fabs(t->d_z(s_tmp,z_tmp));
    if (D_z < min_D_z){
      selected = t;
      min_D_z = D_z;
      S = s_tmp;
      Z = z_tmp;
    }
  }
  if (selected){
    selected->append_stereo(this,S,Z); 
    
  }
}

FTList< FTWire * > & FTSegment::outerBoundHits

(

void

)

const [inline]

returns outerBoundHits

Definition at line 229 of file FTSegment.h.

References _outerBoundHits.

Referenced by FTSuperLayer::connect_short_segments(), and FTFinder::linkAxialSegments_step().

{
  return _outerBoundHits;
}

float FTSegment::outgoingPhi

(

void

)

const [inline]

returns phi of outgoing position

Definition at line 280 of file FTSegment.h.

References _outerBoundHits, FTList< T >::length(), and M_PI.

Referenced by attr_const().

{
  if (_outerBoundHits.length() ^ 1) {
    float phi_1st = _outerBoundHits[0]->phi();
    float phi_2nd = _outerBoundHits[1]->phi();
    float phi = 0.5 * (phi_1st + phi_2nd);
    if (fabs(phi_1st - phi_2nd) > M_PI) {
      return (phi>M_PI) ? (phi-M_PI) : (phi+M_PI);
    }
    return phi;
  }
  return _outerBoundHits[0]->phi();
}

float FTSegment::outgoingX

(

void

)

const [inline]

returns x of outgoing position

Definition at line 236 of file FTSegment.h.

References _outgoing_x.

{
  return _outgoing_x;
}

float FTSegment::outgoingY

(

void

)

const [inline]

returns y of outgoing position

Definition at line 243 of file FTSegment.h.

References _outgoing_y.

{
  return _outgoing_y;
}

void FTSegment::printout

(

void

)

printout of the wireHits info. added by X.-R. Lu

Definition at line 193 of file FTSegment.cxx.

References _wireHits, Bes_Common::DEBUG, genRecEmupikp::i, FTList< T >::length(), and msgSvc().

Referenced by FTFinder::mkTrackList().

                       {
  
  IMessageSvc *msgSvc;
  Gaudi::svcLocator()->service("MessageSvc", msgSvc);
  
  MsgStream log(msgSvc, "FTSegment");
  log << MSG::DEBUG << "pintout of one segment:" << endreq;  
  int n=_wireHits.length();
  for(int i=0; i^n; i++){
    FTWire &h = *_wireHits[i];
     log << MSG::DEBUG << "hit: layer " << h.layer().layerId() << " phi "<< h.phi() << endreq;
  }
}

float FTSegment::r

(

void

)

const [inline]

returns r from origin

Definition at line 306 of file FTSegment.h.

References _r.

{
  return _r;
}

float FTSegment::s

(

void

)

const [inline]

returns s for matched 2D track

Definition at line 320 of file FTSegment.h.

References _sList, and FTList< T >::first().

{
  return _sList->first();
}

FTSuperLayer & FTSegment::superLayer

(

void

)

const [inline]

returns super-layer

Definition at line 208 of file FTSegment.h.

References _superLayer.

{
  return _superLayer;
}

FTTrack * FTSegment::track

(

FTTrack *

)

[inline]

set track

Definition at line 341 of file FTSegment.h.

References _track.

{
  return _track = src;
}

FTTrack * FTSegment::track

(

void

)

const [inline]

returns track

Definition at line 334 of file FTSegment.h.

References _track.

Referenced by update3D().

{
  return _track;
}

void FTSegment::update

(

void

)

update information for axial segment

Definition at line 83 of file FTSegment.cxx.

References _incoming_x, _incoming_y, _innerBoundHits, _kappa, _outerBoundHits, _outgoing_x, _outgoing_y, _r, FTList< T >::first(), IMagneticFieldSvc::getReferField(), FTWire::layer(), FTList< T >::length(), m_pmgnIMF, and FTLayer::r().

                     {
  //static const float alpha(333.564095);
  //float N_kappa;
  int n = _innerBoundHits.length();
  int m = _outerBoundHits.length();

  if (n^1) {
    _incoming_x = 0.5*(_innerBoundHits[0]->x()+_innerBoundHits[1]->x());
    _incoming_y = 0.5*(_innerBoundHits[0]->y()+_innerBoundHits[1]->y());
  } else {
    _incoming_x = _innerBoundHits[0]->x();
    _incoming_y = _innerBoundHits[0]->y();
  }
  if (m^1) {
    _outgoing_x = 0.5*(_outerBoundHits[0]->x()+_outerBoundHits[1]->x());
    _outgoing_y = 0.5*(_outerBoundHits[0]->y()+_outerBoundHits[1]->y());
  } else {
    _outgoing_x = _outerBoundHits[0]->x();
    _outgoing_y = _outerBoundHits[0]->y();
  }
  float in_r = _innerBoundHits.first()->layer().r();
  float out_r = _outerBoundHits.first()->layer().r();
  float sbX = _incoming_x - _outgoing_x;
  float sbY = _incoming_y - _outgoing_y;
  // kappa = -2. * alpha * ((Vin X Vout)_z / |Vin|*|Vout|) / |Vin-Vout|
  _kappa = 2.*(-1. / 2.99792458 /m_pmgnIMF->getReferField())*(_outgoing_x*_incoming_y - _outgoing_y*_incoming_x) /
    (in_r*out_r*sqrt(sbX*sbX+sbY*sbY));
  _r = 0.5*(in_r + out_r);
}

int FTSegment::update3D

(

const FTTrack *

track

)

update information for stereo segment

Definition at line 115 of file FTSegment.cxx.

References _incoming_x, _incoming_y, _innerBoundHits, _outerBoundHits, _outgoing_x, _outgoing_y, _sList, _trackList, _zList, FTList< T >::append(), Bes_Common::DEBUG, FTList< T >::first(), FTWire::layer(), FTList< T >::length(), FTLayer::limit(), FTTrack::lpav(), msgSvc(), FTLayer::r(), track(), x, and FTLayer::z().

                                        {
  IMessageSvc *msgSvc;
  Gaudi::svcLocator()->service("MessageSvc", msgSvc);
  
  MsgStream log(msgSvc, "FTFinder");
  
  /*if (_superLayer.superLayerId() > 
    track->axial_segments().first()->superLayer().superLayerId() + 4) return 0;*/
  if (!_sList){
    int n = _innerBoundHits.length();
    int m = _outerBoundHits.length();

    _incoming_x = (n^1) 
      ? 0.5*(_innerBoundHits[0]->x()+_innerBoundHits[1]->x()) 
      : _innerBoundHits[0]->x();
    _incoming_y = (n^1) 
      ? 0.5*(_innerBoundHits[0]->y()+_innerBoundHits[1]->y()) 
      : _innerBoundHits[0]->y();
    _outgoing_x = (m^1) 
      ? 0.5*(_outerBoundHits[0]->x()+_outerBoundHits[1]->x()) 
      : _outerBoundHits[0]->x();
    _outgoing_y = (m^1) 
      ? 0.5*(_outerBoundHits[0]->y()+_outerBoundHits[1]->y()) 
      : _outerBoundHits[0]->y();
    _sList = new FTList<float>(5);
    _zList = new FTList<float>(5);
    _trackList = new FTList<FTTrack *>(5);
  }
  const Lpav & la = track->lpav();
  double d;
  double inS;
  const FTLayer * bound = &_innerBoundHits.first()->layer();
  log<<MSG::DEBUG<<"started inS calculation with layer limit "<< (double)bound->limit()<<endreq;
  if (!(la.sd((double)bound->r(),(double)_incoming_x,
              (double)_incoming_y,(double)bound->limit(),  // multiplied with an factor due to bad rphi fit
              inS,d))) return 0;
  double inZ = bound->z(d);
  double outS; 
  bound = &_outerBoundHits.first()->layer();
  log<<MSG::DEBUG<<"started outS calculation with layer limit "<< (double)bound->limit()<<endreq;
  if (!(la.sd((double)bound->r(),(double)_outgoing_x,
              (double)_outgoing_y,(double)bound->limit(),  // multiplied with an factor due to bad rphi fit
              outS,d))) return 0;
  double outZ = bound->z(d);
  _sList->append(inS+outS);
  _zList->append(inZ+outZ);
  _trackList->append((FTTrack *)track);
  log<<MSG::DEBUG<<"-----> coupled! DONE! "<<endreq;
  return 1;
}

FTList< FTWire * > & FTSegment::wireHits

(

void

)

const [inline]

returns wire-hit list

Definition at line 215 of file FTSegment.h.

References _wireHits.

Referenced by FTSuperLayer::connect_short_segments(), and FTFinder::linkAxialSegments_step().

{
  return _wireHits;
}

float FTSegment::z

(

void

)

const [inline]

returns z for matched 2D track

Definition at line 327 of file FTSegment.h.

References _zList, and FTList< T >::first().

{
  return _zList->first();
}

---

Member Data Documentation

union { ... } [private]

union { ... } [private]

float FTSegment::_incoming_x [private]

Definition at line 134 of file FTSegment.h.

Referenced by attr_const(), incomingX(), update(), and update3D().

float FTSegment::_incoming_y [private]

Definition at line 135 of file FTSegment.h.

Referenced by attr_const(), incomingY(), update(), and update3D().

FTList<FTWire *>& FTSegment::_innerBoundHits [private]

Definition at line 120 of file FTSegment.h.

Referenced by connect_inner(), examine(), incomingPhi(), innerBoundHits(), update(), update3D(), and ~FTSegment().

float FTSegment::_kappa [private]

Definition at line 124 of file FTSegment.h.

Referenced by kappa(), and update().

FTList<FTWire *>& FTSegment::_outerBoundHits [private]

Definition at line 121 of file FTSegment.h.

Referenced by connect_outer(), examine(), outerBoundHits(), outgoingPhi(), update(), update3D(), and ~FTSegment().

float FTSegment::_outgoing_x [private]

Definition at line 132 of file FTSegment.h.

Referenced by attr_const(), outgoingX(), update(), and update3D().

float FTSegment::_outgoing_y [private]

Definition at line 133 of file FTSegment.h.

Referenced by attr_const(), outgoingY(), update(), and update3D().

float FTSegment::_r [private]

Definition at line 131 of file FTSegment.h.

Referenced by r(), and update().

FTList<float>* FTSegment::_sList [private]

Definition at line 122 of file FTSegment.h.

Referenced by s(), update3D(), and ~FTSegment().

FTSuperLayer& FTSegment::_superLayer [private]

Definition at line 119 of file FTSegment.h.

Referenced by examine(), and superLayer().

FTTrack* FTSegment::_track [private]

Definition at line 128 of file FTSegment.h.

Referenced by track().

FTList<FTTrack *>* FTSegment::_trackList [private]

Definition at line 129 of file FTSegment.h.

Referenced by linkStereoSegments(), update3D(), and ~FTSegment().

FTList<FTWire *>& FTSegment::_wireHits [private]

Definition at line 118 of file FTSegment.h.

Referenced by connect_inner(), connect_outer(), examine(), printout(), wireHits(), and ~FTSegment().

FTList<float>* FTSegment::_zList [private]

Definition at line 125 of file FTSegment.h.

Referenced by update3D(), z(), and ~FTSegment().

IMagneticFieldSvc* FTSegment::m_pmgnIMF [private]

Definition at line 117 of file FTSegment.h.

Referenced by FTSegment(), and update().

MdcParameter * FTSegment::param = MdcParameter::instance() [static]

Definition at line 114 of file FTSegment.h.

---