EwkMuDQM Class Reference

#include <EwkMuDQM.h>

Inheritance diagram for EwkMuDQM:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	beginRun (const edm::Run &, const edm::EventSetup &)
virtual void	endJob ()
virtual void	endRun (const edm::Run &, const edm::EventSetup &)
	EwkMuDQM (const edm::ParameterSet &)
void	init_histograms ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Attributes
MonitorElement *	acop_after_
MonitorElement *	acop_before_
double	acopCut_
MonitorElement *	chi2_after_
MonitorElement *	chi2_before_
MonitorElement *	dimuonmass_after_
MonitorElement *	dimuonmass_before_
MonitorElement *	dimuonSAmass_after_
MonitorElement *	dimuonSAmass_before_
MonitorElement *	dimuonSASAmass_after_
MonitorElement *	dimuonSASAmass_before_
MonitorElement *	dxy_after_
MonitorElement *	dxy_before_
double	dxyCut_
double	eJetMin_
MonitorElement *	eta_after_
MonitorElement *	eta_before_
double	etaCut_
MonitorElement *	goodewkmuon_after_
MonitorElement *	goodewkmuon_before_
HLTConfigProvider	hltConfigProvider_
bool	isAlsoTrackerMuon_
bool	isCombinedIso_
MonitorElement *	iso_after_
MonitorElement *	iso_before_
double	isoCut03_
bool	isRelativeIso_
bool	isValidHltConfig_
edm::InputTag	jetTag_
MonitorElement *	met_after_
MonitorElement *	met_before_
double	metMax_
double	metMin_
edm::InputTag	metTag_
MonitorElement *	mt_after_
MonitorElement *	mt_before_
double	mtMax_
double	mtMin_
MonitorElement *	muoncharge_after_
MonitorElement *	muoncharge_before_
MonitorElement *	muonhits_after_
MonitorElement *	muonhits_before_
int	muonHitsCut_
edm::InputTag	muonTag_
unsigned int	nall
MonitorElement *	nhits_after_
MonitorElement *	nhits_before_
unsigned int	nhlt
unsigned int	niso
int	nJetMax_
MonitorElement *	njets_after_
MonitorElement *	njets_before_
int	nMatchesCut_
unsigned int	nmet
double	normalizedChi2Cut_
MonitorElement *	npvs_after_
MonitorElement *	npvs_before_
unsigned int	nrec
unsigned int	nsel
MonitorElement *	nz1_after_
MonitorElement *	nz1_before_
MonitorElement *	nz2_after_
MonitorElement *	nz2_before_
int	pixelHitsCut_
MonitorElement *	pt_after_
MonitorElement *	pt_before_
double	ptCut_
MonitorElement *	ptmuonZ_after_
double	ptThrForZ1_
double	ptThrForZ2_
DQMStore *	theDbe
MonitorElement *	tkmu_after_
MonitorElement *	tkmu_before_
int	trackerHitsCut_
MonitorElement *	trig_after_
MonitorElement *	trig_before_
edm::InputTag	trigTag_
edm::InputTag	vertexTag_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

DQM offline for EWKMu

Definition at line 19 of file EwkMuDQM.h.

Constructor & Destructor Documentation

EwkMuDQM::EwkMuDQM

(

const edm::ParameterSet &

cfg

)

Definition at line 36 of file EwkMuDQM.cc.

References isValidHltConfig_.

                                             :
      // Input collections
      trigTag_(cfg.getUntrackedParameter<edm::InputTag> ("TrigTag", edm::InputTag("TriggerResults::HLT"))),
      muonTag_(cfg.getUntrackedParameter<edm::InputTag> ("MuonTag", edm::InputTag("muons"))),
      metTag_(cfg.getUntrackedParameter<edm::InputTag> ("METTag", edm::InputTag("pfmet"))),
      jetTag_(cfg.getUntrackedParameter<edm::InputTag> ("JetTag", edm::InputTag("ak5PFJets"))),
      vertexTag_(cfg.getUntrackedParameter<edm::InputTag> ("VertexTag", edm::InputTag("offlinePrimaryVertices"))),

      // Main cuts 
      ptCut_(cfg.getUntrackedParameter<double>("PtCut", 25.)),
      etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.1)),
      isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso", true)),
      isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso", false)),
      isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03", 0.1)),
      mtMin_(cfg.getUntrackedParameter<double>("MtMin", 50.)),
      mtMax_(cfg.getUntrackedParameter<double>("MtMax", 200.)),
      metMin_(cfg.getUntrackedParameter<double>("MetMin", -999999.)),
      metMax_(cfg.getUntrackedParameter<double>("MetMax", 999999.)),
      acopCut_(cfg.getUntrackedParameter<double>("AcopCut", 999.)),

      // Muon quality cuts
      dxyCut_(cfg.getUntrackedParameter<double>("DxyCut", 0.2)),   // dxy < 0.2 cm 
      normalizedChi2Cut_(cfg.getUntrackedParameter<double>("NormalizedChi2Cut", 10.)), // chi2/ndof (of global fit) <10.0
      trackerHitsCut_(cfg.getUntrackedParameter<int>("TrackerHitsCut", 11)),  // Tracker Hits >10 
      pixelHitsCut_(cfg.getUntrackedParameter<int>("PixelHitsCut", 1)), // Pixel Hits >0
      muonHitsCut_(cfg.getUntrackedParameter<int>("MuonHitsCut", 1)),  // Valid Muon Hits >0 
      isAlsoTrackerMuon_(cfg.getUntrackedParameter<bool>("IsAlsoTrackerMuon", true)),
      nMatchesCut_(cfg.getUntrackedParameter<int>("NMatchesCut", 2)), // At least 2 Chambers with matches 

      // 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))
{
  isValidHltConfig_ = false;
}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 192 of file EwkMuDQM.cc.

References acop_after_, acop_before_, acopCut_, DeDxDiscriminatorTools::charge(), reco::LeafCandidate::charge(), dimuonmass_after_, dimuonmass_before_, dxy_after_, dxy_before_, dxyCut_, eJetMin_, reco::MuonIsolation::emEt, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), eta(), eta_after_, eta_before_, etaCut_, MonitorElement::Fill(), newFWLiteAna::found, edm::Event::getByLabel(), reco::Muon::globalTrack(), goodewkmuon_after_, goodewkmuon_before_, reco::MuonIsolation::hadEt, hltConfigProvider_, i, reco::Muon::innerTrack(), isCombinedIso_, reco::Muon::isGlobalMuon(), edm::Ref< C, T, F >::isNull(), iso_after_, iso_before_, isoCut03_, reco::Muon::isolationR03(), isRelativeIso_, reco::Muon::isTrackerMuon(), reco::Vertex::isValid(), j, metsig::jet, jetTag_, LogTrace, M_PI, CaloMET_cfi::met, met_after_, met_before_, metMax_, metMin_, metTag_, mt_after_, mt_before_, mtMax_, mtMin_, muoncharge_after_, muoncharge_before_, muonHitsCut_, muonTag_, nall, nJetMax_, njets_after_, njets_before_, nMatchesCut_, normalizedChi2Cut_, npvs_after_, npvs_before_, reco::Muon::numberOfMatches(), nz1_after_, nz1_before_, nz2_after_, nz2_before_, reco::LeafCandidate::p(), reco::LeafCandidate::phi(), pixelHitsCut_, HLTConfigProvider::prescaleSize(), HLTConfigProvider::prescaleValue(), reco::LeafCandidate::pt(), pt_after_, pt_before_, ptCut_, ptmuonZ_after_, ptThrForZ1_, ptThrForZ2_, reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), mathSSE::sqrt(), reco::MuonIsolation::sumPt, trackerHitsCut_, trig_after_, trig_before_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), trigNames, trigTag_, Geom::Phi< T >::value(), GoodVertex_cfg::vertexCollection, and vertexTag_.

                                                                 {
      
      // 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.getByLabel(muonTag_, muonCollection)) {
    //LogWarning("") << ">>> Muon collection does not exist !!!";
    return;
      }
      unsigned int muonCollectionSize = muonCollection->size();

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


      // Loop to reject/control Z->mumu is done separately
      unsigned int nmuonsForZ1 = 0;
      unsigned int nmuonsForZ2 = 0;
      bool cosmic = false;
      for (unsigned int i=0; i<muonCollectionSize; i++) {
            const Muon& mu = muonCollection->at(i);
            if (!mu.isGlobalMuon()) continue;
            double pt = mu.pt();
            double dxy = mu.innerTrack()->dxy(beamSpotHandle->position());

            if (fabs(dxy)>1) { cosmic=true; break;} 

            if (pt>ptThrForZ1_) nmuonsForZ1++;
            if (pt>ptThrForZ2_) nmuonsForZ2++;

            for (unsigned int j=i+1; j<muonCollectionSize; j++) {
                 const Muon& mu2 = muonCollection->at(j);
                 if (mu2.isGlobalMuon() && (mu.charge()*mu2.charge()==-1) ){
                         const math::XYZTLorentzVector ZRecoGlb (mu.px()+mu2.px(), mu.py()+mu2.py() , mu.pz()+mu2.pz(), mu.p()+mu2.p());
                         dimuonmass_before_->Fill(ZRecoGlb.mass());
                 }
            }
      }
     if(cosmic) return;

      LogTrace("") << "> Z rejection: muons above " << ptThrForZ1_ << " [GeV]: " << nmuonsForZ1;
      LogTrace("") << "> Z rejection: muons above " << ptThrForZ2_ << " [GeV]: " << nmuonsForZ2;
      nz1_before_->Fill(nmuonsForZ1);
      nz2_before_->Fill(nmuonsForZ2);
  
      // MET
      Handle<View<MET> > metCollection;
      if (!ev.getByLabel(metTag_, metCollection)) {
    //LogWarning("") << ">>> MET collection does not exist !!!";
    return;
      }
      const MET& met = metCollection->at(0);
      double met_et = met.pt();
      LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met.px() << ", " << met.py() << " [GeV]";
      met_before_->Fill(met_et);

      // Vertices in the event
      Handle<View<reco::Vertex> > vertexCollection;
           if (!ev.getByLabel(vertexTag_, vertexCollection)) {
                 LogError("") << ">>> Vertex collection does not exist !!!";
                 return;
            }
      unsigned int vertexCollectionSize = vertexCollection->size();

      

      int nvvertex = 0;
      for (unsigned int i=0; i<vertexCollectionSize; i++) {
            const Vertex& vertex = vertexCollection->at(i);
            if (vertex.isValid()) nvvertex++;
      }

      npvs_before_->Fill(nvvertex);

      bool trigger_fired = false;
      Handle<TriggerResults> triggerResults;
      if (!ev.getByLabel(trigTag_, triggerResults)) {
    //LogWarning("") << ">>> TRIGGER collection does not exist !!!";
    return;
      }
      const edm::TriggerNames & trigNames = ev.triggerNames(*triggerResults);

      for (unsigned int i=0; i<triggerResults->size(); i++)
      {
              const std::string trigName = trigNames.triggerName(i);
              size_t found = trigName.find("HLT_Mu");
              if ( found == std::string::npos) continue;

              bool prescaled=false;    
              for (unsigned int ps= 0; ps<  hltConfigProvider_.prescaleSize(); ps++){
                  const unsigned int prescaleValue = hltConfigProvider_.prescaleValue(ps, trigName) ;
                  if (prescaleValue != 1) prescaled =true;
              }
              if(prescaled) continue;    

              if( triggerResults->accept(i) )   trigger_fired=true;
      }     
      trig_before_->Fill(trigger_fired);






      // Jet collection
      Handle<View<Jet> > jetCollection;
      if (!ev.getByLabel(jetTag_, jetCollection)) {
    //LogError("") << ">>> JET collection does not exist !!!";
    return;
      }
      unsigned int jetCollectionSize = jetCollection->size();
      int njets = 0;
      for (unsigned int i=0; i<jetCollectionSize; i++) {
            const Jet& jet = jetCollection->at(i);
                  double minDistance=99999; // This is in order to use PFJets
                  for (unsigned int i=0; i<muonCollectionSize; i++) {
                        const Muon& mu = muonCollection->at(i);
                        double distance = sqrt( (mu.eta()-jet.eta())*(mu.eta()-jet.eta()) +(mu.phi()-jet.phi())*(mu.phi()-jet.phi()) );      
                        if (minDistance>distance) minDistance=distance;
                  }
                  if (minDistance<0.3) continue; // 0.3 is the isolation cone around the muon
            if (jet.et()>eJetMin_) njets++;
      }
      LogTrace("") << ">>> Total number of jets: " << jetCollectionSize;
      LogTrace("") << ">>> Number of jets above " << eJetMin_ << " [GeV]: " << njets;
      njets_before_->Fill(njets);

      // Start counting
      nall++;

      // 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;
      bool pv_hist_done = false;
      bool charge_hist_done = false;

      // Central W->mu nu selection criteria
      const int NFLAGS = 11;
      bool muon_sel[NFLAGS];
      bool muon4Z=false;

      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; 
            if (fabs(eta)<etaCut_) muon_sel[1] = true; 

            double charge=mu.charge();

            // d0, chi2, nhits quality cuts
            double dxy = gm->dxy(beamSpotHandle->position());
            double normalizedChi2 = gm->normalizedChi2();
            double trackerHits = tk->hitPattern().numberOfValidTrackerHits();
            int pixelHits = tk->hitPattern().numberOfValidPixelHits();
            int muonHits = gm->hitPattern().numberOfValidMuonHits();
            int nMatches = mu.numberOfMatches();

            LogTrace("") << "\t... dxy, normalizedChi2, trackerHits, isTrackerMuon?: " << dxy << " [cm], " << normalizedChi2 << ", " << trackerHits << ", " << mu.isTrackerMuon();
            if (fabs(dxy)<dxyCut_) muon_sel[2] = true; 

            bool quality=true;
            
            if (normalizedChi2>normalizedChi2Cut_) quality =false; 
            if (trackerHits<trackerHitsCut_) quality =false;
            if (pixelHits<pixelHitsCut_) quality =false;
            if (muonHits<muonHitsCut_) quality=false;;
            if (!mu.isTrackerMuon()) quality=false;
            if (nMatches<nMatchesCut_) quality=false;
            muon_sel[3]=quality;


            pt_before_->Fill(pt);
            eta_before_->Fill(eta);
            dxy_before_->Fill(dxy);
            muoncharge_before_->Fill(charge);
            goodewkmuon_before_->Fill(quality);

            // 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[4] = true; 

            LogTrace("") << "\t... isolation value" << isovar <<", isolated? " << muon_sel[6];
            iso_before_->Fill(isovar);


            // HLT (not mtched to muon for the time being)
            if (trigger_fired) muon_sel[5] = true; 

            // For Z:
            if (pt>ptThrForZ1_ && fabs(eta)<etaCut_ && fabs(dxy)<dxyCut_ && quality && trigger_fired && isovar<isoCut03_) { muon4Z = true;}


            // 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[6] = true; 
            mt_before_->Fill(massT);
            if (met_et>metMin_ && met_et<metMax_) muon_sel[7] = true; 

            // 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[8] = true; 
            acop_before_->Fill(acop);

            // Remaining flags (from global event information)
            if (nmuonsForZ1<1 || nmuonsForZ2<2) muon_sel[9] = true; 
            if (njets<=nJetMax_) muon_sel[10] = true; 

            // 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<4 && !muon_sel[j]) rec_sel_this = false;
                  if (j<5 && !muon_sel[j]) iso_sel_this = false;
                  if (j<6 && !muon_sel[j]) hlt_sel_this = false;
                  if (j<8 && !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) 
                        pt_after_->Fill(pt);
                  if (!muon_sel[1] || flags_passed==NFLAGS) 
                        eta_after_->Fill(eta);
                  if (!muon_sel[2] || flags_passed==NFLAGS) 
                        dxy_after_->Fill(dxy);
                  if (!muon_sel[3] || flags_passed==NFLAGS)
                        goodewkmuon_after_->Fill(quality);
                  if (!muon_sel[4] || flags_passed==NFLAGS) 
                        iso_after_->Fill(isovar);
                  if (!muon_sel[5] || flags_passed==NFLAGS) 
                        if (!hlt_hist_done) trig_after_->Fill(trigger_fired);
                        hlt_hist_done = true;
                  if (!muon_sel[6] || flags_passed==NFLAGS) 
                        mt_after_->Fill(massT);
                  if (!muon_sel[7] || flags_passed==NFLAGS) 
                        if (!met_hist_done) met_after_->Fill(met_et);
                        met_hist_done = true;
                  if (!muon_sel[8] || flags_passed==NFLAGS) 
                        acop_after_->Fill(acop);
                  if (!muon_sel[9] || flags_passed==NFLAGS) 
                        if (!nz1_hist_done) nz1_after_->Fill(nmuonsForZ1);
                        nz1_hist_done = true;
                  if (!muon_sel[9] || flags_passed==NFLAGS) 
                        if (!nz2_hist_done) nz2_after_->Fill(nmuonsForZ2);
                        nz2_hist_done = true;
                  if (!muon_sel[10] || flags_passed==NFLAGS) { 
                        if (!njets_hist_done) njets_after_->Fill(njets);
                        njets_hist_done = true;
                  if( flags_passed==NFLAGS ) {
                        if (!pv_hist_done) npvs_after_->Fill(nvvertex);
                        if (!charge_hist_done) muoncharge_after_->Fill(charge);
                     }
                        pv_hist_done=true; charge_hist_done=true;   
                  }  
            }


            // The cases in which the event is rejected as a Z are considered independently:
            if ( muon4Z &&  !muon_sel[9]){
                   // Plots for 2 muons       
                   bool usedMuon=false;
                   for (unsigned int j=i+1; j<muonCollectionSize; j++) {
                         const Muon& mu2 = muonCollection->at(j);
                              if (!mu2.isGlobalMuon()) continue;
                              if (mu2.charge() * charge != -1 ) continue;
                                    double pt2 = mu2.pt(); if (pt2<=ptThrForZ1_) continue;
                                    double eta2=mu2.eta(); if (fabs(eta2)>etaCut_) continue;
                                    double isovar2 = mu2.isolationR03().sumPt; 
                                    if (isCombinedIso_) {
                                          isovar2 += mu2.isolationR03().emEt;
                                          isovar2 += mu2.isolationR03().hadEt;
                                    }
                                    if (isRelativeIso_) isovar2 /= pt2;
                                    if (isovar2>=isoCut03_) continue;
                               const math::XYZTLorentzVector ZRecoGlb (mu.px()+mu2.px(), mu.py()+mu2.py() , mu.pz()+mu2.pz(), mu.p()+mu2.p());
                               dimuonmass_after_->Fill(ZRecoGlb.mass());
                               if(!usedMuon){ptmuonZ_after_->Fill(mu.pt()); usedMuon=true;}
                  }
            }



      }




      

      return;

}

void EwkMuDQM::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 94 of file EwkMuDQM.cc.

References init_histograms(), cmsCodeRules.cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), and theDbe.

                        {
      theDbe = Service<DQMStore>().operator->();
      theDbe->setCurrentFolder("Physics/EwkMuDQM");
      init_histograms();
}

void EwkMuDQM::beginRun ( const edm::Run &  r, const edm::EventSetup &  iSet  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 76 of file EwkMuDQM.cc.

References hltConfigProvider_, HLTConfigProvider::init(), isValidHltConfig_, nall, nhlt, niso, nmet, nrec, and nsel.

                                                            {
      nall = 0;
      nsel = 0;

      nrec = 0; 
      niso = 0; 
      nhlt = 0; 
      nmet = 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, iSet, "HLT", isConfigChanged );

}

void EwkMuDQM::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 185 of file EwkMuDQM.cc.

                      {
}

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 188 of file EwkMuDQM.cc.

                                                          {

}

void EwkMuDQM::init_histograms

(

)

Definition at line 100 of file EwkMuDQM.cc.

References acop_after_, acop_before_, DQMStore::book1D(), dimuonmass_after_, dimuonmass_before_, dxy_after_, dxy_before_, eJetMin_, eta_after_, eta_before_, goodewkmuon_after_, goodewkmuon_before_, i, isCombinedIso_, iso_after_, iso_before_, isRelativeIso_, jetTag_, edm::InputTag::label(), M_PI, met_after_, met_before_, metTag_, mt_after_, mt_before_, muoncharge_after_, muoncharge_before_, njets_after_, njets_before_, npvs_after_, npvs_before_, nz1_after_, nz1_before_, nz2_after_, nz2_before_, pt_after_, pt_before_, ptmuonZ_after_, ptThrForZ1_, ptThrForZ2_, theDbe, trig_after_, and trig_before_.

Referenced by beginJob().

                               {

      char chtitle[256] = "";
      for (int i=0; i<2; ++i) {
            snprintf(chtitle, 255, "Muon transverse momentum (global muon) [GeV]");
            pt_before_ = theDbe->book1D("PT_BEFORECUTS",chtitle,100,0.,100.);
            pt_after_ = theDbe->book1D("PT_LASTCUT",chtitle,100,0.,100.);

            snprintf(chtitle, 255, "Muon pseudo-rapidity");
            eta_before_ = theDbe->book1D("ETA_BEFORECUTS",chtitle,50,-2.5,2.5);
            eta_after_ = theDbe->book1D("ETA_LASTCUT",chtitle,50,-2.5,2.5);

            snprintf(chtitle, 255, "Muon transverse distance to beam spot [cm]");
            dxy_before_ = theDbe->book1D("DXY_BEFORECUTS",chtitle,100,-0.5,0.5);
            dxy_after_ = theDbe->book1D("DXY_LASTCUT",chtitle,100,-0.5,0.5);

            snprintf(chtitle, 255, "Quality-muon flag");
            goodewkmuon_before_ = theDbe->book1D("GOODEWKMUON_BEFORECUTS",chtitle,2,-0.5,1.5);
            goodewkmuon_after_ = theDbe->book1D("GOODEWKMUON_LASTCUT",chtitle,2,-0.5,1.5);

            if (isRelativeIso_) {
                  if (isCombinedIso_) {
                        snprintf(chtitle, 255, "Relative (combined) isolation variable");
                  } else {
                        snprintf(chtitle, 255, "Relative (tracker) isolation variable");
                  }
                  iso_before_ = theDbe->book1D("ISO_BEFORECUTS",chtitle,100, 0., 1.);
                  iso_after_ = theDbe->book1D("ISO_LASTCUT",chtitle,100, 0., 1.);
            } else {
                  if (isCombinedIso_) {
                        snprintf(chtitle, 255, "Absolute (combined) isolation variable [GeV]");
                  } else {
                        snprintf(chtitle, 255, "Absolute (tracker) isolation variable [GeV]");
                  }
                  iso_before_ = theDbe->book1D("ISO_BEFORECUTS",chtitle,100, 0., 20.);
                  iso_after_ = theDbe->book1D("ISO_LASTCUT",chtitle,100, 0., 20.);
            }

            snprintf(chtitle, 255, "HLT_Mu* Trigger response");
            trig_before_ = theDbe->book1D("TRIG_BEFORECUTS",chtitle,2,-0.5,1.5);
            trig_after_ = theDbe->book1D("TRIG_LASTCUT",chtitle,2,-0.5,1.5);

            snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
            mt_before_ = theDbe->book1D("MT_BEFORECUTS",chtitle,150,0.,300.);
            mt_after_ = theDbe->book1D("MT_LASTCUT",chtitle,150,0.,300.);

            snprintf(chtitle, 255, "Missing transverse energy (%s) [GeV]", metTag_.label().data());
            met_before_ = theDbe->book1D("MET_BEFORECUTS",chtitle,100,0.,200.);
            met_after_ = theDbe->book1D("MET_LASTCUT",chtitle,100,0.,200.);

            snprintf(chtitle, 255, "MU-MET (%s) acoplanarity", metTag_.label().data());
            acop_before_ = theDbe->book1D("ACOP_BEFORECUTS",chtitle,50,0.,M_PI);
            acop_after_ = theDbe->book1D("ACOP_LASTCUT",chtitle,50,0.,M_PI);

            snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ1_);
            nz1_before_ = theDbe->book1D("NZ1_BEFORECUTS",chtitle,10,-0.5,9.5);
            nz1_after_ = theDbe->book1D("NZ1_LASTCUT",chtitle,10,-0.5,9.5);

            snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ2_);
            nz2_before_ = theDbe->book1D("NZ2_BEFORECUTS",chtitle,10,-0.5,9.5);
            nz2_after_ = theDbe->book1D("NZ2_LASTCUT",chtitle,10,-0.5,9.5);

            snprintf(chtitle, 255, "Number of jets (%s) above %.2f GeV", jetTag_.label().data(), eJetMin_);
            njets_before_ = theDbe->book1D("NJETS_BEFORECUTS",chtitle,10,-0.5,9.5);
            njets_after_ = theDbe->book1D("NJETS_LASTCUT",chtitle,10,-0.5,9.5);

            snprintf(chtitle, 255, "DiMuonMass (2 globals)");
            dimuonmass_before_= theDbe->book1D("DIMUONMASS_BEFORECUTS",chtitle,100,0,200);
            dimuonmass_after_= theDbe->book1D("DIMUONMASS_AFTERZCUTS",chtitle,100,0,200);

            snprintf(chtitle, 255, "Global pt for Muons in Z");
            ptmuonZ_after_= theDbe->book1D("PT_AFTERZCUT",chtitle,100,0.,100.);

            snprintf(chtitle, 255, "Number of Valid Primary Vertices");
            npvs_before_ = theDbe->book1D("NPVs_BEFORECUTS",chtitle,10,-0.5,9.5);
            npvs_after_ = theDbe->book1D("NPVs_LASTCUT",chtitle,10,-0.5,9.5);

            snprintf(chtitle, 255, "Muon Charge");
            muoncharge_before_ = theDbe->book1D("MUONCHARGE_BEFORECUTS",chtitle,3,-1.5,1.5);
            muoncharge_after_ = theDbe->book1D("MUONCHARGE_LASTCUT",chtitle,3,-1.5,1.5);

      }
}

Member Data Documentation

MonitorElement* EwkMuDQM::acop_after_

private

Definition at line 114 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::acop_before_

private

Definition at line 113 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::acopCut_

private

Definition at line 47 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::chi2_after_

private

Definition at line 87 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::chi2_before_

private

Definition at line 86 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::dimuonmass_after_

private

Definition at line 126 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::dimuonmass_before_

private

Definition at line 125 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::dimuonSAmass_after_

private

Definition at line 129 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::dimuonSAmass_before_

private

Definition at line 128 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::dimuonSASAmass_after_

private

Definition at line 132 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::dimuonSASAmass_before_

private

Definition at line 131 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::dxy_after_

private

Definition at line 84 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::dxy_before_

private

Definition at line 83 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::dxyCut_

private

Definition at line 49 of file EwkMuDQM.h.

Referenced by analyze().

double EwkMuDQM::eJetMin_

private

Definition at line 60 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::eta_after_

private

Definition at line 81 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::eta_before_

private

Definition at line 80 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::etaCut_

private

Definition at line 39 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::goodewkmuon_after_

private

Definition at line 96 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::goodewkmuon_before_

private

Definition at line 95 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

HLTConfigProvider EwkMuDQM::hltConfigProvider_

private

Definition at line 64 of file EwkMuDQM.h.

Referenced by analyze(), and beginRun().

bool EwkMuDQM::isAlsoTrackerMuon_

private

Definition at line 54 of file EwkMuDQM.h.

bool EwkMuDQM::isCombinedIso_

private

Definition at line 41 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::iso_after_

private

Definition at line 102 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::iso_before_

private

Definition at line 101 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::isoCut03_

private

Definition at line 42 of file EwkMuDQM.h.

Referenced by analyze().

bool EwkMuDQM::isRelativeIso_

private

Definition at line 40 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

bool EwkMuDQM::isValidHltConfig_

private

Definition at line 63 of file EwkMuDQM.h.

Referenced by beginRun(), and EwkMuDQM().

edm::InputTag EwkMuDQM::jetTag_

private

Definition at line 35 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::met_after_

private

Definition at line 111 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::met_before_

private

Definition at line 110 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::metMax_

private

Definition at line 46 of file EwkMuDQM.h.

Referenced by analyze().

double EwkMuDQM::metMin_

private

Definition at line 45 of file EwkMuDQM.h.

Referenced by analyze().

edm::InputTag EwkMuDQM::metTag_

private

Definition at line 33 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::mt_after_

private

Definition at line 108 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::mt_before_

private

Definition at line 107 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::mtMax_

private

Definition at line 44 of file EwkMuDQM.h.

Referenced by analyze().

double EwkMuDQM::mtMin_

private

Definition at line 43 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::muoncharge_after_

private

Definition at line 138 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::muoncharge_before_

private

Definition at line 137 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::muonhits_after_

private

Definition at line 93 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::muonhits_before_

private

Definition at line 92 of file EwkMuDQM.h.

int EwkMuDQM::muonHitsCut_

private

Definition at line 53 of file EwkMuDQM.h.

Referenced by analyze().

edm::InputTag EwkMuDQM::muonTag_

private

Definition at line 32 of file EwkMuDQM.h.

Referenced by analyze().

unsigned int EwkMuDQM::nall

private

Definition at line 67 of file EwkMuDQM.h.

Referenced by analyze(), and beginRun().

MonitorElement* EwkMuDQM::nhits_after_

private

Definition at line 90 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::nhits_before_

private

Definition at line 89 of file EwkMuDQM.h.

unsigned int EwkMuDQM::nhlt

private

Definition at line 70 of file EwkMuDQM.h.

Referenced by beginRun().

unsigned int EwkMuDQM::niso

private

Definition at line 69 of file EwkMuDQM.h.

Referenced by beginRun().

int EwkMuDQM::nJetMax_

private

Definition at line 61 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::njets_after_

private

Definition at line 123 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::njets_before_

private

Definition at line 122 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

int EwkMuDQM::nMatchesCut_

private

Definition at line 55 of file EwkMuDQM.h.

Referenced by analyze().

unsigned int EwkMuDQM::nmet

private

Definition at line 71 of file EwkMuDQM.h.

Referenced by beginRun().

double EwkMuDQM::normalizedChi2Cut_

private

Definition at line 50 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::npvs_after_

private

Definition at line 135 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::npvs_before_

private

Definition at line 134 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

unsigned int EwkMuDQM::nrec

private

Definition at line 68 of file EwkMuDQM.h.

Referenced by beginRun().

unsigned int EwkMuDQM::nsel

private

Definition at line 72 of file EwkMuDQM.h.

Referenced by beginRun().

MonitorElement* EwkMuDQM::nz1_after_

private

Definition at line 117 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::nz1_before_

private

Definition at line 116 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::nz2_after_

private

Definition at line 120 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::nz2_before_

private

Definition at line 119 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

int EwkMuDQM::pixelHitsCut_

private

Definition at line 52 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::pt_after_

private

Definition at line 78 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::pt_before_

private

Definition at line 77 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::ptCut_

private

Definition at line 38 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::ptmuonZ_after_

private

Definition at line 140 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::ptThrForZ1_

private

Definition at line 57 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

double EwkMuDQM::ptThrForZ2_

private

Definition at line 58 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

DQMStore* EwkMuDQM::theDbe

private

Definition at line 75 of file EwkMuDQM.h.

Referenced by beginJob(), and init_histograms().

MonitorElement* EwkMuDQM::tkmu_after_

private

Definition at line 99 of file EwkMuDQM.h.

MonitorElement* EwkMuDQM::tkmu_before_

private

Definition at line 98 of file EwkMuDQM.h.

int EwkMuDQM::trackerHitsCut_

private

Definition at line 51 of file EwkMuDQM.h.

Referenced by analyze().

MonitorElement* EwkMuDQM::trig_after_

private

Definition at line 105 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

MonitorElement* EwkMuDQM::trig_before_

private

Definition at line 104 of file EwkMuDQM.h.

Referenced by analyze(), and init_histograms().

edm::InputTag EwkMuDQM::trigTag_

private

Definition at line 31 of file EwkMuDQM.h.

Referenced by analyze().

edm::InputTag EwkMuDQM::vertexTag_

private

Definition at line 36 of file EwkMuDQM.h.

Referenced by analyze().