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 40 of file PixelInactiveAreaFinder.h.

Member Typedef Documentation

template<typename T >

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

Definition at line 54 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 42 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 370 of file PixelInactiveAreaFinder.cc.

References edm::VecArray< T, N >::emplace_back(), B2GTnPMonitor_cfi::item, eostools::move(), runTheMatrix::ret, 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(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 399 of file PixelInactiveAreaFinder.cc.

References edm::VecArray< T, N >::emplace_back(), dataset::end, Exception, mps_fire::i, PixelInactiveAreaFinder::layerSetIndexInactiveToActive_, LogDebug, LogTrace, eostools::move(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, sistrip::SpyUtilities::range(), runTheMatrix::ret, contentValuesCheck::ss, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by areasAndLayerSets().

                                                                                                    {
  // 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.

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

private

Definition at line 61 of file PixelInactiveAreaFinder.h.

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

private

Definition at line 62 of file PixelInactiveAreaFinder.h.

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

private

Definition at line 66 of file PixelInactiveAreaFinder.h.