WMuNuValidator Class Reference

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
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
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)
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	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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 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 dataset::name.

Referenced by filter().

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

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

override

Definition at line 304 of file WMuNuValidator.cc.

References edm::HLTGlobalStatus::accept(), acopCut_, beamSpotToken_, DEFINE_FWK_MODULE, PVValHelper::dxy, dxyCut_, eJetMin_, reco::MuonIsolation::emEt, reco::LeafCandidate::et(), PVValHelper::eta, reco::LeafCandidate::eta(), etaCut_, fastOption_, fill_histogram(), edm::Event::getByToken(), muon::GlobalMuonPromptTight, reco::Muon::globalTrack(), reco::MuonIsolation::hadEt, mps_fire::i, reco::Muon::innerTrack(), isCombinedIso_, reco::Muon::isGlobalMuon(), muon::isGoodMuon(), edm::Ref< C, T, F >::isNull(), isoCut03_, reco::Muon::isolationR03(), isRelativeIso_, reco::Muon::isTrackerMuon(), metsig::jet, jetToken_, LogTrace, M_PI, RazorAnalyzer::met, MT2Analyzer::metCollection, metIncludesMuons_, metMax_, metMin_, metToken_, mtMax_, mtMin_, RPCpg::mu, HiRecoMuon_cff::muonCollection, muonToken_, muonTrig_, nall, nhlt, niso, nJetMax_, nmet, nrec, nsel, reco::LeafCandidate::phi(), reco::BeamSpot::position(), 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 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;

}

void WMuNuValidator::init_histograms

(

)

Definition at line 153 of file WMuNuValidator.cc.

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

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 badGlobalMuonTaggerFWLite.BadGlobalMuonTagger::badMuons(), and 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().