TCurlFinder.cxx File Reference

#include <fstream>
#include <iostream>
#include <math.h>
#include <stdlib.h>
#include "TrkReco/TCircle.h"
#include "TrkReco/TCurlFinder.h"
#include "TrkReco/TMDCWire.h"
#include "TrkReco/TMDCWireHit.h"
#include "TrkReco/TMLink.h"
#include "TrkReco/TSegmentCurl.h"
#include "TrkReco/TTrack.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 2594 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 179 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 163 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 1520 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;
}