TTrackManager Class Reference

A manager of TTrack information to make outputs as MdcRec_trk. More...

#include <TTrackManager.h>

Inheritance diagram for TTrackManager:

Public Member Functions
	TTrackManager ()
	Constructor.
virtual	~TTrackManager ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
std::string	name (void) const
	returns name.
std::string	version (void) const
	returns version.
const AList< TTrack > &	allTracks (void) const
	returns a list of all reconstructed tracks including junks.
const AList< TTrack > &	tracks (void) const
	returns a list of reconstructed tracks.
const AList< TTrack > &	tracks2D (void) const
	returns a list of 2D tracks.
const AList< TTrack > &	tracksFinal (void) const
	returns a list of tracks writen to MdcRec_trk.
void	fittingFlag (unsigned)
	sets fitting flag.
double	maxMomentum (double)
	sets the max. momentum.
int	debugLevel (void) const
	returns/sets debug level.
int	debugLevel (int)
float	paraT0 (void) const
void	clear (void)
	clears all internal information.
void	append (AList< TTrack > &list)
	appends (2D) tracks. 'list' will be cleaned up.
void	append2D (AList< TTrack > &list)
void	refit (void)
	refits tracks.
void	finish (void)
	finishes tracks.
void	movePivot (void)
	moves pivot of tracks.
StatusCode	determineT0 (unsigned level, unsigned nMaxTracks)
	determines T0 and refit all tracks.
void	setCurlerFlags (void)
	tests for curlers.
void	sortTracksByQuality (void)
	sorts tracks.
void	sortTracksByPt (void)
void	clearTables (void) const
	clears tables.
StatusCode	makeTds (RecMdcTrackCol *trackList, RecMdcHitCol *hitList, int tkStat=1, int runge=0, int cal=0)
	stores track info. into TDS. by Zang Shilei
void	saveTables (void)
	stores track info. into Panther table.
void	saveMCTables (void) const
	stores MC track info. into Panther table.
void	sortBanksByPt (void) const
	sorts RECMDC_TRK tables.
void	treatCurler (MdcTrk &curl, MdcRec_trk_add &cdc, unsigned flag) const
	final decision for a curler.
void	determineIP (void)
	determines IP.
void	maskCurlHits (const AList< TMDCWireHit > &axial, const AList< TMDCWireHit > &stereo, const AList< TTrack > &tracks) const
	masks hits on found curl tracks.
void	salvageAssociateHits (const AList< TMDCWireHit > &, float maxSigma2)
	salvages hits for dE/dx(not for track fitting).
TMLink &	divide (const TTrack &t, AList< TMLink > *l) const
TMLink &	divideByIp (const TTrack &t, AList< TMLink > *l) const
void	removeHitsAcrossOverIp (AList< TMLink > &) const
TTrack *	closest (const AList< TTrack > &, const TMDCWireHit &) const
	returns a track which is the closest to a hit.
void	salvage (const AList< TMDCWireHit > &) const
	salvages remaining hits.
void	mask (void) const
	masks hits out which are in tail of curly tracks.
void	maskNormal (TTrack &) const
void	maskCurl (TTrack &) const
void	maskOut (TTrack &, const AList< TMLink > &) const
void	maskMultiHits (TTrack &) const
void	merge (void)
void	sett0bes (double t0)
	TTrackManager ()
	Constructor.
virtual	~TTrackManager ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
std::string	name (void) const
	returns name.
std::string	version (void) const
	returns version.
const AList< TTrack > &	allTracks (void) const
	returns a list of all reconstructed tracks including junks.
const AList< TTrack > &	tracks (void) const
	returns a list of reconstructed tracks.
const AList< TTrack > &	tracks2D (void) const
	returns a list of 2D tracks.
const AList< TTrack > &	tracksFinal (void) const
	returns a list of tracks writen to MdcRec_trk.
void	fittingFlag (unsigned)
	sets fitting flag.
double	maxMomentum (double)
	sets the max. momentum.
int	debugLevel (void) const
	returns/sets debug level.
int	debugLevel (int)
float	paraT0 (void) const
void	clear (void)
	clears all internal information.
void	append (AList< TTrack > &list)
	appends (2D) tracks. 'list' will be cleaned up.
void	append2D (AList< TTrack > &list)
void	refit (void)
	refits tracks.
void	finish (void)
	finishes tracks.
void	movePivot (void)
	moves pivot of tracks.
StatusCode	determineT0 (unsigned level, unsigned nMaxTracks)
	determines T0 and refit all tracks.
void	setCurlerFlags (void)
	tests for curlers.
void	sortTracksByQuality (void)
	sorts tracks.
void	sortTracksByPt (void)
void	clearTables (void) const
	clears tables.
StatusCode	makeTds (RecMdcTrackCol *trackList, RecMdcHitCol *hitList, int tkStat=1, int runge=0, int cal=0)
	stores track info. into TDS. by Zang Shilei
void	saveTables (void)
	stores track info. into Panther table.
void	saveMCTables (void) const
	stores MC track info. into Panther table.
void	sortBanksByPt (void) const
	sorts RECMDC_TRK tables.
void	treatCurler (MdcTrk &curl, MdcRec_trk_add &cdc, unsigned flag) const
	final decision for a curler.
void	determineIP (void)
	determines IP.
void	maskCurlHits (const AList< TMDCWireHit > &axial, const AList< TMDCWireHit > &stereo, const AList< TTrack > &tracks) const
	masks hits on found curl tracks.
void	salvageAssociateHits (const AList< TMDCWireHit > &, float maxSigma2)
	salvages hits for dE/dx(not for track fitting).
TMLink &	divide (const TTrack &t, AList< TMLink > *l) const
TMLink &	divideByIp (const TTrack &t, AList< TMLink > *l) const
void	removeHitsAcrossOverIp (AList< TMLink > &) const
TTrack *	closest (const AList< TTrack > &, const TMDCWireHit &) const
	returns a track which is the closest to a hit.
void	salvage (const AList< TMDCWireHit > &) const
	salvages remaining hits.
void	mask (void) const
	masks hits out which are in tail of curly tracks.
void	maskNormal (TTrack &) const
void	maskCurl (TTrack &) const
void	maskOut (TTrack &, const AList< TMLink > &) const
void	maskMultiHits (TTrack &) const
void	merge (void)
void	sett0bes (double t0)
	TTrackManager ()
	Constructor.
virtual	~TTrackManager ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
std::string	name (void) const
	returns name.
std::string	version (void) const
	returns version.
const AList< TTrack > &	allTracks (void) const
	returns a list of all reconstructed tracks including junks.
const AList< TTrack > &	tracks (void) const
	returns a list of reconstructed tracks.
const AList< TTrack > &	tracks2D (void) const
	returns a list of 2D tracks.
const AList< TTrack > &	tracksFinal (void) const
	returns a list of tracks writen to MdcRec_trk.
void	fittingFlag (unsigned)
	sets fitting flag.
double	maxMomentum (double)
	sets the max. momentum.
int	debugLevel (void) const
	returns/sets debug level.
int	debugLevel (int)
float	paraT0 (void) const
void	clear (void)
	clears all internal information.
void	append (AList< TTrack > &list)
	appends (2D) tracks. 'list' will be cleaned up.
void	append2D (AList< TTrack > &list)
void	refit (void)
	refits tracks.
void	finish (void)
	finishes tracks.
void	movePivot (void)
	moves pivot of tracks.
StatusCode	determineT0 (unsigned level, unsigned nMaxTracks)
	determines T0 and refit all tracks.
void	setCurlerFlags (void)
	tests for curlers.
void	sortTracksByQuality (void)
	sorts tracks.
void	sortTracksByPt (void)
void	clearTables (void) const
	clears tables.
StatusCode	makeTds (RecMdcTrackCol *trackList, RecMdcHitCol *hitList, int tkStat=1, int runge=0, int cal=0)
	stores track info. into TDS. by Zang Shilei
void	saveTables (void)
	stores track info. into Panther table.
void	saveMCTables (void) const
	stores MC track info. into Panther table.
void	sortBanksByPt (void) const
	sorts RECMDC_TRK tables.
void	treatCurler (MdcTrk &curl, MdcRec_trk_add &cdc, unsigned flag) const
	final decision for a curler.
void	determineIP (void)
	determines IP.
void	maskCurlHits (const AList< TMDCWireHit > &axial, const AList< TMDCWireHit > &stereo, const AList< TTrack > &tracks) const
	masks hits on found curl tracks.
void	salvageAssociateHits (const AList< TMDCWireHit > &, float maxSigma2)
	salvages hits for dE/dx(not for track fitting).
TMLink &	divide (const TTrack &t, AList< TMLink > *l) const
TMLink &	divideByIp (const TTrack &t, AList< TMLink > *l) const
void	removeHitsAcrossOverIp (AList< TMLink > &) const
TTrack *	closest (const AList< TTrack > &, const TMDCWireHit &) const
	returns a track which is the closest to a hit.
void	salvage (const AList< TMDCWireHit > &) const
	salvages remaining hits.
void	mask (void) const
	masks hits out which are in tail of curly tracks.
void	maskNormal (TTrack &) const
void	maskCurl (TTrack &) const
void	maskOut (TTrack &, const AList< TMLink > &) const
void	maskMultiHits (TTrack &) const
void	merge (void)
void	sett0bes (double t0)
Public Member Functions inherited from TUpdater
	TUpdater ()
	Constructor.
virtual	~TUpdater ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual bool	updated (void) const
	returns true if an object is updated.
virtual void	update (void)
	updates an object.
virtual bool	initialized (void) const
	returns true if an object is initialized.
virtual void	initialize (void)
	initializes an object.
	TUpdater ()
	Constructor.
virtual	~TUpdater ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual bool	updated (void) const
	returns true if an object is updated.
virtual void	update (void)
	updates an object.
virtual bool	initialized (void) const
	returns true if an object is initialized.
virtual void	initialize (void)
	initializes an object.
	TUpdater ()
	Constructor.
virtual	~TUpdater ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual bool	updated (void) const
	returns true if an object is updated.
virtual void	update (void)
	updates an object.
virtual bool	initialized (void) const
	returns true if an object is initialized.
virtual void	initialize (void)
	initializes an object.

Static Public Member Functions
static bool	goodTrack (const TTrack &, bool track2D=false)
	checks goodness of a track.
static void	maskBadHits (const AList< TTrack > &, float maxSigma2)
	masks hits with large chisq as associated hits. Pull in TMLink is used.
static bool	goodTrack (const TTrack &, bool track2D=false)
	checks goodness of a track.
static void	maskBadHits (const AList< TTrack > &, float maxSigma2)
	masks hits with large chisq as associated hits. Pull in TMLink is used.
static bool	goodTrack (const TTrack &, bool track2D=false)
	checks goodness of a track.
static void	maskBadHits (const AList< TTrack > &, float maxSigma2)
	masks hits with large chisq as associated hits. Pull in TMLink is used.

Detailed Description

A manager of TTrack information to make outputs as MdcRec_trk.

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

Constructor & Destructor Documentation

TTrackManager() [1/3]

TTrackManager::TTrackManager

(

)

Constructor.

Definition at line 75 of file TTrackManager.cxx.

    : _maxMomentum( 10. )
    , _sigmaCurlerMergeTest( sqrt( 100. ) )
    , _nCurlerMergeTest( 4 )
    , _debugLevel( 0 )
    , _fitter( "TTrackManager Fitter" )
    , _cFitter( "TTrackManager 2D Fitter" )
    , _s( 0 )
    , m_pmgnIMF( nullptr )
    , m_rawDataProviderSvc( nullptr ) {}

~TTrackManager() [1/3]

TTrackManager::~TTrackManager ( )

virtual

Destructor.

Definition at line 86 of file TTrackManager.cxx.

{}

TTrackManager() [2/3]

TTrackManager::TTrackManager

(

)

Constructor.

~TTrackManager() [2/3]

virtual TTrackManager::~TTrackManager ( )

virtual

Destructor.

TTrackManager() [3/3]

TTrackManager::TTrackManager

(

)

Constructor.

~TTrackManager() [3/3]

virtual TTrackManager::~TTrackManager ( )

virtual

Destructor.

Member Function Documentation

allTracks() [1/3]

const AList< TTrack > & TTrackManager::allTracks ( void ) const

inline

returns a list of all reconstructed tracks including junks.

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

{ return _tracksAll; }

allTracks() [2/3]

const AList< TTrack > & TTrackManager::allTracks

(

void

)

const

returns a list of all reconstructed tracks including junks.

allTracks() [3/3]

const AList< TTrack > & TTrackManager::allTracks

(

void

)

const

returns a list of all reconstructed tracks including junks.

append() [1/3]

void TTrackManager::append

(

AList< TTrack > &

list

)

appends (2D) tracks. 'list' will be cleaned up.

Definition at line 862 of file TTrackManager.cxx.

                                                {
  _tracksAll.append( list );
  _tracks.append( selectGoodTracks( list ) );
  list.removeAll();
}

append() [2/3]

void TTrackManager::append

(

AList< TTrack > &

list

)

appends (2D) tracks. 'list' will be cleaned up.

append() [3/3]

void TTrackManager::append

(

AList< TTrack > &

list

)

appends (2D) tracks. 'list' will be cleaned up.

append2D() [1/3]

void TTrackManager::append2D

(

AList< TTrack > &

list

)

Definition at line 868 of file TTrackManager.cxx.

                                                  {
  _tracksAll.append( list );
  _tracks2D.append( selectGoodTracks( list, true ) );
  _tracks2D.sort( SortByPt );
  list.removeAll();
}

append2D() [2/3]

void TTrackManager::append2D

(

AList< TTrack > &

list

)

append2D() [3/3]

void TTrackManager::append2D

(

AList< TTrack > &

list

)

clear() [1/3]

void TTrackManager::clear ( void )

virtual

clears all internal information.

Reimplemented from TUpdater.

Definition at line 826 of file TTrackManager.cxx.

                                {
  HepAListDeleteAll( _tracksAll );
  _tracks.removeAll();
  _tracks2D.removeAll();
  _tracksFinal.removeAll();
  HepAListDeleteAll( _associateHits );
  static bool first = true;
  if ( first )
  {
    first = false;
    int size;
    //  _s = (struct summary *)
    //      BASF_Sharedmem->get_pointer(BASF_Sharedmem->get_id(),
    //                  "TrkMgrSum",
    //                  & size);
  }
}

clear() [2/3]

void TTrackManager::clear ( void )

virtual

clears all internal information.

Reimplemented from TUpdater.

clear() [3/3]

void TTrackManager::clear ( void )

virtual

clears all internal information.

Reimplemented from TUpdater.

clearTables() [1/3]

void TTrackManager::clearTables

(

void

)

const

clears tables.

Definition at line 2521 of file TTrackManager.cxx.

                                            {
  //    BsDelEnt(RECMDC_TRK, BBS_ID_ALL);
  //    BsDelEnt(RECMDC_TRK_ADD, BBS_ID_ALL);
  //    BsDelEnt(RECMDC_MCTRK, BBS_ID_ALL);
  //    BsDelEnt(RECMDC_MCTRK2HEP, BBS_ID_ALL);
  unsigned n = MdcRecTrkCol::getMdcRecTrkCol()->size();
  for ( unsigned i = 0; i < n; i++ ) delete &( *MdcRecTrkCol::getMdcRecTrkCol() )[i];
 
  n = MdcRecTrkAddCol::getMdcRecTrkAddCol()->size();
  for ( unsigned i = 0; i < n; i++ ) delete &( *MdcRecTrkAddCol::getMdcRecTrkAddCol() )[i];
 
  n = MdcRecMctrkCol::getMdcRecMctrkCol()->size();
  for ( unsigned i = 0; i < n; i++ ) delete &( *MdcRecMctrkCol::getMdcRecMctrkCol() )[i];
 
  n = MdcRecMctrk2hepCol::getMdcRecMctrk2hepCol()->size();
  for ( unsigned i = 0; i < n; i++ )
    delete &( *MdcRecMctrk2hepCol::getMdcRecMctrk2hepCol() )[i];
 
  //...Clear track association...
  //    unsigned n = BsCouTab(RECMDC_WIRHIT);
  n = MdcRecWirhitCol::getMdcRecWirhitCol()->size();
  for ( unsigned i = 0; i < n; i++ )
  {
    //  reccdc_wirhit & h = * (reccdc_wirhit *)
    //      BsGetEnt(RECMDC_WIRHIT, i + 1, BBS_No_Index);
    //  h.m_trk = 0;
    MdcRec_wirhit& h = ( *MdcRecWirhitCol::getMdcRecWirhitCol() )[i];
    h.trk            = 0;
  }
  //    n = BsCouTab(DATMDC_MCWIRHIT);
  n = MdcDatMcwirhitCol::getMdcDatMcwirhitCol()->size();
  for ( unsigned i = 0; i < n; i++ )
  {
    //  datcdc_mcwirhit & h = * (datcdc_mcwirhit *)
    //      BsGetEnt(DATMDC_MCWIRHIT, i + 1, BBS_No_Index);
    //  h.m_trk = 0;
    MdcDat_mcwirhit& h = ( *MdcDatMcwirhitCol::getMdcDatMcwirhitCol() )[i];
    h.trk              = 0;
  }
}

clearTables() [2/3]

void TTrackManager::clearTables

(

void

)

const

clears tables.

clearTables() [3/3]

void TTrackManager::clearTables

(

void

)

const

clears tables.

closest() [1/3]

TTrack * TTrackManager::closest	(	const AList< TTrack > &	tracks,
		const TMDCWireHit &	hit ) const

returns a track which is the closest to a hit.

Definition at line 229 of file TTrackManager.cxx.

                                                                                          {
 
  TMLink t;
  t.hit( &hit );
  unsigned n           = tracks.length();
  double   minDistance = MAXDOUBLE;
  TTrack*  minTrk      = NULL;
 
  //...Loop over all tracks...
  for ( unsigned i = 0; i < n; i++ )
  {
    TTrack& trk = *tracks[i];
    int     err = trk.approach( t );
    if ( err < 0 ) continue;
    if ( minDistance > t.distance() )
    {
      minDistance = t.distance();
      minTrk      = &trk;
    }
  }
 
  return minTrk;
}

Referenced by salvage().

closest() [2/3]

TTrack * TTrackManager::closest	(	const AList< TTrack > &	,
		const TMDCWireHit &	) const

returns a track which is the closest to a hit.

closest() [3/3]

TTrack * TTrackManager::closest	(	const AList< TTrack > &	,
		const TMDCWireHit &	) const

returns a track which is the closest to a hit.

debugLevel() [1/6]

int TTrackManager::debugLevel ( int a )

inline

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

{ return _debugLevel = a; }

debugLevel() [2/6]

int TTrackManager::debugLevel

(

int

)

debugLevel() [3/6]

int TTrackManager::debugLevel

(

int

)

debugLevel() [4/6]

int TTrackManager::debugLevel ( void ) const

inline

returns/sets debug level.

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

{ return _debugLevel; }

debugLevel() [5/6]

int TTrackManager::debugLevel

(

void

)

const

returns/sets debug level.

debugLevel() [6/6]

int TTrackManager::debugLevel

(

void

)

const

returns/sets debug level.

determineIP() [1/3]

void TTrackManager::determineIP

(

void

)

determines IP.

Definition at line 2624 of file TTrackManager.cxx.

                                      {
  static const HepVector3D InitialVertex( 0., 0., 0. );
 
  //...Track selection...
  //    unsigned n = BsCouTab(RECTRK);
  unsigned             n = MdcTrkCol::getMdcTrkCol()->size();
  AList<MdcTrk_localz> zList;
  for ( unsigned i = 0; i < n; i++ )
  {
    //  const rectrk & t = * (rectrk *) BsGetEnt(RECTRK, i + 1, BBS_No_Index);
    const MdcTrk& t = ( *MdcTrkCol::getMdcTrkCol() )[i];
 
    if ( t.prekal == 0 ) continue;
    //  const reccdc_trk_add & c = * (reccdc_trk_add *)
    //      BsGetEnt(RECMDC_TRK_ADD, t.m_prekal, BBS_No_Index);
    const MdcRec_trk_add& c = *t.prekal->add;
 
    //...Only good tracks...
    if ( c.quality ) continue;
 
    //...Require SVD hits...
    //  const rectrk_global & g = * (rectrk_global *) BsGetEnt(RECTRK_GLOBAL,
    //                                 t.m_glob[2],
    //                                 BBS_No_Index);
    const MdcTrk_global& g = *t.glob[2];
 
    if ( !&g ) continue;
    if ( g.nhits[3] < 2 ) continue;
    if ( g.nhits[4] < 2 ) continue;
 
    //...OK...
    //  const rectrk_localz & z = * (rectrk_localz *) BsGetEnt(RECTRK_LOCALZ,
    //                                 t.m_zero[2],
    //                                 BBS_No_Index);
    MdcTrk_localz& z = *t.zero[2];
 
    if ( !&z ) continue;
    //  zList.append((rectrk_localz &) z);
    zList.append( z );
  }
  unsigned nZ = zList.length();
  if ( nZ < 2 ) return;
 
  //...Fitting...
  //      kvertexfitter kvf;
  //      kvf.initialVertex(initialVertex);
  //      for (unsigned i = 0; i < nZ; i++) {
  //    kvf.addTrack();
  //      }
}

determineIP() [2/3]

void TTrackManager::determineIP

(

void

)

determines IP.

determineIP() [3/3]

void TTrackManager::determineIP

(

void

)

determines IP.

determineT0() [1/3]

StatusCode TTrackManager::determineT0	(	unsigned	level,
		unsigned	nMaxTracks )

determines T0 and refit all tracks.

Definition at line 1272 of file TTrackManager.cxx.

                                                                     {
#ifdef TRKRECO_DEBUG_DETAIL
  if ( level == 0 )
  {
    std::cout << "TTrackManager::determineT0 !!! called with level = 0";
    std::cout << std::endl;
  }
#endif
 
  IMessageSvc* msgSvc;
  Gaudi::svcLocator()->service( "MessageSvc", msgSvc );
  MsgStream log( msgSvc, "TTrackManager" );
 
  IDataProviderSvc* eventSvc = NULL;
  Gaudi::svcLocator()->service( "EventDataSvc", eventSvc );
 
  static bool     first   = true;
  static unsigned methode = 0;
  if ( first )
  {
    first = false;
 
    if ( level == 1 ) { _cFitter.fit2D( true ); }
    else if ( level == 2 )
    {
      // default setting
    }
    else if ( level == 3 )
    {
      //        _cFitter.sag(true);  //Liuqg
    }
    else if ( level == 4 )
    {
      //        _cFitter.sag(true);  //Liuqg
      _cFitter.propagation( true );
    }
    else if ( level == 5 )
    {
      //        _cFitter.sag(true);  //Liuqg
      _cFitter.propagation( true );
      _cFitter.tof( true );
    }
    else if ( level == 6 )
    {
      methode = 1;
      //        _cFitter.sag(true); //Liuqg
      _cFitter.propagation( true );
      _cFitter.tof( true );
    }
    else if ( level == 7 )
    {
      methode = 2;
      //        _cFitter.sag(true);  //Liuqg
      _cFitter.propagation( true );
      _cFitter.tof( true );
    }
  }
 
  unsigned n = _tracks.length();
  if ( !n ) return StatusCode::SUCCESS;
 
  if ( nMax == 0 ) nMax = n;
  if ( n > nMax ) n = nMax;
 
  //    float t0 = 0.;
  _t0 = 999.;
 
  // read t0 from TDS
  float                      t0_sta = -1;
  float                      tev    = 0;
  SmartDataPtr<RecEsTimeCol> aevtimeCol( eventSvc, "/Event/Recon/RecEsTimeCol" );
  if ( aevtimeCol )
  {
    RecEsTimeCol::iterator iter_evt = aevtimeCol->begin();
    t0_sta                          = ( *iter_evt )->getStat();
    tev                             = ( *iter_evt )->getTest();
    //  cout<<"t0_sta: "<<t0_sta<<"  tev: "<<tev<<endl;
  }
  else { log << MSG::WARNING << "Could not find RecEsTimeCol" << endmsg; }
 
  if ( methode == 0 ) _t0 = T0( n );
 
  else if ( methode == 1 ) _t0 = T0Fit( n );
 
  else if ( methode == 2 )
  { // revise method == 2 to BESIII environment. Liuqg
    if ( t0_sta != 1 )
    { // 1: tof  11:tof after reco  //no tof
      _t0 = T0Fit( n );
      // cout<<"t0: "<<_t0<<endl;
    }
  }
 
  // std::cout << "reccdc_timing=" << BsCouTab(RECMDC_TIMING) << std::endl;
  /*    else if (methode == 2 && BsCouTab(RECMDC_TIMING) != 0) {
        struct reccdc_timing * r0 = (struct reccdc_timing *)
        BsGetEnt(RECMDC_TIMING, BsCouTab(RECMDC_TIMING), BBS_No_Index);
        if (r0->m_quality == 102) {
        if (BsCouTab(BELLE_EVENT)) {
        struct belle_event * b0 = (struct belle_event *)
        BsGetEnt(BELLE_EVENT, 1, BBS_No_Index);
        if(1==b0->m_ExpMC)                    t0 = T0Fit(n);
        if(2==b0->m_ExpMC && r0->m_time !=0.) t0 = T0Fit(n);
        }
        }
        else if (r0->m_quality == 100) t0 = T0Fit(n);
  //      std::cout << "quality=" << r0->m_quality << std::endl;
  }
  */
 
  //...For debug...
  if ( _debugLevel )
  {
    std::cout << "TTrackManager::determineT0 ... methode=" << methode;
    std::cout << ", T0 offset=" << -_t0;
    std::cout << ", # of tracks used=" << n << std::endl;
  }
 
  //...store them... Liuqg
  float t0_rec    = 999.;
  int   t0_recSta = 8;
  if ( fabs( _t0 + tev ) < 4 ) t0_rec = 0;
  if ( fabs( _t0 + tev - 8 ) < 4 ) t0_rec = 8;
  if ( fabs( _t0 + tev - 16 ) < 4 ) t0_rec = 16;
  log << MSG::INFO << "beginning to make RecEsTimeCol" << endmsg;
  IDataManagerSvc* dataManSvc = dynamic_cast<IDataManagerSvc*>( eventSvc );
  DataObject*      aEvTime;
  eventSvc->findObject( "/Event/Recon/RecEsTimeCol", aEvTime );
  if ( aEvTime != NULL && t0_rec < 25 )
  {
    dataManSvc->clearSubTree( "/Event/Recon/RecEsTimeCol" );
    eventSvc->unregisterObject( "/Event/Recon/RecEsTimeCol" );
  }
  if ( t0_rec < 25 )
  {
 
    RecEsTimeCol* aEvTimeCol = new RecEsTimeCol;
    RecEsTime*    aevtime    = new RecEsTime;
    aevtime->setTest( t0_rec );
    aevtime->setStat( t0_recSta );
    aEvTimeCol->push_back( aevtime );
 
    // register event time  to TDS
    // check whether the t0  has been already registered
    StatusCode evtime = eventSvc->registerObject( "/Event/Recon/RecEsTimeCol", aEvTimeCol );
 
    if ( evtime != StatusCode::SUCCESS )
    {
      log << MSG::FATAL << "Could not register Event Start Time" << endmsg;
      return ( StatusCode::FAILURE );
    }
  }
  return ( StatusCode::SUCCESS );
}

determineT0() [2/3]

StatusCode TTrackManager::determineT0	(	unsigned	level,
		unsigned	nMaxTracks )

determines T0 and refit all tracks.

determineT0() [3/3]

StatusCode TTrackManager::determineT0	(	unsigned	level,
		unsigned	nMaxTracks )

determines T0 and refit all tracks.

divide() [1/3]

TMLink & TTrackManager::divide	(	const TTrack &	t,
		AList< TMLink > *	l ) const

Definition at line 1006 of file TTrackManager.cxx.

                                                                       {
  TMLink&           start  = *OuterMost( t.links() );
  const HepPoint3D& center = t.helix().center();
  const HepVector3D a      = start.positionOnTrack() - center;
  for ( unsigned j = 0; j < t.nLinks(); j++ )
  {
    if ( t[j] == &start ) continue;
    TMLink&           k = *t[j];
    const HepVector3D b = k.positionOnTrack() - center;
    if ( a.cross( b ).z() >= 0. ) l[0].append( k );
    else l[1].append( k );
  }
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    outer link = " << start.hit()->wire()->name() << std::endl;
  std::cout << "        nLinks of 0 = " << l[0].length() << std::endl;
  std::cout << "        nLinks of 1 = " << l[1].length() << std::endl;
#endif
 
  if ( l[0].length() == 0 || l[1].length() == 0 ) return divideByIp( t, l );
 
  return start;
}

divide() [2/3]

TMLink & TTrackManager::divide	(	const TTrack &	t,
		AList< TMLink > *	l ) const

divide() [3/3]

TMLink & TTrackManager::divide	(	const TTrack &	t,
		AList< TMLink > *	l ) const

divideByIp() [1/3]

TMLink & TTrackManager::divideByIp	(	const TTrack &	t,
		AList< TMLink > *	l ) const

Definition at line 1030 of file TTrackManager.cxx.

                                                                           {
  l[0].removeAll();
  l[1].removeAll();
 
  const HepPoint3D& center = t.helix().center();
  const HepVector3D a      = ORIGIN - center;
  for ( unsigned j = 0; j < t.nLinks(); j++ )
  {
    TMLink&           k = *t[j];
    const HepVector3D b = k.positionOnTrack() - center;
    if ( a.cross( b ).z() >= 0. ) l[0].append( k );
    else l[1].append( k );
  }
 
  //...This is a dummy...
  TMLink& start = *OuterMost( t.links() );
  return start;
}

Referenced by divide(), maskCurl(), and maskNormal().

divideByIp() [2/3]

TMLink & TTrackManager::divideByIp	(	const TTrack &	t,
		AList< TMLink > *	l ) const

divideByIp() [3/3]

TMLink & TTrackManager::divideByIp	(	const TTrack &	t,
		AList< TMLink > *	l ) const

dump() [1/3]

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

dumps debug information.

Definition at line 90 of file TTrackManager.cxx.

                                                                           {
  bool def = ( msg == "" ) ? true : false;
  /*
      if (msg.find("summary") != std::string::npos || msg.find("detail") != std::string::npos)
     { struct summary s;
          //        bzero((char*)& s, sizeof(struct summary));
          memset((char*)& s, 0, sizeof(struct summary));
          for (int i = 0; i < BASF_Sharedmem->nprocess(); i++) {
              int size;
              struct summary & r = * (struct summary *)
                  BASF_Sharedmem->get_pointer(i, "TrkMgrSum", & size);
              s._nEvents += r._nEvents;
              s._nTracks += r._nTracks;
              s._nTracksAll += r._nTracksAll;
              s._nTracks2D += r._nTracks2D;
              s._nTracksFinal += r._nTracksFinal;
              s._nSuperMoms += r._nSuperMoms;
              s._nToBeMerged += r._nToBeMerged;
              s._nToBeMergedMoreThanTwo += r._nToBeMergedMoreThanTwo;
          }
 
          std::cout << pre;
          std::cout << "all events : " << s._nEvents << std::endl;
          std::cout << pre;
          std::cout << "all tracks : " << s._nTracksAll << std::endl;
          std::cout << pre;
          std::cout << "    good tracks  : " << s._nTracks << std::endl;
          std::cout << pre;
          std::cout << "    2D tracks    : " << s._nTracks2D << std::endl;
          std::cout << pre;
          std::cout << "    final tracks : " << s._nTracksFinal << std::endl;
          std::cout << pre;
          std::cout << "    super mom.   : " << s._nSuperMoms << std::endl;
          std::cout << pre;
          std::cout << "    to be mreged : " << s._nToBeMerged << std::endl;
          std::cout << pre;
          std::cout << "    to be mreged2: " << s._nToBeMergedMoreThanTwo
                    << std::endl;
      }
  */
  if ( def || msg.find( "eventSummary" ) != std::string::npos ||
       msg.find( "detail" ) != std::string::npos )
  {
    std::cout << pre << "tracks reconstructed : " << _tracksAll.length() << std::endl;
    std::cout << pre << "good tracks          : " << _tracks.length() << std::endl;
    std::cout << pre << "2D tracks            : " << _tracks2D.length() << std::endl;
    std::cout << pre << "Track list:" << std::endl;
 
    std::string tab = pre;
    std::string spc = " ";
    for ( unsigned i = 0; i < _tracksAll.length(); i++ )
    {
      std::cout << tab << TrackDump( *_tracksAll[i] ) << std::endl;
      if ( msg.find( "hits" ) != std::string::npos ||
           msg.find( "detail" ) != std::string::npos )
        Dump( _tracksAll[i]->links(), "hits sort flag" );
    }
  }
}

dump() [2/3]

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

dumps debug information.

dump() [3/3]

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

dumps debug information.

finish() [1/3]

void TTrackManager::finish

(

void

)

finishes tracks.

Definition at line 844 of file TTrackManager.cxx.

                                 {
  refit();
  movePivot();
  if ( _debugLevel > 1 )
  {
    std::cout << name() << " ... finishing" << std::endl;
    //  unsigned n = _tracksAll.length();
    unsigned n = _tracks.length();
    for ( unsigned i = 0; i < n; i++ )
    {
      //        TTrack & t = * _tracksAll[i];
      TTrack& t = *_tracks[i];
      std::cout << "    " << t.name() << std::endl;
      t.dump( "hits mc track flag sort", "    " );
    }
  }
}

finish() [2/3]

void TTrackManager::finish

(

void

)

finishes tracks.

finish() [3/3]

void TTrackManager::finish

(

void

)

finishes tracks.

fittingFlag() [1/3]

void TTrackManager::fittingFlag ( unsigned a )

inline

sets fitting flag.

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

                                                   {
  if ( a & 1 ) _fitter.sag( true );
  if ( a & 2 ) _fitter.propagation( true );
  if ( a & 4 ) _fitter.tof( true );
  if ( a & 8 ) _fitter.freeT0( true );
}

fittingFlag() [2/3]

void TTrackManager::fittingFlag

(

unsigned

)

sets fitting flag.

fittingFlag() [3/3]

void TTrackManager::fittingFlag

(

unsigned

)

sets fitting flag.

goodTrack() [1/3]

bool TTrackManager::goodTrack ( const TTrack & t, bool track2D = false )

static

checks goodness of a track.

Definition at line 3025 of file TTrackManager.cxx.

                                                             {
 
  //...Check number of hits...
  if ( !checkNumberOfHits( t, track2D ) ) return false;
 
  //...Check helix parameter...
  if ( HelixHasNan( t.helix() ) ) return false;
 
  return true;
}

goodTrack() [2/3]

bool TTrackManager::goodTrack ( const TTrack & , bool track2D = false )

static

checks goodness of a track.

goodTrack() [3/3]

bool TTrackManager::goodTrack ( const TTrack & , bool track2D = false )

static

checks goodness of a track.

makeTds() [1/3]

StatusCode TTrackManager::makeTds	(	RecMdcTrackCol *	trackList,
		RecMdcHitCol *	hitList,
		int	tkStat = 1,
		int	runge = 0,
		int	cal = 0 )

stores track info. into TDS. by Zang Shilei

added by Jike Wang

add by jialk: returns driftTime prop time correction & entra angle

check the result:MdcTrackCol

Definition at line 255 of file TTrackManager.cxx.

                                  { // yzhang change interface 2010-05-14
  IMessageSvc* msgSvc;
  Gaudi::svcLocator()->service( "MessageSvc", msgSvc );
 
  MsgStream log( msgSvc, "TrkReco" );
 
  IDataProviderSvc* eventSvc = NULL;
  Gaudi::svcLocator()->service( "EventDataSvc", eventSvc );
 
#ifdef TRKRECO_DEBUG
  if ( eventSvc ) { log << MSG::INFO << "makeTds:event Svc has been found" << endmsg; }
  else
  {
    log << MSG::FATAL << "makeTds:Could not find eventSvc" << endmsg;
    return StatusCode::FAILURE;
  }
#endif
  // check whether  Recon  already registered
  //    DataObject *aReconEvent;
  //    eventSvc->findObject("/Event/Recon",aReconEvent);
  //    if(!aReconEvent) {
  //       ReconEvent* recevt = new ReconEvent;
  //       StatusCode sc = eventSvc->registerObject("/Event/Recon",recevt );
  //       if(sc!=StatusCode::SUCCESS) {
  //        log << MSG::FATAL << "Could not register ReconEvent" <<endmsg;
  //        return( StatusCode::FAILURE);
  //      }
  //   }
  //
  //   /// Unregister Tracks
  //   IDataManagerSvc *dataManSvc;
  //   if(doClear){//yzhang add, do NOT clear Tds when associate rec
  //     dataManSvc= dynamic_cast<IDataManagerSvc*>(eventSvc);
  //     DataObject *aTrackCol;
  //     eventSvc->findObject("/Event/Recon/RecMdcTrackCol",aTrackCol);
  //     if(aTrackCol != NULL) {
  //       dataManSvc->clearSubTree("/Event/Recon/RecMdcTrackCol");
  //       eventSvc->unregisterObject("/Event/Recon/RecMdcTrackCol");
  //     }
  //     DataObject *aRecHitCol;
  //     eventSvc->findObject("/Event/Recon/RecMdcHitCol",aRecHitCol);
  //     if(aRecHitCol != NULL) {
  //       dataManSvc->clearSubTree("/Event/Recon/RecMdcHitCol");
  //       eventSvc->unregisterObject("/Event/Recon/RecMdcHitCol");
  //     }
  //   }
  //   /// Writing
  //   RecMdcTrackCol* trkcol = new RecMdcTrackCol;
  //   RecMdcHitCol* hitcol = new RecMdcHitCol;
 
  unsigned n      = _tracks.length();
  int      nTdsTk = trkcol->size();
  for ( int i = 0; i < n; i++ )
  {
    RecMdcTrack* trk        = new RecMdcTrack;
    TTrack*      t          = _tracks[i];
    int          trackindex = i + nTdsTk; // for combination
 
    HitRefVec     hitrefvec;
    AList<TMLink> hits  = t->links();
    float         chisq = 0;
 
    HepPoint3D pos    = t->helix().pivot();
    int        charge = -1;
    HepVector  m_a( 5, 0 );
    m_a            = t->helix().a();
    int statfinder = t->getFinderType();
    if ( abs( statfinder ) > 1000 && brunge ) statfinder = 103;
    if ( abs( statfinder ) > 1000 && !brunge ) statfinder = 3;
    // be cautious
    if ( !brunge ) m_a[2] = -m_a[2];
    if ( abs( m_a[1] ) > 999 ) m_a[1] = 0;
    while ( m_a[1] < 0 ) { m_a[1] += 2 * pi; }
    while ( m_a[1] > 2 * pi ) { m_a[1] -= 2 * pi; }
    /// added by Jike Wang
    const double x0 = t->helix().pivot().x();
    const double y0 = t->helix().pivot().y();
 
    const double dr    = m_a[0];
    const double phi0  = m_a[1];
    const double kappa = m_a[2];
    const double dz    = m_a[3];
    const double tanl  = m_a[4];
 
    const double Bz    = -1000 * ( getPmgnIMF()->getReferField() );
    const double alpha = 333.564095 / Bz;
 
    const double cox = x0 + dr * cos( phi0 ) - alpha * cos( phi0 ) / kappa;
    const double coy = y0 + dr * sin( phi0 ) - alpha * sin( phi0 ) / kappa;
 
    unsigned nHits    = t->links().length();
    unsigned nStereos = 0;
    unsigned firstLyr = 44;
    unsigned lastLyr  = 0;
    for ( unsigned j = 0; j < nHits; j++ )
    {
 
      TMLink* l = hits[j];
 
      HepPoint3D ontrack = l->positionOnTrack();
      HepPoint3D onwire  = l->positionOnWire();
      HepPoint3D dir( ontrack.y() - coy, cox - ontrack.x(), 0 );
      double     pos_phi = onwire.phi();
      double     dir_phi = dir.phi();
      while ( pos_phi > pi ) { pos_phi -= pi; }
      while ( pos_phi < 0 ) { pos_phi += pi; }
      while ( dir_phi > pi ) { dir_phi -= pi; }
      while ( dir_phi < 0 ) { dir_phi += pi; }
      double entrangle = dir_phi - pos_phi;
      while ( entrangle > pi / 2 ) entrangle -= pi;
      while ( entrangle < ( -1 ) * pi / 2 ) entrangle += pi;
 
      // jialk setFltLen to tds
      int    imass   = 3;
      float  tl      = t->helix().a()[4];
      float  f       = sqrt( 1. + tl * tl );
      float  s       = fabs( t->helix().curv() ) * fabs( l->dPhi() ) * f;
      float  p       = f / fabs( t->helix().a()[2] );
      float  Mass[5] = { 0.000511, 0.105, 0.140, 0.493677, 0.98327 };
      double tof = s / 29.98 * sqrt( 1.0 + ( Mass[imass - 1] / p ) * ( Mass[imass - 1] / p ) );
 
      RecMdcHit* hit = new RecMdcHit;
      hit->setId( l->hit()->reccdc()->id );
      //       hit->setTrkId(i);
      hit->setTrkId( trackindex ); // for combination
      hit->setDriftDistLeft( l->drift( 0 ) );
      hit->setDriftDistRight( l->drift( 1 ) );
      hit->setErrDriftDistLeft( l->dDrift( 0 ) );
      hit->setErrDriftDistRight( l->dDrift( 1 ) );
      hit->setChisqAdd( l->pull() );
      hit->setFlagLR( l->leftRight() );
      //  hit->setStat(l->hit()->state());
      hit->setStat( 1 );
      //  std::cout<<"state :"<< l->hit()->state() << std::endl;
      hit->setAdc( l->hit()->reccdc()->adc );
      //           hit->setTdc((l->hit()->reccdc()->tdc + t0bes)/0.09375); //jialk
      hit->setTdc( l->hit()->reccdc()->timechannel );
      hit->setFltLen( tof * 30 ); // jialk
      // std::cout<<"tdc :"<< l->hit()->reccdc()->tdc <<"  "<<"t0bes : "<< t0bes <<std::endl;
      if ( cal ) { hit->setDoca( l->distancenew() ); }
      else { hit->setDoca( l->distance() ); }
      hit->setZhit( l->positionOnTrack().z() );
 
      /// add by jialk: returns driftTime prop time correction & entra angle
      hit->setEntra( entrangle );
      hit->setDriftT( l->getDriftTime() );
 
      Identifier hitIdf = MdcID::wire_id( l->wire()->layerId(), l->wire()->localId() );
      hit->setMdcId( hitIdf );
      hitcol->push_back( hit );
 
      SmartRef<RecMdcHit> refhit( hit );
      hitrefvec.push_back( refhit );
      chisq += l->pull();
      if ( l->wire()->stereo() ) ++nStereos;
      if ( firstLyr > l->wire()->layerId() ) firstLyr = l->wire()->layerId();
      if ( lastLyr < l->wire()->layerId() ) lastLyr = l->wire()->layerId();
    }
 
    HepSymMatrix m_ea = t->helix().Ea();
    double       errorMat[15];
    int          k = 0;
    for ( int ie = 0; ie < 5; ie++ )
    {
      for ( int je = ie; je < 5; je++ )
      {
        errorMat[k] = m_ea[ie][je];
        k++;
      }
    }
 
    trk->setTrackId( trackindex );
    trk->setHelix( m_a );
    trk->setPxy( t->pt() );
    trk->setPx( t->pt() * ( -sin( t->helix().phi0() ) ) );
    trk->setPy( t->pt() * cos( t->helix().phi0() ) );
    trk->setPz( t->pz() );
    trk->setP( t->ptot() );
 
    // jialk
    double theta = acos( ( t->pz() ) / t->ptot() );
    trk->setTheta( theta );
 
    double poca_dr   = t->helix().dr();
    double poca_phi0 = t->helix().phi0();
    trk->setPhi( poca_phi0 );
    HepPoint3D poca( pos.x() + poca_dr * cos( poca_phi0 ),
                     pos.y() + poca_dr * sin( poca_phi0 ), pos.z() + t->helix().dz() );
    trk->setPoca( poca );
    trk->setX( poca.x() );
    trk->setY( poca.y() );
    trk->setZ( poca.z() );
    trk->setR( sqrt( poca.x() * poca.x() + poca.y() * poca.y() ) );
    trk->setPivot( pos );
    trk->setVX0( pos.x() );
    trk->setVY0( pos.y() );
    trk->setVZ0( pos.z() );
 
    trk->setError( m_ea );
    //    trk->setError(errorMat);   //...............2
 
    //    double poca_dr = t->helix().dr();
    //    double poca_phi0 = t->helix().phi0();
    //    HepPoint3D poca(poca_dr*cos(poca_phi0),poca_dr*sin(poca_phi0),t->helix().dz());
    //    trk->setVX0(pos.x()+poca_dr*cos(poca_phi0));
    //    trk->setVY0(pos.y()+poca_dr*sin(poca_phi0));
    //    trk->setVZ0(pos.z()+t->helix().dz());
    //    cout<<"pivot:("<<pos.x()<<","<<pos.y()<<","<<pos.z()<<")"<<endl;
    //    cout<<"poca:("<<trk->getVX0()<<","<<trk->getVY0()<<","<<trk->getVZ0()<<")"<<endl;
 
    trk->setChi2( chisq );
    trk->setNdof( nHits - 5 );
    if ( t->quality() & TrackQuality2D ) trk->setNdof( nHits - 3 );
 
    TMLink* last = OuterMost( hits );
    t->approach( *last );
    trk->setFiTerm( last->dPhi() );
 
    trk->setNhits( nHits );
    trk->setNster( nStereos );
    trk->setStat( statfinder ); // yzhang add stat: ConformalFinder=2, CurlFiner=3
    // xuqn 2013-02-26
    // trk->setCharge(int(t->charge())*(-1));
    if ( !brunge ) trk->setCharge( int( t->charge() ) * ( -1 ) );
    else trk->setCharge( t->charge() );
    trk->setVecHits( hitrefvec );
    trk->setFirstLayer( firstLyr );
    trk->setLastLayer( lastLyr );
    // cout<<"first: "<<firstLyr<<"  last: "<<lastLyr<<endl;
    trkcol->push_back( trk );
  }
 
  // StatusCode trksc;
  // trksc = eventSvc->registerObject("/Event/Recon/RecMdcTrackCol", trkcol);
  // if( trksc.isFailure() ) {
  //   log << MSG::FATAL << "Could not register MdcTrack" << endmsg;
  //   return trksc;
  // }
  //
  // StatusCode hitsc = eventSvc->registerObject("/Event/Recon/RecMdcHitCol", hitcol);
  // if ( hitsc.isFailure() ) {
  //   log << MSG::FATAL << "Could not register MdcRecHit" << endmsg;
  //   return hitsc;
  // }
  // log << MSG::INFO << "MdcTrackCol registered successfully!" <<endmsg;
 
  /// check the result:MdcTrackCol
#ifdef TRKRECO_DEBUG
  SmartDataPtr<RecMdcTrackCol> newtrkCol( eventSvc, "/Event/Recon/RecMdcTrackCol" );
  if ( !newtrkCol )
  {
    log << MSG::FATAL << "Could not find RecMdcTrackCol" << endmsg;
    return ( StatusCode::FAILURE );
  }
  log << MSG::INFO << "Begin to check RecMdcTrackCol" << endmsg;
  std::cout << std::fixed << std::setprecision( 6 ); // yzhang debug for 720
  RecMdcTrackCol::iterator iter_trk = newtrkCol->begin();
  for ( ; iter_trk != newtrkCol->end(); iter_trk++ )
  {
    const RecMdcTrack* tk = *iter_trk;
    std::cout << "//==== " << name() << " print RecMdcTrack No." << tk->trackId() << " :"
              << std::endl;
    cout << " dr " << tk->helix( 0 ) << " phi0 " << tk->helix( 1 ) << " cpa " << tk->helix( 2 )
         << " z0 " << tk->helix( 3 ) << " tanl " << tk->helix( 4 ) << endl;
    bool printTrackDetail = true;
    if ( printTrackDetail )
    {
      std::cout << " q " << tk->charge() << " theta " << tk->theta() << " phi " << tk->phi()
                << " x0 " << tk->x() << " y0 " << tk->y() << " z0 " << tk->z() << " r0 "
                << tk->r() << endl;
      std::cout << " p " << tk->p() << " pt " << tk->pxy() << " pxyz(" << tk->px() << ","
                << tk->py() << "," << tk->pz() << ")" << endl;
      std::cout << " tkStat " << tk->stat() << " chi2 " << tk->chi2() << " ndof " << tk->ndof()
                << " nhit " << tk->getNhits() << " nst " << tk->nster() << endl;
      bool printErrMat = true;
      if ( printErrMat )
      {
        std::cout << "errmat   " << std::endl;
        for ( int i = 0; i < 15; i++ ) { std::cout << " " << tk->err( i ); }
        std::cout << "   " << std::endl;
      }
    }
 
    bool printHit = true;
    if ( printHit )
    {
      int  haveDigi[43][288];
      bool printMcTk = true;
      if ( printMcTk )
      {
        for ( int ii = 0; ii < 43; ii++ )
        {
          for ( int jj = 0; jj < 43; jj++ ) { haveDigi[ii][jj] = -9999; }
        }
        MdcDigiVec           mdcDigiVec = getRawDataProviderSvc()->getMdcDigiVec();
        MdcDigiCol::iterator iter       = mdcDigiVec.begin();
        std::cout << name() << "//==== " << name()
                  << " print MdcDigiVec, nDigi=" << mdcDigiVec.size() << " :" << std::endl;
        for ( int iDigi = 0; iter != mdcDigiVec.end(); iter++, iDigi++ )
        {
          long l         = MdcID::layer( ( *iter )->identify() );
          long w         = MdcID::wire( ( *iter )->identify() );
          haveDigi[l][w] = ( *iter )->getTrackIndex();
        }
      }
 
      int nhits = tk->getVecHits().size();
      std::cout << "nHits=" << nhits << std::endl;
      cout << "No. ";
      if ( printMcTk ) cout << "mcTkId ";
      cout << "(layer,wire,lr) stat z" << endl;
      for ( int ii = 0; ii < nhits; ii++ )
      {
        Identifier id( tk->getVecHits()[ii]->getMdcId() );
        int        layer = MdcID::layer( id );
        int        wire  = MdcID::wire( id );
        cout << ii << ":";
        if ( printMcTk ) { cout << haveDigi[layer][wire]; }
        cout << " (" << layer << "," << wire << "," << tk->getVecHits()[ii]->getFlagLR()
             << ")  " << tk->getVecHits()[ii]->getStat() << "  "
             << tk->getVecHits()[ii]->getZhit() << " " << std::endl;
      } // end of hit list
      std::cout << "  " << std::endl;
    }
    /*
       std::cout << "retrieved MDC track:"
       << "Track Id: " << (*iter_trk)->getTrkId()
       << "   Pivot: " << (*iter_trk)->getVX0() << " "
       << (*iter_trk)->getVY0() << " " << (*iter_trk)->getVZ0()
       << endl
       << "Phi0: " << (*iter_trk)->getFi0() << "   Error of Phi0 "
       << (*iter_trk)->getError()(2,2) << endmsg
       << "kappa: " << (*iter_trk)->getCpa() << "   Error of kappa "
       << (*iter_trk)->getError()(3,3) << endmsg
       << "Tanl: " << (*iter_trk)->getTanl() << "   Error of Tanl "
       << (*iter_trk)->getError()(5,5) << endmsg
       << "Chisq of fit: "<< (*iter_trk)->getChisq()
       << "   Phi terminal: "<< (*iter_trk)->getFiTerm()
       << endl
       << "Number of hits: "<< (*iter_trk)->getNhits()
       << "   Number of stereo hits " << (*iter_trk)->getNster()
       << endmsg;
 
       log << MSG::DEBUG << "hitList of this  track:" << endmsg;
       std::vector<MdcRecHit*> gothits = (*iter_trk)->getVecHits();
       std::vector<MdcRecHit*>::iterator it_gothit = gothits.begin();
       for( ; it_gothit != gothits.end(); it_gothit++){
       log << MSG::DEBUG << "hits Id: "<<(*it_gothit)->getId()
       << "   hits DDL&DDR " <<(*it_gothit)->getDriftDistLeft()
       << "   hits MDC IDentifier " <<(*it_gothit)->getMdcId()
       << endmsg
       << "   hits TDC " <<(*it_gothit)->getTdc()
       << "   hits ADC " <<(*it_gothit)->getAdc() << endmsg;
       }
       */
  }
  std::cout << std::defaultfloat;
#endif
  return StatusCode::SUCCESS;
}

makeTds() [2/3]

StatusCode TTrackManager::makeTds	(	RecMdcTrackCol *	trackList,
		RecMdcHitCol *	hitList,
		int	tkStat = 1,
		int	runge = 0,
		int	cal = 0 )

stores track info. into TDS. by Zang Shilei

makeTds() [3/3]

StatusCode TTrackManager::makeTds	(	RecMdcTrackCol *	trackList,
		RecMdcHitCol *	hitList,
		int	tkStat = 1,
		int	runge = 0,
		int	cal = 0 )

stores track info. into TDS. by Zang Shilei

mask() [1/3]

void TTrackManager::mask

(

void

)

const

masks hits out which are in tail of curly tracks.

Definition at line 921 of file TTrackManager.cxx.

                                     {
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << name() << " ... masking" << std::endl;
#endif
 
  unsigned n = _tracks.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    TTrack& t = *_tracks[i];
 
    //...Skip if no core...
    //   This should not be happend...
    if ( !t.cores().length() ) continue;
 
    //...Counts # of hits per layer...
    unsigned nHits[43];
    NHits( t.cores(), nHits );
 
    //...Check each layer...
    bool needMask = false;
    for ( unsigned j = 0; j < 43; j++ )
    {
      if ( nHits[j] > 1 )
      {
        AList<TMLink> linksInLayer = SameLayer( t.links(), j );
        if ( Width( linksInLayer ) > 2 )
        {
          needMask = true;
 
#ifdef TRKRECO_DEBUG_DETAIL
          Dump( linksInLayer, "sort", "    -->" );
#endif
          break;
        }
      }
    }
    if ( !needMask ) continue;
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    trk" << i << "(id is tmp) needs mask" << std::endl;
    std::cout << "        type = " << t.type() << std::endl;
#endif
 
    //...Switch by track type...
    switch ( t.type() )
    {
    case TrackTypeNormal:
      maskNormal( t );
      maskMultiHits( t );
      break;
    case TrackTypeCurl:
      maskCurl( t );
      maskMultiHits( t );
      break;
    default: break;
    }
 
    //...Refit...
    // refit() ???
    _fitter.fit( t );
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    masking result : ";
    t.dump( "detail sort", "    " );
#endif
  }
}

mask() [2/3]

void TTrackManager::mask

(

void

)

const

masks hits out which are in tail of curly tracks.

mask() [3/3]

void TTrackManager::mask

(

void

)

const

masks hits out which are in tail of curly tracks.

maskBadHits() [1/3]

void TTrackManager::maskBadHits ( const AList< TTrack > & tracks, float maxSigma2 )

static

masks hits with large chisq as associated hits. Pull in TMLink is used.

Definition at line 2492 of file TTrackManager.cxx.

                                                                              {
#ifdef TRKRECO_DEBUG
  std::cout << "... trkmgr::maskBadHits" << std::endl;
#endif
 
  unsigned n = tracks.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    TTrack&             t           = *tracks[i];
    bool                toBeUpdated = false;
    const AList<TMLink> links       = t.links();
    unsigned            nHits       = links.length();
    for ( unsigned j = 0; j < nHits; j++ )
    {
      if ( links[j]->pull() > maxSigma2 )
      {
        links[j]->hit()->state( links[j]->hit()->state() | WireHitInvalidForFit );
        toBeUpdated = true;
#ifdef TRKRECO_DEBUG
        std::cout << "    " << t.name() << " : ";
        std::cout << links[j]->wire()->name() << "(pull=";
        std::cout << links[j]->pull() << ") is masked" << std::endl;
#endif
      }
    }
    if ( toBeUpdated ) t.update();
  }
}

Referenced by TConformalFinder::doit().

maskBadHits() [2/3]

void TTrackManager::maskBadHits ( const AList< TTrack > & , float maxSigma2 )

static

masks hits with large chisq as associated hits. Pull in TMLink is used.

maskBadHits() [3/3]

void TTrackManager::maskBadHits ( const AList< TTrack > & , float maxSigma2 )

static

masks hits with large chisq as associated hits. Pull in TMLink is used.

maskCurl() [1/3]

void TTrackManager::maskCurl

(

TTrack &

t

)

const

Definition at line 1216 of file TTrackManager.cxx.

                                              {
 
  //...Divide into two tracks...
  AList<TMLink> l[2];
  TMLink&       start = divideByIp( t, l );
  if ( l[0].length() == 0 ) return;
  if ( l[1].length() == 0 ) return;
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    curl : divided by IP" << std::endl;
  std::cout << "    0:";
  Dump( l[0], "flag sort" );
  std::cout << "    1:";
  Dump( l[1], "flag sort" );
  std::cout << std::endl;
#endif
 
  //...Which should be masked out ?...
  unsigned maskSide = 2;
 
  //...1. Check by # of super layers...
  if ( NSuperLayers( l[0] ) < NSuperLayers( l[1] ) ) maskSide = 0;
  else if ( NSuperLayers( l[0] ) > NSuperLayers( l[1] ) ) maskSide = 1;
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    NSuperLayers 0, 1 = " << NSuperLayers( l[0] ) << ", ";
  std::cout << NSuperLayers( l[1] ) << std::endl;
#endif
  if ( maskSide != 2 )
  {
    maskOut( t, l[maskSide] );
    return;
  }
 
  //...Make two tracks...
  TTrack* tt[2];
  tt[0] = new TTrack( t );
  tt[1] = new TTrack( t );
  tt[0]->remove( l[1] );
  tt[1]->remove( l[0] );
  _fitter.fit( *tt[0] );
  _fitter.fit( *tt[1] );
  Helix h0 = Helix( tt[0]->helix() );
  Helix h1 = Helix( tt[1]->helix() );
 
  //...Check by z...
  h0.pivot( ORIGIN );
  h1.pivot( ORIGIN );
  if ( fabs( h0.dz() ) < fabs( h1.dz() ) ) maskSide = 1;
  else maskSide = 0;
 
  delete tt[0];
  delete tt[1];
  maskOut( t, l[maskSide] );
  return;
}

Referenced by mask().

maskCurl() [2/3]

void TTrackManager::maskCurl

(

TTrack &

)

const

maskCurl() [3/3]

void TTrackManager::maskCurl

(

TTrack &

)

const

maskCurlHits() [1/3]

void TTrackManager::maskCurlHits	(	const AList< TMDCWireHit > &	axial,
		const AList< TMDCWireHit > &	stereo,
		const AList< TTrack > &	tracks ) const

masks hits on found curl tracks.

Definition at line 150 of file TTrackManager.cxx.

                                                                           {
  //...Coded by jtanaka...
 
  int i = 0;
  while ( const TTrack* t = tracks[i++] )
  {
    int j = 0;
    while ( TMDCWireHit* a = axial[j++] )
    {
      double x = t->helix().center().x() - a->xyPosition().x();
      double y = t->helix().center().y() - a->xyPosition().y();
      double r = sqrt( x * x + y * y );
      double R = fabs( t->helix().radius() );
      double q = t->helix().center().x() * a->xyPosition().y() -
                 t->helix().center().y() * a->xyPosition().x();
      double qq = q * t->charge();
      if ( R - 2. < r && r < R + 2. && qq > 0. ) { a->state( a->state() | WireHitUsed ); }
    }
    j = 0;
    while ( TMDCWireHit* s = stereo[j++] )
    {
      double x = t->helix().center().x() - s->xyPosition().x();
      double y = t->helix().center().y() - s->xyPosition().y();
      double r = sqrt( x * x + y * y );
      double R = fabs( t->helix().radius() );
      double q = t->helix().center().x() * s->xyPosition().y() -
                 t->helix().center().y() * s->xyPosition().x();
      double qq = q * t->charge();
      if ( R - 2.5 < r && r < R + 2.5 && qq > 0. ) { s->state( s->state() | WireHitUsed ); }
    }
  }
}

maskCurlHits() [2/3]

void TTrackManager::maskCurlHits	(	const AList< TMDCWireHit > &	axial,
		const AList< TMDCWireHit > &	stereo,
		const AList< TTrack > &	tracks ) const

masks hits on found curl tracks.

maskCurlHits() [3/3]

void TTrackManager::maskCurlHits	(	const AList< TMDCWireHit > &	axial,
		const AList< TMDCWireHit > &	stereo,
		const AList< TTrack > &	tracks ) const

masks hits on found curl tracks.

maskMultiHits() [1/3]

void TTrackManager::maskMultiHits

(

TTrack &

t

)

const

Definition at line 1103 of file TTrackManager.cxx.

                                                   {
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "... masking multi-hits" << std::endl;
#endif
 
  if ( !t.cores().length() ) return;
  AList<TMLink> cores        = t.cores();
  unsigned      n            = cores.length();
  bool          layerLimited = false;
  AList<TMLink> bads;
 
  cores.sort( SortByWireId );
  for ( unsigned i = 0; i < n; i++ )
  {
    if ( layerLimited )
    {
      bads.append( cores[i] );
      continue;
    }
    AList<TMLink> linksInLayer = SameLayer( cores, cores[i]->wire()->layerId() );
    if ( linksInLayer.length() > 3 )
    {
      bads.append( cores[i] );
      layerLimited = true;
    }
  }
  maskOut( t, bads );
}

Referenced by mask().

maskMultiHits() [2/3]

void TTrackManager::maskMultiHits

(

TTrack &

)

const

maskMultiHits() [3/3]

void TTrackManager::maskMultiHits

(

TTrack &

)

const

maskNormal() [1/3]

void TTrackManager::maskNormal

(

TTrack &

t

)

const

Definition at line 1132 of file TTrackManager.cxx.

                                                {
 
  //...Divide into two tracks...
  AList<TMLink> l[2];
  TMLink&       start = divideByIp( t, l );
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    normal : divided by IP" << std::endl;
  std::cout << "    0=";
  for ( unsigned j = 0; j < l[0].length(); j++ )
  { std::cout << "," << l[0][j]->wire()->name(); }
  std::cout << std::endl;
  std::cout << "    1=";
  for ( unsigned j = 0; j < l[1].length(); j++ )
  { std::cout << "," << l[1][j]->wire()->name(); }
  std::cout << std::endl;
#endif
 
  //...Which should be masked out ?...
  unsigned maskSide = 2;
 
  //...1. Check by # of super layers...
  if ( NSuperLayers( l[0] ) < NSuperLayers( l[1] ) ) maskSide = 0;
  else if ( NSuperLayers( l[0] ) > NSuperLayers( l[1] ) ) maskSide = 1;
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    NSuperLayers 0, 1 = " << NSuperLayers( l[0] ) << ", ";
  std::cout << NSuperLayers( l[1] ) << std::endl;
#endif
  if ( maskSide != 2 )
  {
    maskOut( t, l[maskSide] );
    return;
  }
 
  //...2. Check by the inner-most layer...
  unsigned i0 = InnerMost( l[0] )->wire()->layerId();
  unsigned i1 = InnerMost( l[1] )->wire()->layerId();
  if ( i0 < i1 ) maskSide = 1;
  else if ( i0 > i1 ) maskSide = 0;
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    i0, i1 = " << i0 << ", " << i1 << std::endl;
#endif
  if ( maskSide != 2 )
  {
    maskOut( t, l[maskSide] );
    return;
  }
 
  //...3. Check by # of layers...
  if ( NLayers( l[0] ) < NLayers( l[1] ) ) maskSide = 0;
  else if ( NLayers( l[0] ) > NLayers( l[1] ) ) maskSide = 1;
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    NLayers 0, 1 = " << NLayers( l[0] ) << ", ";
  std::cout << NLayers( l[1] ) << std::endl;
#endif
  if ( maskSide != 2 )
  {
    maskOut( t, l[maskSide] );
    return;
  }
 
  //...4. Check by pt...
  if ( maskSide == 2 )
  {
    TTrack* tt[2];
    for ( unsigned j = 0; j < 2; j++ )
    {
      tt[j] = new TTrack( t );
      tt[j]->remove( l[j] );
      _fitter.fit( *tt[j] );
    }
    if ( tt[1]->pt() > tt[0]->pt() ) maskSide = 1;
    else maskSide = 0;
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    pt 0 = " << tt[1]->pt() << std::endl;
    std::cout << "    pt 1 = " << tt[0]->pt() << std::endl;
#endif
    delete tt[0];
    delete tt[1];
  }
  maskOut( t, l[maskSide] );
  return;
}

Referenced by mask().

maskNormal() [2/3]

void TTrackManager::maskNormal

(

TTrack &

)

const

maskNormal() [3/3]

void TTrackManager::maskNormal

(

TTrack &

)

const

maskOut() [1/3]

void TTrackManager::maskOut	(	TTrack &	t,
		const AList< TMLink > &	links ) const

Definition at line 1088 of file TTrackManager.cxx.

                                                                         {
  unsigned n = links.length();
  if ( !n ) return;
  for ( unsigned i = 0; i < n; i++ )
  {
    const TMDCWireHit& hit = *links[i]->hit();
    hit.state( hit.state() | WireHitInvalidForFit );
  }
  t._fitted = false;
 
#ifdef TRKRECO_DEBUG_DETAIL
  Dump( links, "detail", "    TTrackManager::maskOut ... masking " );
#endif
}

Referenced by maskCurl(), maskMultiHits(), and maskNormal().

maskOut() [2/3]

void TTrackManager::maskOut	(	TTrack &	,
		const AList< TMLink > &	) const

maskOut() [3/3]

void TTrackManager::maskOut	(	TTrack &	,
		const AList< TMLink > &	) const

maxMomentum() [1/3]

double TTrackManager::maxMomentum ( double a )

inline

sets the max. momentum.

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

{ return _maxMomentum = a; }

maxMomentum() [2/3]

double TTrackManager::maxMomentum

(

double

)

sets the max. momentum.

maxMomentum() [3/3]

double TTrackManager::maxMomentum

(

double

)

sets the max. momentum.

merge() [1/3]

void TTrackManager::merge

(

void

)

Definition at line 1542 of file TTrackManager.cxx.

                                {
 
  //...Merging...
  unsigned n = _tracks.length();
  // cout<<"tracks: "<<n<<endl;
  AList<TTrack> bads;
  unsigned*     flagTrk = (unsigned*)malloc( n * sizeof( unsigned ) );
  for ( unsigned i = 0; i < n; i++ ) flagTrk[i] = 0;
 
  //...Search a track to be merged...
  for ( unsigned i0 = 0; i0 < n; i0++ )
  {
 
    if ( flagTrk[i0] != 0 ) continue;
    TTrack& t0 = *_tracks[i0];
    if ( !( t0.pt() < 0.13 ) ) continue;
 
    unsigned Noverlap( 0 ), Nall( 0 );
    float    OverlapRatioMax( -999. );
    unsigned MaxID( 0 );
 
    for ( unsigned i1 = 0; i1 < n; i1++ )
    {
 
      if ( i0 == i1 || flagTrk[i1] != 0 ) continue;
      TTrack& t1 = *_tracks[i1];
      if ( !( t1.pt() < 0.13 ) ) continue;
      Nall = t1.links().length();
      if ( !Nall ) continue;
 
      Noverlap = 0;
      for ( unsigned j = 0; j < Nall; j++ )
      {
        TMLink&            l    = *t1.links()[j];
        const TMDCWireHit& whit = *l.hit();
        double             load( 3. );                 // jialk original is 2
        if ( whit.state() & WireHitStereo ) load = 4.; // jialk original is 3
 
        double x = t0.helix().center().x() - l.positionOnTrack().x();
        double y = t0.helix().center().y() - l.positionOnTrack().y();
        double r = sqrt( x * x + y * y );
        double R = fabs( t0.helix().radius() );
 
        if ( ( R - load ) < r && r < ( R + load ) ) Noverlap++;
      }
 
      if ( !Noverlap ) continue;
      float tmpRatio = float( Noverlap ) / float( Nall );
 
      if ( tmpRatio > OverlapRatioMax )
      {
        OverlapRatioMax = tmpRatio;
        MaxID           = i1;
      }
    }
    // jialk caution original is 0.8
    if ( OverlapRatioMax < 0.7 ) continue;
 
    //...Mask should be done...
    unsigned      MaskID[2] = { MaxID, i0 };
    AList<TMLink> l[2];
 
    for ( unsigned j0 = 0; j0 < 2; j0++ )
    {
      for ( unsigned j1 = 0; j1 < _tracks[MaskID[j0]]->nLinks(); j1++ )
      {
        TMLink& k = *_tracks[MaskID[j0]]->links()[j1];
        l[j0].append( k );
      }
    }
    // _tracks[i0]->links().append( _tracks[MaxID]->links() );
    // _tracks[MaxID]->links().append( _tracks[i0]->links ());
    _tracks[i0]->append( _tracks[MaxID]->links() );
    _tracks[MaxID]->append( _tracks[i0]->links() );
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    mask & merge " << std::endl;
    std::cout << "    0:";
    Dump( l[0], "flag sort" );
    std::cout << "    1:";
    Dump( l[1], "flag sort" );
    std::cout << std::endl;
#endif
 
    //...Which should be masked out ?...
    unsigned maskSide = 2;
 
#if 0
    //...0. Check by # of super layers... ( not applied now )
    unsigned super0 = NSuperLayers(l[0]);
    unsigned super1 = NSuperLayers(l[1]);
 
    if( super0 < super1 ) maskSide = 0;
    else if ( super0 > super1 ) maskSide = 1;
 
#  ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    NSuperLayers 0, 1 = " << NSuperLayers(l[0]) << ", ";
    std::cout << NSuperLayers(l[1]) << std::endl;
#  endif
 
    if (maskSide == 2) {
#endif
 
    //...1. Check by the inner-most layer...
    unsigned inner0 = InnerMost( l[0] )->wire()->layerId();
    unsigned inner1 = InnerMost( l[1] )->wire()->layerId();
    if ( inner0 < inner1 ) maskSide = 1;
    else if ( inner0 > inner1 ) maskSide = 0;
 
    if ( maskSide == 2 )
    {
 
      //...2. Check by dz
 
      //...Make two tracks...
      TTrack* tt[2];
      tt[0] = new TTrack( *( _tracks[MaskID[0]] ) );
      tt[1] = new TTrack( *( _tracks[MaskID[1]] ) );
      _fitter.fit( *tt[0] );
      _fitter.fit( *tt[1] );
      Helix h0 = Helix( tt[0]->helix() );
      Helix h1 = Helix( tt[1]->helix() );
 
      //...Check dz...
      h0.pivot( ORIGIN );
      h1.pivot( ORIGIN );
      if ( fabs( h0.dz() ) < fabs( h1.dz() ) ) maskSide = 1;
      else maskSide = 0;
 
      delete tt[0];
      delete tt[1];
    }
#if 0
    }
#endif
    bads.append( _tracks[MaskID[maskSide]] );
    flagTrk[MaskID[maskSide]] = 1;
  }
 
  // cout<<"bads: "<<bads.length()<<endl;
  _tracks.remove( bads );
 
  //*****  Masking *****
  n = _tracks.length();
 
  for ( unsigned i = 0; i < n; i++ )
  {
    TTrack& t = *_tracks[i];
    for ( unsigned j = 0; j < t.links().length(); j++ )
    {
      TMLink&            l    = *t.links()[j];
      const TMDCWireHit& whit = *l.hit();
 
      if ( !( whit.state() & WireHitFittingValid ) ) continue;
 
      // within half circle or not?
      double q = t.helix().center().x() * l.positionOnTrack().y() -
                 t.helix().center().y() * l.positionOnTrack().x();
      double qq = q * t.charge();
 
      if ( qq > 0 ) whit.state( whit.state() & ~WireHitInvalidForFit );
      else whit.state( whit.state() | WireHitInvalidForFit );
    }
  }
 
  free( flagTrk );
}

merge() [2/3]

void TTrackManager::merge

(

void

)

merge() [3/3]

void TTrackManager::merge

(

void

)

movePivot() [1/3]

void TTrackManager::movePivot

(

void

)

moves pivot of tracks.

Definition at line 815 of file TTrackManager.cxx.

                                    {
  unsigned n = _tracksAll.length();
  //  unsigned n = _tracks.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    if ( _tracksAll[i]->nLinks() ) _tracksAll[i]->movePivot();
    //  if (_tracks[i]->nLinks()) _tracks[i]->movePivot();
  }
  nameTracks();
}

Referenced by finish().

movePivot() [2/3]

void TTrackManager::movePivot

(

void

)

moves pivot of tracks.

movePivot() [3/3]

void TTrackManager::movePivot

(

void

)

moves pivot of tracks.

name() [1/3]

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

inline

returns name.

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

{ return std::string( "Track Manager" ); }

Referenced by finish(), makeTds(), mask(), refit(), and salvage().

name() [2/3]

std::string TTrackManager::name

(

void

)

const

returns name.

name() [3/3]

std::string TTrackManager::name

(

void

)

const

returns name.

paraT0() [1/3]

float TTrackManager::paraT0 ( void ) const

inline

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

{ return _t0; }

paraT0() [2/3]

float TTrackManager::paraT0

(

void

)

const

paraT0() [3/3]

float TTrackManager::paraT0

(

void

)

const

refit() [1/3]

void TTrackManager::refit

(

void

)

refits tracks.

Definition at line 875 of file TTrackManager.cxx.

                                {
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << name() << " ... refitting" << std::endl;
#endif
  unsigned      n = _tracks.length();
  AList<TTrack> bads;
  for ( unsigned i = 0; i < n; i++ )
  {
    TTrack& t = *_tracks[i];
    int     err;
    err = _fitter.fit( t );
    if ( err < 0 )
    {
      bads.append( t );
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "    " << t.name();
      std::cout << " rejected because of fitting failure" << std::endl;
#endif
      continue;
    }
    t.refine( 30. * 10. );
    err = _fitter.fit( t );
    if ( err < 0 )
    {
      bads.append( t );
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "    " << t.name();
      std::cout << " rejected because of fitting failure" << std::endl;
#endif
      continue;
    }
    t.refine( 30. * 1. );
    err = _fitter.fit( t );
    if ( err < 0 )
    {
      bads.append( t );
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "    " << t.name();
      std::cout << " rejected because of fitting failure" << std::endl;
#endif
      continue;
    }
  }
  _tracks.remove( bads );
}

Referenced by finish().

refit() [2/3]

void TTrackManager::refit

(

void

)

refits tracks.

refit() [3/3]

void TTrackManager::refit

(

void

)

refits tracks.

removeHitsAcrossOverIp() [1/3]

void TTrackManager::removeHitsAcrossOverIp

(

AList< TMLink > &

l

)

const

Definition at line 1049 of file TTrackManager.cxx.

                                                                   {
 
  //...Calculate average phi...
  unsigned n      = l.length();
  float    phiSum = 0.;
  for ( unsigned i = 0; i < n; i++ )
  {
    const TMDCWire& w     = *l[i]->hit()->wire();
    unsigned        j     = w.localId();
    unsigned        nWire = w.layer()->nWires();
 
    float phi = (float)j / (float)nWire;
    phiSum += phi;
  }
  float average = phiSum / (float)n;
 
  AList<TMLink> cross;
  for ( unsigned i = 0; i < n; i++ )
  {
    const TMDCWire& w     = *l[i]->hit()->wire();
    unsigned        j     = w.localId();
    unsigned        nWire = w.layer()->nWires();
 
    float phi = (float)j / (float)nWire;
    float dif = fabs( phi - average );
    if ( dif > 0.5 ) dif = 1. - dif;
 
    if ( dif > 0.3 ) cross.append( l[i] );
  }
  l.remove( cross );
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    Cross over IP reduction : ";
  for ( unsigned i = 0; i < cross.length(); i++ )
  { std::cout << cross[i]->wire()->name() << ","; }
  std::cout << std::endl;
#endif
}

removeHitsAcrossOverIp() [2/3]

void TTrackManager::removeHitsAcrossOverIp

(

AList< TMLink > &

)

const

removeHitsAcrossOverIp() [3/3]

void TTrackManager::removeHitsAcrossOverIp

(

AList< TMLink > &

)

const

salvage() [1/3]

void TTrackManager::salvage

(

const AList< TMDCWireHit > &

hits

)

const

salvages remaining hits.

Definition at line 185 of file TTrackManager.cxx.

                                                                  {
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << name() << " ... salvaging" << std::endl;
  std::cout << "    # of given hits=" << hits.length() << std::endl;
#endif
 
  //...Check arguments...
  unsigned nTracks = _tracks.length();
  if ( nTracks == 0 ) return;
  unsigned nHits = hits.length();
  if ( nHits == 0 ) return;
 
  //...Hit loop...
  for ( unsigned i = 0; i < nHits; i++ )
  {
    TMDCWireHit& h = *hits[i];
 
    //...Already used?...
    if ( h.state() & WireHitUsed ) continue;
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    checking " << h.wire()->name() << std::endl;
#endif
 
    //...Select the closest track to a hit...
    TTrack* best = closest( _tracks, h );
#ifdef TRKRECO_DEBUG_DETAIL
    if ( !best )
    {
      std::cout << "        no track candidate returned";
      std::cout << "by TTrackManager::closest" << std::endl;
    }
#endif
    if ( !best ) continue;
 
    //...Try to append this hit...
    AList<TMLink> link;
    link.append( new TMLink( 0, &h ) );
    best->appendByApproach( link, 30. );
    // best->assign(WireHitConformalFinder);
    best->finder( TrackTrackManager );
  }
}

salvage() [2/3]

void TTrackManager::salvage

(

const AList< TMDCWireHit > &

)

const

salvages remaining hits.

salvage() [3/3]

void TTrackManager::salvage

(

const AList< TMDCWireHit > &

)

const

salvages remaining hits.

salvageAssociateHits() [1/3]

void TTrackManager::salvageAssociateHits	(	const AList< TMDCWireHit > &	hits,
		float	maxSigma2 )

salvages hits for dE/dx(not for track fitting).

Definition at line 2381 of file TTrackManager.cxx.

                                                                                          {
  // #define TRKRECO_DEBUG
  // #define TRKRECO_DEBUG_DETAIL
 
#ifdef TRKRECO_DEBUG
  std::cout << "... trkmgr::salvaging associate hits" << std::endl;
  std::cout << "    # of given hits=" << hits.length() << std::endl;
#endif
 
  //...Check arguments...
  unsigned nTracks = _tracks.length();
  if ( nTracks == 0 ) return;
  unsigned nHits = hits.length();
  if ( nHits == 0 ) return;
 
  static const TPoint2D o( 0., 0. );
 
  //...Hit loop...
  for ( unsigned i = 0; i < nHits; i++ )
  {
    TMDCWireHit& h = *hits[i];
 
    //...Already used?...
    if ( h.state() & WireHitUsed ) continue;
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    checking " << h.wire()->name();
#endif
 
    //...Track loop...
    AList<TMLink> toBeDeleted;
    TMLink*       best      = NULL;
    TTrack*       bestTrack = NULL;
    for ( unsigned j = 0; j < nTracks; j++ )
    {
      TTrack& t = *_tracks[j];
 
      //...Pre-selection...
      TPoint2D c  = t.center();
      TPoint2D co = -c;
      TPoint2D x  = h.wire()->xyPosition();
 
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << " : c= " << co.cross( x - c ) * t.charge();
      std::cout << ", d=" << fabs( ( x - c ).mag() - fabs( t.radius() ) );
#endif
 
      if ( co.cross( x - c ) * t.charge() > 0. ) continue;
      if ( fabs( ( x - c ).mag() - fabs( t.radius() ) ) > 5. ) continue;
 
      //...Try to append this hit...
      TMLink& link = *new TMLink( 0, &h );
      int     err  = t.approach( link );
      if ( err < 0 )
      {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << " : " << t.name() << " approach failure";
#endif
        toBeDeleted.append( link );
        continue;
      }
 
      //...Calculate sigma...
      float distance = link.distance();
      float diff     = fabs( distance - link.hit()->drift() );
      float sigma    = diff / link.hit()->dDrift();
      link.pull( sigma * sigma );
 
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << " : " << t.name() << " pull = " << link.pull();
#endif
      if ( link.pull() > maxSigma2 )
      {
        toBeDeleted.append( link );
        continue;
      }
 
      if ( best )
      {
        if ( best->pull() > link.pull() )
        {
          toBeDeleted.append( best );
          best      = &link;
          bestTrack = &t;
        }
        else { toBeDeleted.append( link ); }
      }
      else
      {
        best      = &link;
        bestTrack = &t;
      }
    }
 
    if ( best )
    {
      bestTrack->append( *best );
      best->hit()->state( best->hit()->state() | WireHitInvalidForFit );
      _associateHits.append( best );
#ifdef TRKRECO_DEBUG
      std::cout << "    " << best->hit()->wire()->name();
      std::cout << " -> " << bestTrack->name() << std::endl;
#endif
    }
    HepAListDeleteAll( toBeDeleted );
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << std::endl;
#endif
  }
}

salvageAssociateHits() [2/3]

void TTrackManager::salvageAssociateHits	(	const AList< TMDCWireHit > &	,
		float	maxSigma2 )

salvages hits for dE/dx(not for track fitting).

salvageAssociateHits() [3/3]

void TTrackManager::salvageAssociateHits	(	const AList< TMDCWireHit > &	,
		float	maxSigma2 )

salvages hits for dE/dx(not for track fitting).

saveMCTables() [1/3]

void TTrackManager::saveMCTables

(

void

)

const

stores MC track info. into Panther table.

Definition at line 768 of file TTrackManager.cxx.

                                             {
  unsigned n = _tracksFinal.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    const TTrack& t = *_tracksFinal[i];
 
    //  struct reccdc_trk * r;
    //  r = (struct reccdc_trk *) BsGetEnt(RECMDC_TRK, i + 1, BBS_No_Index);
    MdcRec_trk* r = &( *MdcRecTrkCol::getMdcRecTrkCol() )[i];
 
    //...Set type...
 
    //...Hit loop...
    const AList<TMLink>& hits  = t.finalHits();
    unsigned             nHits = hits.length();
    for ( unsigned j = 0; j < nHits; j++ )
    {
      TMLink*          l = hits[j];
      MdcRec_wirhit*   h = l->hit()->reccdc();
      MdcDat_mcwirhit* m = l->hit()->mc()->datcdc();
      m->trk             = r;
      //        struct reccdc_mctrk2hep * c;
      //        c = (struct reccdc_mctrk2hep *) BsNewEnt(RECMDC_MCTRK2HEP);
      MdcRec_mctrk2hep* c = new MdcRec_mctrk2hep;
      MdcRecMctrk2hepCol::getMdcRecMctrk2hepCol()->push_back( *c );
      c->wir = h;
      c->trk = r;
      c->hep = l->hit()->mc()->hep()->gen();
    }
 
    const TTrackMC* const mc = t.mc();
    //  struct reccdc_mctrk * m;
    //  m = (struct reccdc_mctrk *) BsNewEnt(RECMDC_MCTRK);
    // //   MdcRec_mctrk* m = &(*MdcRecMctrkCol::getMdcRecMctrkCol())[0];
    MdcRec_mctrk* m = new MdcRec_mctrk;
    MdcRecMctrkCol::getMdcRecMctrkCol()->push_back( *m );
    m->wirFrac    = mc->wireFraction();
    m->wirFracHep = mc->wireFractionHEP();
    m->charge     = mc->charge();
    m->ptFrac     = mc->ptFraction();
    m->pzFrac     = mc->pzFraction();
    m->quality    = mc->quality();
    if ( mc->hep() ) m->hep = mc->hep()->gen();
    else m->hep = 0;
  }
}

saveMCTables() [2/3]

void TTrackManager::saveMCTables

(

void

)

const

stores MC track info. into Panther table.

saveMCTables() [3/3]

void TTrackManager::saveMCTables

(

void

)

const

stores MC track info. into Panther table.

saveTables() [1/3]

void TTrackManager::saveTables

(

void

)

stores track info. into Panther table.

Definition at line 617 of file TTrackManager.cxx.

                                     {
#ifdef TRKRECO_DEBUG
  std::cout << "TTrackManager::saveTables ... # 3D tracks=" << _tracks.length()
            << ", # 2D tracks=" << _tracks2D.length() << ", all tracks=" << _tracksAll.length()
            << std::endl;
#endif
 
  //...For 3D tracks...
  AList<TTrack> badTracks;
  unsigned      n  = _tracks.length();
  unsigned*     id = (unsigned*)malloc( n * sizeof( unsigned ) );
  //    bzero((char *) id, n * sizeof(unsigned));
  memset( (char*)id, 0, n * sizeof( unsigned ) );
  for ( unsigned i = 0; i < n; i++ )
  {
    TTrack& t = *_tracks[i];
    if ( !t.nLinks() )
    {
      badTracks.append( (TTrack&)t );
      continue;
    }
 
    //...Copy track parameters...
    MdcRec_trk*     r   = 0;
    MdcRec_trk_add* a   = 0;
    int             err = copyTrack( t, &r, &a );
    if ( err )
    {
      badTracks.append( t );
      continue;
    }
    _tracksFinal.append( t );
 
    //...Type and quality...
    id[i]       = r->id;
    r->stat     = t.state();
    a->kind     = t.type();
    a->decision = t.finder();
    a->stat     = t.fitting();
    //      if (a->m_kind == TrackTypeCosmic) {
    //          a->m_quality = TrackQualityCosmic;
    //      }
    //      if (t.daughter() && (_tracks.index(t.daughter()) >= 0))
    //          a->m_daughter = _tracks.index(t.daughter()) + 1;
  }
 
  //...Daughter treatment...
  for ( unsigned i = 0; i < n; i++ )
  {
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "id[" << i << "]=" << id[i] << std::endl;
#endif
    if ( !( id[i] ) ) continue;
    if ( !( _tracks[i]->daughter() ) ) continue;
 
    int dId = _tracks.index( _tracks[i]->daughter() );
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    dId=" << dId;
    if ( dId >= 0 ) std::cout << ", id[dId]=" << id[dId];
    std::cout << std::endl;
#endif
 
    if ( dId >= 0 )
    {
      if ( id[dId] )
      {
        //      reccdc_trk_add * a = (reccdc_trk_add *)
        //          BsGetEnt(RECMDC_TRK_ADD, id[i], BBS_No_Index);
        //      a->m_daughter = id[dId];
        MdcRec_trk_add* a = &( *MdcRecTrkAddCol::getMdcRecTrkAddCol() )[i];
      }
    }
  }
  free( id );
 
  //...Remove bad tracks...
  _tracks.remove( badTracks );
  badTracks.removeAll();
 
  //...For 2D tracks...
  n = _tracks2D.length();
  for ( unsigned i = 0; i < n; i++ )
  {
    TTrack& t = *_tracks2D[i];
 
    //...Copy track parameters...
    MdcRec_trk*     r   = 0;
    MdcRec_trk_add* a   = 0;
    int             err = copyTrack( t, &r, &a );
    if ( err )
    {
#ifdef TRKRECO_DEBUG
      std::cout << "TTrackManager::saveTables !!! bad 2D tracks found"
                << " : err=" << err << std::endl
                << TrackDump( t ) << std::endl;
#endif
      badTracks.append( t );
      continue;
    }
    _tracksFinal.append( t );
 
    //...Reset helix parameter...
    //      r->m_helix[3] = 0.;
    //      r->m_helix[4] = 0.;
    //      r->m_nhits -= r->m_nster;
    //      r->m_nster = 0;
 
    //...Table filling...
    r->stat     = t.state();
    a->kind     = t.type();
    a->decision = t.finder();
    // a->m_quality = t.quality();
    a->quality = TrackQuality2D;
    a->stat    = t.fitting();
 
#ifdef TRKRECO_DEBUG
    if ( ( r->ndf == 0 ) && ( r->chiSq > 0. ) )
    {
      std::cout << "TTrackManager::saveTables !!! chisq>0 with ndf=0." << std::endl
                << "    Here is a track dump"
                << "    " << TrackDump( t ) << std::endl;
      t.dump( "detail" );
    }
    if ( ( r->ndf > 0 ) && ( r->chiSq == 0. ) )
    {
      std::cout << "TTrackManager::saveTables !!! chisq=0 with ndf>0." << std::endl
                << "    Here is a track dump"
                << "    " << TrackDump( t ) << std::endl;
      t.dump( "detail" );
    }
 
    if ( r->ndf == 0 )
      std::cout << "TTrackManager::saveTables ... ndf = 0" << std::endl
                << "    " << TrackDump( t ) << std::endl;
    if ( r->chiSq == 0. )
      std::cout << "TTrackManager::saveTables ... chisq = 0" << std::endl
                << "    " << TrackDump( t ) << std::endl;
#endif
  }
  _tracks2D.remove( badTracks );
 
  //...Statistics...
  ++_s->_nEvents;
  _s->_nTracks += _tracks.length();
  _s->_nTracksAll += _tracksAll.length();
  _s->_nTracks2D += _tracks2D.length();
  _s->_nTracksFinal += _tracksFinal.length();
}

saveTables() [2/3]

void TTrackManager::saveTables

(

void

)

stores track info. into Panther table.

saveTables() [3/3]

void TTrackManager::saveTables

(

void

)

stores track info. into Panther table.

setCurlerFlags() [1/3]

void TTrackManager::setCurlerFlags

(

void

)

tests for curlers.

Definition at line 2344 of file TTrackManager.cxx.

                                         {
  unsigned n = _tracks.length();
  if ( n < 2 ) return;
 
  for ( unsigned i = 0; i < n - 1; i++ )
  {
    TTrack& t0 = *_tracks[i];
    if ( t0.type() != TrackTypeCurl ) continue;
    float c0 = t0.charge();
 
    for ( unsigned j = i + 1; j < n; j++ )
    {
      TTrack& t1 = *_tracks[j];
      float   c1 = t1.charge();
      if ( c0 * c1 > 0. ) continue;
      if ( t1.type() != TrackTypeCurl ) continue;
 
      bool     toBeMerged = false;
      unsigned n0         = t0.testByApproach( t1.cores(), _sigmaCurlerMergeTest );
      if ( n0 > _nCurlerMergeTest ) toBeMerged = true;
      if ( !toBeMerged )
      {
        unsigned n1 = t1.testByApproach( t0.cores(), _sigmaCurlerMergeTest );
        if ( n1 > _nCurlerMergeTest ) toBeMerged = true;
      }
 
      if ( toBeMerged )
      {
        //      ++_s->_nToBeMerged;
        if ( ( t0.daughter() ) || ( t1.daughter() ) ) ++_s->_nToBeMergedMoreThanTwo;
        t0.daughter( &t1 );
        t1.daughter( &t0 );
      }
    }
  }
}

setCurlerFlags() [2/3]

void TTrackManager::setCurlerFlags

(

void

)

tests for curlers.

setCurlerFlags() [3/3]

void TTrackManager::setCurlerFlags

(

void

)

tests for curlers.

sett0bes() [1/3]

void TTrackManager::sett0bes ( double t0 )

inline

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

{ t0bes = t0; }

sett0bes() [2/3]

void TTrackManager::sett0bes ( double t0 )

inline

Definition at line 201 of file InstallArea/x86_64-el9-gcc13-opt/include/TrkReco/TTrackManager.h.

{ t0bes = t0; }

sett0bes() [3/3]

void TTrackManager::sett0bes ( double t0 )

inline

Definition at line 201 of file Reconstruction/TrkReco/include/TrkReco/TTrackManager.h.

{ t0bes = t0; }

sortBanksByPt() [1/3]

void TTrackManager::sortBanksByPt

(

void

)

const

sorts RECMDC_TRK tables.

Definition at line 2019 of file TTrackManager.cxx.

                                              {
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "trkmgr::sortBanksByPt : # of tracks="
            //   << BsCouTab(RECMDC_TRK_ADD) << std::endl;
            << MdcRecTrkAddCol::getMdcRecTrkAddCol()->size() << std::endl;
#endif
 
  //    unsigned n = BsCouTab(RECMDC_TRK_ADD);
  unsigned n = MdcRecTrkAddCol::getMdcRecTrkAddCol()->size();
  if ( n < 2 ) return;
 
  //...Sort RECMDC...
  unsigned* id = (unsigned*)malloc( n * sizeof( unsigned ) );
  for ( unsigned i = 0; i < n; i++ ) id[i] = i;
  for ( unsigned i = 0; i < n - 1; i++ )
  {
    //  reccdc_trk & cdc0 =
    //      * (reccdc_trk *) BsGetEnt(RECMDC_TRK,
    //                    i + 1,
    //                    BBS_No_Index);
    //  reccdc_trk_add & add0 =
    //      * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD,
    //                    i + 1,
    //                    BBS_No_Index);
    //  reccdc_mctrk & mc0 =
    //      * (reccdc_mctrk *) BsGetEnt(RECMDC_MCTRK,
    //                  i + 1,
    //                  BBS_No_Index);
    MdcRec_trk&     cdc0 = ( *MdcRecTrkCol::getMdcRecTrkCol() )[i];
    MdcRec_trk_add& add0 = ( *MdcRecTrkAddCol::getMdcRecTrkAddCol() )[i];
    MdcRec_mctrk&   mc0  = ( *MdcRecMctrkCol::getMdcRecMctrkCol() )[i];
 
    float    bestPt      = 1. / fabs( cdc0.helix[2] );
    unsigned bestQuality = add0.quality;
    for ( unsigned j = i + 1; j < n; j++ )
    {
      //        reccdc_trk & cdc1 =
      //        * (reccdc_trk *) BsGetEnt(RECMDC_TRK,
      //                      j + 1,
      //                      BBS_No_Index);
      //        reccdc_trk_add & add1 =
      //        * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD,
      //                          j + 1,
      //                          BBS_No_Index);
      //        reccdc_mctrk & mc1 =
      //        * (reccdc_mctrk *) BsGetEnt(RECMDC_MCTRK,
      //                        j + 1,
      //                        BBS_No_Index);
      MdcRec_trk&     cdc1 = ( *MdcRecTrkCol::getMdcRecTrkCol() )[j];
      MdcRec_trk_add& add1 = ( *MdcRecTrkAddCol::getMdcRecTrkAddCol() )[j];
      MdcRec_mctrk&   mc1  = ( *MdcRecMctrkCol::getMdcRecMctrkCol() )[j];
 
      float pt = 1. / fabs( cdc1.helix[2] );
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "i,j=" << i << "," << j << " : quality i,j=" << bestQuality << ","
                << add1.quality << std::endl;
#endif
      unsigned quality = add1.quality;
      if ( quality > bestQuality ) continue;
      else if ( quality < bestQuality )
      {
        bestQuality = quality;
        bestPt      = pt;
        swapReccdc( cdc0, add0, mc0, cdc1, add1, mc1 );
        unsigned tmp = id[i];
        id[i]        = id[j];
        id[j]        = tmp;
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "swapped" << std::endl;
#endif
        continue;
      }
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "i,j=" << i << "," << j << " : pt i,j=" << bestPt << "," << pt << std::endl;
#endif
      if ( pt > bestPt )
      {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "swapping ... " << &cdc0 << "," << &add0 << "," << &mc0 << " <-> "
                  << &cdc1 << "," << &add1 << "," << &mc1 << std::endl;
#endif
        bestQuality = quality;
        bestPt      = pt;
        swapReccdc( cdc0, add0, mc0, cdc1, add1, mc1 );
        unsigned tmp = id[i];
        id[i]        = id[j];
        id[j]        = tmp;
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "swapped" << std::endl;
#endif
      }
    }
  }
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "trkmgr::sortBanksByPt : first phase finished" << std::endl;
#endif
 
  tagReccdc( id, n );
  free( id );
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "trkmgr::sortBanksByPt : second phase finished" << std::endl;
#endif
 
#if 0
  //...Sort RECTRK...
  n = BsCouTab(RECTRK);
  id = (unsigned *) malloc(n * sizeof(unsigned));
  for (unsigned i = 0; i < n; i++) id[i] = i;
  if (n > 1) {
    unsigned i = 0;
    while (i < n - 1) {
      rectrk & t = * (rectrk *) BsGetEnt(RECTRK, i + 1, BBS_No_Index);
      if (t.m_prekal == (i + 1)) {
    ++i;
    continue;
      }
 
      rectrk & s = * (rectrk *) BsGetEnt(RECTRK,
      t.m_prekal,
      BBS_No_Index);
 
      swapRectrk(t, s);
      unsigned tmp = id[i];
      id[i] = id[s.m_ID - 1];
      id[s.m_ID - 1] = tmp;
 
      //std::cout << "swap " << i + 1 << " and " << s.m_ID << std::endl;
 
      reccdc_trk_add & a =
    * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD,
        t.m_prekal,
        BBS_No_Index);
      reccdc_trk_add & b =
    * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD,
        s.m_prekal,
        BBS_No_Index);
      a.m_rectrk = t.m_ID;
      b.m_rectrk = s.m_ID;
    }
  }
#else
  /*
  // jtanaka 000925 -->
  n = BsCouTab(RECTRK);
  id = (unsigned *) malloc(n * sizeof(unsigned));
  for (unsigned i = 0; i < n; i++) id[i] = i;
  int foundId = 0;
  while(foundId != n){
  rectrk & t = * (rectrk *) BsGetEnt(RECTRK, foundId + 1, BBS_No_Index);
  int minPrekal = t.m_prekal;
  int exchangeId = foundId;
  for(int i=foundId+1;i<n;++i){
  rectrk & s = * (rectrk *) BsGetEnt(RECTRK, i + 1, BBS_No_Index);
  if(s.m_prekal < minPrekal){
  minPrekal = s.m_prekal;
  exchangeId = i;
  }
  }
  if(exchangeId != foundId){
  rectrk & s = * (rectrk *) BsGetEnt(RECTRK,
  exchangeId + 1,
  BBS_No_Index);
 
  swapRectrk(t, s);
  unsigned tmp = id[t.m_ID - 1];
  id[t.m_ID - 1] = id[s.m_ID - 1];
  id[s.m_ID - 1] = tmp;
 
  reccdc_trk_add & a =
   * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD,
   t.m_prekal,
   BBS_No_Index);
   reccdc_trk_add & b =
   * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD,
   s.m_prekal,
   BBS_No_Index);
   a.m_rectrk = t.m_ID;
   b.m_rectrk = s.m_ID;
   }
   ++foundId;
   }
  // <-- jtanaka 000925
  */
#endif
 
  tagRectrk( id, n );
  free( id );
}

sortBanksByPt() [2/3]

void TTrackManager::sortBanksByPt

(

void

)

const

sorts RECMDC_TRK tables.

sortBanksByPt() [3/3]

void TTrackManager::sortBanksByPt

(

void

)

const

sorts RECMDC_TRK tables.

sortTracksByPt() [1/3]

void TTrackManager::sortTracksByPt

(

void

)

Definition at line 1888 of file TTrackManager.cxx.

                                         {
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "trkmgr::sortTracksByPt : # of tracks=" << _tracks.length() << std::endl;
#endif
 
  unsigned n = _tracks.length();
  if ( n < 2 ) return;
 
  for ( unsigned i = 0; i < n - 1; i++ )
  {
    TTrack& t0     = *_tracks[i];
    float   bestPt = t0.pt();
    for ( unsigned j = i + 1; j < n; j++ )
    {
      TTrack& t1 = *_tracks[j];
      float   pt = t1.pt();
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "i,j=" << i << "," << j << " : pt i,j=" << bestPt << "," << pt << std::endl;
#endif
      if ( pt > bestPt )
      {
        bestPt = pt;
        _tracks.swap( i, j );
      }
    }
  }
}

sortTracksByPt() [2/3]

void TTrackManager::sortTracksByPt

(

void

)

sortTracksByPt() [3/3]

void TTrackManager::sortTracksByPt

(

void

)

sortTracksByQuality() [1/3]

void TTrackManager::sortTracksByQuality

(

void

)

sorts tracks.

Definition at line 1865 of file TTrackManager.cxx.

                                              {
  unsigned n = _tracks.length();
  if ( n < 2 ) return;
 
  for ( unsigned i = 0; i < n - 1; i++ )
  {
    TTrack& t0         = *_tracks[i];
    float   bestRChisq = HUGE_VAL;
    if ( t0.ndf() > 0 ) bestRChisq = t0.chi2() / t0.ndf();
    for ( unsigned j = i + 1; j < n; j++ )
    {
      TTrack& t1     = *_tracks[j];
      float   rChisq = HUGE_VAL;
      if ( t1.ndf() > 0 ) rChisq = t1.chi2() / t1.ndf();
      if ( rChisq < bestRChisq )
      {
        bestRChisq = rChisq;
        _tracks.swap( i, j );
      }
    }
  }
}

sortTracksByQuality() [2/3]

void TTrackManager::sortTracksByQuality

(

void

)

sorts tracks.

sortTracksByQuality() [3/3]

void TTrackManager::sortTracksByQuality

(

void

)

sorts tracks.

tracks() [1/3]

const AList< TTrack > & TTrackManager::tracks ( void ) const

inline

returns a list of reconstructed tracks.

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

{ return _tracks; }

Referenced by closest(), maskBadHits(), and maskCurlHits().

tracks() [2/3]

const AList< TTrack > & TTrackManager::tracks

(

void

)

const

returns a list of reconstructed tracks.

tracks() [3/3]

const AList< TTrack > & TTrackManager::tracks

(

void

)

const

returns a list of reconstructed tracks.

tracks2D() [1/3]

const AList< TTrack > & TTrackManager::tracks2D ( void ) const

inline

returns a list of 2D tracks.

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

{ return _tracks2D; }

tracks2D() [2/3]

const AList< TTrack > & TTrackManager::tracks2D

(

void

)

const

returns a list of 2D tracks.

tracks2D() [3/3]

const AList< TTrack > & TTrackManager::tracks2D

(

void

)

const

returns a list of 2D tracks.

tracksFinal() [1/3]

const AList< TTrack > & TTrackManager::tracksFinal ( void ) const

inline

returns a list of tracks writen to MdcRec_trk.

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

{ return _tracksFinal; }

tracksFinal() [2/3]

const AList< TTrack > & TTrackManager::tracksFinal

(

void

)

const

returns a list of tracks writen to MdcRec_trk.

tracksFinal() [3/3]

const AList< TTrack > & TTrackManager::tracksFinal

(

void

)

const

returns a list of tracks writen to MdcRec_trk.

treatCurler() [1/3]

void TTrackManager::treatCurler	(	MdcTrk &	curl,
		MdcRec_trk_add &	cdc,
		unsigned	flag ) const

final decision for a curler.

Definition at line 1916 of file TTrackManager.cxx.

                                                                                         {
  //...Originally coded by j.tanaka...
 
  //...Check inputs...
  if ( trk1.zero[2] == 0 ) return;
  if ( cdc1.daughter == 0 ) return;
 
  //...The other side...
  //    reccdc_trk_add & cdc2 = * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD,
  //                              cdc1.m_daughter,
  //                              BBS_No_Index);
  //    MdcRec_trk_add & cdc2 = (*MdcRecTrkAddCol::getMdcRecTrkAddCol())[cdc1.daughter.id];
  MdcRec_trk_add& cdc2 = *cdc1.daughter->add;
 
  if ( cdc2.daughter == 0 ) return;
  if ( cdc2.rectrk == 0 ) return;
  //    rectrk & trk2 = * (rectrk *) BsGetEnt(RECTRK, cdc2.m_rectrk, BBS_No_Index);
  MdcTrk& trk2 = *cdc2.rectrk;
 
  if ( trk2.zero[2] == 0 ) return;
 
  //...Obtain RECTRK_LOCALZ...
  //    rectrk_localz & z1 = * (rectrk_localz *) BsGetEnt(RECTRK_LOCALZ,
  //                              trk1.m_zero[2],
  //                              BBS_No_Index);
  //    rectrk_localz & z2 = * (rectrk_localz *) BsGetEnt(RECTRK_LOCALZ,
  //                              trk2.m_zero[2],
  //                              BBS_No_Index);
  MdcTrk_localz& z1 = *trk1.zero[2];
  MdcTrk_localz& z2 = *trk2.zero[2];
 
  //...Pointer to mother and daughter...
  MdcRec_trk_add* mother   = &cdc1;
  MdcRec_trk_add* daughter = &cdc2;
 
  //      //...By dr...
  //      if (flag == 1) {
  //    float dr1 = fabs(z1.m_helix[0]);
  //    float dr2 = fabs(z2.m_helix[0]);
  //    if (dr1 > dr2) {
  //        mother = & cdc2;
  //        daughter = & cdc1;
  //    }
  //      }
 
  //      //...By dz...
  //      else {
  //    float dz1 = fabs(z1.m_helix[3]);
  //    float dz2 = fabs(z2.m_helix[3]);
  //    if (dz1 > dz2) {
  //        mother = & cdc2;
  //        daughter = & cdc1;
  //    }
  //      }
 
  //...By dz + dr...
  if ( flag == 3 )
  {
    float dz1 = fabs( z1.helix[3] );
    float dz2 = fabs( z2.helix[3] );
    if ( fabs( dz1 - dz2 ) < 2. ) flag = 1;
    else flag = 2;
  }
 
  //...By dr...
  if ( flag == 1 )
  {
    float dr1 = fabs( z1.helix[0] );
    float dr2 = fabs( z2.helix[0] );
    if ( dr1 > dr2 )
    {
      mother   = &cdc2;
      daughter = &cdc1;
    }
  }
 
  //...By dz...
  else if ( flag == 2 )
  {
    float dz1 = fabs( z1.helix[3] );
    float dz2 = fabs( z2.helix[3] );
    if ( dz1 > dz2 )
    {
      mother   = &cdc2;
      daughter = &cdc1;
    }
  }
 
  //...Update information...
  mother->quality &= ( ~TrackQualityOutsideCurler );
  mother->likelihood[0] = 1.;
  mother->decision |= TrackTrackManager;
  // zsl
  mother->daughter = daughter->body;
  mother->mother   = 0;
  // zsl end;
  daughter->quality |= TrackQualityOutsideCurler;
  daughter->likelihood[0] = 0.;
  daughter->mother        = mother->body;
  daughter->daughter      = 0;
  daughter->decision |= TrackTrackManager;
}

treatCurler() [2/3]

void TTrackManager::treatCurler	(	MdcTrk &	curl,
		MdcRec_trk_add &	cdc,
		unsigned	flag ) const

final decision for a curler.

treatCurler() [3/3]

void TTrackManager::treatCurler	(	MdcTrk &	curl,
		MdcRec_trk_add &	cdc,
		unsigned	flag ) const

final decision for a curler.

version() [1/3]

std::string TTrackManager::version

(

void

)

const

returns version.

Definition at line 88 of file TTrackManager.cxx.

{ return std::string( "2.27" ); }

version() [2/3]

std::string TTrackManager::version

(

void

)

const

returns version.

version() [3/3]

std::string TTrackManager::version

(

void

)

const

returns version.

---

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

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