Inheritance diagram for ClusterTPAssociationProducer:

Public Types
typedef std::vector< OmniClusterRef >	OmniClusterCollection

Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

	~ClusterTPAssociationProducer () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
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

void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override

Additional Inherited Members
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 32 of file ClusterTPAssociationProducer.cc.

Member Typedef Documentation

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

Definition at line 34 of file ClusterTPAssociationProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 59 of file ClusterTPAssociationProducer.cc.

     : sipixelSimLinksToken_(
           consumes<edm::DetSetVector<PixelDigiSimLink> >(cfg.getParameter<edm::InputTag>("pixelSimLinkSrc"))),
       sistripSimLinksToken_(
           consumes<edm::DetSetVector<StripDigiSimLink> >(cfg.getParameter<edm::InputTag>("stripSimLinkSrc"))),
       siphase2OTSimLinksToken_(
           consumes<edm::DetSetVector<PixelDigiSimLink> >(cfg.getParameter<edm::InputTag>("phase2OTSimLinkSrc"))),
       pixelClustersToken_(
           consumes<edmNew::DetSetVector<SiPixelCluster> >(cfg.getParameter<edm::InputTag>("pixelClusterSrc"))),
       stripClustersToken_(
           consumes<edmNew::DetSetVector<SiStripCluster> >(cfg.getParameter<edm::InputTag>("stripClusterSrc"))),
       phase2OTClustersToken_(consumes<edmNew::DetSetVector<Phase2TrackerCluster1D> >(
           cfg.getParameter<edm::InputTag>("phase2OTClusterSrc"))),
       trackingParticleToken_(
           consumes<TrackingParticleCollection>(cfg.getParameter<edm::InputTag>("trackingParticleSrc"))),
       throwOnMissingCollections_(cfg.getParameter<bool>("throwOnMissingCollections")) {
   produces<ClusterTPAssociation>();
 }

ClusterTPAssociationProducer::~ClusterTPAssociationProducer ( )

override

Definition at line 78 of file ClusterTPAssociationProducer.cc.

78 {}

Member Function Documentation

void ClusterTPAssociationProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 80 of file ClusterTPAssociationProducer.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                               {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("simTrackSrc", edm::InputTag("g4SimHits"));
   desc.add<edm::InputTag>("pixelSimLinkSrc", edm::InputTag("simSiPixelDigis"));
   desc.add<edm::InputTag>("stripSimLinkSrc", edm::InputTag("simSiStripDigis"));
   desc.add<edm::InputTag>("phase2OTSimLinkSrc", edm::InputTag("simSiPixelDigis", "Tracker"));
   desc.add<edm::InputTag>("pixelClusterSrc", edm::InputTag("siPixelClusters"));
   desc.add<edm::InputTag>("stripClusterSrc", edm::InputTag("siStripClusters"));
   desc.add<edm::InputTag>("phase2OTClusterSrc", edm::InputTag("siPhase2Clusters"));
   desc.add<edm::InputTag>("trackingParticleSrc", edm::InputTag("mix", "MergedTrackTruth"));
   desc.add<bool>("throwOnMissingCollections", true);
   descriptions.add("tpClusterProducerDefault", desc);
 }

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 289 of file ClusterTPAssociationProducer.cc.

References edm::DetSet< T >::data, DEFINE_FWK_MODULE, 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::StreamID	,
		edm::Event &	iEvent,
		const edm::EventSetup &	es
	)		const

overrideprivate

Definition at line 94 of file ClusterTPAssociationProducer.cc.

References SiStripCluster::amplitudes(), edmNew::DetSet< T >::begin(), Phase2TrackerCluster1D::column(), runTauDisplay::eid, edmNew::DetSet< T >::end(), plotBeamSpotDB::first, Phase2TrackerCluster1D::firstRow(), SiStripCluster::firstStrip(), edm::Event::getByToken(), edm::HandleBase::id(), edm::HandleBase::isValid(), plotBeamSpotDB::last, edmNew::makeRefTo(), taus_updatedMVAIds_cff::mapping, SiPixelCluster::maxPixelCol(), SiPixelCluster::maxPixelRow(), SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), eostools::move(), seedClusterRemoverPhase2_cfi::phase2OTClusters, phase2OTClustersToken_, HiRegitMuonDetachedTripletStep_cff::pixelClusters, pixelClustersToken_, PixelChannelIdentifier::pixelToChannel(), Phase2TrackerDigi::pixelToChannel(), edm::Handle< T >::product(), edm::Event::put(), siphase2OTSimLinksToken_, sipixelSimLinksToken_, sistripSimLinksToken_, Phase2TrackerCluster1D::size(), HiRegitMuonDetachedTripletStep_cff::stripClusters, stripClustersToken_, throwOnMissingCollections_, and trackingParticleToken_.

                                                                                                        {
   // 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);
 
   // Phase2 Cluster
   edm::Handle<edmNew::DetSetVector<Phase2TrackerCluster1D> > phase2OTClusters;
   bool foundPhase2OTClusters = iEvent.getByToken(phase2OTClustersToken_, phase2OTClusters);
 
   // Pixel DigiSimLink
   edm::Handle<edm::DetSetVector<PixelDigiSimLink> > sipixelSimLinks;
   //  iEvent.getByLabel(_pixelSimLinkSrc, sipixelSimLinks);
   auto pixelSimLinksFound = iEvent.getByToken(sipixelSimLinksToken_, sipixelSimLinks);
   if (not throwOnMissingCollections_ and foundPixelClusters and not pixelSimLinksFound) {
     auto clusterTPList = std::make_unique<ClusterTPAssociation>();
     iEvent.put(std::move(clusterTPList));
     return;
   }
 
   // SiStrip DigiSimLink
   edm::Handle<edm::DetSetVector<StripDigiSimLink> > sistripSimLinks;
   auto stripSimLinksFound = iEvent.getByToken(sistripSimLinksToken_, sistripSimLinks);
   if (not throwOnMissingCollections_ and foundStripClusters and not stripSimLinksFound) {
     auto clusterTPList = std::make_unique<ClusterTPAssociation>();
     iEvent.put(std::move(clusterTPList));
     return;
   }
 
   // Phase2 OT DigiSimLink
   edm::Handle<edm::DetSetVector<PixelDigiSimLink> > siphase2OTSimLinks;
   auto phase2OTSimLinksFound = iEvent.getByToken(siphase2OTSimLinksToken_, siphase2OTSimLinks);
   if (not throwOnMissingCollections_ and foundPhase2OTClusters and not phase2OTSimLinksFound) {
     auto clusterTPList = std::make_unique<ClusterTPAssociation>();
     iEvent.put(std::move(clusterTPList));
     return;
   }
 
   // TrackingParticle
   edm::Handle<TrackingParticleCollection> TPCollectionH;
   auto tpFound = iEvent.getByToken(trackingParticleToken_, TPCollectionH);
   if (not throwOnMissingCollections_ and not tpFound) {
     auto clusterTPList = std::make_unique<ClusterTPAssociation>();
     iEvent.put(std::move(clusterTPList));
     return;
   }
 
   auto clusterTPList = std::make_unique<ClusterTPAssociation>(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
     clusterTPList->addKeyID(pixelClusters.id());
     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.empty())
               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->emplace_back(OmniClusterRef(c_ref), ipos->second);
           }
         }
       }
     }
   }
 
   if (foundStripClusters) {
     // Strip Clusters
     clusterTPList->addKeyID(stripClusters.id());
     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.empty())
             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->emplace_back(OmniClusterRef(c_ref), ipos->second);
           }
         }
       }
     }
   }
 
   if (foundPhase2OTClusters) {
     // Phase2 Clusters
     clusterTPList->addKeyID(phase2OTClusters.id());
     if (phase2OTClusters.isValid()) {
       for (edmNew::DetSetVector<Phase2TrackerCluster1D>::const_iterator iter = phase2OTClusters->begin(false),
                                                                         eter = phase2OTClusters->end(false);
            iter != eter;
            ++iter) {
         if (!(*iter).isValid())
           continue;
         uint32_t detid = iter->id();
         DetId detId(detid);
         edmNew::DetSet<Phase2TrackerCluster1D> link_phase2 = (*iter);
         for (edmNew::DetSet<Phase2TrackerCluster1D>::const_iterator di = link_phase2.begin(); di != link_phase2.end();
              di++) {
           const Phase2TrackerCluster1D& cluster = (*di);
           edm::Ref<edmNew::DetSetVector<Phase2TrackerCluster1D>, Phase2TrackerCluster1D> c_ref =
               edmNew::makeRefTo(phase2OTClusters, di);
 
           std::set<std::pair<uint32_t, EncodedEventId> > simTkIds;
 
           for (unsigned int istr(0); istr < cluster.size(); ++istr) {
             uint32_t channel = Phase2TrackerDigi::pixelToChannel(cluster.firstRow() + istr, cluster.column());
             std::vector<std::pair<uint32_t, EncodedEventId> > trkid(
                 getSimTrackId<PixelDigiSimLink>(siphase2OTSimLinks, detId, channel));
             if (trkid.empty())
               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()) {
               clusterTPList->emplace_back(OmniClusterRef(c_ref), ipos->second);
             }
           }
         }
       }
     }
   }
   clusterTPList->sortAndUnique();
   iEvent.put(std::move(clusterTPList));
 }

Member Data Documentation

edm::EDGetTokenT<edmNew::DetSetVector<Phase2TrackerCluster1D> > ClusterTPAssociationProducer::phase2OTClustersToken_

private

Definition at line 54 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

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

private

Definition at line 52 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

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

private

Definition at line 51 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

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

private

Definition at line 49 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

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

private

Definition at line 50 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

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

private

Definition at line 53 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

bool ClusterTPAssociationProducer::throwOnMissingCollections_

private

Definition at line 56 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

edm::EDGetTokenT<TrackingParticleCollection> ClusterTPAssociationProducer::trackingParticleToken_

private

Definition at line 55 of file ClusterTPAssociationProducer.cc.

Referenced by produce().

Private Attributes
edm::EDGetTokenT< edmNew::DetSetVector< Phase2TrackerCluster1D > >	phase2OTClustersToken_

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

edm::EDGetTokenT< edm::DetSetVector< PixelDigiSimLink > >	siphase2OTSimLinksToken_

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

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

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

bool	throwOnMissingCollections_

edm::EDGetTokenT< TrackingParticleCollection >	trackingParticleToken_

Public Types

Public Member Functions

Static 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