#include <FTFinder.h>

Public Member Functions
void	setAlgorithmPointer (Algorithm *)
	returns FTFinder pointer
	FTFinder ()
	Constructors and destructor.
void	init ()
	initializer(creates geometry)
void	term ()
	terminator
void	begin_run ()
	begin run function(reads constants)
void	event ()
	track finder core
FTSuperLayer *	superLayer (int id) const
	returns superlayer
FTList< FTTrack * > &	tracks ()
	returns track list
CLHEP::Hep3Vector	vertex () const
	returns event primary vertex
int	getWireId (FTWire *) const
	returns wire ID for given FTWire object
float	t2x (const FTLayer *l, const float t) const
	convert t to x
float	x2t (const FTLayer *l, const float x) const
	convert x to t
void	setBunchtime (double bunchtime)
	set bunchtime for MC events
void	setT0cal (double t0cal)
	set crude mdc time calibconst (~200 ns) for t0 calculation

Public Attributes
int	eventNo
int	runNo
int	expflag
float	t0Estime
float	tOffSet
float	evtTiming
const HepPoint3D	pivot
float	Testime
int	i_rPhiFit
int	tEstFlag
FTList< FTList< float > >	tEstime

Detailed Description

Definition at line 25 of file FTFinder.h.

Constructor & Destructor Documentation

◆ FTFinder()

FTFinder::FTFinder ( )

Constructors and destructor.

Definition at line 111 of file FTFinder.cxx.

    : tOffSet( 0 )
    , t0Estime( -999 )
    , tEstFlag( 0 )
    , tEstime( 10 )
    , _tracks( 0 )
    , _linkedSegments( 0 )
    , _axialSegCollect( 0 )
    , _vx( 0 )
    , _vy( 0 )
    , _vz( 0 )
    , _ExpNo( 0 )
    , _RunNo( 0 )
    , m_total_trk( 0 )
    , pivot( 0, 0, 0 ) {
  StatusCode scmgn = Gaudi::svcLocator()->service( "MagneticFieldSvc", m_pmgnIMF );
  if ( scmgn != StatusCode::SUCCESS )
  { std::cout << "Unable to open Magnetic field service" << std::endl; }
}

Member Function Documentation

◆ begin_run()

void FTFinder::begin_run ( )

begin run function(reads constants)

Definition at line 154 of file FTFinder.cxx.

                         {
  eventNo = 0;
  runNo   = 0;
  expflag = 0;
  if ( !param->_doIt ) return;
  _ExpNo = 0;
  _RunNo = 0;
 
  IMdcGeomSvc* mdcGeomSvc;
  StatusCode   sc = Gaudi::svcLocator()->service( "MdcGeomSvc", mdcGeomSvc );
 
  if ( !sc.isSuccess() ) return;
 
  m_geom = FTGeom::instance();
}

◆ event()

void FTFinder::event ( )

track finder core

– vertex(r-phi) fit and event timing correction

– 2D track re-fitting

Definition at line 170 of file FTFinder.cxx.

                     {
  if ( !param->_doIt ) return;
 
  IMessageSvc* msgSvc;
  Gaudi::svcLocator()->service( "MessageSvc", msgSvc ).ignore();
 
  MsgStream log( msgSvc, "FTFinder" );
 
  IDataProviderSvc* eventSvc;
  Gaudi::svcLocator()->service( "EventDataSvc", eventSvc ).ignore();
 
  //--
  // clear old information
  //--
  clear();
 
  // check whether  Recon  already registered
  DataObject* aReconEvent;
  eventSvc->findObject( "/Event/Recon", aReconEvent ).ignore();
  if ( !aReconEvent )
  {
    // register ReconEvent Data Object to TDS;
    ReconEvent* recevt = new ReconEvent;
    StatusCode  sc     = eventSvc->registerObject( "/Event/Recon", recevt );
    if ( sc != StatusCode::SUCCESS )
    {
      log << MSG::FATAL << "Could not register ReconEvent" << endmsg;
      return;
    }
  }
 
  // register Event start time
  IDataManagerSvc* dataManSvc;
  Gaudi::svcLocator()->service( "EventDataSvc", dataManSvc ).ignore();
 
  DataObject* aEsTimeEvent;
  eventSvc->findObject( "/Event/Recon/RecEsTimeCol", aEsTimeEvent ).ignore();
  if ( aEsTimeEvent )
  {
    dataManSvc->clearSubTree( "/Event/Recon/RecEsTimeCol" ).ignore();
    eventSvc->unregisterObject( "/Event/Recon/RecEsTimeCol" ).ignore();
  }
 
  RecEsTimeCol* aRecEsTimeCol = new RecEsTimeCol;
  StatusCode    est;
  est = eventSvc->registerObject( "/Event/Recon/RecEsTimeCol", aRecEsTimeCol );
  if ( est.isFailure() )
  {
    log << MSG::FATAL << "Could not register RecEsTimeCol" << endmsg;
    return;
  }
  log << MSG::DEBUG << "RecEsTimeCol registered successfully!" << endmsg;
 
  //--
  // update wirehit information
  //--
  updateMdc();
 
  //--
  // segment finding
  //--
  for ( int i = 0; i < 11; i++ ) { m_geom->getSuperLayer( i )->mkSegmentList(); }
 
  //--
  // reduce noise and get start time from segment fit
  //--
  for ( int i = 0; i < 10; i++ ) { m_geom->getSuperLayer( i )->reduce_noise( tEstime ); }
 
  getTestime();
 
  //--
  // axial segment linking
  //--
  mkTrackList();
 
  //--
  // 2D track fitting
  //--
  //(*_superLayer).reAppendSalvage();
  log << MSG::DEBUG << "number of 2D tracks: " << _tracks.size() << endmsg;
 
#ifndef OnlineMode
  num_2Dtrk += _tracks.size();
#endif
 
  for ( auto it = _tracks.begin(); it != _tracks.end(); )
  {
    if ( !( *it )->r_phiFit() )
    {
      delete *it;
      it            = _tracks.erase( it );
      int tmp_index = std::distance( _tracks.begin(), it );
      log << MSG::DEBUG << "===========> deleted 2D track :  " << tmp_index << endmsg;
    }
    else it++;
  }
 
  if ( t0Estime != -999 )
  {
    // begin to make Event start time
    RecEsTime* arecestime = new RecEsTime;
    arecestime->setTest( t0Estime );
    arecestime->setStat( tEstFlag );
    aRecEsTimeCol->push_back( arecestime );
  }
 
  if ( !_tracks.size() )
  {
    makeTds().ignore();
    return;
  }
 
  ///--
  /// vertex(r-phi) fit and event timing correction
  //--
  int vtx_flag = VertexFit( 0 );
  evtTiming    = ( param->_evtTimeCorr ) ? CorrectEvtTiming() : 0;
 
  ///--
  /// 2D track re-fitting
  //--
  if ( param->_hitscut == 1 )
  {
    for ( auto trk_i : _tracks ) { trk_i->r_phiReFit( _vx, _vy, vtx_flag ); }
  }
 
  //--
  // cut bad hits from 2D track re-fitting
  if ( param->_hitscut == 2 )
  {
    for ( auto trk_i : _tracks )
    {
      for ( int j = 0; j < 2; j++ )
      {
        i_rPhiFit = trk_i->r_phi2Fit( _vx, _vy, vtx_flag );
        if ( i_rPhiFit != -99 ) trk_i->r_phi3Fit( i_rPhiFit, _vx, _vy, vtx_flag );
      }
      trk_i->r_phi4Fit( _vx, _vy, vtx_flag );
    }
  }
 
  //--
  // stereo segment linking
  //--
  mkTrack3D();
 
  //--
  // final track fittng
  //--
  log << MSG::DEBUG << "number of 3D tracks: " << _tracks.size() << endmsg;
 
#ifndef OnlineMode
  num_3Dtrk += _tracks.size();
#endif
 
  for ( auto it = _tracks.begin(); it != _tracks.end(); )
  {
 
    int trk_index = std::distance( _tracks.begin(), it ) + 1;
 
#ifndef OnlineMode
    log << MSG::DEBUG << "=======> 3D track: " << trk_index << endmsg;
    ( *it )->printout();
#endif
 
    if ( ( *it )->get_nhits() < 18 )
    {
      delete ( *it );
      log << MSG::DEBUG << "================> deleted 3D track :  " << trk_index << endmsg;
      it = _tracks.erase( it );
    }
    else it++;
  }
 
  if ( !_tracks.size() )
  {
    makeTds().ignore();
    return;
  }
 
  for ( auto trk_i : _tracks ) trk_i->s_zFit();
 
  //--
  // find primary event vertex
  //--
  if ( param->_findEventVertex ) { VertexFit( 1 ); }
 
  log << MSG::DEBUG << "final number of tracks: " << _tracks.size() << endmsg;
 
#ifndef OnlineMode
  num_finaltrk += _tracks.size();
  if ( param->_mkMdst ) makeMdst();
#endif
 
  //--
  // output tracking result into TDS
  // by wangdy
  //--
  if ( param->_mkTds ) makeTds().ignore();
}

◆ getWireId()

int FTFinder::getWireId ( FTWire * w ) const

returns wire ID for given FTWire object

Definition at line 89 of file FTFinder.cxx.

89{ return w->getWireId(); }

w

double w

Definition EvtPhokharaDef.hh:28

◆ init()

void FTFinder::init ( )

initializer(creates geometry)

Definition at line 131 of file FTFinder.cxx.

                    {
  if ( !param->_doIt ) return;
 
  // Get the Particle Properties Service
  IPartPropSvc* p_PartPropSvc;
  StatusCode    PartPropStatus = Gaudi::svcLocator()->service( "PartPropSvc", p_PartPropSvc );
 
  if ( !PartPropStatus.isSuccess() )
  {
    std::cerr << "Could not initialize Particle Properties Service" << std::endl;
    return;
  }
  m_particleTable = p_PartPropSvc->PDT();
 
  StatusCode RawData =
      Gaudi::svcLocator()->service( "RawDataProviderSvc", m_rawDataProviderSvc );
  if ( !RawData.isSuccess() )
  {
    std::cerr << "Could not load RawDataProviderSvc!" << m_rawDataProviderSvc << endmsg;
    return;
  }
}

◆ setAlgorithmPointer()

void FTFinder::setAlgorithmPointer ( Algorithm * alg )

returns FTFinder pointer

Definition at line 87 of file FTFinder.cxx.

87{ m_algorithm = alg; }

◆ setBunchtime()

void FTFinder::setBunchtime ( double bunchtime )

inline

set bunchtime for MC events

Definition at line 68 of file FTFinder.h.

68{ _bunchtime = bunchtime; };

◆ setT0cal()

void FTFinder::setT0cal ( double t0cal )

inline

set crude mdc time calibconst (~200 ns) for t0 calculation

Definition at line 70 of file FTFinder.h.

70{ _t0cal = t0cal; };

◆ superLayer()

FTSuperLayer * FTFinder::superLayer ( int id ) const

returns superlayer

Definition at line 85 of file FTFinder.cxx.

85{ return m_geom->getSuperLayer( id ); }

◆ t2x()

float FTFinder::t2x	(	const FTLayer *	l,
		const float	t ) const

convert t to x

Definition at line 95 of file FTFinder.cxx.

                                                           {
  float x = ( t > 0.0f ) ? (( param->_xtCoEff ))*sqrt( t * l->csize() ) : 0.0f;
  if ( x > 0.47f * l->csize() )
  {
    x = 0.0004f * t +
        0.47f * l->csize() *
            ( 1.0f - 0.0004f * 0.47f / ( ( param->_xtCoEff ) * ( param->_xtCoEff ) ) );
  }
 
  return x;
}

◆ term()

void FTFinder::term ( )

terminator

Definition at line 1783 of file FTFinder.cxx.

                    {
  for ( auto t : _tracks ) delete t;
  for ( auto s : _linkedSegments ) delete s;
 
  if ( param->_doIt ) clear();
}

◆ tracks()

FTList< FTTrack * > & FTFinder::tracks ( )

returns track list

Definition at line 86 of file FTFinder.cxx.

86{ return _tracks; }

◆ vertex()

Hep3Vector FTFinder::vertex ( ) const

returns event primary vertex

Definition at line 107 of file FTFinder.cxx.

107{ return Hep3Vector( _vx, _vy, _vz ); }

◆ x2t()

float FTFinder::x2t	(	const FTLayer *	l,
		const float	x ) const

convert x to t

Definition at line 91 of file FTFinder.cxx.

                                                           {
  return ( x * x ) / ( param->_xtCoEff * param->_xtCoEff * l->csize() );
}

Member Data Documentation

◆ eventNo

int FTFinder::eventNo

Definition at line 128 of file FTFinder.h.

Referenced by begin_run().

◆ evtTiming

float FTFinder::evtTiming

Definition at line 133 of file FTFinder.h.

Referenced by event().

◆ expflag

int FTFinder::expflag

Definition at line 130 of file FTFinder.h.

Referenced by begin_run().

◆ i_rPhiFit

int FTFinder::i_rPhiFit

Definition at line 136 of file FTFinder.h.

Referenced by event().

◆ pivot

const HepPoint3D FTFinder::pivot

Definition at line 134 of file FTFinder.h.

Referenced by FTFinder().

◆ runNo

int FTFinder::runNo

Definition at line 129 of file FTFinder.h.

Referenced by begin_run().

◆ t0Estime

float FTFinder::t0Estime

Definition at line 131 of file FTFinder.h.

Referenced by event(), and FTFinder().

◆ tEstFlag

int FTFinder::tEstFlag

Definition at line 137 of file FTFinder.h.

Referenced by event(), and FTFinder().

◆ Testime

float FTFinder::Testime

Definition at line 135 of file FTFinder.h.

◆ tEstime

FTList<FTList<float> > FTFinder::tEstime

Definition at line 138 of file FTFinder.h.

Referenced by event(), and FTFinder().

◆ tOffSet

float FTFinder::tOffSet

Definition at line 132 of file FTFinder.h.

Referenced by FTFinder().

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

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ FTFinder()

Member Function Documentation

◆ begin_run()

◆ event()

◆ getWireId()

◆ init()

◆ setAlgorithmPointer()

◆ setBunchtime()

◆ setT0cal()

◆ superLayer()

◆ t2x()

◆ term()

◆ tracks()

◆ vertex()

◆ x2t()

Member Data Documentation

◆ eventNo

◆ evtTiming

◆ expflag

◆ i_rPhiFit

◆ pivot

◆ runNo

◆ t0Estime

◆ tEstFlag

◆ Testime

◆ tEstime

◆ tOffSet