GeometryInterface Class Reference

#include <GeometryInterface.h>

Classes
struct	InterestingQuantities

Public Types
typedef ID	Column
typedef int	ID
typedef double	Value
typedef std::vector< std::pair< Column, Value > >	Values

Public Member Functions
std::vector< InterestingQuantities > const &	allModules ()
Value	binWidth (ID id)
std::pair< Column, Value >	extract (Column const &col, InterestingQuantities const &iq)
Value	extract (ID id, DetId did, edm::Event *ev=nullptr, int16_t col=0, int16_t row=0)
void	extractColumns (std::vector< Column > const &names, InterestingQuantities const &iq, Values &out)
std::string	formatValue (Column, Value)
	GeometryInterface (const edm::ParameterSet &, edm::ConsumesCollector &&, edm::Transition transition=edm::Transition::BeginRun)
ID	intern (std::string const &id)
void	load (edm::EventSetup const &iSetup)
bool	loaded ()
Value	maxValue (ID id)
Value	minValue (ID id)
std::string	pretty (Column col)
std::string	unintern (ID id)

Static Public Attributes
static const Value	UNDEFINED = 999999999.9f

Private Member Functions
void	addExtractor (ID id, std::function< Value(InterestingQuantities const &iq)> func, Value min=UNDEFINED, Value max=UNDEFINED, Value binwidth=1)
void	loadFEDCabling (const SiPixelFedCablingMap *)
void	loadFromSiPixelCoordinates (const TrackerGeometry &, const TrackerTopology &, const SiPixelFedCablingMap &, const edm::ParameterSet &)
void	loadFromTopology (const TrackerGeometry &, const TrackerTopology &, const edm::ParameterSet &)
void	loadModuleLevel (const edm::ParameterSet &iConfig)
void	loadTimebased (const edm::ParameterSet &iConfig)

Private Attributes
std::vector< InterestingQuantities >	all_modules
std::map< ID, Value >	bin_width
std::string	cablingMapLabel_
std::vector< std::function< Value(InterestingQuantities const &iq)> >	extractors
std::map< std::pair< Column, Value >, std::string >	format_value
const edm::ParameterSet	iConfig
std::map< std::string, ID >	ids {std::make_pair(std::string("INVALID"), ID(0))}
bool	is_loaded = false
edm::ESGetToken< SiPixelFedCablingMap, SiPixelFedCablingMapRcd >	labeledSiPixelFedCablingMapToken_
ID	max_id = 0
std::map< ID, Value >	max_value
std::map< ID, Value >	min_value
edm::ESGetToken< SiPixelFedCablingMap, SiPixelFedCablingMapRcd >	siPixelFedCablingMapToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopologyToken_

Detailed Description

Definition at line 34 of file GeometryInterface.h.

Member Typedef Documentation

Column

typedef ID GeometryInterface::Column

Definition at line 40 of file GeometryInterface.h.

ID

typedef int GeometryInterface::ID

Definition at line 37 of file GeometryInterface.h.

Value

typedef double GeometryInterface::Value

Definition at line 41 of file GeometryInterface.h.

Values

typedef std::vector<std::pair<Column, Value> > GeometryInterface::Values

Definition at line 46 of file GeometryInterface.h.

Constructor & Destructor Documentation

GeometryInterface()

GeometryInterface::GeometryInterface	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC,
		edm::Transition	transition = edm::Transition::BeginRun
	)

Definition at line 34 of file GeometryInterface.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

    : iConfig(conf), cablingMapLabel_{conf.getParameter<std::string>("CablingMapLabel")} {
  if (transition == edm::Transition::BeginRun) {
    trackerGeometryToken_ = iC.esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>();
    trackerTopologyToken_ = iC.esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>();
    siPixelFedCablingMapToken_ =
        iC.esConsumes<SiPixelFedCablingMap, SiPixelFedCablingMapRcd, edm::Transition::BeginRun>();
    if (!cablingMapLabel_.empty()) {
      labeledSiPixelFedCablingMapToken_ =
          iC.esConsumes<SiPixelFedCablingMap, SiPixelFedCablingMapRcd, edm::Transition::BeginRun>(
              edm::ESInputTag("", cablingMapLabel_));
    }
  } else {
    trackerGeometryToken_ =
        iC.esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::EndLuminosityBlock>();
    trackerTopologyToken_ = iC.esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::EndLuminosityBlock>();
    siPixelFedCablingMapToken_ =
        iC.esConsumes<SiPixelFedCablingMap, SiPixelFedCablingMapRcd, edm::Transition::EndLuminosityBlock>();
    if (!cablingMapLabel_.empty()) {
      labeledSiPixelFedCablingMapToken_ =
          iC.esConsumes<SiPixelFedCablingMap, SiPixelFedCablingMapRcd, edm::Transition::EndLuminosityBlock>(
              edm::ESInputTag("", cablingMapLabel_));
    }
  }
}

Member Function Documentation

addExtractor()

void GeometryInterface::addExtractor ( ID  id, std::function< Value(InterestingQuantities const &iq)>  func, Value  min = UNDEFINED, Value  max = UNDEFINED, Value  binwidth = 1  )

inlineprivate

Definition at line 169 of file GeometryInterface.h.

References bin_width, extractors, EcalMonitorTask_cff::func, EcalPhiSymFlatTableProducers_cfi::id, WZElectronSkims53X_cff::max, max_value, SiStripPI::min, and min_value.

Referenced by loadFEDCabling(), loadFromSiPixelCoordinates(), loadFromTopology(), loadModuleLevel(), and loadTimebased().

                                        {
    max_value[id] = max;
    min_value[id] = min;
    bin_width[id] = binwidth;
    extractors[id] = func;
  }

allModules()

std::vector<InterestingQuantities> const& GeometryInterface::allModules ( )

inline

Definition at line 99 of file GeometryInterface.h.

References all_modules.

Referenced by SiPixelPhase1GeometryDebug::analyze(), HistogramManager::book(), and HistogramManager::loadFromDQMStore().

{ return all_modules; }

binWidth()

Value GeometryInterface::binWidth ( ID  id )

inline

Definition at line 103 of file GeometryInterface.h.

References bin_width, and EcalPhiSymFlatTableProducers_cfi::id.

Referenced by HistogramManager::book().

{ return bin_width[id]; };

extract() [1/2]

std::pair<Column, Value> GeometryInterface::extract ( Column const &  col, InterestingQuantities const &  iq  )

inline

Definition at line 76 of file GeometryInterface.h.

References cms::cuda::assert(), cuy::col, extractors, EcalPhiSymFlatTableProducers_cfi::id, UNDEFINED, unintern(), and heppy_batch::val.

Referenced by SiPixelPhase1GeometryDebug::analyze(), HistogramManager::book(), extractColumns(), and HistogramManager::fillInternal().

                                                                                     {
    assert(col != 0 || !"Extracting invalid column.");
    ID id = col;
    assert(ID(extractors.size()) > id || !"extractors vector too small!");
    auto& ex = extractors[id];
    if (!ex) {  // we have never heard about this. This is a typo for sure.
      edm::LogError("GeometryInterface") << "Undefined column used: " << unintern(id) << ". Check your spelling.\n";
    } else {
      auto val = ex(iq);
      if (val != UNDEFINED) {
        return std::make_pair(Column{id}, val);
      }
    }
    return std::make_pair(col, UNDEFINED);
  }

extract() [2/2]

Value GeometryInterface::extract ( ID  id, DetId  did, edm::Event *  ev = nullptr, int16_t  col = 0, int16_t  row = 0  )

inline

Definition at line 92 of file GeometryInterface.h.

References cuy::col, makeMEIFBenchmarkPlots::ev, extractors, and EcalPhiSymFlatTableProducers_cfi::id.

                                                                                            {
    InterestingQuantities iq = {ev, did, col, row};
    return extractors[id](iq);
  }

extractColumns()

void GeometryInterface::extractColumns ( std::vector< Column > const &  names, InterestingQuantities const &  iq, Values &  out  )

inline

Definition at line 66 of file GeometryInterface.h.

References cuy::col, extract(), names, MillePedeFileConverter_cfg::out, and heppy_batch::val.

Referenced by HistogramManager::book(), HistogramManager::executeExtend(), HistogramManager::executeGroupBy(), HistogramManager::executePerEventHarvesting(), HistogramManager::fill(), and HistogramManager::loadFromDQMStore().

                                                                                                    {
    out.clear();
    for (Column const& col : names) {
      auto val = extract(col, iq);
      out.push_back(val);
    }
  }

formatValue()

std::string GeometryInterface::formatValue	(	Column	col,
		Value	val
	)

Definition at line 485 of file GeometryInterface.cc.

References cuy::col, format_value, createfilelist::int, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, mergeVDriftHistosByStation::name, pretty(), AlCaHLTBitMon_QueryRunRegistry::string, to_string(), UNDEFINED, heppy_batch::val, and relativeConstraints::value.

Referenced by HistogramManager::makePathName().

                                                              {
  auto it = format_value.find(std::make_pair(col, val));
  if (it != format_value.end())
    return it->second;

  // non-number output names (_pO etc.) are hardwired here.
  std::string name = pretty(col);
  std::string value = "_" + std::to_string(int(val));
  if (val == 0)
    value = "";                     // hide Barrel_0 etc.
  if (name == "PXDisk" && val > 0)  // +/- sign for disk num
    value = "_+" + std::to_string(int(val));
  // pretty (legacy?) names for Shells and HalfCylinders
  std::map<int, std::string> shellname{{11, "_mI"}, {12, "_mO"}, {21, "_pI"}, {22, "_pO"}};
  if (name == "HalfCylinder" || name == "Shell")
    value = shellname[int(val)];
  if (val == UNDEFINED)
    value = "_UNDEFINED";
  return format_value[std::make_pair(col, val)] = name + value;
}

intern()

ID GeometryInterface::intern ( std::string const &  id )

inline

Definition at line 106 of file GeometryInterface.h.

References extractors, EcalPhiSymFlatTableProducers_cfi::id, ids, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and max_id.

Referenced by SiPixelPhase1GeometryDebug::analyze(), loadFEDCabling(), loadFromSiPixelCoordinates(), loadFromTopology(), loadModuleLevel(), loadTimebased(), and SummationSpecification::parse_columns().

                                 {
    auto it = ids.find(id);
    if (it == ids.end()) {
      ids[id] = ++max_id;
      extractors.resize(max_id + 1);
    }
    return ids[id];
  };

load()

void GeometryInterface::load

(

edm::EventSetup const &

iSetup

)

Definition at line 62 of file GeometryInterface.cc.

References nano_mu_local_reco_cff::bool, cablingMapLabel_, MillePedeFileConverter_cfg::e, extractors, edm::EventSetup::getData(), iConfig, ids, is_loaded, labeledSiPixelFedCablingMapToken_, loadFEDCabling(), loadFromSiPixelCoordinates(), loadFromTopology(), loadModuleLevel(), loadTimebased(), CrabHelper::log, max_value, min_value, siPixelFedCablingMapToken_, trackerGeometryToken_, and trackerTopologyToken_.

Referenced by HistogramManager::book(), and HistogramManager::executePerLumiHarvesting().

                                                        {
  const TrackerGeometry& trackerGeometry = iSetup.getData(trackerGeometryToken_);
  const TrackerTopology& trackerTopology = iSetup.getData(trackerTopologyToken_);
  const SiPixelFedCablingMap& siPixelFedCablingMap = iSetup.getData(siPixelFedCablingMapToken_);

  const SiPixelFedCablingMap* labeledSiPixelFedCablingMap = &siPixelFedCablingMap;
  if (!cablingMapLabel_.empty()) {
    labeledSiPixelFedCablingMap = &iSetup.getData(labeledSiPixelFedCablingMapToken_);
  }

  //loadFromAlignment(iSetup, iConfig);
  loadFromTopology(trackerGeometry, trackerTopology, iConfig);
  loadTimebased(iConfig);
  loadModuleLevel(iConfig);
  loadFEDCabling(labeledSiPixelFedCablingMap);
  loadFromSiPixelCoordinates(trackerGeometry, trackerTopology, siPixelFedCablingMap, iConfig);
  edm::LogInfo log("GeometryInterface");
  log << "Known colum names:\n";
  for (const auto& e : ids)
    log << "+++ column: " << e.first << " ok " << bool(extractors[e.second]) << " min " << min_value[e.second]
        << " max " << max_value[e.second] << "\n";
  is_loaded = true;
}

loaded()

bool GeometryInterface::loaded ( )

inline

Definition at line 52 of file GeometryInterface.h.

References is_loaded.

Referenced by HistogramManager::book(), and HistogramManager::executePerLumiHarvesting().

{ return is_loaded; };

loadFEDCabling()

void GeometryInterface::loadFEDCabling ( const SiPixelFedCablingMap *  labeledSiPixelFedCablingMap )

private

Definition at line 453 of file GeometryInterface.cc.

References addExtractor(), GeometryInterface::InterestingQuantities::col, runTauDisplay::gp, intern(), DetId::rawId(), GeometryInterface::InterestingQuantities::row, and GeometryInterface::InterestingQuantities::sourceModule.

Referenced by load().

                                                                                              {
  std::shared_ptr<SiPixelFrameReverter> siPixelFrameReverter =
      // I think passing the bare pointer here is safe, but who knows...
      std::make_shared<SiPixelFrameReverter>(labeledSiPixelFedCablingMap);

  addExtractor(intern("FED"), [siPixelFrameReverter](InterestingQuantities const& iq) {
    if (iq.sourceModule == 0xFFFFFFFF)
      return Value(iq.col);  // hijacked for the raw data plugin
    return Value(siPixelFrameReverter->findFedId(iq.sourceModule.rawId()));
  });

  // TODO: ranges should be set manually below, since booking probably cannot
  // infer them correctly (no ROC-level granularity)
  // PERF: this is slow. Prefer SiPixelCordinates versions here.
  addExtractor(intern("LinkInFed"), [siPixelFrameReverter](InterestingQuantities const& iq) {
    if (iq.sourceModule == 0xFFFFFFFF)
      return Value(iq.row);  // hijacked for the raw data plugin
    sipixelobjects::GlobalPixel gp = {iq.row, iq.col};
    return Value(siPixelFrameReverter->findLinkInFed(iq.sourceModule.rawId(), gp));
  });
  // not sure if this is useful anywhere.
  addExtractor(intern("RocInLink"), [siPixelFrameReverter](InterestingQuantities const& iq) {
    sipixelobjects::GlobalPixel gp = {iq.row, iq.col};
    return Value(siPixelFrameReverter->findRocInLink(iq.sourceModule.rawId(), gp));
  });
  // This might be equivalent to our ROC numbering.
  addExtractor(intern("RocInDet"), [siPixelFrameReverter](InterestingQuantities const& iq) {
    sipixelobjects::GlobalPixel gp = {iq.row, iq.col};
    return Value(siPixelFrameReverter->findRocInDet(iq.sourceModule.rawId(), gp));
  });
}

loadFromSiPixelCoordinates()

void GeometryInterface::loadFromSiPixelCoordinates ( const TrackerGeometry &  trackerGeometry, const TrackerTopology &  trackerTopology, const SiPixelFedCablingMap &  siPixelFedCablingMap, const edm::ParameterSet &  iConfig  )

private

Definition at line 197 of file GeometryInterface.cc.

References addExtractor(), GeometryInterface::InterestingQuantities::col, extractors, edm::ParameterSet::getParameter(), iConfig, recoMuon::in, intern(), SiPixelPI::quadrant(), GeometryInterface::InterestingQuantities::row, GeometryInterface::InterestingQuantities::sourceModule, and UNDEFINED.

Referenced by load().

                                                                                   {
  // TODO: SiPixelCoordinates has a large overlap with theis GeometryInterface
  // in general.
  // Rough convention is to use own code for things that are easy and fast to
  // determine, and use SiPixelCoordinates for complicated things.
  // SiPixelCoordinates uses lookup maps for everything, so it is faster than
  // most other code, but still slow on DQM scales.
  int phase = iConfig.getParameter<int>("upgradePhase");

  // this shared pointer is kept alive by the references in the lambdas that follow.
  // That is a bit less obvious than keeping it as a member but more correct.
  auto coord = std::make_shared<SiPixelCoordinates>(phase);

  // note that we should reeinit for each event. But this probably won't explode
  // thanks to the massive memoization in SiPixelCoordinates which is completely
  // initialized while booking.
  coord->init(&trackerTopology, &trackerGeometry, &siPixelFedCablingMap);

  // SiPixelCoordinates uses a different convention for UNDEFINED:
  auto from_coord = [](double in) { return (in == -9999.0) ? UNDEFINED : Value(in); };

  // Rings are a concept that cannot be derived from bitmasks.
  addExtractor(intern("PXRing"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->ring(iq.sourceModule));
  });

  // Quadrant names.
  auto pxbarrel = extractors[intern("PXBarrel")];
  addExtractor(
      intern("HalfCylinder"),
      [coord, pxbarrel](InterestingQuantities const& iq) {
        if (pxbarrel(iq) != UNDEFINED)
          return UNDEFINED;
        int quadrant = coord->quadrant(iq.sourceModule);
        switch (quadrant) {
          case 1:
            return Value(12);  // mO
          case 2:
            return Value(11);  // mI
          case 3:
            return Value(22);  // pO
          case 4:
            return Value(21);  // pI
          default:
            return UNDEFINED;
        }
      },
      0,
      0  // N/A
  );
  addExtractor(
      intern("Shell"),
      [coord, pxbarrel](InterestingQuantities const& iq) {
        if (pxbarrel(iq) == UNDEFINED)
          return UNDEFINED;
        int quadrant = coord->quadrant(iq.sourceModule);
        switch (quadrant) {
          case 1:
            return Value(12);  // mO
          case 2:
            return Value(11);  // mI
          case 3:
            return Value(22);  // pO
          case 4:
            return Value(21);  // pI
          default:
            return UNDEFINED;
        }
      },
      0,
      0  // N/A
  );

  // Online Numbering.
  addExtractor(intern("SignedLadder"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->signed_ladder(iq.sourceModule()));
  });
  addExtractor(intern("SignedModule"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->signed_module(iq.sourceModule()));
  });
  addExtractor(intern("SignedBlade"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->signed_blade(iq.sourceModule()));
  });

  // Flipped and Outer ladders
  addExtractor(intern("OuterLadder"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->outer(iq.sourceModule()));
  });
  addExtractor(intern("FlippedLadder"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->flipped(iq.sourceModule()));
  });

  // Pixel Map axis.
  // TODO: automatic range and binning for phase0 are incorrect.
  // Should be set manually here.
  addExtractor(
      intern("SignedModuleCoord"),  // BPIX x
      [coord, from_coord](InterestingQuantities const& iq) {
        return from_coord(coord->signed_module_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
      },
      UNDEFINED,
      UNDEFINED,
      1.0 / 8.0);
  addExtractor(
      intern("SignedLadderCoord"),  // BPIX y
      [coord, from_coord](InterestingQuantities const& iq) {
        return from_coord(coord->signed_ladder_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
      },
      UNDEFINED,
      UNDEFINED,
      1.0 / 2.0);
  addExtractor(
      intern("SignedDiskCoord"),  // FPIX x (per-ring)
      [coord, from_coord](InterestingQuantities const& iq) {
        return from_coord(coord->signed_disk_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
      },
      UNDEFINED,
      UNDEFINED,
      1.0 / 8.0);
  addExtractor(
      intern("SignedDiskRingCoord"),  // FPIX x (FPIX-as-one-plot)
      [coord, from_coord](InterestingQuantities const& iq) {
        return from_coord(coord->signed_disk_ring_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
      },
      UNDEFINED,
      UNDEFINED,
      1.0 / 16.0);
  addExtractor(
      intern("SignedBladePanelCoord"),  // FPIX y
      [coord, from_coord, phase](InterestingQuantities const& iq) {
        if (phase == 0) {
          return from_coord(coord->signed_blade_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
        } else if (phase == 1) {
          return from_coord(
              coord->signed_blade_panel_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
        } else {
          return UNDEFINED;  // TODO: phase2
        }
      },
      UNDEFINED,
      UNDEFINED,
      phase == 1 ? 0.25 : 0.2);
  addExtractor(
      intern("SignedShiftedBladePanelCoord"),  // FPIX-as-one y
      [coord, from_coord, phase](InterestingQuantities const& iq) {
        if (phase == 0) {
          return from_coord(coord->signed_blade_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
        } else if (phase == 1) {
          return from_coord(
              coord->signed_shifted_blade_panel_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
        } else {
          return UNDEFINED;  // TODO: phase2
        }
      },
      UNDEFINED,
      UNDEFINED,
      phase == 1 ? 0.25 : 0.1  // half-roc for phase0
  );
  addExtractor(
      intern("SignedBladePanel"),  // per-module FPIX y
      [coord, from_coord](InterestingQuantities const& iq) {
        return from_coord(coord->signed_blade_panel_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
      },
      UNDEFINED,
      UNDEFINED,
      1.0 / 2.0);

  addExtractor(
      intern("SignedBladePanel"),
      [coord, from_coord](InterestingQuantities const& iq) {
        return from_coord(coord->signed_blade_panel_coord(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
      },
      UNDEFINED,
      UNDEFINED,
      1.0 / 2.0);

  // more readout-related things.
  addExtractor(intern("ROC"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->roc(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
  });
  addExtractor(intern("Sector"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->sector(iq.sourceModule()));
  });
  addExtractor(intern("Channel"), [coord, from_coord](InterestingQuantities const& iq) {
    return from_coord(coord->channel(iq.sourceModule(), std::make_pair(int(iq.row), int(iq.col))));
  });
}

loadFromTopology()

void GeometryInterface::loadFromTopology ( const TrackerGeometry &  trackerGeometry, const TrackerTopology &  trackerTopology, const edm::ParameterSet &  iConfig  )

private

Definition at line 86 of file GeometryInterface.cc.

References addExtractor(), all_modules, cms::cuda::assert(), ALCARECOPPSCalTrackBasedSel_cff::detid, TrackerGeometry::detIds(), MillePedeFileConverter_cfg::e, makeMuonMisalignmentScenario::endcap, extractors, format_value, edm::ParameterSet::getParameter(), iConfig, EcalPhiSymFlatTableProducers_cfi::id, TrackerGeometry::idToDetUnit(), intern(), SiPixelOfflineDQM_client_cff::isUpgrade, mod(), mergeVDriftHistosByStation::name, PixelTopology::ncolumns(), PixelTopology::nrows(), TrackerTopology::PBLadder, TrackerTopology::PBLayer, TrackerTopology::PBModule, TrackerTopology::PFBlade, TrackerTopology::PFDisk, TrackerTopology::PFModule, TrackerTopology::PFPanel, TrackerTopology::PFSide, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, nano_mu_digi_cff::rawId, GeometryInterface::InterestingQuantities::sourceModule, PixelGeomDetUnit::specificTopology(), AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::tt, and UNDEFINED.

Referenced by load().

                                                                         {
  std::vector<ID> geomquantities;

  struct TTField {
    const TrackerTopology* tt;
    TrackerTopology::DetIdFields field;
    Value operator()(InterestingQuantities const& iq) {
      if (tt->hasField(iq.sourceModule, field))
        return tt->getField(iq.sourceModule, field);
      else
        return UNDEFINED;
    };
  };

  const TrackerTopology* tt = &trackerTopology;

  std::vector<std::pair<std::string, TTField>> namedPartitions{
      {"PXEndcap", {tt, TrackerTopology::PFSide}},

      {"PXLayer", {tt, TrackerTopology::PBLayer}},
      {"PXLadder", {tt, TrackerTopology::PBLadder}},
      {"PXBModule", {tt, TrackerTopology::PBModule}},

      {"PXBlade", {tt, TrackerTopology::PFBlade}},
      {"PXDisk", {tt, TrackerTopology::PFDisk}},
      {"PXPanel", {tt, TrackerTopology::PFPanel}},
      {"PXFModule", {tt, TrackerTopology::PFModule}},
  };

  for (auto& e : namedPartitions) {
    geomquantities.push_back(intern(e.first));
    addExtractor(intern(e.first), e.second, UNDEFINED, UNDEFINED);
  }

  auto pxbarrel = [](InterestingQuantities const& iq) {
    return iq.sourceModule.subdetId() == PixelSubdetector::PixelBarrel ? 0 : UNDEFINED;
  };
  auto pxforward = [](InterestingQuantities const& iq) {
    return iq.sourceModule.subdetId() == PixelSubdetector::PixelEndcap ? 0 : UNDEFINED;
  };
  auto pxall = [](InterestingQuantities const& iq) { return 0; };
  addExtractor(intern("PXBarrel"), pxbarrel, 0, 0);
  addExtractor(intern("PXForward"), pxforward, 0, 0);
  addExtractor(intern("PXAll"), pxall, 0, 0);

  // Redefine the disk numbering to use the sign
  auto pxendcap = extractors[intern("PXEndcap")];
  auto diskid = intern("PXDisk");
  auto pxdisk = extractors[diskid];
  extractors[diskid] = [pxdisk, pxendcap](InterestingQuantities const& iq) {
    auto disk = pxdisk(iq);
    if (disk == UNDEFINED)
      return UNDEFINED;
    auto endcap = pxendcap(iq);
    return endcap == 1 ? -disk : disk;
  };

  // DetId and module names
  auto detid = [](InterestingQuantities const& iq) {
    uint32_t id = iq.sourceModule.rawId();
    return Value(id);
  };
  addExtractor(intern("DetId"), detid, 0, 0  // No sane value possible here.
  );
  // these are just aliases with special handling in formatting
  // the names are created with PixelBarrelName et. al. later
  addExtractor(intern("PXModuleName"), detid, 0, 0);

  int phase = iConfig.getParameter<int>("upgradePhase");
  bool isUpgrade = phase == 1;

  // Now traverse the detector and collect whatever we need.
  auto detids = trackerGeometry.detIds();
  for (DetId id : detids) {
    if (id.subdetId() != PixelSubdetector::PixelBarrel && id.subdetId() != PixelSubdetector::PixelEndcap)
      continue;
    auto iq = InterestingQuantities{nullptr, id, 0, 0};

    // prepare pretty names
    std::string name = "";
    if (id.subdetId() == PixelSubdetector::PixelBarrel) {  // Barrel
      PixelBarrelName mod(id, tt, isUpgrade);
      name = mod.name();
    } else {  // assume Endcap
      PixelEndcapName mod(id, tt, isUpgrade);
      name = mod.name();
    }
    format_value[std::make_pair(intern("PXModuleName"), Value(id.rawId()))] = name;

    // we record each module 4 times, one for each corner, so we also get ROCs
    // in booking (at least for the ranges)
    const PixelGeomDetUnit* detUnit = dynamic_cast<const PixelGeomDetUnit*>(trackerGeometry.idToDetUnit(id));
    assert(detUnit);
    const PixelTopology* topo = &detUnit->specificTopology();
    iq.row = 0;
    iq.col = 0;
    all_modules.push_back(iq);
    iq.row = topo->nrows() - 1;
    iq.col = 0;
    all_modules.push_back(iq);
    iq.row = 0;
    iq.col = topo->ncolumns() - 1;
    all_modules.push_back(iq);
    iq.row = topo->nrows() - 1;
    iq.col = topo->ncolumns() - 1;
    all_modules.push_back(iq);
  }
}

loadModuleLevel()

void GeometryInterface::loadModuleLevel ( const edm::ParameterSet &  iConfig )

private

Definition at line 439 of file GeometryInterface.cc.

References addExtractor(), GeometryInterface::InterestingQuantities::col, edm::ParameterSet::getParameter(), iConfig, intern(), and GeometryInterface::InterestingQuantities::row.

Referenced by load().

                                                                      {
  // stuff that is within modules. Might require some phase0/phase1/strip switching later
  addExtractor(
      intern("row"),
      [](InterestingQuantities const& iq) { return Value(iq.row); },
      0,
      iConfig.getParameter<int>("module_rows") - 1);
  addExtractor(
      intern("col"),
      [](InterestingQuantities const& iq) { return Value(iq.col); },
      0,
      iConfig.getParameter<int>("module_cols") - 1);
}

loadTimebased()

void GeometryInterface::loadTimebased ( const edm::ParameterSet &  iConfig )

private

Definition at line 388 of file GeometryInterface.cc.

References addExtractor(), edm::EventBase::bunchCrossing(), edm::ParameterSet::getParameter(), iConfig, intern(), HistogramManager_cfi::lumiblock, edm::EventBase::luminosityBlock(), HistogramManager_cfi::n_onlineblocks, HistogramManager_cfi::onlineblock, GeometryInterface::InterestingQuantities::sourceEvent, and UNDEFINED.

Referenced by load().

                                                                    {
  // extractors for quantities that are roughly time-based. We cannot book plots based on these; they have to
  // be grouped away in step1.
  addExtractor(
      intern("Lumisection"),
      [](InterestingQuantities const& iq) {
        if (!iq.sourceEvent)
          return UNDEFINED;
        return Value(iq.sourceEvent->luminosityBlock());
      },
      1,
      iConfig.getParameter<int>("max_lumisection"));

  int onlineblock = iConfig.getParameter<int>("onlineblock");
  int n_onlineblocks = iConfig.getParameter<int>("n_onlineblocks");
  addExtractor(
      intern("OnlineBlock"),
      [onlineblock](InterestingQuantities const& iq) {
        if (!iq.sourceEvent)
          return UNDEFINED;
        return Value(onlineblock + iq.sourceEvent->luminosityBlock() / onlineblock);
      },
      // note: this range is not visible anywhere (if the RenderPlugin does its job),
      // but the strange range allows the RenderPlugin to know the block size.
      onlineblock,
      onlineblock + n_onlineblocks - 1);

  int lumiblock = iConfig.getParameter<int>("lumiblock");
  addExtractor(
      intern("LumiBlock"),
      [lumiblock](InterestingQuantities const& iq) {
        if (!iq.sourceEvent)
          return UNDEFINED;
        // The '-1' is for making 1-10 the same block rather than 0-9
        // The '+0.5' makes the block span an integer range rather n.5-m.5
        return Value(((iq.sourceEvent->luminosityBlock() - 1) / lumiblock) + 0.5);
      },
      -0.5,
      iConfig.getParameter<int>("max_lumisection") / lumiblock);

  addExtractor(
      intern("BX"),
      [](InterestingQuantities const& iq) {
        if (!iq.sourceEvent)
          return UNDEFINED;
        return Value(iq.sourceEvent->bunchCrossing());
      },
      1,
      iConfig.getParameter<int>("max_bunchcrossing"));
}

maxValue()

Value GeometryInterface::maxValue ( ID  id )

inline

Definition at line 101 of file GeometryInterface.h.

References EcalPhiSymFlatTableProducers_cfi::id, and max_value.

Referenced by HistogramManager::book().

{ return max_value[id]; };

minValue()

Value GeometryInterface::minValue ( ID  id )

inline

Definition at line 102 of file GeometryInterface.h.

References EcalPhiSymFlatTableProducers_cfi::id, and min_value.

Referenced by HistogramManager::book().

{ return min_value[id]; };

pretty()

std::string GeometryInterface::pretty ( Column  col )

inline

Definition at line 125 of file GeometryInterface.h.

References cuy::col, and unintern().

Referenced by HistogramManager::book(), HistogramManager::executeExtend(), formatValue(), and HistogramManager::makePathName().

{ return unintern(col); }

unintern()

std::string GeometryInterface::unintern ( ID  id )

inline

Definition at line 118 of file GeometryInterface.h.

References MillePedeFileConverter_cfg::e, and ids.

Referenced by extract(), and pretty().

                            {
    for (auto& e : ids)
      if (e.second == id)
        return e.first;
    return "INVALID";
  }

Member Data Documentation

all_modules

std::vector<InterestingQuantities> GeometryInterface::all_modules

private

Definition at line 180 of file GeometryInterface.h.

Referenced by allModules(), and loadFromTopology().

bin_width

std::map<ID, Value> GeometryInterface::bin_width

private

Definition at line 163 of file GeometryInterface.h.

Referenced by addExtractor(), and binWidth().

cablingMapLabel_

std::string GeometryInterface::cablingMapLabel_

private

Definition at line 152 of file GeometryInterface.h.

Referenced by load().

extractors

std::vector<std::function<Value(InterestingQuantities const& iq)> > GeometryInterface::extractors

private

Definition at line 158 of file GeometryInterface.h.

Referenced by addExtractor(), extract(), intern(), load(), loadFromSiPixelCoordinates(), and loadFromTopology().

format_value

std::map<std::pair<Column, Value>, std::string> GeometryInterface::format_value

private

Definition at line 167 of file GeometryInterface.h.

Referenced by formatValue(), and loadFromTopology().

iConfig

const edm::ParameterSet GeometryInterface::iConfig

private

Definition at line 139 of file GeometryInterface.h.

Referenced by load(), loadFromSiPixelCoordinates(), loadFromTopology(), loadModuleLevel(), and loadTimebased().

ids

std::map<std::string, ID> GeometryInterface::ids {std::make_pair(std::string("INVALID"), ID(0))}

private

Definition at line 183 of file GeometryInterface.h.

Referenced by intern(), load(), and unintern().

is_loaded

bool GeometryInterface::is_loaded = false

private

Definition at line 154 of file GeometryInterface.h.

Referenced by load(), and loaded().

labeledSiPixelFedCablingMapToken_

edm::ESGetToken<SiPixelFedCablingMap, SiPixelFedCablingMapRcd> GeometryInterface::labeledSiPixelFedCablingMapToken_

private

Definition at line 151 of file GeometryInterface.h.

Referenced by load().

max_id

ID GeometryInterface::max_id = 0

private

Definition at line 184 of file GeometryInterface.h.

Referenced by intern().

max_value

std::map<ID, Value> GeometryInterface::max_value

private

Definition at line 161 of file GeometryInterface.h.

Referenced by addExtractor(), load(), and maxValue().

min_value

std::map<ID, Value> GeometryInterface::min_value

private

Definition at line 162 of file GeometryInterface.h.

Referenced by addExtractor(), load(), and minValue().

siPixelFedCablingMapToken_

edm::ESGetToken<SiPixelFedCablingMap, SiPixelFedCablingMapRcd> GeometryInterface::siPixelFedCablingMapToken_

private

Definition at line 150 of file GeometryInterface.h.

Referenced by load().

trackerGeometryToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> GeometryInterface::trackerGeometryToken_

private

Definition at line 141 of file GeometryInterface.h.

Referenced by load().

trackerTopologyToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> GeometryInterface::trackerTopologyToken_

private

Definition at line 142 of file GeometryInterface.h.

Referenced by load().

UNDEFINED

const GeometryInterface::Value GeometryInterface::UNDEFINED = 999999999.9f

static

Definition at line 42 of file GeometryInterface.h.

Referenced by SiPixelPhase1GeometryDebug::analyze(), HistogramManager::book(), extract(), formatValue(), loadFromSiPixelCoordinates(), loadFromTopology(), and loadTimebased().