#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 63 of file ClusterSummaryProducer.h.

Member Typedef Documentation

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

private

Definition at line 72 of file ClusterSummaryProducer.h.

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

private

Definition at line 73 of file ClusterSummaryProducer.h.

Constructor & Destructor Documentation

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, and ClusterSummary::subDetSelections.

     : 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 ( )

inlineoverride

Definition at line 66 of file ClusterSummaryProducer.h.

References beginStream(), and produce().

66 {};

Member Function Documentation

void ClusterSummaryProducer::beginStream ( edm::StreamID )

overrideprivate

Definition at line 127 of file ClusterSummaryProducer.cc.

References DEFINE_FWK_MODULE, and moduleNames.

Referenced by ~ClusterSummaryProducer().

                                                     {
   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;
 }

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(), gather_cfg::cout, doPixels, doStrips, edmNew::DetSetVector< T >::end(), dqmdumpme::first, dqmMemoryStats::float, edm::Event::getByToken(), ClusterSummary::getClusChargeByIndex(), ClusterSummary::getClusSizeByIndex(), ClusterSummary::getNClusByIndex(), edmNew::DetSetVector< T >::id(), moduleNames, eostools::move(), LumiMonitor_cff::pixelClusters, pixelClusters_, edm::Event::put(), ClusterSummary::reset(), mps_fire::result, HLT_2018_cff::stripClusters, and stripClusters_.

Referenced by ~ClusterSummaryProducer().

                                                                                   {
   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;
       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;
     }
     std::cout << "-------------------------------------------------------" << 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));
 }