#include <MuonLinksProducerForHLT.h>

Inheritance diagram for MuonLinksProducerForHLT:

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

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

virtual	~MuonLinksProducerForHLT ()

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

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

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

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::EDGetTokenT< reco::MuonTrackLinksCollection >	linkToken_

double	pMin

double	ptMin

double	shareHitFraction

edm::InputTag	theInclusiveTrackCollectionInInput

edm::InputTag	theLinkCollectionInInput

edm::EDGetTokenT< reco::TrackCollection >	trackToken_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::global::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

Author: R. Bellan - UCSB ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Simple producer to make reco::MuonTrackLinks collection out of the global muons from "muons" collection to restore dropped links used as input for MuonIdProducer.

Author: R. Bellan - UCSB ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 25 of file MuonLinksProducerForHLT.h.

Constructor & Destructor Documentation

MuonLinksProducerForHLT::MuonLinksProducerForHLT ( const edm::ParameterSet & iConfig )

explicit

Definition at line 23 of file MuonLinksProducerForHLT.cc.

References edm::ParameterSet::getParameter(), linkToken_, pMin, ptMin, shareHitFraction, theInclusiveTrackCollectionInInput, theLinkCollectionInInput, and trackToken_.

 {
    produces<reco::MuonTrackLinksCollection>();
    theLinkCollectionInInput = iConfig.getParameter<edm::InputTag>("LinkCollection");
    theInclusiveTrackCollectionInInput = iConfig.getParameter<edm::InputTag>("InclusiveTrackerTrackCollection");
    ptMin = iConfig.getParameter<double>("ptMin");
    pMin = iConfig.getParameter<double>("pMin");
    shareHitFraction = iConfig.getParameter<double>("shareHitFraction");
 
    linkToken_ = consumes<reco::MuonTrackLinksCollection>(theLinkCollectionInInput);
    trackToken_ = consumes<reco::TrackCollection>(theInclusiveTrackCollectionInInput);
    
 }

MuonLinksProducerForHLT::~MuonLinksProducerForHLT ( )

virtual

Definition at line 37 of file MuonLinksProducerForHLT.cc.

38 {

39 }

Member Function Documentation

void MuonLinksProducerForHLT::produce	(	edm::StreamID	,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

overridevirtual

Definition at line 41 of file MuonLinksProducerForHLT.cc.

References runEdmFileComparison::found, dedxEstimators_cff::fraction, edm::Event::getByToken(), linkToken_, eostools::move(), convertSQLitetoXML_cfg::output, pMin, ptMin, edm::Event::put(), shareHitFraction, TrackingRecHit::some, HiIsolationCommonParameters_cff::track, and trackToken_.

 {
    auto output = std::make_unique<reco::MuonTrackLinksCollection>();
 
    edm::Handle<reco::MuonTrackLinksCollection> links; 
    iEvent.getByToken(linkToken_, links);
 
    edm::Handle<reco::TrackCollection> incTracks; 
    iEvent.getByToken(trackToken_, incTracks);
 
    for(reco::MuonTrackLinksCollection::const_iterator link = links->begin(); 
        link != links->end(); ++link){
      bool found = false;
      unsigned int trackIndex = 0;
      unsigned int muonTrackHits = link->trackerTrack()->extra()->recHitsSize();
      for(reco::TrackCollection::const_iterator track = incTracks->begin();
      track != incTracks->end(); ++track, ++trackIndex){      
        if ( track->pt() < ptMin ) continue;
        if ( track->p() < pMin ) continue;
        //std::cout << "pt (muon/track) " << link->trackerTrack()->pt() << " " << track->pt() << std::endl;
        unsigned trackHits = track->extra()->recHitsSize();
        //std::cout << "hits (muon/track) " << muonTrackHits  << " " << trackHits() << std::endl;
        unsigned int smallestNumberOfHits = trackHits < muonTrackHits ? trackHits : muonTrackHits;
        int numberOfCommonDetIds = 0;
        for ( auto hit = track->extra()->recHitsBegin();
          hit != track->extra()->recHitsEnd(); ++hit ) {
      for ( auto mit = link->trackerTrack()->extra()->recHitsBegin();
          mit != link->trackerTrack()->extra()->recHitsEnd(); ++mit ) {
        if ( (*hit)->geographicalId() == (*mit)->geographicalId() && 
         (*hit)->sharesInput((*mit),TrackingRecHit::some) ) { 
          numberOfCommonDetIds++;
          break;
        }
      }
        }
        double fraction = (double)numberOfCommonDetIds/smallestNumberOfHits;
        // std::cout << "Overlap/Smallest/fraction = " << numberOfCommonDetIds << " " << smallestNumberOfHits << " " << fraction << std::endl;
        if( fraction > shareHitFraction ) { 
      output->push_back(reco::MuonTrackLinks(reco::TrackRef(incTracks,trackIndex), 
                         link->standAloneTrack(), 
                         link->globalTrack() ) );
      found = true;
      break;
        }
      }
      if (!found) 
        output->push_back(reco::MuonTrackLinks(link->trackerTrack(), 
                           link->standAloneTrack(), 
                           link->globalTrack() ) );
    }  
    iEvent.put(std::move(output));
 }