#include <ClusterTPAssociationProducer.h>

Inheritance diagram for ClusterTPAssociationProducer:

Public Types
typedef std::vector < OmniClusterRef >	OmniClusterCollection

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	ClusterTPAssociationProducer (const edm::ParameterSet &)

	~ClusterTPAssociationProducer ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

template<typename T >
std::vector< std::pair < uint32_t, EncodedEventId > >	getSimTrackId (const edm::Handle< edm::DetSetVector< T > > &simLinks, const DetId &detId, uint32_t channel) const

virtual void	produce (edm::Event &, const edm::EventSetup &)

Additional Inherited Members
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

Detailed Description

Definition at line 24 of file ClusterTPAssociationProducer.h.

Member Typedef Documentation

typedef std::vector<OmniClusterRef> ClusterTPAssociationProducer::OmniClusterCollection

Definition at line 28 of file ClusterTPAssociationProducer.h.

Constructor & Destructor Documentation

ClusterTPAssociationProducer::ClusterTPAssociationProducer ( const edm::ParameterSet & cfg )

explicit

Definition at line 24 of file ClusterTPAssociationProducer.cc.

References edm::ParameterSet::getParameter(), pixelClustersToken_, sipixelSimLinksToken_, sistripSimLinksToken_, stripClustersToken_, and trackingParticleToken_.

   : _verbose(cfg.getParameter<bool>("verbose")),
     _pixelSimLinkSrc(cfg.getParameter<edm::InputTag>("pixelSimLinkSrc")),
     _stripSimLinkSrc(cfg.getParameter<edm::InputTag>("stripSimLinkSrc")),
     _pixelClusterSrc(cfg.getParameter<edm::InputTag>("pixelClusterSrc")),
     _stripClusterSrc(cfg.getParameter<edm::InputTag>("stripClusterSrc")),
     _trackingParticleSrc(cfg.getParameter<edm::InputTag>("trackingParticleSrc"))
 {
 
   sipixelSimLinksToken_ = consumes<edm::DetSetVector<PixelDigiSimLink> >(cfg.getParameter<edm::InputTag>("pixelSimLinkSrc"));
   sistripSimLinksToken_ = consumes<edm::DetSetVector<StripDigiSimLink> >(cfg.getParameter<edm::InputTag>("stripSimLinkSrc"));
   pixelClustersToken_ = consumes<edmNew::DetSetVector<SiPixelCluster> >(cfg.getParameter<edm::InputTag>("pixelClusterSrc"));
   stripClustersToken_ = consumes<edmNew::DetSetVector<SiStripCluster> >(cfg.getParameter<edm::InputTag>("stripClusterSrc"));
   trackingParticleToken_ = consumes<TrackingParticleCollection>(cfg.getParameter<edm::InputTag>("trackingParticleSrc"));
 
   produces<ClusterTPAssociationList>();
 }

ClusterTPAssociationProducer::~ClusterTPAssociationProducer ( )

Definition at line 42 of file ClusterTPAssociationProducer.cc.

42 {

43 }

Member Function Documentation

virtual void ClusterTPAssociationProducer::beginJob ( void )

inlineprivatevirtual

Definition at line 34 of file ClusterTPAssociationProducer.h.

34 {}

virtual void ClusterTPAssociationProducer::endJob ( void )

inlineprivatevirtual

Definition at line 36 of file ClusterTPAssociationProducer.h.

36 {}

template<typename T >

std::vector< std::pair< uint32_t, EncodedEventId > > ClusterTPAssociationProducer::getSimTrackId	(	const edm::Handle< edm::DetSetVector< T > > &	simLinks,
		const DetId &	detId,
		uint32_t	channel
	)		const

private

Definition at line 163 of file ClusterTPAssociationProducer.cc.

References edm::DetSet< T >::data, and edm::DetSet< T >::push_back().

 {
   //std::pair<uint32_t, EncodedEventId> simTrkId;
   std::vector<std::pair<uint32_t, EncodedEventId> > simTrkId;
   auto isearch = simLinks->find(detId);
   if (isearch != simLinks->end()) {
     // Loop over DigiSimLink in this det unit
     edm::DetSet<T> link_detset = (*isearch);
     for (typename edm::DetSet<T>::const_iterator it  = link_detset.data.begin(); 
                                                  it != link_detset.data.end(); ++it) {
       if (channel == it->channel()) {
         simTrkId.push_back(std::make_pair(it->SimTrackId(), it->eventId()));
       } 
     }
   }
   return simTrkId;
 }

void ClusterTPAssociationProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	es
	)

privatevirtual

Implements edm::stream::EDProducerBase.

Definition at line 45 of file ClusterTPAssociationProducer.cc.

References SiStripCluster::amplitudes(), edmNew::DetSet< T >::begin(), clusterTPAssociationListGreater(), cond::rpcobgas::detid, edmNew::DetSet< T >::end(), plotBeamSpotDB::first, SiStripCluster::firstStrip(), edm::Event::getByToken(), prof2calltree::last, edmNew::makeRefTo(), SiPixelCluster::maxPixelCol(), SiPixelCluster::maxPixelRow(), SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), HLT_25ns14e33_v1_cff::pixelClusters, pixelClustersToken_, PixelChannelIdentifier::pixelToChannel(), edm::Handle< T >::product(), edm::Event::put(), sipixelSimLinksToken_, sistripSimLinksToken_, python.multivaluedict::sort(), HLT_25ns14e33_v1_cff::stripClusters, stripClustersToken_, and trackingParticleToken_.

                                                                                     {
   auto clusterTPList = std::make_unique<ClusterTPAssociationList>();
  
   // Pixel DigiSimLink
   edm::Handle<edm::DetSetVector<PixelDigiSimLink> > sipixelSimLinks;
   //  iEvent.getByLabel(_pixelSimLinkSrc, sipixelSimLinks);
   iEvent.getByToken(sipixelSimLinksToken_,sipixelSimLinks);
 
   // SiStrip DigiSimLink
   edm::Handle<edm::DetSetVector<StripDigiSimLink> > sistripSimLinks;
   iEvent.getByToken(sistripSimLinksToken_,sistripSimLinks);
 
   // Pixel Cluster
   edm::Handle<edmNew::DetSetVector<SiPixelCluster> > pixelClusters;
   bool foundPixelClusters = iEvent.getByToken(pixelClustersToken_,pixelClusters);
 
   // Strip Cluster
   edm::Handle<edmNew::DetSetVector<SiStripCluster> > stripClusters;
   bool foundStripClusters = iEvent.getByToken(stripClustersToken_,stripClusters);
 
 
   // TrackingParticle
   edm::Handle<TrackingParticleCollection>  TPCollectionH;
   iEvent.getByToken(trackingParticleToken_,TPCollectionH);
 
   // prepare temporary map between SimTrackId and TrackingParticle index
   std::map<std::pair<size_t, EncodedEventId>, TrackingParticleRef> mapping;
   for (TrackingParticleCollection::size_type itp = 0;
                                              itp < TPCollectionH.product()->size(); ++itp) {
     TrackingParticleRef trackingParticle(TPCollectionH, itp);
 
     // SimTracks inside TrackingParticle
     EncodedEventId eid(trackingParticle->eventId());
     //size_t index = 0;
     for (std::vector<SimTrack>::const_iterator itrk  = trackingParticle->g4Track_begin(); 
                                                itrk != trackingParticle->g4Track_end(); ++itrk) {
       std::pair<uint32_t, EncodedEventId> trkid(itrk->trackId(), eid);
       //std::cout << "creating map for id: " << trkid.first << " with tp: " << trackingParticle.key() << std::endl;
       mapping.insert(std::make_pair(trkid, trackingParticle));
     }
   }
 
   if ( foundPixelClusters ) {
     // Pixel Clusters 
     for (edmNew::DetSetVector<SiPixelCluster>::const_iterator iter  = pixelClusters->begin(); 
                                                             iter != pixelClusters->end(); ++iter) {
       uint32_t detid = iter->id(); 
       DetId detId(detid);
       edmNew::DetSet<SiPixelCluster> link_pixel = (*iter);
       for (edmNew::DetSet<SiPixelCluster>::const_iterator di  = link_pixel.begin(); 
        di != link_pixel.end(); ++di) {
     const SiPixelCluster& cluster = (*di);
     edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> c_ref = 
       edmNew::makeRefTo(pixelClusters, di);
     
     std::set<std::pair<uint32_t, EncodedEventId> > simTkIds; 
     for (int irow = cluster.minPixelRow(); irow <= cluster.maxPixelRow(); ++irow) {
       for (int icol = cluster.minPixelCol(); icol <= cluster.maxPixelCol(); ++icol) {
         uint32_t channel = PixelChannelIdentifier::pixelToChannel(irow, icol);
         std::vector<std::pair<uint32_t, EncodedEventId> > trkid(getSimTrackId<PixelDigiSimLink>(sipixelSimLinks, detId, channel));
         if (trkid.size()==0) continue; 
         simTkIds.insert(trkid.begin(),trkid.end());
       }
     }
     for (std::set<std::pair<uint32_t, EncodedEventId> >::const_iterator iset  = simTkIds.begin(); 
          iset != simTkIds.end(); iset++) {
       auto ipos = mapping.find(*iset);
       if (ipos != mapping.end()) {
         //std::cout << "cluster in detid: " << detid << " from tp: " << ipos->second.key() << " " << iset->first << std::endl;
         clusterTPList->push_back(std::make_pair(OmniClusterRef(c_ref), ipos->second));
       }
     }
       }
     }
   }
 
   if ( foundStripClusters ) {
     // Strip Clusters
     for (edmNew::DetSetVector<SiStripCluster>::const_iterator iter  = stripClusters->begin(false), eter = stripClusters->end(false); 
      iter != eter; ++iter) {
       if (!(*iter).isValid()) continue;
       uint32_t detid = iter->id();  
       DetId detId(detid);
       edmNew::DetSet<SiStripCluster> link_strip = (*iter);
       for (edmNew::DetSet<SiStripCluster>::const_iterator di  = link_strip.begin(); 
        di != link_strip.end(); di++) {
     const SiStripCluster& cluster = (*di);
     edm::Ref<edmNew::DetSetVector<SiStripCluster>, SiStripCluster> c_ref = 
       edmNew::makeRefTo(stripClusters, di);
     
     std::set<std::pair<uint32_t, EncodedEventId> > simTkIds; 
     int first  = cluster.firstStrip();     
     int last   = first + cluster.amplitudes().size();
     
     for (int istr = first; istr < last; ++istr) {
       std::vector<std::pair<uint32_t, EncodedEventId> > trkid(getSimTrackId<StripDigiSimLink>(sistripSimLinks, detId, istr));
       if (trkid.size()==0) continue; 
       simTkIds.insert(trkid.begin(),trkid.end());
     }
     for (std::set<std::pair<uint32_t, EncodedEventId> >::const_iterator iset  = simTkIds.begin(); 
          iset != simTkIds.end(); iset++) {
       auto ipos = mapping.find(*iset);
       if (ipos != mapping.end()) {
         //std::cout << "cluster in detid: " << detid << " from tp: " << ipos->second.key() << " " << iset->first << std::endl;
         clusterTPList->push_back(std::make_pair(OmniClusterRef(c_ref), ipos->second));
       } 
     }
       } 
     }
   }
 
   std::sort(clusterTPList->begin(), clusterTPList->end(), clusterTPAssociationListGreater);
   iEvent.put(std::move(clusterTPList));
 }

Member Data Documentation

edm::InputTag ClusterTPAssociationProducer::_pixelClusterSrc

private

Definition at line 45 of file ClusterTPAssociationProducer.h.

edm::InputTag ClusterTPAssociationProducer::_pixelSimLinkSrc

private

Definition at line 43 of file ClusterTPAssociationProducer.h.

edm::InputTag ClusterTPAssociationProducer::_stripClusterSrc

private

Definition at line 46 of file ClusterTPAssociationProducer.h.

edm::InputTag ClusterTPAssociationProducer::_stripSimLinkSrc

private

Definition at line 44 of file ClusterTPAssociationProducer.h.

edm::InputTag ClusterTPAssociationProducer::_trackingParticleSrc

private

Definition at line 47 of file ClusterTPAssociationProducer.h.

bool ClusterTPAssociationProducer::_verbose

private

Definition at line 42 of file ClusterTPAssociationProducer.h.

Referenced by cmsPerfSuite.PerfSuite::mkCandleDir(), cmsPerfSuite.PerfSuite::optionParse(), cmsPerfSuite.PerfSuite::printFlush(), cmsPerfSuite.PerfSuite::runcmd(), and cmsPerfSuite.PerfSuite::runCmdSet().

edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > ClusterTPAssociationProducer::pixelClustersToken_

private

Definition at line 51 of file ClusterTPAssociationProducer.h.

Referenced by ClusterTPAssociationProducer(), and produce().

edm::EDGetTokenT<edm::DetSetVector<PixelDigiSimLink> > ClusterTPAssociationProducer::sipixelSimLinksToken_

private

Definition at line 49 of file ClusterTPAssociationProducer.h.

Referenced by ClusterTPAssociationProducer(), and produce().

edm::EDGetTokenT<edm::DetSetVector<StripDigiSimLink> > ClusterTPAssociationProducer::sistripSimLinksToken_

private

Definition at line 50 of file ClusterTPAssociationProducer.h.

Referenced by ClusterTPAssociationProducer(), and produce().

edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > ClusterTPAssociationProducer::stripClustersToken_

private

Definition at line 52 of file ClusterTPAssociationProducer.h.

Referenced by ClusterTPAssociationProducer(), and produce().

edm::EDGetTokenT<TrackingParticleCollection> ClusterTPAssociationProducer::trackingParticleToken_

private

Definition at line 53 of file ClusterTPAssociationProducer.h.

Referenced by ClusterTPAssociationProducer(), and produce().

Private Attributes
edm::InputTag	_pixelClusterSrc

edm::InputTag	_pixelSimLinkSrc

edm::InputTag	_stripClusterSrc

edm::InputTag	_stripSimLinkSrc

edm::InputTag	_trackingParticleSrc

bool	_verbose

edm::EDGetTokenT < edmNew::DetSetVector < SiPixelCluster > >	pixelClustersToken_

edm::EDGetTokenT < edm::DetSetVector < PixelDigiSimLink > >	sipixelSimLinksToken_

edm::EDGetTokenT < edm::DetSetVector < StripDigiSimLink > >	sistripSimLinksToken_

edm::EDGetTokenT < edmNew::DetSetVector < SiStripCluster > >	stripClustersToken_

edm::EDGetTokenT < TrackingParticleCollection >	trackingParticleToken_

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation