CosmicNavigationSchool Class Reference

Inheritance diagram for CosmicNavigationSchool:

Classes
class	CosmicNavigationSchoolConfiguration

Public Member Functions
void	build (const GeometricSearchTracker *theTracker, const MagneticField *field, const CosmicNavigationSchoolConfiguration conf)
	CosmicNavigationSchool (const GeometricSearchTracker *theTracker, const MagneticField *field)
	~CosmicNavigationSchool () override
Public Member Functions inherited from SimpleNavigationSchool
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 Member Functions
	CosmicNavigationSchool ()
Protected Member Functions inherited from SimpleNavigationSchool
virtual float	barrelLength ()
void	cleanMemory ()
virtual void	establishInverseRelations ()
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)

Private Member Functions
void	buildAdditionalBarrelLinks ()
void	buildAdditionalForwardLinks (SymmetricLayerFinder &symFinder)
virtual void	establishInverseRelations ()
void	establishInverseRelations (SymmetricLayerFinder &symFinder)
void	linkBarrelLayers (SymmetricLayerFinder &symFinder) override

Additional Inherited Members
Public Types inherited from NavigationSchool
typedef std::vector< NavigableLayer * >	StateType
Protected Types inherited from SimpleNavigationSchool
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 Attributes inherited from SimpleNavigationSchool
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

Detailed Description

Concrete navigation school for cosmics in the Tracker

Definition at line 15 of file CosmicNavigationSchool.cc.

Constructor & Destructor Documentation

CosmicNavigationSchool() [1/2]

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

Definition at line 90 of file CosmicNavigationSchool.cc.

                                                                           {
  build(theInputTracker, field, CosmicNavigationSchoolConfiguration());
}

References newFWLiteAna::build.

~CosmicNavigationSchool()

CosmicNavigationSchool::~CosmicNavigationSchool ( )

inlineoverride

Definition at line 18 of file CosmicNavigationSchool.cc.

{ cleanMemory(); }

References SimpleNavigationSchool::cleanMemory().

CosmicNavigationSchool() [2/2]

CosmicNavigationSchool::CosmicNavigationSchool ( )

inlineprotected

Definition at line 48 of file CosmicNavigationSchool.cc.

{}

Member Function Documentation

build()

void CosmicNavigationSchool::build	(	const GeometricSearchTracker *	theTracker,
		const MagneticField *	field,
		const CosmicNavigationSchoolConfiguration	conf
	)

Definition at line 95 of file CosmicNavigationSchool.cc.

                                                                                   {
  LogTrace("CosmicNavigationSchool") << "*********Running CosmicNavigationSchool***********";
  theBarrelLength = 0;
  theField = field;
  theTracker = theInputTracker;
 
  theAllDetLayersInSystem = &theInputTracker->allLayers();
  theAllNavigableLayer.resize(theInputTracker->allLayers().size(), nullptr);
 
  // Get barrel layers
  vector<BarrelDetLayer const*> const& blc = theTracker->barrelLayers();
  for (auto i = blc.begin(); i != blc.end(); i++) {
    if (conf.noPXB && GeomDetEnumerators::isTrackerPixel((*i)->subDetector()))
      continue;
    if (conf.noTOB && (*i)->subDetector() == GeomDetEnumerators::TOB)
      continue;
    if (conf.noTIB && (*i)->subDetector() == GeomDetEnumerators::TIB)
      continue;
    theBarrelLayers.push_back((*i));
  }
 
  // get forward layers
  vector<ForwardDetLayer const*> const& flc = theTracker->forwardLayers();
  for (auto i = flc.begin(); i != flc.end(); i++) {
    if (conf.noPXF && GeomDetEnumerators::isTrackerPixel((*i)->subDetector()))
      continue;
    if (conf.noTEC && (*i)->subDetector() == GeomDetEnumerators::TEC)
      continue;
    if (conf.noTID && (*i)->subDetector() == GeomDetEnumerators::TID)
      continue;
    theForwardLayers.push_back((*i));
  }
 
  FDLI middle =
      find_if(theForwardLayers.begin(), theForwardLayers.end(), [](auto const* 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(symFinder);
 
  if (conf.self) {
    // set the self search by hand
    //   NavigationSetter setter(*this);
 
    //add TOB1->TOB1 inward link
    const std::vector<const BarrelDetLayer*>& tobL = theInputTracker->tobLayers();
    if (!tobL.empty()) {
      if (conf.allSelf) {
        LogDebug("CosmicNavigationSchool") << " adding all TOB self search.";
        for (auto lIt = tobL.begin(); lIt != tobL.end(); ++lIt)
          dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[(*lIt)->seqNum()])->theSelfSearch = true;
      } else {
        SimpleNavigableLayer* navigableLayer =
            dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[tobL.front()->seqNum()]);
        LogDebug("CosmicNavigationSchool") << " adding TOB1 to TOB1.";
        navigableLayer->theSelfSearch = true;
      }
    }
    const std::vector<const BarrelDetLayer*>& tibL = theInputTracker->tibLayers();
    if (!tibL.empty()) {
      if (conf.allSelf) {
        LogDebug("CosmicNavigationSchool") << " adding all TIB self search.";
        for (auto lIt = tibL.begin(); lIt != tibL.end(); ++lIt)
          dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[(*lIt)->seqNum()])->theSelfSearch = true;
      } else {
        SimpleNavigableLayer* navigableLayer =
            dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[tibL.front()->seqNum()]);
        LogDebug("CosmicNavigationSchool") << " adding tib1 to tib1.";
        navigableLayer->theSelfSearch = true;
      }
    }
    const std::vector<const BarrelDetLayer*>& pxbL = theInputTracker->pixelBarrelLayers();
    if (!pxbL.empty()) {
      if (conf.allSelf) {
        LogDebug("CosmicNavigationSchool") << " adding all PXB self search.";
        for (auto lIt = pxbL.begin(); lIt != pxbL.end(); ++lIt)
          dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[(*lIt)->seqNum()])->theSelfSearch = true;
      } else {
        SimpleNavigableLayer* navigableLayer =
            dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[pxbL.front()->seqNum()]);
        LogDebug("CosmicNavigationSchool") << " adding pxb1 to pxb1.";
        navigableLayer->theSelfSearch = true;
      }
    }
  }
}

References a, GeometricSearchTracker::allLayers(), CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::allSelf, mps_fire::i, GeomDetEnumerators::isTrackerPixel(), LogDebug, LogTrace, CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::noPXB, CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::noPXF, CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::noTEC, CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::noTIB, CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::noTID, CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::noTOB, GeometricSearchTracker::pixelBarrelLayers(), CosmicNavigationSchool::CosmicNavigationSchoolConfiguration::self, GeomDetEnumerators::TEC, SimpleNavigableLayer::theSelfSearch, GeomDetEnumerators::TIB, GeometricSearchTracker::tibLayers(), GeomDetEnumerators::TID, GeomDetEnumerators::TOB, and GeometricSearchTracker::tobLayers().

buildAdditionalBarrelLinks()

void CosmicNavigationSchool::buildAdditionalBarrelLinks ( )

private

Definition at line 261 of file CosmicNavigationSchool.cc.

                                                        {
  for (auto i = theBarrelLayers.begin(); i != theBarrelLayers.end(); i++) {
    SimpleNavigableLayer* navigableLayer = dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[(*i)->seqNum()]);
    if (i + 1 != theBarrelLayers.end())
      navigableLayer->setAdditionalLink(*(i + 1), outsideIn);
  }
}

References mps_fire::i, outsideIn, and SimpleNavigableLayer::setAdditionalLink().

buildAdditionalForwardLinks()

void CosmicNavigationSchool::buildAdditionalForwardLinks ( SymmetricLayerFinder &  symFinder )

private

Definition at line 269 of file CosmicNavigationSchool.cc.

                                                                                        {
  //the first layer of FPIX should not check the crossing side (since there are no inner layers to be tryed first)
  SimpleNavigableLayer* firstR =
      dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[theRightLayers.front()->seqNum()]);
  SimpleNavigableLayer* firstL =
      dynamic_cast<SimpleNavigableLayer*>(theAllNavigableLayer[theLeftLayers.front()->seqNum()]);
  firstR->setCheckCrossingSide(false);
  firstL->setCheckCrossingSide(false);
 
  for (auto i : theRightLayers) {
    //look for first bigger barrel layer and link to it outsideIn
    SimpleForwardNavigableLayer* nfl = dynamic_cast<SimpleForwardNavigableLayer*>(theAllNavigableLayer[(i)->seqNum()]);
    SimpleForwardNavigableLayer* mnfl =
        dynamic_cast<SimpleForwardNavigableLayer*>(theAllNavigableLayer[symFinder.mirror(i)->seqNum()]);
    for (auto j : theBarrelLayers) {
      if ((i)->specificSurface().outerRadius() < (j)->specificSurface().radius() &&
          fabs((i)->specificSurface().position().z()) < (j)->surface().bounds().length() / 2.) {
        nfl->setAdditionalLink(j, outsideIn);
        mnfl->setAdditionalLink(j, outsideIn);
        break;
      }
    }
  }
}

References mps_fire::i, dqmiolumiharvest::j, SymmetricLayerFinder::mirror(), outsideIn, position, DetLayer::seqNum(), SimpleForwardNavigableLayer::setAdditionalLink(), and SimpleNavigableLayer::setCheckCrossingSide().

establishInverseRelations() [1/2]

void SimpleNavigationSchool::establishInverseRelations

private

Definition at line 324 of file SimpleNavigationSchool.cc.

                                                       {
  // 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]);
  }
}

establishInverseRelations() [2/2]

void CosmicNavigationSchool::establishInverseRelations ( SymmetricLayerFinder &  symFinder )

private

Definition at line 219 of file CosmicNavigationSchool.cc.

                                                                                      {
  //again: standard part is identical to SimpleNavigationSchool one.
  //After the standard link, special outsideIn links are added
 
  // 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);
    }
  }
 
  for (auto nl : theAllNavigableLayer) {
    if (!nl)
      continue;
    auto navigableLayer = static_cast<SimpleNavigableLayer*>(nl);
    auto dl = nl->detLayer();
    navigableLayer->setInwardLinks(reachedBarrelLayersMap[dl], reachedForwardLayersMap[dl]);
  }
 
  //buildAdditionalBarrelLinks();
  buildAdditionalForwardLinks(symFinder);
}

References mps_fire::i, and insideOut.

linkBarrelLayers()

void CosmicNavigationSchool::linkBarrelLayers ( SymmetricLayerFinder &  symFinder )

overrideprivatevirtual

Reimplemented from SimpleNavigationSchool.

Definition at line 189 of file CosmicNavigationSchool.cc.

                                                                             {
  //identical to the SimpleNavigationSchool one, but it allows crossing over the tracker
  //is some non-standard link is needed, it should probably be added here
 
  // 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., false));
  }
}

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