A tracking module. More...

#include <TrkReco.h>

Inheritance diagram for TrkReco:

Public Member Functions
	TrkReco (const std::string &name, ISvcLocator *pSvcLocator)
	returns TrkReco.

StatusCode	initialize ()

StatusCode	execute ()

StatusCode	finalize ()

StatusCode	beginRun ()

StatusCode	endRun ()

void	disp_stat (const char *)
	initializes TrkReco.

void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	returns name.

const TFinderBase *	confFinder (void) const
	returns a pointer to the rphi finder.

const TFinderBase *	curlFinder (void) const
	returns a pointer to the curl finder.

const TTrackManager &	trackManager (void) const
	returns a pointer to TTrackManager.

void	clear (void)
	clears all TMDC information.

void	fastClear (void)
	clears TMDC information.

const AList< TTrack > &	tracks (void) const
	returns a list of reconstructed tracks.

Public Attributes
int	b_mcPar

int	b_mcHit

int	b_tuple

int	b_goodTrk

int	b_timeTest

bool	m_dropHot

bool	m_combineTracking

bool	m_keepBadTdc

bool	m_keepUnmatch

int	b_RungeKuttaCorrection

float	b_cdcVersion

std::string	_cdcVersion

float	b_fudgeFactor

int	useESTime

int	b_debugLevel

int	b_useAllHits

int	b_doT0Reset

int	b_nT0ResetMax

int	b_doMCAnalysis

float	b_helixFitterChisqMax

int	b_doPerfectFinder

int	b_perfectFitting

int	b_conformalFinder

int	b_doConformalFinder

int	b_doConformalFastFinder

int	b_doConformalSlowFinder

int	b_conformalPerfectSegmentFinding

int	b_conformalFittingFlag

float	b_conformalMaxSigma

int	b_conformalMinNLinksForSegment

int	b_conformalMinNCores

int	b_conformalMinNSegments

float	b_salvageLevel

int	b_conformalSalvageLoadWidth

int	b_conformalStereoMode

int	b_conformalStereoLoadWidth

float	b_conformalStereoMaxSigma

float	b_conformalStereoSzSegmentDistance

float	b_conformalStereoSzLinkDistance

int	b_doConformalFinderStereo

int	b_doConformalFinderCosmic

float	b_conformalFraction

float	b_conformalStereoZ3

float	b_conformalStereoZ4

float	b_conformalStereoChisq3

float	b_conformalStereoChisq4

int	b_conformalFittingCorrections

int	m_CalibFlag

double	b_momentumCut

int	b_fittingFlag

int	b_doSalvage

int	b_doMerge

int	b_doT0Determination

int	b_nTracksForT0

int	b_sortMode

int	b_doAssociation

float	b_associateSigma

int	b_test

int	b_doCurlFinder

int	min_segment

int	min_salvage

double	bad_distance_for_salvage

double	good_distance_for_salvage

int	min_sequence

int	min_fullwire

double	range_for_axial_search

double	range_for_stereo_search

int	superlayer_for_stereo_search

double	range_for_axial_last2d_search

double	range_for_stereo_last2d_search

double	trace2d_distance

double	trace2d_first_distance

double	trace3d_distance

int	determine_one_track

double	selector_max_impact

double	selector_max_sigma

double	selector_strange_pz

double	selector_replace_dz

int	stereo_2dfind

int	merge_exe

double	merge_ratio

double	merge_z_diff

double	mask_distance

double	ratio_used_wire

double	range_for_stereo1

double	range_for_stereo2

double	range_for_stereo3

double	range_for_stereo4

double	range_for_stereo5

double	range_for_stereo6

double	z_cut

double	z_diff_for_last_attend

int	svd_reconstruction

double	min_svd_electrons

int	on_correction

int	output_2dtracks

int	curl_version

double	minimum_seedLength

double	minimum_2DTrackLength

double	minimum_3DTrackLength

double	minimum_closeHitsLength

double	MIN_RADIUS_OF_STRANGE_TRACK

double	ULTIMATE_MIN_RADIUS_OF_STRANGE_TRACK

Detailed Description

A tracking module.

Definition at line 46 of file TrkReco.h.

Constructor & Destructor Documentation

◆ TrkReco()

TrkReco::TrkReco	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

returns TrkReco.

Initiate Parameters for TrkReco

Definition at line 80 of file TrkReco.cxx.

                                                                :
  Algorithm(name, pSvcLocator),
  _cdc(0),
  _perfectFinder(0),
  _rkfitter("range fitter",
      false,
      0,
      true),
  _confFinder(0),
  _curlFinder(0),
  _nEvents(0) {
 
    /// Initiate Parameters for TrkReco
    initPara();
    
    TrkRecoHelixFitterChisqMax = b_helixFitterChisqMax;
}

Member Function Documentation

◆ beginRun()

StatusCode TrkReco::beginRun ( )

Get MdcGeomSvc

Definition at line 302 of file TrkReco.cxx.

                            {
  /// Get MdcGeomSvc
    IMdcGeomSvc* imdcGeomSvc;
    StatusCode sc1 =Gaudi::svcLocator()->service("MdcGeomSvc", imdcGeomSvc);
    _mdcGeomSvc =  dynamic_cast<MdcGeomSvc*> (imdcGeomSvc);
    if (sc1.isFailure()) {
      return( StatusCode::FAILURE);
    }
     _cdc = cdcInit();
 
    return StatusCode::SUCCESS;  
}

◆ clear()

void TrkReco::clear ( void )

clears all TMDC information.

Definition at line 956 of file TrkReco.cxx.

                   {
 
  //...Clear track candidates of the last event...
  HepAListDeleteAll(_mcTracks);
 
  //...Clear finders...
  if (b_doPerfectFinder) _perfectFinder->clear();
  if (b_doConformalFinder) _confFinder->clear();
  if (b_doCurlFinder) _curlFinder->clear();
  _trackManager.clear();
  _trackManager.clearTables();
}

Referenced by execute(), and fastClear().

◆ confFinder()

const TFinderBase * TrkReco::confFinder ( void ) const

inline

returns a pointer to the rphi finder.

Definition at line 320 of file TrkReco.h.

                              {
    return (TFinderBase *) _confFinder;
}

◆ curlFinder()

const TFinderBase * TrkReco::curlFinder ( void ) const

inline

returns a pointer to the curl finder.

Definition at line 326 of file TrkReco.h.

                              {
    return (TFinderBase *) _curlFinder;
}

◆ disp_stat()

void TrkReco::disp_stat ( const char * m )

initializes TrkReco.

terminates TrkReco. processes an event. temporarily required

Definition at line 315 of file TrkReco.cxx.

                                 {
    std::string msg = m;
//    dump(msg);
}

◆ dump()

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

returns name.

returns version. dumps debug information.

◆ endRun()

StatusCode TrkReco::endRun ( )

inline

Definition at line 56 of file TrkReco.h.

56{};

◆ execute()

StatusCode TrkReco::execute ( )

Definition at line 321 of file TrkReco.cxx.

                            {
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in execute()" << endreq;
  StatusCode sc;   
  if (b_timeTest && b_tuple) m_timer[1]->start();
 
  // Initiate state of all sense wire
  if (b_goodTrk && b_tuple)
    for (int ii=0;ii<43;ii++){
      for (int jj=0;jj<288;jj++){
    havedigi[ii][jj]= -99;//no hit/noise
      }
    }
  
 //------------------------------------
  // Initialize track collection in TDS
  //------------------------------------
  //yzhang add 2010-04-30 
  //Clear TDS tracks and hits
  if(!b_RungeKuttaCorrection){
  IDataManagerSvc *dataManSvc;
  DataObject *aTrackCol;
  DataObject *aRecHitCol;
  if(!m_combineTracking){
    dataManSvc = dynamic_cast<IDataManagerSvc*> (eventSvc().get());
    eventSvc()->findObject("/Event/Recon/RecMdcTrackCol",aTrackCol);
    if(aTrackCol != NULL) {
      dataManSvc->clearSubTree("/Event/Recon/RecMdcTrackCol");
      eventSvc()->unregisterObject("/Event/Recon/RecMdcTrackCol");
    }
    eventSvc()->findObject("/Event/Recon/RecMdcHitCol",aRecHitCol);
    if(aRecHitCol != NULL) {
      dataManSvc->clearSubTree("/Event/Recon/RecMdcHitCol");
      eventSvc()->unregisterObject("/Event/Recon/RecMdcHitCol");
    }
  }
  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" <<endreq;
      return( StatusCode::FAILURE);
    }
  }
  RecMdcTrackCol* trkcol;
  eventSvc()->findObject("/Event/Recon/RecMdcTrackCol",aTrackCol);
  if (aTrackCol) {
    trkcol = dynamic_cast<RecMdcTrackCol*> (aTrackCol);
  }else{
    trkcol= new RecMdcTrackCol;
    sc =  eventSvc()->registerObject(EventModel::Recon::RecMdcTrackCol, trkcol);
    if(!sc.isSuccess()) {
      log << MSG::FATAL << " Could not register RecMdcTrack collection" <<endreq;
      return StatusCode::FAILURE;
    }
  }
  RecMdcHitCol* hitcol;
  eventSvc()->findObject("/Event/Recon/RecMdcHitCol",aRecHitCol);
  if (aRecHitCol) {
    hitcol = dynamic_cast<RecMdcHitCol*> (aRecHitCol);
  }else{
    hitcol= new RecMdcHitCol;
    sc =  eventSvc()->registerObject(EventModel::Recon::RecMdcHitCol, hitcol);
    if(!sc.isSuccess()) {
     log << MSG::FATAL << " Could not register RecMdcHit collection" <<endreq;
      return StatusCode::FAILURE;
    }
  }  
  // Part 1: Get Tev
  t0_bes = 0.;
  t0Sta  = -1;
  if (useESTime) { 
    SmartDataPtr<RecEsTimeCol> aevtimeCol(eventSvc(),"/Event/Recon/RecEsTimeCol");
//bugcheck jialk    
    if (aevtimeCol  && aevtimeCol->size() ) {
      RecEsTimeCol::iterator iter_evt = aevtimeCol->begin();
      t0_bes =  (*iter_evt)->getTest();
      t0Sta  =  (*iter_evt)->getStat();
    }else{
      log << MSG::WARNING << "Could not find EsTimeCol" << endreq;
      return StatusCode::SUCCESS;
    }
  }
 
  // Part 2: Get the event header
  SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
  if (!eventHeader) {
    log << MSG::FATAL << "Could not find Event Header" << endreq;
    return( StatusCode::FAILURE);
  }
  _nEvents = eventHeader->eventNumber();
  if (b_tuple) std::cout << "TrkReco ... processing ev# " << _nEvents << std::endl;
 
  // Part 3: Retrieve MDC digi 
  int digiId;
  uint32_t getDigiFlag = 0;
  if(m_dropHot || m_combineTracking)getDigiFlag |= MdcRawDataProvider::b_dropHot;
  if(m_keepBadTdc)  getDigiFlag |= MdcRawDataProvider::b_keepBadTdc;
  if(m_keepUnmatch) getDigiFlag |= MdcRawDataProvider::b_keepUnmatch;
  MdcDigiVec mdcDigiVec = _rawDataProviderSvc->getMdcDigiVec(getDigiFlag);
  if (0 == mdcDigiVec.size()){
    log << MSG::WARNING << " No hits in MdcDigiVec" << endreq;
    return StatusCode::SUCCESS;
  }
/*  //jialk in order to reject events with exceeding # of total hits   2009/03/31   
if (mdcDigiVec.size() > 2000){   
     log << MSG::WARNING << " Too many hits in MdcDigiVec" << endreq;   
     return StatusCode::SUCCESS;   
   }*/ 
  MdcDigiVec::iterator iter1 = mdcDigiVec.begin();
  MC_DIGI_SIZE = mdcDigiVec.size();
  
  digiId = 0; 
  Identifier mdcId;
  int layerId;
  int wireId;
 
  //Clear the old MdcRec_wirhit tables and create the hits' info for the new event.
  unsigned nt = MdcRecWirhitCol::getMdcRecWirhitCol()->size();
  //cout<<"Col size of last Event's WirHit = "<<nt<<endl;
  MdcRecWirhitCol::getMdcRecWirhitCol()->clear();
 
  for (;iter1 != mdcDigiVec.end(); iter1++, digiId++) {
    //log << MSG::INFO << "MDC digit No: " << digiId << endreq;
    mdcId = (*iter1)->identify();
    layerId = MdcID::layer(mdcId);
    wireId = MdcID::wire(mdcId);
    /*log << MSG::INFO
      << " time_channel = " << (*iter1)->getTimeChannel()
      <<" time = "<<(*iter1)->getTimeChannel()* 40./10000/1000000
      << " charge_channel = " << (*iter1)->getChargeChannel()
      << " layerId = " << layerId 
      << " wireId = " << wireId
      << endreq;*/
 
    if (b_goodTrk && b_tuple) havedigi[layerId][wireId] = (*iter1)->getTrackIndex();  //-1:noise  0-n:tracks
 
        MdcRec_wirhit mhit;
 
        mhit.id = digiId;
    mhit.geo = _mdcGeomSvc->Wire(layerId,wireId);
    //Apply crude TOF, Belle: tof=1.074*radius/c, here we use 1.28 instead of 1.074. 
    double tof;
    tof = 1.28 * mhit.geo->Lyr()->Radius()/299.8;  //unit of Radius is mm.  
    mhit.tdc = RawDataUtil::MdcTime((*iter1)->getTimeChannel());
 
    //cout<<" ..  mhit.tdc = "<<mhit.tdc<<"    t0_bes = "<<t0_bes<<endl;
    mhit.tdc -= t0_bes;
 
    //jialk 
    _trackManager.sett0bes(t0_bes);
    //mhit.adc = RawDataUtil::MdcCharge((*iter1)->getChargeChannel());
    //mhit.tdc = (*iter1)->getTimeChannel();
    mhit.adc = (*iter1)->getChargeChannel();
    //cout<<"raw tdc(ns): "<<mhit.tdc<<"; tof(ns): "<<tof<<endl;
//  double dist2 = _mdcCalibFunSvc->rawTimeNoTOFToDist(mhit.tdc-tof, layerId, wireId, 2, 0.0);
 //       mhit.erddl = _mdcCalibFunSvc->getSigma(layerId, 2, dist2, 0.0);
        double timewalk=_mdcCalibFunSvc->getTimeWalk(layerId,mhit.adc); 
        double T0 = _mdcCalibFunSvc->getT0(layerId,wireId);
        double drifttime = mhit.tdc-tof-timewalk-T0;
        double dist2 = _mdcCalibFunSvc->driftTimeToDist(drifttime,layerId, wireId, 2, 0);//by liucy 2010/05/12
        mhit.erddl = _mdcCalibFunSvc->getSigma(layerId, 2 , dist2, 0.0,0.0,0.0,mhit.adc);
//cout<<"getSigma: "<<mhit.erddl<<endl;
    mhit.ddl = dist2/10.;  //mm->cm
        mhit.ddr = mhit.ddl;
    mhit.erddl = mhit.erddl/10.;  //mm->cm
        mhit.erddr = mhit.erddl;
 
        mhit.lr = 2;
        mhit.stat = 0;
        mhit.stat = mhit.stat |= 1048576; //bit20  WireHitTimeValid
        mhit.stat = mhit.stat |= 2097152; //bit21  WireHitChargeValid
        mhit.stat = mhit.stat |= 4194304; //bit22  WireHitFindingValid
        mhit.stat = mhit.stat |= 1073741824; //bit30 
    //cout<<"layerNo = "<<mhit.geo->Layer()<<"; "<<mdigi.digi[i].layerNo<<endl;
    //cout<<"cellNo = "<<mhit.geo->Cell()<<"; "<<mdigi.digi[i].cellNo<<endl;
    //cout<<"NCell of this layer = "<<mhit.geo->Lyr()->NCell()<<endl;
    //cout<<"NCell of this layer = "<<mhit.geo->Lyr()->NCell()<<endl;
        MdcRecWirhitCol::getMdcRecWirhitCol()->push_back(mhit);
 
    if (b_tuple) {
    t10_tdc = mhit.tdc;
    t10_adc = mhit.adc;
    t10_drift = mhit.ddl;
    t10_dDrift = mhit.erddl;
    t10_lyrId = layerId;
    t10_localId = wireId;
    m_tuple10->write();
    }
   }
 
    unsigned nT0Reset = 0;
 
    //...Starting point...
TrkReco_start:
 
    //...Clear myself...
    clear();
    
    //jialk in order to reject events with exceeding # of total hits   2009/03/31   
    if (mdcDigiVec.size() > 2000){   
     log << MSG::WARNING << " Too many hits in MdcDigiVec" << endreq;   
     return StatusCode::SUCCESS;   
    }
 
    //...Update TMDC...
    _cdc->update(b_doMCAnalysis);
 
    //...Get lists of hits...
    unsigned mask = 0;
    if (! b_useAllHits) mask = WireHitFindingValid;
    const AList<TMDCWireHit> & axialHits = _cdc->axialHits(mask);
    const AList<TMDCWireHit> & stereoHits = _cdc->stereoHits(mask);
    const AList<TMDCWireHit> & allHits = _cdc->hits(mask);
    //cout<<"axial: "<<axialHits.length()<<"  stereo: "<<stereoHits.length()<<endl; 
    
    //...Storage for tracks...
    AList<TTrack> tracks;
    AList<TTrack> tracks2D;
 
    //...Perfect finder...
    if (b_doPerfectFinder) {
      _perfectFinder->doit(axialHits, stereoHits, tracks, tracks2D);
      _trackManager.append(tracks);
    }
 
    else {
    //...Conformal finder...
    if (b_doConformalFinder) {
 
        //...T0 reset option...
        if (b_doT0Reset) {
        if (b_nT0ResetMax > nT0Reset)
            ((TConformalFinder *) _confFinder)->doT0Reset(true);
        else
            ((TConformalFinder *) _confFinder)->doT0Reset(false);
        }           
 
        _confFinder->doit(axialHits, stereoHits, tracks, tracks2D);
 
        //...T0 reset...
        if (b_doT0Reset) {
        ++nT0Reset;
        if (((TConformalFinder *) _confFinder)->T0ResetDone())
            goto TrkReco_start;
        }
 
        //cout<<"tracks: "<<tracks.length()<<endl;
        
        //...Stores tracks...
        _trackManager.append(tracks);
        _trackManager.append2D(tracks2D);
        if (b_conformalFinder == 0) {
        if (b_doSalvage == 1) _trackManager.salvage(allHits);
        if (b_doSalvage) _trackManager.mask();
        }
    }
    
    //...Curl finder...
        if (b_doCurlFinder) {
        if ((! b_doSalvage) && (b_conformalFinder == 0))
        _trackManager.maskCurlHits(axialHits,
                       stereoHits,
                       _trackManager.tracks());
        AList<TTrack> confTracks = _trackManager.tracks();
        tracks.append(confTracks);
        _curlFinder->doit(axialHits, stereoHits, tracks, tracks2D);
            tracks.remove(confTracks);
        //_trackManager.append(tracks);
    }
 
    //...Finishes tracks...
//  if ((b_doSalvage) && (b_conformalFinder == 0)) _trackManager.refit();   
 
        //...Appends tracks which are reconstructed by CurlFinder...
    if (b_doCurlFinder) {
      _trackManager.append(tracks);
      _trackManager.append2D(tracks2D);
    }
       
    //...Merge & Mask ...
    //if ((b_doMerge) && (b_conformalFinder == 0)) _trackManager.merge();
    _trackManager.merge();  //Liuqg
    //if (b_conformalFinder != 0) _trackManager.setCurlerFlags();
 
    //...Salvage for dE/dx...
    //if (b_doAssociation)
    //  _trackManager.salvageAssociateHits(allHits, b_associateSigma);
    }
 
    //...Move a pivot... //move to the innermost hit, remove this step now 2005/10/17 zang shilei
    //_trackManager.movePivot();
 
    //...Save Panther tables...
    // _trackManager.checkNumberOfHits();
    //if (b_sortMode == 0) _trackManager.sortTracksByQuality();
    //else                 _trackManager.sortTracksByPt();
    /*if (b_doMCAnalysis) {
      if (mcEvent()) {
    mcInformation();
    _trackManager.saveMCTables();
      }
    }*/
 
    if (b_tuple) FillTuple();
 
    //write Tds 
    int t_tkStat = 2;
    if (!b_doConformalFinder && b_doCurlFinder) t_tkStat = 3;//FIXME
    StatusCode scTds = _trackManager.makeTds(trkcol, hitcol, t_tkStat,b_RungeKuttaCorrection,m_CalibFlag);
    if (scTds !=  StatusCode::SUCCESS) return( StatusCode::FAILURE);
 
    //...T0 correction...
    if (b_doT0Determination) {
      _trackManager.determineT0(b_doT0Determination, b_nTracksForT0);
      if (b_tuple) t3_t0Rec = _trackManager.paraT0();
    }
  }
  else if(b_RungeKuttaCorrection==1){
    if (useESTime) { 
      SmartDataPtr<RecEsTimeCol> aevtimeCol(eventSvc(),"/Event/Recon/RecEsTimeCol");
      if (aevtimeCol  && aevtimeCol->size() ) {
    RecEsTimeCol::iterator iter_evt = aevtimeCol->begin();
    t0_bes =  (*iter_evt)->getTest();
    t0Sta  =  (*iter_evt)->getStat();
      }else{
    log << MSG::WARNING << "Could not find EsTimeCol" << endreq;
    return StatusCode::SUCCESS;
      }
    }
 
    SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
    if (!eventHeader) {
      log << MSG::FATAL << "Could not find Event Header" << endreq;
      return( StatusCode::FAILURE);
    }
    _nEvents = eventHeader->eventNumber();
    if (b_tuple) std::cout << "TrkReco ... processing ev# " << _nEvents << std::endl;
    AList<TTrack> rktracks;
    int digiId;
    uint32_t getDigiFlag = 0;
    if(m_dropHot || m_combineTracking)getDigiFlag |= MdcRawDataProvider::b_dropHot;
    if(m_keepBadTdc)  getDigiFlag |= MdcRawDataProvider::b_keepBadTdc;
    if(m_keepUnmatch) getDigiFlag |= MdcRawDataProvider::b_keepUnmatch;
    MdcDigiVec mdcDigiVec = _rawDataProviderSvc->getMdcDigiVec(getDigiFlag);
    if (0 == mdcDigiVec.size()){
      log << MSG::WARNING << " No hits in MdcDigiVec" << endreq;
      return StatusCode::SUCCESS;
    }
    MdcDigiVec::iterator iter1 = mdcDigiVec.begin();
    MC_DIGI_SIZE = mdcDigiVec.size();
 
    digiId = 0;
    Identifier mdcId;
    int layerId;
    int wireId;
 
    unsigned nt = MdcRecWirhitCol::getMdcRecWirhitCol()->size();
    MdcRecWirhitCol::getMdcRecWirhitCol()->clear();
 
    for (;iter1 != mdcDigiVec.end(); iter1++, digiId++) {
      mdcId = (*iter1)->identify();
      layerId = MdcID::layer(mdcId);
      wireId = MdcID::wire(mdcId);
      if (b_goodTrk && b_tuple) havedigi[layerId][wireId] = (*iter1)->getTrackIndex();
 
      MdcRec_wirhit mhit;
 
      mhit.id = digiId;
      mhit.geo = _mdcGeomSvc->Wire(layerId,wireId);
      double tof;
      tof = 1.28 * mhit.geo->Lyr()->Radius()/299.8;
      mhit.tdc = RawDataUtil::MdcTime((*iter1)->getTimeChannel());
      mhit.timechannel=(*iter1)->getTimeChannel();
      mhit.tdc -= t0_bes;
 
      _trackManager.sett0bes(t0_bes);
      mhit.adc = (*iter1)->getChargeChannel();
      double timewalk=_mdcCalibFunSvc->getTimeWalk(layerId,mhit.adc);
      double T0 = _mdcCalibFunSvc->getT0(layerId,wireId);
      double drifttime = mhit.tdc-tof-timewalk-T0;
      double dist2 = _mdcCalibFunSvc->driftTimeToDist(drifttime,layerId, wireId, 2, 0);
      mhit.erddl = _mdcCalibFunSvc->getSigma(layerId, 2 , dist2, 0.0,0.0,0.0,mhit.adc);
      mhit.ddl = dist2/10.;
      mhit.ddr = mhit.ddl;
      mhit.erddl = mhit.erddl/10.;
      mhit.erddr = mhit.erddl;
 
      mhit.lr = 2;
      mhit.stat = 0;
      mhit.stat = mhit.stat |= 1048576;
      mhit.stat = mhit.stat |= 2097152;
      mhit.stat = mhit.stat |= 4194304;
      mhit.stat = mhit.stat |= 1073741824;
      MdcRecWirhitCol::getMdcRecWirhitCol()->push_back(mhit);
 
      if (b_tuple) {
    t10_tdc = mhit.tdc;
    t10_adc = mhit.adc;
    t10_drift = mhit.ddl;
    t10_dDrift = mhit.erddl;
    t10_lyrId = layerId;
    t10_localId = wireId;
    m_tuple10->write();
      }
    }
 
    unsigned nT0Reset = 0;
 
    clear();
    if (mdcDigiVec.size() > 2000){
      log << MSG::WARNING << " Too many hits in MdcDigiVec" << endreq;
      return StatusCode::SUCCESS;
    }
    int counter=0;
    _cdc->update(b_doMCAnalysis);
 
    unsigned mask = 0;
    if (! b_useAllHits) mask = WireHitFindingValid;
    const AList<TMDCWireHit> & axialHits = _cdc->axialHits(mask);
    const AList<TMDCWireHit> & stereoHits = _cdc->stereoHits(mask);
    const AList<TMDCWireHit> & allHits = _cdc->hits(mask);
    AList<TMLink> _allHits[3];
    TConformalFinder0::conformalTransformationDriftCircle(ORIGIN, axialHits, _allHits[0]);
    TConformalFinder0::conformalTransformationDriftCircle(ORIGIN, stereoHits, _allHits[1]);
    _allHits[2].append(_allHits[0]);
    _allHits[2].append(_allHits[1]);
 
 
   SmartDataPtr<RecMdcTrackCol> mdcTracks(eventSvc(),EventModel::Recon::RecMdcTrackCol);
   if( ! mdcTracks )
   {
     log << MSG::ERROR << "Unable to retrieve RecMdcTrackCol" << endreq;
     return  StatusCode::FAILURE;
   } else {
     log << MSG::DEBUG << "RecMdcTrackCol retrieved of size "<< mdcTracks->size() << endreq;
     int ntrk=0;
     for(RecMdcTrackCol::iterator it=mdcTracks->begin(); it!=mdcTracks->end(); it++)
     {ntrk++;
       TRunge *r = new TRunge(**it);
       TTrack *t1 = new TTrack(*r);
       HitRefVec gothits = (*it)->getVecHits();
       HitRefVec::iterator it_gothit = gothits.begin();
       unsigned stat=(*it)->stat();
       int nhit=0;
       for( ; it_gothit != gothits.end(); it_gothit++){
         for(int i=0;i<_allHits[2].length();i++){
           int lyrraw=_allHits[2][i]->wire()->layerId();
           int wireraw=_allHits[2][i]->wire()->localId();
           int g_layer       = MdcID::layer((*it_gothit)->getMdcId());
           int g_wire        = MdcID::wire((*it_gothit)->getMdcId());
           if(lyrraw==g_layer&&g_wire==wireraw){
             nhit++;
             t1->append(*_allHits[2][i]);
           }
         }
       }
       TRunge *rr = new TRunge(*t1);
       TRunge *tt = new TRunge(*t1);
       int err=_rkfitter.fit(*rr);
       int nhits=rr->links().length();
       int ndrop=0;
       int nmax=0;
       if(nhits<=10)nmax=0;
       if(nhits>10)nmax=(int)nhit*0.3;
       for(int ii=0;ii<nmax;ii++){ 
    
        ndrop=rr->DropWorst();
        if(ndrop)err=_rkfitter.fit(*rr);
        if(err)break;
       }  
       if(err==-2) counter++;
       if(m_CalibFlag==1) {//This is for Calibration 
         tt->removeLinks();
         tt->append(rr->links());
         for(int i=0;i<43;i++){
            err= _rkfitter.fit(*tt,0,i);
         }
        rr->removeLinks();
        rr->append(tt->links());
        TTrack *t =new TTrack(*rr);
        if(err==0){
            rktracks.append(*t); 
            t->setFinderType(100+stat);
        }
       }
       if(m_CalibFlag==0){
       TTrack *t =new TTrack(*rr);
       if(err==0){
           
        t->setFinderType(100+stat);
         rktracks.append(*t);
       }
       }
       delete r ;
       delete t1;
       delete rr;
       delete tt;
     }
   }
 
  _trackManager.append(rktracks);
  IDataManagerSvc *dataManSvc;
  DataObject *aTrackCol;
  DataObject *aRecHitCol;
  if(!m_combineTracking){
    dataManSvc = dynamic_cast<IDataManagerSvc*> (eventSvc().get());
    eventSvc()->findObject("/Event/Recon/RecMdcTrackCol",aTrackCol);
    if(aTrackCol != NULL) {
      dataManSvc->clearSubTree("/Event/Recon/RecMdcTrackCol");
      eventSvc()->unregisterObject("/Event/Recon/RecMdcTrackCol");
    }
    eventSvc()->findObject("/Event/Recon/RecMdcHitCol",aRecHitCol);
    if(aRecHitCol != NULL) {
      dataManSvc->clearSubTree("/Event/Recon/RecMdcHitCol");
      eventSvc()->unregisterObject("/Event/Recon/RecMdcHitCol");
    }
  }
  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" <<endreq;
      return( StatusCode::FAILURE);
    }
  }
  RecMdcTrackCol* trkcol;
  eventSvc()->findObject("/Event/Recon/RecMdcTrackCol",aTrackCol);
  if (aTrackCol) {
    trkcol = dynamic_cast<RecMdcTrackCol*> (aTrackCol);
  }else{
    trkcol= new RecMdcTrackCol;
    sc =  eventSvc()->registerObject(EventModel::Recon::RecMdcTrackCol, trkcol);
    if(!sc.isSuccess()) {
      log << MSG::FATAL << " Could not register RecMdcTrack collection" <<endreq;
      return StatusCode::FAILURE;
    }
  }
  RecMdcHitCol* hitcol;
  eventSvc()->findObject("/Event/Recon/RecMdcHitCol",aRecHitCol);
  if (aRecHitCol) {
    hitcol = dynamic_cast<RecMdcHitCol*> (aRecHitCol);
  }else{
    hitcol= new RecMdcHitCol;
    sc =  eventSvc()->registerObject(EventModel::Recon::RecMdcHitCol, hitcol);
    if(!sc.isSuccess()) {
      log << MSG::FATAL << " Could not register RecMdcHit collection" <<endreq;
      return StatusCode::FAILURE;
    }
  }
  if (b_tuple) FillTuple();
 
  int t_tkStat = 2;
  if (!b_doConformalFinder && b_doCurlFinder) t_tkStat = 3;//FIXME
  StatusCode scTds = _trackManager.makeTds(trkcol, hitcol, t_tkStat,b_RungeKuttaCorrection,m_CalibFlag);
  if (b_doT0Determination) {
    _trackManager.determineT0(b_doT0Determination, b_nTracksForT0);
    if (b_tuple) t3_t0Rec = _trackManager.paraT0();
  }
   for (unsigned i = 0; i < 3; i++) {
        if (i == 2)
            HepAListDeleteAll(_allHits[i]);
        else
            _allHits[i].removeAll();
    }
  
  rktracks.removeAll(); 
  }
  //...For debug...
  /*if (b_debugLevel) {
    std::cout << "TrkReco ... ev# " << _nEvents << " processed,"
    << " #tracks found=" << _trackManager.allTracks().length()
    << ", #good tracks=" << _trackManager.tracks().length()
    << ", #2D tracks=" << _trackManager.tracks2D().length()
    << std::endl;
    if (b_debugLevel > 1) _trackManager.dump("eventSummary hits");
    else _trackManager.dump("eventSummary");
    }*/
 
  //TUpdater::update();
  return StatusCode::SUCCESS;
}

◆ fastClear()

void TrkReco::fastClear ( void )

clears TMDC information.

Definition at line 970 of file TrkReco.cxx.

                       {
  clear();
}

◆ finalize()

StatusCode TrkReco::finalize ( )

Definition at line 284 of file TrkReco.cxx.

                             {
 
    if (b_timeTest && b_tuple) m_timersvc->print();
 
    //...Clear TrkReco objects...
    if (_cdc) _cdc->fastClear();
    if (b_doConformalFinder && _confFinder) _confFinder->clear();
    if (b_doCurlFinder && _curlFinder) _curlFinder->clear();
 
    //...Summary...
//    dump("summary");
 
    //...Delete TrkReco objects...
    if (_confFinder) delete _confFinder;
    if (_curlFinder) delete _curlFinder;
 
    return StatusCode::SUCCESS;
}

◆ initialize()

StatusCode TrkReco::initialize ( )

Get RawDataProviderSvc

Get MdcGeomSvc

Get MdcCalibFunSvc

Initialize NTuple

Definition at line 99 of file TrkReco.cxx.

                              {
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "in initialize()" << endreq;
 
    StatusCode sc;
 
    /// Get RawDataProviderSvc
    IRawDataProviderSvc* irawDataProviderSvc;
    sc = service ("RawDataProviderSvc", irawDataProviderSvc);
    _rawDataProviderSvc = dynamic_cast<RawDataProviderSvc*> (irawDataProviderSvc);
    if ( sc.isFailure() ){
      log << MSG::FATAL << "Could not load RawDataProviderSvc!" << endreq;
      return StatusCode::FAILURE;
    }
 
    /// Get MdcGeomSvc
 //   IMdcGeomSvc* imdcGeomSvc;
 //   sc = service("MdcGeomSvc", imdcGeomSvc);
 //   _mdcGeomSvc =  dynamic_cast<MdcGeomSvc*> (imdcGeomSvc);
 //   if (sc.isFailure()) {
 //     return( StatusCode::FAILURE);  
 //   }
 
    /// Get MdcCalibFunSvc
    IMdcCalibFunSvc* imdcCalibSvc;
    sc = service ("MdcCalibFunSvc",imdcCalibSvc);
    _mdcCalibFunSvc = dynamic_cast<MdcCalibFunSvc*> (imdcCalibSvc);
    if ( sc.isFailure() ){
      log << MSG::FATAL << "Could not load MdcCalibFunSvc!" << endreq;
    }
 
    /// Initialize NTuple
    if (b_tuple) InitTuple();
 
#ifdef TRKRECO_DEBUG
    b_debugLevel = 2;
#endif
    
    //...Create TMDC...
//    _cdc = cdcInit();
 
    //...TrkReco...
    TrkRecoHelixFitterChisqMax = b_helixFitterChisqMax;
 
    //...Create rphi finder...
/*    _perfectFinder = new TPerfectFinder(b_perfectFitting,
                    b_conformalMaxSigma,
                    b_conformalStereoMaxSigma,
                    b_conformalFittingCorrections);
*/
    if ((! _confFinder) && (b_conformalFinder == 0)) {
       _confFinder = new TConformalFinder0(b_conformalMaxSigma,
                        b_conformalFraction,
                        b_conformalStereoZ3,
                        b_conformalStereoZ4,
                        b_conformalStereoChisq3,
                        b_conformalStereoChisq4,
                        b_conformalStereoMaxSigma,
                        b_conformalFittingCorrections,
                        b_salvageLevel,
                        b_doConformalFinderCosmic);
    }
    else if (! _confFinder) {
    _confFinder = new TConformalFinder(b_doConformalFastFinder,
                       b_doConformalSlowFinder,
                       b_conformalPerfectSegmentFinding,
                       b_conformalMaxSigma,
                       b_conformalStereoMaxSigma,
                       b_salvageLevel,
                       b_conformalMinNLinksForSegment,
                       b_conformalMinNCores,
                       b_conformalMinNSegments,
                       b_conformalSalvageLoadWidth,
                       b_conformalStereoMode,
                       b_conformalStereoLoadWidth,
                       b_conformalStereoSzSegmentDistance,
                       b_conformalStereoSzLinkDistance,
                       b_conformalFittingFlag);
    }
    _confFinder->debugLevel(b_debugLevel);
    _confFinder->doStereo(b_doConformalFinderStereo);
 
    //...Salvage flag is ignored in new conf. finder...
    if (b_doSalvage == 2) _confFinder->doSalvage(true);
 
    //...Create curl finder...
    if (! _curlFinder)
        _curlFinder = new TCurlFinder((unsigned)min_segment,
                      (unsigned)min_salvage,
                      bad_distance_for_salvage,
                      good_distance_for_salvage,
                      (unsigned)min_sequence,
                      (unsigned)min_fullwire,
                      range_for_axial_search,
                      range_for_stereo_search,
                      (unsigned)superlayer_for_stereo_search,
                      range_for_axial_last2d_search,
                      range_for_stereo_last2d_search,
                      trace2d_distance,
                      trace2d_first_distance,
                      trace3d_distance,
                      (unsigned)determine_one_track,
                      selector_max_impact,
                      selector_max_sigma,
                      selector_strange_pz,
                      selector_replace_dz,
                      (unsigned)stereo_2dfind,
                      (unsigned)merge_exe,
                      merge_ratio,
                      merge_z_diff,
                      mask_distance,
                      ratio_used_wire,
                      range_for_stereo1,
                      range_for_stereo2,
                      range_for_stereo3,
                      range_for_stereo4,
                      range_for_stereo5,
                      range_for_stereo6,
                      z_cut,
                      z_diff_for_last_attend,
                      (unsigned)svd_reconstruction,
                      min_svd_electrons,
                      (unsigned)on_correction,
                      (unsigned)output_2dtracks,
                      (unsigned)curl_version,
                      //jialk
                      minimum_seedLength, 
                      minimum_2DTrackLength,
                      minimum_3DTrackLength,
                      minimum_closeHitsLength,
                      MIN_RADIUS_OF_STRANGE_TRACK, 
                      ULTIMATE_MIN_RADIUS_OF_STRANGE_TRACK);
    _curlFinder->debugLevel(b_debugLevel);
 
    //...Track manager setup...
    _trackManager.maxMomentum(b_momentumCut);
    _trackManager.debugLevel(b_debugLevel);
    _trackManager.fittingFlag(b_fittingFlag);
 
    //...Initialize...
//zsl    TUpdater::initialize();
 
    //...For debug...
//    dump("parameter");
 
    TrkRecoTest = b_test;
 
    if (b_timeTest && b_tuple) {
      StatusCode sc = service( "BesTimerSvc", m_timersvc);
      if( sc.isFailure() ) {
    log << MSG::WARNING << " Unable to locate BesTimerSvc" << endreq;
    return StatusCode::FAILURE;
      } 
      m_timer[1] = m_timersvc->addItem("Execution");
      m_timer[1]->propName("nExecution");
    }
    
    return StatusCode::SUCCESS;   
}

◆ trackManager()

const TTrackManager & TrkReco::trackManager ( void ) const

inline

returns a pointer to TTrackManager.

Definition at line 346 of file TrkReco.h.

                                {
    return _trackManager;
}

◆ tracks()

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

inline

returns a list of reconstructed tracks.

Definition at line 332 of file TrkReco.h.

                          {
    return _trackManager.tracksFinal();
}

Referenced by execute().

Member Data Documentation

◆ _cdcVersion

std::string TrkReco::_cdcVersion

Definition at line 138 of file TrkReco.h.

◆ b_associateSigma

float TrkReco::b_associateSigma

Definition at line 195 of file TrkReco.h.

◆ b_cdcVersion

float TrkReco::b_cdcVersion

Definition at line 137 of file TrkReco.h.

◆ b_conformalFinder

int TrkReco::b_conformalFinder

Definition at line 157 of file TrkReco.h.

Referenced by execute(), and initialize().

◆ b_conformalFittingCorrections

int TrkReco::b_conformalFittingCorrections

Definition at line 183 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalFittingFlag

int TrkReco::b_conformalFittingFlag

Definition at line 162 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalFraction

float TrkReco::b_conformalFraction

Definition at line 178 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalMaxSigma

float TrkReco::b_conformalMaxSigma

Definition at line 163 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalMinNCores

int TrkReco::b_conformalMinNCores

Definition at line 165 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalMinNLinksForSegment

int TrkReco::b_conformalMinNLinksForSegment

Definition at line 164 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalMinNSegments

int TrkReco::b_conformalMinNSegments

Definition at line 166 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalPerfectSegmentFinding

int TrkReco::b_conformalPerfectSegmentFinding

Definition at line 161 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalSalvageLoadWidth

int TrkReco::b_conformalSalvageLoadWidth

Definition at line 168 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoChisq3

float TrkReco::b_conformalStereoChisq3

Definition at line 181 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoChisq4

float TrkReco::b_conformalStereoChisq4

Definition at line 182 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoLoadWidth

int TrkReco::b_conformalStereoLoadWidth

Definition at line 170 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoMaxSigma

float TrkReco::b_conformalStereoMaxSigma

Definition at line 171 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoMode

int TrkReco::b_conformalStereoMode

Definition at line 169 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoSzLinkDistance

float TrkReco::b_conformalStereoSzLinkDistance

Definition at line 173 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoSzSegmentDistance

float TrkReco::b_conformalStereoSzSegmentDistance

Definition at line 172 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoZ3

float TrkReco::b_conformalStereoZ3

Definition at line 179 of file TrkReco.h.

Referenced by initialize().

◆ b_conformalStereoZ4

float TrkReco::b_conformalStereoZ4

Definition at line 180 of file TrkReco.h.

Referenced by initialize().

◆ b_debugLevel

int TrkReco::b_debugLevel

Definition at line 145 of file TrkReco.h.

Referenced by initialize().

◆ b_doAssociation

int TrkReco::b_doAssociation

Definition at line 194 of file TrkReco.h.

◆ b_doConformalFastFinder

int TrkReco::b_doConformalFastFinder

Definition at line 159 of file TrkReco.h.

Referenced by initialize().

◆ b_doConformalFinder

int TrkReco::b_doConformalFinder

Definition at line 158 of file TrkReco.h.

Referenced by clear(), execute(), and finalize().

◆ b_doConformalFinderCosmic

int TrkReco::b_doConformalFinderCosmic

Definition at line 177 of file TrkReco.h.

Referenced by initialize().

◆ b_doConformalFinderStereo

int TrkReco::b_doConformalFinderStereo

Definition at line 176 of file TrkReco.h.

Referenced by initialize().

◆ b_doConformalSlowFinder

int TrkReco::b_doConformalSlowFinder

Definition at line 160 of file TrkReco.h.

Referenced by initialize().

◆ b_doCurlFinder

int TrkReco::b_doCurlFinder

Definition at line 199 of file TrkReco.h.

Referenced by clear(), execute(), and finalize().

◆ b_doMCAnalysis

int TrkReco::b_doMCAnalysis

Definition at line 149 of file TrkReco.h.

Referenced by execute().

◆ b_doMerge

int TrkReco::b_doMerge

Definition at line 190 of file TrkReco.h.

◆ b_doPerfectFinder

int TrkReco::b_doPerfectFinder

Definition at line 153 of file TrkReco.h.

Referenced by clear(), and execute().

◆ b_doSalvage

int TrkReco::b_doSalvage

Definition at line 189 of file TrkReco.h.

Referenced by execute(), and initialize().

◆ b_doT0Determination

int TrkReco::b_doT0Determination

Definition at line 191 of file TrkReco.h.

Referenced by execute().

◆ b_doT0Reset

int TrkReco::b_doT0Reset

Definition at line 147 of file TrkReco.h.

Referenced by execute().

◆ b_fittingFlag

int TrkReco::b_fittingFlag

Definition at line 188 of file TrkReco.h.

Referenced by initialize().

◆ b_fudgeFactor

float TrkReco::b_fudgeFactor

Definition at line 139 of file TrkReco.h.

◆ b_goodTrk

int TrkReco::b_goodTrk

Definition at line 129 of file TrkReco.h.

Referenced by execute().

◆ b_helixFitterChisqMax

float TrkReco::b_helixFitterChisqMax

Definition at line 150 of file TrkReco.h.

Referenced by initialize().

◆ b_mcHit

int TrkReco::b_mcHit

Definition at line 127 of file TrkReco.h.

◆ b_mcPar

int TrkReco::b_mcPar

Definition at line 126 of file TrkReco.h.

◆ b_momentumCut

double TrkReco::b_momentumCut

Definition at line 187 of file TrkReco.h.

Referenced by initialize().

◆ b_nT0ResetMax

int TrkReco::b_nT0ResetMax

Definition at line 148 of file TrkReco.h.

Referenced by execute().

◆ b_nTracksForT0

int TrkReco::b_nTracksForT0

Definition at line 192 of file TrkReco.h.

Referenced by execute().

◆ b_perfectFitting

int TrkReco::b_perfectFitting

Definition at line 154 of file TrkReco.h.

◆ b_RungeKuttaCorrection

int TrkReco::b_RungeKuttaCorrection

Definition at line 135 of file TrkReco.h.

Referenced by execute().

◆ b_salvageLevel

float TrkReco::b_salvageLevel

Definition at line 167 of file TrkReco.h.

Referenced by initialize().

◆ b_sortMode

int TrkReco::b_sortMode

Definition at line 193 of file TrkReco.h.

◆ b_test

int TrkReco::b_test

Definition at line 196 of file TrkReco.h.

Referenced by initialize().

◆ b_timeTest

int TrkReco::b_timeTest

Definition at line 130 of file TrkReco.h.

Referenced by execute(), finalize(), and initialize().

◆ b_tuple

int TrkReco::b_tuple

Definition at line 128 of file TrkReco.h.

Referenced by execute(), finalize(), and initialize().

◆ b_useAllHits

int TrkReco::b_useAllHits

Definition at line 146 of file TrkReco.h.

Referenced by execute().

◆ bad_distance_for_salvage

double TrkReco::bad_distance_for_salvage

Definition at line 202 of file TrkReco.h.

Referenced by initialize().

◆ curl_version

int TrkReco::curl_version

Definition at line 237 of file TrkReco.h.

Referenced by initialize().

◆ determine_one_track

int TrkReco::determine_one_track

Definition at line 214 of file TrkReco.h.

Referenced by initialize().

◆ good_distance_for_salvage

double TrkReco::good_distance_for_salvage

Definition at line 203 of file TrkReco.h.

Referenced by initialize().

◆ m_CalibFlag

int TrkReco::m_CalibFlag

Definition at line 184 of file TrkReco.h.

Referenced by execute().

◆ m_combineTracking

bool TrkReco::m_combineTracking

Definition at line 132 of file TrkReco.h.

Referenced by execute().

◆ m_dropHot

bool TrkReco::m_dropHot

Definition at line 131 of file TrkReco.h.

Referenced by execute().

◆ m_keepBadTdc

bool TrkReco::m_keepBadTdc

Definition at line 133 of file TrkReco.h.

Referenced by execute().

◆ m_keepUnmatch

bool TrkReco::m_keepUnmatch

Definition at line 134 of file TrkReco.h.

Referenced by execute().

◆ mask_distance

double TrkReco::mask_distance

Definition at line 223 of file TrkReco.h.

Referenced by initialize().

◆ merge_exe

int TrkReco::merge_exe

Definition at line 220 of file TrkReco.h.

Referenced by initialize().

◆ merge_ratio

double TrkReco::merge_ratio

Definition at line 221 of file TrkReco.h.

Referenced by initialize().

◆ merge_z_diff

double TrkReco::merge_z_diff

Definition at line 222 of file TrkReco.h.

Referenced by initialize().

◆ min_fullwire

int TrkReco::min_fullwire

Definition at line 205 of file TrkReco.h.

Referenced by initialize().

◆ MIN_RADIUS_OF_STRANGE_TRACK

double TrkReco::MIN_RADIUS_OF_STRANGE_TRACK

Definition at line 244 of file TrkReco.h.

Referenced by initialize().

◆ min_salvage

int TrkReco::min_salvage

Definition at line 201 of file TrkReco.h.

Referenced by initialize().

◆ min_segment

int TrkReco::min_segment

Definition at line 200 of file TrkReco.h.

Referenced by initialize().

◆ min_sequence

int TrkReco::min_sequence

Definition at line 204 of file TrkReco.h.

Referenced by initialize().

◆ min_svd_electrons

double TrkReco::min_svd_electrons

Definition at line 234 of file TrkReco.h.

Referenced by initialize().

◆ minimum_2DTrackLength

double TrkReco::minimum_2DTrackLength

Definition at line 241 of file TrkReco.h.

Referenced by initialize().

◆ minimum_3DTrackLength

double TrkReco::minimum_3DTrackLength

Definition at line 242 of file TrkReco.h.

Referenced by initialize().

◆ minimum_closeHitsLength

double TrkReco::minimum_closeHitsLength

Definition at line 243 of file TrkReco.h.

Referenced by initialize().

◆ minimum_seedLength

double TrkReco::minimum_seedLength

Definition at line 240 of file TrkReco.h.

Referenced by initialize().

◆ on_correction

int TrkReco::on_correction

Definition at line 235 of file TrkReco.h.

Referenced by initialize().

◆ output_2dtracks

int TrkReco::output_2dtracks

Definition at line 236 of file TrkReco.h.

Referenced by initialize().

◆ range_for_axial_last2d_search

double TrkReco::range_for_axial_last2d_search

Definition at line 209 of file TrkReco.h.

Referenced by initialize().

◆ range_for_axial_search

double TrkReco::range_for_axial_search

Definition at line 206 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo1

double TrkReco::range_for_stereo1

Definition at line 225 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo2

double TrkReco::range_for_stereo2

Definition at line 226 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo3

double TrkReco::range_for_stereo3

Definition at line 227 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo4

double TrkReco::range_for_stereo4

Definition at line 228 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo5

double TrkReco::range_for_stereo5

Definition at line 229 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo6

double TrkReco::range_for_stereo6

Definition at line 230 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo_last2d_search

double TrkReco::range_for_stereo_last2d_search

Definition at line 210 of file TrkReco.h.

Referenced by initialize().

◆ range_for_stereo_search

double TrkReco::range_for_stereo_search

Definition at line 207 of file TrkReco.h.

Referenced by initialize().

◆ ratio_used_wire

double TrkReco::ratio_used_wire

Definition at line 224 of file TrkReco.h.

Referenced by initialize().

◆ selector_max_impact

double TrkReco::selector_max_impact

Definition at line 215 of file TrkReco.h.

Referenced by initialize().

◆ selector_max_sigma

double TrkReco::selector_max_sigma

Definition at line 216 of file TrkReco.h.

Referenced by initialize().

◆ selector_replace_dz

double TrkReco::selector_replace_dz

Definition at line 218 of file TrkReco.h.

Referenced by initialize().

◆ selector_strange_pz

double TrkReco::selector_strange_pz

Definition at line 217 of file TrkReco.h.

Referenced by initialize().

◆ stereo_2dfind

int TrkReco::stereo_2dfind

Definition at line 219 of file TrkReco.h.

Referenced by initialize().

◆ superlayer_for_stereo_search

int TrkReco::superlayer_for_stereo_search

Definition at line 208 of file TrkReco.h.

Referenced by initialize().

◆ svd_reconstruction

int TrkReco::svd_reconstruction

Definition at line 233 of file TrkReco.h.

Referenced by initialize().

◆ trace2d_distance

double TrkReco::trace2d_distance

Definition at line 211 of file TrkReco.h.

Referenced by initialize().

◆ trace2d_first_distance

double TrkReco::trace2d_first_distance

Definition at line 212 of file TrkReco.h.

Referenced by initialize().

◆ trace3d_distance

double TrkReco::trace3d_distance

Definition at line 213 of file TrkReco.h.

Referenced by initialize().

◆ ULTIMATE_MIN_RADIUS_OF_STRANGE_TRACK

double TrkReco::ULTIMATE_MIN_RADIUS_OF_STRANGE_TRACK

Definition at line 245 of file TrkReco.h.

Referenced by initialize().

◆ useESTime

int TrkReco::useESTime

Definition at line 142 of file TrkReco.h.

Referenced by execute().

◆ z_cut

double TrkReco::z_cut

Definition at line 231 of file TrkReco.h.

Referenced by initialize().

◆ z_diff_for_last_attend

double TrkReco::z_diff_for_last_attend

Definition at line 232 of file TrkReco.h.

Referenced by initialize().

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

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TrkReco()

Member Function Documentation

◆ beginRun()

◆ clear()

◆ confFinder()

◆ curlFinder()

◆ disp_stat()

◆ dump()

◆ endRun()

◆ execute()

◆ fastClear()

◆ finalize()

◆ initialize()

◆ trackManager()

◆ tracks()

Member Data Documentation

◆ _cdcVersion

◆ b_associateSigma

◆ b_cdcVersion

◆ b_conformalFinder

◆ b_conformalFittingCorrections

◆ b_conformalFittingFlag

◆ b_conformalFraction

◆ b_conformalMaxSigma

◆ b_conformalMinNCores

◆ b_conformalMinNLinksForSegment

◆ b_conformalMinNSegments

◆ b_conformalPerfectSegmentFinding

◆ b_conformalSalvageLoadWidth

◆ b_conformalStereoChisq3

◆ b_conformalStereoChisq4

◆ b_conformalStereoLoadWidth

◆ b_conformalStereoMaxSigma

◆ b_conformalStereoMode

◆ b_conformalStereoSzLinkDistance

◆ b_conformalStereoSzSegmentDistance

◆ b_conformalStereoZ3

◆ b_conformalStereoZ4

◆ b_debugLevel

◆ b_doAssociation

◆ b_doConformalFastFinder

◆ b_doConformalFinder

◆ b_doConformalFinderCosmic

◆ b_doConformalFinderStereo

◆ b_doConformalSlowFinder

◆ b_doCurlFinder

◆ b_doMCAnalysis

◆ b_doMerge

◆ b_doPerfectFinder

◆ b_doSalvage

◆ b_doT0Determination

◆ b_doT0Reset

◆ b_fittingFlag

◆ b_fudgeFactor

◆ b_goodTrk

◆ b_helixFitterChisqMax

◆ b_mcHit

◆ b_mcPar

◆ b_momentumCut

◆ b_nT0ResetMax

◆ b_nTracksForT0

◆ b_perfectFitting

◆ b_RungeKuttaCorrection

◆ b_salvageLevel

◆ b_sortMode

◆ b_test

◆ b_timeTest

◆ b_tuple

◆ b_useAllHits

◆ bad_distance_for_salvage

◆ curl_version

◆ determine_one_track

◆ good_distance_for_salvage

◆ m_CalibFlag

◆ m_combineTracking

◆ m_dropHot