#include <ttbarEventSelector.h>

Inheritance diagram for ttbarEventSelector:

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

	ttbarEventSelector (const edm::ParameterSet &)

Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default

	EDFilter (const EDFilter &)=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

const EDFilter &	operator= (const EDFilter &)=delete

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

	EDFilterBase (const EDFilterBase &)=delete

ModuleDescription const &	moduleDescription () const

const EDFilterBase &	operator= (const EDFilterBase &)=delete

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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

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

std::vector< ESResolverIndex > 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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
int	EventCategory (int &nEle, int &nMu, int &nJets, int &nbJets)

double	getMt (const TLorentzVector &vlep, const reco::PFMET &obj)

Private Attributes
const edm::InputTag	bjetsTag_

const edm::EDGetTokenT< reco::JetTagCollection >	bjetsToken_

const edm::InputTag	bsTag_

const edm::EDGetTokenT< reco::BeamSpot >	bsToken_

const double	btagFactor_

const edm::InputTag	electronTag_

const edm::EDGetTokenT< reco::GsfElectronCollection >	electronToken_

const edm::InputTag	jetsTag_

const edm::EDGetTokenT< reco::PFJetCollection >	jetsToken_

const double	maxD0_

const double	maxDeltaEtaInEB_

const double	maxDeltaEtaInEE_

const double	maxDeltaPhiInEB_

const double	maxDeltaPhiInEE_

const double	maxDz_

const double	maxEta_Jets_

const double	maxEtaEle_

const double	maxEtaHighest_Jets_

const double	maxEtaMu_

const double	maxHOEEB_

const double	maxHOEEE_

const double	maxInvMass_

const double	maxIsoEle_

const double	maxIsoMu_

const double	maxMet_

const double	maxNormChi2_

const double	maxSigmaiEiEEB_

const double	maxSigmaiEiEEE_

const double	maxWmass_

const double	minChambers_

const double	minInvMass_

const double	minMatchedStations_

const double	minMatches_

const double	minMet_

const int	minPixelHits_

const double	minPt_

const double	minPt_Jets_

const double	minPtHighest_Jets_

const double	minPtHighestEle_

const double	minPtHighestMu_

const int	minStripHits_

const double	minWmass_

const edm::InputTag	muonTag_

const edm::EDGetTokenT< reco::MuonCollection >	muonToken_

const edm::InputTag	pfmetTag_

const edm::EDGetTokenT< reco::PFMETCollection >	pfmetToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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::stream::EDFilterBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

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 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 38 of file ttbarEventSelector.h.

Constructor & Destructor Documentation

◆ ttbarEventSelector()

ttbarEventSelector::ttbarEventSelector ( const edm::ParameterSet & ps )

explicit

Definition at line 15 of file ttbarEventSelector.cc.

     : electronTag_(ps.getUntrackedParameter<edm::InputTag>("electronInputTag", edm::InputTag("gedGsfElectrons"))),
       jetsTag_(ps.getUntrackedParameter<edm::InputTag>("jetsInputTag", edm::InputTag("ak4PFJetsCHS"))),
       bjetsTag_(ps.getUntrackedParameter<edm::InputTag>("bjetsInputTag", edm::InputTag("pfDeepCSVJetTags", "probb"))),
       pfmetTag_(ps.getUntrackedParameter<edm::InputTag>("pfmetTag", edm::InputTag("pfMet"))),
       muonTag_(ps.getUntrackedParameter<edm::InputTag>("muonInputTag", edm::InputTag("muons"))),
       bsTag_(ps.getUntrackedParameter<edm::InputTag>("offlineBeamSpot", edm::InputTag("offlineBeamSpot"))),
       electronToken_(consumes<reco::GsfElectronCollection>(electronTag_)),
       jetsToken_(consumes<reco::PFJetCollection>(jetsTag_)),
       bjetsToken_(consumes<reco::JetTagCollection>(bjetsTag_)),
       pfmetToken_(consumes<reco::PFMETCollection>(pfmetTag_)),
       muonToken_(consumes<reco::MuonCollection>(muonTag_)),
       bsToken_(consumes<reco::BeamSpot>(bsTag_)),
       maxEtaEle_(ps.getUntrackedParameter<double>("maxEta", 2.4)),
       maxEtaMu_(ps.getUntrackedParameter<double>("maxEta", 2.4)),
       minPt_(ps.getUntrackedParameter<double>("minPt", 5)),
 
       // for Electron only
       maxDeltaPhiInEB_(ps.getUntrackedParameter<double>("maxDeltaPhiInEB", .15)),
       maxDeltaEtaInEB_(ps.getUntrackedParameter<double>("maxDeltaEtaInEB", .007)),
       maxHOEEB_(ps.getUntrackedParameter<double>("maxHOEEB", .12)),
       maxSigmaiEiEEB_(ps.getUntrackedParameter<double>("maxSigmaiEiEEB", .01)),
       maxDeltaPhiInEE_(ps.getUntrackedParameter<double>("maxDeltaPhiInEE", .1)),
       maxDeltaEtaInEE_(ps.getUntrackedParameter<double>("maxDeltaEtaInEE", .009)),
       maxHOEEE_(ps.getUntrackedParameter<double>("maxHOEEB_", .10)),
       maxSigmaiEiEEE_(ps.getUntrackedParameter<double>("maxSigmaiEiEEE", .03)),
 
       // for Muon only
       minChambers_(ps.getUntrackedParameter<uint32_t>("minChambers", 2)),
       minMatches_(ps.getUntrackedParameter<uint32_t>("minMatches", 2)),
       minMatchedStations_(ps.getUntrackedParameter<double>("minMatchedStations", 2)),
 
       // for Jets only
       maxEtaHighest_Jets_(ps.getUntrackedParameter<double>("maxEtaHighest_Jets", 2.4)),
       maxEta_Jets_(ps.getUntrackedParameter<double>("maxEta_Jets", 3.0)),
 
       // for b-tag only
       btagFactor_(ps.getUntrackedParameter<double>("btagFactor", 0.6)),
 
       maxNormChi2_(ps.getUntrackedParameter<double>("maxNormChi2", 10)),
       maxD0_(ps.getUntrackedParameter<double>("maxD0", 0.02)),
       maxDz_(ps.getUntrackedParameter<double>("maxDz", 20.)),
       minPixelHits_(ps.getUntrackedParameter<uint32_t>("minPixelHits", 1)),
       minStripHits_(ps.getUntrackedParameter<uint32_t>("minStripHits", 8)),
       maxIsoEle_(ps.getUntrackedParameter<double>("maxIso", 0.5)),
       maxIsoMu_(ps.getUntrackedParameter<double>("maxIso", 0.3)),
       minPtHighestMu_(ps.getUntrackedParameter<double>("minPtHighestMu", 24)),
       minPtHighestEle_(ps.getUntrackedParameter<double>("minPtHighestEle", 32)),
       minPtHighest_Jets_(ps.getUntrackedParameter<double>("minPtHighest_Jets", 30)),
       minPt_Jets_(ps.getUntrackedParameter<double>("minPt_Jets", 20)),
       minInvMass_(ps.getUntrackedParameter<double>("minInvMass", 140)),
       maxInvMass_(ps.getUntrackedParameter<double>("maxInvMass", 200)),
       minMet_(ps.getUntrackedParameter<double>("minMet", 50)),
       maxMet_(ps.getUntrackedParameter<double>("maxMet", 80)),
       minWmass_(ps.getUntrackedParameter<double>("minWmass", 50)),
       maxWmass_(ps.getUntrackedParameter<double>("maxWmass", 130)) {}

Member Function Documentation

◆ EventCategory()

int ttbarEventSelector::EventCategory	(	int &	nEle,
		int &	nMu,
		int &	nJets,
		int &	nbJets
	)

private

Definition at line 362 of file ttbarEventSelector.cc.

References eostools::cat(), l1tJetFileWriter_cfi::nJets, and l1ctLayer1_cff::nMu.

Referenced by filter().

                                                                                   {
   int cat = 0;
   if ((nEle >= 2 || nMu >= 2) && nJets > 1 && nbJets > 1) {  // di-lepton
     if (nEle >= 2 && nJets < 2)
       cat = 11;
     else if (nMu >= 2 && nJets < 2)
       cat = 12;
     else if (nEle >= 1 && nMu >= 1 && nJets < 2)
       cat = 13;
   } else if ((nEle >= 1 || nMu >= 1) && nJets > 3 && nbJets > 2) {  //semi-lepton
     if (nEle >= 1 && nJets > 1)
       cat = 21;
     else if (nMu >= 1 && nJets > 1)
       cat = 22;
   } else if ((nEle < 1 && nMu < 1) && nJets > 5 && nbJets > 1)
     cat = 30;
   else
     cat = 50;
   return cat;
 }

◆ filter()

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

overridevirtual

Implements edm::stream::EDFilterBase.

Definition at line 72 of file ttbarEventSelector.cc.

References pwdgSkimBPark_cfi::beamSpot, bjetsToken_, bsTag_, bsToken_, btagFactor_, TtFullHadKinFitProducer_cfi::bTags, nano_mu_local_reco_cff::chi2, electronIdCutBasedExt_cfi::deltaEtaIn, electronIdCutBasedExt_cfi::deltaPhiIn, electronTag_, electronToken_, EventCategory(), dqmdumpme::first, getMt(), muon::GlobalMuonPromptTight, EgHLTOffHistBins_cfi::hOverE, mps_fire::i, iEvent, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), dqmiolumiharvest::j, PDWG_EXODelayedJetMET_cff::jets, jetsToken_, SiStripPI::max, maxD0_, maxDeltaEtaInEB_, maxDeltaEtaInEE_, maxDeltaPhiInEB_, maxDeltaPhiInEE_, maxDz_, maxEta_Jets_, maxEtaMu_, maxHOEEB_, maxHOEEE_, maxIsoEle_, maxIsoMu_, maxMet_, maxNormChi2_, maxSigmaiEiEEB_, maxSigmaiEiEEE_, maxWmass_, minChambers_, minMatchedStations_, minMatches_, minPixelHits_, minPt_, minPt_Jets_, minPtHighest_Jets_, minPtHighestEle_, minPtHighestMu_, minStripHits_, minWmass_, TtSemiLepEvtBuilder_cfi::mt, MT2Analyzer::mt2, amptDefaultParameters_cff::mu, muonTag_, muonToken_, ndof, reco::HitPattern::numberOfValidPixelHits(), reco::HitPattern::numberOfValidStripHits(), nano_mu_reco_cff::pfIso04, pfmetTag_, pfmetToken_, edm::Handle< T >::product(), edm::second(), reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, and reco::GsfElectron::PflowIsolationVariables::sumPUPt.

                                                                              {
   // beamspot
   edm::Handle<reco::BeamSpot> beamSpot;
   iEvent.getByToken(bsToken_, beamSpot);
 
   int le = 0, lm = 0, lj = 0, lbj = 0;
 
   // Read Electron Collection
   edm::Handle<reco::GsfElectronCollection> electronColl;
   iEvent.getByToken(electronToken_, electronColl);
   std::vector<TLorentzVector> list_ele;
   std::vector<int> chrgeList_ele;
 
   if (electronColl.isValid()) {
     for (auto const& ele : *electronColl) {
       if (!ele.ecalDriven())
         continue;
       if (ele.pt() < minPt_)
         continue;
       // set a max Eta cut
       if (!(ele.isEB() || ele.isEE()))
         continue;
 
       double hOverE = ele.hadronicOverEm();
       double sigmaee = ele.sigmaIetaIeta();
       double deltaPhiIn = ele.deltaPhiSuperClusterTrackAtVtx();
       double deltaEtaIn = ele.deltaEtaSuperClusterTrackAtVtx();
 
       // separate cut for barrel and endcap
       if (ele.isEB()) {
         if (fabs(deltaPhiIn) >= maxDeltaPhiInEB_ && fabs(deltaEtaIn) >= maxDeltaEtaInEB_ && hOverE >= maxHOEEB_ &&
             sigmaee >= maxSigmaiEiEEB_)
           continue;
       } else if (ele.isEE()) {
         if (fabs(deltaPhiIn) >= maxDeltaPhiInEE_ && fabs(deltaEtaIn) >= maxDeltaEtaInEE_ && hOverE >= maxHOEEE_ &&
             sigmaee >= maxSigmaiEiEEE_)
           continue;
       }
 
       reco::GsfTrackRef trk = ele.gsfTrack();
       if (!trk.isNonnull())
         continue;  // only electrons with tracks
       double chi2 = trk->chi2();
       double ndof = trk->ndof();
       double chbyndof = (ndof > 0) ? chi2 / ndof : 0;
       if (chbyndof >= maxNormChi2_)
         continue;
 
       double trkd0 = trk->d0();
       if (beamSpot.isValid()) {
         trkd0 = -(trk->dxy(beamSpot->position()));
       } else {
         edm::LogError("ElectronSelector") << "Error >> Failed to get BeamSpot for label: " << bsTag_;
       }
       if (std::fabs(trkd0) >= maxD0_)
         continue;
 
       const reco::HitPattern& hitp = trk->hitPattern();
       int nPixelHits = hitp.numberOfValidPixelHits();
       if (nPixelHits < minPixelHits_)
         continue;
 
       int nStripHits = hitp.numberOfValidStripHits();
       if (nStripHits < minStripHits_)
         continue;
 
       // PF Isolation
       reco::GsfElectron::PflowIsolationVariables pfIso = ele.pfIsolationVariables();
       const float eiso =
           pfIso.sumChargedHadronPt + std::max(0.0, pfIso.sumNeutralHadronEt + pfIso.sumPhotonEt - 0.5 * pfIso.sumPUPt);
       if (eiso > maxIsoEle_ * ele.pt())
         continue;
 
       TLorentzVector lv_ele;
       lv_ele.SetPtEtaPhiE(ele.pt(), ele.eta(), ele.phi(), ele.energy());
       list_ele.push_back(lv_ele);
       chrgeList_ele.push_back(ele.charge());
     }
     le = list_ele.size();
   } else {
     edm::LogError("ElectronSelector") << "Error >> Failed to get ElectronCollection for label: " << electronTag_;
   }
 
   // Read Muon Collection
   edm::Handle<reco::MuonCollection> muonColl;
   iEvent.getByToken(muonToken_, muonColl);
 
   std::vector<TLorentzVector> list_mu;
   std::vector<int> chrgeList_mu;
   if (muonColl.isValid()) {
     for (auto const& mu : *muonColl) {
       if (!mu.isGlobalMuon())
         continue;
       if (!mu.isPFMuon())
         continue;
       if (std::fabs(mu.eta()) >= maxEtaMu_)
         continue;
       if (mu.pt() < minPt_)
         continue;
 
       reco::TrackRef gtk = mu.globalTrack();
       double chi2 = gtk->chi2();
       double ndof = gtk->ndof();
       double chbyndof = (ndof > 0) ? chi2 / ndof : 0;
       if (chbyndof >= maxNormChi2_)
         continue;
 
       reco::TrackRef tk = mu.innerTrack();
       if (beamSpot.isValid()) {
         double trkd0 = -(tk->dxy(beamSpot->position()));
         if (std::fabs(trkd0) >= maxD0_)
           continue;
         double trkdz = tk->dz(beamSpot->position());
         if (std::fabs(trkdz) >= maxDz_)
           continue;
       } else {
         edm::LogError("MuonSelector") << "Error >> Failed to get BeamSpot for label: " << bsTag_;
       }
 
       const reco::HitPattern& hitp = gtk->hitPattern();
       if (hitp.numberOfValidPixelHits() < minPixelHits_)
         continue;
       if (hitp.numberOfValidStripHits() < minStripHits_)
         continue;
 
       // Hits/section in the muon chamber
       if (mu.numberOfChambers() < minChambers_)
         continue;
       if (mu.numberOfMatches() < minMatches_)
         continue;
       if (mu.numberOfMatchedStations() < minMatchedStations_)
         continue;
       if (!muon::isGoodMuon(mu, muon::GlobalMuonPromptTight))
         continue;
 
       // PF Isolation
       const reco::MuonPFIsolation& pfIso04 = mu.pfIsolationR04();
       double absiso = pfIso04.sumChargedHadronPt +
                       std::max(0.0, pfIso04.sumNeutralHadronEt + pfIso04.sumPhotonEt - 0.5 * pfIso04.sumPUPt);
       if (absiso > maxIsoMu_ * mu.pt())
         continue;
 
       TLorentzVector lv_mu;
       lv_mu.SetPtEtaPhiE(mu.pt(), mu.eta(), mu.phi(), mu.energy());
       list_mu.push_back(lv_mu);
       chrgeList_mu.push_back(mu.charge());
     }
     lm = list_mu.size();
   } else {
     edm::LogError("MuonSelector") << "Error >> Failed to get MuonCollection for label: " << muonTag_;
     return false;
   }
 
   // for Track jet collections
   edm::Handle<reco::PFJetCollection> jetColl;
   iEvent.getByToken(jetsToken_, jetColl);
   std::vector<TLorentzVector> list_jets;
 
   if (jetColl.isValid()) {
     for (const auto& jets : *jetColl) {
       if (jets.pt() < minPt_Jets_)
         continue;
       if (std::fabs(jets.eta()) > maxEta_Jets_)
         continue;
       TLorentzVector lv_jets;  // lv_bJets;
       lv_jets.SetPtEtaPhiE(jets.pt(), jets.eta(), jets.phi(), jets.energy());
       list_jets.push_back(lv_jets);
     }
     lj = list_jets.size();
   }
 
   edm::Handle<reco::JetTagCollection> bTagHandle;
   iEvent.getByToken(bjetsToken_, bTagHandle);
   const reco::JetTagCollection& bTags = *(bTagHandle.product());
   std::vector<TLorentzVector> list_bjets;
 
   if (!bTags.empty()) {
     for (unsigned bj = 0; bj != bTags.size(); ++bj) {
       TLorentzVector lv_bjets;
       lv_bjets.SetPtEtaPhiE(
           bTags[bj].first->pt(), bTags[bj].first->eta(), bTags[bj].first->phi(), bTags[bj].first->energy());
       if ((bTags[bj].second > btagFactor_) && (lv_bjets.Pt() > minPt_Jets_))
         list_bjets.push_back(lv_bjets);
     }
     lbj = list_bjets.size();
   }
 
   // for MET collection
   edm::Handle<reco::PFMETCollection> pfColl;
   iEvent.getByToken(pfmetToken_, pfColl);
   if (EventCategory(le, lm, lj, lbj) == 11) {  // dilepton- ele ele
     if (list_ele[0].Pt() < minPtHighestEle_)
       return false;
     if ((list_ele[0].Pt() < list_mu[0].Pt()) || (list_ele[1].Pt() < list_mu[0].Pt()))
       return false;
     if (chrgeList_ele[0] + chrgeList_ele[1] != 0)
       return false;
     if (pfColl.isValid()) {
       double mt1 = getMt(list_ele[0], pfColl->front());
       double mt2 = getMt(list_ele[1], pfColl->front());
       double mt = mt1 + mt2;
       if (mt < 2 * minMet_ || mt > 2 * maxMet_)
         return false;
     } else {
       edm::LogError("ttbarEventSelector")
           << "Error >> Failed to get PFMETCollection in dilepton ele-ele channel for label: " << pfmetTag_;
       return false;
     }
   } else if (EventCategory(le, lm, lj, lbj) == 12) {  // dilepton - mu mu
     if (list_mu[0].Pt() < minPtHighestMu_)
       return false;
     if ((list_mu[0].Pt() < list_ele[0].Pt()) || (list_mu[1].Pt() < list_ele[0].Pt()))
       return false;
     if (chrgeList_mu[0] + chrgeList_mu[1] != 0)
       return false;
     if (pfColl.isValid()) {
       double mt1 = getMt(list_mu[0], pfColl->front());
       double mt2 = getMt(list_mu[1], pfColl->front());
       double mt = mt1 + mt2;
       if (mt < 2 * minMet_ || mt > 2 * maxMet_)
         return false;
     } else {
       edm::LogError("ttbarEventSelector")
           << "Error >> Failed to get PFMETCollection in dilepton mu-mu channel for label: " << pfmetTag_;
       return false;
     }
   } else if (EventCategory(le, lm, lj, lbj) == 13) {  // dilepton - ele mu
     if ((list_mu[0].Pt() < list_ele[1].Pt()) || (list_ele[0].Pt() < list_mu[1].Pt()))
       return false;
     if ((list_mu[0].Pt() < minPtHighestMu_) || (list_ele[0].Pt() < minPtHighestEle_))
       return false;
     if (chrgeList_mu[0] + chrgeList_ele[0] != 0)
       return false;
     if (pfColl.isValid()) {
       double mt1 = getMt(list_mu[0], pfColl->front());
       double mt2 = getMt(list_ele[0], pfColl->front());
       double mt = mt1 + mt2;
       if (mt < 2 * minMet_ || mt > 2 * maxMet_)
         return false;
     } else {
       edm::LogError("ttbarEventSelector")
           << "Error >> Failed to get PFMETCollection in dilepton ele-mu channel for label: " << pfmetTag_;
       return false;
     }
   }
   if (EventCategory(le, lm, lj, lbj) == 21) {  // semilepton - ele or mu
     if (list_jets[0].Pt() < minPtHighest_Jets_)
       return false;
 
     if (list_ele[0].Pt() < minPtHighestEle_)
       return false;
     // Both should not be present at the same time
     if ((!list_ele.empty() && list_ele[0].Pt() > minPtHighestEle_) &&
         (!list_mu.empty() && list_mu[0].Pt() > minPtHighestMu_))
       return false;
 
     return true;
   }
   if (EventCategory(le, lm, lj, lbj) == 22) {  // semilepton - ele or mu
     if (list_jets[0].Pt() < minPtHighest_Jets_)
       return false;
 
     if (list_mu[0].Pt() < minPtHighestMu_)
       return false;
     // Both should not be present at the same time
     if ((!list_ele.empty() && list_ele[0].Pt() > minPtHighestEle_) &&
         (!list_mu.empty() && list_mu[0].Pt() > minPtHighestMu_))
       return false;
 
     return true;
   } else if (EventCategory(le, lm, lj, lbj) == 30) {
     if (list_jets[0].Pt() < minPtHighest_Jets_)
       return false;
     for (int i = 0; i < 4; i++) {
       TLorentzVector vjet1;
       for (int j = i + 1; j < 4; j++) {
         TLorentzVector vjet2;
         vjet1 = list_jets[i];
         vjet2 = list_jets[j];
         TLorentzVector vjet = vjet1 + vjet2;
         if (vjet.M() < minWmass_ || vjet.M() > maxWmass_)
           return false;
       }
     }
   } else if (EventCategory(le, lm, lj, lbj) == 50)
     return false;
 
   return false;
 }

◆ getMt()

double ttbarEventSelector::getMt	(	const TLorentzVector &	vlep,
		const reco::PFMET &	obj
	)

private

Definition at line 383 of file ttbarEventSelector.cc.

References funct::cos(), SiPixelRawToDigiRegional_cfi::deltaPhi, BTaggingMonitor_cfi::met, getGTfromDQMFile::obj, phi, funct::sin(), and mathSSE::sqrt().

Referenced by filter().

                                                                                  {
   double met = obj.et();
   double phi = obj.phi();
 
   TLorentzVector vmet;
   double metx = met * std::cos(phi);
   double mety = met * std::sin(phi);
   vmet.SetPxPyPzE(metx, mety, 0.0, met);
 
   // transverse mass
   TLorentzVector vw = vlep + vmet;
 
   return std::sqrt(2 * vlep.Et() * met * (1 - std::cos(deltaPhi(vlep.Phi(), phi))));
 }

Member Data Documentation

◆ bjetsTag_

const edm::InputTag ttbarEventSelector::bjetsTag_

private

Definition at line 48 of file ttbarEventSelector.h.

◆ bjetsToken_

const edm::EDGetTokenT<reco::JetTagCollection> ttbarEventSelector::bjetsToken_

private

Definition at line 54 of file ttbarEventSelector.h.

Referenced by filter().

◆ bsTag_

const edm::InputTag ttbarEventSelector::bsTag_

private

Definition at line 51 of file ttbarEventSelector.h.

Referenced by filter().

◆ bsToken_

const edm::EDGetTokenT<reco::BeamSpot> ttbarEventSelector::bsToken_

private

Definition at line 57 of file ttbarEventSelector.h.

Referenced by filter().

◆ btagFactor_

const double ttbarEventSelector::btagFactor_

private

Definition at line 78 of file ttbarEventSelector.h.

Referenced by filter().

◆ electronTag_

const edm::InputTag ttbarEventSelector::electronTag_

private

Definition at line 46 of file ttbarEventSelector.h.

Referenced by filter().

◆ electronToken_

const edm::EDGetTokenT<reco::GsfElectronCollection> ttbarEventSelector::electronToken_

private

Definition at line 52 of file ttbarEventSelector.h.

Referenced by filter().

◆ jetsTag_

const edm::InputTag ttbarEventSelector::jetsTag_

private

Definition at line 47 of file ttbarEventSelector.h.

◆ jetsToken_

const edm::EDGetTokenT<reco::PFJetCollection> ttbarEventSelector::jetsToken_

private

Definition at line 53 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxD0_

const double ttbarEventSelector::maxD0_

private

Definition at line 81 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxDeltaEtaInEB_

const double ttbarEventSelector::maxDeltaEtaInEB_

private

Definition at line 63 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxDeltaEtaInEE_

const double ttbarEventSelector::maxDeltaEtaInEE_

private

Definition at line 67 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxDeltaPhiInEB_

const double ttbarEventSelector::maxDeltaPhiInEB_

private

Definition at line 62 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxDeltaPhiInEE_

const double ttbarEventSelector::maxDeltaPhiInEE_

private

Definition at line 66 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxDz_

const double ttbarEventSelector::maxDz_

private

Definition at line 82 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxEta_Jets_

const double ttbarEventSelector::maxEta_Jets_

private

Definition at line 76 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxEtaEle_

const double ttbarEventSelector::maxEtaEle_

private

Definition at line 59 of file ttbarEventSelector.h.

◆ maxEtaHighest_Jets_

const double ttbarEventSelector::maxEtaHighest_Jets_

private

Definition at line 75 of file ttbarEventSelector.h.

◆ maxEtaMu_

const double ttbarEventSelector::maxEtaMu_

private

Definition at line 60 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxHOEEB_

const double ttbarEventSelector::maxHOEEB_

private

Definition at line 64 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxHOEEE_

const double ttbarEventSelector::maxHOEEE_

private

Definition at line 68 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxInvMass_

const double ttbarEventSelector::maxInvMass_

private

Definition at line 92 of file ttbarEventSelector.h.

◆ maxIsoEle_

const double ttbarEventSelector::maxIsoEle_

private

Definition at line 85 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxIsoMu_

const double ttbarEventSelector::maxIsoMu_

private

Definition at line 86 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxMet_

const double ttbarEventSelector::maxMet_

private

Definition at line 94 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxNormChi2_

const double ttbarEventSelector::maxNormChi2_

private

Definition at line 80 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxSigmaiEiEEB_

const double ttbarEventSelector::maxSigmaiEiEEB_

private

Definition at line 65 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxSigmaiEiEEE_

const double ttbarEventSelector::maxSigmaiEiEEE_

private

Definition at line 69 of file ttbarEventSelector.h.

Referenced by filter().

◆ maxWmass_

const double ttbarEventSelector::maxWmass_

private

Definition at line 96 of file ttbarEventSelector.h.

Referenced by filter().

◆ minChambers_

const double ttbarEventSelector::minChambers_

private

Definition at line 71 of file ttbarEventSelector.h.

Referenced by filter().

◆ minInvMass_

const double ttbarEventSelector::minInvMass_

private

Definition at line 91 of file ttbarEventSelector.h.

◆ minMatchedStations_

const double ttbarEventSelector::minMatchedStations_

private

Definition at line 73 of file ttbarEventSelector.h.

Referenced by filter().

◆ minMatches_

const double ttbarEventSelector::minMatches_

private

Definition at line 72 of file ttbarEventSelector.h.

Referenced by filter().

◆ minMet_

const double ttbarEventSelector::minMet_

private

Definition at line 93 of file ttbarEventSelector.h.

◆ minPixelHits_

const int ttbarEventSelector::minPixelHits_

private

Definition at line 83 of file ttbarEventSelector.h.

Referenced by filter().

◆ minPt_

const double ttbarEventSelector::minPt_

private

Definition at line 61 of file ttbarEventSelector.h.

Referenced by filter().

◆ minPt_Jets_

const double ttbarEventSelector::minPt_Jets_

private

Definition at line 90 of file ttbarEventSelector.h.

Referenced by filter().

◆ minPtHighest_Jets_

const double ttbarEventSelector::minPtHighest_Jets_

private

Definition at line 89 of file ttbarEventSelector.h.

Referenced by filter().

◆ minPtHighestEle_

const double ttbarEventSelector::minPtHighestEle_

private

Definition at line 88 of file ttbarEventSelector.h.

Referenced by filter().

◆ minPtHighestMu_

const double ttbarEventSelector::minPtHighestMu_

private

Definition at line 87 of file ttbarEventSelector.h.

Referenced by filter().

◆ minStripHits_

const int ttbarEventSelector::minStripHits_

private

Definition at line 84 of file ttbarEventSelector.h.

Referenced by filter().

◆ minWmass_

const double ttbarEventSelector::minWmass_

private

Definition at line 95 of file ttbarEventSelector.h.

Referenced by filter().

◆ muonTag_

const edm::InputTag ttbarEventSelector::muonTag_

private

Definition at line 50 of file ttbarEventSelector.h.

Referenced by filter().

◆ muonToken_

const edm::EDGetTokenT<reco::MuonCollection> ttbarEventSelector::muonToken_

private

Definition at line 56 of file ttbarEventSelector.h.

Referenced by filter().

◆ pfmetTag_

const edm::InputTag ttbarEventSelector::pfmetTag_

private

Definition at line 49 of file ttbarEventSelector.h.

Referenced by filter().

◆ pfmetToken_

const edm::EDGetTokenT<reco::PFMETCollection> ttbarEventSelector::pfmetToken_

private

Definition at line 55 of file ttbarEventSelector.h.

Referenced by filter().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ttbarEventSelector()

Member Function Documentation

◆ EventCategory()

◆ filter()

◆ getMt()

Member Data Documentation

◆ bjetsTag_

◆ bjetsToken_

◆ bsTag_

◆ bsToken_

◆ btagFactor_

◆ electronTag_

◆ electronToken_

◆ jetsTag_

◆ jetsToken_

◆ maxD0_

◆ maxDeltaEtaInEB_

◆ maxDeltaEtaInEE_

◆ maxDeltaPhiInEB_

◆ maxDeltaPhiInEE_

◆ maxDz_

◆ maxEta_Jets_

◆ maxEtaEle_

◆ maxEtaHighest_Jets_

◆ maxEtaMu_

◆ maxHOEEB_

◆ maxHOEEE_

◆ maxInvMass_

◆ maxIsoEle_

◆ maxIsoMu_

◆ maxMet_

◆ maxNormChi2_

◆ maxSigmaiEiEEB_

◆ maxSigmaiEiEEE_

◆ maxWmass_

◆ minChambers_

◆ minInvMass_

◆ minMatchedStations_

◆ minMatches_

◆ minMet_

◆ minPixelHits_

◆ minPt_

◆ minPt_Jets_

◆ minPtHighest_Jets_

◆ minPtHighestEle_

◆ minPtHighestMu_

◆ minStripHits_

◆ minWmass_

◆ muonTag_

◆ muonToken_

◆ pfmetTag_

◆ pfmetToken_