EwkMuLumiMonitorDQM Class Reference

#include <EwkMuLumiMonitorDQM.h>

Inheritance diagram for EwkMuLumiMonitorDQM:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
virtual void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
virtual void	endRun (const edm::Run &, const edm::EventSetup &) override
	EwkMuLumiMonitorDQM (const edm::ParameterSet &)
void	init_histograms ()
bool	IsMuMatchedToHLTMu (const reco::Muon &, const std::vector< reco::Particle > &, double, double)
double	muIso (const reco::Muon &)
double	tkIso (const reco::Track &, edm::Handle< reco::TrackCollection >, edm::Handle< CaloTowerCollection >)
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
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
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
double	acopCut_
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
edm::EDGetTokenT< CaloTowerCollection >	caloTowerToken_
double	deltaRTrk_
double	deltaRVetoTrk_
double	dxyCut_
double	etaMuCut_
MonitorElement *	highMass1HLT_
MonitorElement *	highMass2HLT_
MonitorElement *	highMassGlbSta_
MonitorElement *	highMassGlbTrk_
MonitorElement *	highMassIsBothGlbTrkThanW_
MonitorElement *	highMassNotIso_
HLTConfigProvider	hltConfigProvider_
bool	isCombinedIso_
double	isoCut03_
bool	isRelativeIso_
bool	isValidHltConfig_
bool	isW_
bool	isZGlbSta_
bool	isZGlbTrk_
bool	isZGolden1HLT_
bool	isZGolden2HLT_
bool	isZGoldenNoIso_
MonitorElement *	mass1HLT_
MonitorElement *	mass2HLT_
MonitorElement *	massGlbSta_
MonitorElement *	massGlbTrk_
MonitorElement *	massIsBothGlbTrkThanW_
MonitorElement *	massNotIso_
double	maxDeltaR_
double	maxDPtRel_
bool	metIncludesMuons_
edm::EDGetTokenT< edm::View< reco::MET > >	metToken_
double	mtMax_
double	mtMin_
edm::EDGetTokenT< edm::View< reco::Muon > >	muonToken_
unsigned int	n1hlt
unsigned int	n2hlt
unsigned int	nall
unsigned int	nEvWithHighPtMu
unsigned int	nGlbSta
unsigned int	nGlbTrk
unsigned int	nhlt
unsigned int	nInKinRange
unsigned int	niso
unsigned int	nNotIso
unsigned int	nsel
unsigned int	nTMass
unsigned int	nW
double	ptMuCut_
double	ptThreshold_
MonitorElement *	TMass_
edm::EDGetTokenT< reco::TrackCollection >	trackToken_
edm::EDGetTokenT< trigger::TriggerEvent >	trigEvToken_
edm::InputTag	trigTag_
edm::EDGetTokenT< edm::TriggerResults >	trigToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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

DQM offline for EWK MuLumiMonitor: intended for luminosity purposes using Z/W

Authors

: Michele de Gruttola, INFN Naples - Maria Cepeda, CIEMAT on behalf EWK-Muon group

Definition at line 37 of file EwkMuLumiMonitorDQM.h.

Constructor & Destructor Documentation

EwkMuLumiMonitorDQM::EwkMuLumiMonitorDQM

(

const edm::ParameterSet &

cfg

)

Definition at line 36 of file EwkMuLumiMonitorDQM.cc.

References isValidHltConfig_.

    :
      // Input collections
      trigTag_(cfg.getUntrackedParameter<edm::InputTag>(
          "TrigTag", edm::InputTag("TriggerResults::HLT"))),
      trigToken_(consumes<edm::TriggerResults>(trigTag_)),
      trigEvToken_(consumes<trigger::TriggerEvent>(
          cfg.getUntrackedParameter<edm::InputTag>("triggerEvent"))),
      beamSpotToken_(
          consumes<reco::BeamSpot>(cfg.getUntrackedParameter<edm::InputTag>(
              "offlineBeamSpot", edm::InputTag("offlineBeamSpot")))),
      muonToken_(consumes<edm::View<reco::Muon> >(
          cfg.getUntrackedParameter<edm::InputTag>("muons"))),
      trackToken_(consumes<reco::TrackCollection>(
          cfg.getUntrackedParameter<edm::InputTag>("tracks"))),
      caloTowerToken_(consumes<CaloTowerCollection>(
          cfg.getUntrackedParameter<edm::InputTag>("calotower"))),
      metToken_(consumes<edm::View<reco::MET> >(
          cfg.getUntrackedParameter<edm::InputTag>("metTag"))),
      metIncludesMuons_(cfg.getUntrackedParameter<bool>("METIncludesMuons")),
      // Main cuts
      // massMin_(cfg.getUntrackedParameter<double>("MtMin", 20.)),
      //   massMax_(cfg.getUntrackedParameter<double>("MtMax", 2000.))
      //  hltPath_(cfg.getUntrackedParameter<std::string> ("hltPath")) ,
      //  L3FilterName_(cfg.getUntrackedParameter<std::string>
      // ("L3FilterName")),
      ptMuCut_(cfg.getUntrackedParameter<double>("ptMuCut")),
      etaMuCut_(cfg.getUntrackedParameter<double>("etaMuCut")),
      isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso")),
      isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso")),
      isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03")),
      //  deltaRTrk_(cfg.getUntrackedParameter<double>("deltaRTrk")),
      ptThreshold_(cfg.getUntrackedParameter<double>("ptThreshold")),
      // deltaRVetoTrk_(cfg.getUntrackedParameter<double>("deltaRVetoTrk")),
      maxDPtRel_(cfg.getUntrackedParameter<double>("maxDPtRel")),
      maxDeltaR_(cfg.getUntrackedParameter<double>("maxDeltaR")),
      mtMin_(cfg.getUntrackedParameter<double>("mtMin")),
      mtMax_(cfg.getUntrackedParameter<double>("mtMax")),
      acopCut_(cfg.getUntrackedParameter<double>("acopCut")),
      dxyCut_(cfg.getUntrackedParameter<double>("DxyCut")) {
  // just to initialize
  isValidHltConfig_ = false;

}

Member Function Documentation

void EwkMuLumiMonitorDQM::analyze ( const edm::Event &  ev, const edm::EventSetup &    )

overridevirtual

number of prescale sets available

Definition at line 216 of file EwkMuLumiMonitorDQM.cc.

References edm::HLTGlobalStatus::accept(), acopCut_, beamSpotToken_, caloTowerToken_, reco::LeafCandidate::charge(), reco::TrackBase::charge(), PVValHelper::dxy, reco::TrackBase::dxy(), dxyCut_, edm::InputTag::encode(), PVValHelper::eta, reco::LeafCandidate::eta(), reco::TrackBase::eta(), etaMuCut_, f, MonitorElement::Fill(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), highMass1HLT_, highMass2HLT_, highMassGlbSta_, highMassGlbTrk_, highMassIsBothGlbTrkThanW_, highMassNotIso_, hltConfigProvider_, mps_fire::i, reco::Muon::innerTrack(), reco::Muon::isGlobalMuon(), IsMuMatchedToHLTMu(), isoCut03_, reco::Muon::isStandAloneMuon(), isW_, isZGlbSta_, isZGlbTrk_, isZGolden1HLT_, isZGolden2HLT_, isZGoldenNoIso_, gen::k, LogTrace, M_PI, ResonanceBuilder::mass, mass1HLT_, mass2HLT_, massGlbSta_, massGlbTrk_, massIsBothGlbTrkThanW_, massNotIso_, maxDeltaR_, maxDPtRel_, RazorAnalyzer::met, MT2Analyzer::metCollection, metIncludesMuons_, metToken_, HLTConfigProvider::moduleLabels(), mtMax_, RPCpg::mu, muIso(), electronCleaner_cfi::muons, muonToken_, gen::n, n1hlt, n2hlt, nall, dataset::name, nEvWithHighPtMu, nGlbSta, nGlbTrk, nhlt, niso, nNotIso, nsel, nTMass, nW, MillePedeFileConverter_cfg::out, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), reco::BeamSpot::position(), HLTConfigProvider::prescaleSize(), HLTConfigProvider::prescaleValue(), EnergyCorrector::pt, reco::LeafCandidate::pt(), reco::TrackBase::pt(), hiDetachedQuadStep_cff::pt2, ptMuCut_, reco::LeafCandidate::px(), reco::TrackBase::px(), reco::LeafCandidate::py(), reco::TrackBase::py(), reco::TrackBase::pz(), alignCSCRings::s, findQualityFiles::size, trigger::TriggerEvent::sizeFilters(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, tkIso(), TMass_, l1t::tracks, trackToken_, trigEvToken_, HLTConfigProvider::triggerIndex(), HLTConfigProvider::triggerNames(), triggerResults, trigToken_, and findQualityFiles::v.

                                                                    {
  nall++;
  bool hlt_sel = false;
  double iso1 = -1;
  double iso2 = -1;
  bool isMu1Iso = false;
  bool isMu2Iso = false;
  bool singleTrigFlag1 = false;
  bool singleTrigFlag2 = false;
  isZGolden1HLT_ = false;
  isZGolden2HLT_ = false;
  isZGoldenNoIso_ = false;
  isZGlbSta_ = false;
  isZGlbTrk_ = false;
  isW_ = false;
  // Trigger
  bool trigger_fired = false;

  Handle<TriggerResults> triggerResults;
  if (!ev.getByToken(trigToken_, triggerResults)) {
    // LogWarning("") << ">>> TRIGGER collection does not exist !!!";
    return;
  }

  ev.getByToken(trigToken_, triggerResults);
  /*
    const edm::TriggerNames & trigNames = ev.triggerNames(*triggerResults);


  for (size_t i=0; i<triggerResults->size(); i++) {
  std::string trigName = trigNames.triggerName(i);
  //std::cout << " trigName == " << trigName << std::endl;
    if ( trigName == hltPath_ && triggerResults->accept(i)) {
    trigger_fired = true;
    hlt_sel=true;
    nhlt++;
    }
    }
  */
  // see the trigger single muon which are present
  string lowestMuonUnprescaledTrig = "";
  bool lowestMuonUnprescaledTrigFound = false;
  const std::vector<std::string>& triggerNames =
      hltConfigProvider_.triggerNames();
  for (size_t ts = 0; ts < triggerNames.size(); ts++) {
    string trig = triggerNames[ts];
    size_t f = trig.find("HLT_Mu");
    if ((f != std::string::npos)) {
      // std::cout << "single muon trigger present: " << trig << std::endl;
      // See if the trigger is prescaled;
      bool prescaled = false;
      const unsigned int prescaleSize = hltConfigProvider_.prescaleSize();
      for (unsigned int ps = 0; ps < prescaleSize; ps++) {
        const unsigned int prescaleValue =
            hltConfigProvider_.prescaleValue(ps, trig);
        if (prescaleValue != 1) prescaled = true;
        //  std::cout<< " prescaleValue[" << ps << "] =" << prescaleValue
        //<<std::endl;
      }
      if (!prescaled) {
        // looking now for the lowest hlt path not prescaled, with name of the
        // form HLT_MuX or HLTMuX_vY
        for (unsigned int n = 9; n < 100; n++) {
          string lowestTrig = "HLT_Mu";
          string lowestTrigv0 = "copy";
          std::stringstream out;
          out << n;
          std::string s = out.str();
          lowestTrig.append(s);
          lowestTrigv0 = lowestTrig;
          for (unsigned int v = 1; v < 10; v++) {
            lowestTrig.append("_v");
            std::stringstream oout;
            oout << v;
            std::string ss = oout.str();
            lowestTrig.append(ss);
            if (trig == lowestTrig) lowestMuonUnprescaledTrig = trig;
            if (trig == lowestTrig) lowestMuonUnprescaledTrigFound = true;
            if (trig == lowestTrig) break;
          }
          if (lowestMuonUnprescaledTrigFound) break;

          lowestTrig = lowestTrigv0;
          if (trig == lowestTrig) lowestMuonUnprescaledTrig = trig;
          //      if (trig==lowestTrig) {std::cout << " before break, lowestTrig
          // lowest single muon trigger present unprescaled: " << lowestTrig <<
          // std::endl; }
          if (trig == lowestTrig) lowestMuonUnprescaledTrigFound = true;
          if (trig == lowestTrig) break;
        }
        if (lowestMuonUnprescaledTrigFound) break;
      }
    }
  }
  //  std::cout << "after break, lowest single muon trigger present unprescaled:
  // " << lowestMuonUnprescaledTrig << std::endl;
  unsigned int triggerIndex;  // index of trigger path

  // See if event passed trigger paths
  std::string hltPath_ = lowestMuonUnprescaledTrig;

  triggerIndex = hltConfigProvider_.triggerIndex(hltPath_);
  if (triggerIndex < triggerResults->size())
    trigger_fired = triggerResults->accept(triggerIndex);
  std::string L3FilterName_ = "";
  if (trigger_fired) {
    const std::vector<std::string>& moduleLabs =
        hltConfigProvider_.moduleLabels(hltPath_);
    /*for (size_t k =0; k < moduleLabs.size()-1 ; k++){
      std::cout << "moduleLabs[" << k << "] == " << moduleLabs[k] << std::endl;
    }
    */
    // the l3 filter name is just the last module....
    size_t moduleLabsSizeMinus2 = moduleLabs.size() - 2;
    //  std::cout<<"moduleLabs[" << moduleLabsSizeMinus2 << "]== "<<
    // moduleLabs[moduleLabsSizeMinus2] << std::endl;

    L3FilterName_ = moduleLabs[moduleLabsSizeMinus2];
  }

  edm::Handle<trigger::TriggerEvent> handleTriggerEvent;
  LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << hltPath_
               << ")";
  if (!ev.getByToken(trigEvToken_, handleTriggerEvent)) {
    // LogWarning( "errorTriggerEventValid" ) << "trigger::TriggerEvent product
    // with InputTag " << trigEv_.encode() << " not in event";
    return;
  }
  ev.getByToken(trigEvToken_, handleTriggerEvent);
  const trigger::TriggerObjectCollection& toc(handleTriggerEvent->getObjects());
  size_t nMuHLT = 0;
  std::vector<reco::Particle> HLTMuMatched;
  for (size_t ia = 0; ia < handleTriggerEvent->sizeFilters(); ++ia) {
    std::string fullname = handleTriggerEvent->filterTag(ia).encode();
    // std::cout<< "fullname::== " << fullname<< std::endl;
    std::string name;
    size_t p = fullname.find_first_of(':');
    if (p != std::string::npos) {
      name = fullname.substr(0, p);
    } else {
      name = fullname;
    }
    if (!toc.empty()) {
      const trigger::Keys& k = handleTriggerEvent->filterKeys(ia);
      for (trigger::Keys::const_iterator ki = k.begin(); ki != k.end(); ++ki) {
        // looking at all the single muon l3 trigger present, for example
        // hltSingleMu15L3Filtered15.....
        if (name == L3FilterName_) {
          // trigger_fired = true;
          hlt_sel = true;
          nhlt++;
          HLTMuMatched.push_back(toc[*ki].particle());
          nMuHLT++;
        }
      }
    }
  }

  // Beam spot
  Handle<reco::BeamSpot> beamSpotHandle;
  if (!ev.getByToken(beamSpotToken_, beamSpotHandle)) {
    // LogWarning("") << ">>> No beam spot found !!!";
    return;
  }

  //  looping on muon....
  Handle<View<Muon> > muons;
  if (!ev.getByToken(muonToken_, muons)) {
    // LogError("") << ">>> muon collection does not exist !!!";
    return;
  }

  ev.getByToken(muonToken_, muons);
  // saving only muons with pt> ptMuCut and eta<etaMuCut, and dxy<dxyCut
  std::vector<reco::Muon> highPtGlbMuons;
  std::vector<reco::Muon> highPtStaMuons;

  for (size_t i = 0; i < muons->size(); i++) {
    const reco::Muon& mu = muons->at(i);
    double pt = mu.pt();
    double eta = mu.eta();
    if (pt > ptMuCut_ && fabs(eta) < etaMuCut_) {
      if (mu.isGlobalMuon()) {
        // check the dxy....
        double dxy = mu.innerTrack()->dxy(beamSpotHandle->position());
        if (fabs(dxy) > dxyCut_) continue;
        highPtGlbMuons.push_back(mu);
      }
      if (mu.isGlobalMuon()) continue;
      // if is not, look if it is a standalone....
      if (mu.isStandAloneMuon()) highPtStaMuons.push_back(mu);
      nEvWithHighPtMu++;
    }
  }
  size_t nHighPtGlbMu = highPtGlbMuons.size();
  size_t nHighPtStaMu = highPtStaMuons.size();
  if (hlt_sel && (nHighPtGlbMu > 0)) {
    // loop on high pt muons if there's at least two with opposite charge build
    // a Z, more then one z candidate is foreseen.........
    // stop the loop after 10 cicles....
    (nHighPtGlbMu > 10) ? nHighPtGlbMu = 10 : 1;
    // Z selection
    if (nHighPtGlbMu > 1) {
      for (unsigned int i = 0; i < nHighPtGlbMu; i++) {
        reco::Muon muon1 = highPtGlbMuons[i];
        math::XYZTLorentzVector mu1(muon1.p4());
        // double pt1= muon1.pt();
        for (unsigned int j = i + 1; j < nHighPtGlbMu; ++j) {
          reco::Muon muon2 = highPtGlbMuons[j];
          math::XYZTLorentzVector mu2(muon2.p4());
          // double pt2= muon2.pt();
          if (muon1.charge() == muon2.charge()) continue;
          math::XYZTLorentzVector pair = mu1 + mu2;
          double mass = pair.M();
          // checking isolation and hlt maching
          iso1 = muIso(muon1);
          iso2 = muIso(muon2);
          isMu1Iso = (iso1 < isoCut03_);
          isMu2Iso = (iso2 < isoCut03_);
          singleTrigFlag1 =
              IsMuMatchedToHLTMu(muon1, HLTMuMatched, maxDeltaR_, maxDPtRel_);
          singleTrigFlag2 =
              IsMuMatchedToHLTMu(muon2, HLTMuMatched, maxDeltaR_, maxDPtRel_);
          if (singleTrigFlag1 && singleTrigFlag2) isZGolden2HLT_ = true;
          if ((singleTrigFlag1 && !singleTrigFlag2) ||
              (!singleTrigFlag1 && singleTrigFlag2))
            isZGolden1HLT_ = true;
          // Z Golden passing all criteria, with 2 or 1 muon matched to an HLT
          // muon. Using the two cathegories as a control sample for the HLT
          // matching efficiency
          if (isZGolden2HLT_ || isZGolden1HLT_) {
            // a Z golden has been found, let's remove the two muons from the
            // list, dome for avoiding resolution effect enter muons in the
            // standalone and tracker collections.........
            highPtGlbMuons.erase(highPtGlbMuons.begin() + i);
            highPtGlbMuons.erase(highPtGlbMuons.begin() + j);
            // and updating the number of high pt muons....
            nHighPtGlbMu = highPtGlbMuons.size();

            if (isMu1Iso && isMu2Iso) {
              niso++;
              if (isZGolden2HLT_) {
                n2hlt++;
                mass2HLT_->Fill(mass);
                highMass2HLT_->Fill(mass);
                /*
                  if (pt1 > pt2) {
                  highest_mupt2HLT_ -> Fill (pt1);
                  lowest_mupt2HLT_ -> Fill (pt2);
                  } else {
                  highest_mupt2HLT_ -> Fill (pt2);
                  lowest_mupt2HLT_ -> Fill (pt1);
                  }
                */
              }
              if (isZGolden1HLT_) {
                n1hlt++;
                mass1HLT_->Fill(mass);
                highMass1HLT_->Fill(mass);
                /*
                  if (pt1 >pt2) {
                  highest_mupt1HLT_ -> Fill (pt1);
                  lowest_mupt1HLT_ -> Fill (pt2);
                } else {
                  highest_mupt1HLT_ -> Fill (pt2);
                  lowest_mupt1HLT_ -> Fill (pt1);
                }
                */
              }
            } else {
              // ZGlbGlb when at least one of the two muon is not isolated and
              // at least one HLT matched, used as control sample for the
              // isolation efficiency
              // filling events with one muon not isolated and both hlt mathced
              if (((isMu1Iso && !isMu2Iso) || (!isMu1Iso && isMu2Iso)) &&
                  (isZGolden2HLT_ && isZGolden1HLT_)) {
                isZGoldenNoIso_ = true;
                nNotIso++;
                massNotIso_->Fill(mass);
                highMassNotIso_->Fill(mass);
              }
              /*
                if (pt1 > pt2) {
                highest_muptNotIso_ -> Fill (pt1);
                lowest_muptNotIso_ -> Fill (pt2);
                } else {
                highest_muptNotIso_ -> Fill (pt2);
                lowest_muptNotIso_ -> Fill (pt1);
                }
              */
            }
          }
        }
      }
    }
    // W->MuNu selection (if at least one muon has been selected)
    // looking for a W if a Z is not found.... let's think if we prefer to
    // exclude zMuMuNotIso or zMuSta....
    if (!(isZGolden2HLT_ || isZGolden1HLT_)) {
      Handle<View<MET> > metCollection;
      if (!ev.getByToken(metToken_, metCollection)) {
        // LogError("") << ">>> MET collection does not exist !!!";
        return;
      }
      const MET& met = metCollection->at(0);
      nW = 0;
      for (unsigned int i = 0; i < nHighPtGlbMu; i++) {
        reco::Muon muon1 = highPtGlbMuons[i];
        /*
               quality cut not implemented
            Quality Cuts           double dxy =
gm->dxy(beamSpotHandle->position());
            Cut in 0.2           double trackerHits =
gm->hitPattern().numberOfValidTrackerHits();
            Cut in 11           bool quality = fabs(dxy)<dxyCut_  &&
muon::isGoodMuon(mu,muon::GlobalMuonPromptTight) && trackerHits>=trackerHitsCut_
&& mu.isTrackerMuon() ;
if (!quality) continue;
        */
        // isolation cut and hlt maching
        iso1 = muIso(muon1);
        isMu1Iso = (iso1 < isoCut03_);
        if (!isMu1Iso) continue;
        // checking if muon is matched to any HLT muon....
        singleTrigFlag1 =
            IsMuMatchedToHLTMu(muon1, HLTMuMatched, maxDeltaR_, maxDPtRel_);
        if (!singleTrigFlag1) continue;
        // std::cout << " is GlobMu macthecd to HLT: " << singleTrigFlag1 <<
        // std::endl;
        // MT cuts
        double met_px = met.px();
        double met_py = met.py();

        if (!metIncludesMuons_) {
          for (unsigned int i = 0; i < nHighPtGlbMu; i++) {
            reco::Muon muon1 = highPtGlbMuons[i];
            met_px -= muon1.px();
            met_py -= muon1.py();
          }
        }
        double met_et = met.pt();  // sqrt(met_px*met_px+met_py*met_py);
        LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met_px
                     << ", " << met_py << " [GeV]";
        double w_et = met_et + muon1.pt();
        double w_px = met_px + muon1.px();
        double w_py = met_py + muon1.py();
        double massT = w_et * w_et - w_px * w_px - w_py * w_py;
        massT = (massT > 0) ? sqrt(massT) : 0;
        // Acoplanarity cuts
        Geom::Phi<double> deltaphi(muon1.phi() - atan2(met_py, met_px));
        double acop = M_PI - fabs(deltaphi.value());
        // std::cout << " is acp of W candidate less then cut: " << (acop<
        // acopCut_) << std::endl;
        if (acop > acopCut_) continue;  // Cut in 2.0
        TMass_->Fill(massT);
        if (massT < mtMin_ || massT > mtMax_) continue;  // Cut in (50,200) GeV
        // we give the number of W only in the tMass selected but we have a look
        // at mass tails to check the QCD background
        isW_ = true;
        nW++;
        nTMass++;
      }
    }
    // if a zGlobal is not selected, look at the dimuonGlobalOneStandAlone and
    // dimuonGlobalOneTrack...., used as a control sample for the tracking
    // efficiency  and muon system efficiency
    if (!(isZGolden2HLT_ || isZGolden1HLT_ || isZGoldenNoIso_)) {
      for (unsigned int i = 0; i < nHighPtGlbMu; ++i) {
        reco::Muon glbMuon = highPtGlbMuons[i];
        math::XYZTLorentzVector mu1(glbMuon.p4());
        // double pt1= glbMuon.pt();
        // checking that the global muon is hlt matched otherwise skip the event
        singleTrigFlag1 =
            IsMuMatchedToHLTMu(glbMuon, HLTMuMatched, maxDeltaR_, maxDPtRel_);
        if (!singleTrigFlag1) continue;
        // checking that the global muon is isolated matched otherwise skip the
        // event
        iso1 = muIso(glbMuon);
        isMu1Iso = (iso1 < isoCut03_);
        if (!isMu1Iso) continue;
        // look at the standalone muon ....
        // stop the loop after 10 cicles....
        (nHighPtStaMu > 10) ? nHighPtStaMu = 10 : 1;
        for (unsigned int j = 0; j < nHighPtStaMu; ++j) {
          reco::Muon staMuon = highPtStaMuons[j];
          math::XYZTLorentzVector mu2(staMuon.p4());
          // double pt2= staMuon.pt();
          if (glbMuon.charge() == staMuon.charge()) continue;
          math::XYZTLorentzVector pair = mu1 + mu2;
          double mass = pair.M();
          iso2 = muIso(staMuon);
          isMu2Iso = (iso2 < isoCut03_);
          LogTrace("") << "\t... isolation value" << iso1 << ", isolated? "
                       << isMu1Iso;
          LogTrace("") << "\t... isolation value" << iso2 << ", isolated? "
                       << isMu2Iso;
          // requiring theat the global mu is mathed to the HLT  and both are
          // isolated
          if (isMu2Iso) {
            isZGlbSta_ = true;
            nGlbSta++;
            massGlbSta_->Fill(mass);
            highMassGlbSta_->Fill(mass);
            /*
              if (pt1 > pt2) {
              highest_muptGlbSta_ -> Fill (pt1);
              lowest_muptGlbSta_ -> Fill (pt2);
            } else {
              highest_muptGlbSta_ -> Fill (pt2);
              lowest_muptGlbSta_ -> Fill (pt1);
            }
            */
          }
        }
        // look at the tracks....
        Handle<TrackCollection> tracks;
        if (!ev.getByToken(trackToken_, tracks)) {
          // LogError("") << ">>> track collection does not exist !!!";
          return;
        }
        ev.getByToken(trackToken_, tracks);
        Handle<CaloTowerCollection> calotower;
        if (!ev.getByToken(caloTowerToken_, calotower)) {
          // LogError("") << ">>> calotower collection does not exist !!!";
          return;
        }
        ev.getByToken(caloTowerToken_, calotower);
        // avoid to loop on more than 5000 trks
        size_t nTrk = tracks->size();
        (nTrk > 5000) ? nTrk = 5000 : 1;
        for (unsigned int j = 0; j < nTrk; j++) {
          const reco::Track& tk = tracks->at(j);
          if (glbMuon.charge() == tk.charge()) continue;
          double pt2 = tk.pt();
          double eta = tk.eta();
          double dxy = tk.dxy(beamSpotHandle->position());
          if (pt2 < ptMuCut_ || fabs(eta) > etaMuCut_ || fabs(dxy) > dxyCut_)
            continue;
          // assuming that the track is a mu....
          math::XYZTLorentzVector mu2(
              tk.px(), tk.py(), tk.pz(),
              sqrt((tk.p() * tk.p()) + (0.10566 * 0.10566)));
          math::XYZTLorentzVector pair = mu1 + mu2;
          double mass = pair.M();
          // now requiring that the tracks is isolated.......
          iso2 = tkIso(tk, tracks, calotower);
          isMu2Iso = (iso2 < isoCut03_);
          //          std::cout << "found a track with rel/comb iso: " << iso2
          //<< std::endl;
          if (isMu2Iso) {
            isZGlbTrk_ = true;
            nGlbTrk++;
            massGlbTrk_->Fill(mass);
            highMassGlbTrk_->Fill(mass);
            if (!isW_) continue;
            massIsBothGlbTrkThanW_->Fill(mass);
            highMassIsBothGlbTrkThanW_->Fill(mass);
            /*
              if (pt1 > pt2) {
              highest_muptGlbTrk_ -> Fill (pt1);
              lowest_muptGlbTrk_ -> Fill (pt2);
              } else {
              highest_muptGlbTrk_ -> Fill (pt2);
              lowest_muptGlbTrk_ -> Fill (pt1);
            }
            */
          }
        }
      }
    }
  }
  if ((hlt_sel || isZGolden2HLT_ || isZGolden1HLT_ || isZGoldenNoIso_ ||
       isZGlbSta_ || isZGlbTrk_ || isW_)) {
    nsel++;
    LogTrace("") << ">>>> Event ACCEPTED";
  } else {
    LogTrace("") << ">>>> Event REJECTED";
  }
  return;
}

void EwkMuLumiMonitorDQM::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 106 of file EwkMuLumiMonitorDQM.cc.

References DQMStore::IBooker::book1D(), highMass1HLT_, highMass2HLT_, highMassGlbSta_, highMassGlbTrk_, highMassIsBothGlbTrkThanW_, highMassNotIso_, mass1HLT_, mass2HLT_, massGlbSta_, massGlbTrk_, massIsBothGlbTrkThanW_, massNotIso_, DQMStore::IBooker::setCurrentFolder(), and TMass_.

                                         {

  ibooker.setCurrentFolder("Physics/EwkMuLumiMonitorDQM");

  mass2HLT_ = ibooker.book1D("Z_2HLT_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
  mass1HLT_ = ibooker.book1D("Z_1HLT_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
  massNotIso_ = ibooker.book1D("Z_NOTISO_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
  massGlbSta_ = ibooker.book1D("Z_GLBSTA_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
  massGlbTrk_ = ibooker.book1D("Z_GLBTRK_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
  massIsBothGlbTrkThanW_ = ibooker.book1D("Z_ISBOTHGLBTRKTHANW_MASS",
      "Z mass [GeV/c^{2}]", 200, 0., 200.);

  highMass2HLT_ = ibooker.book1D("Z_2HLT_HIGHMASS",
      "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
  highMass1HLT_ = ibooker.book1D("Z_1HLT_HIGHMASS",
      "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
  highMassNotIso_ = ibooker.book1D("Z_NOTISO_HIGHMASS",
      "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
  highMassGlbSta_ = ibooker.book1D("Z_GLBSTA_HIGHMASS",
      "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
  highMassGlbTrk_ = ibooker.book1D("Z_GLBTRK_HIGHMASS",
      "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
  highMassIsBothGlbTrkThanW_ = ibooker.book1D("Z_ISBOTHGLBTRKTHANW_HIGHMASS",
      "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);

  TMass_ = ibooker.book1D("TMASS", "Transverse mass [GeV]", 300, 0., 300.);
}

void EwkMuLumiMonitorDQM::dqmBeginRun ( const edm::Run &  r, const edm::EventSetup &  iSetup  )

overridevirtual

Definition at line 81 of file EwkMuLumiMonitorDQM.cc.

References hltConfigProvider_, HLTConfigProvider::init(), isValidHltConfig_, n1hlt, n2hlt, nall, nEvWithHighPtMu, nGlbSta, nGlbTrk, nhlt, nInKinRange, niso, nNotIso, nsel, nTMass, nW, edm::InputTag::process(), and trigTag_.

                                                                            {
  nall = 0;
  nEvWithHighPtMu = 0;
  nInKinRange = 0;
  nsel = 0;
  niso = 0;
  nhlt = 0;
  n1hlt = 0;
  n2hlt = 0;
  nNotIso = 0;
  nGlbSta = 0;
  nGlbTrk = 0;
  nTMass = 0;
  nW = 0;

  // passed as parameter to HLTConfigProvider::init(), not yet used
  bool isConfigChanged = false;

  // isValidHltConfig_ used to short-circuit analyze() in case of problems
  isValidHltConfig_ =
      hltConfigProvider_.init(r, iSetup, trigTag_.process(), isConfigChanged);
  // std::cout << "hlt config trigger is valid??" << isValidHltConfig_ <<
  // std::endl;
}

void EwkMuLumiMonitorDQM::endRun ( const edm::Run &  r, const edm::EventSetup &    )

overridevirtual

Definition at line 197 of file EwkMuLumiMonitorDQM.cc.

References n1hlt, n2hlt, nall, nGlbSta, nhlt, nNotIso, nsel, and nTMass.

                                                                {

  LogVerbatim("") << "\n>>>>>> Z/W SELECTION SUMMARY BEGIN >>>>>>>>>>>>>>>";
  LogVerbatim("") << "Total numer of events analyzed: " << nall << " [events]";
  LogVerbatim("") << "Total numer of events selected: " << nsel << " [events]";

  LogVerbatim("") << "Passing HLT criteria:           " << nhlt << " [events] ";
  LogVerbatim("") << "Passing 2 HLT match criteria:           " << n2hlt
                  << " [events] ";
  LogVerbatim("") << "Passing 1 HLT match criteria:           " << n1hlt
                  << " [events] ";
  LogVerbatim("") << "Passing Not Iso criteria:           " << nNotIso
                  << " [events] ";
  LogVerbatim("") << "Passing GlbSta criteria:           " << nGlbSta
                  << " [events] ";
  LogVerbatim("") << "Passing W criteria:           " << nTMass << " [events] ";
  LogVerbatim("") << ">>>>>> Z/W SELECTION SUMMARY END   >>>>>>>>>>>>>>>\n";
}

void EwkMuLumiMonitorDQM::init_histograms

(

)

bool EwkMuLumiMonitorDQM::IsMuMatchedToHLTMu	(	const reco::Muon &	mu,
		const std::vector< reco::Particle > &	HLTMu,
		double	DR,
		double	DPtRel
	)

Definition at line 182 of file EwkMuLumiMonitorDQM.cc.

References deltaR(), gen::k, nJettinessAdder_cfi::nPass, EnergyCorrector::pt, and reco::LeafCandidate::pt().

Referenced by analyze().

                   {
  size_t dim = HLTMu.size();
  size_t nPass = 0;
  if (dim == 0) return false;
  for (size_t k = 0; k < dim; k++) {
    if ((deltaR(HLTMu[k], mu) < DR) &&
        (fabs(HLTMu[k].pt() - mu.pt()) / HLTMu[k].pt() < DPtRel)) {
      nPass++;
    }
  }
  return (nPass > 0);
}

double EwkMuLumiMonitorDQM::muIso

(

const reco::Muon &

mu

)

Definition at line 135 of file EwkMuLumiMonitorDQM.cc.

References reco::MuonIsolation::emEt, reco::MuonIsolation::hadEt, isCombinedIso_, reco::Muon::isolationR03(), isRelativeIso_, reco::LeafCandidate::pt(), and reco::MuonIsolation::sumPt.

Referenced by analyze().

                                                    {
  double isovar = mu.isolationR03().sumPt;
  if (isCombinedIso_) {
    isovar += mu.isolationR03().emEt;
    isovar += mu.isolationR03().hadEt;
  }
  if (isRelativeIso_) isovar /= mu.pt();
  return isovar;
}

double EwkMuLumiMonitorDQM::tkIso	(	const reco::Track &	tk,
		edm::Handle< reco::TrackCollection >	tracks,
		edm::Handle< CaloTowerCollection >	calotower
	)

Definition at line 145 of file EwkMuLumiMonitorDQM.cc.

References edm::SortedCollection< T, SORT >::begin(), deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, PVValHelper::dz, HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), reco::TrackBase::eta(), mps_fire::i, isCombinedIso_, isRelativeIso_, reco::Track::outerEta(), reco::Track::outerPhi(), reco::TrackBase::phi(), reco::TrackBase::pt(), ptThreshold_, mathSSE::sqrt(), reco::TrackBase::vx(), reco::TrackBase::vy(), and reco::TrackBase::vz().

Referenced by analyze().

                                                                         {
  double ptSum = 0;
  for (size_t i = 0; i < tracks->size(); ++i) {
    const reco::Track& elem = tracks->at(i);
    double elemPt = elem.pt();
    // same parameter used for muIsolation: dR [0.01, IsoCut03_], |dZ|<0.2,
    // |d_r(xy)|<0.1
    double elemVx = elem.vx();
    double elemVy = elem.vy();
    double elemD0 = sqrt(elemVx * elemVx + elemVy * elemVy);
    if (elemD0 > 0.2) continue;
    double dz = fabs(elem.vz() - tk.vz());
    if (dz > 0.1) continue;
    // evaluate only for tracks with pt>ptTreshold
    if (elemPt < ptThreshold_) continue;
    double dR = deltaR(elem.eta(), elem.phi(), tk.eta(), tk.phi());
    // isolation in a cone with dR=0.3, and vetoing the track itself
    if ((dR < 0.01) || (dR > 0.3)) continue;
    ptSum += elemPt;
  }
  if (isCombinedIso_) {
    // loop on clusters....
    for (CaloTowerCollection::const_iterator it = calotower->begin();
         it != calotower->end(); ++it) {
      double dR = deltaR(it->eta(), it->phi(), tk.outerEta(), tk.outerPhi());
      // veto value is 0.1 for towers....
      if ((dR < 0.1) || (dR > 0.3)) continue;
      ptSum += it->emEnergy();
      ptSum += it->hadEnergy();
    }
  }
  if (isRelativeIso_) ptSum /= tk.pt();
  return ptSum;
}

Member Data Documentation

double EwkMuLumiMonitorDQM::acopCut_

private

Definition at line 82 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

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

private

Definition at line 58 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

edm::EDGetTokenT<CaloTowerCollection> EwkMuLumiMonitorDQM::caloTowerToken_

private

Definition at line 61 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

double EwkMuLumiMonitorDQM::deltaRTrk_

private

Definition at line 75 of file EwkMuLumiMonitorDQM.h.

double EwkMuLumiMonitorDQM::deltaRVetoTrk_

private

Definition at line 77 of file EwkMuLumiMonitorDQM.h.

double EwkMuLumiMonitorDQM::dxyCut_

private

Definition at line 83 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

double EwkMuLumiMonitorDQM::etaMuCut_

private

Definition at line 69 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

MonitorElement* EwkMuLumiMonitorDQM::highMass1HLT_

private

Definition at line 91 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::highMass2HLT_

private

Definition at line 86 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::highMassGlbSta_

private

Definition at line 101 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::highMassGlbTrk_

private

Definition at line 106 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::highMassIsBothGlbTrkThanW_

private

Definition at line 113 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::highMassNotIso_

private

Definition at line 96 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

HLTConfigProvider EwkMuLumiMonitorDQM::hltConfigProvider_

private

Definition at line 137 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

bool EwkMuLumiMonitorDQM::isCombinedIso_

private

Definition at line 72 of file EwkMuLumiMonitorDQM.h.

Referenced by muIso(), and tkIso().

double EwkMuLumiMonitorDQM::isoCut03_

private

Definition at line 73 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

bool EwkMuLumiMonitorDQM::isRelativeIso_

private

Definition at line 71 of file EwkMuLumiMonitorDQM.h.

Referenced by muIso(), and tkIso().

bool EwkMuLumiMonitorDQM::isValidHltConfig_

private

Definition at line 136 of file EwkMuLumiMonitorDQM.h.

Referenced by dqmBeginRun(), and EwkMuLumiMonitorDQM().

bool EwkMuLumiMonitorDQM::isW_

private

Definition at line 134 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

bool EwkMuLumiMonitorDQM::isZGlbSta_

private

Definition at line 132 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

bool EwkMuLumiMonitorDQM::isZGlbTrk_

private

Definition at line 133 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

bool EwkMuLumiMonitorDQM::isZGolden1HLT_

private

Definition at line 129 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

bool EwkMuLumiMonitorDQM::isZGolden2HLT_

private

Definition at line 130 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

bool EwkMuLumiMonitorDQM::isZGoldenNoIso_

private

Definition at line 131 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

MonitorElement* EwkMuLumiMonitorDQM::mass1HLT_

private

Definition at line 90 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::mass2HLT_

private

Definition at line 85 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::massGlbSta_

private

Definition at line 100 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::massGlbTrk_

private

Definition at line 105 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::massIsBothGlbTrkThanW_

private

Definition at line 112 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EwkMuLumiMonitorDQM::massNotIso_

private

Definition at line 95 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

double EwkMuLumiMonitorDQM::maxDeltaR_

private

Definition at line 79 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

double EwkMuLumiMonitorDQM::maxDPtRel_

private

Definition at line 78 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

bool EwkMuLumiMonitorDQM::metIncludesMuons_

private

Definition at line 63 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

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

private

Definition at line 62 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

double EwkMuLumiMonitorDQM::mtMax_

private

Definition at line 81 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

double EwkMuLumiMonitorDQM::mtMin_

private

Definition at line 80 of file EwkMuLumiMonitorDQM.h.

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

private

Definition at line 59 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

unsigned int EwkMuLumiMonitorDQM::n1hlt

private

Definition at line 121 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::n2hlt

private

Definition at line 122 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::nall

private

Definition at line 115 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::nEvWithHighPtMu

private

Definition at line 116 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

unsigned int EwkMuLumiMonitorDQM::nGlbSta

private

Definition at line 124 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::nGlbTrk

private

Definition at line 125 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

unsigned int EwkMuLumiMonitorDQM::nhlt

private

Definition at line 120 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::nInKinRange

private

Definition at line 117 of file EwkMuLumiMonitorDQM.h.

Referenced by dqmBeginRun().

unsigned int EwkMuLumiMonitorDQM::niso

private

Definition at line 119 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

unsigned int EwkMuLumiMonitorDQM::nNotIso

private

Definition at line 123 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::nsel

private

Definition at line 118 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::nTMass

private

Definition at line 126 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), dqmBeginRun(), and endRun().

unsigned int EwkMuLumiMonitorDQM::nW

private

Definition at line 127 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

double EwkMuLumiMonitorDQM::ptMuCut_

private

Definition at line 68 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

double EwkMuLumiMonitorDQM::ptThreshold_

private

Definition at line 76 of file EwkMuLumiMonitorDQM.h.

Referenced by tkIso().

MonitorElement* EwkMuLumiMonitorDQM::TMass_

private

Definition at line 110 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<reco::TrackCollection> EwkMuLumiMonitorDQM::trackToken_

private

Definition at line 60 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

edm::EDGetTokenT<trigger::TriggerEvent> EwkMuLumiMonitorDQM::trigEvToken_

private

Definition at line 57 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

edm::InputTag EwkMuLumiMonitorDQM::trigTag_

private

Definition at line 55 of file EwkMuLumiMonitorDQM.h.

Referenced by dqmBeginRun().

edm::EDGetTokenT<edm::TriggerResults> EwkMuLumiMonitorDQM::trigToken_

private

Definition at line 56 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().