HICSimpleNavigationSchool Class Reference

#include <HICSimpleNavigationSchool.h>

Inheritance diagram for HICSimpleNavigationSchool:

Public Member Functions
int	getExcludedBarrelLayer ()
	HICSimpleNavigationSchool ()
	HICSimpleNavigationSchool (const GeometricSearchTracker *theTracker, const MagneticField *field)
	HICSimpleNavigationSchool (const GeometricSearchTracker *theTracker, const MagneticField *field, int j, int l)
virtual StateType	navigableLayers () const
void	setExcludedBarrelLayer (int &j)
Public Member Functions inherited from NavigationSchool
const std::vector< DetLayer * > &	allLayersInSystem () const
	NavigationSchool ()
virtual	~NavigationSchool ()

Protected Types
typedef std::vector < BarrelDetLayer * >	BDLC
typedef BDLC::iterator	BDLI
typedef std::vector < SimpleBarrelNavigableLayer * >	BNLCType
typedef BDLC::const_iterator	ConstBDLI
typedef FDLC::const_iterator	ConstFDLI
typedef std::vector< const DetLayer * >	DLC
typedef DLC::iterator	DLI
typedef std::vector < ForwardDetLayer * >	FDLC
typedef FDLC::iterator	FDLI
typedef std::vector < SimpleForwardNavigableLayer * >	FNLCType

Protected Member Functions
virtual float	barrelLength ()
virtual void	establishInverseRelations ()
virtual void	linkBarrelLayers (SymmetricLayerFinder &symFinder)
virtual void	linkForwardLayers (SymmetricLayerFinder &symFinder)
virtual void	linkNextBarrelLayer (ForwardDetLayer *fl, BDLC &)
virtual void	linkNextForwardLayer (BarrelDetLayer *, FDLC &)
virtual void	linkNextLargerLayer (BDLI, BDLI, BDLC &)
virtual void	linkNextLayerInGroup (FDLI fli, const FDLC &group, FDLC &reachableFL)
virtual void	linkOuterGroup (ForwardDetLayer *fl, const FDLC &group, FDLC &reachableFL)
virtual void	linkWithinGroup (FDLI fl, const FDLC &group, FDLC &reachableFL)
virtual ConstFDLI	outerRadiusIncrease (FDLI fl, const FDLC &group)
virtual std::vector< FDLC >	splitForwardLayers ()

Protected Attributes
int	excludedBarrelLayer
int	excludedForwardLayer
BDLC	theBarrelLayers
float	theBarrelLength
BNLCType	theBarrelNLC
std::vector< DetLayer * >	theDetLayers
const MagneticField *	theField
FDLC	theForwardLayers
FNLCType	theForwardNLC
FDLC	theLeftLayers
FDLC	theRightLayers
const GeometricSearchTracker *	theTracker
Protected Attributes inherited from NavigationSchool
const std::vector< DetLayer * > *	theAllDetLayersInSystem

Additional Inherited Members
Public Types inherited from NavigationSchool
typedef std::vector < NavigableLayer * >	StateType

Detailed Description

Concrete navigation school for the Tracker

Definition at line 20 of file HICSimpleNavigationSchool.h.

Member Typedef Documentation

typedef std::vector<BarrelDetLayer*> HICSimpleNavigationSchool::BDLC

protected

Definition at line 42 of file HICSimpleNavigationSchool.h.

typedef BDLC::iterator HICSimpleNavigationSchool::BDLI

protected

Definition at line 45 of file HICSimpleNavigationSchool.h.

typedef std::vector< SimpleBarrelNavigableLayer*> HICSimpleNavigationSchool::BNLCType

protected

Definition at line 56 of file HICSimpleNavigationSchool.h.

typedef BDLC::const_iterator HICSimpleNavigationSchool::ConstBDLI

protected

Definition at line 47 of file HICSimpleNavigationSchool.h.

typedef FDLC::const_iterator HICSimpleNavigationSchool::ConstFDLI

protected

Definition at line 48 of file HICSimpleNavigationSchool.h.

typedef std::vector<const DetLayer*> HICSimpleNavigationSchool::DLC

protected

Definition at line 41 of file HICSimpleNavigationSchool.h.

typedef DLC::iterator HICSimpleNavigationSchool::DLI

protected

Definition at line 44 of file HICSimpleNavigationSchool.h.

typedef std::vector<ForwardDetLayer*> HICSimpleNavigationSchool::FDLC

protected

Definition at line 43 of file HICSimpleNavigationSchool.h.

typedef FDLC::iterator HICSimpleNavigationSchool::FDLI

protected

Definition at line 46 of file HICSimpleNavigationSchool.h.

typedef std::vector< SimpleForwardNavigableLayer*> HICSimpleNavigationSchool::FNLCType

protected

Definition at line 57 of file HICSimpleNavigationSchool.h.

Constructor & Destructor Documentation

HICSimpleNavigationSchool::HICSimpleNavigationSchool ( )

inline

Definition at line 23 of file HICSimpleNavigationSchool.h.

: theField(0),theTracker(0){};

HICSimpleNavigationSchool::HICSimpleNavigationSchool	(	const GeometricSearchTracker *	theTracker,
		const MagneticField *	field
	)

Definition at line 29 of file HICSimpleNavigationSchool.cc.

References GeometricSearchTracker::allLayers(), GeometricSearchTracker::barrelLayers(), establishInverseRelations(), GeometricSearchTracker::forwardLayers(), i, linkBarrelLayers(), linkForwardLayers(), NavigationSchool::theAllDetLayersInSystem, theBarrelLayers, theDetLayers, theForwardLayers, theLeftLayers, theRightLayers, and theTracker.

                                                             : 
  theBarrelLength(0),theField(field), theTracker(theInputTracker)
{

  theAllDetLayersInSystem=&theInputTracker->allLayers();
//  vector<DetLayer*> theDetLayers;
  // Get barrel layers
  vector<BarrelDetLayer*> blc = theTracker->barrelLayers(); 
//  int kk=-1;
  for ( vector<BarrelDetLayer*>::iterator i = blc.begin(); i != blc.end(); i++) {
    if((*i)->specificSurface().radius()<20. || (*i)->specificSurface().radius()>65.) {
//     cout<<"Barrel surface "<<(*i)->specificSurface().radius()<<endl;
//     kk++;
//     if(kk==excludedBarrelLayer) continue; 
     theBarrelLayers.push_back( (*i) );
     theDetLayers.push_back( (*i) );
    }
  }

  // get forward layers

  vector<ForwardDetLayer*> flc = theTracker->forwardLayers(); 
  for ( vector<ForwardDetLayer*>::iterator i = flc.begin(); i != flc.end(); i++) {
    if(((*i)->specificSurface().outerRadius()>60.&& (*i)->specificSurface().innerRadius()>20.)|| (*i)->specificSurface().innerRadius()<10.) {
 //     cout<<" Endcap surface "<<(*i)->specificSurface().outerRadius()<<" "<<(*i)->specificSurface().innerRadius()<<endl;
      theForwardLayers.push_back( (*i) );
      theDetLayers.push_back( (*i) );
    }
  }
 
//  cout<<" Number of detectors"<< theDetLayers.size()<<endl;
 
  FDLI middle = find_if( theForwardLayers.begin(), theForwardLayers.end(),
             not1(DetBelowZ(0)));
 // cout<<" Point 0 "<<endl;
  theLeftLayers  = FDLC( theForwardLayers.begin(), middle);
 //  cout<<" Point 1 "<<endl;
  theRightLayers = FDLC( middle, theForwardLayers.end());
 //  cout<<" Point 2 "<<endl; 
  SymmetricLayerFinder symFinder( theForwardLayers);
 //  cout<<" Point 3 "<<endl;
  // only work on positive Z side; negative by mirror symmetry later
  linkBarrelLayers( symFinder);
  //  cout<<" Point 4 "<<endl;
  linkForwardLayers( symFinder);
   // cout<<" Point 5 "<<endl;
  establishInverseRelations();
  // cout<<" Point 6 "<<endl;
}

HICSimpleNavigationSchool::HICSimpleNavigationSchool	(	const GeometricSearchTracker *	theTracker,
		const MagneticField *	field,
		int	j,
		int	l
	)

Definition at line 80 of file HICSimpleNavigationSchool.cc.

References GeometricSearchTracker::allLayers(), GeometricSearchTracker::barrelLayers(), establishInverseRelations(), excludedBarrelLayer, excludedForwardLayer, GeometricSearchTracker::forwardLayers(), i, GetRecoTauVFromDQM_MC_cff::kk, linkBarrelLayers(), linkForwardLayers(), NavigationSchool::theAllDetLayersInSystem, theBarrelLayers, theDetLayers, theForwardLayers, theLeftLayers, theRightLayers, and theTracker.

                                                                             : 
  theBarrelLength(0),theField(field), theTracker(theInputTracker)
{

  excludedBarrelLayer = ll;
  excludedForwardLayer = nn;

  theAllDetLayersInSystem=&theInputTracker->allLayers();
//  vector<DetLayer*> theDetLayers;
  // Get barrel layers
  vector<BarrelDetLayer*> blc = theTracker->barrelLayers(); 
  int kk = 0;

  //cout<<" kk= "<<kk<<" excludedBarrel "<<excludedBarrelLayer<<" excludedForward "<<excludedForwardLayer<<endl;

  for ( vector<BarrelDetLayer*>::iterator i = blc.begin(); i != blc.end(); i++) {
    if((*i)->specificSurface().radius()<20. || (*i)->specificSurface().radius()>65.) {
//     cout<<"Barrel surface "<<kk<<" "<<(*i)->specificSurface().radius()<<" "<<excludedBarrelLayer<<endl;
     if(kk != excludedBarrelLayer ) {
      //cout<<"Included layer "<<endl;  
      theBarrelLayers.push_back( (*i) );
      theDetLayers.push_back( (*i) );
     }
       else {
  //      cout<<"Excluded layer "<<endl;
       }
    }
    kk++; 
  }

  // get forward layers
  if( excludedForwardLayer > -1 ) {
  kk=0; 
  vector<ForwardDetLayer*> flc = theTracker->forwardLayers(); 
  for ( vector<ForwardDetLayer*>::iterator i = flc.begin(); i != flc.end(); i++) 
  {
    if(((*i)->specificSurface().outerRadius()>60.&& 
                  (*i)->specificSurface().innerRadius()>20.)|| (*i)->specificSurface().innerRadius()<10.) 
    {
//      cout<<" Endcap surface "<<kk<<" "<<(*i)->specificSurface().position().z()<<endl;
       if(kk != excludedForwardLayer && kk != excludedForwardLayer + 14) 
       {
        //cout<<"Excluded layer "<<endl;
        theForwardLayers.push_back( (*i) );
        theDetLayers.push_back( (*i) );
       } else {
  //       cout<<"Excluded layer "<<endl; 
       }
     } // radius
      kk++;
   } // forward layers
   } else {
  vector<ForwardDetLayer*> flc = theTracker->forwardLayers(); 
  for ( vector<ForwardDetLayer*>::iterator i = flc.begin(); i != flc.end(); i++) 
  {
    if(((*i)->specificSurface().outerRadius()>60.&& 
                  (*i)->specificSurface().innerRadius()>20.)|| (*i)->specificSurface().innerRadius()<10.) 
    {
//      cout<<" Endcap surface "<<kk<<" "<<(*i)->specificSurface().position().z()<<endl;
        //cout<<"Excluded layer "<<endl;
        theForwardLayers.push_back( (*i) );
        theDetLayers.push_back( (*i) );
     } // radius
      kk++;
   } // forward layers
   } // No Excluded layer
 
//  cout<<" Number of detectors"<< theDetLayers.size()<<endl;
 
  FDLI middle = find_if( theForwardLayers.begin(), theForwardLayers.end(),
             not1(DetBelowZ(0)));
  //cout<<" Point 0 "<<endl;
  theLeftLayers  = FDLC( theForwardLayers.begin(), middle);
  // cout<<" Point 1 "<<endl;
  theRightLayers = FDLC( middle, theForwardLayers.end());
  // cout<<" Point 2 "<<endl; 
  SymmetricLayerFinder symFinder( theForwardLayers);
  // cout<<" Point 3 "<<endl;
  // only work on positive Z side; negative by mirror symmetry later
  linkBarrelLayers( symFinder);
  //  cout<<" Point 4 "<<endl;
  linkForwardLayers( symFinder);
  //  cout<<" Point 5 "<<endl;
  establishInverseRelations();
  // cout<<" Point 6 "<<endl;
}

Member Function Documentation

float HICSimpleNavigationSchool::barrelLength ( )

protectedvirtual

Definition at line 450 of file HICSimpleNavigationSchool.cc.

References i, LogDebug, max(), theBarrelLayers, and theBarrelLength.

Referenced by linkNextBarrelLayer().

{
  if ( theBarrelLength < 1.) {
    for (BDLI i=theBarrelLayers.begin(); i!=theBarrelLayers.end(); i++) {
      theBarrelLength = max( theBarrelLength,
                 (**i).surface().bounds().length() / 2.f);
    }

    LogDebug("TkNavigation") << "The barrel length is " << theBarrelLength ;
  }
  return theBarrelLength;
}

void HICSimpleNavigationSchool::establishInverseRelations ( )

protectedvirtual

Definition at line 463 of file HICSimpleNavigationSchool.cc.

References i, insideOut, SimpleNavigableLayer::setInwardLinks(), theBarrelLayers, theDetLayers, and theForwardLayers.

Referenced by HICSimpleNavigationSchool().

                                                          {

  NavigationSetter setter(*this);

    // find for each layer which are the barrel and forward
    // layers that point to it
    typedef map<const DetLayer*, vector<BarrelDetLayer*>, less<const DetLayer*> > BarrelMapType;
    typedef map<const DetLayer*, vector<ForwardDetLayer*>, less<const DetLayer*> > ForwardMapType;

   // cout<<"SimpleNavigationSchool::establishInverseRelations::point 0 "<<endl;
    BarrelMapType reachedBarrelLayersMap;
    ForwardMapType reachedForwardLayersMap;


    for ( BDLI bli = theBarrelLayers.begin();
        bli!=theBarrelLayers.end(); bli++) {
      DLC reachedLC = (**bli).nextLayers( insideOut);
      for ( DLI i = reachedLC.begin(); i != reachedLC.end(); i++) {
        reachedBarrelLayersMap[*i].push_back( *bli);
      }
    }
  //  cout<<"SimpleNavigationSchool::establishInverseRelations::point 1 "<<endl;
    for ( FDLI fli = theForwardLayers.begin();
        fli!=theForwardLayers.end(); fli++) {
      DLC reachedLC = (**fli).nextLayers( insideOut);
      for ( DLI i = reachedLC.begin(); i != reachedLC.end(); i++) {
        reachedForwardLayersMap[*i].push_back( *fli);
      }
    }

  //  cout<<"SimpleNavigationSchool::establishInverseRelations::point 2:: size od detlayers "<<theDetLayers.size()<<endl;
//    vector<DetLayer*> lc = theTracker->allLayers();
    for ( vector<DetLayer*>::iterator i = theDetLayers.begin(); i != theDetLayers.end(); i++) {
      SimpleNavigableLayer* navigableLayer =
    dynamic_cast<SimpleNavigableLayer*>((**i).navigableLayer());
      navigableLayer->setInwardLinks( reachedBarrelLayersMap[*i],reachedForwardLayersMap[*i] );
    }
  // cout<<"SimpleNavigationSchool::establishInverseRelations::point 3 "<<endl; 
}

int HICSimpleNavigationSchool::getExcludedBarrelLayer ( )

inline

Definition at line 31 of file HICSimpleNavigationSchool.h.

References excludedBarrelLayer.

{return excludedBarrelLayer;};

void HICSimpleNavigationSchool::linkBarrelLayers ( SymmetricLayerFinder &  symFinder )

protectedvirtual

Definition at line 184 of file HICSimpleNavigationSchool.cc.

References i, linkNextForwardLayer(), linkNextLargerLayer(), SymmetricLayerFinder::mirror(), theBarrelLayers, theBarrelNLC, and theField.

Referenced by HICSimpleNavigationSchool().

{
  // Link barrel layers outwards
  for ( BDLI i = theBarrelLayers.begin(); i != theBarrelLayers.end(); i++) {
    BDLC reachableBL;
    FDLC leftFL;
    FDLC rightFL;

    // always add next barrel layer first
    if ( i+1 != theBarrelLayers.end()) reachableBL.push_back(*(i+1));
 
    // Add closest reachable forward layer (except for last BarrelLayer)
    if (i != theBarrelLayers.end() - 1) {
      linkNextForwardLayer( *i, rightFL);
    }

    // Add next BarrelLayer with length larger than the current BL
    if ( i+2 < theBarrelLayers.end()) {
      linkNextLargerLayer( i, theBarrelLayers.end(), reachableBL);
    }

    theBarrelNLC.push_back( new 
       SimpleBarrelNavigableLayer( *i, reachableBL,
                   symFinder.mirror(rightFL),
                   rightFL,theField, 5.));
  }
}

void HICSimpleNavigationSchool::linkForwardLayers ( SymmetricLayerFinder &  symFinder )

protectedvirtual

Definition at line 249 of file HICSimpleNavigationSchool.cc.

References g, i, linkNextBarrelLayer(), linkNextLayerInGroup(), linkOuterGroup(), linkWithinGroup(), LogDebug, SymmetricLayerFinder::mirror(), splitForwardLayers(), theField, and theForwardNLC.

Referenced by HICSimpleNavigationSchool().

{

  // handle right side first, groups are only on the right 
  vector<FDLC> groups = splitForwardLayers();

  LogDebug("TkNavigation") << "SimpleNavigationSchool, Forward groups size = " << groups.size() ;
  for (vector<FDLC>::iterator g = groups.begin(); g != groups.end(); g++) {
    LogDebug("TkNavigation") << "group " << g - groups.begin() << " has " 
                 << g->size() << " layers " ;
  }

  for ( vector<FDLC>::iterator group = groups.begin();
    group != groups.end(); group++) {

    for ( FDLI i = group->begin(); i != group->end(); i++) {

      BDLC reachableBL;
      FDLC reachableFL;
 
      // Always connect to next barrel layer first, if exists
      linkNextBarrelLayer( *i, reachableBL);

      // Then always connect to next forward layer of "same" size, 
      // and layers of larger inner Radius
      linkNextLayerInGroup( i, *group, reachableFL);

      // Then connect to next N fw layers of next size
      if ( group+1 != groups.end()) {
    linkOuterGroup( *i, *(group+1), reachableFL);
      }

      // or connect within the group if outer radius increases
      linkWithinGroup( i, *group, reachableFL);

      theForwardNLC.push_back( new SimpleForwardNavigableLayer( *i,reachableBL,
                                reachableFL,
                                theField,
                                5.));
      theForwardNLC.push_back( new SimpleForwardNavigableLayer( symFinder.mirror(*i),
                                reachableBL,
                                symFinder.mirror(reachableFL),
                                theField,
                                5.));

    }
  }

//    // now the left side by symmetry
//    for ( FDLI ileft = theLeftLayers.begin(); 
//      ileft != theLeftLayers.end(); ileft++) {
//      ForwardDetLayer* right = symFinder.mirror( *ileft);
    
//      theForwardNLC.push_back( new 
//         SimpleForwardNavigableLayer( *ileft , right->nextBarrelLayers(),
//                            symFinder.mirror(right->nextForwardLayers())));
//    }
}

void HICSimpleNavigationSchool::linkNextBarrelLayer ( ForwardDetLayer *  fl, BDLC &  reachableBL  )

protectedvirtual

Definition at line 308 of file HICSimpleNavigationSchool.cc.

References barrelLength(), GeometricSearchDet::position(), theBarrelLayers, and PV3DBase< T, PVType, FrameType >::z().

Referenced by linkForwardLayers().

{
  if ( fl->position().z() > barrelLength()) return;

 // float outerRadius = fl->specificSurface().outerRadius();
  float zpos        = fl->position().z();
  for ( BDLI bli = theBarrelLayers.begin(); bli != theBarrelLayers.end(); bli++) {
//    if ( outerRadius < (**bli).specificSurface().radius() &&
//   zpos        < (**bli).surface().bounds().length() / 2.)
   if( fabs(zpos) > 130. && fabs(zpos) < 132. ) 
    {
      //cout<<" Forward layer "<<fl->position().z()<<" to Barrel layer "<<(**bli).specificSurface().radius()<<endl;
      reachableBL.push_back( *bli);
      return;
    }
  }
}

void HICSimpleNavigationSchool::linkNextForwardLayer ( BarrelDetLayer *  bl, FDLC &  rightFL  )

protectedvirtual

Definition at line 212 of file HICSimpleNavigationSchool.cc.

References theRightLayers.

Referenced by linkBarrelLayers().

{
  // find first forward layer with larger Z and larger outer radius
 // float length = bl->surface().bounds().length() / 2.;
//  float radius = bl->specificSurface().radius();
  for ( FDLI fli = theRightLayers.begin();
    fli != theRightLayers.end(); fli++) {
   // if ( length < (**fli).position().z() &&
   //    radius < (**fli).specificSurface().outerRadius()) 
  if(fabs((**fli).position().z())>130. && fabs((**fli).position().z())<132.)
    {
#ifdef DEBUG
//      cout<<" Add to barrel layer "<<radius<<" Forward layer "<<(**fli).position().z()<<endl;
#endif
      rightFL.push_back( *fli);
      return;
    }
  }
}

void HICSimpleNavigationSchool::linkNextLargerLayer ( BDLI  bli, BDLI  end, BDLC &  reachableBL  )

protectedvirtual

Definition at line 233 of file HICSimpleNavigationSchool.cc.

References end, epsilon, and i.

Referenced by linkBarrelLayers().

{
  // compare length of next layer with length of following ones
  float length = (**(bli+1)).surface().bounds().length();
  float epsilon = 0.1;

  for ( BDLI i = bli+2; i < end; i++) {
    if ( length + epsilon < (**i).surface().bounds().length()) {
      reachableBL.push_back( *i);
      return;
    }
  }
}

void HICSimpleNavigationSchool::linkNextLayerInGroup ( FDLI  fli, const FDLC &  group, FDLC &  reachableFL  )

protectedvirtual

Definition at line 328 of file HICSimpleNavigationSchool.cc.

References epsilon, and i.

Referenced by linkForwardLayers().

{
  // Always connect to next forward layer of "same" size, if exists
  if ( fli+1 != group.end()) {
    reachableFL.push_back( *(fli+1));
    // If that layer has an inner radius larger then the current one
    // also connect ALL next disks of same radius.
    float innerRThis = (**fli).specificSurface().innerRadius();
    float innerRNext =  (**(fli+1)).specificSurface().innerRadius();
    const float epsilon = 2.f;

    if (innerRNext > innerRThis + epsilon) {
      // next disk is smaller, so it doesn't cover fully subsequent ones
      // of same radius

      int i = 2;
      while ( (fli+i) != group.end()) {
    if ( (**(fli+i)).specificSurface().innerRadius() < 
         innerRNext + epsilon) {
      // following disk has not increased in ineer radius 
      reachableFL.push_back( *(fli+i));
      i++;
    } else {
      break;
    }
      }
    }
  }
}

void HICSimpleNavigationSchool::linkOuterGroup ( ForwardDetLayer *  fl, const FDLC &  group, FDLC &  reachableFL  )

protectedvirtual

Definition at line 361 of file HICSimpleNavigationSchool.cc.

References first, prof2calltree::last, min, GeometricSearchDet::position(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by linkForwardLayers().

{

  // insert N layers with Z grater than fl

  ConstFDLI first = find_if( group.begin(), group.end(), 
                 not1( DetBelowZ( fl->position().z())));
  if ( first != group.end()) {

    // Hard-wired constant!!!!!!
    ConstFDLI last = min( first + 7, group.end());

    reachableFL.insert( reachableFL.end(), first, last);
  }
}

void HICSimpleNavigationSchool::linkWithinGroup ( FDLI  fl, const FDLC &  group, FDLC &  reachableFL  )

protectedvirtual

Definition at line 379 of file HICSimpleNavigationSchool.cc.

References outerRadiusIncrease().

Referenced by linkForwardLayers().

{
  ConstFDLI biggerLayer = outerRadiusIncrease( fl, group);
  if ( biggerLayer != group.end() && biggerLayer != fl+1) {
    reachableFL.push_back( *biggerLayer);
  }
}

HICSimpleNavigationSchool::StateType HICSimpleNavigationSchool::navigableLayers ( ) const

virtual

Implements NavigationSchool.

Definition at line 169 of file HICSimpleNavigationSchool.cc.

References query::result, theBarrelNLC, and theForwardNLC.

{
  StateType result;
  for ( vector< SimpleBarrelNavigableLayer*>::const_iterator 
      ib = theBarrelNLC.begin(); ib != theBarrelNLC.end(); ib++) {
    result.push_back( *ib);
  }
  for ( vector< SimpleForwardNavigableLayer*>::const_iterator 
      ifl = theForwardNLC.begin(); ifl != theForwardNLC.end(); ifl++) {
    result.push_back( *ifl);
  }
  return result;
}

HICSimpleNavigationSchool::ConstFDLI HICSimpleNavigationSchool::outerRadiusIncrease ( FDLI  fl, const FDLC &  group  )

protectedvirtual

Definition at line 390 of file HICSimpleNavigationSchool.cc.

References epsilon.

Referenced by linkWithinGroup().

{
  const float epsilon = 5.f;
  float outerRadius = (**fl).specificSurface().outerRadius();
  while ( ++fl != group.end()) {
    if ( (**fl).specificSurface().outerRadius() > outerRadius + epsilon) {
      return fl;
    }
  }
  return fl;
}

void HICSimpleNavigationSchool::setExcludedBarrelLayer ( int &  j )

inline

Definition at line 29 of file HICSimpleNavigationSchool.h.

References excludedBarrelLayer, and j.

{excludedBarrelLayer=j;}; 

vector< HICSimpleNavigationSchool::FDLC > HICSimpleNavigationSchool::splitForwardLayers ( )

protectedvirtual

Definition at line 403 of file HICSimpleNavigationSchool.cc.

References begin, cond::rpcobimon::current, end, i, LogDebug, query::result, python.multivaluedict::sort(), and theRightLayers.

Referenced by linkForwardLayers().

{
  // only work on positive Z side; negative by mirror symmetry later

  FDLC myRightLayers( theRightLayers);
  FDLI begin = myRightLayers.begin();
  FDLI end   = myRightLayers.end();

  // sort according to inner radius
  sort ( begin, end, DiskLessInnerRadius()); 

  // partition in cylinders
  vector<FDLC> result;
  FDLC current;
  current.push_back( *begin);
  for ( FDLI i = begin+1; i != end; i++) {

    LogDebug("TkNavigation") << "(**i).specificSurface().innerRadius()      = "
                 << (**i).specificSurface().innerRadius() << endl
                 << "(**(i-1)).specificSurface().outerRadius()) = "
                 << (**(i-1)).specificSurface().outerRadius() ;

    // if inner radius of i is larger than outer radius of i-1 then split!
    if ( (**i).specificSurface().innerRadius() > 
     (**(i-1)).specificSurface().outerRadius()) {

      LogDebug("TkNavigation") << "found break between groups" ;

      // sort layers in group along Z
      sort ( current.begin(), current.end(), DetLessZ());

      result.push_back(current);
      current.clear();
    }
    current.push_back(*i);
  }
  result.push_back(current); // save last one too 

  // now sort subsets in Z
  for ( vector<FDLC>::iterator ivec = result.begin();
    ivec != result.end(); ivec++) {
    sort( ivec->begin(), ivec->end(), DetLessZ());
  }

  return result;
}

Member Data Documentation

int HICSimpleNavigationSchool::excludedBarrelLayer

protected

Definition at line 39 of file HICSimpleNavigationSchool.h.

Referenced by getExcludedBarrelLayer(), HICSimpleNavigationSchool(), and setExcludedBarrelLayer().

int HICSimpleNavigationSchool::excludedForwardLayer

protected

Definition at line 40 of file HICSimpleNavigationSchool.h.

Referenced by HICSimpleNavigationSchool().

BDLC HICSimpleNavigationSchool::theBarrelLayers

protected

Definition at line 50 of file HICSimpleNavigationSchool.h.

Referenced by barrelLength(), establishInverseRelations(), HICSimpleNavigationSchool(), linkBarrelLayers(), and linkNextBarrelLayer().

float HICSimpleNavigationSchool::theBarrelLength

protected

Definition at line 54 of file HICSimpleNavigationSchool.h.

Referenced by barrelLength().

BNLCType HICSimpleNavigationSchool::theBarrelNLC

protected

Definition at line 58 of file HICSimpleNavigationSchool.h.

Referenced by linkBarrelLayers(), and navigableLayers().

std::vector<DetLayer*> HICSimpleNavigationSchool::theDetLayers

protected

Definition at line 86 of file HICSimpleNavigationSchool.h.

Referenced by establishInverseRelations(), and HICSimpleNavigationSchool().

const MagneticField* HICSimpleNavigationSchool::theField

protected

Definition at line 87 of file HICSimpleNavigationSchool.h.

Referenced by linkBarrelLayers(), and linkForwardLayers().

FDLC HICSimpleNavigationSchool::theForwardLayers

protected

Definition at line 51 of file HICSimpleNavigationSchool.h.

Referenced by establishInverseRelations(), and HICSimpleNavigationSchool().

FNLCType HICSimpleNavigationSchool::theForwardNLC

protected

Definition at line 59 of file HICSimpleNavigationSchool.h.

Referenced by linkForwardLayers(), and navigableLayers().

FDLC HICSimpleNavigationSchool::theLeftLayers

protected

Definition at line 53 of file HICSimpleNavigationSchool.h.

Referenced by HICSimpleNavigationSchool().

FDLC HICSimpleNavigationSchool::theRightLayers

protected

Definition at line 52 of file HICSimpleNavigationSchool.h.

Referenced by HICSimpleNavigationSchool(), linkNextForwardLayer(), and splitForwardLayers().

const GeometricSearchTracker* HICSimpleNavigationSchool::theTracker

protected

Definition at line 88 of file HICSimpleNavigationSchool.h.

Referenced by HICSimpleNavigationSchool().