#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

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, selectors, 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 ( )

inlineoverride

Definition at line 66 of file ClusterSummaryProducer.h.

66 {};

Member Function Documentation

void ClusterSummaryProducer::beginStream ( edm::StreamID )

overrideprivate

Definition at line 127 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;
 }

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(), first, edm::Event::getByToken(), ClusterSummary::getClusChargeByIndex(), ClusterSummary::getClusSizeByIndex(), ClusterSummary::getNClusByIndex(), edmNew::DetSetVector< T >::id(), callgraph::module, moduleNames, eostools::move(), HLT_FULL_cff::pixelClusters, pixelClusters_, edm::Event::put(), ClusterSummary::reset(), mps_fire::result, selectors, HLT_FULL_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;
       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));
 }