#include <HLTTauMCProducer.h>

Inheritance diagram for HLTTauMCProducer:

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

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

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

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

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
enum	tauDecayModes { kElectron, kMuon, kOneProng0pi0, kOneProng1pi0, kOneProng2pi0, kThreeProng0pi0, kThreeProng1pi0, kOther, kUndefined }

Private Member Functions
void	getGenDecayProducts (const reco::GenParticleRef &, reco::GenParticleRefVector &, int status=1, int pdgId=0) const

Private Attributes
const double	etaMax_

const double	etaMin_

const std::vector< int >	m_PDG_

const edm::EDGetTokenT < reco::GenParticleCollection >	MC_

const edm::EDGetTokenT < reco::GenMETCollection >	MCMET_

const double	phiMax_

const double	phiMin_

const double	ptMinMCElectron_

const double	ptMinMCMuon_

const double	ptMinMCTau_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 33 of file HLTTauMCProducer.h.

Member Enumeration Documentation

enum HLTTauMCProducer::tauDecayModes

private

Enumerator
kElectron
kMuon
kOneProng0pi0
kOneProng1pi0
kOneProng2pi0
kThreeProng0pi0
kThreeProng1pi0
kOther
kUndefined

Definition at line 45 of file HLTTauMCProducer.h.

                      {
     kElectron,
     kMuon,
     kOneProng0pi0,
     kOneProng1pi0,
     kOneProng2pi0,
     kThreeProng0pi0,
     kThreeProng1pi0,
     kOther,
     kUndefined
   };

Constructor & Destructor Documentation

HLTTauMCProducer::HLTTauMCProducer ( const edm::ParameterSet & mc )

explicit

Definition at line 7 of file HLTTauMCProducer.cc.

References edm::ParameterSet::getUntrackedParameter().

     : MC_{consumes<GenParticleCollection>(mc.getUntrackedParameter<edm::InputTag>("GenParticles"))},
       MCMET_{consumes<GenMETCollection>(mc.getUntrackedParameter<edm::InputTag>("GenMET"))},
       ptMinMCTau_{mc.getUntrackedParameter<double>("ptMinTau", 5.)},
       ptMinMCElectron_{mc.getUntrackedParameter<double>("ptMinElectron", 5.)},
       ptMinMCMuon_{mc.getUntrackedParameter<double>("ptMinMuon", 2.)},
       m_PDG_{mc.getUntrackedParameter<std::vector<int>>("BosonID")},
       etaMin_{mc.getUntrackedParameter<double>("EtaMin", -2.5)},
       etaMax_{mc.getUntrackedParameter<double>("EtaMax", 2.5)},
       phiMin_{mc.getUntrackedParameter<double>("PhiMin", -3.15)},
       phiMax_{mc.getUntrackedParameter<double>("PhiMax", 3.15)} {
   // One Parameter Set per Collection
 
   produces<LorentzVectorCollection>("LeptonicTauLeptons");
   produces<LorentzVectorCollection>("LeptonicTauElectrons");
   produces<LorentzVectorCollection>("LeptonicTauMuons");
   produces<LorentzVectorCollection>("HadronicTauOneProng");
   produces<LorentzVectorCollection>("HadronicTauThreeProng");
   produces<LorentzVectorCollection>("HadronicTauOneAndThreeProng");
   produces<LorentzVectorCollection>("TauOther");
   produces<LorentzVectorCollection>("Neutrina");
   produces<LorentzVectorCollection>("MET");
   produces<std::vector<int>>("Mothers");
 }

Member Function Documentation

void HLTTauMCProducer::getGenDecayProducts	(	const reco::GenParticleRef &	mother,
		reco::GenParticleRefVector &	products,
		int	status = `1`,
		int	pdgId = `0`
	)		const

private

Definition at line 262 of file HLTTauMCProducer.cc.

References funct::abs(), edm::RefVector< C, T, F >::begin(), ztail::d, edm::RefVector< C, T, F >::end(), edm::RefVector< C, T, F >::push_back(), and mps_update::status.

Referenced by produce().

                                                             {
   const GenParticleRefVector &daughterRefs = mother->daughterRefVector();
 
   for (GenParticleRefVector::const_iterator d = daughterRefs.begin(); d != daughterRefs.end(); ++d) {
     if ((status == 0 || (*d)->status() == status) && (pdgId == 0 || std::abs((*d)->pdgId()) == pdgId)) {
       products.push_back(*d);
     } else
       getGenDecayProducts(*d, products, status, pdgId);
   }
 }

void HLTTauMCProducer::produce	(	edm::StreamID	,
		edm::Event &	iEvent,
		const edm::EventSetup &	iES
	)		const

overridevirtual

Implements edm::global::EDProducerBase.

Definition at line 32 of file HLTTauMCProducer.cc.

References funct::abs(), edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::empty(), edm::RefVector< C, T, F >::end(), etaMax_, etaMin_, HepMCValidationHelper::genMet(), genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByToken(), getGenDecayProducts(), reco::GenParticle::isPromptDecayed(), edm::HandleBase::isValid(), kElectron, kMuon, kOneProng0pi0, kOneProng1pi0, kOneProng2pi0, kOther, kThreeProng0pi0, kThreeProng1pi0, m_PDG_, MC_, MCMET_, HLT_FULL_cff::MET, eostools::move(), patZpeak::numMuons, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::pdgId(), phiMax_, phiMin_, pi, ptMinMCElectron_, ptMinMCMuon_, ptMinMCTau_, edm::RefVector< C, T, F >::push_back(), edm::Event::put(), and submitPVValidationJobs::t.

                                                                                             {
   // All the code from HLTTauMCInfo is here :-)
 
   unique_ptr<LorentzVectorCollection> product_Electrons(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_Muons(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_Leptons(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_OneProng(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_ThreeProng(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_OneAndThreeProng(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_Other(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_Neutrina(new LorentzVectorCollection);
   unique_ptr<LorentzVectorCollection> product_MET(new LorentzVectorCollection);
   unique_ptr<std::vector<int>> product_Mothers(new std::vector<int>);
 
   edm::Handle<GenParticleCollection> genParticles;
   iEvent.getByToken(MC_, genParticles);
 
   if (!genParticles.isValid())
     return;
 
   // Look for MET
   edm::Handle<reco::GenMETCollection> genMet;
   iEvent.getByToken(MCMET_, genMet);
   LorentzVector MET(0., 0., 0., 0.);
   if (genMet.isValid()) {
     MET = LorentzVector(genMet->front().px(), genMet->front().py(), 0, genMet->front().pt());
   }
   product_MET->push_back(MET);
 
   // Look for primary bosons
   // It is not guaranteed that primary bosons are stored in event history.
   // Is it really needed when check if taus from the boson is removed?
   // Kept for backward compatibility
   for (GenParticleCollection::const_iterator p = genParticles->begin(); p != genParticles->end(); ++p) {
     // Check the PDG ID
     bool pdg_ok = false;
     for (size_t pi = 0; pi < m_PDG_.size(); ++pi) {
       if (abs((*p).pdgId()) == m_PDG_[pi] && ((*p).isHardProcess() || (*p).status() == 3)) {
         pdg_ok = true;
         // cout<<" Bsoson particles: "<< (*p).pdgId()<< " " <<(*p).status() << "
         // "<< pdg_ok<<endl;
         break;
       }
     }
 
     // Check if the boson is one of interest and if there is a valid vertex
     if (pdg_ok) {
       product_Mothers->push_back((*p).pdgId());
 
       TLorentzVector Boson((*p).px(), (*p).py(), (*p).pz(), (*p).energy());
     }
   }  // End of search for the bosons
 
   // Look for taus
   GenParticleRefVector allTaus;
   unsigned index = 0;
   for (GenParticleCollection::const_iterator p = genParticles->begin(); p != genParticles->end(); ++p, ++index) {
     const GenParticle &genP = *p;
     // accept only isPromptDecayed() particles
     if (!genP.isPromptDecayed())
       continue;
     // check if it is tau, i.e. if |pdgId|=15
     if (std::abs(genP.pdgId()) == 15) {
       GenParticleRef genRef(genParticles, index);
       // check if it is the last tau in decay/radiation chain
       GenParticleRefVector daugTaus;
       getGenDecayProducts(genRef, daugTaus, 0, 15);
       if (daugTaus.empty())
         allTaus.push_back(genRef);
     }
   }
 
   // Find stable tau decay products and build visible taus
   for (GenParticleRefVector::const_iterator t = allTaus.begin(); t != allTaus.end(); ++t) {
     // look for all stable (status=1) decay products
     GenParticleRefVector decayProducts;
     getGenDecayProducts(*t, decayProducts, 1);
 
     // build visible taus and recognize decay mode
     if (!decayProducts.empty()) {
       LorentzVector Visible_Taus(0., 0., 0., 0.);
       LorentzVector TauDecayProduct(0., 0., 0., 0.);
       LorentzVector Neutrino(0., 0., 0., 0.);
 
       int numElectrons = 0;
       int numMuons = 0;
       int numChargedPions = 0;
       int numNeutralPions = 0;
       int numPhotons = 0;
       int numNeutrinos = 0;
       int numOtherParticles = 0;
 
       for (GenParticleRefVector::const_iterator pit = decayProducts.begin(); pit != decayProducts.end(); ++pit) {
         int pdg_id = abs((*pit)->pdgId());
         if (pdg_id == 11)
           numElectrons++;
         else if (pdg_id == 13)
           numMuons++;
         else if (pdg_id == 211 || pdg_id == 321)
           numChargedPions++;  // Count both pi+ and K+
         else if (pdg_id == 111 || pdg_id == 130 || pdg_id == 310)
           numNeutralPions++;  // Count both pi0 and K0_L/S
         else if (pdg_id == 12 || pdg_id == 14 || pdg_id == 16) {
           numNeutrinos++;
           if (pdg_id == 16) {
             Neutrino.SetPxPyPzE((*pit)->px(), (*pit)->py(), (*pit)->pz(), (*pit)->energy());
           }
         } else if (pdg_id == 22)
           numPhotons++;
         else {
           numOtherParticles++;
         }
 
         if (pdg_id != 12 && pdg_id != 14 && pdg_id != 16) {
           TauDecayProduct.SetPxPyPzE((*pit)->px(), (*pit)->py(), (*pit)->pz(), (*pit)->energy());
           Visible_Taus += TauDecayProduct;
         }
         //        cout<< "This has to be the same: " << (*pit)->pdgId()
         //<< " "<< (*pit)->status()<< " mother: "<< (*pit)->mother()->pdgId() <<
         // endl;
       }
 
       int tauDecayMode = kOther;
 
       if (numOtherParticles == 0) {
         if (numElectrons == 1) {
           //--- tau decays into electrons
           tauDecayMode = kElectron;
         } else if (numMuons == 1) {
           //--- tau decays into muons
           tauDecayMode = kMuon;
         } else {
           //--- hadronic tau decays
           switch (numChargedPions) {
             case 1:
               if (numNeutralPions != 0) {
                 tauDecayMode = kOther;
                 break;
               }
               switch (numPhotons) {
                 case 0:
                   tauDecayMode = kOneProng0pi0;
                   break;
                 case 2:
                   tauDecayMode = kOneProng1pi0;
                   break;
                 case 4:
                   tauDecayMode = kOneProng2pi0;
                   break;
                 default:
                   tauDecayMode = kOther;
                   break;
               }
               break;
             case 3:
               if (numNeutralPions != 0) {
                 tauDecayMode = kOther;
                 break;
               }
               switch (numPhotons) {
                 case 0:
                   tauDecayMode = kThreeProng0pi0;
                   break;
                 case 2:
                   tauDecayMode = kThreeProng1pi0;
                   break;
                 default:
                   tauDecayMode = kOther;
                   break;
               }
               break;
           }
         }
       }
 
       //          cout<< "So we have a: " << tauDecayMode <<endl;
       if (tauDecayMode == kElectron) {
         if ((Visible_Taus.eta() > etaMin_ && Visible_Taus.eta() < etaMax_ && Visible_Taus.phi() > phiMin_ &&
              Visible_Taus.phi() < phiMax_) &&
             (Visible_Taus.pt() > ptMinMCElectron_)) {
           product_Electrons->push_back(Visible_Taus);
           product_Leptons->push_back(Visible_Taus);
         }
       } else if (tauDecayMode == kMuon) {
         if ((Visible_Taus.eta() > etaMin_ && Visible_Taus.eta() < etaMax_ && Visible_Taus.phi() > phiMin_ &&
              Visible_Taus.phi() < phiMax_) &&
             (Visible_Taus.pt() > ptMinMCMuon_)) {
           product_Muons->push_back(Visible_Taus);
           product_Leptons->push_back(Visible_Taus);
         }
       } else if (tauDecayMode == kOneProng0pi0 || tauDecayMode == kOneProng1pi0 || tauDecayMode == kOneProng2pi0) {
         if ((Visible_Taus.eta() > etaMin_ && Visible_Taus.eta() < etaMax_ && Visible_Taus.phi() > phiMin_ &&
              Visible_Taus.phi() < phiMax_) &&
             (Visible_Taus.pt() > ptMinMCTau_)) {
           product_OneProng->push_back(Visible_Taus);
           product_OneAndThreeProng->push_back(Visible_Taus);
           product_Neutrina->push_back(Neutrino);
         }
       } else if (tauDecayMode == kThreeProng0pi0 || tauDecayMode == kThreeProng1pi0) {
         if ((Visible_Taus.eta() > etaMin_ && Visible_Taus.eta() < etaMax_ && Visible_Taus.phi() > phiMin_ &&
              Visible_Taus.phi() < phiMax_) &&
             (Visible_Taus.pt() > ptMinMCTau_)) {
           product_ThreeProng->push_back(Visible_Taus);
           product_OneAndThreeProng->push_back(Visible_Taus);
           product_Neutrina->push_back(Neutrino);
         }
       } else if (tauDecayMode == kOther) {
         if ((Visible_Taus.eta() > etaMin_ && Visible_Taus.eta() < etaMax_ && Visible_Taus.phi() > phiMin_ &&
              Visible_Taus.phi() < phiMax_) &&
             (Visible_Taus.pt() > ptMinMCTau_)) {
           product_Other->push_back(Visible_Taus);
         }
       }
     }
   }
 
   iEvent.put(std::move(product_Leptons), "LeptonicTauLeptons");
   iEvent.put(std::move(product_Electrons), "LeptonicTauElectrons");
   iEvent.put(std::move(product_Muons), "LeptonicTauMuons");
   iEvent.put(std::move(product_OneProng), "HadronicTauOneProng");
   iEvent.put(std::move(product_ThreeProng), "HadronicTauThreeProng");
   iEvent.put(std::move(product_OneAndThreeProng), "HadronicTauOneAndThreeProng");
   iEvent.put(std::move(product_Other), "TauOther");
   iEvent.put(std::move(product_Neutrina), "Neutrina");
   iEvent.put(std::move(product_MET), "MET");
   iEvent.put(std::move(product_Mothers), "Mothers");
 }

Member Data Documentation

const double HLTTauMCProducer::etaMax_

private

Definition at line 63 of file HLTTauMCProducer.h.

Referenced by produce().

const double HLTTauMCProducer::etaMin_

private

Definition at line 63 of file HLTTauMCProducer.h.

Referenced by produce().

const std::vector<int> HLTTauMCProducer::m_PDG_

private

Definition at line 62 of file HLTTauMCProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::GenParticleCollection> HLTTauMCProducer::MC_

private

Definition at line 57 of file HLTTauMCProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::GenMETCollection> HLTTauMCProducer::MCMET_

private

Definition at line 58 of file HLTTauMCProducer.h.

Referenced by produce().

const double HLTTauMCProducer::phiMax_

private

Definition at line 63 of file HLTTauMCProducer.h.

Referenced by produce().

const double HLTTauMCProducer::phiMin_

private

Definition at line 63 of file HLTTauMCProducer.h.

Referenced by produce().

const double HLTTauMCProducer::ptMinMCElectron_

private

Definition at line 60 of file HLTTauMCProducer.h.

Referenced by produce().

const double HLTTauMCProducer::ptMinMCMuon_

private

Definition at line 61 of file HLTTauMCProducer.h.

Referenced by produce().

const double HLTTauMCProducer::ptMinMCTau_

private

Definition at line 59 of file HLTTauMCProducer.h.

Referenced by produce().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation