SimpleForwardNavigableLayer Class Reference

#include <SimpleForwardNavigableLayer.h>

Inheritance diagram for SimpleForwardNavigableLayer:

Public Member Functions
std::vector< const DetLayer *>	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection, int &counter) const final
std::vector< const DetLayer * >	compatibleLayers (NavigationDirection direction) const override
std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection dir) const override
const DetLayer *	detLayer () const override
std::vector< const DetLayer * >	nextLayers (NavigationDirection direction) const override
std::vector< const DetLayer * >	nextLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection) const override
void	setAdditionalLink (const DetLayer *, NavigationDirection direction=insideOut) override
void	setDetLayer (const DetLayer *dl) override
void	setInwardLinks (const BDLC &, const FDLC &, TkLayerLess sorter=TkLayerLess(outsideIn)) override
	SimpleForwardNavigableLayer (const ForwardDetLayer *detLayer, const BDLC &outerBL, const FDLC &outerFL, const MagneticField *field, float epsilon, bool checkCrossingSide=true)
Public Member Functions inherited from SimpleNavigableLayer
virtual std::vector< const DetLayer *>	compatibleLayers (NavigationDirection direction) const=0
virtual std::vector< const DetLayer *>	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection, int &counter) const
virtual std::vector< const DetLayer *>	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection) const
std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection, int &counter) const final
void	setCheckCrossingSide (bool docheck)
	SimpleNavigableLayer (const MagneticField *field, float eps, bool checkCrossingSide=true)
Public Member Functions inherited from NavigableLayer
void	setSchool (NavigationSchool const *sh)
virtual	~NavigableLayer ()

Private Attributes
const ForwardDetLayer *	theDetLayer
BDLC	theInnerBarrelLayers
FDLC	theInnerForwardLayers
DLC	theInnerLayers
BDLC	theOuterBarrelLayers
FDLC	theOuterForwardLayers
DLC	theOuterLayers

Additional Inherited Members
Public Types inherited from SimpleNavigableLayer
typedef std::vector< const BarrelDetLayer * >	BDLC
typedef std::vector< const DetLayer * >	DLC
typedef std::vector< const ForwardDetLayer * >	FDLC
Public Attributes inherited from SimpleNavigableLayer
bool	theSelfSearch
Protected Types inherited from SimpleNavigableLayer
typedef BDLC::iterator	BDLI
typedef BDLC::const_iterator	ConstBDLI
typedef FDLC::const_iterator	ConstFDLI
typedef FDLC::iterator	FDLI
typedef TrajectoryStateOnSurface	TSOS
Protected Member Functions inherited from SimpleNavigableLayer
TSOS	crossingState (const FreeTrajectoryState &fts, PropagationDirection dir) const
AnalyticalPropagator	propagator (PropagationDirection dir) const
bool	wellInside (const FreeTrajectoryState &fts, PropagationDirection dir, const BarrelDetLayer *bl, DLC &result) const
bool	wellInside (const FreeTrajectoryState &fts, PropagationDirection dir, const ForwardDetLayer *bl, DLC &result) const
bool	wellInside (const FreeTrajectoryState &fts, PropagationDirection dir, ConstBDLI begin, ConstBDLI end, DLC &result) const
bool	wellInside (const FreeTrajectoryState &fts, PropagationDirection dir, const DLC &layers, DLC &result) const
bool	wellInside (const FreeTrajectoryState &fts, PropagationDirection dir, ConstFDLI begin, ConstFDLI end, DLC &result) const
Protected Attributes inherited from SimpleNavigableLayer
bool	theCheckCrossingSide
float	theEpsilon
const MagneticField *	theField
Protected Attributes inherited from NavigableLayer
NavigationSchool const *	school = nullptr

Detailed Description

A concrete NavigableLayer for the forward

Definition at line 10 of file SimpleForwardNavigableLayer.h.

Constructor & Destructor Documentation

SimpleForwardNavigableLayer()

SimpleForwardNavigableLayer::SimpleForwardNavigableLayer	(	const ForwardDetLayer *	detLayer,
		const BDLC &	outerBL,
		const FDLC &	outerFL,
		const MagneticField *	field,
		float	epsilon,
		bool	checkCrossingSide = true
	)

Definition at line 16 of file SimpleForwardNavigableLayer.cc.

References jetUpdater_cfi::sort, theOuterBarrelLayers, theOuterForwardLayers, and theOuterLayers.

    : SimpleNavigableLayer(field, epsilon, checkCrossingSide),
      theDetLayer(detLayer),
      theOuterBarrelLayers(outerBL),
      theInnerBarrelLayers(0),
      theOuterForwardLayers(outerFL),
      theInnerForwardLayers(0),
      theOuterLayers(0),
      theInnerLayers(0) {
  // put barrel and forward layers together
  theOuterLayers.reserve(outerBL.size() + outerFL.size());
  for (ConstBDLI bl = outerBL.begin(); bl != outerBL.end(); bl++)
    theOuterLayers.push_back(*bl);
  for (ConstFDLI fl = outerFL.begin(); fl != outerFL.end(); fl++)
    theOuterLayers.push_back(*fl);

  // sort the outer layers
  sort(theOuterLayers.begin(), theOuterLayers.end(), TkLayerLess());
  sort(theOuterForwardLayers.begin(), theOuterForwardLayers.end(), TkLayerLess());
  sort(theOuterBarrelLayers.begin(), theOuterBarrelLayers.end(), TkLayerLess());
}

Member Function Documentation

compatibleLayers() [1/3]

std::vector< const DetLayer * > SimpleNavigableLayer::compatibleLayers

final

Definition at line 214 of file SimpleNavigableLayer.cc.

                                                                                        {
  typedef std::vector<const DetLayer*> Lvect;
  typedef std::set<const DetLayer*> Lset;

  //initiate the first iteration
  const Lvect& someLayers = nextLayers(fts, timeDirection);
  if (someLayers.empty()) {
    LogDebug("SimpleNavigableLayer") << "Number of compatible layers: " << 0;
    return someLayers;
  }

  Lset collect;                     //a container of unique instances. to avoid duplicates
  Lset layerToTry, nextLayerToTry;  //set used for iterations
  layerToTry.insert(someLayers.begin(), someLayers.end());

  while (!layerToTry.empty() && (counter++) <= 150) {
    LogDebug("SimpleNavigableLayer") << counter << "] going to check on : " << layerToTry.size() << " next layers.";
    //clear this set first, it will be swaped with layerToTry
    nextLayerToTry.clear();
    for (auto toTry : layerToTry) {
      //add the layer you tried.
      LogDebug("SimpleNavigableLayer") << counter << "] adding layer with pointer: " << (toTry)
                                       << " first detid: " << (toTry)->basicComponents().front()->geographicalId();
      if (!collect.insert(toTry).second)
        continue;

      //find the next layers from it
      Lvect&& nextLayers = school->nextLayers(*toTry, fts, timeDirection);
      LogDebug("SimpleNavigableLayer") << counter << "] this layer has : " << nextLayers.size() << " next layers.";
      nextLayerToTry.insert(nextLayers.begin(), nextLayers.end());
    }  // layerToTry
    //swap now that you where to go next.
    layerToTry.swap(nextLayerToTry);
  }
  if (counter >= 150) {
    edm::LogWarning("SimpleNavigableLayer") << "WARNING: compatibleLayers() more than 150 iterations!!! Bailing out..";
    counter = -1;
    return Lvect();
  }

  LogDebug("SimpleNavigableLayer") << "Number of compatible layers: " << collect.size();

  return Lvect(collect.begin(), collect.end());
}

compatibleLayers() [2/3]

vector< const DetLayer * > SimpleForwardNavigableLayer::compatibleLayers ( NavigationDirection  direction ) const

overridevirtual

Implements NavigableLayer.

Definition at line 101 of file SimpleForwardNavigableLayer.cc.

                                                                                                   {
  edm::LogError("TkNavigation") << "ERROR: compatibleLayers() method used without all reachableLayers are set";
  throw DetLayerException("compatibleLayers() method used without all reachableLayers are set");
  return vector<const DetLayer*>();
}

compatibleLayers() [3/3]

std::vector<const DetLayer*> SimpleForwardNavigableLayer::compatibleLayers ( const FreeTrajectoryState &  fts, PropagationDirection  dir  ) const

inlineoverridevirtual

Reimplemented from NavigableLayer.

Definition at line 29 of file SimpleForwardNavigableLayer.h.

References SimpleNavigableLayer::compatibleLayers(), and DeadROC_duringRun::dir.

                                                                                         {
    int counter = 0;
    return SimpleNavigableLayer::compatibleLayers(fts, dir, counter);
  }

detLayer()

const DetLayer* SimpleForwardNavigableLayer::detLayer ( ) const

inlineoverridevirtual

Implements NavigableLayer.

Definition at line 37 of file SimpleForwardNavigableLayer.h.

{ return theDetLayer; }

nextLayers() [1/2]

vector< const DetLayer * > SimpleForwardNavigableLayer::nextLayers ( NavigationDirection  direction ) const

overridevirtual

Implements NavigableLayer.

Definition at line 43 of file SimpleForwardNavigableLayer.cc.

References DeadROC_duringRun::dir, insideOut, mps_fire::result, theInnerLayers, and theOuterLayers.

                                                                                             {
  vector<const DetLayer*> result;

  // the order is the one in which layers
  // should be checked for a reasonable trajectory

  if (dir == insideOut) {
    return theOuterLayers;
  } else {
    return theInnerLayers;
  }

  return result;
}

nextLayers() [2/2]

vector< const DetLayer * > SimpleForwardNavigableLayer::nextLayers ( const FreeTrajectoryState &  fts, PropagationDirection  timeDirection  ) const

overridevirtual

Implements NavigableLayer.

Definition at line 58 of file SimpleForwardNavigableLayer.cc.

References DeadROC_duringRun::dir, FreeTrajectoryState::hasError(), LIKELY, FreeTrajectoryState::momentum(), oppositeToMomentum, FreeTrajectoryState::parameters(), FreeTrajectoryState::position(), position, mps_fire::result, theInnerBarrelLayers, theInnerForwardLayers, theInnerLayers, theOuterBarrelLayers, theOuterForwardLayers, theOuterLayers, and SimpleNavigableLayer::wellInside().

                                                                                                {
  // This method contains the sequence in which the layers are tested.
  // The iteration stops as soon as a layer contains the propagated state
  // within epsilon

  vector<const DetLayer*> result;

  FreeTrajectoryState ftsWithoutErrors = (fts.hasError()) ? FreeTrajectoryState(fts.parameters()) : fts;

  auto const position = fts.position();
  auto const momentum = fts.momentum();

  //establish whether the tracks is crossing the tracker from outer layers to inner ones
  //or from inner to outer
  float zpos = position.z();
  bool isInOutTrackFWD = momentum.z() * zpos > 0;
  GlobalVector transversePosition(position.x(), position.y(), 0);
  bool isInOutTrackBarrel = (transversePosition.dot(momentum) > 0);

  //establish whether inner or outer layers are crossed after propagation, according
  //to BOTH propagationDirection AND track momentum
  bool dirOppositeXORisInOutTrackBarrel =
      (!(dir == oppositeToMomentum) && isInOutTrackBarrel) || ((dir == oppositeToMomentum) && !isInOutTrackBarrel);
  bool dirOppositeXORisInOutTrackFWD =
      (!(dir == oppositeToMomentum) && isInOutTrackFWD) || ((dir == oppositeToMomentum) && !isInOutTrackFWD);
  //bool dirOppositeXORisInOutTrack = ( !(dir == oppositeToMomentum) && isInOutTrack) || ( (dir == oppositeToMomentum) && !isInOutTrack);

  if LIKELY (dirOppositeXORisInOutTrackFWD && dirOppositeXORisInOutTrackBarrel) {  //standard tracks
    wellInside(ftsWithoutErrors, dir, theOuterLayers, result);
  } else if (!dirOppositeXORisInOutTrackFWD && !dirOppositeXORisInOutTrackBarrel) {  // !dirOppositeXORisInOutTrack
    wellInside(ftsWithoutErrors, dir, theInnerLayers, result);
  } else if (!dirOppositeXORisInOutTrackFWD && dirOppositeXORisInOutTrackBarrel) {
    wellInside(ftsWithoutErrors, dir, theInnerForwardLayers.begin(), theInnerForwardLayers.end(), result);
    wellInside(ftsWithoutErrors, dir, theOuterBarrelLayers.begin(), theOuterBarrelLayers.end(), result);
  } else {
    wellInside(ftsWithoutErrors, dir, theInnerBarrelLayers.begin(), theInnerBarrelLayers.end(), result);
    wellInside(ftsWithoutErrors, dir, theOuterForwardLayers.begin(), theOuterForwardLayers.end(), result);
  }

  return result;
}

setAdditionalLink()

void SimpleForwardNavigableLayer::setAdditionalLink ( const DetLayer *  additional, NavigationDirection  direction = insideOut  )

overridevirtual

Implements SimpleNavigableLayer.

Definition at line 128 of file SimpleForwardNavigableLayer.cc.

References insideOut, theInnerBarrelLayers, theInnerForwardLayers, theInnerLayers, theOuterBarrelLayers, theOuterForwardLayers, and theOuterLayers.

Referenced by CosmicNavigationSchool::buildAdditionalForwardLinks().

                                                                                                             {
  const ForwardDetLayer* fadditional = dynamic_cast<const ForwardDetLayer*>(additional);
  const BarrelDetLayer* badditional = dynamic_cast<const BarrelDetLayer*>(additional);
  if (badditional) {
    if (direction == insideOut) {
      theOuterBarrelLayers.push_back(badditional);
      theOuterLayers.push_back(badditional);
      return;
    }
    theInnerBarrelLayers.push_back(badditional);
    theInnerLayers.push_back(badditional);
    return;
  } else if (fadditional) {
    if (direction == insideOut) {
      theOuterForwardLayers.push_back(fadditional);
      theOuterLayers.push_back(badditional);
      return;
    }
    theInnerForwardLayers.push_back(fadditional);
    theInnerLayers.push_back(badditional);
    return;
  }
  edm::LogError("TkNavigation") << "trying to add neither a ForwardDetLayer nor a BarrelDetLayer";
  return;
}

setDetLayer()

void SimpleForwardNavigableLayer::setDetLayer ( const DetLayer *  dl )

overridevirtual

Implements NavigableLayer.

Definition at line 107 of file SimpleForwardNavigableLayer.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr.

                                                                {
  cerr << "Warning: SimpleForwardNavigableLayer::setDetLayer called." << endl << "This should never happen!" << endl;
}

setInwardLinks()

void SimpleForwardNavigableLayer::setInwardLinks ( const BDLC &  innerBL, const FDLC &  innerFL, TkLayerLess  sorter = TkLayerLess(outsideIn)  )

overridevirtual

Implements SimpleNavigableLayer.

Definition at line 111 of file SimpleForwardNavigableLayer.cc.

References jetUpdater_cfi::sort, MCScenario_CRAFT1_22X::sorter(), theInnerBarrelLayers, theInnerForwardLayers, and theInnerLayers.

                                                                                                             {
  theInnerBarrelLayers = innerBL;
  theInnerForwardLayers = innerFL;

  theInnerLayers.clear();
  theInnerLayers.reserve(innerBL.size() + innerFL.size());
  for (ConstBDLI bl = innerBL.begin(); bl != innerBL.end(); bl++)
    theInnerLayers.push_back(*bl);
  for (ConstFDLI fl = innerFL.begin(); fl != innerFL.end(); fl++)
    theInnerLayers.push_back(*fl);

  // sort the inner layers
  sort(theInnerLayers.begin(), theInnerLayers.end(), sorter);
  sort(theInnerForwardLayers.begin(), theInnerForwardLayers.end(), sorter);
  sort(theInnerBarrelLayers.begin(), theInnerBarrelLayers.end(), sorter);
}

Member Data Documentation

theDetLayer

const ForwardDetLayer* SimpleForwardNavigableLayer::theDetLayer

private

Definition at line 43 of file SimpleForwardNavigableLayer.h.

theInnerBarrelLayers

BDLC SimpleForwardNavigableLayer::theInnerBarrelLayers

private

Definition at line 45 of file SimpleForwardNavigableLayer.h.

Referenced by nextLayers(), setAdditionalLink(), and setInwardLinks().

theInnerForwardLayers

FDLC SimpleForwardNavigableLayer::theInnerForwardLayers

private

Definition at line 48 of file SimpleForwardNavigableLayer.h.

Referenced by nextLayers(), setAdditionalLink(), and setInwardLinks().

theInnerLayers

DLC SimpleForwardNavigableLayer::theInnerLayers

private

Definition at line 51 of file SimpleForwardNavigableLayer.h.

Referenced by nextLayers(), setAdditionalLink(), and setInwardLinks().

theOuterBarrelLayers

BDLC SimpleForwardNavigableLayer::theOuterBarrelLayers

private

Definition at line 44 of file SimpleForwardNavigableLayer.h.

Referenced by nextLayers(), setAdditionalLink(), and SimpleForwardNavigableLayer().

theOuterForwardLayers

FDLC SimpleForwardNavigableLayer::theOuterForwardLayers

private

Definition at line 47 of file SimpleForwardNavigableLayer.h.

Referenced by nextLayers(), setAdditionalLink(), and SimpleForwardNavigableLayer().

theOuterLayers

DLC SimpleForwardNavigableLayer::theOuterLayers

private

Definition at line 50 of file SimpleForwardNavigableLayer.h.

Referenced by nextLayers(), setAdditionalLink(), and SimpleForwardNavigableLayer().