Inheritance diagram for WMuNuSelector:

Public Member Functions
void	beginJob () override

void	endJob () override

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

void	init_histograms ()

	WMuNuSelector (const edm::ParameterSet &)

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 Attributes
double	acopCut_

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

double	dxyCut_

double	eJetMin_

double	etaCut_

std::map< std::string, TH1D * >	h1_

std::map< std::string, TH2D * >	h2_

bool	isAlsoTrackerMuon_

bool	isCombinedIso_

double	isoCut03_

bool	isRelativeIso_

edm::EDGetTokenT< edm::View< reco::Jet > >	jetToken_

double	metMax_

double	metMin_

double	mtMax_

double	mtMin_

edm::EDGetTokenT< edm::View< reco::Muon > >	muonToken_

const std::string	muonTrig_

double	nacop

double	nall

double	ncharge

double	nid

double	niso

int	nJetMax_

double	nkin

double	nmass

double	normalizedChi2Cut_

double	npresel

double	nsel

double	ntrig

bool	plotHistograms_

double	ptCut_

double	ptThrForZ1_

double	ptThrForZ2_

int	selectByCharge_

int	trackerHitsCut_

edm::EDGetTokenT< edm::TriggerResults >	trigToken_

edm::EDGetTokenT< reco::WMuNuCandidateCollection >	WMuNuCollectionToken_

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::ProducerBase
ProducesCollector	producesCollector ()

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 35 of file WMuNuSelector.cc.

Constructor & Destructor Documentation

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

Definition at line 102 of file WMuNuSelector.cc.

     :  // Fast selection (no histograms)
       plotHistograms_(cfg.getUntrackedParameter<bool>("plotHistograms", true)),
 
       // Input collections
       trigToken_(consumes<TriggerResults>(
           cfg.getUntrackedParameter<edm::InputTag>("TrigTag", edm::InputTag("TriggerResults::HLT")))),
       muonToken_(consumes<View<Muon> >(cfg.getUntrackedParameter<edm::InputTag>("MuonTag", edm::InputTag("muons")))),
       jetToken_(
           consumes<View<Jet> >(cfg.getUntrackedParameter<edm::InputTag>("JetTag", edm::InputTag("sisCone5CaloJets")))),
       beamSpotToken_(consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"))),
       WMuNuCollectionToken_(consumes<reco::WMuNuCandidateCollection>(
           cfg.getUntrackedParameter<edm::InputTag>("WMuNuCollectionTag", edm::InputTag("WMuNus")))),
 
       // Preselection cuts
       muonTrig_(cfg.getUntrackedParameter<std::string>("MuonTrig", "HLT_Mu9")),
       ptThrForZ1_(cfg.getUntrackedParameter<double>("PtThrForZ1", 20.)),
       ptThrForZ2_(cfg.getUntrackedParameter<double>("PtThrForZ2", 10.)),
       eJetMin_(cfg.getUntrackedParameter<double>("EJetMin", 999999.)),
       nJetMax_(cfg.getUntrackedParameter<int>("NJetMax", 999999)),
 
       // Main cuts
       ptCut_(cfg.getUntrackedParameter<double>("PtCut", 25.)),
       etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.1)),
       isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso", true)),
       isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso", false)),
       isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03", 0.1)),
       mtMin_(cfg.getUntrackedParameter<double>("MtMin", 50.)),
       mtMax_(cfg.getUntrackedParameter<double>("MtMax", 200.)),
       metMin_(cfg.getUntrackedParameter<double>("MetMin", -999999.)),
       metMax_(cfg.getUntrackedParameter<double>("MetMax", 999999.)),
       acopCut_(cfg.getUntrackedParameter<double>("AcopCut", 2.)),
 
       // Muon quality cuts
       dxyCut_(cfg.getUntrackedParameter<double>("DxyCut", 0.2)),
       normalizedChi2Cut_(cfg.getUntrackedParameter<double>("NormalizedChi2Cut", 10.)),
       trackerHitsCut_(cfg.getUntrackedParameter<int>("TrackerHitsCut", 11)),
       isAlsoTrackerMuon_(cfg.getUntrackedParameter<bool>("IsAlsoTrackerMuon", true)),
 
       // W+/W- Selection
       selectByCharge_(cfg.getUntrackedParameter<int>("SelectByCharge", 0)) {}

Member Function Documentation

void WMuNuSelector::beginJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 144 of file WMuNuSelector.cc.

References h1_, h2_, M_PI, TFileService::make(), nacop, nall, ncharge, nid, niso, nkin, npresel, nsel, ntrig, and plotHistograms_.

                              {
   nall = 0;
   ntrig = 0;
   npresel = 0;
   ncharge = 0;
   nkin = 0;
   nid = 0;
   nacop = 0;
   niso = 0;
   nsel = 0;
 
   if (plotHistograms_) {
     edm::Service<TFileService> fs;
     h1_["hNWCand_sel"] =
         fs->make<TH1D>("NWCand_sel", "Nb. of WCandidates passing pre-selection (ordered by pt)", 10, 0., 10.);
     h1_["hPtMu_sel"] = fs->make<TH1D>("ptMu_sel", "Pt mu", 100, 0., 100.);
     h1_["hEtaMu_sel"] = fs->make<TH1D>("etaMu_sel", "Eta mu", 50, -2.5, 2.5);
     h1_["hd0_sel"] = fs->make<TH1D>("d0_sel", "Impact parameter", 1000, -1., 1.);
     h1_["hNHits_sel"] = fs->make<TH1D>("NumberOfValidHits_sel", "Number of Hits in Silicon", 100, 0., 100.);
     h1_["hNormChi2_sel"] = fs->make<TH1D>("NormChi2_sel", "Chi2/ndof of global track", 1000, 0., 50.);
     h1_["hTracker_sel"] = fs->make<TH1D>("isTrackerMuon_sel", "is Tracker Muon?", 2, 0., 2.);
     h1_["hMET_sel"] = fs->make<TH1D>("MET_sel", "Missing Transverse Energy (GeV)", 300, 0., 300.);
     h1_["hTMass_sel"] = fs->make<TH1D>("TMass_sel", "Rec. Transverse Mass (GeV)", 300, 0., 300.);
     h1_["hAcop_sel"] = fs->make<TH1D>("Acop_sel", "Mu-MET acoplanarity", 50, 0., M_PI);
     h1_["hPtSum_sel"] = fs->make<TH1D>("ptSum_sel", "Track Isolation, Sum pT (GeV)", 200, 0., 100.);
     h1_["hPtSumN_sel"] = fs->make<TH1D>("ptSumN_sel", "Track Isolation, Sum pT/pT", 1000, 0., 10.);
     h1_["hCal_sel"] = fs->make<TH1D>("Cal_sel", "Calorimetric isolation, HCAL+ECAL (GeV)", 200, 0., 100.);
     h1_["hIsoTotN_sel"] = fs->make<TH1D>("isoTotN_sel", "(Sum pT + Cal)/pT", 1000, 0., 10.);
     h1_["hIsoTot_sel"] = fs->make<TH1D>("isoTot_sel", "(Sum pT + Cal)", 200, 0., 100.);
     h2_["hTMass_PtSum_inclusive"] = fs->make<TH2D>(
         "TMass_PtSum_inclusive", "Rec. Transverse Mass (GeV) vs Sum pT (GeV)", 200, 0., 100., 300, 0., 300.);
     h2_["hTMass_PtSumNorm_inclusive"] = fs->make<TH2D>(
         "TMass_PtSumNorm_inclusive", "Rec. Transverse Mass (GeV) vs Sum Pt / Pt", 1000, 0, 10, 300, 0, 300);
     h2_["hTMass_TotIsoNorm_inclusive"] = fs->make<TH2D>(
         "TMass_TotIsoNorm_inclusive", "Rec. Transverse Mass (GeV) vs (Sum Pt + Cal)/Pt", 10000, 0, 10, 200, 0, 200);
     h2_["hMET_PtSum_inclusive"] = fs->make<TH2D>(
         "MET_PtSum_inclusive", "Missing Transverse Energy (GeV) vs Sum Pt (GeV)", 200, 0., 100., 300, 0., 300.);
     h2_["hMET_PtSumNorm_inclusive"] = fs->make<TH2D>(
         "MET_PtSumNorm_inclusive", "Missing Transverse Energy (GeV) vs Sum Pt/Pt", 1000, 0, 10, 300, 0, 300);
     h2_["hMET_TotIsoNorm_inclusive"] = fs->make<TH2D>(
         "MET_TotIsoNorm_inclusive", "Missing Transverse Energy (GeV) vs (SumPt + Cal)/Pt", 10000, 0, 10, 200, 0, 200);
   }
 }

void WMuNuSelector::endJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 188 of file WMuNuSelector.cc.

References Vispa.Plugins.EdmBrowser.EdmDataAccessor::all(), runTauDisplay::eid, nacop, nall, ncharge, nid, niso, nkin, npresel, nsel, ntrig, selectByCharge_, and mathSSE::sqrt().

                            {
   double all = nall;
   double epresel = npresel / all;
   double etrig = ntrig / all;
   double ekin = nkin / all;
   double eid = nid / all;
   double eacop = nacop / all;
   double eiso = niso / all;
   double esel = nsel / all;
 
   LogVerbatim("") << "\n>>>>>> W SELECTION SUMMARY BEGIN >>>>>>>>>>>>>>>";
   LogVerbatim("") << "Total number of events analyzed: " << nall << " [events]";
   LogVerbatim("") << "Total number of events triggered: " << ntrig << " [events]";
   LogVerbatim("") << "Total number of events pre-selected: " << npresel << " [events]";
   LogVerbatim("") << "Total number of events after kinematic cuts: " << nkin << " [events]";
   LogVerbatim("") << "Total number of events after Muon ID cuts: " << nid << " [events]";
   LogVerbatim("") << "Total number of events after Acop cut: " << nacop << " [events]";
   LogVerbatim("") << "Total number of events after iso cut: " << niso << " [events]";
   LogVerbatim("") << "Total number of events selected: " << nsel << " [events]";
   LogVerbatim("") << "Efficiencies:";
   LogVerbatim("") << "Trigger Efficiency:                   "
                   << "(" << setprecision(4) << etrig * 100. << " +/- " << setprecision(2)
                   << sqrt(etrig * (1 - etrig) / all) * 100. << ")%";
   LogVerbatim("") << "Pre-Selection Efficiency:             "
                   << "(" << setprecision(4) << epresel * 100. << " +/- " << setprecision(2)
                   << sqrt(epresel * (1 - epresel) / all) * 100. << ")%";
   LogVerbatim("") << "Pt, Eta Selection Efficiency:         "
                   << "(" << setprecision(4) << ekin * 100. << " +/- " << setprecision(2)
                   << sqrt(ekin * (1 - ekin) / all) * 100. << ")%";
   LogVerbatim("") << "MuonID Efficiency:                    "
                   << "(" << setprecision(4) << eid * 100. << " +/- " << setprecision(2)
                   << sqrt(eid * (1 - eid) / all) * 100. << ")%";
   LogVerbatim("") << "Acop Efficiency:                      "
                   << "(" << setprecision(4) << eacop * 100. << " +/- " << setprecision(2)
                   << sqrt(eacop * (1 - eacop) / all) * 100. << ")%";
   LogVerbatim("") << "Iso Efficiency:                       "
                   << "(" << setprecision(4) << eiso * 100. << " +/- " << setprecision(2)
                   << sqrt(eiso * (1 - eiso) / all) * 100. << ")%";
   LogVerbatim("") << "Selection Efficiency:             "
                   << "(" << setprecision(4) << esel * 100. << " +/- " << setprecision(2)
                   << sqrt(esel * (1 - esel) / nall) * 100. << ")%";
 
   if (fabs(selectByCharge_) == 1) {
     esel = nsel / ncharge;
     LogVerbatim("") << "\n>>>>>> W+(-) SELECTION >>>>>>>>>>>>>>>";
     LogVerbatim("") << "Total number of W+(-) events pre-selected: " << ncharge << " [events]";
     LogVerbatim("") << "Total number of events selected: " << nsel << " [events]";
     LogVerbatim("") << "Selection Efficiency only W+(-): "
                     << "(" << setprecision(4) << esel * 100. << " +/- " << setprecision(2)
                     << sqrt(esel * (1 - esel) / ncharge) * 100. << ")%";
   }
   LogVerbatim("") << ">>>>>> W SELECTION SUMMARY END   >>>>>>>>>>>>>>>\n";
 }

bool WMuNuSelector::filter	(	edm::Event &	ev,
		const edm::EventSetup &
	)

override

Definition at line 242 of file WMuNuSelector.cc.

References edm::HLTGlobalStatus::accept(), reco::WMuNuCandidate::acop(), acopCut_, beamSpotToken_, reco::LeafCandidate::charge(), DEFINE_FWK_MODULE, PVValHelper::dxy, eJetMin_, reco::MuonIsolation::emEt, reco::WMuNuCandidate::eT(), reco::LeafCandidate::et(), PVValHelper::eta, reco::LeafCandidate::eta(), etaCut_, edm::Event::getByToken(), reco::WMuNuCandidate::getMuon(), reco::WMuNuCandidate::getNeutrino(), muon::GlobalMuonPromptTight, reco::Muon::globalTrack(), h1_, h2_, reco::MuonIsolation::hadEt, mps_fire::i, isCombinedIso_, reco::Muon::isGlobalMuon(), muon::isGoodMuon(), isoCut03_, reco::Muon::isolationR03(), isRelativeIso_, reco::Muon::isTrackerMuon(), metsig::jet, jetfilter_cfi::jetCollection, jetToken_, LogTrace, reco::WMuNuCandidate::massT(), BTaggingMonitor_cfi::met, metMax_, mtMax_, amptDefaultParameters_cff::mu, pdwgLeptonRecoSkim_cfi::muonCollection, muonToken_, muonTrig_, nacop, nall, ncharge, nid, niso, nJetMax_, BTaggingMonitoring_cff::njets, nkin, TrackCollections2monitor_cff::normalizedChi2, normalizedChi2Cut_, npresel, nsel, ntrig, plotHistograms_, reco::BeamSpot::position(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), ptCut_, ptThrForZ1_, ptThrForZ2_, reco::LeafCandidate::px(), reco::LeafCandidate::py(), selectByCharge_, reco::MuonIsolation::sumPt, trackerHitsCut_, edm::TriggerNames::triggerIndex(), L1TEGammaOffline_cfi::triggerNames, edm::Event::triggerNames(), triggerResults, trigToken_, and WMuNuCollectionToken_.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

                                                        {
   nall++;
 
   // Repeat Pre-Selection Cuts just in case...
   // Muon collection
   Handle<View<Muon> > muonCollection;
   if (!ev.getByToken(muonToken_, muonCollection)) {
     LogError("") << ">>> Muon collection does not exist !!!";
     return false;
   }
   unsigned int muonCollectionSize = muonCollection->size();
 
   // Trigger
   Handle<TriggerResults> triggerResults;
   if (!ev.getByToken(trigToken_, triggerResults)) {
     LogError("") << ">>> TRIGGER collection does not exist !!!";
     return false;
   }
   const edm::TriggerNames& triggerNames = ev.triggerNames(*triggerResults);
   bool trigger_fired = false;
   int itrig1 = triggerNames.triggerIndex(muonTrig_);
   if (triggerResults->accept(itrig1))
     trigger_fired = true;
   LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << muonTrig_ << ")";
 
   // Loop to reject/control Z->mumu is done separately
   unsigned int nmuonsForZ1 = 0;
   unsigned int nmuonsForZ2 = 0;
   for (unsigned int i = 0; i < muonCollectionSize; i++) {
     const Muon& mu = muonCollection->at(i);
     if (!mu.isGlobalMuon())
       continue;
     double pt = mu.pt();
     if (pt > ptThrForZ1_)
       nmuonsForZ1++;
     if (pt > ptThrForZ2_)
       nmuonsForZ2++;
   }
   LogTrace("") << "> Z rejection: muons above " << ptThrForZ1_ << " [GeV]: " << nmuonsForZ1;
   LogTrace("") << "> Z rejection: muons above " << ptThrForZ2_ << " [GeV]: " << nmuonsForZ2;
 
   // Jet collection
   Handle<View<Jet> > jetCollection;
   if (!ev.getByToken(jetToken_, jetCollection)) {
     LogError("") << ">>> JET collection does not exist !!!";
     return false;
   }
   unsigned int jetCollectionSize = jetCollection->size();
   int njets = 0;
   for (unsigned int i = 0; i < jetCollectionSize; i++) {
     const Jet& jet = jetCollection->at(i);
     if (jet.et() > eJetMin_)
       njets++;
   }
   LogTrace("") << ">>> Total number of jets: " << jetCollectionSize;
   LogTrace("") << ">>> Number of jets above " << eJetMin_ << " [GeV]: " << njets;
 
   // Beam spot
   Handle<reco::BeamSpot> beamSpotHandle;
   if (!ev.getByToken(beamSpotToken_, beamSpotHandle)) {
     LogTrace("") << ">>> No beam spot found !!!";
     return false;
   }
 
   // Get WMuNu candidates from file:
 
   Handle<reco::WMuNuCandidateCollection> WMuNuCollection;
   if (!ev.getByToken(WMuNuCollectionToken_, WMuNuCollection)) {
     LogTrace("") << ">>> WMuNu not found !!!";
     return false;
   }
 
   if (WMuNuCollection->empty()) {
     LogTrace("") << "No WMuNu Candidates in the Event!";
     return false;
   }
   if (WMuNuCollection->size() > 1) {
     LogTrace("") << "This event contains more than one W Candidate";
   }
 
   // W->mu nu selection criteria
 
   LogTrace("") << "> WMuNu Candidate with: ";
   const WMuNuCandidate& WMuNu = WMuNuCollection->at(0);
   // WMuNuCandidates are ordered by Pt!
   // The Inclusive Selection WMuNu Candidate is the first one
 
   const reco::Muon& mu = WMuNu.getMuon();
   const reco::MET& met = WMuNu.getNeutrino();
   if (plotHistograms_) {
     h1_["hNWCand_sel"]->Fill(WMuNuCollection->size());
   }
 
   // Preselection cuts:
 
   if (!trigger_fired) {
     LogTrace("") << "Event did not fire the Trigger";
     return false;
   }
   ntrig++;
 
   if (nmuonsForZ1 >= 1 && nmuonsForZ2 >= 2) {
     LogTrace("") << "Z Candidate!!";
     return false;
   }
   if (njets > nJetMax_) {
     LogTrace("") << "NJets > threshold";
     return false;
   }
 
   npresel++;
 
   // Select Ws by charge:
 
   if (selectByCharge_ * WMuNu.charge() == -1) {
     ncharge++;
   }
 
   // W->mu nu selection criteria
 
   if (!mu.isGlobalMuon())
     return false;
 
   reco::TrackRef gm = mu.globalTrack();
   //reco::TrackRef tk = mu.innerTrack();
 
   // Pt,eta cuts
   double pt = mu.pt();
   double eta = mu.eta();
   LogTrace("") << "\t... Muon pt, eta: " << pt << " [GeV], " << eta;
   if (plotHistograms_) {
     h1_["hPtMu_sel"]->Fill(pt);
   }
   if (pt < ptCut_)
     return false;
   if (plotHistograms_) {
     h1_["hEtaMu_sel"]->Fill(eta);
   }
   if (fabs(eta) > etaCut_)
     return false;
 
   nkin++;
 
   // d0, chi2, nhits quality cuts
   double dxy = gm->dxy(beamSpotHandle->position());
   double normalizedChi2 = gm->normalizedChi2();
   LogTrace("") << "Im here" << endl;
   double trackerHits = gm->hitPattern().numberOfValidTrackerHits();
   LogTrace("") << "\t... Muon dxy, normalizedChi2, trackerHits, isTrackerMuon?: " << dxy << " [cm], " << normalizedChi2
                << ", " << trackerHits << ", " << mu.isTrackerMuon();
 
   if (plotHistograms_) {
     h1_["hd0_sel"]->Fill(dxy);
   }
   if (!muon::isGoodMuon(mu, muon::GlobalMuonPromptTight))
     return false;
   if (plotHistograms_) {
     h1_["hNormChi2_sel"]->Fill(normalizedChi2);
   }
   if (normalizedChi2 > normalizedChi2Cut_)
     return false;
   if (plotHistograms_) {
     h1_["hNHits_sel"]->Fill(trackerHits);
   }
   if (trackerHits < trackerHitsCut_)
     return false;
   if (plotHistograms_) {
     h1_["hTracker_sel"]->Fill(mu.isTrackerMuon());
   }
   if (!mu.isTrackerMuon())
     return false;
 
   nid++;
 
   // Acoplanarity cuts
   double acop = WMuNu.acop();
   LogTrace("") << "\t... acoplanarity: " << acop;
 
   // Isolation cuts
   double SumPt = mu.isolationR03().sumPt;
   double isovar = SumPt;
   double Cal = mu.isolationR03().emEt + mu.isolationR03().hadEt;
   if (isCombinedIso_)
     isovar += Cal;
   if (plotHistograms_) {
     h1_["hPtSum_sel"]->Fill(SumPt);
     h1_["hPtSumN_sel"]->Fill(SumPt / pt);
     h1_["hCal_sel"]->Fill(Cal);
     h1_["hIsoTot_sel"]->Fill((SumPt + Cal));
     h1_["hIsoTotN_sel"]->Fill((SumPt + Cal) / pt);
   }
 
   if (isRelativeIso_)
     isovar /= pt;
   bool iso = (isovar <= isoCut03_);
   LogTrace("") << "\t... isolation value" << isovar << ", isolated? " << iso;
 
   double met_et = met.pt();
   LogTrace("") << "\t... Met  pt: " << WMuNu.getNeutrino().pt() << "[GeV]";
 
   double massT = WMuNu.massT();
   double w_et = WMuNu.eT();
 
   LogTrace("") << "\t... W mass, W_et, W_px, W_py: " << massT << ", " << w_et << ", " << WMuNu.px() << ", "
                << WMuNu.py() << " [GeV]";
 
   // Plot 2D Histograms before final cuts
   if (plotHistograms_ && acop < acopCut_) {
     h2_["hTMass_PtSum_inclusive"]->Fill(SumPt, massT);
     h2_["hTMass_PtSumNorm_inclusive"]->Fill(SumPt / pt, massT);
     h2_["hTMass_TotIsoNorm_inclusive"]->Fill((SumPt + Cal) / pt, massT);
     h2_["hMET_PtSum_inclusive"]->Fill(SumPt, met_et);
     h2_["hMET_PtSumNorm_inclusive"]->Fill(SumPt / pt, met_et);
     h2_["hMET_TotIsoNorm_inclusive"]->Fill((SumPt + Cal) / pt, met_et);
   }
 
   if (!iso)
     return false;
 
   niso++;
 
   if (plotHistograms_) {
     h1_["hAcop_sel"]->Fill(acop);
   }
   if (acop >= acopCut_)
     return false;
 
   nacop++;
 
   if (plotHistograms_) {
     h1_["hMET_sel"]->Fill(met_et);
     h1_["hTMass_sel"]->Fill(massT);
   }
 
   if (massT <= mtMin_ || massT >= mtMax_)
     return false;
   if (met_et <= metMin_ || met_et >= metMax_)
     return false;
 
   LogTrace("") << ">>>> Event ACCEPTED";
   nsel++;
 
   // (To be continued ;-) )
 
   return true;
 }

void WMuNuSelector::init_histograms ( )

Member Data Documentation

double WMuNuSelector::acopCut_

private

Definition at line 64 of file WMuNuSelector.cc.

Referenced by filter().

edm::EDGetTokenT<reco::BeamSpot> WMuNuSelector::beamSpotToken_

private

Definition at line 48 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::dxyCut_

private

Definition at line 66 of file WMuNuSelector.cc.

double WMuNuSelector::eJetMin_

private

Definition at line 53 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::etaCut_

private

Definition at line 56 of file WMuNuSelector.cc.

Referenced by filter().

std::map<std::string, TH1D*> WMuNuSelector::h1_

private

Definition at line 79 of file WMuNuSelector.cc.

Referenced by beginJob(), and filter().

std::map<std::string, TH2D*> WMuNuSelector::h2_

private

Definition at line 80 of file WMuNuSelector.cc.

Referenced by beginJob(), and filter().

bool WMuNuSelector::isAlsoTrackerMuon_

private

Definition at line 69 of file WMuNuSelector.cc.

bool WMuNuSelector::isCombinedIso_

private

Definition at line 58 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::isoCut03_

private

Definition at line 59 of file WMuNuSelector.cc.

Referenced by filter().

bool WMuNuSelector::isRelativeIso_

private

Definition at line 57 of file WMuNuSelector.cc.

Referenced by filter().

edm::EDGetTokenT<edm::View<reco::Jet> > WMuNuSelector::jetToken_

private

Definition at line 47 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::metMax_

private

Definition at line 63 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::metMin_

private

Definition at line 62 of file WMuNuSelector.cc.

double WMuNuSelector::mtMax_

private

Definition at line 61 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::mtMin_

private

Definition at line 60 of file WMuNuSelector.cc.

edm::EDGetTokenT<edm::View<reco::Muon> > WMuNuSelector::muonToken_

private

Definition at line 46 of file WMuNuSelector.cc.

Referenced by filter().

const std::string WMuNuSelector::muonTrig_

private

Definition at line 50 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::nacop

private

Definition at line 77 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

double WMuNuSelector::nall

private

Definition at line 73 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

double WMuNuSelector::ncharge

private

Definition at line 76 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

double WMuNuSelector::nid

private

Definition at line 77 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

double WMuNuSelector::niso

private

Definition at line 77 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

int WMuNuSelector::nJetMax_

private

Definition at line 54 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::nkin

private

Definition at line 77 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

double WMuNuSelector::nmass

private

Definition at line 77 of file WMuNuSelector.cc.

double WMuNuSelector::normalizedChi2Cut_

private

Definition at line 67 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::npresel

private

Definition at line 74 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

double WMuNuSelector::nsel

private

Definition at line 75 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

double WMuNuSelector::ntrig

private

Definition at line 74 of file WMuNuSelector.cc.

Referenced by beginJob(), endJob(), and filter().

bool WMuNuSelector::plotHistograms_

private

Definition at line 44 of file WMuNuSelector.cc.

Referenced by beginJob(), and filter().

double WMuNuSelector::ptCut_

private

Definition at line 55 of file WMuNuSelector.cc.

Referenced by badGlobalMuonTaggerFWLite.BadGlobalMuonTagger::badMuons(), and filter().

double WMuNuSelector::ptThrForZ1_

private

Definition at line 51 of file WMuNuSelector.cc.

Referenced by filter().

double WMuNuSelector::ptThrForZ2_

private

Definition at line 52 of file WMuNuSelector.cc.

Referenced by filter().

int WMuNuSelector::selectByCharge_

private

Definition at line 71 of file WMuNuSelector.cc.

Referenced by endJob(), and filter().

int WMuNuSelector::trackerHitsCut_

private

Definition at line 68 of file WMuNuSelector.cc.

Referenced by filter().

edm::EDGetTokenT<edm::TriggerResults> WMuNuSelector::trigToken_

private

Definition at line 45 of file WMuNuSelector.cc.

Referenced by filter().

edm::EDGetTokenT<reco::WMuNuCandidateCollection> WMuNuSelector::WMuNuCollectionToken_

private

Definition at line 49 of file WMuNuSelector.cc.

Referenced by filter().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation