/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/TrkReco/TrkReco-00-08-59-patch4-slc6tag/src/TCurlFinder.cxx File Reference

#include <iostream>
#include <fstream>
#include <math.h>
#include <stdlib.h>
#include "CLHEP/String/Strings.h"
#include "TrkReco/TCurlFinder.h"
#include "TrkReco/TMDCWire.h"
#include "TrkReco/TMDCWireHit.h"
#include "TrkReco/TMLink.h"
#include "TrkReco/TCircle.h"
#include "TrkReco/TTrack.h"
#include "TrkReco/TSegmentCurl.h"
#include "MdcTables/MdcTables.h"

Go to the source code of this file.

Defines
#define	DEBUG_CURL_DUMP   0
#define	DEBUG_CURL_SEGMENT   0
#define	DEBUG_CURL_GNUPLOT   0
#define	DEBUG_CURL_MC   0
Functions
int	sortBySequentialLength (const void *av, const void *bv)
int	sortByArcLength (const void *av, const void *bv)
int	TCurlFinder_doubleCompare (const void *i, const void *j)
TMLink *	findIsolatedCloseHits (TMLink *link)

---

Define Documentation

#define DEBUG_CURL_DUMP   0

Definition at line 26 of file TCurlFinder.cxx.

#define DEBUG_CURL_GNUPLOT   0

Definition at line 28 of file TCurlFinder.cxx.

#define DEBUG_CURL_MC   0

Definition at line 29 of file TCurlFinder.cxx.

#define DEBUG_CURL_SEGMENT   0

Definition at line 27 of file TCurlFinder.cxx.

---

Function Documentation

TMLink * findIsolatedCloseHits

(

TMLink *

link

)

Definition at line 2400 of file TCurlFinder.cxx.

References genRecEmupikp::i, ganga-rec::j, TMDCWire::layerId(), TMDCWire::localId(), TMLink::neighbor(), and TMLink::wire().

Referenced by TCurlFinder::findCloseHits(), TCurlFinder::searchAxialCand(), and TCurlFinder::searchHits().

{
  int nNeighbor = 0;
  int nIsolated = 0;
  TMLink * isolatedLink[2] = {NULL,NULL};
  unsigned layerID = link->wire()->layerId();
  unsigned localID = link->wire()->localId();
  for(int i=0;i<6;++i){
    if(link->neighbor(i)){
      if(link->neighbor(i)->wire()->layerId() == layerID){
        ++nNeighbor;
        int isolated = 1;
        int testEach = 0; // for test
        for(int j=0;j<6;++j){
          if(link->neighbor(i)->neighbor(j)){
            if(link->neighbor(i)->neighbor(j)->wire()->layerId() == layerID &&
               link->neighbor(i)->neighbor(j)->wire()->localId() != localID){
              isolated = 0;
            }
#if 1
            else if(link->neighbor(i)->neighbor(j)->wire()->layerId() == layerID &&
                    link->neighbor(i)->neighbor(j)->wire()->localId() == localID){
              testEach = 1;
            }
#endif
          }else break;
        }
        if(isolated == 1){
          if(nIsolated < 2)
            isolatedLink[nIsolated] = link->neighbor(i);
          ++nIsolated;    
        }
#if 1
        if(testEach == 0){
          std::cout << "Why?? Neighborhood info. dose not exist!!" << std::endl;
        }
#endif
      }
    }else break;
  }
#if 0
  std::cout << "isolated/neighbor # = " << nIsolated << "/" << nNeighbor << std::endl;
  std::cout << "layer ID = " << layerID << " ";
  std::cout << "local ID = " << localID << " --> ";
  if(isolatedLink[0])std::cout << isolatedLink[0]->wire()->localId() << " ";
  if(isolatedLink[1])std::cout << isolatedLink[1]->wire()->localId() << " ";
  std::cout << std::endl;
#endif
  if(nIsolated == 1 &&
     nNeighbor == 1 && isolatedLink[0])return isolatedLink[0];
  else return NULL;
}

int sortByArcLength	(	const void *	av,
		const void *	bv
	)

Definition at line 226 of file TCurlFinder.cxx.

                                                {
  const TMLink **a((const TMLink**)av);
  const TMLink **b((const TMLink**)bv);
  if( (*a)->position().x() > (*b)->position().x() ){
    return 1;
  }else if( (*a)->position().x() == (*b)->position().x() ){
    return 0;
  }else{
    return -1;
  }
}

int sortBySequentialLength	(	const void *	av,
		const void *	bv
	)

Definition at line 200 of file TCurlFinder.cxx.

Referenced by TCurlFinder::checkSortSegments().

                                                         {
  const TSegmentCurl **a((const TSegmentCurl**)av);
  const TSegmentCurl **b((const TSegmentCurl**)bv);
  if( (*a)->maxSeq() < (*b)->maxSeq() ){
    return 1;
  }else if( (*a)->maxSeq() == (*b)->maxSeq() ){
    return 0;
  }else{
    return -1;
  }
}

int TCurlFinder_doubleCompare	(	const void *	i,
		const void *	j
	)

Definition at line 1454 of file TCurlFinder.cxx.

Referenced by TCurlFinder::trace2DTrack(), and TCurlFinder::trace3DTrack().

                                                       {
  if(*(static_cast<const double*>(i)) > *(static_cast<const double*>(j)))return 1;
  if(*(static_cast<const double*>(i)) < *(static_cast<const double*>(j)))return -1;
  return 0;
}

---