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.

Macros
#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)

Macro Definition Documentation

DEBUG_CURL_DUMP

#define DEBUG_CURL_DUMP   0

Definition at line 26 of file TCurlFinder.cxx.

DEBUG_CURL_GNUPLOT

#define DEBUG_CURL_GNUPLOT   0

Definition at line 28 of file TCurlFinder.cxx.

DEBUG_CURL_MC

#define DEBUG_CURL_MC   0

Definition at line 29 of file TCurlFinder.cxx.

DEBUG_CURL_SEGMENT

#define DEBUG_CURL_SEGMENT   0

Definition at line 27 of file TCurlFinder.cxx.

Function Documentation

findIsolatedCloseHits()

TMLink * findIsolatedCloseHits

(

TMLink *

link

)

Definition at line 2400 of file TCurlFinder.cxx.

{
  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;
}

sortByArcLength()

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;
  }
}

sortBySequentialLength()

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

Definition at line 200 of file TCurlFinder.cxx.

                                                         {
  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;
  }
}

TCurlFinder_doubleCompare()

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

Definition at line 1454 of file TCurlFinder.cxx.

                                                       {
  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;
}