TBuilder Class Reference

A class to build a track. More...

#include <TBuilder.h>

Public Member Functions
	TBuilder (const std::string &name, float maxSigma, float maxSigmaStereo, float salvageLevel, float szSegmentDistance, float szLinkDistance, unsigned fittingFlag)
	Constructor with salvage level.
virtual	~TBuilder ()
	Destructor.
const 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.
TTrack *	buildRphi (const AList< TMLink > &) const
	builds a r/phi track from TMLinks or from Segments.
TTrack *	buildRphi (const AList< TSegment > &) const
TTrack *	buildStereoNew (const TTrack &t, AList< TSegment > &goodSegments, AList< TSegment > &badSegments) const
TTrack *	buildStereo (const TTrack &t, AList< TSegment > &) const
TTrack *	buildStereo (TTrack &t, TMLine &l, const AList< TMLink > &links) const
TTrack *	buildStereoOld (TTrack &t, const AList< TMLink > &allLinks, const AList< TMLink > &isolatedLinks) const
TTrack *	build (TTrack &t, const TMLine &l) const
void	salvage (TTrack &t, AList< TMLink > &hits) const
	salvages hits.
void	salvage (TTrack &t, AList< TSegment > &segments) const
TMLine *	initialLine (const TTrack &, AList< TSegment > &) const
	makes a line.
TMLine *	initialLineOld (const TTrack &, AList< TSegment > &) const
TMLine *	initialLine1 (const TTrack &, const AList< TSegment > &, const AList< TMLink > &) const
TMLine *	initialLine2 (const TTrack &, const AList< TMLink > &) const
AList< TMLine >	searchInitialLines (unsigned nSuperLayers) const
AList< TMLine >	searchLines6 (void) const
AList< TMLine >	searchLines5 (void) const
AList< TMLine >	searchLines4 (void) const
AList< TMLine >	searchLines3 (void) const
AList< TMLine >	searchLines2 (void) const
AList< TMLine >	searchLines1 (void) const
TMLine	searchLine (const TMLine &initialLine) const
AList< TSegment >	selectStereoSegment (const TMLine &line, const AList< TSegment > &list, const AList< TMLink > &szList) const
void	removeFarSegment (const TMLine &, AList< TSegment > &, AList< TMLink > &) const
	TBuilder (const std::string &name, float maxSigma, float maxSigmaStereo, float salvageLevel, float szSegmentDistance, float szLinkDistance, unsigned fittingFlag)
	Constructor with salvage level.
virtual	~TBuilder ()
	Destructor.
const 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.
TTrack *	buildRphi (const AList< TMLink > &) const
	builds a r/phi track from TMLinks or from Segments.
TTrack *	buildRphi (const AList< TSegment > &) const
TTrack *	buildStereoNew (const TTrack &t, AList< TSegment > &goodSegments, AList< TSegment > &badSegments) const
TTrack *	buildStereo (const TTrack &t, AList< TSegment > &) const
TTrack *	buildStereo (TTrack &t, TMLine &l, const AList< TMLink > &links) const
TTrack *	buildStereoOld (TTrack &t, const AList< TMLink > &allLinks, const AList< TMLink > &isolatedLinks) const
TTrack *	build (TTrack &t, const TMLine &l) const
void	salvage (TTrack &t, AList< TMLink > &hits) const
	salvages hits.
void	salvage (TTrack &t, AList< TSegment > &segments) const
TMLine *	initialLine (const TTrack &, AList< TSegment > &) const
	makes a line.
TMLine *	initialLineOld (const TTrack &, AList< TSegment > &) const
TMLine *	initialLine1 (const TTrack &, const AList< TSegment > &, const AList< TMLink > &) const
TMLine *	initialLine2 (const TTrack &, const AList< TMLink > &) const
AList< TMLine >	searchInitialLines (unsigned nSuperLayers) const
AList< TMLine >	searchLines6 (void) const
AList< TMLine >	searchLines5 (void) const
AList< TMLine >	searchLines4 (void) const
AList< TMLine >	searchLines3 (void) const
AList< TMLine >	searchLines2 (void) const
AList< TMLine >	searchLines1 (void) const
TMLine	searchLine (const TMLine &initialLine) const
AList< TSegment >	selectStereoSegment (const TMLine &line, const AList< TSegment > &list, const AList< TMLink > &szList) const
void	removeFarSegment (const TMLine &, AList< TSegment > &, AList< TMLink > &) const
	TBuilder (const std::string &name, float maxSigma, float maxSigmaStereo, float salvageLevel, float szSegmentDistance, float szLinkDistance, unsigned fittingFlag)
	Constructor with salvage level.
virtual	~TBuilder ()
	Destructor.
const 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.
TTrack *	buildRphi (const AList< TMLink > &) const
	builds a r/phi track from TMLinks or from Segments.
TTrack *	buildRphi (const AList< TSegment > &) const
TTrack *	buildStereoNew (const TTrack &t, AList< TSegment > &goodSegments, AList< TSegment > &badSegments) const
TTrack *	buildStereo (const TTrack &t, AList< TSegment > &) const
TTrack *	buildStereo (TTrack &t, TMLine &l, const AList< TMLink > &links) const
TTrack *	buildStereoOld (TTrack &t, const AList< TMLink > &allLinks, const AList< TMLink > &isolatedLinks) const
TTrack *	build (TTrack &t, const TMLine &l) const
void	salvage (TTrack &t, AList< TMLink > &hits) const
	salvages hits.
void	salvage (TTrack &t, AList< TSegment > &segments) const
TMLine *	initialLine (const TTrack &, AList< TSegment > &) const
	makes a line.
TMLine *	initialLineOld (const TTrack &, AList< TSegment > &) const
TMLine *	initialLine1 (const TTrack &, const AList< TSegment > &, const AList< TMLink > &) const
TMLine *	initialLine2 (const TTrack &, const AList< TMLink > &) const
AList< TMLine >	searchInitialLines (unsigned nSuperLayers) const
AList< TMLine >	searchLines6 (void) const
AList< TMLine >	searchLines5 (void) const
AList< TMLine >	searchLines4 (void) const
AList< TMLine >	searchLines3 (void) const
AList< TMLine >	searchLines2 (void) const
AList< TMLine >	searchLines1 (void) const
TMLine	searchLine (const TMLine &initialLine) const
AList< TSegment >	selectStereoSegment (const TMLine &line, const AList< TSegment > &list, const AList< TMLink > &szList) const
void	removeFarSegment (const TMLine &, AList< TSegment > &, AList< TMLink > &) const

Detailed Description

A class to build a track.

Definition at line 34 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TBuilder.h.

Constructor & Destructor Documentation

TBuilder() [1/3]

TBuilder::TBuilder	(	const std::string &	name,
		float	maxSigma,
		float	maxSigmaStereo,
		float	salvageLevel,
		float	szSegmentDistance,
		float	szLinkDistance,
		unsigned	fittingFlag )

Constructor with salvage level.

Definition at line 36 of file TBuilder.cxx.

    : _name( a )
    , _fitter( "TBuilder Fitter" )
    , _maxSigma( maxSigma )
    , _maxSigmaStereo( maxSigmaStereo )
    , _salvageLevel( sqrt( salvageLevel ) )
    ,
    //  _minNCores(3), //20060307
    _minNCores( 2 )
    , // liuqg20070312
    _szSegmentDistance( szSegmentDistance )
    , _szLinkDistance( szLinkDistance ) {
  if ( fittingFlag & 1 ) _fitter.sag( true );
  if ( fittingFlag & 2 ) _fitter.propagation( true );
  if ( fittingFlag & 4 ) _fitter.tof( true );
  if ( fittingFlag & 8 ) _fitter.freeT0( true );
}

~TBuilder() [1/3]

TBuilder::~TBuilder ( )

virtual

Destructor.

Definition at line 56 of file TBuilder.cxx.

{}

TBuilder() [2/3]

TBuilder::TBuilder	(	const std::string &	name,
		float	maxSigma,
		float	maxSigmaStereo,
		float	salvageLevel,
		float	szSegmentDistance,
		float	szLinkDistance,
		unsigned	fittingFlag )

Constructor with salvage level.

~TBuilder() [2/3]

virtual TBuilder::~TBuilder ( )

virtual

Destructor.

TBuilder() [3/3]

TBuilder::TBuilder	(	const std::string &	name,
		float	maxSigma,
		float	maxSigmaStereo,
		float	salvageLevel,
		float	szSegmentDistance,
		float	szLinkDistance,
		unsigned	fittingFlag )

Constructor with salvage level.

~TBuilder() [3/3]

virtual TBuilder::~TBuilder ( )

virtual

Destructor.

Member Function Documentation

build() [1/3]

TTrack * TBuilder::build	(	TTrack &	t,
		const TMLine &	l ) const

Definition at line 1188 of file TBuilder.cxx.

                                                          {
  AList<TMLink> links = l.links();
  AList<TMLink> toRemove;
  unsigned      n = links.length();
  for ( unsigned i = 1; i < n; i++ )
  {
    if ( links[i - 1]->link() == links[i]->link() )
    {
      toRemove.append( links[i] );
      continue;
    }
    if ( i < 2 ) continue;
    if ( links[i - 2]->link() == links[i]->link() ) toRemove.append( links[i] );
  }
  links.remove( toRemove );
 
  //...Pick up links...
  n = links.length();
  for ( unsigned i = 0; i < n; i++ ) t.append( *links[i]->link() );
 
  Vector a( 5 );
  a    = t.helix().a();
  a[3] = l.b();
  // a[4] = t.charge() * l.a();//yzhang delete 2012-05-08
  //     //yzhang add 2012-05-08
  if ( _fitter.getMagneticFieldPointer() == NULL )
  {
    std::cout << "   " << __FILE__ << "   " << __LINE__ << " Could not get MagneticFieldSvc  "
              << std::endl;
  }
  const double Bz = -1000 * _fitter.getMagneticFieldPointer()->getReferField();
  if ( Bz > 0 )
  {
    a[4] = t.charge() * l.a(); // yzhang add 2012-05-08 ???
  }
  else
  {
    a[4] = -1. * t.charge() * l.a(); // yzhang add 2012-05-08  ???
  }
  // zhangy
  t._helix->a( a );
 
// #ifdef TRKRECO_WINDOW
//     TTrack * td = new TTrack(t);
//     sz.append(* td, leda_green);
// #endif
#ifdef TRKRECO_DEBUG_DETAIL
  static unsigned nTrk = 0;
  ++nTrk;
  //      for (unsigned i = 0; i < t.nCores(); i++)
  //    if (t.cores()[i]->wire()->name() == "43=220") {
  //        std::cout << "43=220 removed" << std::endl;
  //        t.remove(* t.cores()[i]);
  //    }
  HepVector3D pp;
  if ( links[0]->hit()->mc() ) pp = links[0]->hit()->mc()->hep()->p().vect();
  const HepVector3D p0 = pp;
  //      if (nTrk == 1) {
  //    Helix tmph(Point3D(0, 0, 0), p0, +1);
  //    * t._helix = tmph;
  //      }
  t.dump( "detail", "before fit" );
  std::cout << "Pdif mag=" << ( t.p() - p0 ).mag() << std::endl;
//      t.links().removeAll();
//      for (unsigned i = 0; i < BsCouTab(RECMDC_WIRHIT
//    }
#endif
 
  //...Refine...
  AList<TMLink> bad;
  int           err = _fitter.fit( t );
 
#ifdef TRKRECO_DEBUG_DETAIL
  t.dump( "detail", "after fit" );
  std::cout << "Pdif mag=" << ( t.p() - p0 ).mag() << std::endl;
#endif
 
  t.refine( bad, _maxSigmaStereo * 100. );
  err = _fitter.fit( t );
  t.refine( bad, _maxSigmaStereo * 10. );
  err = _fitter.fit( t );
  t.refine( bad, _maxSigmaStereo );
#ifdef TRKRECO_DEBUG_DETAIL
//    if (nTrk == 2) {
//      Helix tmph(Point3D(0, 0, 0), p0, +1);
//      * t._helix = tmph;
//      std::cout << "initial mom" << p0 << " is set" << std::endl;
//    }
#endif
  err = _fitter.fit( t );
 
#ifdef TRKRECO_DEBUG_DETAIL
  t.dump( "detail", "         " );
  std::cout << "Pdif mag=" << ( t.p() - p0 ).mag() << std::endl;
#endif
 
  //...Termination...
  return &t;
}

Referenced by buildStereoNew().

build() [2/3]

TTrack * TBuilder::build	(	TTrack &	t,
		const TMLine &	l ) const

build() [3/3]

TTrack * TBuilder::build	(	TTrack &	t,
		const TMLine &	l ) const

buildRphi() [1/6]

TTrack * TBuilder::buildRphi

(

const AList< TMLink > &

list

)

const

builds a r/phi track from TMLinks or from Segments.

Definition at line 101 of file TBuilder.cxx.

                                                             {
#ifdef TRKRECO_DEBUG
  std::cout << "... building rphi by links : # of links = ";
  std::cout << list.length() << std::endl;
#endif
 
  //...Classify TMLink's...
  AList<TMLink> cores;
  AList<TMLink> nonCores;
  SeparateCores( list, cores, nonCores );
 
#ifdef TRKRECO_DEBUG
  cout << "  ... cores..." << endl;
  for ( int ii = 0; ii < cores.length(); ++ii )
  {
    cout << "layer: " << cores[ii]->wire()->layerId()
         << "  local: " << cores[ii]->wire()->localId() << endl;
  }
  cout << "  ...noncores..." << endl;
  for ( int ii = 0; ii < nonCores.length(); ++ii )
  {
    cout << "layer: " << nonCores[ii]->wire()->layerId()
         << "  local: " << nonCores[ii]->wire()->localId() << endl;
  }
#endif
  //...Check # of links...
  unsigned nCores = cores.length();
  if ( nCores < _minNCores )
  {
#ifdef TRKRECO_DEBUG
    std::cout << "... building rphi failure : # of cores(=" << nCores;
    std::cout << ") is less then " << _minNCores << std::endl;
#endif
 
    return NULL;
  }
 
  //...Make a circle...
#ifdef TRKRECO_DEBUG
  std::cout << "links in list = " << list.length() << std::endl;
  std::cout << "... making a circle : # cores =" << cores.length() << std::endl;
#endif
  TCircle c( cores );
  int     err = c.fit();
  if ( err < 0 )
  {
#ifdef TRKRECO_DEBUG
    std::cout << "... building rphi failure : circle fit error = ";
    std::cout << err << std::endl;
#endif
    return NULL;
  }
  //    cout<<"radius "<<c.radius()<<" charge "<<c.charge()<<" center  "<<c.center()<<endl;
  //...Make a track...
  TTrack* t = new TTrack( c );
  //   t->dump("hits sort flag", "    ");
  AList<TMLink> bad;
  err = _fitter.fit( *t );
  // cout<<"fit 1 done err: "<<err<<endl;
  //  t->dump("hits sort flag", "    ");
  if ( err < 0 ) goto discard;
  t->refine( bad, _maxSigma * 100. );
  // cout<<"refine 1 done"<<endl;
  // t->dump("hits sort flag", "    ");
  err = _fitter.fit( *t );
  // cout<<"fit 2 done err: "<<err<<endl;
  // t->dump("hits sort flag", "    ");
  t->refine( bad, _maxSigma * 10. );
  // cout<<"refine 2 done"<<endl;
  // t->dump("hits sort flag", "    ");
  err = _fitter.fit( *t );
  // cout<<"fit 3 done err: "<<err<<endl;
  // t->dump("hits sort flag", "    ");
  t->refine( bad, _maxSigma );
  // cout<<"refine 3 done"<<endl;
  // t->dump("hits sort flag", "    ");
  err = _fitter.fit( *t );
  // cout<<"fit 4 done err: "<<err<<endl;
 
  if ( err < 0 ) goto discard;
#ifdef TRKRECO_DEBUG_DETAIL
  c.dump( "detail", "    ccl> " );
  t->dump( "detail", "    1st> " );
#endif
 
  //...Try to append non-core hits...
#ifdef TRKRECO_DEBUG
  std::cout << "... appending non-core hits : # = " << nonCores.length();
  std::cout << std::endl;
#endif
  t->appendByApproach( nonCores, _salvageLevel );
#ifdef TRKRECO_DEBUG
  t->dump( "hits sort flag", "    " );
#endif
  //   t->dump("hits sort flag", "    ");
  return t;
 
  //...Something happened...
discard:
#ifdef TRKRECO_DEBUG
  std::cout << "... building rphi failure : helix fit error = ";
  std::cout << err << std::endl;
#endif
  delete t;
  return NULL;
}

Referenced by buildRphi().

buildRphi() [2/6]

TTrack * TBuilder::buildRphi

(

const AList< TMLink > &

)

const

builds a r/phi track from TMLinks or from Segments.

buildRphi() [3/6]

TTrack * TBuilder::buildRphi

(

const AList< TMLink > &

)

const

builds a r/phi track from TMLinks or from Segments.

buildRphi() [4/6]

TTrack * TBuilder::buildRphi

(

const AList< TSegment > &

list

)

const

Definition at line 60 of file TBuilder.cxx.

                                                               {
 
#ifdef TRKRECO_DEBUG
  std::cout << "... building rphi by segments : # of segments = ";
  std::cout << list.length() << std::endl;
  for ( unsigned i = 0; i < list.length(); i++ ) list[i]->dump( "hits sort flag", "    " );
#endif
  // yuany
  //     for (unsigned i = 0; i < list.length(); i++)
  //     list[i]->dump("", "    ");
 
  //...Pick up links...
  AList<TMLink> links = Links( list );
  // yuany
  //    for(unsigned i=0;i<links.length();i++){
  //      cout<<"wire name "<<links[i]->wire()->name()<<endl;
  //    }
 
  //...Main funtion...
  TTrack* t = buildRphi( links );
 
  //...Check used segments...
  if ( t )
  {
    const AList<TMLink>& usedLinks = t->links();
    unsigned             n         = list.length();
    for ( unsigned i = 0; i < n; i++ )
    {
      TSegment&     segment = *list[i];
      AList<TMLink> used    = Links( segment, *t );
      if ( used.length() )
      {
        t->segments().append( segment );
        segment.tracks().append( t );
      }
    }
  }
 
  return t;
}

buildRphi() [5/6]

TTrack * TBuilder::buildRphi

(

const AList< TSegment > &

)

const

buildRphi() [6/6]

TTrack * TBuilder::buildRphi

(

const AList< TSegment > &

)

const

buildStereo() [1/6]

TTrack * TBuilder::buildStereo	(	const TTrack &	t,
		AList< TSegment > &	segments ) const

Definition at line 657 of file TBuilder.cxx.

                                                                                 {
 
  TTrack& t = *new TTrack( t0 );
 
#ifdef TRKRECO_DEBUG
  std::cout << "... building stereo by segments : # of segments = ";
  std::cout << segments.length() << std::endl;
  for ( unsigned i = 0; i < segments.length(); i++ )
    segments[i]->dump( "hits sort flag", "    " );
#endif
#ifdef TRKRECO_WINDOW
  sz.clear();
  sz.skip( false );
  sz.mode( 2 );
  sz.appendSz( t, segments, leda_black );
  AList<TSegment> tmps = segments;
  std::string     s    = "# of segments = " + std::to_string( int( segments.length() ) );
#endif
 
  //...Find initial line...
  TMLine* line = initialLineOld( t, segments );
  if ( !line )
  {
#ifdef TRKRECO_DEBUG
    std::cout << "... building stereo failure : no initial line found" << std::endl;
#endif
#ifdef TRKRECO_WINDOW
    s = "no initial line found : " + s;
    sz.text( s );
    sz.wait();
#endif
    return NULL;
  }
#ifdef TRKRECO_WINDOW
  sz.append( *line, leda_red );
  sz.text( s );
  sz.wait();
#endif
 
  //...Pick up links...
  AList<TMLink> links;
  unsigned      n = segments.length();
  for ( unsigned i = 0; i < n; i++ ) links.append( segments[i]->links() );
 
  //...Main funtion...
  TTrack* ts = buildStereo( t, *line, links );
 
  //...Check used segments...
  if ( ts )
  {
    AList<TMLink> usedLinks = StereoHits( t.links() );
    for ( unsigned i = 0; i < n; i++ )
    {
      TSegment&     segment = *segments[i];
      AList<TMLink> used    = Links( segment, t );
      if ( used.length() )
      {
        t.segments().append( segment );
        segment.tracks().append( t );
      }
    }
  }
 
  HepAListDeleteAll( (AList<TMLink>&)line->links() );
  delete line;
  return ts;
}

Referenced by buildStereo().

buildStereo() [2/6]

TTrack * TBuilder::buildStereo	(	const TTrack &	t,
		AList< TSegment > &	) const

buildStereo() [3/6]

TTrack * TBuilder::buildStereo	(	const TTrack &	t,
		AList< TSegment > &	) const

buildStereo() [4/6]

TTrack * TBuilder::buildStereo	(	TTrack &	t,
		TMLine &	l,
		const AList< TMLink > &	links ) const

Definition at line 527 of file TBuilder.cxx.

                                                                                            {
#ifdef TRKRECO_DEBUG
  std::cout << "... building stereo by links : # of links = ";
  std::cout << list.length() << std::endl;
#endif
 
  //...Classify TMLink's...
  AList<TMLink> cores;
  AList<TMLink> nonCores;
  SeparateCores( list, cores, nonCores );
  //...Check # of links...
  unsigned nCores = cores.length();
  if ( nCores < _minNCores )
  {
#ifdef TRKRECO_DEBUG
    std::cout << "... stereo building failure : # of cores(=" << nCores;
    std::cout << ") is less then " << _minNCores << std::endl;
#endif
    return NULL;
  }
#ifdef TRKRECO_WINDOW
  sz.appendSz( track, cores, leda_brown );
#endif
 
  //...Cal. left and right position...
  AList<TMLink> forNewLine;
  for ( unsigned i = 0; i < nCores; i++ )
  {
    TMLink& t = *cores[i];
    for ( unsigned i = 0; i < 2; i++ )
    {
      TMLink& tt = *new TMLink( t );
      tt.leftRight( i );
      int err = track.szPosition( tt );
      if ( err ) { delete &tt; }
      else
      {
        if ( line.distance( tt ) < _szLinkDistance )
        {
          tt.link( &t );
          forNewLine.append( tt );
        }
        else { delete &tt; }
      }
    }
  }
 
  //...Create new line...
#ifdef TRKRECO_DEBUG
  std::cout << "    ... creating a new line" << std::endl;
#endif
  unsigned                 nNewLine = forNewLine.length();
  TMLine                   newLine( forNewLine );
  static TRobustLineFitter fitter( "abc" );
  int                      err = fitter.fit( newLine );
  // int err = newLine.fit();
  if ( err < 0 )
  {
    HepAListDeleteAll( forNewLine );
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    ... 2nd linear fit failure. nLinks(";
    std::cout << forNewLine.length() << ")" << std::endl;
#endif
    return NULL;
  }
 
#ifdef TRKRECO_WINDOW
  sz.append( newLine, leda_green );
  sz.wait();
#endif
 
#ifdef TRKRECO_DEBUG_DETAIL
  Dump( forNewLine, "sort hits stereo", "        " );
#endif
 
  //...Decide left/right...
  AList<TMLink> toRemove;
  TMLink*       last = NULL;
  for ( unsigned i = 0; i < nNewLine; i++ )
  {
    if ( last == NULL ) last = forNewLine[i]->link();
    else
    {
      if ( last == forNewLine[i]->link() )
      {
        if ( newLine.distance( *forNewLine[i - 1] ) > newLine.distance( *forNewLine[i] ) )
          toRemove.append( forNewLine[i - 1] );
        else toRemove.append( forNewLine[i] );
        last = NULL;
      }
      else { last = forNewLine[i]->link(); }
    }
  }
  forNewLine.remove( toRemove );
  nNewLine = forNewLine.length();
 
#ifdef TRKRECO_DEBUG_DETAIL
  Dump( toRemove, "sort hits stereo", "      x " );
  Dump( forNewLine, "sort hits stereo", "        " );
#endif
 
  //...3D fit...
  for ( unsigned i = 0; i < nNewLine; i++ ) track.append( *forNewLine[i]->link() );
  Vector a( 5 );
  a    = track.helix().a();
  a[3] = newLine.b();
  a[4] = track.charge() * newLine.a();
  track._helix->a( a );
 
  //...Refine...
  AList<TMLink> bad;
  err = _fitter.fit( track );
  track.refine( bad, _maxSigmaStereo * 100. );
  err = _fitter.fit( track );
  track.refine( bad, _maxSigmaStereo * 10. );
  err = _fitter.fit( track );
  track.refine( bad, _maxSigmaStereo );
  err = _fitter.fit( track );
 
#ifdef TRKRECO_WINDOW
  sz.text( "stereo finished" );
  sz.oneShot( track, leda_blue );
#endif
 
  //...Termination...
  HepAListDeleteAll( toRemove );
  HepAListDeleteAll( forNewLine );
  return &track;
}

buildStereo() [5/6]

TTrack * TBuilder::buildStereo	(	TTrack &	t,
		TMLine &	l,
		const AList< TMLink > &	links ) const

buildStereo() [6/6]

TTrack * TBuilder::buildStereo	(	TTrack &	t,
		TMLine &	l,
		const AList< TMLink > &	links ) const

buildStereoNew() [1/3]

TTrack * TBuilder::buildStereoNew	(	const TTrack &	t,
		AList< TSegment > &	goodSegments,
		AList< TSegment > &	badSegments ) const

Definition at line 1477 of file TBuilder.cxx.

                                                                       {
 
#ifdef TRKRECO_DEBUG
  std::cout << "... building stereo by segments(new) : # of segments = ";
  std::cout << segments.length() << std::endl;
  for ( unsigned i = 0; i < segments.length(); i++ )
    segments[i]->dump( "hits sort flag", "    " );
#endif
#ifdef TRKRECO_WINDOW
  sz.clear();
#endif
  TTrack*              bestCandidate = NULL;
  AList<AList<TMLink>> poorSeeds;
  bool                 ok           = initializeForStereo( t, segments, badSegments );
  unsigned             nSuperLayers = _nSuperLayers + 1;
  //        std::cout<<"nSuperLayers in buildStereoNew:"<<_nSuperLayers<<endl;
  //    cout<<"ok of initializeForStereo: "<<ok<<endl;
 
  if ( !ok ) goto endOfBuilding;
  //...Main loop...
  while ( --nSuperLayers )
  {
 
    //...Initial line search...
    AList<TMLine> initialLines = searchInitialLines( nSuperLayers );
#ifdef TRKRECO_WINDOW
    sz.clear();
    sz.skip( false );
    sz.mode( 2 );
    sz.appendSz( t, segments, leda_black );
    AList<TSegment> tmps = segments;
    std::string     s    = "# of segments = " + std::to_string( int( segments.length() ) );
    sz.appendSz( t, _allLinks, leda_black );
    sz.append( _forLine, leda_brown );
    s = "nSprLyr=" + std::to_string( int( nSuperLayers ) ) +
        ", # of initial lines = " + std::to_string( int( initialLines.length() ) );
    for ( unsigned i = 0; i < initialLines.length(); i++ )
      sz.append( *initialLines[i], leda_red );
    sz.text( s );
    sz.wait();
#endif
    //  cout<<"length of initialLines = "<<initialLines.length()<<endl;
    if ( initialLines.length() == 0 ) continue;
    //...Line loop...
    bool     found         = false;
    unsigned nInitialLines = initialLines.length();
    //        cout<<"nInitialLines:"<<nInitialLines<<endl;
 
    for ( unsigned i = 0; i < nInitialLines; i++ )
    {
 
      //...Linear fit...
      const TMLine& line    = *initialLines[i];
      TMLine        newLine = searchLine( line );
 
      //...Skip if this is a seed of the poor result...
      if ( poorSeeds.length() )
      {
        bool poorCase = false;
        for ( unsigned j = 0; j < poorSeeds.length(); j++ )
        {
          if ( poorSeeds[j]->length() == newLine.links().length() )
          {
            AList<TMLink> tmp = *poorSeeds[j];
            tmp.remove( newLine.links() );
            if ( tmp.length() == 0 )
            {
#ifdef TRKRECO_DEBUG_DETAIL
              std::cout << "    ... This is a poor seed :"
                        << " skipped"
                        << " : # of poor seeds = " << poorSeeds.length() << std::endl;
#endif
              poorCase = true;
              break;
            }
          }
        }
        if ( poorCase ) continue;
      }
 
      //...Is this a good line...
      if ( !stereoQuality( newLine.links() ) ) continue;
 
      //...3D fit...
      TTrack* t3d = new TTrack( t );
      t3d         = build( *t3d, newLine );
      if ( t3d == NULL ) continue;
 
      //...Check quality...
      unsigned quality = stereoQuality( t3d->links() );
      //        cout<<"point3, quality:"<<quality<<endl;
 
      //...Best case...
      if ( quality == 2 )
      {
#ifdef TRKRECO_WINDOW
        sz.text( "stereo finished" );
        sz.oneShot( *t3d, leda_blue );
#endif
        if ( bestCandidate ) delete bestCandidate;
        bestCandidate = t3d;
        found         = true;
        break;
      }
 
      //...Poor case...
      AList<TMLink>* tmpL = new AList<TMLink>();
      tmpL->append( newLine.links() );
      poorSeeds.append( tmpL );
      if ( quality == 0 )
      {
#ifdef TRKRECO_WINDOW
        sz.text( "this candidate discarded" );
        sz.oneShot( *t3d, leda_black );
#endif
        delete t3d;
        continue;
      }
 
      //...Not enough...
      if ( bestCandidate )
      {
        if ( bestCandidate->cores().length() < t3d->cores().length() )
        {
#ifdef TRKRECO_WINDOW
          sz.text( "new candidate" );
          sz.append( *bestCandidate, leda_brown );
          sz.oneShot( *t3d, leda_green );
          sz.remove( *bestCandidate );
#endif
          delete bestCandidate;
          bestCandidate = t3d;
        }
        else
        {
#ifdef TRKRECO_WINDOW
          sz.text( "this candidate discarded" );
          sz.oneShot( *t3d, leda_black );
#endif
          delete t3d;
        }
      }
      else
      {
        bestCandidate = t3d;
#ifdef TRKRECO_WINDOW
        sz.text( "new candidate" );
        sz.oneShot( *bestCandidate, leda_green );
        sz.remove( *t3d );
#endif
      }
    }
 
    //...Termination of a loop...
    HepAListDeleteAll( initialLines );
    if ( found ) break;
  }
 
endOfBuilding:
  _allLinks.removeAll();
  for ( unsigned i = 0; i < 6; i++ ) HepAListDeleteAll( _links[i] );
  HepAListDeleteAll( _forLine );
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    ... # of poor seeds = " << poorSeeds.length() << std::endl;
#endif
#ifdef TRKRECO_WINDOW
  if ( bestCandidate == NULL )
  {
    sz.text( "3D failed" );
    sz.wait();
  }
#endif
 
  //...Check used segments...
  if ( bestCandidate )
  {
    AList<TMLink> usedLinks = StereoHits( bestCandidate->links() );
    for ( unsigned i = 0; i < segments.length(); i++ )
    {
      TSegment&     segment = *segments[i];
      AList<TMLink> used    = Links( segment, *bestCandidate );
      if ( used.length() )
      {
        bestCandidate->segments().append( segment );
        segment.tracks().append( bestCandidate );
      }
    }
  }
 
  if ( poorSeeds.length() ) HepAListDeleteAll( poorSeeds );
 
#ifdef TRASA_DEBUG_DETAIL
  if ( bestCandidate == NULL ) std::cout << "... building stereo(new) failed" << std::endl;
  else std::cout << "... building stereo(new) ok" << std::endl;
#endif
  return bestCandidate;
}

buildStereoNew() [2/3]

TTrack * TBuilder::buildStereoNew	(	const TTrack &	t,
		AList< TSegment > &	goodSegments,
		AList< TSegment > &	badSegments ) const

buildStereoNew() [3/3]

TTrack * TBuilder::buildStereoNew	(	const TTrack &	t,
		AList< TSegment > &	goodSegments,
		AList< TSegment > &	badSegments ) const

buildStereoOld() [1/3]

TTrack * TBuilder::buildStereoOld	(	TTrack &	t,
		const AList< TMLink > &	allLinks,
		const AList< TMLink > &	isolatedLinks ) const

buildStereoOld() [2/3]

TTrack * TBuilder::buildStereoOld	(	TTrack &	t,
		const AList< TMLink > &	allLinks,
		const AList< TMLink > &	isolatedLinks ) const

buildStereoOld() [3/3]

TTrack * TBuilder::buildStereoOld	(	TTrack &	t,
		const AList< TMLink > &	allLinks,
		const AList< TMLink > &	isolatedLinks ) const

dump() [1/3]

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

dumps debug information.

Definition at line 58 of file TBuilder.cxx.

{}

Referenced by buildRphi(), buildStereo(), buildStereoNew(), and salvage().

dump() [2/3]

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

dumps debug information.

dump() [3/3]

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

dumps debug information.

initialLine() [1/3]

TMLine * TBuilder::initialLine	(	const TTrack &	t,
		AList< TSegment > &	list ) const

makes a line.

Definition at line 308 of file TBuilder.cxx.

                                                                            {
 
  //...Select good segments in acceptance...
  AList<TSegment> bad;
  AList<TMLink>   lList;
  unsigned        n = list.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    TMLink& l   = *new TMLink();
    int     err = t.szPosition( *list[i], l );
    if ( err )
    {
      bad.append( list[i] );
      delete &l;
      continue;
    }
 
    const HepPoint3D& p = l.position();
    if ( p.y() > l.wire()->forwardPosition().z() || p.y() < l.wire()->backwardPosition().z() )
    {
      bad.append( list[i] );
      delete &l;
      continue;
    }
    lList.append( l );
  }
  list.remove( bad );
  n = list.length();
 
  //...Cal. expected # of super layers...
  unsigned               nAMaxSL  = 0;
  unsigned               nAMinSL  = 10;
  const AList<TSegment>& segments = t.segments();
  unsigned               nA       = segments.length();
  for ( unsigned i = 0; i < nA; i++ )
  {
    unsigned sl = segments[i]->superLayerId();
    if ( sl > nAMaxSL ) nAMaxSL = sl;
    if ( sl < nAMinSL ) nAMinSL = sl;
  }
  unsigned nExpected = ( nAMaxSL - nAMinSL ) / 2;
#ifdef TRKRECO_DEBUG
  std::cout << "    ... initialLine : axial super layer usage = " << nAMinSL;
  std::cout << " ~ " << nAMaxSL << std::endl;
  std::cout << "                    : expected stereo super layers = ";
  std::cout << nExpected << std::endl;
#endif
 
  TMLine*         line   = NULL;
  AList<TSegment> list0  = list;
  AList<TMLink>   lList0 = lList;
  while ( 1 )
  {
    bool last = ( NSuperLayers( lList ) <= nExpected );
 
    TMLine* line = initialLine2( t, lList );
    if ( line )
    {
      AList<TSegment> tmp = selectStereoSegment( *line, list0, lList0 );
#ifdef TRKRECO_DEBUG
      std::cout << "    ... initialLine2 : # of selected segments = ";
      std::cout << tmp.length() << std::endl;
#endif
      if ( ( tmp.length() >= nExpected ) || ( tmp.length() >= 4 ) || last )
      {
        list = tmp;
        return line;
      }
      delete line;
    }
 
    line = initialLine1( t, list, lList );
    if ( line )
    {
      AList<TSegment> tmp = selectStereoSegment( *line, list0, lList0 );
      bool            ok  = ( tmp.length() >= nExpected );
#ifdef TRKRECO_DEBUG
      std::cout << "    ... initialLine2 : # of selected segments = ";
      std::cout << tmp.length() << " : ok, last = ";
      std::cout << ok << ", " << last << std::endl;
#endif
      if ( ok )
      {
        list = tmp;
        return line;
      }
      removeFarSegment( *line, list, lList );
      delete line;
    }
    else { return NULL; }
  }
}

Referenced by buildStereoOld(), and searchLine().

initialLine() [2/3]

TMLine * TBuilder::initialLine	(	const TTrack &	,
		AList< TSegment > &	) const

makes a line.

initialLine() [3/3]

TMLine * TBuilder::initialLine	(	const TTrack &	,
		AList< TSegment > &	) const

makes a line.

initialLine1() [1/3]

TMLine * TBuilder::initialLine1	(	const TTrack &	t,
		const AList< TSegment > &	list,
		const AList< TMLink > &	lList ) const

Definition at line 240 of file TBuilder.cxx.

                                                                   {
 
  //...Check input...
  if ( list.length() < 2 ) return NULL;
 
  //...Check super layer pattern...
  AList<TSegment> sl[5];
  AList<TMLink>   tl[5];
  unsigned        n = list.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    unsigned j = list[i]->superLayerId() / 2;
    sl[j].append( list[i] );
    tl[j].append( lList[i] );
  }
#ifdef TRKRECO_DEBUG
  std::cout << "        ... initialLine1 : super layer ptn = ";
  for ( unsigned i = 0; i < 5; i++ ) std::cout << sl[i].length();
  std::cout << std::endl;
#endif
 
  //...Count single segment layer...
  unsigned nSingle = 0;
  for ( unsigned i = 0; i < 5; i++ )
    if ( sl[i].length() == 1 ) ++nSingle;
#ifdef TRKRECO_DEBUG
  std::cout << "        ... # of single segment layer = " << nSingle << std::endl;
#endif
  if ( nSingle < 2 ) return NULL;
 
  //...Happy case...
  AList<TSegment> bestCombination;
  AList<TMLink>   forLine;
  for ( unsigned i = 0; i < 5; i++ )
  {
    if ( sl[i].length() != 1 ) continue;
    bestCombination.append( sl[i] );
    forLine.append( tl[i] );
  }
  TMLine& line = *new TMLine( forLine );
  int     err  = line.fit();
  if ( err )
  {
    delete &line;
    return NULL;
  }
  return &line;
}

Referenced by initialLine(), and initialLineOld().

initialLine1() [2/3]

TMLine * TBuilder::initialLine1	(	const TTrack &	,
		const AList< TSegment > &	,
		const AList< TMLink > &	) const

initialLine1() [3/3]

TMLine * TBuilder::initialLine1	(	const TTrack &	,
		const AList< TSegment > &	,
		const AList< TMLink > &	) const

initialLine2() [1/3]

TMLine * TBuilder::initialLine2	(	const TTrack &	t,
		const AList< TMLink > &	lList ) const

Definition at line 290 of file TBuilder.cxx.

                                                                                  {
#ifdef TRKRECO_DEBUG
  std::cout << "        ... initlialLine2 : # of links = " << lList.length();
  std::cout << std::endl;
#endif
 
  //...Static object...
  static TRobustLineFitter fitter( "can you work?" );
  TMLine*                  line = new TMLine( lList );
  int                      err  = fitter.fit( *line );
  if ( err )
  {
    delete line;
    return NULL;
  }
  return line;
}

Referenced by initialLine(), and initialLineOld().

initialLine2() [2/3]

TMLine * TBuilder::initialLine2	(	const TTrack &	,
		const AList< TMLink > &	) const

initialLine2() [3/3]

TMLine * TBuilder::initialLine2	(	const TTrack &	,
		const AList< TMLink > &	) const

initialLineOld() [1/3]

TMLine * TBuilder::initialLineOld	(	const TTrack &	t,
		AList< TSegment > &	list ) const

Definition at line 401 of file TBuilder.cxx.

                                                                               {
 
  //...Select good segments in acceptance...
  AList<TSegment> bad;
  AList<TMLink>   lList;
  unsigned        n = list.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    TMLink& l   = *new TMLink();
    int     err = t.szPosition( *list[i], l );
    if ( err )
    {
      bad.append( list[i] );
      delete &l;
      continue;
    }
 
    const HepPoint3D& p = l.position();
    if ( p.y() > l.wire()->forwardPosition().z() || p.y() < l.wire()->backwardPosition().z() )
    {
      bad.append( list[i] );
      delete &l;
      continue;
    }
    lList.append( l );
  }
  list.remove( bad );
  n = list.length();
 
  //...Cal. expected # of super layers...
  unsigned               nAMaxSL  = 0;
  unsigned               nAMinSL  = 10;
  const AList<TSegment>& segments = t.segments();
  unsigned               nA       = segments.length();
  for ( unsigned i = 0; i < nA; i++ )
  {
    unsigned sl = segments[i]->superLayerId();
    if ( sl > nAMaxSL ) nAMaxSL = sl;
    if ( sl < nAMinSL ) nAMinSL = sl;
  }
  unsigned nExpected = ( nAMaxSL - nAMinSL ) / 2;
#ifdef TRKRECO_DEBUG
  std::cout << "    ... initialLine : axial super layer usage = " << nAMinSL;
  std::cout << " ~ " << nAMaxSL << std::endl;
  std::cout << "                    : expected stereo super layers = ";
  std::cout << nExpected << std::endl;
#endif
 
  TMLine*         line   = NULL;
  AList<TSegment> list0  = list;
  AList<TMLink>   lList0 = lList;
  while ( 1 )
  {
    TMLine* line = initialLine1( t, list, lList );
    if ( line )
    {
      AList<TSegment> tmp = selectStereoSegment( *line, list0, lList0 );
#ifdef TRKRECO_DEBUG
      std::cout << "    ... initialLine1 : # of selected segments = ";
      std::cout << tmp.length() << std::endl;
#endif
      if ( tmp.length() >= nExpected )
      {
        lList0.remove( line->links() );
        HepAListDeleteAll( lList0 );
        list = tmp;
        return line;
      }
      delete line;
    }
 
    line = initialLine2( t, lList );
    if ( line )
    {
      AList<TSegment> tmp  = selectStereoSegment( *line, list0, lList0 );
      bool            ok   = ( tmp.length() >= nExpected );
      bool            last = ( NSuperLayers( lList ) <= nExpected );
#ifdef TRKRECO_DEBUG
      std::cout << "    ... initialLine2 : # of selected segments = ";
      std::cout << tmp.length() << " : ok, last = ";
      std::cout << ok << ", " << last << std::endl;
#endif
      if ( ok || last )
      {
        lList0.remove( line->links() );
        HepAListDeleteAll( lList0 );
        list = tmp;
        return line;
      }
      removeFarSegment( *line, list, lList );
      delete line;
    }
    else
    {
      HepAListDeleteAll( lList0 );
      return NULL;
    }
  }
}

Referenced by buildStereo().

initialLineOld() [2/3]

TMLine * TBuilder::initialLineOld	(	const TTrack &	,
		AList< TSegment > &	) const

initialLineOld() [3/3]

TMLine * TBuilder::initialLineOld	(	const TTrack &	,
		AList< TSegment > &	) const

name() [1/3]

const std::string & TBuilder::name ( void ) const

inline

returns name.

Definition at line 130 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TBuilder.h.

{ return _name; }

Referenced by TBuilder().

name() [2/3]

const std::string & TBuilder::name

(

void

)

const

returns name.

name() [3/3]

const std::string & TBuilder::name

(

void

)

const

returns name.

removeFarSegment() [1/3]

void TBuilder::removeFarSegment	(	const TMLine &	line,
		AList< TSegment > &	list,
		AList< TMLink > &	szList ) const

Definition at line 513 of file TBuilder.cxx.

                                                               {
  unsigned n           = list.length();
  float    maxDistance = 0.;
  unsigned maxId       = 0;
  for ( unsigned i = 0; i < n; i++ )
  {
    float distance = line.distance( *szList[i] );
    if ( distance > maxDistance ) maxId = i;
  }
  list.remove( maxId );
  szList.remove( maxId );
}

Referenced by initialLine(), and initialLineOld().

removeFarSegment() [2/3]

void TBuilder::removeFarSegment	(	const TMLine &	,
		AList< TSegment > &	,
		AList< TMLink > &	) const

removeFarSegment() [3/3]

void TBuilder::removeFarSegment	(	const TMLine &	,
		AList< TSegment > &	,
		AList< TMLink > &	) const

salvage() [1/6]

void TBuilder::salvage	(	TTrack &	t,
		AList< TMLink > &	hits ) const

salvages hits.

Definition at line 208 of file TBuilder.cxx.

                                                             {
 
#ifdef TRKRECO_DEBUG
  std::cout << "... salvaging(TBuilder) : # of given hits=" << hits.length();
  std::cout << ", salvage level=" << _salvageLevel << std::endl;
  Dump( hits, "hits sort flag" );
#endif
 
  unsigned nHits = hits.length();
  if ( nHits == 0 ) return;
 
  //...Try to append this hit...
  t.appendByApproach( hits, _salvageLevel );
  _fitter.fit( t );
}

Referenced by salvage().

salvage() [2/6]

void TBuilder::salvage	(	TTrack &	t,
		AList< TMLink > &	hits ) const

salvages hits.

salvage() [3/6]

void TBuilder::salvage	(	TTrack &	t,
		AList< TMLink > &	hits ) const

salvages hits.

salvage() [4/6]

void TBuilder::salvage	(	TTrack &	t,
		AList< TSegment > &	segments ) const

Definition at line 224 of file TBuilder.cxx.

                                                               {
 
#ifdef TRKRECO_DEBUG
  std::cout << "... salvaging by segments : # of segments = ";
  std::cout << list.length() << std::endl;
  for ( unsigned i = 0; i < list.length(); i++ ) list[i]->dump( "hits sort flag", "    " );
#endif
 
  //...Pick up links...
  AList<TMLink> links;
  unsigned      n = list.length();
  for ( unsigned i = 0; i < n; i++ ) links.append( list[i]->links() );
 
  salvage( t, links );
}

salvage() [5/6]

void TBuilder::salvage	(	TTrack &	t,
		AList< TSegment > &	segments ) const

salvage() [6/6]

void TBuilder::salvage	(	TTrack &	t,
		AList< TSegment > &	segments ) const

searchInitialLines() [1/3]

AList< TMLine > TBuilder::searchInitialLines

(

unsigned

nSuperLayers

)

const

Definition at line 725 of file TBuilder.cxx.

                                                                          {
#ifdef TRKRECO_DEBUG
  std::cout << "    ... searchInitialLines : # of segments in sl : ";
  for ( unsigned i = 0; i < 5; i++ ) std::cout << _links[i].length() << ",";
  std::cout << std::endl
            << "                           : max # of super layers=" << nSuperLayerMax
            << std::endl;
#endif
 
  AList<TMLine> lines;
  if ( nSuperLayerMax > 5 ) lines.append( searchLines6() );
  else if ( nSuperLayerMax > 4 ) lines.append( searchLines5() );
  else if ( nSuperLayerMax > 3 ) lines.append( searchLines4() );
  else if ( nSuperLayerMax > 2 ) lines.append( searchLines3() );
  else if ( nSuperLayerMax > 1 ) lines.append( searchLines2() );
  //      else
  //    lines.append(searchLines1());
 
  lines.sort( SortByB );
  return lines;
}

Referenced by buildStereoNew().

searchInitialLines() [2/3]

AList< TMLine > TBuilder::searchInitialLines

(

unsigned

nSuperLayers

)

const

searchInitialLines() [3/3]

AList< TMLine > TBuilder::searchInitialLines

(

unsigned

nSuperLayers

)

const

searchLine() [1/3]

TMLine TBuilder::searchLine

(

const TMLine &

initialLine

)

const

Definition at line 1109 of file TBuilder.cxx.

                                                             {
#ifdef TRKRECO_WINDOW
  std::string old = sz.text();
#endif
 
  static float _szLinkMaxDistance = 5;
 
  unsigned nLinks = _forLine.length();
  TMLine   line0  = initialLine;
  line0.remove( initialLine.links() );
#ifdef TRKRECO_DEBUG
  line0.fitted( true );
#endif
  unsigned      nGoodLinks0 = 0;
  AList<TMLink> targets;
  unsigned      nLoops = 0;
  while ( 1 )
  {
    ++nLoops;
 
    unsigned nGoodLinks = 0;
    for ( unsigned i = 0; i < nLinks; i++ )
    {
      float distance = line0.distance( *_forLine[i] );
      if ( distance < _szLinkDistance ) targets.append( _forLine[i] );
      if ( distance < _szLinkMaxDistance ) { ++nGoodLinks; }
    }
 
#ifdef TRKRECO_DEBUG
    std::cout << "    ... searchLine : # of close hits(last) = " << nGoodLinks0
              << ", # of close hits = " << nGoodLinks << std::endl;
#endif
 
    //...Check condition...
    if ( nGoodLinks0 == nGoodLinks ) break;
 
    //...Do fit...
    TMLine                   newLine( targets );
    static TRobustLineFitter fitter( "robust fitter" );
    int                      err = fitter.fit( newLine );
    if ( err )
    {
#ifdef TRKRECO_WINDOW
      std::cout << "    ... searchLine : failed to fit" << std::endl;
      std::string s = "line search failed";
      sz.text( s );
      sz.wait();
#endif
      break;
    }
 
    //...To protect infinite loop...
    if ( nLoops > 10 )
    {
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "        reached to max # of loops(10) : break";
#endif
      break;
    }
 
#ifdef TRKRECO_WINDOW
    std::string ss = ", # of close hits(last) = " + std::to_string( int( nGoodLinks0 ) ) +
                     ", # of close hits = " + std::to_string( int( nGoodLinks ) );
    sz.append( line0, leda_brown );
    sz.append( newLine, leda_green );
    sz.text( old + ", nLoops = " + std::to_string( int( nLoops ) ) + ss );
    sz.wait();
    sz.remove( line0 );
    sz.remove( newLine );
#endif
 
    line0       = newLine;
    nGoodLinks0 = nGoodLinks;
    targets.removeAll();
  }
 
  return line0;
}

Referenced by buildStereoNew().

searchLine() [2/3]

TMLine TBuilder::searchLine

(

const TMLine &

initialLine

)

const

searchLine() [3/3]

TMLine TBuilder::searchLine

(

const TMLine &

initialLine

)

const

searchLines1() [1/3]

AList< TMLine > TBuilder::searchLines1

(

void

)

const

Definition at line 1094 of file TBuilder.cxx.

                                                 {
  AList<TMLine> lines;
 
  TMLine& line = *new TMLine( _forLine );
  int     err  = line.fit();
  if ( err ) { delete &line; }
  else { lines.append( line ); }
 
#ifdef TRKRECO_DEBUG
  std::cout << "        ... searchLines1 : # of lines found = " << lines.length() << std::endl;
#endif
 
  return lines;
}

searchLines1() [2/3]

AList< TMLine > TBuilder::searchLines1

(

void

)

const

searchLines1() [3/3]

AList< TMLine > TBuilder::searchLines1

(

void

)

const

searchLines2() [1/3]

AList< TMLine > TBuilder::searchLines2

(

void

)

const

Definition at line 1024 of file TBuilder.cxx.

                                                 {
  AList<TMLine> lines;
 
  const AList<TMLink>* l[2];
  l[0] = l[1] = nullptr; // add initialize 25-05-15
  for ( unsigned skip = 0; skip < 15; skip++ )
  {
    if ( skip == 0 )
    {
      l[0] = &_links[0];
      l[1] = &_links[1];
    }
    else if ( skip == 1 ) { l[1] = &_links[2]; }
    else if ( skip == 2 ) { l[1] = &_links[3]; }
    else if ( skip == 3 ) { l[1] = &_links[4]; }
    else if ( skip == 4 ) { l[1] = &_links[5]; }
    else if ( skip == 5 )
    {
      l[0] = &_links[1];
      l[1] = &_links[2];
    }
    else if ( skip == 6 ) { l[1] = &_links[3]; }
    else if ( skip == 7 ) { l[1] = &_links[4]; }
    else if ( skip == 8 ) { l[1] = &_links[5]; }
    else if ( skip == 9 )
    {
      l[0] = &_links[2];
      l[1] = &_links[3];
    }
    else if ( skip == 10 ) { l[1] = &_links[4]; }
    else if ( skip == 11 ) { l[1] = &_links[5]; }
    else if ( skip == 12 )
    {
      l[0] = &_links[3];
      l[1] = &_links[4];
    }
    else if ( skip == 13 ) { l[1] = &_links[5]; }
    else if ( skip == 14 ) { l[0] = &_links[4]; }
 
    unsigned n[2];
    for ( unsigned i = 0; i < 2; i++ ) n[i] = l[i]->length();
 
    for ( unsigned i0 = 0; i0 < n[0]; i0++ )
    {
      for ( unsigned i1 = 0; i1 < n[1]; i1++ )
      {
        AList<TMLink> forLine;
        forLine.append( ( *l[0] )[i0] );
        forLine.append( ( *l[1] )[i1] );
 
        TMLine& line = *new TMLine( forLine );
        int     err  = line.fit();
        if ( err )
        {
          delete &line;
          continue;
        }
 
        lines.append( line );
      }
    }
  }
 
#ifdef TRKRECO_DEBUG
  std::cout << "        ... searchLines2 : # of lines found = " << lines.length() << std::endl;
#endif
 
  return lines;
}

Referenced by searchInitialLines().

searchLines2() [2/3]

AList< TMLine > TBuilder::searchLines2

(

void

)

const

searchLines2() [3/3]

AList< TMLine > TBuilder::searchLines2

(

void

)

const

searchLines3() [1/3]

AList< TMLine > TBuilder::searchLines3

(

void

)

const

Definition at line 944 of file TBuilder.cxx.

                                                 {
  AList<TMLine> lines;
 
  const AList<TMLink>* l[3];
  l[0] = l[1] = l[2] = nullptr; // add initialize 25-05-15
  for ( unsigned skip = 0; skip < 20; skip++ )
  {
    if ( skip == 0 )
    {
      l[0] = &_links[0];
      l[1] = &_links[1];
      l[2] = &_links[2];
    }
    else if ( skip == 1 ) { l[2] = &_links[3]; }
    else if ( skip == 2 ) { l[2] = &_links[4]; }
    else if ( skip == 3 ) { l[2] = &_links[5]; }
    else if ( skip == 4 ) { l[1] = &_links[2]; }
    else if ( skip == 5 ) { l[2] = &_links[4]; }
    else if ( skip == 6 ) { l[2] = &_links[3]; }
    else if ( skip == 7 )
    {
      l[1] = &_links[3];
      l[2] = &_links[4];
    }
    else if ( skip == 8 ) { l[2] = &_links[5]; }
    else if ( skip == 9 ) { l[1] = &_links[4]; }
    else if ( skip == 10 ) { l[0] = &_links[1]; }
    else if ( skip == 11 ) { l[1] = &_links[3]; }
    else if ( skip == 12 ) { l[2] = &_links[4]; }
    else if ( skip == 13 ) { l[1] = &_links[2]; }
    else if ( skip == 14 ) { l[2] = &_links[3]; }
    else if ( skip == 15 ) { l[2] = &_links[5]; }
    else if ( skip == 16 )
    {
      l[0] = &_links[2];
      l[1] = &_links[3];
    }
    else if ( skip == 17 ) { l[2] = &_links[3]; }
    else if ( skip == 18 )
    {
      l[1] = &_links[4];
      l[2] = &_links[5];
    }
    else if ( skip == 19 ) { l[0] = &_links[3]; }
 
    unsigned n[3];
    for ( unsigned i = 0; i < 3; i++ ) n[i] = l[i]->length();
 
    for ( unsigned i0 = 0; i0 < n[0]; i0++ )
    {
      for ( unsigned i1 = 0; i1 < n[1]; i1++ )
      {
        for ( unsigned i2 = 0; i2 < n[2]; i2++ )
        {
          AList<TMLink> forLine;
          forLine.append( ( *l[0] )[i0] );
          forLine.append( ( *l[1] )[i1] );
          forLine.append( ( *l[2] )[i2] );
 
          TMLine& line = *new TMLine( forLine );
          int     err  = line.fit();
          if ( err )
          {
            delete &line;
            continue;
          }
 
          lines.append( line );
        }
      }
    }
  }
 
#ifdef TRKRECO_DEBUG
  std::cout << "        ... searchLines3 : # of lines found = " << lines.length() << std::endl;
#endif
 
  return lines;
}

Referenced by searchInitialLines().

searchLines3() [2/3]

AList< TMLine > TBuilder::searchLines3

(

void

)

const

searchLines3() [3/3]

AList< TMLine > TBuilder::searchLines3

(

void

)

const

searchLines4() [1/3]

AList< TMLine > TBuilder::searchLines4

(

void

)

const

Definition at line 861 of file TBuilder.cxx.

                                                 {
  AList<TMLine> lines;
 
  const AList<TMLink>* l[4];
  for (int i = 0; i < 4; i++ ) l[i] = nullptr; // add initialize 25-05-15
  for ( unsigned skip = 0; skip < 15; skip++ )
  {
    if ( skip == 0 )
    {
      l[0] = &_links[0];
      l[1] = &_links[1];
      l[2] = &_links[2];
      l[3] = &_links[3];
    }
    else if ( skip == 1 ) { l[3] = &_links[4]; }
    else if ( skip == 2 ) { l[3] = &_links[5]; }
    else if ( skip == 3 )
    {
      l[2] = &_links[3];
      l[3] = &_links[4];
    }
    else if ( skip == 4 ) { l[3] = &_links[5]; }
    else if ( skip == 5 ) { l[2] = &_links[4]; }
    else if ( skip == 6 )
    {
      l[1] = &_links[2];
      l[2] = &_links[3];
      l[3] = &_links[4];
    }
    else if ( skip == 7 ) { l[3] = &_links[5]; }
    else if ( skip == 8 ) { l[2] = &_links[4]; }
    else if ( skip == 9 ) { l[1] = &_links[3]; }
    else if ( skip == 10 )
    {
      l[0] = &_links[1];
      l[1] = &_links[2];
      l[2] = &_links[3];
      l[3] = &_links[4];
    }
    else if ( skip == 11 ) { l[3] = &_links[5]; }
    else if ( skip == 12 ) { l[2] = &_links[4]; }
    else if ( skip == 13 ) { l[1] = &_links[3]; }
    else if ( skip == 14 ) { l[0] = &_links[2]; }
 
    unsigned n[4];
    for ( unsigned i = 0; i < 4; i++ ) n[i] = l[i]->length();
 
    for ( unsigned i0 = 0; i0 < n[0]; i0++ )
    {
      for ( unsigned i1 = 0; i1 < n[1]; i1++ )
      {
        for ( unsigned i2 = 0; i2 < n[2]; i2++ )
        {
          for ( unsigned i3 = 0; i3 < n[3]; i3++ )
          {
            AList<TMLink> forLine;
            forLine.append( ( *l[0] )[i0] );
            forLine.append( ( *l[1] )[i1] );
            forLine.append( ( *l[2] )[i2] );
            forLine.append( ( *l[3] )[i3] );
 
            TMLine& line = *new TMLine( forLine );
            int     err  = line.fit();
            if ( err )
            {
              delete &line;
              continue;
            }
 
            lines.append( line );
          }
        }
      }
    }
  }
 
#ifdef TRKRECO_DEBUG
  std::cout << "        ... searchLines4 : # of lines found = " << lines.length() << std::endl;
#endif
 
  return lines;
}

Referenced by searchInitialLines().

searchLines4() [2/3]

AList< TMLine > TBuilder::searchLines4

(

void

)

const

searchLines4() [3/3]

AList< TMLine > TBuilder::searchLines4

(

void

)

const

searchLines5() [1/3]

AList< TMLine > TBuilder::searchLines5

(

void

)

const

Definition at line 796 of file TBuilder.cxx.

                                                 {
  AList<TMLine> lines;
 
  const AList<TMLink>* l[5];
  for (int i = 0; i < 5; i++ ) l[i] = nullptr; // add initialize 25-05-15
  for ( unsigned skip = 0; skip < 6; skip++ )
  {
    if ( skip == 0 )
    {
      l[0] = &_links[1];
      l[1] = &_links[2];
      l[2] = &_links[3];
      l[3] = &_links[4];
      l[4] = &_links[5];
    }
    else if ( skip == 1 ) { l[0] = &_links[0]; }
    else if ( skip == 2 ) { l[1] = &_links[1]; }
    else if ( skip == 3 ) { l[2] = &_links[2]; }
    else if ( skip == 4 ) { l[3] = &_links[3]; }
    else if ( skip == 5 ) { l[4] = &_links[4]; }
 
    unsigned n[5];
    for ( unsigned i = 0; i < 5; i++ ) n[i] = l[i]->length();
 
    for ( unsigned i0 = 0; i0 < n[0]; i0++ )
    {
      for ( unsigned i1 = 0; i1 < n[1]; i1++ )
      {
        for ( unsigned i2 = 0; i2 < n[2]; i2++ )
        {
          for ( unsigned i3 = 0; i3 < n[3]; i3++ )
          {
            for ( unsigned i4 = 0; i4 < n[4]; i4++ )
            {
 
              AList<TMLink> forLine;
              forLine.append( ( *l[0] )[i0] );
              forLine.append( ( *l[1] )[i1] );
              forLine.append( ( *l[2] )[i2] );
              forLine.append( ( *l[3] )[i3] );
              forLine.append( ( *l[4] )[i4] );
 
              TMLine& line = *new TMLine( forLine );
              int     err  = line.fit();
              if ( err )
              {
                delete &line;
                continue;
              }
 
              lines.append( line );
            }
          }
        }
      }
    }
  }
 
#ifdef TRKRECO_DEBUG
  std::cout << "        ... searchLines5 : # of lines found = " << lines.length() << std::endl;
#endif
 
  return lines;
}

Referenced by searchInitialLines().

searchLines5() [2/3]

AList< TMLine > TBuilder::searchLines5

(

void

)

const

searchLines5() [3/3]

AList< TMLine > TBuilder::searchLines5

(

void

)

const

searchLines6() [1/3]

AList< TMLine > TBuilder::searchLines6

(

void

)

const

Definition at line 747 of file TBuilder.cxx.

                                                 {
  unsigned n[6];
  for ( unsigned i = 0; i < 6; i++ ) n[i] = _links[i].length();
 
  AList<TMLine> lines;
  for ( unsigned i0 = 0; i0 < n[0]; i0++ )
  {
    for ( unsigned i1 = 0; i1 < n[1]; i1++ )
    {
      for ( unsigned i2 = 0; i2 < n[2]; i2++ )
      {
        for ( unsigned i3 = 0; i3 < n[3]; i3++ )
        {
          for ( unsigned i4 = 0; i4 < n[4]; i4++ )
          {
            for ( unsigned i5 = 0; i5 < n[5]; i5++ )
            {
 
              AList<TMLink> forLine;
              forLine.append( _links[0][i0] );
              forLine.append( _links[1][i1] );
              forLine.append( _links[2][i2] );
              forLine.append( _links[3][i3] );
              forLine.append( _links[4][i4] );
              forLine.append( _links[5][i5] );
 
              TMLine& line = *new TMLine( forLine );
              int     err  = line.fit();
              if ( err )
              {
                delete &line;
                continue;
              }
 
              lines.append( line );
            }
          }
        }
      }
    }
  }
 
#ifdef TRKRECO_DEBUG
  std::cout << "        ... searchLines5 : # of lines found = " << lines.length() << std::endl;
#endif
 
  return lines;
}

Referenced by searchInitialLines().

searchLines6() [2/3]

AList< TMLine > TBuilder::searchLines6

(

void

)

const

searchLines6() [3/3]

AList< TMLine > TBuilder::searchLines6

(

void

)

const

selectStereoSegment() [1/3]

AList< TSegment > TBuilder::selectStereoSegment	(	const TMLine &	line,
		const AList< TSegment > &	list,
		const AList< TMLink > &	szList ) const

Definition at line 501 of file TBuilder.cxx.

                                                                                   {
  AList<TSegment> outList;
  unsigned        n = list.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    float distance = line.distance( *szList[i] );
    if ( distance < _szSegmentDistance ) outList.append( list[i] );
  }
  return outList;
}

Referenced by initialLine(), and initialLineOld().

selectStereoSegment() [2/3]

AList< TSegment > TBuilder::selectStereoSegment	(	const TMLine &	line,
		const AList< TSegment > &	list,
		const AList< TMLink > &	szList ) const

selectStereoSegment() [3/3]

AList< TSegment > TBuilder::selectStereoSegment	(	const TMLine &	line,
		const AList< TSegment > &	list,
		const AList< TMLink > &	szList ) const

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TBuilder.h
• InstallArea/x86_64-el9-gcc13-opt/include/TrkReco/TBuilder.h
• Reconstruction/TrkReco/include/TrkReco/TBuilder.h
• TBuilder.cxx