#include <SimpleBarrelNavigableLayer.h>

Inheritance diagram for SimpleBarrelNavigableLayer:

Public Member Functions
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 &theBarrelv, const FDLC &theForwardv, TkLayerLess sorter=TkLayerLess(outsideIn)) override

	SimpleBarrelNavigableLayer (BarrelDetLayer const detLayer, const BDLC &outerBLC, const FDLC &outerLeftFL, const FDLC &outerRightFL, const MagneticField field, float epsilon, bool checkCrossingSide=true)

Public Member Functions inherited from SimpleNavigableLayer
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
virtual void	setDetLayer (DetLayer const *dl)=0

void	setSchool (NavigationSchool const *sh)

virtual	~NavigableLayer ()

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 = 0

Detailed Description

A concrete NavigableLayer for the barrel

Definition at line 12 of file SimpleBarrelNavigableLayer.h.

Constructor & Destructor Documentation

SimpleBarrelNavigableLayer::SimpleBarrelNavigableLayer	(	BarrelDetLayer const *	detLayer,
		const BDLC &	outerBLC,
		const FDLC &	outerLeftFL,
		const FDLC &	outerRightFL,
		const MagneticField *	field,
		float	epsilon,
		bool	checkCrossingSide = `true`
	)

Definition at line 21 of file SimpleBarrelNavigableLayer.cc.

References jetUpdater_cfi::sort, theNegOuterLayers, theOuterBarrelLayers, theOuterLeftForwardLayers, theOuterRightForwardLayers, and thePosOuterLayers.

                                         :
   SimpleNavigableLayer(field,epsilon,checkCrossingSide),
   theDetLayer( detLayer), 
   theOuterBarrelLayers( outerBLC),
   theOuterLeftForwardLayers( outerLeftFL),
   theOuterRightForwardLayers( outerRightFL)
 
 {
   // put barrel and forward layers together
   theNegOuterLayers.reserve( outerBLC.size() + outerLeftFL.size());
   thePosOuterLayers.reserve( outerBLC.size() + outerRightFL.size());
 
   for (ConstBDLI bl=outerBLC.begin(); bl!=outerBLC.end(); bl++) 
     theNegOuterLayers.push_back( *bl);
   thePosOuterLayers = theNegOuterLayers; // barrel part the same
 
   for (ConstFDLI fl=outerLeftFL.begin(); fl!=outerLeftFL.end(); fl++) 
     theNegOuterLayers.push_back( *fl);
   for (ConstFDLI fl=outerRightFL.begin(); fl!=outerRightFL.end(); fl++) 
     thePosOuterLayers.push_back( *fl);
 
   // sort the outer layers 
   sort( theNegOuterLayers.begin(), theNegOuterLayers.end(), TkLayerLess());
   sort( thePosOuterLayers.begin(), thePosOuterLayers.end(), TkLayerLess());
   sort(theOuterBarrelLayers.begin(), theOuterBarrelLayers.end(), TkLayerLess());
   sort(theOuterLeftForwardLayers.begin(), theOuterLeftForwardLayers.end(), TkLayerLess());
   sort(theOuterRightForwardLayers.begin(), theOuterRightForwardLayers.end(), TkLayerLess());
 }

Member Function Documentation

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

overridevirtual

Implements NavigableLayer.

Definition at line 182 of file SimpleBarrelNavigableLayer.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*>();
 }

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

inlineoverridevirtual

Reimplemented from NavigableLayer.

Definition at line 39 of file SimpleBarrelNavigableLayer.h.

References SimpleNavigableLayer::compatibleLayers(), insideOut, and SimpleNavigableLayer::setAdditionalLink().

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

const DetLayer* SimpleBarrelNavigableLayer::detLayer ( ) const

inlineoverridevirtual

Implements NavigableLayer.

Definition at line 47 of file SimpleBarrelNavigableLayer.h.

References outsideIn, NavigableLayer::setDetLayer(), and SimpleNavigableLayer::setInwardLinks().

Referenced by nextLayers().

47 { return theDetLayer;}

SimpleBarrelNavigableLayer::theDetLayer

const BarrelDetLayer * theDetLayer

Definition: SimpleBarrelNavigableLayer.h:53

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

overridevirtual

Implements NavigableLayer.

Definition at line 59 of file SimpleBarrelNavigableLayer.cc.

References GeomDetEnumerators::endcap, mps_fire::i, insideOut, mps_fire::result, theNegInnerLayers, theNegOuterLayers, thePosInnerLayers, and thePosOuterLayers.

 {
   vector<const DetLayer*> result;
   
   // the order is the one in which layers
   // should be checked for a reasonable trajectory
 
   if ( dir == insideOut) {
     result = theNegOuterLayers;
     for ( DLC::const_iterator i=thePosOuterLayers.begin();
       i!=thePosOuterLayers.end(); i++) {
       // avoid duplication of barrel layers
       if ((**i).location() == GeomDetEnumerators::endcap) result.push_back(*i);
     }
   }
   else {
     result = theNegInnerLayers;
     for ( DLC::const_iterator i=thePosInnerLayers.begin();
       i!=thePosInnerLayers.end(); i++) {
       // avoid duplication of barrel layers
       if ((**i).location() == GeomDetEnumerators::endcap) result.push_back(*i);
     }
   }
   return result;
 }

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

overridevirtual

Implements NavigableLayer.

Definition at line 86 of file SimpleBarrelNavigableLayer.cc.

References alongMomentum, detLayer(), FreeTrajectoryState::hasError(), LIKELY, LogDebug, FreeTrajectoryState::momentum(), oppositeToMomentum, FreeTrajectoryState::parameters(), FreeTrajectoryState::position(), position, mps_fire::result, theInnerBarrelLayers, theInnerLeftForwardLayers, theInnerRightForwardLayers, theNegInnerLayers, theNegOuterLayers, theOuterBarrelLayers, theOuterLeftForwardLayers, theOuterRightForwardLayers, thePosInnerLayers, thePosOuterLayers, SimpleNavigableLayer::theSelfSearch, UNLIKELY, 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.
   GlobalVector transversePosition(position.x(), position.y(), 0);
   bool isInOutTrackBarrel  = (transversePosition.dot(momentum)>0);
 
   float zpos = position.z();
   bool isInOutTrackFWD = momentum.z()*zpos>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);
 
   LogDebug("SimpleBarrelNavigableLayer") << "is alongMomentum? " << (dir == alongMomentum) << endl
                      << "isInOutTrackBarrel: " << isInOutTrackBarrel << endl
                      << "isInOutTrackFWD: " << isInOutTrackFWD << endl
                      << "dirOppositeXORisInOutTrackFWD: " << dirOppositeXORisInOutTrackFWD << endl
                      << "dirOppositeXORisInOutTrackBarrel: "<< dirOppositeXORisInOutTrackBarrel << endl;
 
   bool signZmomentumXORdir = (( (momentum.z() > 0) && !(dir == alongMomentum) ) ||
                   (!(momentum.z() > 0) &&  (dir == alongMomentum) )   );
 
 
   if LIKELY( dirOppositeXORisInOutTrackBarrel &&  dirOppositeXORisInOutTrackFWD) {
       if ( signZmomentumXORdir   ) {
     wellInside( ftsWithoutErrors, dir, theNegOuterLayers, result);
       }
       else {
     wellInside( ftsWithoutErrors, dir, thePosOuterLayers, result);
       }
     } else if (!dirOppositeXORisInOutTrackBarrel &&  !dirOppositeXORisInOutTrackFWD){
     if ( signZmomentumXORdir ) {
       wellInside( ftsWithoutErrors, dir, thePosInnerLayers, result);
     }
     else {
       wellInside( ftsWithoutErrors, dir, theNegInnerLayers, result);
     }
   } else if (!dirOppositeXORisInOutTrackBarrel && dirOppositeXORisInOutTrackFWD){
     wellInside(ftsWithoutErrors, dir, theInnerBarrelLayers.begin(), theInnerBarrelLayers.end(), result);    
     
     if (signZmomentumXORdir){   
       wellInside(ftsWithoutErrors, dir, theInnerLeftForwardLayers.begin(), theInnerLeftForwardLayers.end(), result);
       wellInside(ftsWithoutErrors, dir, theOuterLeftForwardLayers.begin(), theOuterLeftForwardLayers.end(), result);    
     }   else {
       wellInside(ftsWithoutErrors, dir, theInnerRightForwardLayers.begin(), theInnerRightForwardLayers.end(), result);
       wellInside(ftsWithoutErrors, dir, theOuterRightForwardLayers.begin(), theOuterRightForwardLayers.end(), result);
     }
   } else {
      if (signZmomentumXORdir){
         wellInside(ftsWithoutErrors, dir, theInnerLeftForwardLayers.begin(), theInnerLeftForwardLayers.end(), result);
      } else {
         wellInside(ftsWithoutErrors, dir, theInnerRightForwardLayers.begin(), theInnerRightForwardLayers.end(), result);
      }  
      wellInside(ftsWithoutErrors, dir, theOuterBarrelLayers.begin(), theOuterBarrelLayers.end(), result);   
   }
 
   bool goingIntoTheBarrel = (!isInOutTrackBarrel && dir==alongMomentum) || (isInOutTrackBarrel && dir==oppositeToMomentum) ;
 
   LogDebug("SimpleBarrelNavigableLayer") << "goingIntoTheBarrel: " << goingIntoTheBarrel;
 
 
   if UNLIKELY(theSelfSearch && result.empty()){
     if (!goingIntoTheBarrel){     
       LogDebug("SimpleBarrelNavigableLayer")<<" state is not going toward the center of the barrel. not adding self search.";}
     else{
       const BarrelDetLayer * bl = reinterpret_cast<const BarrelDetLayer *>(detLayer());      unsigned int before=result.size();
       LogDebug("SimpleBarrelNavigableLayer")<<" I am trying to added myself as a next layer.";
       wellInside(ftsWithoutErrors, dir, bl, result);
       unsigned int after=result.size();
       if (before!=after)
     LogDebug("SimpleBarrelNavigableLayer")<<" I have added myself as a next layer.";
     }
   }
   
   return result;
 }

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

overridevirtual

Implements SimpleNavigableLayer.

Definition at line 237 of file SimpleBarrelNavigableLayer.cc.

References insideOut, GeometricSearchDet::position(), theInnerBarrelLayers, theInnerLeftForwardLayers, theInnerRightForwardLayers, theNegInnerLayers, theNegOuterLayers, theOuterBarrelLayers, theOuterLeftForwardLayers, theOuterRightForwardLayers, thePosInnerLayers, thePosOuterLayers, and PV3DBase< T, PVType, FrameType >::z().

                                                                                                            {
   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);
         theNegOuterLayers.push_back(badditional);
         thePosOuterLayers.push_back(badditional);
         return; 
     }
     theInnerBarrelLayers.push_back(badditional);
     theNegInnerLayers.push_back(badditional);
         thePosInnerLayers.push_back(badditional);   
     return;
   } else if (fadditional){
     double zpos = fadditional->position().z(); 
     if (direction==insideOut){
         if (zpos>0){
             theOuterRightForwardLayers.push_back(fadditional);
             thePosOuterLayers.push_back(fadditional);
             return; 
         }
         theOuterLeftForwardLayers.push_back(fadditional);
         theNegOuterLayers.push_back(fadditional);
         return;
     }
     if (zpos>0){
             theInnerRightForwardLayers.push_back(fadditional);
         thePosInnerLayers.push_back(fadditional);
                 return;
         }
         theInnerLeftForwardLayers.push_back(fadditional);
     theNegInnerLayers.push_back(fadditional);
         return;
   }
   edm::LogError("TkNavigation") << "trying to add neither a ForwardDetLayer nor a BarrelDetLayer";  
   return;   
 }

void SimpleBarrelNavigableLayer::setDetLayer ( const DetLayer * dl )

override

Definition at line 190 of file SimpleBarrelNavigableLayer.cc.

References MessageLogger_cfi::cerr.

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

void SimpleBarrelNavigableLayer::setInwardLinks	(	const BDLC &	theBarrelv,
		const FDLC &	theForwardv,
		TkLayerLess	sorter = `TkLayerLess(outsideIn)`
	)

overridevirtual

Implements SimpleNavigableLayer.

Definition at line 195 of file SimpleBarrelNavigableLayer.cc.

References a, jetUpdater_cfi::sort, MCScenario_CRAFT1_22X::sorter(), theInnerBarrelLayers, theInnerLeftForwardLayers, theInnerRightForwardLayers, theNegInnerLayers, and thePosInnerLayers.

 {
   theInnerBarrelLayers=theBarrelv;
   // sort the inner layers
   sort(theInnerBarrelLayers.begin(), theInnerBarrelLayers.end(),sorter);
 
 
   ConstFDLI middle = find_if( theForwardv.begin(),theForwardv.end(),
           [](const GeometricSearchDet* a){ return a->position().z() >= 0.0; });
   theInnerLeftForwardLayers=FDLC(theForwardv.begin(),middle);
   theInnerRightForwardLayers=FDLC(middle,theForwardv.end());
 
   // sort the inner layers
   sort(theInnerLeftForwardLayers.begin(), theInnerLeftForwardLayers.end(),sorter);
   sort(theInnerRightForwardLayers.begin(), theInnerRightForwardLayers.end(),sorter);
 
 
 
   // put barrel and forward layers together
   theNegInnerLayers.reserve( theInnerBarrelLayers.size() + theInnerLeftForwardLayers.size());
   thePosInnerLayers.reserve( theInnerBarrelLayers.size() + theInnerRightForwardLayers.size());
 
   for (ConstBDLI bl=theInnerBarrelLayers.begin(); bl!=theInnerBarrelLayers.end(); bl++) 
     theNegInnerLayers.push_back( *bl);
   thePosInnerLayers = theNegInnerLayers; // barrel part the same
 
   for (ConstFDLI fl=theInnerLeftForwardLayers.begin(); fl!=theInnerLeftForwardLayers.end(); fl++) 
     theNegInnerLayers.push_back( *fl);
   for (ConstFDLI fl=theInnerRightForwardLayers.begin(); fl!=theInnerRightForwardLayers.end(); fl++) 
     thePosInnerLayers.push_back( *fl);
 
   // sort the inner layers 
   sort( theNegInnerLayers.begin(), theNegInnerLayers.end(), sorter);
   sort( thePosInnerLayers.begin(), thePosInnerLayers.end(), sorter);
   sort(theInnerBarrelLayers.begin(), theInnerBarrelLayers.end(),sorter);
   sort(theInnerLeftForwardLayers.begin(), theInnerLeftForwardLayers.end(),sorter);
   sort(theInnerRightForwardLayers.begin(), theInnerRightForwardLayers.end(),sorter);
 
 }

Member Data Documentation

const BarrelDetLayer* SimpleBarrelNavigableLayer::theDetLayer

private

Definition at line 53 of file SimpleBarrelNavigableLayer.h.

BDLC SimpleBarrelNavigableLayer::theInnerBarrelLayers

private

Definition at line 55 of file SimpleBarrelNavigableLayer.h.

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

FDLC SimpleBarrelNavigableLayer::theInnerLeftForwardLayers

private

Definition at line 60 of file SimpleBarrelNavigableLayer.h.

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

FDLC SimpleBarrelNavigableLayer::theInnerRightForwardLayers

private

Definition at line 61 of file SimpleBarrelNavigableLayer.h.

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

DLC SimpleBarrelNavigableLayer::theNegInnerLayers

private

Definition at line 65 of file SimpleBarrelNavigableLayer.h.

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

DLC SimpleBarrelNavigableLayer::theNegOuterLayers

private

Definition at line 63 of file SimpleBarrelNavigableLayer.h.

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

BDLC SimpleBarrelNavigableLayer::theOuterBarrelLayers

private

Definition at line 54 of file SimpleBarrelNavigableLayer.h.

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

FDLC SimpleBarrelNavigableLayer::theOuterLeftForwardLayers

private

Definition at line 57 of file SimpleBarrelNavigableLayer.h.

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

FDLC SimpleBarrelNavigableLayer::theOuterRightForwardLayers

private

Definition at line 58 of file SimpleBarrelNavigableLayer.h.

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

DLC SimpleBarrelNavigableLayer::thePosInnerLayers

private

Definition at line 66 of file SimpleBarrelNavigableLayer.h.

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

DLC SimpleBarrelNavigableLayer::thePosOuterLayers

private

Definition at line 64 of file SimpleBarrelNavigableLayer.h.

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

Private Attributes
const BarrelDetLayer *	theDetLayer

BDLC	theInnerBarrelLayers

FDLC	theInnerLeftForwardLayers

FDLC	theInnerRightForwardLayers

DLC	theNegInnerLayers

DLC	theNegOuterLayers

BDLC	theOuterBarrelLayers

FDLC	theOuterLeftForwardLayers

FDLC	theOuterRightForwardLayers

DLC	thePosInnerLayers

DLC	thePosOuterLayers

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation