Inheritance diagram for TrajectorySeedFromMuonProducer:

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

	TrajectorySeedFromMuonProducer (const edm::ParameterSet &)

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
edm::InputTag	muonCollectionTag_

bool	skipMatchedMuons_

edm::InputTag	trackCollectionTag_

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

Definition at line 48 of file TrajectorySeedFromMuonProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 61 of file TrajectorySeedFromMuonProducer.cc.

References edm::ParameterSet::getParameter(), muonCollectionTag_, skipMatchedMuons_, and trackCollectionTag_.

 {
   muonCollectionTag_  = iConfig.getParameter<edm::InputTag>("muonCollectionTag");
   trackCollectionTag_ = iConfig.getParameter<edm::InputTag>("trackCollectionTag");
   skipMatchedMuons_   = iConfig.getParameter<bool>("skipMatchedMuons");
   produces<TrajectorySeedCollection>();
 }

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 70 of file TrajectorySeedFromMuonProducer.cc.

References DeDxDiscriminatorTools::charge(), TrajectorySeed::direction(), muonid::findOppositeTrack(), edm::EventSetup::get(), edm::Event::getByLabel(), edm::Ref< C, T, F >::isNonnull(), LogTrace, metsig::muon, muonCollectionTag_, oppositeToMomentum, edm::OwnVector< T, P >::push_back(), edm::Event::put(), DetId::rawId(), dt_dqm_sourceclient_common_cff::reco, query::result, skipMatchedMuons_, and trackCollectionTag_.

 {
   using namespace edm;
   using namespace reco;
   using namespace std;
 
   // Product
   std::auto_ptr<TrajectorySeedCollection> result(new TrajectorySeedCollection());
   
   edm::ESHandle<MagneticField> magneticField;
   iSetup.get<IdealMagneticFieldRecord>().get(magneticField);
   
   edm::ESHandle<TrackerGeometry> trackerGeometry;
   iSetup.get<TrackerDigiGeometryRecord>().get(trackerGeometry);
   
   TrajectoryStateTransform transformer;
   
   edm::Handle<edm::View<Muon> > muonCollectionHandle;
   iEvent.getByLabel(muonCollectionTag_, muonCollectionHandle);
   
   edm::Handle<reco::TrackCollection> trackCollectionHandle;
   iEvent.getByLabel(trackCollectionTag_, trackCollectionHandle);
 
   // Loop over the muon track
   for (edm::View<Muon>::const_iterator muon = muonCollectionHandle->begin();  muon != muonCollectionHandle->end(); ++muon) {
     // muon must have a tracker track
     if( muon->innerTrack().isNull() ) continue;
     edm::RefToBase<reco::Track> track(muon->innerTrack());
     // check if there is a back-to-back track
     if( skipMatchedMuons_ &&
     muonid::findOppositeTrack(trackCollectionHandle, *track).isNonnull() ) continue;
     if( (!track->innerOk()) || (!track->recHit(0)->isValid())) continue;
     GlobalPoint innerPosition(track->innerPosition().x(), track->innerPosition().y(), track->innerPosition().z());
     GlobalVector innerMomentum(track->innerMomentum().x(), track->innerMomentum().y(), track->innerMomentum().z());
     int charge = track->charge();
     const reco::Track::CovarianceMatrix innerStateCovariance = track->innerStateCovariance();
     DetId innerDetId = DetId(track->innerDetId());
     // Construct the information necessary to make a TrajectoryStateOnSurface
     GlobalTrajectoryParameters globalTrajParams(innerPosition, innerMomentum, charge, &(*magneticField));
     CurvilinearTrajectoryError curviError(innerStateCovariance);
     FreeTrajectoryState tracker_state(globalTrajParams, curviError);
     LogTrace("MuonIdentification") << "Track Inner FTS: "<<tracker_state;
     
     TrajectoryStateOnSurface tracker_tsos  = TrajectoryStateOnSurface(globalTrajParams, curviError, trackerGeometry->idToDet(innerDetId)->surface());
     
     // Make Hits, push back the innermost Hit
     edm::OwnVector<TrackingRecHit> trackHits;
     trackHits.push_back(track->recHit(0)->clone() );
     
     // Make TrajectorySeed
     PTrajectoryStateOnDet *PTraj = transformer.persistentState(tracker_tsos, innerDetId.rawId());
     TrajectorySeed trajectorySeed(*PTraj, trackHits, oppositeToMomentum);
     LogTrace("MuonIdentification")<<"Trajectory Seed Direction: "<< trajectorySeed.direction()<<endl;
     result->push_back(trajectorySeed);
   }
   
   iEvent.put(result);
 
 }