PixelInactiveAreaFinder::InactiveAreas Class Reference

#include <PixelInactiveAreaFinder.h>

Public Types
template<typename T >
using	VecArray2 = edm::VecArray< T, 2 >

Public Member Functions
std::vector< std::pair< VecArray2< Area >, std::vector< LayerSetIndex > > >	areasAndLayerSets (const GlobalPoint &point, float zwidth) const
	InactiveAreas (const std::vector< SeedingLayerId > *inactiveLayers, std::vector< DetGroupSpanContainer > &&inactiveSpans, const std::vector< std::pair< unsigned short, unsigned short > > *inactiveLayerPairIndices, const std::vector< std::vector< LayerSetIndex > > *layerSetIndexInactiveToActive)
std::vector< std::pair< VecArray2< DetGroupSpan >, std::vector< LayerSetIndex > > >	spansAndLayerSets (const GlobalPoint &point, float zwidth) const

Private Attributes
const std::vector< std::pair< unsigned short, unsigned short > > *	inactiveLayerPairIndices_
const std::vector< SeedingLayerId > *	inactiveLayers_
std::vector< DetGroupSpanContainer >	inactiveSpans_
const std::vector< std::vector< LayerSetIndex > > *	layerSetIndexInactiveToActive_

Detailed Description

Definition at line 44 of file PixelInactiveAreaFinder.h.

Member Typedef Documentation

template<typename T >

using PixelInactiveAreaFinder::InactiveAreas::VecArray2 = edm::VecArray<T, 2>

Definition at line 57 of file PixelInactiveAreaFinder.h.

Constructor & Destructor Documentation

PixelInactiveAreaFinder::InactiveAreas::InactiveAreas ( const std::vector< SeedingLayerId > *  inactiveLayers, std::vector< DetGroupSpanContainer > &&  inactiveSpans, const std::vector< std::pair< unsigned short, unsigned short > > *  inactiveLayerPairIndices, const std::vector< std::vector< LayerSetIndex > > *  layerSetIndexInactiveToActive  )

inline

Definition at line 46 of file PixelInactiveAreaFinder.h.

                                                                                            :
      inactiveLayers_(inactiveLayers),
      inactiveSpans_(std::move(inactiveSpans)),
      inactiveLayerPairIndices_(inactiveLayerPairIndices),
      layerSetIndexInactiveToActive_(layerSetIndexInactiveToActive)
    {}

Member Function Documentation

std::vector< std::pair< edm::VecArray< PixelInactiveAreaFinder::Area, 2 >, std::vector< PixelInactiveAreaFinder::LayerSetIndex > > > PixelInactiveAreaFinder::InactiveAreas::areasAndLayerSets	(	const GlobalPoint &	point,
		float	zwidth
	)		const

Definition at line 352 of file PixelInactiveAreaFinder.cc.

References edm::VecArray< T, N >::emplace_back(), eostools::move(), and spansAndLayerSets().

Referenced by PixelInactiveAreaTrackingRegionsSeedingLayersProducer::produce().

                                                                                                    {
  auto spansLayerSets = spansAndLayerSets(point, zwidth);

  // TODO: try to remove this conversion...
  std::vector<std::pair<VecArray2<Area>, std::vector<LayerSetIndex> > > ret;
  for(auto& item: spansLayerSets) {
    auto& innerSpan = item.first[0];
    auto& outerSpan = item.first[1];
    VecArray2<Area> areas;
    areas.emplace_back(innerSpan.rSpan.first, innerSpan.rSpan.second,
                       innerSpan.phiSpan.first, innerSpan.phiSpan.second,
                       innerSpan.zSpan.first, innerSpan.zSpan.second);
    areas.emplace_back(outerSpan.rSpan.first, outerSpan.rSpan.second,
                       outerSpan.phiSpan.first, outerSpan.phiSpan.second,
                       outerSpan.zSpan.first, outerSpan.zSpan.second);
    ret.emplace_back(std::move(areas), std::move(item.second));
  }

  return ret;
}

std::vector< std::pair< edm::VecArray< PixelInactiveAreaFinder::DetGroupSpan, 2 >, std::vector< PixelInactiveAreaFinder::LayerSetIndex > > > PixelInactiveAreaFinder::InactiveAreas::spansAndLayerSets	(	const GlobalPoint &	point,
		float	zwidth
	)		const

Definition at line 374 of file PixelInactiveAreaFinder.cc.

References edm::VecArray< T, N >::emplace_back(), end, Exception, mps_fire::i, inactiveLayerPairIndices_, inactiveSpans_, layerSetIndexInactiveToActive_, LogDebug, LogTrace, eostools::move(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by areasAndLayerSets(), and operator<<().

                                                                                                    {
  // TODO: in the future use 2D-r for the origin for the phi overlap check
  const float zmin = point.z()-zwidth;
  const float zmax = point.z()+zwidth;

  std::vector<std::pair<VecArray2<DetGroupSpan>, std::vector<LayerSetIndex> > > ret;

  LogDebug("PixelInactiveAreaFinder") << "Origin at " << point.x() << "," << point.y() << "," << point.z() << " z half width " << zwidth;

  for(LayerSetIndex i=0, end=inactiveLayerPairIndices_->size(); i<end; ++i) {
    const auto& layerIdxPair = (*inactiveLayerPairIndices_)[i];
    const auto& innerSpans = inactiveSpans_[layerIdxPair.first];
    const auto& outerSpans = inactiveSpans_[layerIdxPair.second];

    for(const auto& innerSpan: innerSpans) {
      for(const auto& outerSpan: outerSpans) {

        if(phiRangesOverlap(innerSpan.phiSpan, outerSpan.phiSpan)) {
          std::pair<float,float> range(0,0);

          bool zOverlap = false;
          const auto innerDet = std::get<0>((*inactiveLayers_)[layerIdxPair.first]);
          const auto outerDet = std::get<0>((*inactiveLayers_)[layerIdxPair.second]);
          if(innerDet == GeomDetEnumerators::PixelBarrel) {
            if(outerDet == GeomDetEnumerators::PixelBarrel)
              zOverlap = getZAxisOverlapRangeBarrel(innerSpan, outerSpan, range);
            else
              zOverlap = getZAxisOverlapRangeBarrelEndcap(innerSpan, outerSpan, range);
          }
          else {
            if(outerDet == GeomDetEnumerators::PixelEndcap)
              zOverlap = getZAxisOverlapRangeEndcap(innerSpan, outerSpan, range);
            else
              throw cms::Exception("LogicError") << "Forward->barrel transition is not supported";
          }

          if(zOverlap && zmin <= range.second && range.first <= zmax) {
#ifdef EDM_ML_DEBUG
            Stream ss;
            for(auto ind: (*layerSetIndexInactiveToActive_)[i]) {
              ss << ind << ",";
            }
            ss << "\n  ";
            detGroupSpanInfo(innerSpan, ss);
            ss << "\n  ";
            detGroupSpanInfo(outerSpan, ss);
            LogTrace("PixelInactiveAreaFinder") << " adding areas for active layer sets " << ss.str();
#endif
            VecArray2<DetGroupSpan> vec;
            vec.emplace_back(innerSpan);
            vec.emplace_back(outerSpan);
            ret.emplace_back(std::move(vec), (*layerSetIndexInactiveToActive_)[i]);
          }
        }
      }
    }
  }

  return ret;
}

Member Data Documentation

const std::vector<std::pair<unsigned short, unsigned short> >* PixelInactiveAreaFinder::InactiveAreas::inactiveLayerPairIndices_

private

Definition at line 64 of file PixelInactiveAreaFinder.h.

Referenced by spansAndLayerSets().

const std::vector<SeedingLayerId>* PixelInactiveAreaFinder::InactiveAreas::inactiveLayers_

private

Definition at line 62 of file PixelInactiveAreaFinder.h.

std::vector<DetGroupSpanContainer> PixelInactiveAreaFinder::InactiveAreas::inactiveSpans_

private

Definition at line 63 of file PixelInactiveAreaFinder.h.

Referenced by spansAndLayerSets().

const std::vector<std::vector<LayerSetIndex> >* PixelInactiveAreaFinder::InactiveAreas::layerSetIndexInactiveToActive_

private

Definition at line 65 of file PixelInactiveAreaFinder.h.

Referenced by spansAndLayerSets().