OccupancyPlots Class Reference

#include <myTKAnalyses/DigiInvestigator/src/OccupancyPlots.cc>

Inheritance diagram for OccupancyPlots:

Public Member Functions
	OccupancyPlots (const edm::ParameterSet &)
	~OccupancyPlots () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	beginRun (const edm::Run &, const edm::EventSetup &) override
void	endJob () override
void	endRun (const edm::Run &, const edm::EventSetup &) override

Private Attributes
TProfile **	m_avemultiplicity
TProfile **	m_aveoccupancy
TProfile **	m_averadius
TProfile **	m_avex
TProfile **	m_avey
TProfile **	m_avez
edm::FileInPath	m_fp
std::vector< edm::EDGetTokenT < std::map< unsigned int, int > > >	m_multiplicityMapTokens
TH1F **	m_nchannels_ideal
TH1F **	m_nchannels_real
std::vector< edm::EDGetTokenT < std::map< unsigned int, int > > >	m_occupancyMapTokens
edm::ESGetToken < SiPixelQuality, SiPixelQualityRcd >	m_pixelQualityToken
RunHistogramManager	m_rhm
edm::ESGetToken < SiStripQuality, SiStripQualityRcd >	m_stripQualityToken
edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	m_tkGeomToken
std::map< unsigned int, DetIdSelector >	m_wantedsubdets
TProfile **	m_xavedr
TProfile **	m_xavedrphi
TProfile **	m_xavedz
TProfile **	m_yavedr
TProfile **	m_yavedrphi
TProfile **	m_yavedz
TProfile **	m_zavedr
TProfile **	m_zavedrphi
TProfile **	m_zavedz

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 66 of file OccupancyPlots.cc.

Constructor & Destructor Documentation

OccupancyPlots::OccupancyPlots ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 122 of file OccupancyPlots.cc.

References GlobalPosition_Frontier_DevDB_cff::tag.

    : m_multiplicityMapTokens(edm::vector_transform(
          iConfig.getParameter<std::vector<edm::InputTag> >("multiplicityMaps"),
          [this](edm::InputTag const& tag) { return consumes<std::map<unsigned int, int> >(tag); })),

OccupancyPlots::~OccupancyPlots ( )

override

Definition at line 170 of file OccupancyPlots.cc.

                                {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

void OccupancyPlots::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 180 of file OccupancyPlots.cc.

References edm::Event::getByToken(), m_avemultiplicity, m_aveoccupancy, m_multiplicityMapTokens, m_occupancyMapTokens, VarParsing::mult, and trackerHitRTTI::vector.

                                                                                {
  using namespace edm;

  for (std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > >::const_iterator mapToken =
           m_multiplicityMapTokens.begin();
       mapToken != m_multiplicityMapTokens.end();
       ++mapToken) {
    Handle<std::map<unsigned int, int> > mults;
    iEvent.getByToken(*mapToken, mults);

    for (std::map<unsigned int, int>::const_iterator mult = mults->begin(); mult != mults->end(); mult++) {
      if (m_avemultiplicity && *m_avemultiplicity)
        (*m_avemultiplicity)->Fill(mult->first, mult->second);
    }
  }

  for (std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > >::const_iterator mapToken =
           m_occupancyMapTokens.begin();
       mapToken != m_occupancyMapTokens.end();
       ++mapToken) {
    Handle<std::map<unsigned int, int> > occus;
    iEvent.getByToken(*mapToken, occus);

    for (std::map<unsigned int, int>::const_iterator occu = occus->begin(); occu != occus->end(); occu++) {
      if (m_aveoccupancy && *m_aveoccupancy)
        (*m_aveoccupancy)->Fill(occu->first, occu->second);
    }
  }
}

void OccupancyPlots::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 211 of file OccupancyPlots.cc.

{}

void OccupancyPlots::beginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 213 of file OccupancyPlots.cc.

References RunHistogramManager::beginRun(), and m_rhm.

{ m_rhm.beginRun(iRun); }

void OccupancyPlots::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 369 of file OccupancyPlots.cc.

{}

void OccupancyPlots::endRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 215 of file OccupancyPlots.cc.

References PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, edm::FileInPath::fullPath(), GeomDet::geographicalId(), SiPixelDetInfoFileReader::getAllDetIds(), edm::EventSetup::getData(), if(), 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(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                           {
  const auto& trkgeo = iSetup.getData(m_tkGeomToken);

  // Test new TrackerGeometry features
  LogDebug("IsThereTest") << "Test of TrackerGeometry::isThere";
  LogTrace("IsThereTest") << " is there PixelBarrel: " << trkgeo.isThere(GeomDetEnumerators::PixelBarrel);
  LogTrace("IsThereTest") << " is there PixelEndcap: " << trkgeo.isThere(GeomDetEnumerators::PixelEndcap);
  LogTrace("IsThereTest") << " is there P1PXB: " << trkgeo.isThere(GeomDetEnumerators::P1PXB);
  LogTrace("IsThereTest") << " is there P1PXEC: " << trkgeo.isThere(GeomDetEnumerators::P1PXEC);
  LogTrace("IsThereTest") << " is there P2PXB: " << trkgeo.isThere(GeomDetEnumerators::P2PXB);
  LogTrace("IsThereTest") << " is there P2PXEC: " << trkgeo.isThere(GeomDetEnumerators::P2PXEC);
  LogTrace("IsThereTest") << " is there TIB: " << trkgeo.isThere(GeomDetEnumerators::TIB);
  LogTrace("IsThereTest") << " is there TID: " << trkgeo.isThere(GeomDetEnumerators::TID);
  LogTrace("IsThereTest") << " is there TOB: " << trkgeo.isThere(GeomDetEnumerators::TOB);
  LogTrace("IsThereTest") << " is there TEC: " << trkgeo.isThere(GeomDetEnumerators::TEC);
  LogTrace("IsThereTest") << " is there P2OTB: " << trkgeo.isThere(GeomDetEnumerators::P2OTB);
  LogTrace("IsThereTest") << " is there P2OTEC: " << trkgeo.isThere(GeomDetEnumerators::P2OTEC);

  const Local2DPoint center(0., 0.);
  const Local3DPoint locz(0., 0., 1.);
  const Local3DPoint locx(1., 0., 0.);
  const Local3DPoint locy(0., 1., 0.);
  const GlobalPoint origin(0., 0., 0.);

  TrackingGeometry::DetIdContainer detunits = trkgeo.detUnitIds();

  for (TrackingGeometry::DetIdContainer::const_iterator det = detunits.begin(); det != detunits.end(); ++det) {
    if (det->det() != DetId::Tracker)
      continue;

    edm::LogInfo("DetIdFromGeometry") << det->rawId();

    GlobalPoint position = trkgeo.idToDet(*det)->toGlobal(center);
    GlobalPoint zpos = trkgeo.idToDet(*det)->toGlobal(locz);
    GlobalPoint xpos = trkgeo.idToDet(*det)->toGlobal(locx);
    GlobalPoint ypos = trkgeo.idToDet(*det)->toGlobal(locy);
    GlobalVector posvect = position - origin;
    GlobalVector dz = zpos - position;
    GlobalVector dx = xpos - position;
    GlobalVector dy = ypos - position;

    double dzdr = posvect.perp() > 0 ? (dz.x() * posvect.x() + dz.y() * posvect.y()) / posvect.perp() : 0.;
    double dxdr = posvect.perp() > 0 ? (dx.x() * posvect.x() + dx.y() * posvect.y()) / posvect.perp() : 0.;
    double dydr = posvect.perp() > 0 ? (dy.x() * posvect.x() + dy.y() * posvect.y()) / posvect.perp() : 0.;

    double dzdrphi = posvect.perp() > 0 ? (dz.y() * posvect.x() - dz.x() * posvect.y()) / posvect.perp() : 0.;
    double dxdrphi = posvect.perp() > 0 ? (dx.y() * posvect.x() - dx.x() * posvect.y()) / posvect.perp() : 0.;
    double dydrphi = posvect.perp() > 0 ? (dy.y() * posvect.x() - dy.x() * posvect.y()) / posvect.perp() : 0.;

    for (std::map<unsigned int, DetIdSelector>::const_iterator sel = m_wantedsubdets.begin();
         sel != m_wantedsubdets.end();
         ++sel) {
      if (sel->second.isSelected(*det)) {
        edm::LogInfo("SelectedDetId") << sel->first;
        // average positions
        if (m_averadius && *m_averadius)
          (*m_averadius)->Fill(sel->first, position.perp());
        if (m_avez && *m_avez)
          (*m_avez)->Fill(sel->first, position.z());
        if (m_avex && *m_avex)
          (*m_avex)->Fill(sel->first, position.x());
        if (m_avey && *m_avey)
          (*m_avey)->Fill(sel->first, position.y());
        if (m_zavedr && *m_zavedr)
          (*m_zavedr)->Fill(sel->first, dzdr);
        if (m_zavedz && *m_zavedz)
          (*m_zavedz)->Fill(sel->first, dz.z());
        if (m_zavedrphi && *m_zavedrphi)
          (*m_zavedrphi)->Fill(sel->first, dzdrphi);
        if (m_xavedr && *m_xavedr)
          (*m_xavedr)->Fill(sel->first, dxdr);
        if (m_xavedz && *m_xavedz)
          (*m_xavedz)->Fill(sel->first, dx.z());
        if (m_xavedrphi && *m_xavedrphi)
          (*m_xavedrphi)->Fill(sel->first, dxdrphi);
        if (m_yavedr && *m_yavedr)
          (*m_yavedr)->Fill(sel->first, dydr);
        if (m_yavedz && *m_yavedz)
          (*m_yavedz)->Fill(sel->first, dy.z());
        if (m_yavedrphi && *m_yavedrphi)
          (*m_yavedrphi)->Fill(sel->first, dydrphi);
      }
    }
  }

  // counting the number of channels per module subset

  // the histograms have to be reset to avoid double counting if endRun is called more than once

  if (m_nchannels_ideal && *m_nchannels_ideal)
    (*m_nchannels_ideal)->Reset();
  if (m_nchannels_real && *m_nchannels_real)
    (*m_nchannels_real)->Reset();

  const auto& stripQuality = iSetup.getData(m_stripQualityToken);

  for (const auto det : trkgeo.detUnits()) {
    const StripGeomDetUnit* stripDet = dynamic_cast<const StripGeomDetUnit*>(det);
    if (stripDet != nullptr) {
      const DetId detid = stripDet->geographicalId();

      int nchannideal = stripDet->specificTopology().nstrips();
      //     int nchannreal = stripDet->specificTopology().nstrips();
      int nchannreal = 0;
      for (int strip = 0; strip < nchannideal; ++strip) {
        if (!stripQuality.IsStripBad(detid, strip))
          ++nchannreal;
      }

      for (std::map<unsigned int, DetIdSelector>::const_iterator sel = m_wantedsubdets.begin();
           sel != m_wantedsubdets.end();
           ++sel) {
        if (sel->second.isSelected(detid)) {
          if (m_nchannels_ideal && *m_nchannels_ideal)
            (*m_nchannels_ideal)->Fill(sel->first, nchannideal);
          if (m_nchannels_real && *m_nchannels_real)
            (*m_nchannels_real)->Fill(sel->first, nchannreal);
        }
      }
    }
  }

  const auto& pxlquality = iSetup.getData(m_pixelQualityToken);

  SiPixelDetInfoFileReader pxlreader(m_fp.fullPath());

  const std::vector<uint32_t>& pxldetids = pxlreader.getAllDetIds();

  for (std::vector<uint32_t>::const_iterator detid = pxldetids.begin(); detid != pxldetids.end(); ++detid) {
    int nchannideal = pxlreader.getDetUnitDimensions(*detid).first * pxlreader.getDetUnitDimensions(*detid).second;
    int nchannreal = 0;
    if (!pxlquality.IsModuleBad(*detid)) {
      nchannreal = pxlreader.getDetUnitDimensions(*detid).first * pxlreader.getDetUnitDimensions(*detid).second;
    }
    /*
     int nchannreal = 0;
     for(int strip = 0; strip < nchannideal; ++strip) {
       if(!stripquality.IsStripBad(*detid,strip)) ++nchannreal;
     }
     */

    for (std::map<unsigned int, DetIdSelector>::const_iterator sel = m_wantedsubdets.begin();
         sel != m_wantedsubdets.end();
         ++sel) {
      if (sel->second.isSelected(*detid)) {
        if (m_nchannels_ideal && *m_nchannels_ideal)
          (*m_nchannels_ideal)->Fill(sel->first, nchannideal);
        if (m_nchannels_real && *m_nchannels_real)
          (*m_nchannels_real)->Fill(sel->first, nchannreal);
      }
    }
  }
}

Member Data Documentation

TProfile** OccupancyPlots::m_avemultiplicity

private

Definition at line 90 of file OccupancyPlots.cc.

Referenced by analyze().

TProfile** OccupancyPlots::m_aveoccupancy

private

Definition at line 91 of file OccupancyPlots.cc.

Referenced by analyze().

TProfile** OccupancyPlots::m_averadius

private

Definition at line 96 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_avex

private

Definition at line 98 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_avey

private

Definition at line 99 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_avez

private

Definition at line 97 of file OccupancyPlots.cc.

Referenced by endRun().

edm::FileInPath OccupancyPlots::m_fp

private

Definition at line 82 of file OccupancyPlots.cc.

Referenced by endRun().

std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > > OccupancyPlots::m_multiplicityMapTokens

private

Definition at line 80 of file OccupancyPlots.cc.

Referenced by analyze().

TH1F** OccupancyPlots::m_nchannels_ideal

private

Definition at line 93 of file OccupancyPlots.cc.

Referenced by endRun().

TH1F** OccupancyPlots::m_nchannels_real

private

Definition at line 94 of file OccupancyPlots.cc.

Referenced by endRun().

std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > > OccupancyPlots::m_occupancyMapTokens

private

Definition at line 81 of file OccupancyPlots.cc.

Referenced by analyze().

edm::ESGetToken<SiPixelQuality, SiPixelQualityRcd> OccupancyPlots::m_pixelQualityToken

private

Definition at line 85 of file OccupancyPlots.cc.

Referenced by endRun().

RunHistogramManager OccupancyPlots::m_rhm

private

Definition at line 87 of file OccupancyPlots.cc.

Referenced by beginRun().

edm::ESGetToken<SiStripQuality, SiStripQualityRcd> OccupancyPlots::m_stripQualityToken

private

Definition at line 84 of file OccupancyPlots.cc.

Referenced by endRun().

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> OccupancyPlots::m_tkGeomToken

private

Definition at line 83 of file OccupancyPlots.cc.

Referenced by endRun().

std::map<unsigned int, DetIdSelector> OccupancyPlots::m_wantedsubdets

private

Definition at line 88 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_xavedr

private

Definition at line 106 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_xavedrphi

private

Definition at line 108 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_xavedz

private

Definition at line 107 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_yavedr

private

Definition at line 103 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_yavedrphi

private

Definition at line 105 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_yavedz

private

Definition at line 104 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_zavedr

private

Definition at line 100 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_zavedrphi

private

Definition at line 102 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_zavedz

private

Definition at line 101 of file OccupancyPlots.cc.

Referenced by endRun().