Inheritance diagram for MuonsFromRefitTracksProducer:

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

	~MuonsFromRefitTracksProducer ()

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 Member Functions
virtual void	beginJob ()

reco::Muon *	cloneAndSwitchTrack (const reco::Muon &muon, const reco::TrackRef &newTrack) const

virtual void	endJob ()

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

bool	storeMatchMaps (const edm::Event &event)

Private Attributes
bool	fromCocktail

bool	fromGlobalTrack

bool	fromSigmaSwitch

bool	fromTeVRefit

bool	fromTMR

bool	fromTrackerTrack

double	nSigmaSwitch

double	ptThreshold

edm::InputTag	src

std::string	tevMuonTracks

double	TMRcut

edm::Handle < reco::TrackToTrackMap >	trackMap

edm::Handle < reco::TrackToTrackMap >	trackMapDefault

edm::Handle < reco::TrackToTrackMap >	trackMapFirstHit

edm::Handle < reco::TrackToTrackMap >	trackMapPicky

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 67 of file MuonsFromRefitTracksProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 143 of file MuonsFromRefitTracksProducer.cc.

References fromTeVRefit, and tevMuonTracks.

   : src(cfg.getParameter<edm::InputTag>("src")),
     fromTrackerTrack(cfg.getParameter<bool>("fromTrackerTrack")),
     fromGlobalTrack(cfg.getParameter<bool>("fromGlobalTrack")),
     tevMuonTracks(cfg.getParameter<std::string>("tevMuonTracks")),
     fromCocktail(cfg.getParameter<bool>("fromCocktail")),
     fromTMR(cfg.getParameter<bool>("fromTMR")),
     TMRcut(cfg.getParameter<double>("TMRcut")),
     fromSigmaSwitch(cfg.getParameter<bool>("fromSigmaSwitch")),
     nSigmaSwitch(cfg.getParameter<double>("nSigmaSwitch")),
     ptThreshold(cfg.getParameter<double>("ptThreshold"))
 {
   fromTeVRefit = tevMuonTracks != "none";
   produces<reco::MuonCollection>();
 }

MuonsFromRefitTracksProducer::~MuonsFromRefitTracksProducer ( )

inline

Definition at line 70 of file MuonsFromRefitTracksProducer.cc.

70 {}

Member Function Documentation

virtual void MuonsFromRefitTracksProducer::beginJob ( void )

inlineprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 73 of file MuonsFromRefitTracksProducer.cc.

73 {}

reco::Muon * MuonsFromRefitTracksProducer::cloneAndSwitchTrack	(	const reco::Muon &	muon,
		const reco::TrackRef &	newTrack
	)		const

private

Definition at line 173 of file MuonsFromRefitTracksProducer.cc.

References reco::Muon::clone(), reco::Muon::innerTrack(), reco::Muon::outerTrack(), L1TEmulatorMonitor_cff::p, p4, reco::LeafCandidate::setCharge(), reco::Muon::setGlobalTrack(), reco::Muon::setInnerTrack(), reco::Muon::setOuterTrack(), reco::LeafCandidate::setP4(), reco::LeafCandidate::setVertex(), and mathSSE::sqrt().

Referenced by produce().

                                                                 {
   // Muon mass to make a four-vector out of the new track.
   static const double muMass = 0.10566;
 
   reco::TrackRef tkTrack  = muon.innerTrack();
   reco::TrackRef muTrack  = muon.outerTrack();
       
   // Make up a real Muon from the tracker track.
   reco::Particle::Point vtx(newTrack->vx(), newTrack->vy(), newTrack->vz());
   reco::Particle::LorentzVector p4;
   double p = newTrack->p();
   p4.SetXYZT(newTrack->px(), newTrack->py(), newTrack->pz(),
          sqrt(p*p + muMass*muMass));
 
   reco::Muon* mu = muon.clone();
   mu->setCharge(newTrack->charge());
   mu->setP4(p4);
   mu->setVertex(vtx);
   mu->setGlobalTrack(newTrack);
   mu->setInnerTrack(tkTrack);
   mu->setOuterTrack(muTrack);
   return mu;
 }

virtual void MuonsFromRefitTracksProducer::endJob ( void )

inlineprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 75 of file MuonsFromRefitTracksProducer.cc.

75 {}

void MuonsFromRefitTracksProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	eSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 198 of file MuonsFromRefitTracksProducer.cc.

References edm::ViewBase::clone(), cloneAndSwitchTrack(), edm::HandleBase::failedToGet(), fromCocktail, fromGlobalTrack, fromSigmaSwitch, fromTeVRefit, fromTMR, fromTrackerTrack, edm::Ref< C, T, F >::isNonnull(), prof2calltree::last, metsig::muon, patZpeak::muons, nSigmaSwitch, convertSQLiteXML::ok, ptThreshold, reco::Muon::setGlobalTrack(), reco::Muon::setInnerTrack(), reco::Muon::setOuterTrack(), sigmaSwitch(), src, storeMatchMaps(), tevMuonTracks, muon::tevOptimized(), tevOptimizedTMR(), TMRcut, trackMap, trackMapDefault, trackMapFirstHit, trackMapPicky, and edm::helpers::KeyVal< K, V >::val.

                                                                                        {
   // Get the global muons from the event.
   edm::Handle<edm::View<reco::Muon> > muons;
   event.getByLabel(src, muons);
 
   // If we can't get the global muon collection, or below the
   // track-to-track maps needed, still produce an empty collection of
   // muons so consumers don't throw an exception.
   bool ok = !muons.failedToGet();
 
   // If we're instructed to use the TeV refit tracks in some way, we
   // need the track-to-track maps. If we're making a cocktail muon,
   // get all three track maps (the cocktail ingredients); else just
   // get the map which takes the above global tracks to the desired
   // TeV-muon refitted tracks (firstHit or picky).
   if (ok && fromTeVRefit)
     ok = storeMatchMaps(event);
 
   // Make the output collection.
   std::auto_ptr<reco::MuonCollection> cands(new reco::MuonCollection);
 
   if (ok) {
     edm::View<reco::Muon>::const_iterator muon;
     for (muon = muons->begin(); muon != muons->end(); muon++) {
       // Filter out the so-called trackerMuons and stand-alone muons
       // (and caloMuons, if they were ever to get into the input muons
       // collection).
       if (!muon->isGlobalMuon()) continue;
 
       if (fromTeVRefit || fromSigmaSwitch) {
     // Start out with a null TrackRef.
     reco::TrackRef tevTk;
       
     // If making a cocktail muon, use tevOptimized() to get the track
     // desired. Otherwise, get the refit track from the desired track
     // map.
     if (fromTMR)
       tevTk = tevOptimizedTMR(*muon, *trackMapFirstHit, TMRcut);
     else if (fromCocktail)
       tevTk = muon::tevOptimized(*muon, *trackMapDefault, *trackMapFirstHit,
                      *trackMapPicky);
     else if (fromSigmaSwitch)
       tevTk = sigmaSwitch(*muon, nSigmaSwitch, ptThreshold);
     else {
         reco::TrackToTrackMap::const_iterator tevTkRef =
         trackMap->find(muon->combinedMuon());
       if (tevTkRef != trackMap->end())
         tevTk = tevTkRef->val;
     }
     
     // If the TrackRef is valid, make a new Muon that has the same
     // tracker and stand-alone tracks, but has the refit track as
     // its global track.
     if (tevTk.isNonnull())
       cands->push_back(*cloneAndSwitchTrack(*muon, tevTk));
       }
       else if (fromTrackerTrack)
     cands->push_back(*cloneAndSwitchTrack(*muon, muon->innerTrack()));
       else if (fromGlobalTrack)
     cands->push_back(*cloneAndSwitchTrack(*muon, muon->globalTrack()));
       else {
     cands->push_back(*muon->clone());
 
     // Just cloning does not work in the case of the source being
     // a pat::Muon with embedded track references -- these do not
     // get copied. Explicitly set them.
     reco::Muon& last = cands->at(cands->size()-1);
     if (muon->globalTrack().isTransient())
       last.setGlobalTrack(muon->globalTrack());
     if (muon->innerTrack().isTransient())
       last.setInnerTrack(muon->innerTrack());
     if (muon->outerTrack().isTransient())
       last.setOuterTrack(muon->outerTrack());
       }
     }
   }
   else
     edm::LogWarning("MuonsFromRefitTracksProducer")
       << "either " << src << " or the track map(s) " << tevMuonTracks
       << " not present in the event; producing empty collection";
   
   event.put(cands);
 }

bool MuonsFromRefitTracksProducer::storeMatchMaps ( const edm::Event & event )

private

Definition at line 159 of file MuonsFromRefitTracksProducer.cc.

References edm::HandleBase::failedToGet(), fromCocktail, fromTMR, tevMuonTracks, trackMap, trackMapDefault, trackMapFirstHit, and trackMapPicky.

Referenced by produce().

                                                                        {
   if (fromCocktail || fromTMR) {
     event.getByLabel(tevMuonTracks, "default",  trackMapDefault);
     event.getByLabel(tevMuonTracks, "firstHit", trackMapFirstHit);
     event.getByLabel(tevMuonTracks, "picky",    trackMapPicky);
     return !trackMapDefault.failedToGet() && 
       !trackMapFirstHit.failedToGet() && !trackMapPicky.failedToGet();
   }
   else {
     event.getByLabel(edm::InputTag(tevMuonTracks), trackMap);
     return !trackMap.failedToGet();
   }
 }