Reimplemented from edm::EDAnalyzer.
Definition at line 215 of file OccupancyPlots.cc.
219 LogDebug(
"IsThereTest") <<
"Test of TrackerGeometry::isThere";
241 for (TrackingGeometry::DetIdContainer::const_iterator det = detunits.begin(); det != detunits.end(); ++det) {
248 GlobalPoint zpos = trkgeo.idToDet(*det)->toGlobal(locz);
249 GlobalPoint xpos = trkgeo.idToDet(*det)->toGlobal(locx);
250 GlobalPoint ypos = trkgeo.idToDet(*det)->toGlobal(locy);
256 double dzdr = posvect.
perp() > 0 ? (
dz.x() * posvect.
x() +
dz.y() * posvect.
y()) / posvect.
perp() : 0.;
257 double dxdr = posvect.
perp() > 0 ? (
dx.x() * posvect.
x() +
dx.y() * posvect.
y()) / posvect.
perp() : 0.;
258 double dydr = posvect.
perp() > 0 ? (
dy.x() * posvect.
x() +
dy.y() * posvect.
y()) / posvect.
perp() : 0.;
260 double dzdrphi = posvect.
perp() > 0 ? (
dz.y() * posvect.
x() -
dz.x() * posvect.
y()) / posvect.
perp() : 0.;
261 double dxdrphi = posvect.
perp() > 0 ? (
dx.y() * posvect.
x() -
dx.x() * posvect.
y()) / posvect.
perp() : 0.;
262 double dydrphi = posvect.
perp() > 0 ? (
dy.y() * posvect.
x() -
dy.x() * posvect.
y()) / posvect.
perp() : 0.;
264 for (std::map<unsigned int, DetIdSelector>::const_iterator
sel =
m_wantedsubdets.begin();
267 if (
sel->second.isSelected(*det)) {
279 (*m_zavedr)->Fill(
sel->first, dzdr);
281 (*m_zavedz)->Fill(
sel->first,
dz.z());
283 (*m_zavedrphi)->Fill(
sel->first, dzdrphi);
285 (*m_xavedr)->Fill(
sel->first, dxdr);
287 (*m_xavedz)->Fill(
sel->first,
dx.z());
289 (*m_xavedrphi)->Fill(
sel->first, dxdrphi);
291 (*m_yavedr)->Fill(
sel->first, dydr);
293 (*m_yavedz)->Fill(
sel->first,
dy.z());
295 (*m_yavedrphi)->Fill(
sel->first, dydrphi);
305 (*m_nchannels_ideal)->Reset();
307 (*m_nchannels_real)->Reset();
311 for (
const auto det : trkgeo.detUnits()) {
312 const StripGeomDetUnit* stripDet = dynamic_cast<const StripGeomDetUnit*>(det);
313 if (stripDet !=
nullptr) {
320 if (!stripQuality.IsStripBad(detid,
strip))
324 for (std::map<unsigned int, DetIdSelector>::const_iterator
sel =
m_wantedsubdets.begin();
327 if (
sel->second.isSelected(detid)) {
329 (*m_nchannels_ideal)->Fill(
sel->first, nchannideal);
331 (*m_nchannels_real)->Fill(
sel->first, nchannreal);
341 const std::vector<uint32_t>& pxldetids = pxlreader.
getAllDetIds();
343 for (std::vector<uint32_t>::const_iterator detid = pxldetids.begin(); detid != pxldetids.end(); ++detid) {
344 int nchannideal = pxlreader.getDetUnitDimensions(*detid).first * pxlreader.getDetUnitDimensions(*detid).second;
346 if (!pxlquality.IsModuleBad(*detid)) {
347 nchannreal = pxlreader.getDetUnitDimensions(*detid).first * pxlreader.getDetUnitDimensions(*detid).second;
356 for (std::map<unsigned int, DetIdSelector>::const_iterator
sel =
m_wantedsubdets.begin();
359 if (
sel->second.isSelected(*detid)) {
361 (*m_nchannels_ideal)->Fill(
sel->first, nchannideal);
363 (*m_nchannels_real)->Fill(
sel->first, nchannreal);
References PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, edm::FileInPath::fullPath(), GeomDet::geographicalId(), SiPixelDetInfoFileReader::getAllDetIds(), edm::EventSetup::getData(), LogDebug, LogTrace, m_averadius, m_avex, m_avey, m_avez, m_fp, m_nchannels_ideal, m_nchannels_real, m_pixelQualityToken, m_stripQualityToken, m_tkGeomToken, m_wantedsubdets, m_xavedr, m_xavedrphi, m_xavedz, m_yavedr, m_yavedrphi, m_yavedz, m_zavedr, m_zavedrphi, m_zavedz, StripTopology::nstrips(), GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, GeomDetEnumerators::P2OTB, GeomDetEnumerators::P2OTEC, GeomDetEnumerators::P2PXB, GeomDetEnumerators::P2PXEC, PV3DBase< T, PVType, FrameType >::perp(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, position, EgammaValidation_Wenu_cff::sel, StripGeomDetUnit::specificTopology(), digitizers_cfi::strip, GeomDetEnumerators::TEC, GeomDetEnumerators::TIB, GeomDetEnumerators::TID, GeomDetEnumerators::TOB, DetId::Tracker, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().