#include <MuonLinksProducerForHLT.h>

Inheritance diagram for MuonLinksProducerForHLT:

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

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

virtual	~MuonLinksProducerForHLT ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
double	pMin

double	ptMin

double	shareHitFraction

edm::InputTag	theInclusiveTrackCollectionInInput

edm::InputTag	theLinkCollectionInInput

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Date:: 2011/05/02 16:09:31

Revision:: 1.1

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.

Date:: 2011/05/02 16:09:31

Revision:: 1.1

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

Definition at line 22 of file MuonLinksProducerForHLT.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 31 of file MuonLinksProducerForHLT.cc.

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

 {
    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");
 }

MuonLinksProducerForHLT::~MuonLinksProducerForHLT ( )

virtual

Definition at line 41 of file MuonLinksProducerForHLT.cc.

42 {

43 }

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 45 of file MuonLinksProducerForHLT.cc.

References newFWLiteAna::found, edm::Event::getByLabel(), link(), convertSQLitetoXML_cfg::output, pMin, ptMin, edm::Event::put(), shareHitFraction, TrackingRecHit::some, theInclusiveTrackCollectionInInput, and theLinkCollectionInInput.

 {
    std::auto_ptr<reco::MuonTrackLinksCollection> output(new reco::MuonTrackLinksCollection());
 
    edm::Handle<reco::MuonTrackLinksCollection> links; 
    iEvent.getByLabel(theLinkCollectionInInput, links);
 
    edm::Handle<reco::TrackCollection> incTracks; 
    iEvent.getByLabel(theInclusiveTrackCollectionInInput, 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()->recHits().size();
      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()->recHits().size();
        //std::cout << "hits (muon/track) " << muonTrackHits  << " " << trackHits() << std::endl;
        unsigned int smallestNumberOfHits = trackHits < muonTrackHits ? trackHits : muonTrackHits;
        int numberOfCommonDetIds = 0;
        for ( TrackingRecHitRefVector::const_iterator hit = track->extra()->recHitsBegin();
          hit != track->extra()->recHitsEnd(); ++hit ) {
      for ( TrackingRecHitRefVector::const_iterator mit = link->trackerTrack()->extra()->recHitsBegin();
          mit != link->trackerTrack()->extra()->recHitsEnd(); ++mit ) {
        if ( hit->get()->geographicalId() == mit->get()->geographicalId() && 
         hit->get()->sharesInput(mit->get(),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( output );
 }