MdcTrackListCsmc Class Reference

#include <MdcTrackListCsmc.h>

Inheritance diagram for MdcTrackListCsmc:

Public Member Functions
	MdcTrackListCsmc (const MdcTrackParams &tkPar)
	~MdcTrackListCsmc ()
int	createFromSegs (MdcSegList *segs, const MdcHitMap *, const MdcDetector *, TrkContext &, double trackT0)
int	finish3d (TrkRecoTrk &trk)
void	remove (MdcTrack *atrack)
Public Member Functions inherited from MdcTrackListBase
	MdcTrackListBase (const MdcTrackParams &tkPar)
virtual	~MdcTrackListBase ()
int	nTrack () const
void	setPlot (int plotFlag)
void	newParams (const MdcTrackParams &tkPar)
void	plot () const
void	store (RecMdcTrackCol *, RecMdcHitCol *)
virtual int	createFromSegs (MdcSegList *, const MdcHitMap *, const MdcDetector *, TrkContext &, double bunchTime)=0
int	arbitrateHits ()
void	dropMultiHotInLayer (const MdcTrack *tk)
void	remove (MdcTrack *atrack)
void	setD0Cut (double d0Cut)
void	setZ0Cut (double z0Cut)
void	setPtCut (double ptCut)

Additional Inherited Members
Static Public Attributes inherited from MdcTrackListBase
static double	m_d0Cut = -999.
static double	m_z0Cut = -999.
static double	m_ptCut = -999.
Protected Attributes inherited from MdcTrackListBase
MdcTrackParams	tkParam

Detailed Description

Definition at line 30 of file MdcTrackListCsmc.h.

Constructor & Destructor Documentation

MdcTrackListCsmc()

MdcTrackListCsmc::MdcTrackListCsmc

(

const MdcTrackParams &

tkPar

)

Definition at line 34 of file MdcTrackListCsmc.cxx.

                                                              : 
  MdcTrackListBase(tkPar) { 
// *************************************************************************
    return;
  }

~MdcTrackListCsmc()

MdcTrackListCsmc::~MdcTrackListCsmc

(

)

Definition at line 41 of file MdcTrackListCsmc.cxx.

{} 

Member Function Documentation

createFromSegs()

int MdcTrackListCsmc::createFromSegs ( MdcSegList *  segs, const MdcHitMap *  , const MdcDetector *  gm, TrkContext &  context, double  trackT0  )

virtual

Implements MdcTrackListBase.

Definition at line 45 of file MdcTrackListCsmc.cxx.

                                                     {
// *************************************************************************
  // Forms tracks out of list of superlayer segments
  int nTracks = 0;
  // Create empty list of stereo segs (to save time)
  MdcSegGrouperSt stereoSegs(gm,tkParam.lPrint);
  // *** Create r-phi track
 
  // Create list of axial segments, treated as on straight tracks
  MdcSegGrouperCsmc straightSegs(gm, tkParam.lPrint);
  straightSegs.fillWithSegs(segs);
  //std::cout<< "after straight fill  " << std::endl; 
  //segs->plot();//yzhang debug
  MdcSeg *seed;
  int goodOnly = 1;
  MdcTrack* trialTrack = 0;
  while (1) {
    seed = segs->getSeed(0,goodOnly);
    if (seed == 0 && goodOnly == 1) {
      segs->resetSeed(gm);
      goodOnly = 0;
      continue;
    }
    else if (seed == 0 && goodOnly == 0) {
      break;
    }
    delete trialTrack;
    trialTrack = 0;            
    int success = straightSegs.combineSegs(trialTrack, seed,
                       context, t0, tkParam.maxSegChisqO);
    if (!success) {
      //std::cout<< " MdcTrackListCsmc::combineSegs not successed!" << std::endl;
      continue;
    }
 
    if (tkParam.lPrint > 1) {
      trialTrack->track().printAll(cout);
    }
    // *** End r-phi track-finding
 
    // *** Add stereo to taste
    stereoSegs.fillWithSegs(segs, trialTrack);
 
    success = stereoSegs.combineSegs(trialTrack, 0, context, t0, 
                     tkParam.maxSegChisqO);
    if (success) {
    // Finish 3-d track;
      success = finish3d(trialTrack->track());
      //      success = finish3d(trialTrack->track()); // GSciolla: try to repeat
    }  
    if (!success) {
      //std::cout<< " MdcTrackListCsmc::finish3d not successed!" << std::endl;
      continue;
    }
 
    if (tkParam.lPrint > 1) {
      trialTrack->track().printAll(cout);
    }
 
    nTracks++;
    append(trialTrack);    // Add to list of Tracks
    trialTrack = 0;
    
  }  // end while(1)   
  delete trialTrack;
  
 
  //cout << " Number of Tracks found= " << nTracks  ;
  //cout << " " << endl;
  return nTracks;
}

finish3d()

int MdcTrackListCsmc::finish3d

(

TrkRecoTrk &

trk

)

Definition at line 120 of file MdcTrackListCsmc.cxx.

                                              {
// ************************************************************************
  int success = 0;
  int nParamFit = 0 ;
 
  TrkErrCode fitResult;
 
  // ------------------------
  //   4 param fit (line)
  // ------------------------
  nParamFit = 4;
  TrkLineMaker makeFit;
  makeFit.changeFit(trk);
  makeFit.setFlipAndDrop(trk, true, true);
  //setCosmic(&trk);    // set hot to cosmics OBSOLETE! REMOVE!
  fitResult = trk.hits()->fit();
  makeFit.setFlipAndDrop(trk, false, false);  
  
  // ---------------------------------------------------------------
  // Are there some TDC that can be replaced with later measurements? ( beginning) 
  // ---------------------------------------------------------------
  // Note to Sciolla: multiple hits commented out for now OK! I will put it back later...
  /*  if( _flags.useMultipleHits ) { 
    int NHITS   = trk.nHit(); 
    
    if (fitResult.success()) {
      int ifirst = 0;     
      for (int ihot = 0; ihot < NHITS ; ihot++) {
    
    MdcHitOnTrack* ahot = (MdcHitOnTrack*)trk.hitAList()[ihot];
    double firstTime = ahot->mdcHit()->rawTime(0); 
    
    // Get the list of times for this Digi 
    const MdcDigi* tmpDigi = ahot->mdcHit()->digi() ; 
//  int nTDChits = tmpDigi->getTdcTimeAList()->length() ;
    int nTDChits = tmpDigi->tmdcits() ;
    
    double newTime=0. ; 
    
    for( int il = 1; il<nTDChits ; il++){ 
      
      // get the time corresponding to this doca ---> distance to time 
      newTime = tmpDigi->TdcTime(il); 
      double newDrift = ahot->layer()->timeToDistance(newTime); 
      
      double tmp_doca  = ahot->dcaToWire();
      double new_diff = fabs(fabs(tmp_doca)-fabs(newDrift)); 
      double tmp_drift = ahot->drift();
      double tmp_drift2 = 
        ahot->layer()->timeToDistance(tmpDigi->TdcTime(0)); 
      double old_diff  = fabs(fabs(tmp_doca)-fabs(tmp_drift2)); 
      
      if( ((old_diff-new_diff)>0.1)  // new time better then old of at least 1 mm  
          && (new_diff<=0.8) ) { // and newresid < 2 mm 
        //  if(_flags.debug) {                  
        if(1) { 
          
          if(ifirst==0) { 
        ifirst = 1 ; 
        cout <<" nHits | time1  |fittim | time2 | doca(mm) | drift1 " 
             <<"  drift1(rec) | drift2 | layer | wire | isActive" 
             << " | chi2 | nactive" << endl; 
          }
          
          int isAct = 1 ; 
          if (!ahot->isActive()) { 
        isAct = 0 ; 
        //      cout<< " >> " ;    
        //  cout << "This hit is not active " << endl ;
          }
          
          cout << "   " << nTDChits 
           << "    " << tmpDigi->TdcTime(0) 
           << "  " << ahot->fitTime()
           << "  " <<  newTime  
           << "  " << tmp_doca*10
           << "  " << tmp_drift*10 
           << "  " << tmp_drift2*10 
           << "  " << newDrift*10 
           << "  " << ahot->layer()->layNum() 
           << "  " << ahot->wire()
           << "  " <<  isAct 
           << "  " << trk.chisq()
           << "  " << trk.nActive() 
           << endl;       
          
          cout << " old/new diff (mm) = " << old_diff*10 
           << " | " << new_diff*10
           << endl ; 
        }     
        
        ahot->setTimeIndex(il) ; 
        //  cout << "new rawtime =" 
        //     << ahot->mdcHit()->rawTime(0)<< endl ; // check the new value 
        
        
        // store results in ntuple to see improvements ... before ...       
        HepTuple* tupleMultHits = _manager->ntuple("multiHits");
        
        tupleMultHits->column("ch2Dof1", trk.chisq()/(trk.nActive() - 4));
        tupleMultHits->column("nActiv1", trk.nActive());
        tupleMultHits->column("resid1",ahot->resid());
        tupleMultHits->column("doca1",ahot->resid()+ahot->drift()); 
        tupleMultHits->column("time1",  ahot->fitTime() );   
        tupleMultHits->column("rawtime1",  ahot->mdcHit()->rawTime(il) ); 
        
        // - make the hit active and usable 
        ahot->setActivity(true);
        
        // now refit the track ... 
        fitResult = trk.fit();
        
        // store results in ntuple to see improvements ... after ...        
        tupleMultHits->column("ch2Dof2", trk.chisq()/(trk.nActive() - 4));
        tupleMultHits->column("nActiv2", trk.nActive());
        tupleMultHits->column("resid2",ahot->resid());
        tupleMultHits->column("doca2",ahot->resid()+ahot->drift()); 
        tupleMultHits->column("time2",  ahot->fitTime() );   
        tupleMultHits->column("rawtime2",  ahot->mdcHit()->rawTime(il) );   
        tupleMultHits->dumpData();
      }
    }
      }
    }  
  } // Are there some TDC that can be replaced with later measurements?  (end)
  */
  
  const TrkFit* tkFit = trk.fitResult();
  double chisqperDOF = 0.;
  if (fitResult.success()) {
    int nDOF = tkFit->nActive() - nParamFit;
    if (nDOF > nParamFit) 
      chisqperDOF = tkFit->chisq() / nDOF;
    else 
      chisqperDOF = tkFit->chisq();
    int goodMatch = 1;
    if (chisqperDOF > tkParam.maxChisq) goodMatch = 0;
    if (tkFit->nActive() < tkParam.minHits) goodMatch = 0;
 
    if (goodMatch) success = 1;
  }
  
  return success;
}

Referenced by createFromSegs().

remove()

void MdcTrackListCsmc::remove

(

MdcTrack *

atrack

)

---

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

• MdcTrackListCsmc.h
• MdcTrackListCsmc.cxx