Inheritance diagram for WMuNuValidator:

Public Member Functions
void	beginJob () override

void	endJob () override

void	fill_histogram (const char *, const double &)

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

void	init_histograms ()

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

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

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_

bool	fastOption_

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

bool	isAlsoTrackerMuon_

bool	isCombinedIso_

double	isoCut03_

bool	isRelativeIso_

edm::InputTag	jetTag_

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

bool	metIncludesMuons_

double	metMax_

double	metMin_

edm::InputTag	metTag_

edm::EDGetTokenT< edm::View< reco::MET > >	metToken_

double	mtMax_

double	mtMin_

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

const std::string	muonTrig_

unsigned int	nall

unsigned int	nhlt

unsigned int	niso

int	nJetMax_

unsigned int	nmet

double	normalizedChi2Cut_

unsigned int	nrec

unsigned int	nsel

double	ptCut_

double	ptThrForZ1_

double	ptThrForZ2_

int	trackerHitsCut_

edm::EDGetTokenT< edm::TriggerResults >	trigToken_

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

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

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

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

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)

template<Transition B>
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 (std::string instanceName)

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

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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

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

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 25 of file WMuNuValidator.cc.

Constructor & Destructor Documentation

◆ WMuNuValidator()

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

Definition at line 97 of file WMuNuValidator.cc.

     :  // Fast selection (no histograms or book-keeping)
       fastOption_(cfg.getUntrackedParameter<bool>("FastOption", false)),
  
       // 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")))),
       metTag_(cfg.getUntrackedParameter<edm::InputTag>("METTag", edm::InputTag("met"))),
       metToken_(consumes<View<MET> >(metTag_)),
       metIncludesMuons_(cfg.getUntrackedParameter<bool>("METIncludesMuons", false)),
       jetTag_(cfg.getUntrackedParameter<edm::InputTag>("JetTag", edm::InputTag("sisCone5CaloJets"))),
       jetToken_(consumes<View<Jet> >(jetTag_)),
       beamSpotToken_(consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"))),
  
       // Main cuts
       muonTrig_(cfg.getUntrackedParameter<std::string>("MuonTrig", "HLT_Mu9")),
       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)),
  
       // Z rejection
       ptThrForZ1_(cfg.getUntrackedParameter<double>("PtThrForZ1", 20.)),
       ptThrForZ2_(cfg.getUntrackedParameter<double>("PtThrForZ2", 10.)),
  
       // Top rejection
       eJetMin_(cfg.getUntrackedParameter<double>("EJetMin", 999999.)),
       nJetMax_(cfg.getUntrackedParameter<int>("NJetMax", 999999)) {}

Member Function Documentation

◆ beginJob()

void WMuNuValidator::beginJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 139 of file WMuNuValidator.cc.

                               {
   nall = 0;
   nsel = 0;
  
   if (!fastOption_) {
     nrec = 0;
     niso = 0;
     nhlt = 0;
     nmet = 0;
     init_histograms();
   }
 }

References fastOption_, init_histograms(), nall, nhlt, niso, nmet, nrec, and nsel.

◆ endJob()

void WMuNuValidator::endJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 246 of file WMuNuValidator.cc.

                             {
   double all = nall;
   double esel = nsel / all;
   LogVerbatim("") << "\n>>>>>> W SELECTION SUMMARY BEGIN >>>>>>>>>>>>>>>";
   LogVerbatim("") << "Total numer of events analyzed: " << nall << " [events]";
   LogVerbatim("") << "Total numer of events selected: " << nsel << " [events]";
   LogVerbatim("") << "Overall efficiency:             "
                   << "(" << setprecision(4) << esel * 100. << " +/- " << setprecision(2)
                   << sqrt(esel * (1 - esel) / all) * 100. << ")%";
  
   if (!fastOption_) {
     double erec = nrec / all;
     double eiso = niso / all;
     double ehlt = nhlt / all;
     double emet = nmet / all;
  
     double num = nrec;
     double eff = erec;
     double err = sqrt(eff * (1 - eff) / all);
     LogVerbatim("") << "Passing Pt/Eta/Quality cuts:    " << num << " [events], (" << setprecision(4) << eff * 100.
                     << " +/- " << setprecision(2) << err * 100. << ")%";
  
     num = niso;
     eff = eiso;
     err = sqrt(eff * (1 - eff) / all);
     double effstep = 0.;
     double errstep = 0.;
     if (nrec > 0)
       effstep = eiso / erec;
     if (nrec > 0)
       errstep = sqrt(effstep * (1 - effstep) / nrec);
     LogVerbatim("") << "Passing isolation cuts:         " << num << " [events], (" << setprecision(4) << eff * 100.
                     << " +/- " << setprecision(2) << err * 100. << ")%, to previous step: (" << setprecision(4)
                     << effstep * 100. << " +/- " << setprecision(2) << errstep * 100. << ")%";
  
     num = nhlt;
     eff = ehlt;
     err = sqrt(eff * (1 - eff) / all);
     effstep = 0.;
     errstep = 0.;
     if (niso > 0)
       effstep = ehlt / eiso;
     if (niso > 0)
       errstep = sqrt(effstep * (1 - effstep) / niso);
     LogVerbatim("") << "Passing HLT criteria:           " << num << " [events], (" << setprecision(4) << eff * 100.
                     << " +/- " << setprecision(2) << err * 100. << ")%, to previous step: (" << setprecision(4)
                     << effstep * 100. << " +/- " << setprecision(2) << errstep * 100. << ")%";
  
     num = nmet;
     eff = emet;
     err = sqrt(eff * (1 - eff) / all);
     effstep = 0.;
     errstep = 0.;
     if (nhlt > 0)
       effstep = emet / ehlt;
     if (nhlt > 0)
       errstep = sqrt(effstep * (1 - effstep) / nhlt);
     LogVerbatim("") << "Passing MET/acoplanarity cuts:  " << num << " [events], (" << setprecision(4) << eff * 100.
                     << " +/- " << setprecision(2) << err * 100. << ")%, to previous step: (" << setprecision(4)
                     << effstep * 100. << " +/- " << setprecision(2) << errstep * 100. << ")%";
  
     num = nsel;
     eff = esel;
     err = sqrt(eff * (1 - eff) / all);
     effstep = 0.;
     errstep = 0.;
     if (nmet > 0)
       effstep = esel / emet;
     if (nmet > 0)
       errstep = sqrt(effstep * (1 - effstep) / nmet);
     LogVerbatim("") << "Passing Z/top rejection cuts:   " << num << " [events], (" << setprecision(4) << eff * 100.
                     << " +/- " << setprecision(2) << err * 100. << ")%, to previous step: (" << setprecision(4)
                     << effstep * 100. << " +/- " << setprecision(2) << errstep * 100. << ")%";
   }
  
   LogVerbatim("") << ">>>>>> W SELECTION SUMMARY END   >>>>>>>>>>>>>>>\n";
 }

References python.cmstools::all(), submitPVResolutionJobs::err, fastOption_, nall, nhlt, niso, nmet, nrec, nsel, EgammaValidation_cff::num, and mathSSE::sqrt().

◆ fill_histogram()

void WMuNuValidator::fill_histogram	(	const char *	name,
		const double &	var
	)

Definition at line 240 of file WMuNuValidator.cc.

                                                                        {
   if (fastOption_)
     return;
   h1_[name]->Fill(var);
 }

References fastOption_, h1_, Skims_PA_cff::name, and trigObjTnPSource_cfi::var.

Referenced by filter().

◆ filter()

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

overridevirtual

Implements edm::EDFilter.

Definition at line 324 of file WMuNuValidator.cc.

                                                         {
   // Reset global event selection flags
   bool rec_sel = false;
   bool iso_sel = false;
   bool hlt_sel = false;
   bool met_sel = false;
   bool all_sel = false;
  
   // Muon collection
   Handle<View<Muon> > muonCollection;
   if (!ev.getByToken(muonToken_, muonCollection)) {
     LogError("") << ">>> Muon collection does not exist !!!";
     return false;
   }
   unsigned int muonCollectionSize = muonCollection->size();
  
   // Beam spot
   Handle<reco::BeamSpot> beamSpotHandle;
   if (!ev.getByToken(beamSpotToken_, beamSpotHandle)) {
     LogTrace("") << ">>> No beam spot found !!!";
     return false;
   }
  
   // MET
   double met_px = 0.;
   double met_py = 0.;
   Handle<View<MET> > metCollection;
   if (!ev.getByToken(metToken_, metCollection)) {
     LogError("") << ">>> MET collection does not exist !!!";
     return false;
   }
   const MET& met = metCollection->at(0);
   met_px = met.px();
   met_py = met.py();
   if (!metIncludesMuons_) {
     for (unsigned int i = 0; i < muonCollectionSize; i++) {
       const Muon& mu = muonCollection->at(i);
       if (!mu.isGlobalMuon())
         continue;
       met_px -= mu.px();
       met_py -= mu.py();
     }
   }
   double met_et = sqrt(met_px * met_px + met_py * met_py);
   LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met_px << ", " << met_py << " [GeV]";
   fill_histogram("MET_BEFORECUTS", met_et);
  
   // Trigger
   Handle<TriggerResults> triggerResults;
   if (!ev.getByToken(trigToken_, triggerResults)) {
     LogError("") << ">>> TRIGGER collection does not exist !!!";
     return false;
   }
   const edm::TriggerNames& triggerNames = ev.triggerNames(*triggerResults);
   /*
       for (unsigned int i=0; i<triggerResults->size(); i++) {
             if (triggerResults->accept(i)) {
                   LogTrace("") << "Accept by: " << i << ", Trigger: " << triggerNames.triggerName(i);
             }
       }
     */
   bool trigger_fired = false;
   int itrig1 = triggerNames.triggerIndex(muonTrig_);
   if (triggerResults->accept(itrig1))
     trigger_fired = true;
   LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << muonTrig_ << ")";
   fill_histogram("TRIG_BEFORECUTS", trigger_fired);
  
   // 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;
   fill_histogram("NZ1_BEFORECUTS", nmuonsForZ1);
   fill_histogram("NZ2_BEFORECUTS", 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;
   fill_histogram("NJETS_BEFORECUTS", njets);
  
   // Start counting, reject already events if possible (under FastOption flag)
   nall++;
   if (fastOption_ && !trigger_fired)
     return false;
   if (fastOption_ && nmuonsForZ1 >= 1 && nmuonsForZ2 >= 2)
     return false;
   if (fastOption_ && njets > nJetMax_)
     return false;
  
   // Histograms per event shouldbe done only once, so keep track of them
   bool hlt_hist_done = false;
   bool met_hist_done = false;
   bool nz1_hist_done = false;
   bool nz2_hist_done = false;
   bool njets_hist_done = false;
  
   // Central W->mu nu selection criteria
   const int NFLAGS = 13;
   bool muon_sel[NFLAGS];
   for (unsigned int i = 0; i < muonCollectionSize; i++) {
     for (int j = 0; j < NFLAGS; ++j) {
       muon_sel[j] = false;
     }
  
     const Muon& mu = muonCollection->at(i);
     if (!mu.isGlobalMuon())
       continue;
     if (mu.globalTrack().isNull())
       continue;
     if (mu.innerTrack().isNull())
       continue;
  
     LogTrace("") << "> Wsel: processing muon number " << i << "...";
     reco::TrackRef gm = mu.globalTrack();
     //reco::TrackRef tk = mu.innerTrack();
  
     // Pt,eta cuts
     double pt = mu.pt();
     double eta = mu.eta();
     LogTrace("") << "\t... pt, eta: " << pt << " [GeV], " << eta;
     ;
     if (pt > ptCut_)
       muon_sel[0] = true;
     else if (fastOption_)
       continue;
     if (fabs(eta) < etaCut_)
       muon_sel[1] = true;
     else if (fastOption_)
       continue;
  
     // d0, chi2, nhits quality cuts
     double dxy = gm->dxy(beamSpotHandle->position());
     double normalizedChi2 = gm->normalizedChi2();
     double validmuonhits = gm->hitPattern().numberOfValidMuonHits();
     //double standalonehits=mu.outerTrack()->numberOfValidHits();
     double trackerHits = gm->hitPattern().numberOfValidTrackerHits();
     LogTrace("") << "\t... dxy, normalizedChi2, trackerHits, isTrackerMuon?: " << dxy << " [cm], " << normalizedChi2
                  << ", " << trackerHits << ", " << mu.isTrackerMuon();
     if (fabs(dxy) < dxyCut_)
       muon_sel[2] = true;
     else if (fastOption_)
       continue;
     if (muon::isGoodMuon(mu, muon::GlobalMuonPromptTight))
       muon_sel[3] = true;
     else if (fastOption_)
       continue;
     if (trackerHits >= trackerHitsCut_)
       muon_sel[4] = true;
     else if (fastOption_)
       continue;
     if (mu.isTrackerMuon())
       muon_sel[5] = true;
     else if (fastOption_)
       continue;
  
     fill_histogram("PT_BEFORECUTS", pt);
     fill_histogram("ETA_BEFORECUTS", eta);
     fill_histogram("DXY_BEFORECUTS", dxy);
     fill_histogram("CHI2_BEFORECUTS", normalizedChi2);
     fill_histogram("NHITS_BEFORECUTS", trackerHits);
     fill_histogram("ValidMuonHits_BEFORECUTS", validmuonhits);
     fill_histogram("TKMU_BEFORECUTS", mu.isTrackerMuon());
  
     // Isolation cuts
     double isovar = mu.isolationR03().sumPt;
     if (isCombinedIso_) {
       isovar += mu.isolationR03().emEt;
       isovar += mu.isolationR03().hadEt;
     }
     if (isRelativeIso_)
       isovar /= pt;
     if (isovar < isoCut03_)
       muon_sel[6] = true;
     else if (fastOption_)
       continue;
     LogTrace("") << "\t... isolation value" << isovar << ", isolated? " << muon_sel[6];
     fill_histogram("ISO_BEFORECUTS", isovar);
  
     // HLT (not mtched to muon for the time being)
     if (trigger_fired)
       muon_sel[7] = true;
     else if (fastOption_)
       continue;
  
     // MET/MT cuts
     double w_et = met_et + mu.pt();
     double w_px = met_px + mu.px();
     double w_py = met_py + mu.py();
     double massT = w_et * w_et - w_px * w_px - w_py * w_py;
     massT = (massT > 0) ? sqrt(massT) : 0;
  
     LogTrace("") << "\t... W mass, W_et, W_px, W_py: " << massT << ", " << w_et << ", " << w_px << ", " << w_py
                  << " [GeV]";
     if (massT > mtMin_ && massT < mtMax_)
       muon_sel[8] = true;
     else if (fastOption_)
       continue;
     fill_histogram("MT_BEFORECUTS", massT);
     if (met_et > metMin_ && met_et < metMax_)
       muon_sel[9] = true;
     else if (fastOption_)
       continue;
  
     // Acoplanarity cuts
     Geom::Phi<double> deltaphi(mu.phi() - atan2(met_py, met_px));
     double acop = deltaphi.value();
     if (acop < 0)
       acop = -acop;
     acop = M_PI - acop;
     LogTrace("") << "\t... acoplanarity: " << acop;
     if (acop < acopCut_)
       muon_sel[10] = true;
     else if (fastOption_)
       continue;
     fill_histogram("ACOP_BEFORECUTS", acop);
  
     // Remaining flags (from global event information)
     if (nmuonsForZ1 < 1 || nmuonsForZ2 < 2)
       muon_sel[11] = true;
     else if (fastOption_)
       continue;
     if (njets <= nJetMax_)
       muon_sel[12] = true;
     else if (fastOption_)
       continue;
  
     if (fastOption_) {
       all_sel = true;
       break;
     } else {
       // Collect necessary flags "per muon"
       int flags_passed = 0;
       bool rec_sel_this = true;
       bool iso_sel_this = true;
       bool hlt_sel_this = true;
       bool met_sel_this = true;
       bool all_sel_this = true;
       for (int j = 0; j < NFLAGS; ++j) {
         if (muon_sel[j])
           flags_passed += 1;
         if (j < 6 && !muon_sel[j])
           rec_sel_this = false;
         if (j < 7 && !muon_sel[j])
           iso_sel_this = false;
         if (j < 8 && !muon_sel[j])
           hlt_sel_this = false;
         if (j < 11 && !muon_sel[j])
           met_sel_this = false;
         if (!muon_sel[j])
           all_sel_this = false;
       }
  
       // "rec" => pt,eta and quality cuts are satisfied
       if (rec_sel_this)
         rec_sel = true;
       // "iso" => "rec" AND "muon is isolated"
       if (iso_sel_this)
         iso_sel = true;
       // "hlt" => "iso" AND "event is triggered"
       if (hlt_sel_this)
         hlt_sel = true;
       // "met" => "hlt" AND "MET/MT and acoplanarity cuts"
       if (met_sel_this)
         met_sel = true;
       // "all" => "met" AND "Z/top rejection cuts"
       if (all_sel_this)
         all_sel = true;
  
       // Do N-1 histograms now (and only once for global event quantities)
       if (flags_passed >= (NFLAGS - 1)) {
         if (!muon_sel[0] || flags_passed == NFLAGS)
           fill_histogram("PT_LASTCUT", pt);
         if (!muon_sel[1] || flags_passed == NFLAGS)
           fill_histogram("ETA_LASTCUT", eta);
         if (!muon_sel[2] || flags_passed == NFLAGS)
           fill_histogram("DXY_LASTCUT", dxy);
         if (!muon_sel[3] || flags_passed == NFLAGS) {
           fill_histogram("CHI2_LASTCUT", normalizedChi2);
           fill_histogram("ValidMuonHits_LASTCUT", validmuonhits);
         }
         if (!muon_sel[4] || flags_passed == NFLAGS)
           fill_histogram("NHITS_LASTCUT", trackerHits);
         if (!muon_sel[5] || flags_passed == NFLAGS)
           fill_histogram("TKMU_LASTCUT", mu.isTrackerMuon());
         if (!muon_sel[6] || flags_passed == NFLAGS)
           fill_histogram("ISO_LASTCUT", isovar);
         if (!muon_sel[7] || flags_passed == NFLAGS) {
           if (!hlt_hist_done) {
             fill_histogram("TRIG_LASTCUT", trigger_fired);
             hlt_hist_done = true;
           }
         }
         if (!muon_sel[8] || flags_passed == NFLAGS)
           fill_histogram("MT_LASTCUT", massT);
         if (!muon_sel[9] || flags_passed == NFLAGS) {
           if (!met_hist_done) {
             fill_histogram("MET_LASTCUT", met_et);
             met_hist_done = true;
           }
         }
         if (!muon_sel[10] || flags_passed == NFLAGS)
           fill_histogram("ACOP_LASTCUT", acop);
         if (!muon_sel[11] || flags_passed == NFLAGS) {
           if (!nz1_hist_done) {
             fill_histogram("NZ1_LASTCUT", nmuonsForZ1);
             nz1_hist_done = true;
           }
         }
         if (!muon_sel[11] || flags_passed == NFLAGS) {
           if (!nz2_hist_done) {
             fill_histogram("NZ2_LASTCUT", nmuonsForZ2);
             nz2_hist_done = true;
           }
         }
         if (!muon_sel[12] || flags_passed == NFLAGS) {
           if (!njets_hist_done) {
             fill_histogram("NJETS_LASTCUT", njets);
             njets_hist_done = true;
           }
         }
       }
     }
   }
  
   // Collect final flags
   if (!fastOption_) {
     if (rec_sel)
       nrec++;
     if (iso_sel)
       niso++;
     if (hlt_sel)
       nhlt++;
     if (met_sel)
       nmet++;
   }
  
   if (all_sel) {
     nsel++;
     LogTrace("") << ">>>> Event ACCEPTED";
   } else {
     LogTrace("") << ">>>> Event REJECTED";
   }
  
   return all_sel;
 }

References acopCut_, beamSpotToken_, PVValHelper::dxy, dxyCut_, eJetMin_, PVValHelper::eta, etaCut_, ev, fastOption_, fill_histogram(), muon::GlobalMuonPromptTight, mps_fire::i, isCombinedIso_, muon::isGoodMuon(), isoCut03_, isRelativeIso_, dqmiolumiharvest::j, metsig::jet, jetfilter_cfi::jetCollection, jetToken_, LogTrace, M_PI, BTaggingMonitor_cfi::met, susyDQM_cfi::metCollection, metIncludesMuons_, metMax_, metMin_, metToken_, mtMax_, mtMin_, amptDefaultParameters_cff::mu, pdwgLeptonRecoSkim_cfi::muonCollection, muonToken_, muonTrig_, nall, nhlt, niso, nJetMax_, BTaggingMonitoring_cff::njets, nmet, TrackCollections2monitor_cff::normalizedChi2, nrec, nsel, reco::BeamSpot::position(), DiDispStaMuonMonitor_cfi::pt, ptCut_, ptThrForZ1_, ptThrForZ2_, mathSSE::sqrt(), trackerHitsCut_, L1TEGammaOffline_cfi::triggerNames, triggerResults, trigToken_, and Geom::Phi< T1, Range >::value().

◆ init_histograms()

void WMuNuValidator::init_histograms ( )

Definition at line 152 of file WMuNuValidator.cc.

                                      {
   edm::Service<TFileService> fs;
   TFileDirectory subDir0 = fs->mkdir("BeforeCuts");
   TFileDirectory subDir1 = fs->mkdir("LastCut");
   TFileDirectory* subDir[2];
   subDir[0] = &subDir0;
   subDir[1] = &subDir1;
  
   char chname[256] = "";
   char chtitle[256] = "";
   std::string chsuffix[2] = {"_BEFORECUTS", "_LASTCUT"};
  
   for (int i = 0; i < 2; ++i) {
     snprintf(chname, 255, "PT%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Muon transverse momentum [GeV]");
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 100, 0., 100.);
  
     snprintf(chname, 255, "ETA%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Muon pseudo-rapidity");
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 50, -2.5, 2.5);
  
     snprintf(chname, 255, "DXY%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Muon transverse distance to beam spot [cm]");
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 100, -0.5, 0.5);
  
     snprintf(chname, 255, "CHI2%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Normalized Chi2, inner track fit");
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 100, 0., 100.);
  
     snprintf(chname, 255, "NHITS%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Number of hits, inner track");
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 40, -0.5, 39.5);
  
     snprintf(chname, 255, "ValidMuonHits%s", chsuffix[i].data());
     snprintf(chtitle, 255, "number Of Valid Muon Hits");
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 40, -0.5, 39.5);
  
     snprintf(chname, 255, "TKMU%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Tracker-muon flag (for global muons)");
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 2, -0.5, 1.5);
  
     snprintf(chname, 255, "ISO%s", chsuffix[i].data());
     if (isRelativeIso_) {
       if (isCombinedIso_) {
         snprintf(chtitle, 255, "Relative (combined) isolation variable");
       } else {
         snprintf(chtitle, 255, "Relative (tracker) isolation variable");
       }
       h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 100, 0., 1.);
     } else {
       if (isCombinedIso_) {
         snprintf(chtitle, 255, "Absolute (combined) isolation variable [GeV]");
       } else {
         snprintf(chtitle, 255, "Absolute (tracker) isolation variable [GeV]");
       }
       h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 100, 0., 20.);
     }
  
     snprintf(chname, 255, "TRIG%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Trigger response (bit %s)", muonTrig_.data());
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 2, -0.5, 1.5);
  
     snprintf(chname, 255, "MT%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 150, 0., 300.);
  
     snprintf(chname, 255, "MET%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Missing transverse energy (%s) [GeV]", metTag_.label().data());
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 100, 0., 200.);
  
     snprintf(chname, 255, "ACOP%s", chsuffix[i].data());
     snprintf(chtitle, 255, "MU-MET (%s) acoplanarity", metTag_.label().data());
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 50, 0., M_PI);
  
     snprintf(chname, 255, "NZ1%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ1_);
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 10, -0.5, 9.5);
  
     snprintf(chname, 255, "NZ2%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ2_);
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 10, -0.5, 9.5);
  
     snprintf(chname, 255, "NJETS%s", chsuffix[i].data());
     snprintf(chtitle, 255, "Number of jets (%s) above %.2f GeV", jetTag_.label().data(), eJetMin_);
     h1_[chname] = subDir[i]->make<TH1D>(chname, chtitle, 10, -0.5, 9.5);
   }
 }

References data, eJetMin_, h1_, mps_fire::i, isCombinedIso_, isRelativeIso_, jetTag_, edm::InputTag::label(), M_PI, metTag_, TFileService::mkdir(), muonTrig_, ptThrForZ1_, ptThrForZ2_, AlCaHLTBitMon_QueryRunRegistry::string, and PostProcessorHcalIsoTrack_cfi::subDir.

Referenced by beginJob().

Member Data Documentation

◆ acopCut_

double WMuNuValidator::acopCut_

private

Definition at line 54 of file WMuNuValidator.cc.

Referenced by filter().

◆ beamSpotToken_

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

private

Definition at line 43 of file WMuNuValidator.cc.

Referenced by filter().

◆ dxyCut_

double WMuNuValidator::dxyCut_

private

Definition at line 56 of file WMuNuValidator.cc.

Referenced by filter().

◆ eJetMin_

double WMuNuValidator::eJetMin_

private

Definition at line 64 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

◆ etaCut_

double WMuNuValidator::etaCut_

private

Definition at line 46 of file WMuNuValidator.cc.

Referenced by filter().

◆ fastOption_

bool WMuNuValidator::fastOption_

private

Definition at line 35 of file WMuNuValidator.cc.

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

◆ h1_

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

private

Definition at line 74 of file WMuNuValidator.cc.

Referenced by fill_histogram(), and init_histograms().

◆ isAlsoTrackerMuon_

bool WMuNuValidator::isAlsoTrackerMuon_

private

Definition at line 59 of file WMuNuValidator.cc.

◆ isCombinedIso_

bool WMuNuValidator::isCombinedIso_

private

Definition at line 48 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

◆ isoCut03_

double WMuNuValidator::isoCut03_

private

Definition at line 49 of file WMuNuValidator.cc.

Referenced by filter().

◆ isRelativeIso_

bool WMuNuValidator::isRelativeIso_

private

Definition at line 47 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

◆ jetTag_

edm::InputTag WMuNuValidator::jetTag_

private

Definition at line 41 of file WMuNuValidator.cc.

Referenced by init_histograms().

◆ jetToken_

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

private

Definition at line 42 of file WMuNuValidator.cc.

Referenced by filter().

◆ metIncludesMuons_

bool WMuNuValidator::metIncludesMuons_

private

Definition at line 40 of file WMuNuValidator.cc.

Referenced by filter().

◆ metMax_

double WMuNuValidator::metMax_

private

Definition at line 53 of file WMuNuValidator.cc.

Referenced by filter().

◆ metMin_

double WMuNuValidator::metMin_

private

Definition at line 52 of file WMuNuValidator.cc.

Referenced by filter().

◆ metTag_

edm::InputTag WMuNuValidator::metTag_

private

Definition at line 38 of file WMuNuValidator.cc.

Referenced by init_histograms().

◆ metToken_

edm::EDGetTokenT<edm::View<reco::MET> > WMuNuValidator::metToken_

private

Definition at line 39 of file WMuNuValidator.cc.

Referenced by filter().

◆ mtMax_

double WMuNuValidator::mtMax_

private

Definition at line 51 of file WMuNuValidator.cc.

Referenced by filter().

◆ mtMin_

double WMuNuValidator::mtMin_

private

Definition at line 50 of file WMuNuValidator.cc.

Referenced by filter().

◆ muonToken_

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

private

Definition at line 37 of file WMuNuValidator.cc.

Referenced by filter().

◆ muonTrig_

const std::string WMuNuValidator::muonTrig_

private

Definition at line 44 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

◆ nall

unsigned int WMuNuValidator::nall

private

Definition at line 67 of file WMuNuValidator.cc.

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

◆ nhlt

unsigned int WMuNuValidator::nhlt

private

Definition at line 70 of file WMuNuValidator.cc.

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

◆ niso

unsigned int WMuNuValidator::niso

private

Definition at line 69 of file WMuNuValidator.cc.

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

◆ nJetMax_

int WMuNuValidator::nJetMax_

private

Definition at line 65 of file WMuNuValidator.cc.

Referenced by filter().

◆ nmet

unsigned int WMuNuValidator::nmet

private

Definition at line 71 of file WMuNuValidator.cc.

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

◆ normalizedChi2Cut_

double WMuNuValidator::normalizedChi2Cut_

private

Definition at line 57 of file WMuNuValidator.cc.

◆ nrec

unsigned int WMuNuValidator::nrec

private

Definition at line 68 of file WMuNuValidator.cc.

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

◆ nsel

unsigned int WMuNuValidator::nsel

private

Definition at line 72 of file WMuNuValidator.cc.

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

◆ ptCut_

double WMuNuValidator::ptCut_

private

Definition at line 45 of file WMuNuValidator.cc.

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

◆ ptThrForZ1_

double WMuNuValidator::ptThrForZ1_

private

Definition at line 61 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

◆ ptThrForZ2_

double WMuNuValidator::ptThrForZ2_

private

Definition at line 62 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

◆ trackerHitsCut_

int WMuNuValidator::trackerHitsCut_

private

Definition at line 58 of file WMuNuValidator.cc.

Referenced by filter().

◆ trigToken_

edm::EDGetTokenT<edm::TriggerResults> WMuNuValidator::trigToken_

private

Definition at line 36 of file WMuNuValidator.cc.

Referenced by filter().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ WMuNuValidator()

Member Function Documentation

◆ beginJob()

◆ endJob()

◆ fill_histogram()

◆ filter()

◆ init_histograms()

Member Data Documentation

◆ acopCut_

◆ beamSpotToken_

◆ dxyCut_

◆ eJetMin_

◆ etaCut_

◆ fastOption_

◆ h1_

◆ isAlsoTrackerMuon_

◆ isCombinedIso_

◆ isoCut03_

◆ isRelativeIso_

◆ jetTag_

◆ jetToken_

◆ metIncludesMuons_

◆ metMax_

◆ metMin_

◆ metTag_

◆ metToken_

◆ mtMax_

◆ mtMin_

◆ muonToken_

◆ muonTrig_

◆ nall

◆ nhlt

◆ niso

◆ nJetMax_

◆ nmet

◆ normalizedChi2Cut_

◆ nrec

◆ nsel

◆ ptCut_

◆ ptThrForZ1_

◆ ptThrForZ2_

◆ trackerHitsCut_

◆ trigToken_