Produces reco::GenParticle from SimTracks. More...

Inheritance diagram for pat::PATGenCandsFromSimTracksProducer:

Classes
struct	GlobalContext
	Global context for all recursive methods. More...

struct	LessById

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

	~PATGenCandsFromSimTracksProducer ()

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

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef StringCutObjectSelector < reco::GenParticle >	StrFilter

Private Member Functions
virtual void	endJob ()

const SimTrack *	findGeantMother (const SimTrack &tk, const GlobalContext &g) const
	Find the mother of a given GEANT track (or NULL if it can't be found). More...

edm::Ref < reco::GenParticleCollection >	findRef (const SimTrack &tk, GlobalContext &g) const

edm::Ref < reco::GenParticleCollection >	generatorRef_ (const SimTrack &tk, const GlobalContext &g) const
	Used by findRef if the track is a PYTHIA particle. More...

reco::GenParticle	makeGenParticle_ (const SimTrack &tk, const edm::Ref< reco::GenParticleCollection > &mother, const GlobalContext &g) const
	Make a GenParticle for this SimTrack, with a given mother. More...

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

Private Attributes
std::auto_ptr< StrFilter >	filter_

bool	firstEvent_

edm::InputTag	genParticles_
	Collection of GenParticles I need to make refs to. It must also have its associated vector<int> of barcodes, aligned with them. More...

bool	makeMotherLink_
	If true, I'll try to make a link from the GEANT particle to a GenParticle. More...

std::set< int >	motherPdgIds_

std::vector< PdtEntry >	motherPdts_

std::set< int >	pdgIds_

std::vector< PdtEntry >	pdts_

int	setStatus_

edm::InputTag	src_

bool	writeAncestors_
	If true, I'll save GenParticles corresponding to the ancestors of this GEANT particle. Common ancestors are only written once. More...

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)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

Produces reco::GenParticle from SimTracks.

The PATGenCandsFromSimTracksProducer produces GenParticles from SimTracks, so they can be used for MC matching.

Author: Jordan Tucker (original module), Giovanni Petrucciani (PAT integration)

Version

Id:: PATGenCandsFromSimTracksProducer.cc,v 1.9 2013/02/27 23:26:56 wmtan Exp

Definition at line 29 of file PATGenCandsFromSimTracksProducer.cc.

Member Typedef Documentation

typedef StringCutObjectSelector<reco::GenParticle> pat::PATGenCandsFromSimTracksProducer::StrFilter

private

Definition at line 46 of file PATGenCandsFromSimTracksProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 116 of file PATGenCandsFromSimTracksProducer.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::existsAs(), alcazmumu_cfi::filter, filter_, edm::ParameterSet::getParameter(), makeMotherLink_, motherPdts_, pdts_, and writeAncestors_.

                                                                                           :
   firstEvent_(true),
   src_(cfg.getParameter<InputTag>("src")),            // source sim tracks & vertices
   setStatus_(cfg.getParameter<int32_t>("setStatus")), // set status of GenParticle to this code
   makeMotherLink_(cfg.existsAs<bool>("makeMotherLink") ? cfg.getParameter<bool>("makeMotherLink") : false),
   writeAncestors_(cfg.existsAs<bool>("writeAncestors") ? cfg.getParameter<bool>("writeAncestors") : false),
   genParticles_(makeMotherLink_ ? cfg.getParameter<InputTag>("genParticles") : edm::InputTag())
 {
     // Possibly allow a list of particle types
     if (cfg.exists("particleTypes")) {
         pdts_ = cfg.getParameter<vector<PdtEntry> >("particleTypes");
     }
     if (cfg.exists("motherTypes")) {
         motherPdts_ = cfg.getParameter<vector<PdtEntry> >("motherTypes");
     }
 
     // Possibly allow a string cut
     if (cfg.existsAs<string>("filter")) {
         string filter = cfg.getParameter<string>("filter");
         if (!filter.empty()) {
             filter_ = auto_ptr<StrFilter>(new StrFilter(filter));
         }
     }
 
     if (writeAncestors_ && !makeMotherLink_) {
         edm::LogWarning("Configuration") << "PATGenCandsFromSimTracksProducer: " << 
             "you have set 'writeAncestors' to 'true' and 'makeMotherLink' to false;" <<
             "GEANT particles with generator level (e.g.PYHIA) mothers won't have mother links.\n";
     }
     produces<GenParticleCollection>();
 }

pat::PATGenCandsFromSimTracksProducer::~PATGenCandsFromSimTracksProducer ( )

inline

Definition at line 32 of file PATGenCandsFromSimTracksProducer.cc.

32 {}

Member Function Documentation

virtual void pat::PATGenCandsFromSimTracksProducer::endJob ( void )

inlineprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 36 of file PATGenCandsFromSimTracksProducer.cc.

36 {}

const SimTrack * PATGenCandsFromSimTracksProducer::findGeantMother	(	const SimTrack &	tk,
		const GlobalContext &	g
	)		const

private

Find the mother of a given GEANT track (or NULL if it can't be found).

Definition at line 149 of file PATGenCandsFromSimTracksProducer.cc.

References SimTrack::genpartIndex(), customizeTrackingMonitorSeedNumber::idx, SimVertex::noParent(), SimTrack::noVertex(), SimVertex::parentIndex(), pat::PATGenCandsFromSimTracksProducer::GlobalContext::simtks, pat::PATGenCandsFromSimTracksProducer::GlobalContext::simvtxs, and SimTrack::vertIndex().

Referenced by findRef(), and produce().

                                                                                                   {
    assert(tk.genpartIndex() == -1); // MUST NOT be called with a PYTHIA track
    if (!tk.noVertex()) {
        const SimVertex &vtx = g.simvtxs[tk.vertIndex()];
        if (!vtx.noParent()) {
            unsigned int idx = vtx.parentIndex();
            SimTrackContainer::const_iterator it = std::lower_bound(g.simtks.begin(), g.simtks.end(), idx, LessById());
            if ((it != g.simtks.end()) && (it->trackId() == idx)) {
                 return &*it;
            }
        }
    }
    return 0;
 }

edm::Ref< reco::GenParticleCollection > PATGenCandsFromSimTracksProducer::findRef	(	const SimTrack &	tk,
		GlobalContext &	g
	)		const

private

Find the GenParticle reference for a given GEANT or PYTHIA track.

if the track corresponds to a PYTHIA particle, return a ref to that particle
otherwise, if this simtrack has no mother simtrack, return a null ref
otherwise, if writeAncestors is true, make a GenParticle for it and return a ref to it
otherwise, if writeAncestors is false, return the ref to the GEANT mother of this track

Definition at line 165 of file PATGenCandsFromSimTracksProducer.cc.

References findGeantMother(), generatorRef_(), SimTrack::genpartIndex(), makeGenParticle_(), makeMotherLink_, pat::PATGenCandsFromSimTracksProducer::GlobalContext::output, AlCaHLTBitMon_ParallelJobs::p, pat::PATGenCandsFromSimTracksProducer::GlobalContext::refprod, pat::PATGenCandsFromSimTracksProducer::GlobalContext::simTksProcessed, CoreSimTrack::trackId(), and writeAncestors_.

Referenced by produce().

                                                                                     {
     if (tk.genpartIndex() != -1) return makeMotherLink_ ? generatorRef_(tk, g) : edm::Ref<reco::GenParticleCollection>();
     const SimTrack * simMother = findGeantMother(tk, g);
 
     edm::Ref<reco::GenParticleCollection> motherRef;
     if (simMother != 0) motherRef = findRef(*simMother,g);
 
     if (writeAncestors_) {
         // If writing ancestors, I need to serialize myself, and then to return a ref to me
         // But first check if I've already been serialized
         std::map<unsigned int,int>::const_iterator it = g.simTksProcessed.find(tk.trackId());
         if (it != g.simTksProcessed.end()) {
             // just return a ref to it
             assert(it->second > 0);
             return edm::Ref<reco::GenParticleCollection>(g.refprod, (it->second) - 1);
         } else {
             // make genParticle, save, update the map, and return ref to myself
             reco::GenParticle p = makeGenParticle_(tk, motherRef, g); 
             g.output.push_back(p);
             g.simTksProcessed[tk.trackId()] = g.output.size();
             return edm::Ref<reco::GenParticleCollection>(g.refprod, g.output.size()-1 );
         }
     } else {
         // Otherwise, I just return a ref to my mum
         return motherRef; 
     }
 }

edm::Ref< reco::GenParticleCollection > PATGenCandsFromSimTracksProducer::generatorRef_	(	const SimTrack &	tk,
		const GlobalContext &	g
	)		const

private

Used by findRef if the track is a PYTHIA particle.

Definition at line 194 of file PATGenCandsFromSimTracksProducer.cc.

References pat::PATGenCandsFromSimTracksProducer::GlobalContext::barcodesAreSorted, spr::find(), pat::PATGenCandsFromSimTracksProducer::GlobalContext::genBarcodes, SimTrack::genpartIndex(), and pat::PATGenCandsFromSimTracksProducer::GlobalContext::gens.

Referenced by findRef().

                                                                                                 {
     assert(st.genpartIndex() != -1);        
     // Note that st.genpartIndex() is the barcode, not the index within GenParticleCollection, so I have to search the particle
     std::vector<int>::const_iterator it;
     if (g.barcodesAreSorted) {
         it = std::lower_bound(g.genBarcodes->begin(), g.genBarcodes->end(), st.genpartIndex());
     } else {
         it = std::find(       g.genBarcodes->begin(), g.genBarcodes->end(), st.genpartIndex());
     }
 
     // Check that I found something
     // I need to check '*it == st.genpartIndex()' because lower_bound just finds the right spot for an item in a sorted list, not the item
     if ((it != g.genBarcodes->end()) && (*it == st.genpartIndex())) {
         return reco::GenParticleRef(g.gens, it - g.genBarcodes->begin());
     } else {
         return reco::GenParticleRef();
     }
 }

reco::GenParticle PATGenCandsFromSimTracksProducer::makeGenParticle_	(	const SimTrack &	tk,
		const edm::Ref< reco::GenParticleCollection > &	mother,
		const GlobalContext &	g
	)		const

private

Make a GenParticle for this SimTrack, with a given mother.

Definition at line 214 of file PATGenCandsFromSimTracksProducer.cc.

References reco::CompositeRefCandidateT< D >::addMother(), DeDxDiscriminatorTools::charge(), CoreSimTrack::charge(), edm::Ref< C, T, F >::isNonnull(), CoreSimTrack::momentum(), SimTrack::noVertex(), p4, setStatus_, pat::PATGenCandsFromSimTracksProducer::GlobalContext::simvtxs, CoreSimTrack::type(), and SimTrack::vertIndex().

Referenced by findRef(), and produce().

                                                                                                                                                      {
     // Make up a GenParticleCandidate from the GEANT track info.
     int charge = static_cast<int>(tk.charge());
     Particle::LorentzVector p4 = tk.momentum();
     Particle::Point vtx; // = (0,0,0) by default
     if (!tk.noVertex()) vtx = g.simvtxs[tk.vertIndex()].position();
     GenParticle gp(charge, p4, vtx, tk.type(), setStatus_, true);
     if (mother.isNonnull()) gp.addMother(mother);
     return gp;
 }

void PATGenCandsFromSimTracksProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 226 of file PATGenCandsFromSimTracksProducer.cc.

References abs, reco::CompositeRefCandidateT< D >::addMother(), gather_cfg::cout, edm::hlt::Exception, filter_, findGeantMother(), findRef(), firstEvent_, EgammaValidation_cff::genp, genParticles_, i, edm::OrphanHandleBase::id(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::key(), makeGenParticle_(), makeMotherLink_, motherPdgIds_, motherPdts_, pdgIds_, pdts_, python.multivaluedict::sort(), src_, CoreSimTrack::type(), and writeAncestors_.

                                                   {
 
   if (firstEvent_){
     if (!pdts_.empty()) {
       pdgIds_.clear();
       for (vector<PdtEntry>::iterator itp = pdts_.begin(), edp = pdts_.end(); itp != edp; ++itp) {
     itp->setup(iSetup); // decode string->pdgId and vice-versa                                                                                                
     pdgIds_.insert(std::abs(itp->pdgId()));
       }
       pdts_.clear();
     }
     if (!motherPdts_.empty()) {
       motherPdgIds_.clear();
       for (vector<PdtEntry>::iterator itp = motherPdts_.begin(), edp = motherPdts_.end(); itp != edp; ++itp) {
     itp->setup(iSetup); // decode string->pdgId and vice-versa                                                                                                
     motherPdgIds_.insert(std::abs(itp->pdgId()));
       }
       motherPdts_.clear();
     }
     firstEvent_ = false;
   }
 
   // Simulated tracks (i.e. GEANT particles).
   Handle<SimTrackContainer> simtracks;
   event.getByLabel(src_, simtracks);
 
   // Need to check that SimTrackContainer is sorted; otherwise, copy and sort :-(
   std::auto_ptr<SimTrackContainer> simtracksTmp;
   const SimTrackContainer * simtracksSorted = &* simtracks;
   if (makeMotherLink_ || writeAncestors_) {
       if (!__gnu_cxx::is_sorted(simtracks->begin(), simtracks->end(), LessById())) {
           simtracksTmp.reset(new SimTrackContainer(*simtracks));
           std::sort(simtracksTmp->begin(), simtracksTmp->end(), LessById());
           simtracksSorted = &* simtracksTmp;
       }
   }
 
   // Get the associated vertices
   Handle<SimVertexContainer> simvertices;
   event.getByLabel(src_, simvertices);
 
   // Get the GenParticles and barcodes, if needed to set mother links
   Handle<GenParticleCollection> gens;
   Handle<std::vector<int> > genBarcodes;
   bool barcodesAreSorted = true;
   if (makeMotherLink_) {
       event.getByLabel(genParticles_, gens);
       event.getByLabel(genParticles_, genBarcodes);
       if (gens->size() != genBarcodes->size()) throw cms::Exception("Corrupt data") << "Barcodes not of the same size as GenParticles!\n";
       barcodesAreSorted = __gnu_cxx::is_sorted(genBarcodes->begin(), genBarcodes->end());
   }    
 
 
   // make the output collection
   auto_ptr<GenParticleCollection> cands(new GenParticleCollection);
   edm::RefProd<GenParticleCollection> refprod = event.getRefBeforePut<GenParticleCollection>();
 
   GlobalContext globals(*simtracksSorted, *simvertices, gens, genBarcodes, barcodesAreSorted, *cands, refprod);
 
   for (SimTrackContainer::const_iterator isimtrk = simtracks->begin();
           isimtrk != simtracks->end(); ++isimtrk) {
 
       // Skip PYTHIA tracks.
       if (isimtrk->genpartIndex() != -1) continue; 
 
       // Maybe apply the PdgId filter
       if (!pdgIds_.empty()) { // if we have a filter on pdg ids
            if (pdgIds_.find(std::abs(isimtrk->type())) == pdgIds_.end()) continue;
       }
 
       GenParticle genp = makeGenParticle_(*isimtrk, Ref<GenParticleCollection>(), globals);
 
       // Maybe apply filter on the particle
       if (filter_.get() != 0) {
         if (!(*filter_)(genp)) continue;
       }
 
       if (!motherPdgIds_.empty()) {
            const SimTrack *motherSimTk = findGeantMother(*isimtrk, globals);
            if (motherSimTk == 0) continue;
            if (motherPdgIds_.find(std::abs(motherSimTk->type())) == motherPdgIds_.end()) continue;
       }
 
       if (makeMotherLink_ || writeAncestors_) {
           Ref<GenParticleCollection> motherRef;
           const SimTrack * mother = findGeantMother(*isimtrk, globals);
           if (mother != 0) motherRef = findRef(*mother, globals); 
           if (motherRef.isNonnull()) genp.addMother(motherRef);
       }
 
       cands->push_back(genp);
   }
 
   // Write to the Event, and get back a handle (which can be useful for debugging)
   edm::OrphanHandle<reco::GenParticleCollection> orphans = event.put(cands);
 
 #ifdef DEBUG_PATGenCandsFromSimTracksProducer
   std::cout << "Produced a list of " << orphans->size() << " genParticles." << std::endl;
   for (GenParticleCollection::const_iterator it = orphans->begin(), ed = orphans->end(); it != ed; ++it) {
       std::cout << "    ";
       std::cout << "GenParticle #" << (it - orphans->begin()) << ": pdgId " << it->pdgId()
                 << ", pt = " << it->pt() << ", eta = " << it->eta() << ", phi = " << it->phi() 
                 << ", rho = " << it->vertex().Rho() << ", z = " << it->vertex().Z() << std::endl;
       edm::Ref<GenParticleCollection> mom = it->motherRef();
       size_t depth = 2;
       while (mom.isNonnull()) {
           if (mom.id() == orphans.id()) {
               // I need to re-make the ref because they are not working until this module returns.
               mom = edm::Ref<GenParticleCollection>(orphans, mom.key());
           }
           for (size_t i = 0; i < depth; ++i) std::cout << "    ";
           std::cout << "GenParticleRef [" << mom.id() << "/" << mom.key() << "]: pdgId " << mom->pdgId() << ", status = " << mom->status() 
                     << ", pt = " << mom->pt() << ", eta = " << mom->eta() << ", phi = " << mom->phi() 
                     << ", rho = " << mom->vertex().Rho() << ", z = " << mom->vertex().Z() << std::endl;
           if (mom.id() != orphans.id()) break;
           if ((mom->motherRef().id() == mom.id()) && (mom->motherRef().key() == mom.key())) {
               throw cms::Exception("Corrupt Data") << "A particle is it's own mother.\n";
           }
           mom = mom->motherRef(); 
           depth++; 
       }
   }
   std::cout << std::endl;
 #endif
 
 }