#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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 68 of file ClusterSummaryProducer.h.

Member Typedef Documentation

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

private

Definition at line 78 of file ClusterSummaryProducer.h.

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

private

Definition at line 79 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, 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 71 of file ClusterSummaryProducer.h.

References beginStream(), and produce().

71 {};

Member Function Documentation

void ClusterSummaryProducer::beginStream ( edm::StreamID )

overrideprivate

Definition at line 123 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 48 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(), plotBeamSpotDB::first, objects.autophobj::float, edm::Event::getByToken(), ClusterSummary::getClusChargeByIndex(), ClusterSummary::getClusSizeByIndex(), ClusterSummary::getNClusByIndex(), edmNew::DetSetVector< T >::id(), python.rootplot.argparse::module, moduleNames, eostools::move(), HiRegitMuonDetachedTripletStep_cff::pixelClusters, pixelClusters_, edm::Event::put(), ClusterSummary::reset(), mps_fire::result, selectors, HiRegitMuonDetachedTripletStep_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));
 }