ClusterSummaryProducer Class Reference

#include <msegala/ClusterSummaryProducer/src/ClusterSummaryProducer.cc>

Inheritance diagram for ClusterSummaryProducer:

Public Member Functions
	ClusterSummaryProducer (const edm::ParameterSet &)
	~ClusterSummaryProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Types
typedef std::pair< DetIdSelector, ClusterSummary::CMSTracker >	ModuleSelection
typedef std::vector< ModuleSelection >	ModuleSelections

Private Member Functions
void	beginStream (edm::StreamID) override
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
ClusterSummary	cCluster
bool	doPixels
bool	doStrips
std::map< std::string, int >	EnumMap
std::vector< std::string >	moduleNames
edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > >	pixelClusters_
ModuleSelections	selectors
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	stripClusters_
bool	verbose

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 63 of file ClusterSummaryProducer.h.

Member Typedef Documentation

ModuleSelection

typedef std::pair<DetIdSelector, ClusterSummary::CMSTracker> ClusterSummaryProducer::ModuleSelection

private

Definition at line 72 of file ClusterSummaryProducer.h.

ModuleSelections

typedef std::vector<ModuleSelection> ClusterSummaryProducer::ModuleSelections

private

Definition at line 73 of file ClusterSummaryProducer.h.

Constructor & Destructor Documentation

ClusterSummaryProducer()

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

explicit

Definition at line 3 of file ClusterSummaryProducer.cc.

References cCluster, Exception, edm::ParameterSet::getParameter(), moduleNames, ClusterSummary::NTRACKERENUMS, ClusterSummary::NVALIDENUMS, pixelClusters_, corrVsCorr::selection, AlCaHLTBitMon_QueryRunRegistry::string, ClusterSummary::STRIP, stripClusters_, ClusterSummary::subDetNames, ClusterSummary::subDetSelections, and verbose.

    : doStrips(iConfig.getParameter<bool>("doStrips")),
      doPixels(iConfig.getParameter<bool>("doPixels")),
      verbose(iConfig.getParameter<bool>("verbose")) {
  pixelClusters_ =
      consumes<edmNew::DetSetVector<SiPixelCluster> >(iConfig.getParameter<edm::InputTag>("pixelClusters"));
  stripClusters_ =
      consumes<edmNew::DetSetVector<SiStripCluster> >(iConfig.getParameter<edm::InputTag>("stripClusters"));

  ClusterSummary::CMSTracker maxEnum = ClusterSummary::STRIP;

  std::vector<std::string> wantedsubdets = iConfig.getParameter<std::vector<std::string> >("wantedSubDets");
  for (const auto& iS : wantedsubdets) {
    ClusterSummary::CMSTracker subdet = ClusterSummary::NVALIDENUMS;
    for (int iN = 0; iN < ClusterSummary::NVALIDENUMS; ++iN)
      if (ClusterSummary::subDetNames[iN] == iS)
        subdet = ClusterSummary::CMSTracker(iN);
    if (subdet == ClusterSummary::NVALIDENUMS)
      throw cms::Exception("No standard selection: ") << iS;

    selectors.push_back(ModuleSelection(DetIdSelector(ClusterSummary::subDetSelections[subdet]), subdet));
    if (subdet > maxEnum)
      maxEnum = subdet;
    if (verbose)
      moduleNames.push_back(ClusterSummary::subDetNames[subdet]);
  }

  std::vector<edm::ParameterSet> wantedusersubdets_ps =
      iConfig.getParameter<std::vector<edm::ParameterSet> >("wantedUserSubDets");
  for (const auto& iS : wantedusersubdets_ps) {
    ClusterSummary::CMSTracker subdet = (ClusterSummary::CMSTracker)iS.getParameter<unsigned int>("detSelection");
    std::string detname = iS.getParameter<std::string>("detLabel");
    std::vector<std::string> selection = iS.getParameter<std::vector<std::string> >("selection");

    if (subdet <= ClusterSummary::NVALIDENUMS)
      throw cms::Exception("Already predefined selection: ") << subdet;
    if (subdet >= ClusterSummary::NTRACKERENUMS)
      throw cms::Exception("Selection is out of range: ") << subdet;

    selectors.push_back(ModuleSelection(DetIdSelector(selection), subdet));
    if (subdet > maxEnum)
      maxEnum = subdet;
    if (verbose)
      moduleNames.push_back(detname);
  }

  cCluster = ClusterSummary(maxEnum + 1);
  produces<ClusterSummary>().setBranchAlias("trackerClusterSummary");
}

~ClusterSummaryProducer()

ClusterSummaryProducer::~ClusterSummaryProducer ( )

inlineoverride

Definition at line 66 of file ClusterSummaryProducer.h.

{}

Member Function Documentation

beginStream()

void ClusterSummaryProducer::beginStream ( edm::StreamID  )

overrideprivate

Definition at line 135 of file ClusterSummaryProducer.cc.

References moduleNames, and verbose.

                                                    {
  if (!verbose)
    return;
  edm::LogInfo("ClusterSummaryProducer") << "+++++++++++++++++++++++++++++++ " << std::endl << "Getting info on ";
  for (const auto& iS : moduleNames) {
    edm::LogInfo("ClusterSummaryProducer") << iS << " ";
  }
  edm::LogInfo("ClusterSummaryProducer") << std::endl;
}

produce()

void ClusterSummaryProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 53 of file ClusterSummaryProducer.cc.

References ClusterSummary::addClusChargeByIndex(), ClusterSummary::addClusSizeByIndex(), ClusterSummary::addNClusByIndex(), edmNew::DetSetVector< T >::begin(), cCluster, ClusterVariables::charge(), ClusterVariables::clusterSize(), nano_mu_local_reco_cff::clusterSize, ALCARECOPPSCalTrackBasedSel_cff::detid, doPixels, doStrips, edmNew::DetSetVector< T >::end(), dqmdumpme::first, ALCARECOEcalPhiSym_cff::float, ClusterSummary::getClusChargeByIndex(), ClusterSummary::getClusSizeByIndex(), ClusterSummary::getNClusByIndex(), edmNew::DetSetVector< T >::id(), iEvent, callgraph::module, moduleNames, eostools::move(), DQMOffline_LumiMonitoring_cff::pixelClusters, pixelClusters_, ClusterSummary::reset(), mps_fire::result, HLT_2024v14_cff::stripClusters, stripClusters_, and verbose.

                                                                                  {
  using namespace edm;
  cCluster.reset();
  std::vector<bool> selectedVector(selectors.size(), false);

  auto getSelections = [&](const uint32_t detid) {
    for (unsigned int iS = 0; iS < selectors.size(); ++iS)
      selectedVector[iS] = selectors[iS].first.isSelected(detid);
  };
  auto fillSelections = [&](const int clusterSize, const float clusterCharge) {
    for (unsigned int iS = 0; iS < selectors.size(); ++iS) {
      if (!selectedVector[iS])
        continue;
      const ClusterSummary::CMSTracker module = selectors[iS].second;
      cCluster.addNClusByIndex(module, 1);
      cCluster.addClusSizeByIndex(module, clusterSize);
      cCluster.addClusChargeByIndex(module, clusterCharge);
    }
  };

  //===================++++++++++++========================
  //                   For SiStrips
  //===================++++++++++++========================
  if (doStrips) {
    edm::Handle<edmNew::DetSetVector<SiStripCluster> > stripClusters;
    iEvent.getByToken(stripClusters_, stripClusters);
    edmNew::DetSetVector<SiStripCluster>::const_iterator itClusters = stripClusters->begin();
    for (; itClusters != stripClusters->end(); ++itClusters) {
      getSelections(itClusters->id());
      for (edmNew::DetSet<SiStripCluster>::const_iterator cluster = itClusters->begin(); cluster != itClusters->end();
           ++cluster) {
        const ClusterVariables Summaryinfo(*cluster);
        fillSelections(Summaryinfo.clusterSize(), Summaryinfo.charge());
      }
    }
  }

  //===================++++++++++++========================
  //                   For SiPixels
  //===================++++++++++++========================
  if (doPixels) {
    edm::Handle<edmNew::DetSetVector<SiPixelCluster> > pixelClusters;
    iEvent.getByToken(pixelClusters_, pixelClusters);
    edmNew::DetSetVector<SiPixelCluster>::const_iterator itClusters = pixelClusters->begin();
    for (; itClusters != pixelClusters->end(); ++itClusters) {
      getSelections(itClusters->id());
      for (edmNew::DetSet<SiPixelCluster>::const_iterator cluster = itClusters->begin(); cluster != itClusters->end();
           ++cluster) {
        fillSelections(cluster->size(), float(cluster->charge()) / 1000.);
      }
    }
  }

  //===================+++++++++++++========================
  //                   Fill Producer
  //===================+++++++++++++========================
  if (verbose) {
    for (const auto& iS : selectors) {
      const ClusterSummary::CMSTracker module = iS.second;
      if (cCluster.getNClusByIndex(module) != 0) {
        edm::LogInfo("ClusterSummaryProducer")
            << "n" << moduleNames[module] << ", avg size, avg charge = " << cCluster.getNClusByIndex(module) << ", "
            << cCluster.getClusSizeByIndex(module) / cCluster.getNClusByIndex(module) << ", "
            << cCluster.getClusChargeByIndex(module) / cCluster.getNClusByIndex(module) << std::endl;
      } else {
        edm::LogInfo("ClusterSummaryProducer")
            << "n" << moduleNames[module] << ", avg size, avg charge = " << cCluster.getNClusByIndex(module) << ", "
            << "0"
            << ", "
            << "0" << std::endl;
      }
      edm::LogInfo("ClusterSummaryProducer") << "-------------------------------------------------------" << std::endl;
    }
  }

  //Put the filled class into the producer
  auto result = std::make_unique<ClusterSummary>();
  //Cleanup empty selections
  result->copyNonEmpty(cCluster);
  iEvent.put(std::move(result));
}

Member Data Documentation

cCluster

ClusterSummary ClusterSummaryProducer::cCluster

private

Definition at line 83 of file ClusterSummaryProducer.h.

Referenced by ClusterSummaryProducer(), and produce().

doPixels

bool ClusterSummaryProducer::doPixels

private

Definition at line 87 of file ClusterSummaryProducer.h.

Referenced by produce().

doStrips

bool ClusterSummaryProducer::doStrips

private

Definition at line 86 of file ClusterSummaryProducer.h.

Referenced by produce().

EnumMap

std::map<std::string, int> ClusterSummaryProducer::EnumMap

private

Definition at line 84 of file ClusterSummaryProducer.h.

moduleNames

std::vector<std::string> ClusterSummaryProducer::moduleNames

private

Definition at line 81 of file ClusterSummaryProducer.h.

Referenced by beginStream(), ClusterSummaryProducer(), and produce().

pixelClusters_

edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > ClusterSummaryProducer::pixelClusters_

private

Definition at line 77 of file ClusterSummaryProducer.h.

Referenced by ClusterSummaryProducer(), and produce().

selectors

ModuleSelections ClusterSummaryProducer::selectors

private

Definition at line 80 of file ClusterSummaryProducer.h.

stripClusters_

edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > ClusterSummaryProducer::stripClusters_

private

Definition at line 78 of file ClusterSummaryProducer.h.

Referenced by ClusterSummaryProducer(), and produce().

verbose

bool ClusterSummaryProducer::verbose

private

Definition at line 88 of file ClusterSummaryProducer.h.

Referenced by beginStream(), ClusterSummaryProducer(), and produce().