#include <MuonNavigationSchool.h>

Inheritance diagram for MuonNavigationSchool:

Classes
struct	delete_layer

Public Member Functions
	MuonNavigationSchool (const MuonDetLayerGeometry *, bool enableRPC=true)
	Constructor. More...

virtual StateType	navigableLayers () override
	return navigable layers, from base class More...

	~MuonNavigationSchool ()
	Destructor. More...

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
	NavigationDirection. More...

virtual	~NavigationSchool ()

Private Member Functions
void	addBarrelLayer (const BarrelDetLayer *)
	add barrel layer More...

void	addEndcapLayer (const ForwardDetLayer *)
	add endcap layer (backward and forward) More...

float	calculateEta (const float &, const float &) const
	calculate pseudorapidity from r and z More...

void	createInverseLinks ()
	establish inward links More...

void	linkBarrelLayers ()
	link barrel layers More...

void	linkEndcapLayers (const MapE &, std::vector< MuonForwardNavigableLayer * > &)
	link endcap layers More...

Private Attributes
MapE	theBackwardLayers
	+z endcap More...

std::vector < MuonForwardNavigableLayer * >	theBackwardNLC

MapB	theBarrelLayers

std::vector < MuonBarrelNavigableLayer * >	theBarrelNLC
	-z endcap More...

MapE	theForwardLayers
	barrel More...

std::vector < MuonForwardNavigableLayer * >	theForwardNLC

const MuonDetLayerGeometry *	theMuonDetLayerGeometry

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

Protected Member Functions inherited from NavigationSchool
void	setState (const StateType &state)

Protected Attributes inherited from NavigationSchool
const std::vector< const DetLayer * > *	theAllDetLayersInSystem

StateType	theAllNavigableLayer

Detailed Description

Description: Navigation school for the muon system This class defines which DetLayers are reacheable from each Muon DetLayer (DT, CSC and RPC). The reacheableness is based on an eta range criteria.

Author: : Stefano Lacaprara - INFN Padova stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t

Modification:

Chang Liu: The class links maps for nextLayers and compatibleLayers in the same time.

Definition at line 34 of file MuonNavigationSchool.h.

Constructor & Destructor Documentation

MuonNavigationSchool::MuonNavigationSchool	(	const MuonDetLayerGeometry *	muonLayout,
		bool	enableRPC = `true`
	)

Constructor.

Definition at line 38 of file MuonNavigationSchool.cc.

References addBarrelLayer(), addEndcapLayer(), MuonDetLayerGeometry::allBarrelLayers(), MuonDetLayerGeometry::allCSCLayers(), MuonDetLayerGeometry::allDTLayers(), MuonDetLayerGeometry::allEndcapLayers(), MuonDetLayerGeometry::allLayers(), Reference_intrackfit_cff::barrel, createInverseLinks(), Reference_intrackfit_cff::endcap, Exception, i, linkBarrelLayers(), linkEndcapLayers(), NavigationSchool::theAllDetLayersInSystem, NavigationSchool::theAllNavigableLayer, theBackwardLayers, theBackwardNLC, theForwardLayers, and theForwardNLC.

                                                                                                    : theMuonDetLayerGeometry(muonLayout) {
 
   theAllDetLayersInSystem=&muonLayout->allLayers(); 
   theAllNavigableLayer.resize(muonLayout->allLayers().size(),nullptr);
 
 
 
   // get all barrel DetLayers (DT + optional RPC) 
   vector<const DetLayer*> barrel;
   if ( enableRPC ) barrel = muonLayout->allBarrelLayers();
   else barrel = muonLayout->allDTLayers();
 
   for ( auto i = barrel.begin(); i != barrel.end(); i++ ) {
     const BarrelDetLayer* mbp = dynamic_cast<const BarrelDetLayer*>(*i);
     if ( mbp == 0 ) throw cms::Exception("MuonNavigationSchool", "Bad BarrelDetLayer");
     addBarrelLayer(mbp);
   }
 
   // get all endcap DetLayers (CSC + optional RPC)
   vector<const DetLayer*> endcap;
   if ( enableRPC ) endcap = muonLayout->allEndcapLayers();
   else endcap = muonLayout->allCSCLayers();
 
   for ( auto i = endcap.begin(); i != endcap.end(); i++ ) {
     const ForwardDetLayer* mep = dynamic_cast<const ForwardDetLayer*>(*i);
     if ( mep == 0 ) throw cms::Exception("MuonNavigationSchool", "Bad ForwardDetLayer");
     addEndcapLayer(mep);
   }
 
   // create outward links for all DetLayers
   linkBarrelLayers();
   linkEndcapLayers(theForwardLayers,theForwardNLC);
   linkEndcapLayers(theBackwardLayers,theBackwardNLC);
 
   // create inverse links
   createInverseLinks();
 
 }

MuonNavigationSchool::~MuonNavigationSchool ( )

Destructor.

Definition at line 79 of file MuonNavigationSchool.cc.

References theBackwardNLC, theBarrelNLC, and theForwardNLC.

                                             {
 
    for_each(theBarrelNLC.begin(),theBarrelNLC.end(), delete_layer());
    for_each(theForwardNLC.begin(),theForwardNLC.end(), delete_layer());
    for_each(theBackwardNLC.begin(),theBackwardNLC.end(), delete_layer());
 
 }

Member Function Documentation

void MuonNavigationSchool::addBarrelLayer ( const BarrelDetLayer * mbp )

private

add barrel layer

create barrel layer map

Definition at line 115 of file MuonNavigationSchool.cc.

References calculateEta(), CosmicsPD_Skims::radius, BarrelDetLayer::specificSurface(), and theBarrelLayers.

Referenced by MuonNavigationSchool().

                                                                    {
 
   const BoundCylinder& bc = mbp->specificSurface();
   float radius = bc.radius();
   float length = bc.bounds().length()/2.;
 
   float eta_max = calculateEta(radius, length);
   float eta_min = -eta_max;
 
   theBarrelLayers[mbp] = MuonEtaRange(eta_max, eta_min);
 
 }

void MuonNavigationSchool::addEndcapLayer ( const ForwardDetLayer * mep )

private

add endcap layer (backward and forward)

create forwrad/backward layer maps

Definition at line 130 of file MuonNavigationSchool.cc.

References calculateEta(), ForwardDetLayer::specificSurface(), theBackwardLayers, theForwardLayers, and detailsBasic3DVector::z.

Referenced by MuonNavigationSchool().

                                                                     {
 
   const BoundDisk& bd = mep->specificSurface();
   float outRadius = bd.outerRadius();
   float inRadius = bd.innerRadius();
   float thick = bd.bounds().length()/2.;
   float z = bd.position().z();
 
   if ( z > 0. ) {
     float eta_min = calculateEta(outRadius, z-thick);
     float eta_max = calculateEta(inRadius, z+thick);
     theForwardLayers[mep] = MuonEtaRange(eta_max, eta_min);
   } else {
     float eta_max = calculateEta(outRadius, z+thick);
     float eta_min = calculateEta(inRadius, z-thick);
     theBackwardLayers[mep] = MuonEtaRange(eta_max, eta_min);
   }
 
 }

float MuonNavigationSchool::calculateEta	(	const float &	r,
		const float &	z
	)		const

private

calculate pseudorapidity from r and z

Definition at line 152 of file MuonNavigationSchool.cc.

References cmsBatch::log, and funct::tan().

Referenced by addBarrelLayer(), and addEndcapLayer().

                                                                              {
 
   if ( z > 0 ) return -log((tan(atan(r/z)/2.)));
   return log(-(tan(atan(r/z)/2.)));
 
 }

void MuonNavigationSchool::createInverseLinks ( )

private

establish inward links

create inverse links (i.e. inwards)

Definition at line 263 of file MuonNavigationSchool.cc.

Referenced by MuonNavigationSchool().

                                                {
 
   // set outward link
   // NavigationSetter setter(*this);
 
   setState(navigableLayers());
 
 
   // find for each layer which are the layers pointing to it
   typedef map<const DetLayer*, MapB, less<const DetLayer*> > BarrelMapType;
   typedef map<const DetLayer*, MapE, less<const DetLayer*> > ForwardMapType;
 
   // map of all DetLayers which can reach a specific DetLayer
   BarrelMapType reachedBarrelLayersMap;
   ForwardMapType reachedForwardLayersMap;
 
   // map of all DetLayers which is compatible with a specific DetLayer
   BarrelMapType compatibleBarrelLayersMap;
   ForwardMapType compatibleForwardLayersMap;
 
   // collect all reacheable layers starting from a barrel layer
   for ( MapBI bli  = theBarrelLayers.begin(); 
               bli != theBarrelLayers.end(); bli++ ) {
     // barrel
     MuonBarrelNavigableLayer* mbnl =
       dynamic_cast<MuonBarrelNavigableLayer*>(theAllNavigableLayer[((*bli).first)->seqNum()]);
     MapB reacheableB = mbnl->getOuterBarrelLayers();
     for (MapBI i = reacheableB.begin(); i != reacheableB.end(); i++ ) {
       reachedBarrelLayersMap[(*i).first].insert(*bli);
     }
     MapB compatibleB = mbnl->getAllOuterBarrelLayers();
     for (MapBI i = compatibleB.begin(); i != compatibleB.end(); i++ ) {
       compatibleBarrelLayersMap[(*i).first].insert(*bli);
     }
     MapE reacheableE = mbnl->getOuterBackwardLayers();
     for (MapEI i = reacheableE.begin(); i != reacheableE.end(); i++ ) {
       reachedBarrelLayersMap[(*i).first].insert(*bli);
     }
     reacheableE = mbnl->getOuterForwardLayers();
     for (MapEI i = reacheableE.begin(); i != reacheableE.end(); i++ ) {
       reachedBarrelLayersMap[(*i).first].insert(*bli);
     }
     MapE compatibleE = mbnl->getAllOuterBackwardLayers();
     for (MapEI i = compatibleE.begin(); i != compatibleE.end(); i++ ) {
       compatibleBarrelLayersMap[(*i).first].insert(*bli);
     }
     compatibleE = mbnl->getAllOuterForwardLayers();
     for (MapEI i = compatibleE.begin(); i != compatibleE.end(); i++ ) {
       compatibleBarrelLayersMap[(*i).first].insert(*bli);
     }
 
   }
 
   // collect all reacheable layer starting from a backward layer
   for ( MapEI eli  = theBackwardLayers.begin(); 
               eli != theBackwardLayers.end(); eli++ ) {
     MapE reacheableE =
       dynamic_cast<MuonForwardNavigableLayer*>(theAllNavigableLayer[((*eli).first)->seqNum()])->getOuterEndcapLayers();
     for (MapEI i = reacheableE.begin(); i != reacheableE.end(); i++ ) {
       reachedForwardLayersMap[(*i).first].insert(*eli);
     }
   // collect all compatible layer starting from a backward layer
     MapE compatibleE =
       dynamic_cast<MuonForwardNavigableLayer*>(theAllNavigableLayer[((*eli).first)->seqNum()])->getAllOuterEndcapLayers();
     for (MapEI i = compatibleE.begin(); i != compatibleE.end(); i++ ) {
       compatibleForwardLayersMap[(*i).first].insert(*eli);
     }
   }
 
   for ( MapEI eli  = theForwardLayers.begin(); 
               eli != theForwardLayers.end(); eli++ ) {
   // collect all reacheable layer starting from a forward layer
     MapE reacheableE =
       dynamic_cast<MuonForwardNavigableLayer*>(theAllNavigableLayer[((*eli).first)->seqNum()])->getOuterEndcapLayers();
     for (MapEI i = reacheableE.begin(); i != reacheableE.end(); i++ ) {
       reachedForwardLayersMap[(*i).first].insert(*eli);
     }
   // collect all compatible layer starting from a forward layer
     MapE compatibleE =
       dynamic_cast<MuonForwardNavigableLayer*>(theAllNavigableLayer[((*eli).first)->seqNum()])->getAllOuterEndcapLayers();
     for (MapEI i = compatibleE.begin(); i != compatibleE.end(); i++ ) {
       compatibleForwardLayersMap[(*i).first].insert(*eli);
     }
   }
 
   // now set inverse link for barrel layers
   for ( MapBI bli  = theBarrelLayers.begin(); 
               bli != theBarrelLayers.end(); bli++ ) {
     MuonBarrelNavigableLayer* mbnl =
       dynamic_cast<MuonBarrelNavigableLayer*>(theAllNavigableLayer[((*bli).first)->seqNum()]);
     mbnl->setInwardLinks(reachedBarrelLayersMap[(*bli).first]);
     mbnl->setInwardCompatibleLinks(compatibleBarrelLayersMap[(*bli).first]);
 
   }
   //BACKWARD
   for ( MapEI eli  = theBackwardLayers.begin(); 
               eli != theBackwardLayers.end(); eli++ ) {
     MuonForwardNavigableLayer* mfnl =      
       dynamic_cast<MuonForwardNavigableLayer*>(theAllNavigableLayer[((*eli).first)->seqNum()]);
     // for backward next layers
     mfnl->setInwardLinks(reachedBarrelLayersMap[(*eli).first],
                          reachedForwardLayersMap[(*eli).first]);
   // for backward compatible layers
     mfnl->setInwardCompatibleLinks(compatibleBarrelLayersMap[(*eli).first],
                          compatibleForwardLayersMap[(*eli).first]);
   }
   //FORWARD
   for ( MapEI eli  = theForwardLayers.begin(); 
               eli != theForwardLayers.end(); eli++ ) {
     MuonForwardNavigableLayer* mfnl = 
       dynamic_cast<MuonForwardNavigableLayer*>(theAllNavigableLayer[((*eli).first)->seqNum()]);
     // and for forward next layers
     mfnl->setInwardLinks(reachedBarrelLayersMap[(*eli).first],
                          reachedForwardLayersMap[(*eli).first]);
   // and for forward compatible layers
     mfnl->setInwardCompatibleLinks(compatibleBarrelLayersMap[(*eli).first],
                          compatibleForwardLayersMap[(*eli).first]);
   }
 
 }

void MuonNavigationSchool::linkBarrelLayers ( )

private

link barrel layers

linking barrel layers outwards

Definition at line 160 of file MuonNavigationSchool.cc.

References theBackwardLayers, theBarrelLayers, theBarrelNLC, and theForwardLayers.

Referenced by MuonNavigationSchool().

                                             {
 
   for (MapBI bl  = theBarrelLayers.begin();
              bl != theBarrelLayers.end(); bl++) {
 
     MuonEtaRange range = (*bl).second;
 
     // first add next barrel layer
     MapBI plusOne(bl);
     plusOne++;
     MapB outerBarrel;
     MapB allOuterBarrel;
     if ( plusOne != theBarrelLayers.end() ) { outerBarrel.insert(*plusOne);}
     // add all outer barrel layers
     for ( MapBI iMBI = plusOne; iMBI!= theBarrelLayers.end(); iMBI++){
       allOuterBarrel.insert(*iMBI);
     }
     // then add all compatible backward layers with an eta criteria
     MapE allOuterBackward;
     for (MapEI el  = theBackwardLayers.begin();
                el != theBackwardLayers.end(); el++) {
       if ( (*el).second.isCompatible(range) ) {
         allOuterBackward.insert(*el);
       }
     }
     //add the backward next layer with an eta criteria
     MapE outerBackward;
     for (MapEI el  = theBackwardLayers.begin();
                el != theBackwardLayers.end(); el++) {
       if ( (*el).second.isCompatible(range) ) {
         outerBackward.insert(*el);
         break;
       }
     }
 
     // then add all compatible forward layers with an eta criteria
     MapE allOuterForward;
     for (MapEI el  = theForwardLayers.begin();
                el != theForwardLayers.end(); el++) {
       if ( (*el).second.isCompatible(range) ) {
         allOuterForward.insert(*el);
       }
     }
 
     // then add forward next layer with an eta criteria
     MapE outerForward;
     for (MapEI el  = theForwardLayers.begin();
                el != theForwardLayers.end(); el++) {
       if ( (*el).second.isCompatible(range) ) {
         outerForward.insert(*el);
         break;
       }
     }
 
     theBarrelNLC.push_back(new MuonBarrelNavigableLayer(
                        (*bl).first,outerBarrel, outerBackward, outerForward,
                        allOuterBarrel,allOuterBackward,allOuterForward));
 
   }
 
 }

void MuonNavigationSchool::linkEndcapLayers	(	const MapE &	layers,
		std::vector< MuonForwardNavigableLayer * > &	result
	)

private

link endcap layers

linking forward/backward layers outwards

Definition at line 222 of file MuonNavigationSchool.cc.

References MuonEtaRange::add(), MuonEtaRange::isCompatible(), MuonEtaRange::isInside(), prof2calltree::l, bookConverter::max, GetRecoTauVFromDQM_MC_cff::next, and MuonEtaRange::subtract().

Referenced by MuonNavigationSchool().

                                                                                         {
 
   for (MapEI el = layers.begin(); el != layers.end(); el++) {
 
     MuonEtaRange range = (*el).second;
     // first add next endcap layer (if compatible)
     MapEI plusOne(el); 
     plusOne++;
     MapE outerLayers;
     if ( plusOne != layers.end() && (*plusOne).second.isCompatible(range) ) {
         outerLayers.insert(*plusOne);
       if ( !range.isInside((*plusOne).second) ) {
         // then look if the next layer has a wider eta range, if so add it
         MapEI tmpel(plusOne);
         tmpel++;
         MuonEtaRange max((*plusOne).second);
         for ( MapEI l = tmpel; l != layers.end(); l++ ) {
           MuonEtaRange next = (*l).second;
           if ( next.isCompatible(max) && !range.isInside(next) &&
                !next.isInside(max) && next.subtract(max).isInside(range) ) {
             max = max.add(next);
             outerLayers.insert(*l);
           }
         }
       }
     }
 
     MapE allOuterLayers;
     for (MapEI iMEI = plusOne; iMEI!=layers.end(); iMEI++){
       if ((*iMEI).second.isCompatible(range)) allOuterLayers.insert(*iMEI);
     }
     
     result.push_back(new MuonForwardNavigableLayer(
                    (*el).first,outerLayers, allOuterLayers));
   }
 
 }

MuonNavigationSchool::StateType MuonNavigationSchool::navigableLayers ( )

overridevirtual

return navigable layers, from base class

return all Navigable layers

Implements NavigationSchool.

Definition at line 90 of file MuonNavigationSchool.cc.

References cuy::ib, query::result, theBackwardNLC, theBarrelNLC, and theForwardNLC.

Referenced by createInverseLinks().

                                       {
 
   StateType result;
   
   vector<MuonBarrelNavigableLayer*>::const_iterator ib;
   vector<MuonForwardNavigableLayer*>::const_iterator ie;
 
   for ( ib = theBarrelNLC.begin(); ib != theBarrelNLC.end(); ib++ ) {
     result.push_back(*ib);
   }
 
   for ( ie = theForwardNLC.begin(); ie != theForwardNLC.end(); ie++ ) {
     result.push_back(*ie);
   }
 
   for ( ie = theBackwardNLC.begin(); ie != theBackwardNLC.end(); ie++ ) {
     result.push_back(*ie);
   }
   
   return result;
 
 }