SimpleNavigationSchool Class Reference

#include <SimpleNavigationSchool.h>

Inheritance diagram for SimpleNavigationSchool:

Public Member Functions
StateType	navigableLayers () override
	SimpleNavigationSchool ()
	SimpleNavigationSchool (const GeometricSearchTracker *tracker, const MagneticField *field)
	~SimpleNavigationSchool () override
Public Member Functions inherited from TkNavigationSchool
const MagneticField &	field () const
const GeometricSearchTracker &	searchTracker () const
	TkNavigationSchool (const GeometricSearchTracker *tracker, const MagneticField *field)
Public Member Functions inherited from NavigationSchool
const std::vector< const DetLayer * > &	allLayersInSystem () const
template<typename... Args>
std::vector< const DetLayer * >	compatibleLayers (const DetLayer &detLayer, Args &&...args) const
	Returns all layers compatible. More...
	NavigationSchool ()
template<typename... Args>
std::vector< const DetLayer * >	nextLayers (const DetLayer &detLayer, Args &&...args) const
virtual	~NavigationSchool ()

Protected Types
typedef std::vector< const 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< const ForwardDetLayer * >	FDLC
typedef FDLC::iterator	FDLI
typedef std::vector < SimpleForwardNavigableLayer * >	FNLCType

Protected Member Functions
virtual float	barrelLength ()
void	cleanMemory ()
virtual void	establishInverseRelations ()
virtual void	linkBarrelLayers (SymmetricLayerFinder &symFinder)
virtual void	linkForwardLayers (SymmetricLayerFinder &symFinder)
virtual void	linkNextBarrelLayer (ForwardDetLayer const *fl, BDLC &)
virtual void	linkNextForwardLayer (BarrelDetLayer const *, FDLC &)
virtual void	linkNextLargerLayer (BDLI, BDLI, BDLC &)
virtual void	linkNextLayerInGroup (FDLI fli, const FDLC &group, FDLC &reachableFL)
virtual void	linkOuterGroup (ForwardDetLayer const *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 Member Functions inherited from NavigationSchool
void	setState (const StateType &state)

Protected Attributes
BDLC	theBarrelLayers
float	theBarrelLength = 0
BNLCType	theBarrelNLC
FDLC	theForwardLayers
FNLCType	theForwardNLC
FDLC	theLeftLayers
FDLC	theRightLayers
Protected Attributes inherited from TkNavigationSchool
const MagneticField *	theField
const GeometricSearchTracker *	theTracker
Protected Attributes inherited from NavigationSchool
const std::vector< const DetLayer * > *	theAllDetLayersInSystem
StateType	theAllNavigableLayer

Private Member Functions
void	init ()

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

Detailed Description

Concrete navigation school for the Tracker

Definition at line 19 of file SimpleNavigationSchool.h.

Member Typedef Documentation

typedef std::vector<const BarrelDetLayer*> SimpleNavigationSchool::BDLC

protected

Definition at line 33 of file SimpleNavigationSchool.h.

typedef BDLC::iterator SimpleNavigationSchool::BDLI

protected

Definition at line 36 of file SimpleNavigationSchool.h.

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

protected

Definition at line 47 of file SimpleNavigationSchool.h.

typedef BDLC::const_iterator SimpleNavigationSchool::ConstBDLI

protected

Definition at line 38 of file SimpleNavigationSchool.h.

typedef FDLC::const_iterator SimpleNavigationSchool::ConstFDLI

protected

Definition at line 39 of file SimpleNavigationSchool.h.

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

protected

Definition at line 32 of file SimpleNavigationSchool.h.

typedef DLC::iterator SimpleNavigationSchool::DLI

protected

Definition at line 35 of file SimpleNavigationSchool.h.

typedef std::vector<const ForwardDetLayer*> SimpleNavigationSchool::FDLC

protected

Definition at line 34 of file SimpleNavigationSchool.h.

typedef FDLC::iterator SimpleNavigationSchool::FDLI

protected

Definition at line 37 of file SimpleNavigationSchool.h.

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

protected

Definition at line 48 of file SimpleNavigationSchool.h.

Constructor & Destructor Documentation

SimpleNavigationSchool::SimpleNavigationSchool ( )

inline

Definition at line 21 of file SimpleNavigationSchool.h.

: TkNavigationSchool(nullptr, nullptr) {}

SimpleNavigationSchool::SimpleNavigationSchool ( const GeometricSearchTracker *  tracker, const MagneticField *  field  )

inline

Definition at line 22 of file SimpleNavigationSchool.h.

References init.

      : TkNavigationSchool(tracker, field) {
    init();
  }

SimpleNavigationSchool::~SimpleNavigationSchool ( )

inlineoverride

Definition at line 26 of file SimpleNavigationSchool.h.

{ cleanMemory(); }

Member Function Documentation

float SimpleNavigationSchool::barrelLength ( )

protectedvirtual

Definition at line 313 of file SimpleNavigationSchool.cc.

References mps_fire::i, LogDebug, and SiStripPI::max.

                                           {
  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 SimpleNavigationSchool::cleanMemory ( )

protected

Definition at line 53 of file SimpleNavigationSchool.cc.

References cuy::ib.

Referenced by BeamHaloNavigationSchool::~BeamHaloNavigationSchool(), and CosmicNavigationSchool::~CosmicNavigationSchool().

                                         {
  // free the memory allocated to the SimpleNavigableLayers
  for (vector<SimpleBarrelNavigableLayer*>::const_iterator ib = theBarrelNLC.begin(); ib != theBarrelNLC.end(); ib++) {
    delete (*ib);
  }
  theBarrelNLC.clear();
  for (vector<SimpleForwardNavigableLayer*>::const_iterator ifl = theForwardNLC.begin(); ifl != theForwardNLC.end();
       ifl++) {
    delete (*ifl);
  }
  theForwardNLC.clear();
}

void SimpleNavigationSchool::establishInverseRelations ( )

protectedvirtual

Reimplemented in BeamHaloNavigationSchool.

Definition at line 324 of file SimpleNavigationSchool.cc.

References NavigableLayer::detLayer(), mps_fire::i, and insideOut.

                                                       {
  // NavigationSetter setter(*this);

  setState(navigableLayers());

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

  BarrelMapType reachedBarrelLayersMap;
  ForwardMapType reachedForwardLayersMap;

  for (auto bli : theBarrelLayers) {
    auto reachedLC = nextLayers(*bli, insideOut);
    for (auto i : reachedLC) {
      reachedBarrelLayersMap[i].push_back(bli);
    }
  }

  for (auto fli : theForwardLayers) {
    auto reachedLC = nextLayers(*fli, insideOut);
    for (auto i : reachedLC) {
      reachedForwardLayersMap[i].push_back(fli);
    }
  }

  /*
    vector<DetLayer*> lc = theTracker->allLayers();
    for ( vector<DetLayer*>::iterator i = lc.begin(); i != lc.end(); i++) {
      SimpleNavigableLayer* navigableLayer =
    dynamic_cast<SimpleNavigableLayer*>((**i).navigableLayer());
      navigableLayer->setInwardLinks( reachedBarrelLayersMap[*i],reachedForwardLayersMap[*i] );
    }
    */

  for (auto nl : theAllNavigableLayer) {
    if (!nl)
      continue;
    auto navigableLayer = static_cast<SimpleNavigableLayer*>(nl);
    auto dl = nl->detLayer();
    navigableLayer->setInwardLinks(reachedBarrelLayersMap[dl], reachedForwardLayersMap[dl]);
  }
}

void SimpleNavigationSchool::init ( void  )

private

Definition at line 25 of file SimpleNavigationSchool.cc.

References a, and mps_fire::i.

                                  {
  theAllDetLayersInSystem = &theTracker->allLayers();
  theAllNavigableLayer.resize(theTracker->allLayers().size(), nullptr);

  // Get barrel layers
  for (auto i : theTracker->barrelLayers()) {
    theBarrelLayers.push_back(i);
  }

  // get forward layers
  for (auto i : theTracker->forwardLayers()) {
    theForwardLayers.push_back(i);
  }

  FDLI middle = find_if(theForwardLayers.begin(), theForwardLayers.end(), [](const GeometricSearchDet* a) {
    return a->position().z() >= 0.0;
  });
  theLeftLayers = FDLC(theForwardLayers.begin(), middle);
  theRightLayers = FDLC(middle, theForwardLayers.end());

  SymmetricLayerFinder symFinder(theForwardLayers);

  // only work on positive Z side; negative by mirror symmetry later
  linkBarrelLayers(symFinder);
  linkForwardLayers(symFinder);
  establishInverseRelations();
}

void SimpleNavigationSchool::linkBarrelLayers ( SymmetricLayerFinder &  symFinder )

protectedvirtual

Reimplemented in CosmicNavigationSchool.

Definition at line 78 of file SimpleNavigationSchool.cc.

References mps_fire::i, and SymmetricLayerFinder::mirror().

                                                                             {
  // 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 SimpleNavigationSchool::linkForwardLayers ( SymmetricLayerFinder &  symFinder )

protectedvirtual

Definition at line 138 of file SimpleNavigationSchool.cc.

References g, watchdog::group, mps_fire::i, LogDebug, and SymmetricLayerFinder::mirror().

                                                                              {
  // 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 SimpleNavigationSchool::linkNextBarrelLayer ( ForwardDetLayer const *  fl, BDLC &  reachableBL  )

protectedvirtual

Definition at line 184 of file SimpleNavigationSchool.cc.

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

                                                                                             {
  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.) {
      reachableBL.push_back(*bli);
      return;
    }
  }
}

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

protectedvirtual

Definition at line 105 of file SimpleNavigationSchool.cc.

References Surface::bounds(), Bounds::length(), CosmicsPD_Skims::radius, BarrelDetLayer::specificSurface(), and BarrelDetLayer::surface().

                                                                                         {
  // 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()) {
      //search if there are any sovrapposition between forward layers
      for (FDLI fliNext = fli; fliNext != theRightLayers.end(); fliNext++) {
        if ((**fliNext).position().z() < (**fli).position().z() &&
            (**fliNext).specificSurface().innerRadius() < (**fli).specificSurface().outerRadius()) {
          rightFL.push_back(*fliNext);
          return;
        }
      }
      rightFL.push_back(*fli);
      return;
    }
  }
}

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

protectedvirtual

Definition at line 125 of file SimpleNavigationSchool.cc.

References dataset::end, geometryDiff::epsilon, and mps_fire::i.

                                                                                      {
  // 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 SimpleNavigationSchool::linkNextLayerInGroup ( FDLI  fli, const FDLC &  group, FDLC &  reachableFL  )

protectedvirtual

Definition at line 198 of file SimpleNavigationSchool.cc.

References geometryDiff::epsilon, and mps_fire::i.

                                                                                                {
  // 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 SimpleNavigationSchool::linkOuterGroup ( ForwardDetLayer const *  fl, const FDLC &  group, FDLC &  reachableFL  )

protectedvirtual

Definition at line 226 of file SimpleNavigationSchool.cc.

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

                                                                                                           {
  // insert N layers with Z grater than fl

  ConstFDLI first = find_if(group.begin(), group.end(), [fl](const GeometricSearchDet* a) {
    return a->position().z() >= fl->position().z();
  });
  if (first != group.end()) {
    // Hard-wired constant!!!!!!
    ConstFDLI last = min(first + 7, group.end());

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

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

protectedvirtual

Definition at line 240 of file SimpleNavigationSchool.cc.

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

SimpleNavigationSchool::StateType SimpleNavigationSchool::navigableLayers ( )

overridevirtual

Implements NavigationSchool.

Definition at line 66 of file SimpleNavigationSchool.cc.

References cuy::ib, and mps_fire::result.

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

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

protectedvirtual

Definition at line 247 of file SimpleNavigationSchool.cc.

References geometryDiff::epsilon.

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

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

protectedvirtual

Definition at line 258 of file SimpleNavigationSchool.cc.

References a, b, SplitLinear::begin, dataset::end, mps_fire::i, isDetLessZ(), LogDebug, mps_fire::result, and ForwardDetLayer::specificSurface().

                                                                              {
  // 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, but keeping the ordering in z!
  std::stable_sort(begin, end, [](const ForwardDetLayer* a, const ForwardDetLayer* b) {
    return a->specificSurface().innerRadius() < b->specificSurface().innerRadius();
  });

  // partition in cylinders
  vector<FDLC> result;
  FDLC current;
  current.push_back(*begin);
  for (FDLI i = begin + 1; i != end; i++) {
#ifdef EDM_ML_DEBUG
    LogDebug("TkNavigation") << "(**i).specificSurface().innerRadius()      = " << (**i).specificSurface().innerRadius()
                             << endl
                             << "(**(i-1)).specificSurface().outerRadius()) = "
                             << (**(i - 1)).specificSurface().outerRadius();
    LogDebug("TkNavigation") << "(**i).specificSurface().position().z()      = "
                             << (**i).specificSurface().position().z() << endl
                             << "(**(i-1)).specificSurface().position().z() = "
                             << (**(i - 1)).specificSurface().position().z();
#endif

    // if inner radius of i is larger than outer radius of i-1 then split!
    // FIXME: The solution found for phase2 is a bit dirty, we can do better.
    // For phase2 we compare the EXTENDED pixel with the TID to get the assignment right!
    if ((**i).specificSurface().innerRadius() > (**(i - 1)).specificSurface().outerRadius() ||
        (theTracker->posPixelForwardLayers().back()->specificSurface().position().z() >
             theTracker->posTidLayers().front()->specificSurface().position().z() &&
         (**i).specificSurface().position().z() < (**(i - 1)).specificSurface().position().z())) {
      LogDebug("TkNavigation") << "found break between groups";

      // sort layers in group along Z
      std::stable_sort(current.begin(), current.end(), isDetLessZ);

      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++) {
    std::stable_sort(ivec->begin(), ivec->end(), isDetLessZ);
  }

  return result;
}

Member Data Documentation

BDLC SimpleNavigationSchool::theBarrelLayers

protected

Definition at line 41 of file SimpleNavigationSchool.h.

float SimpleNavigationSchool::theBarrelLength = 0

protected

Definition at line 45 of file SimpleNavigationSchool.h.

BNLCType SimpleNavigationSchool::theBarrelNLC

protected

Definition at line 49 of file SimpleNavigationSchool.h.

FDLC SimpleNavigationSchool::theForwardLayers

protected

Definition at line 42 of file SimpleNavigationSchool.h.

FNLCType SimpleNavigationSchool::theForwardNLC

protected

Definition at line 50 of file SimpleNavigationSchool.h.

FDLC SimpleNavigationSchool::theLeftLayers

protected

Definition at line 44 of file SimpleNavigationSchool.h.

FDLC SimpleNavigationSchool::theRightLayers

protected

Definition at line 43 of file SimpleNavigationSchool.h.