Inheritance diagram for WMuNuValidator:

Public Member Functions
virtual void	beginJob () override

virtual void	endJob () override

void	fill_histogram (const char *, const double &)

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

void	init_histograms ()

	WMuNuValidator (const edm::ParameterSet &)

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

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

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

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 &	itemsToGetFromEvent () 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

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

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

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 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 26 of file WMuNuValidator.cc.

Constructor & Destructor Documentation

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

void WMuNuValidator::beginJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 140 of file WMuNuValidator.cc.

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

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

void WMuNuValidator::endJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 245 of file WMuNuValidator.cc.

References Vispa.Plugins.EdmBrowser.EdmDataAccessor::all(), fastOption_, nall, nhlt, niso, nmet, nrec, nsel, pileupDistInMC::num, and mathSSE::sqrt().

                             {
       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";
 }

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

Definition at line 240 of file WMuNuValidator.cc.

References fastOption_, h1_, and mergeVDriftHistosByStation::name.

Referenced by filter().

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

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

overridevirtual

Implements edm::EDFilter.

Definition at line 304 of file WMuNuValidator.cc.

References acopCut_, beamSpotToken_, dxyCut_, eJetMin_, reco::MuonIsolation::emEt, reco::LeafCandidate::et(), eta, reco::LeafCandidate::eta(), etaCut_, fastOption_, fill_histogram(), edm::Event::getByToken(), muon::GlobalMuonPromptTight, reco::Muon::globalTrack(), reco::MuonIsolation::hadEt, i, reco::Muon::innerTrack(), isCombinedIso_, reco::Muon::isGlobalMuon(), muon::isGoodMuon(), edm::Ref< C, T, F >::isNull(), isoCut03_, reco::Muon::isolationR03(), isRelativeIso_, reco::Muon::isTrackerMuon(), j, metsig::jet, jetToken_, LogTrace, M_PI, objects.METAnalyzer::met, MT2Analyzer::metCollection, metIncludesMuons_, metMax_, metMin_, metToken_, mtMax_, mtMin_, RPCpg::mu, muonToken_, muonTrig_, nall, nhlt, niso, nJetMax_, HLT_25ns10e33_v2_cff::njets, nmet, nrec, nsel, reco::LeafCandidate::phi(), EnergyCorrector::pt, reco::LeafCandidate::pt(), ptCut_, ptThrForZ1_, ptThrForZ2_, reco::LeafCandidate::px(), reco::LeafCandidate::py(), mathSSE::sqrt(), reco::MuonIsolation::sumPt, trackerHitsCut_, edm::TriggerNames::triggerIndex(), edm::Event::triggerNames(), triggerResults, trigToken_, and Geom::Phi< T >::value().

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

                                                            {
 
       // 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 0;
       }
       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 0;
       }
       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 0;
       }
       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;
 
 }

void WMuNuValidator::init_histograms ( )

Definition at line 153 of file WMuNuValidator.cc.

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

Referenced by beginJob().

                                      {
       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);
 
       }
 }

Member Data Documentation

double WMuNuValidator::acopCut_

private

Definition at line 54 of file WMuNuValidator.cc.

Referenced by filter().

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

private

Definition at line 43 of file WMuNuValidator.cc.

Referenced by filter().

double WMuNuValidator::dxyCut_

private

Definition at line 56 of file WMuNuValidator.cc.

Referenced by filter().

double WMuNuValidator::eJetMin_

private

Definition at line 64 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

double WMuNuValidator::etaCut_

private

Definition at line 46 of file WMuNuValidator.cc.

Referenced by filter().

bool WMuNuValidator::fastOption_

private

Definition at line 35 of file WMuNuValidator.cc.

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

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

private

Definition at line 74 of file WMuNuValidator.cc.

Referenced by fill_histogram(), and init_histograms().

bool WMuNuValidator::isAlsoTrackerMuon_

private

Definition at line 59 of file WMuNuValidator.cc.

bool WMuNuValidator::isCombinedIso_

private

Definition at line 48 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

double WMuNuValidator::isoCut03_

private

Definition at line 49 of file WMuNuValidator.cc.

Referenced by filter().

bool WMuNuValidator::isRelativeIso_

private

Definition at line 47 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

edm::InputTag WMuNuValidator::jetTag_

private

Definition at line 41 of file WMuNuValidator.cc.

Referenced by init_histograms().

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

private

Definition at line 42 of file WMuNuValidator.cc.

Referenced by filter().

bool WMuNuValidator::metIncludesMuons_

private

Definition at line 40 of file WMuNuValidator.cc.

Referenced by filter().

double WMuNuValidator::metMax_

private

Definition at line 53 of file WMuNuValidator.cc.

Referenced by filter().

double WMuNuValidator::metMin_

private

Definition at line 52 of file WMuNuValidator.cc.

Referenced by filter().

edm::InputTag WMuNuValidator::metTag_

private

Definition at line 38 of file WMuNuValidator.cc.

Referenced by init_histograms().

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

private

Definition at line 39 of file WMuNuValidator.cc.

Referenced by filter().

double WMuNuValidator::mtMax_

private

Definition at line 51 of file WMuNuValidator.cc.

Referenced by filter().

double WMuNuValidator::mtMin_

private

Definition at line 50 of file WMuNuValidator.cc.

Referenced by filter().

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

private

Definition at line 37 of file WMuNuValidator.cc.

Referenced by filter().

const std::string WMuNuValidator::muonTrig_

private

Definition at line 44 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

unsigned int WMuNuValidator::nall

private

Definition at line 67 of file WMuNuValidator.cc.

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

unsigned int WMuNuValidator::nhlt

private

Definition at line 70 of file WMuNuValidator.cc.

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

unsigned int WMuNuValidator::niso

private

Definition at line 69 of file WMuNuValidator.cc.

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

int WMuNuValidator::nJetMax_

private

Definition at line 65 of file WMuNuValidator.cc.

Referenced by filter().

unsigned int WMuNuValidator::nmet

private

Definition at line 71 of file WMuNuValidator.cc.

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

double WMuNuValidator::normalizedChi2Cut_

private

Definition at line 57 of file WMuNuValidator.cc.

unsigned int WMuNuValidator::nrec

private

Definition at line 68 of file WMuNuValidator.cc.

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

unsigned int WMuNuValidator::nsel

private

Definition at line 72 of file WMuNuValidator.cc.

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

double WMuNuValidator::ptCut_

private

Definition at line 45 of file WMuNuValidator.cc.

Referenced by filter().

double WMuNuValidator::ptThrForZ1_

private

Definition at line 61 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

double WMuNuValidator::ptThrForZ2_

private

Definition at line 62 of file WMuNuValidator.cc.

Referenced by filter(), and init_histograms().

int WMuNuValidator::trackerHitsCut_

private

Definition at line 58 of file WMuNuValidator.cc.

Referenced by filter().

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

Member Function Documentation

Member Data Documentation