PixelInactiveAreaFinder.cc File Reference

#include "PixelInactiveAreaFinder.h"
#include "FWCore/Utilities/interface/VecArray.h"
#include "FWCore/Utilities/interface/transform.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelQuality.h"
#include "DataFormats/Math/interface/deltaPhi.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "TrackingTools/TransientTrackingRecHit/interface/SeedingLayerSetsLooper.h"
#include <fstream>
#include <queue>
#include <algorithm>

Go to the source code of this file.

Functions
	badPixelFEDChannelsTokens_ (edm::vector_transform(iConfig.getParameter< std::vector< edm::InputTag >>("badPixelFEDChannelCollectionLabels"), [&](const auto &tag){return iC.consumes< PixelFEDChannelCollection >(tag);}))
std::ostream &	operator<< (std::ostream &os, SeedingLayerSetsBuilder::SeedingLayerId layer)
	pixelQualityToken_ (iC.esConsumes())
	trackerGeometryToken_ (iC.esConsumes())
	trackerTopologyToken_ (iC.esConsumes())

Function Documentation

badPixelFEDChannelsTokens_

(

edm::

vector_transformiConfig.getParameter< std::vector< edm::InputTag >>("badPixelFEDChannelCollectionLabels"),[&](const auto &tag){return iC.consumes< PixelFEDChannelCollection >(tag);}

)

std::ostream& operator<<	(	std::ostream &	os,
		SeedingLayerSetsBuilder::SeedingLayerId	layer
	)

Definition at line 16 of file PixelInactiveAreaFinder.cc.

References phase1PixelTopology::layer, NegEndcap, GeomDetEnumerators::PixelBarrel, and PosEndcap.

                                                                                    {
  if (std::get<0>(layer) == GeomDetEnumerators::PixelBarrel) {
    os << "BPix";
  } else {
    os << "FPix";
  }
  os << std::get<2>(layer);
  if (std::get<1>(layer) == TrackerDetSide::PosEndcap) {
    os << "_pos";
  } else if (std::get<1>(layer) == TrackerDetSide::NegEndcap) {
    os << "_neg";
  }
  return os;
}

pixelQualityToken_

(

iC.

esConsumes()

)

Definition at line 475 of file PixelInactiveAreaFinder.cc.

References cms::cuda::assert(), HLT_FULL_cff::distance, Exception, spr::find(), newFWLiteAna::found, mps_fire::i, phase1PixelTopology::layer, LogDebug, LogTrace, runTheMatrix::ret, and HLT_FULL_cff::seedingLayers.

Referenced by MeasurementTrackerESProducer::MeasurementTrackerESProducer(), and MeasurementTrackerESProducer::produce().

                                          {
#ifdef EDM_ML_DEBUG
  for (const auto& layer : seedingLayers) {
    LogTrace("PixelInactiveAreaFinder") << "Input layer subdet " << std::get<0>(layer) << " side "
                                        << static_cast<unsigned int>(std::get<1>(layer)) << " layer "
                                        << std::get<2>(layer);
  }
#endif

  auto findOrAdd = [&](SeedingLayerId layer) -> unsigned short {
    auto found = std::find(inactiveLayers_.cbegin(), inactiveLayers_.cend(), layer);
    if (found == inactiveLayers_.cend()) {
      auto ret = inactiveLayers_.size();
      inactiveLayers_.push_back(layer);
      return ret;
    }
    return std::distance(inactiveLayers_.cbegin(), found);
  };

  // mapping from active layer pairs to inactive layer pairs
  const auto activeToInactiveMap = createActiveToInactiveMap();

  // convert input layer pairs (that are for active layers) to layer
  // pairs to look for inactive areas
  LayerSetIndex i = 0;
  for (const auto& layerSet : seedingLayerSetsLooper.makeRange(seedingLayers)) {
    assert(layerSet.size() == 2);
    auto found = activeToInactiveMap.find(std::make_pair(layerSet[0], layerSet[1]));
    if (found == activeToInactiveMap.end()) {
      throw cms::Exception("Configuration")
          << "Encountered layer pair " << layerSet[0] << "+" << layerSet[1]
          << " not found from the internal 'active layer pairs' to 'inactive layer pairs' mapping; either fix the "
             "input or the mapping (in PixelInactiveAreaFinder.cc)";
    }

    LogTrace("PixelInactiveAreaFinder") << "Input layer set " << layerSet[0] << "+" << layerSet[1];
    for (const auto& inactiveLayerSet : found->second) {
      auto innerInd = findOrAdd(inactiveLayerSet.first);
      auto outerInd = findOrAdd(inactiveLayerSet.second);

      auto found = std::find(
          inactiveLayerSetIndices_.cbegin(), inactiveLayerSetIndices_.cend(), std::make_pair(innerInd, outerInd));
      if (found == inactiveLayerSetIndices_.end()) {
        inactiveLayerSetIndices_.emplace_back(innerInd, outerInd);
        layerSetIndexInactiveToActive_.push_back(std::vector<LayerSetIndex>{i});
      } else {
        layerSetIndexInactiveToActive_.at(std::distance(inactiveLayerSetIndices_.cbegin(), found))
            .push_back(i);  // TODO: move to operator[] once finished
      }

      LogTrace("PixelInactiveAreaFinder")
          << " inactive layer set " << inactiveLayerSet.first << "+" << inactiveLayerSet.second;
    }

    ++i;
  }

#ifdef EDM_ML_DEBUG
  LogDebug("PixelInactiveAreaFinder") << "All inactive layer sets";
  for (const auto& idxPair : inactiveLayerSetIndices_) {
    LogTrace("PixelInactiveAreaFinder") << " " << inactiveLayers_[idxPair.first] << "+"
                                        << inactiveLayers_[idxPair.second];
  }
#endif
}

trackerGeometryToken_

(

iC.

esConsumes()

)

Referenced by SelectedElectronFEDListProducer< TEle, TCand >::produce().

trackerTopologyToken_

(

iC.

esConsumes()

)

Referenced by MeasurementTrackerESProducer::MeasurementTrackerESProducer(), VertexProducer::produce(), and MeasurementTrackerESProducer::produce().