#include <VBFGenJetFilter.h>

Inheritance diagram for VBFGenJetFilter:

Public Member Functions
bool	filter (edm::Event &, const edm::EventSetup &) override

	VBFGenJetFilter (const edm::ParameterSet &)

	~VBFGenJetFilter () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
int	charge (const int &Id)

std::vector< const reco::GenJet * >	filterGenJets (const std::vector< reco::GenJet > *jets)

std::vector< const reco::GenParticle * >	filterGenLeptons (const std::vector< reco::GenParticle > *particles)

std::vector< HepMC::GenParticle * >	getNu (const HepMC::GenEvent *particles)

std::vector< HepMC::GenParticle * >	getSt3 (const HepMC::GenEvent *particles)

std::vector< HepMC::GenParticle * >	getVisibleDecayProducts (HepMC::GenParticle *particle)

double	nuMET (std::vector< HepMC::GenParticle * > vNu)

void	printGenVector (std::vector< HepMC::GenParticle * > vec)

Private Attributes
double	deltaRJetLep

double	etaMax

double	etaMin

bool	leadJetsNoLepMass

edm::EDGetTokenT< reco::GenJetCollection >	m_inputTag_GenJetCollection

edm::EDGetTokenT< reco::GenParticleCollection >	m_inputTag_GenParticleCollection

double	maxDeltaEta

double	maxDeltaPhi

double	maxInvMass

double	maxLeadingJetsInvMass

double	minDeltaEta

double	minDeltaPhi

double	minInvMass

double	minLeadingJetsInvMass

bool	oppositeHemisphere

double	ptMin

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	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::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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

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

Definition at line 29 of file VBFGenJetFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file VBFGenJetFilter.cc.

References edm::ParameterSet::getUntrackedParameter(), leadJetsNoLepMass, m_inputTag_GenJetCollection, and m_inputTag_GenParticleCollection.

                                                                :
 oppositeHemisphere(iConfig.getUntrackedParameter<bool>  ("oppositeHemisphere",false)),
 leadJetsNoLepMass (iConfig.getUntrackedParameter<bool>  ("leadJetsNoLepMass", false)),
 ptMin             (iConfig.getUntrackedParameter<double>("minPt",                20)),
 etaMin            (iConfig.getUntrackedParameter<double>("minEta",             -5.0)),
 etaMax            (iConfig.getUntrackedParameter<double>("maxEta",              5.0)),
 minInvMass        (iConfig.getUntrackedParameter<double>("minInvMass",          0.0)),
 maxInvMass        (iConfig.getUntrackedParameter<double>("maxInvMass",      99999.0)),
 minDeltaPhi       (iConfig.getUntrackedParameter<double>("minDeltaPhi",        -1.0)),
 maxDeltaPhi       (iConfig.getUntrackedParameter<double>("maxDeltaPhi",     99999.0)),
 minDeltaEta       (iConfig.getUntrackedParameter<double>("minDeltaEta",        -1.0)),
 maxDeltaEta       (iConfig.getUntrackedParameter<double>("maxDeltaEta",     99999.0)),
 minLeadingJetsInvMass (iConfig.getUntrackedParameter<double>("minLeadingJetsInvMass",          0.0)),
 maxLeadingJetsInvMass (iConfig.getUntrackedParameter<double>("maxLeadingJetsInvMass",      99999.0)),
 deltaRJetLep          (iConfig.getUntrackedParameter<double>("deltaRJetLep",                   0.3))
 {
   
   m_inputTag_GenJetCollection = consumes<reco::GenJetCollection>(iConfig.getUntrackedParameter<edm::InputTag>("inputTag_GenJetCollection",edm::InputTag("ak5GenJetsNoNu")));
   if (leadJetsNoLepMass) m_inputTag_GenParticleCollection  = consumes<reco::GenParticleCollection>(iConfig.getUntrackedParameter<edm::InputTag>("genParticles",edm::InputTag("genParticles")));
   
 }

VBFGenJetFilter::~VBFGenJetFilter ( )

override

Definition at line 40 of file VBFGenJetFilter.cc.

                                  {
   
 }

Member Function Documentation

int VBFGenJetFilter::charge ( const int & Id )

private

bool VBFGenJetFilter::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 84 of file VBFGenJetFilter.cc.

References a, b, DEFINE_FWK_MODULE, reco::deltaPhi(), reco::deltaR2(), deltaRJetLep, particleFlow_cfi::dEta, particleFlow_cfi::dPhi, filterGenJets(), filterGenLeptons(), ttbarCategorization_cff::genJets, GenHFHadronMatcher_cfi::genParticles, edm::Event::getByToken(), mps_fire::i, boostedTaus_cff::jetIdx, leadJetsNoLepMass, m_inputTag_GenJetCollection, m_inputTag_GenParticleCollection, maxDeltaEta, maxDeltaPhi, maxInvMass, maxLeadingJetsInvMass, minLeadingJetsInvMass, oppositeHemisphere, p4, reco::LeafCandidate::p4(), and edm::Handle< T >::product().

 {
   using namespace edm;
   
   Handle< vector<reco::GenJet> > handleGenJets;
   iEvent.getByToken(m_inputTag_GenJetCollection, handleGenJets);
   const vector<reco::GenJet>* genJets = handleGenJets.product();
   
   // Getting filtered generator jets
   vector<const reco::GenJet*> filGenJets = filterGenJets(genJets);
   
   // If we do not find at least 2 jets veto the event
   if(filGenJets.size()<2){return false;}
   
   
   // Testing dijet mass   
   if(leadJetsNoLepMass) { 
 
       Handle<reco::GenParticleCollection> genParticelesCollection;
       iEvent.getByToken(m_inputTag_GenParticleCollection, genParticelesCollection);
       const vector<reco::GenParticle>*  genParticles = genParticelesCollection.product();
 
       // Getting filtered generator muons
       vector<const reco::GenParticle*> filGenLep = filterGenLeptons(genParticles);
 
       // Getting p4 of jet with no lepton
       vector<math::XYZTLorentzVector> genJetsWithoutLeptonsP4;
       unsigned int jetIdx = 0;
 
       while(genJetsWithoutLeptonsP4.size()<2 && jetIdx < filGenJets.size()) {
           bool jetWhitoutLep = true;
     
           const math::XYZTLorentzVector & p4J= (filGenJets[jetIdx])->p4();
           for(unsigned int i = 0; i < filGenLep.size() && jetWhitoutLep; ++i) {
               if(reco::deltaR2((filGenLep[i])->p4(), p4J) < deltaRJetLep*deltaRJetLep)
                   jetWhitoutLep = false;
           }
           if (jetWhitoutLep)  genJetsWithoutLeptonsP4.push_back(p4J);
           ++jetIdx;
       }
 
       // Checking the invariant mass of the leading jets
       if (genJetsWithoutLeptonsP4.size() < 2) return false;
       float invMassLeadingJet = (genJetsWithoutLeptonsP4[0] + genJetsWithoutLeptonsP4[1]).M();
       if ( invMassLeadingJet > minLeadingJetsInvMass  && invMassLeadingJet < maxLeadingJetsInvMass) return true;
       else return false;
   }
   
   
   for(unsigned a=0; a<filGenJets.size(); a++){
     for(unsigned b=a+1; b<filGenJets.size(); b++){    
       
       const reco::GenJet* pA = filGenJets[a];
       const reco::GenJet* pB = filGenJets[b];
       
       // Getting the dijet vector
       math::XYZTLorentzVector diJet = pA->p4() + pB->p4();
       
       // Testing opposite hemispheres
       double dijetProd = pA->p4().eta()*pB->p4().eta();
       if(oppositeHemisphere && dijetProd>=0){continue;}
       
       // Testing dijet mass
       double invMass = diJet.mass();
       if(invMass<=minInvMass || invMass>maxInvMass){continue;}
       
       // Testing dijet delta eta
       double dEta = fabs(pA->p4().eta()-pB->p4().eta());
       if(dEta<=minDeltaEta || dEta>maxDeltaEta){continue;}
 
       // Testing dijet delta phi
       double dPhi = fabs(reco::deltaPhi(pA->p4().phi(),pB->p4().phi()));
       if(dPhi<=minDeltaPhi || dPhi>maxDeltaPhi){continue;}
       
       return true;
     }
   }
 
   return false;
 }

vector< const reco::GenJet * > VBFGenJetFilter::filterGenJets ( const std::vector< reco::GenJet > * jets )

private

Definition at line 65 of file VBFGenJetFilter.cc.

References etaMax, etaMin, mps_fire::i, MillePedeFileConverter_cfg::out, reco::LeafCandidate::p4(), and ptMin.

Referenced by filter().

                                                                                       {
   
   vector<const reco::GenJet*> out;
   
   for(unsigned i=0; i<jets->size(); i++){
     
     const reco::GenJet* j = &((*jets)[i]);
     
     if(j->p4().pt() >ptMin &&  j->p4().eta()>etaMin && j->p4().eta()<etaMax)
     {
       out.push_back(j);
     }
   }
   
   return out;
 }

vector< const reco::GenParticle * > VBFGenJetFilter::filterGenLeptons ( const std::vector< reco::GenParticle > * particles )

private

Definition at line 45 of file VBFGenJetFilter.cc.

References funct::abs(), MillePedeFileConverter_cfg::out, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by filter().

                                                                                                         {
   vector<const reco::GenParticle*> out;
 
   for(const auto & p : *particles){
       
       int absPdgId = std::abs(p.pdgId());
       
       if(((absPdgId == 11) || (absPdgId == 13) || (absPdgId == 15)) && p.isHardProcess()) {
           out.push_back(&p);              
       }
       
   }      
   return out;
 }

std::vector<HepMC::GenParticle*> VBFGenJetFilter::getNu ( const HepMC::GenEvent * particles )

private

std::vector<HepMC::GenParticle*> VBFGenJetFilter::getSt3 ( const HepMC::GenEvent * particles )

private

std::vector<HepMC::GenParticle*> VBFGenJetFilter::getVisibleDecayProducts ( HepMC::GenParticle * particle )

private

double VBFGenJetFilter::nuMET ( std::vector< HepMC::GenParticle * > vNu )

private

void VBFGenJetFilter::printGenVector ( std::vector< HepMC::GenParticle * > vec )

private

Member Data Documentation

double VBFGenJetFilter::deltaRJetLep

private

Definition at line 68 of file VBFGenJetFilter.h.

Referenced by filter().

double VBFGenJetFilter::etaMax

private

Definition at line 59 of file VBFGenJetFilter.h.

Referenced by filterGenJets().

double VBFGenJetFilter::etaMin

private

Definition at line 58 of file VBFGenJetFilter.h.

Referenced by filterGenJets().

bool VBFGenJetFilter::leadJetsNoLepMass

private

Definition at line 56 of file VBFGenJetFilter.h.

Referenced by filter(), and VBFGenJetFilter().

edm::EDGetTokenT< reco::GenJetCollection > VBFGenJetFilter::m_inputTag_GenJetCollection

private

Definition at line 71 of file VBFGenJetFilter.h.

Referenced by filter(), and VBFGenJetFilter().

edm::EDGetTokenT< reco::GenParticleCollection > VBFGenJetFilter::m_inputTag_GenParticleCollection

private

Definition at line 72 of file VBFGenJetFilter.h.

Referenced by filter(), and VBFGenJetFilter().

double VBFGenJetFilter::maxDeltaEta

private

Definition at line 65 of file VBFGenJetFilter.h.

Referenced by filter().

double VBFGenJetFilter::maxDeltaPhi

private

Definition at line 63 of file VBFGenJetFilter.h.

Referenced by filter().

double VBFGenJetFilter::maxInvMass

private

Definition at line 61 of file VBFGenJetFilter.h.

Referenced by filter().

double VBFGenJetFilter::maxLeadingJetsInvMass

private

Definition at line 67 of file VBFGenJetFilter.h.

Referenced by filter().

double VBFGenJetFilter::minDeltaEta

private

Definition at line 64 of file VBFGenJetFilter.h.

double VBFGenJetFilter::minDeltaPhi

private

Definition at line 62 of file VBFGenJetFilter.h.

double VBFGenJetFilter::minInvMass

private

Definition at line 60 of file VBFGenJetFilter.h.

double VBFGenJetFilter::minLeadingJetsInvMass

private

Definition at line 66 of file VBFGenJetFilter.h.

Referenced by filter().

bool VBFGenJetFilter::oppositeHemisphere

private

Definition at line 55 of file VBFGenJetFilter.h.

Referenced by filter().

double VBFGenJetFilter::ptMin

private

Definition at line 57 of file VBFGenJetFilter.h.

Referenced by filterGenJets().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation