#include <L3TkMuonProducer.h>

Inheritance diagram for L3TkMuonProducer:

Public Types
typedef edm::Ref < L3MuonTrajectorySeedCollection >	SeedRef

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

Public Member Functions
	L3TkMuonProducer (const edm::ParameterSet &)
	constructor with config More...

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

virtual	~L3TkMuonProducer ()
	destructor More...

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef std::map< pseudoRef, std::pair< reco::TrackRef, SeedRef > >	LXtoL3sMap

typedef std::pair< unsigned int, unsigned int >	pseudoRef

Private Member Functions
pseudoRef	makePseudoRef (const L3MuonTrajectorySeed &s)

bool	sharedSeed (const L3MuonTrajectorySeed &s1, const L3MuonTrajectorySeed &s2)

Static Private Member Functions
static bool	seedRefBypT (const SeedRef &s1, const SeedRef &s2)

static bool	trackRefBypT (const reco::TrackRef &t1, const reco::TrackRef &t2)

Private Attributes
edm::InputTag	theL3CollectionLabel

edm::EDGetTokenT < reco::TrackCollection >	trackToken_

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

This module creates a skimed list of reco::Track (pointing to the original TrackExtra and TrackingRecHitOwnedVector One highest pT track per L1/L2 is selected, requiring some quality.

Author: J-R Vlimant.

Definition at line 25 of file L3TkMuonProducer.h.

Member Typedef Documentation

typedef std::map<pseudoRef, std::pair<reco::TrackRef,SeedRef> > L3TkMuonProducer::LXtoL3sMap

private

Definition at line 50 of file L3TkMuonProducer.h.

typedef std::pair<unsigned int,unsigned int> L3TkMuonProducer::pseudoRef

private

Definition at line 49 of file L3TkMuonProducer.h.

typedef edm::Ref<L3MuonTrajectorySeedCollection> L3TkMuonProducer::SeedRef

Definition at line 38 of file L3TkMuonProducer.h.

Constructor & Destructor Documentation

L3TkMuonProducer::L3TkMuonProducer ( const edm::ParameterSet & parameterSet )

constructor with config

Definition at line 27 of file L3TkMuonProducer.cc.

References edm::ParameterSet::getParameter(), LogTrace, edm::BranchDescription::moduleLabel(), edm::BranchDescription::processName(), edm::BranchDescription::productInstanceName(), and edm::BranchDescription::unwrappedTypeID().

                                                                   {
   LogTrace("Muon|RecoMuon|L3TkMuonProducer")<<" constructor called";
 
   // StandAlone Collection Label
   theL3CollectionLabel = parameterSet.getParameter<InputTag>("InputObjects");
   trackToken_ = consumes<reco::TrackCollection>(theL3CollectionLabel); 
   produces<TrackCollection>();
   produces<TrackExtraCollection>();
   produces<TrackingRecHitCollection>();
 
 
 
   callWhenNewProductsRegistered( [this](const edm::BranchDescription& iBD) {
                    edm::TypeID id(typeid(L3MuonTrajectorySeedCollection));
                    if(iBD.unwrappedTypeID() == id) {
                      this->mayConsume<L3MuonTrajectorySeedCollection>(edm::InputTag{iBD.moduleLabel(), iBD.productInstanceName(),iBD.processName()} );
                    }
                  });
 
 
 
 }

L3TkMuonProducer::~L3TkMuonProducer ( )

virtual

destructor

Definition at line 51 of file L3TkMuonProducer.cc.

References LogTrace.

                                    {
   LogTrace("Muon|RecoMuon|L3TkMuonProducer")<<" L3TkMuonProducer destructor called";
 }

Member Function Documentation

pseudoRef L3TkMuonProducer::makePseudoRef ( const L3MuonTrajectorySeed & s )

inlineprivate

Definition at line 52 of file L3TkMuonProducer.h.

References edm::ProductID::id(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isNull(), edm::Ref< C, T, F >::key(), L3MuonTrajectorySeed::l1Particle(), and L3MuonTrajectorySeed::l2Track().

                                                         {
     reco::TrackRef l2ref = s.l2Track();
     if (l2ref.isNull()){
       l1extra::L1MuonParticleRef l1ref = s.l1Particle();
       return std::make_pair(l1ref.id().id(),l1ref.key());
     }else return std::make_pair(l2ref.id().id(),l2ref.key());
   }

void L3TkMuonProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

virtual

produce candidates

reconstruct muons

Implements edm::EDProducer.

Definition at line 109 of file L3TkMuonProducer.cc.

References TrajectorySeed::direction(), f, i, edm::ProductID::id(), edm::Ref< C, T, F >::id(), reco::Track::innerDetId(), reco::Track::innerMomentum(), reco::Track::innerOk(), reco::Track::innerPosition(), reco::Track::innerStateCovariance(), edm::Ref< C, T, F >::key(), TSGFromL2_cfi::l3seeds(), LogDebug, findQualityFiles::maxI, metname, edm::Provenance::moduleLabel(), reco::Track::outerDetId(), reco::Track::outerMomentum(), reco::Track::outerOk(), reco::Track::outerPosition(), reco::Track::outerStateCovariance(), printseed(), printvector(), edm::Provenance::processName(), edm::Provenance::productInstanceName(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), python.multivaluedict::sort(), and testEve_cfg::tracks.

                                                                         {
   const string metname = "Muon|RecoMuon|L3TkMuonProducer";
   
   // Take the L3 container
   LogDebug(metname)<<" Taking the L3/GLB muons: "<<theL3CollectionLabel.label();
   Handle<TrackCollection> tracks; 
   event.getByToken(trackToken_,tracks);
 
   //make the LX->L3s pools
   LXtoL3sMap LXtoL3s;
 
   unsigned int maxI = tracks->size();
   bool gotL3seeds=false;
   edm::Handle<L3MuonTrajectorySeedCollection> l3seeds;
 
   //make a list of reference to tracker tracks
   vector<TrackRef> orderedTrackTracks(maxI);
   for (unsigned int i=0;i!=maxI;i++) orderedTrackTracks[i]=TrackRef(tracks,i);
   LogDebug(metname)<<"vector of L3 tracks before ordering:\n"<<printvector(orderedTrackTracks);
   //order them in pT
   sort(orderedTrackTracks.begin(),orderedTrackTracks.end(),trackRefBypT);
   LogDebug(metname)<<"vector of L3 tracks after ordering:\n"<<printvector(orderedTrackTracks);
   //loop over then
   for (unsigned int i=0;i!=maxI;i++) {
     TrackRef & tk=orderedTrackTracks[i];  
     SeedRef l3seedRef = tk->seedRef().castTo<SeedRef>();
 
     vector<SeedRef> allPossibleOrderedLx; // with identical hit-set
     //add the direct relation
     allPossibleOrderedLx.push_back(l3seedRef);
     LogDebug(metname)<<"adding the seed ref: "<<l3seedRef.id().id()<<":"<<l3seedRef.key()<<" for this tracker track: "<<tk.id().id()<<":"<<tk.key();
 
     //add the relations due to shared seeds
     //check whether there is a "shared" seed in addition
     if (!gotL3seeds){
       //need to fetch the handle from the ref
       const edm::Provenance & seedsProv=event.getProvenance(l3seedRef.id());
       edm::InputTag l3seedsTag(seedsProv.moduleLabel(), seedsProv.productInstanceName(), seedsProv.processName());
       event.getByLabel(l3seedsTag, l3seeds);
       gotL3seeds=true;
       LogDebug(metname)<<"got seeds handle from: "<<l3seedsTag;
     }
     //loop the other seeds in the collection
     for (unsigned int iS=0;iS!=l3seeds->size();++iS){
       const L3MuonTrajectorySeed & seed = (*l3seeds)[iS];
       const L3MuonTrajectorySeed & thisSeed = *l3seedRef;
       if (l3seedRef.key()==iS) continue; //not care about this one
       //compare this seed with the seed in the collection
       if (sharedSeed(seed,thisSeed)){
     SeedRef thisSharedSeedRef(l3seeds,iS);
     LogDebug(metname)<<"shared seeds: \n"<<printseed(l3seedRef)<<" and: \n"<<printseed(thisSharedSeedRef)
              <<"\nadding ANOTHER seed ref: "<<thisSharedSeedRef.id().id()<<":"<<thisSharedSeedRef.key()<<" for this tracker track: "<<tk.id().id()<<":"<<tk.key();
     //  edm::LogError(metname)<<" we have a shared seed right there.";
     allPossibleOrderedLx.push_back(thisSharedSeedRef);
       }//seed is shared
     }//loop all other existing seed for overlaps
 
     //now you have the full list of Lx objects that have seed this tracker track.
     // order the list in pT of Lx objects
     LogDebug(metname)<<"list of possible Lx objects for tracker track: "<<tk.id().id()<<":"<<tk.key()<<" before ordering\n"<<printvector(allPossibleOrderedLx);
     sort(allPossibleOrderedLx.begin(),allPossibleOrderedLx.end(),seedRefBypT);
     LogDebug(metname)<<"list of possible Lx objects for tracker track: "<<tk.id().id()<<":"<<tk.key()<<" after ordering\n"<<printvector(allPossibleOrderedLx);
     // assign this tracker track to the highest pT Lx.
     for (unsigned int iL=0;iL!=allPossibleOrderedLx.size();++iL){
       SeedRef thisRef=allPossibleOrderedLx[iL];
       pseudoRef ref = makePseudoRef(*thisRef);
       LogDebug(metname)<<"seed ref: "<<thisRef.id().id()<<":"<<thisRef.key()<<" transcribe to pseudoref: "<<ref.first<<":"<<ref.second;
       LXtoL3sMap::iterator f=LXtoL3s.find(ref);
       if (f!=LXtoL3s.end()){
     //there's already an entry. because of the prior ordering in pT of the tracker track refs
     // the track ref already there *has* a higher pT: this one cannot compete and should be assigned to the next Lx;
     LogDebug(metname)<<"this tracker track: "<<tk.id().id()<<":"<<tk.key()<<" ("<< tk->pt()<<")"
              <<"\n cannot compete in pT with track: "<<f->second.first.id().id()<<":"<<f->second.first.key()<<" ("<<f->second.first->pt()<<")"
              <<"\n already assigned to pseudo ref: "<<ref.first<<":"<<ref.second<<" which corresponds to seedRef: "<<f->second.second.id().id()<<":"<<f->second.second.key();
     continue;
       }else{
     //there was no entry yet. make the assignement
     LogDebug(metname)<<"this tracker track: "<<tk.id().id()<<":"<<tk.key()
              <<" is assigned to pseudo ref: "<<ref.first<<":"<<ref.second<<" which corresponds to seedRef: "<<thisRef.id().id()<<":"<<thisRef.key();
     LXtoL3s[ref] = std::make_pair(tk,thisRef);
     //once assigned. break
     break;
       }
     }//loop possible Lx for possible assignement
   }//loop over ordered list of tracker track refs
 
 
   //prepare the output
   std::auto_ptr<TrackCollection> outTracks( new TrackCollection(LXtoL3s.size()));
   std::auto_ptr<TrackExtraCollection> outTrackExtras( new TrackExtraCollection(LXtoL3s.size()));
   reco::TrackExtraRefProd rTrackExtras = event.getRefBeforePut<TrackExtraCollection>();
   std::auto_ptr<TrackingRecHitCollection> outRecHits( new TrackingRecHitCollection());
   TrackingRecHitRefProd rHits = event.getRefBeforePut<TrackingRecHitCollection>();
 
   LogDebug(metname)<<"reading the map to make "<< LXtoL3s.size()<<"products.";
   //fill the collection from the map
   LXtoL3sMap::iterator f=LXtoL3s.begin();
   unsigned int i=0;
   for (;f!=LXtoL3s.end();++f,++i){
 
     LogDebug(metname)<<"copy the track over, and make ref to extra";
     const Track & trk = *(f->second.first);
     (*outTracks)[i] = Track(trk);
     (*outTracks)[i].setExtra( TrackExtraRef(rTrackExtras,i));
 
     LogDebug(metname)<<"copy the trackExtra too, and change the seedref";
     edm::RefToBase<TrajectorySeed> seedRef(f->second.second);
     //do not use the copy constructor, otherwise the hit Ref are still the same
     (*outTrackExtras)[i] = TrackExtra(
                       trk.outerPosition(), trk.outerMomentum(), trk.outerOk(),
                       trk.innerPosition(), trk.innerMomentum(), trk.innerOk(),
                       trk.outerStateCovariance(), trk.outerDetId(),
                       trk.innerStateCovariance(), trk.innerDetId(),
                       seedRef->direction(),seedRef
                       );
 
     LogDebug(metname)<<"copy the hits too";
     unsigned int iRH=0;
     for( trackingRecHit_iterator hit = trk.recHitsBegin(); hit != trk.recHitsEnd(); ++ hit,++iRH ) {
       outRecHits->push_back((*hit)->clone());
       (*outTrackExtras)[i].add( TrackingRecHitRef( rHits, iRH));
     }
   }
   
   LogDebug(metname)<<"made: "<<outTracks->size()<<" tracks, "<<outTrackExtras->size()<<" extras and "<<outRecHits->size()<<" rechits.";
 
   //put the collection in the event
   LogDebug(metname)<<"loading...";
   event.put(outTracks);
   event.put(outTrackExtras);
   event.put(outRecHits);
   LogDebug(metname)<<" Event loaded"
            <<"================================";
 }

static bool L3TkMuonProducer::seedRefBypT	(	const SeedRef &	s1,
		const SeedRef &	s2
	)

inlinestaticprivate

Definition at line 64 of file L3TkMuonProducer.h.

References edm::Ref< C, T, F >::isNull().

                                                                   {
     double pt1,pt2;
     reco::TrackRef l2ref1 = s1->l2Track();
     if (l2ref1.isNull()) pt1=s1->l1Particle()->pt();
     else pt1=l2ref1->pt();
     reco::TrackRef l2ref2 = s2->l2Track();
     if (l2ref2.isNull()) pt2=s2->l1Particle()->pt();
     else pt2=l2ref2->pt();
     return (pt1>pt2);
   }

bool L3TkMuonProducer::sharedSeed	(	const L3MuonTrajectorySeed &	s1,
		const L3MuonTrajectorySeed &	s2
	)

private

Definition at line 55 of file L3TkMuonProducer.cc.

References TrackingRecHit::all, TrajectorySeed::nHits(), diffTwoXMLs::r1, diffTwoXMLs::r2, and TrajectorySeed::recHits().

                                                                                               {
   //quit right away on nH=0
   if (s1.nHits()==0 || s2.nHits()==0) return false;
   //quit right away if not the same number of hits
   if (s1.nHits()!=s2.nHits()) return false;
   TrajectorySeed::range r1=s1.recHits();
   TrajectorySeed::range r2=s2.recHits();
   TrajectorySeed::const_iterator i1,i2;
   TrajectorySeed::const_iterator & i1_e=r1.second,&i2_e=r2.second;
   TrajectorySeed::const_iterator & i1_b=r1.first,&i2_b=r2.first;
   //quit right away if first detId does not match. front exist because of ==0 ->quit test
   if(i1_b->geographicalId() != i2_b->geographicalId()) return false;
   //then check hit by hit if they are the same
   for (i1=i1_b,i2=i2_b;i1!=i1_e,i2!=i2_e;++i1,++i2){
     if (!i1->sharesInput(&(*i2),TrackingRecHit::all)) return false;
   }
   return true;
 }

static bool L3TkMuonProducer::trackRefBypT	(	const reco::TrackRef &	t1,
		const reco::TrackRef &	t2
	)

inlinestaticprivate

Definition at line 75 of file L3TkMuonProducer.h.

                                                                            {
     return (t1->pt()>t2->pt());
   }

Member Data Documentation

edm::InputTag L3TkMuonProducer::theL3CollectionLabel

private

Definition at line 43 of file L3TkMuonProducer.h.

edm::EDGetTokenT<reco::TrackCollection> L3TkMuonProducer::trackToken_

private

Definition at line 44 of file L3TkMuonProducer.h.

Public Types

Public Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation