CMS 3D CMS Logo

Public Member Functions | Protected Member Functions | Private Attributes

TopHLTDiMuonDQM Class Reference

#include <TopHLTDiMuonDQM.h>

Inheritance diagram for TopHLTDiMuonDQM:
edm::EDAnalyzer

List of all members.

Public Member Functions

 TopHLTDiMuonDQM (const edm::ParameterSet &)
 ~TopHLTDiMuonDQM ()

Protected Member Functions

void analyze (const edm::Event &, const edm::EventSetup &)
void beginJob ()
void beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)
void beginRun (const edm::Run &, const edm::EventSetup &)
void endJob ()
void endLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)
void endRun (const edm::Run &, const edm::EventSetup &)

Private Attributes

MonitorElementCombRelIso03
DQMStoredbe_
MonitorElementDeltaEtaMuonsRC
MonitorElementDeltaEtaMuonsWC
MonitorElementDeltaPhiMuonsRC
MonitorElementDeltaPhiMuonsWC
MonitorElementDeltaR_Match
MonitorElementDeltaR_Reco
MonitorElementDeltaR_Trig
MonitorElementDiMuonMassRC
MonitorElementDiMuonMassRC_LOGX
MonitorElementDiMuonMassWC
MonitorElementDiMuonMassWC_LOGX
MonitorElementEtaMuons
std::vector< std::string > hltPaths_L1_
std::vector< std::string > hltPaths_L3_
std::vector< std::string > hltPaths_sig_
std::vector< std::string > hltPaths_trig_
double MassWindow_down_
double MassWindow_up_
std::string monitorName_
double muon_eta_cut_
double muon_iso_cut_
double muon_pT_cut_
MonitorElementMuonEfficiency_eta
MonitorElementMuonEfficiency_eta_sig
MonitorElementMuonEfficiency_eta_trig
MonitorElementMuonEfficiency_phi
MonitorElementMuonEfficiency_phi_sig
MonitorElementMuonEfficiency_phi_trig
MonitorElementMuonEfficiency_pT
MonitorElementMuonEfficiency_pT_LOGX
MonitorElementMuonEfficiency_pT_LOGX_sig
MonitorElementMuonEfficiency_pT_LOGX_trig
MonitorElementMuonEfficiency_pT_sig
MonitorElementMuonEfficiency_pT_trig
edm::InputTag muons_
MonitorElementNMuons
MonitorElementNMuons_charge
MonitorElementNMuons_iso
MonitorElementNTracks
MonitorElementPhiMuons
MonitorElementPixelHits_muons
MonitorElementPtMuons
MonitorElementPtMuons_LOGX
MonitorElementTrackerHits_muons
MonitorElementTrigger_Match
MonitorElementTriggerEfficiencies
MonitorElementTriggerEfficiencies_sig
MonitorElementTriggerEfficiencies_trig
edm::InputTag triggerEvent_
edm::InputTag triggerFilter_
edm::InputTag triggerResults_
MonitorElementTrigs
edm::InputTag vertex_
double vertex_X_cut_
double vertex_Y_cut_
double vertex_Z_cut_
MonitorElementVxVy_muons
MonitorElementVz_muons

Detailed Description

Definition at line 51 of file TopHLTDiMuonDQM.h.


Constructor & Destructor Documentation

TopHLTDiMuonDQM::TopHLTDiMuonDQM ( const edm::ParameterSet ps)

Definition at line 15 of file TopHLTDiMuonDQM.cc.

References edm::ParameterSet::getParameter().

                                                            {

  monitorName_ = ps.getParameter<string>("monitorName");

  triggerResults_ = ps.getParameter<edm::InputTag>("TriggerResults");
  triggerEvent_   = ps.getParameter<edm::InputTag>("TriggerEvent");
  triggerFilter_  = ps.getParameter<edm::InputTag>("TriggerFilter");

  hltPaths_L1_   = ps.getParameter<vector<string> >("hltPaths_L1");
  hltPaths_L3_   = ps.getParameter<vector<string> >("hltPaths_L3");
  hltPaths_sig_  = ps.getParameter<vector<string> >("hltPaths_sig");
  hltPaths_trig_ = ps.getParameter<vector<string> >("hltPaths_trig");

  vertex_       = ps.getParameter<edm::InputTag>("vertexCollection");
  vertex_X_cut_ = ps.getParameter<double>("vertex_X_cut");
  vertex_Y_cut_ = ps.getParameter<double>("vertex_Y_cut");
  vertex_Z_cut_ = ps.getParameter<double>("vertex_Z_cut");

  muons_        = ps.getParameter<edm::InputTag>("muonCollection");
  muon_pT_cut_  = ps.getParameter<double>("muon_pT_cut");
  muon_eta_cut_ = ps.getParameter<double>("muon_eta_cut");
  muon_iso_cut_ = ps.getParameter<double>("muon_iso_cut");

  MassWindow_up_   = ps.getParameter<double>("MassWindow_up");
  MassWindow_down_ = ps.getParameter<double>("MassWindow_down");

}
TopHLTDiMuonDQM::~TopHLTDiMuonDQM ( )

Definition at line 44 of file TopHLTDiMuonDQM.cc.

                                  {

}

Member Function Documentation

void TopHLTDiMuonDQM::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
) [protected, virtual]

Implements edm::EDAnalyzer.

Definition at line 180 of file TopHLTDiMuonDQM.cc.

References abs, SiPixelRawToDigiRegional_cfi::deltaPhi, deltaR(), reco::MuonIsolation::emEt, reco::Particle::eta(), reco::LeafCandidate::eta(), trigger::TriggerObject::eta(), edm::HandleBase::failedToGet(), edm::Event::getByLabel(), reco::MuonIsolation::hadEt, CreateSqliteForCondDB_cfg::hltPaths, reco::MuonIsolation::hoEt, i, reco::Vertex::isFake(), j, gen::k, relativeConstraints::keys, metsig::muon, patZpeak::muons, trigger::TriggerObject::particle(), reco::Particle::pdgId(), reco::LeafCandidate::phi(), reco::Particle::phi(), trigger::TriggerObject::pt(), mathSSE::sqrt(), reco::MuonIsolation::sumPt, reco::Vertex::tracksSize(), edm::TriggerNames::triggerName(), edm::Event::triggerNames(), reco::Vertex::x(), reco::Vertex::y(), and reco::Vertex::z().

                                                                                  {

  // ------------------------
  //  Global Event Variables
  // ------------------------

  vector<string> hltPaths = hltPaths_L3_;

  vector<reco::Particle> Triggered_muons;
  reco::MuonCollection    Isolated_muons;
  reco::MuonCollection     Matched_muons;

  const int N_TriggerPaths = hltPaths.size();
  const int N_SignalPaths  = hltPaths_sig_.size();
  const int N_ControlPaths = hltPaths_trig_.size();

  bool Fired_Signal_Trigger[ 10] = {false};
  bool Fired_Control_Trigger[10] = {false};

  double DilepMass = 0.;

  double deltaR_Trig  = 1000.;
  double deltaR_Reco  =    0.;
  double deltaR_Match =    0.;

  int N_iso_mu = 0;

  double vertex_X = 100.;
  double vertex_Y = 100.;
  double vertex_Z = 100.;

  // -------------------------
  //  Analyze Trigger Results
  // -------------------------

  edm::Handle<TriggerResults> trigResults;
  iEvent.getByLabel(triggerResults_, trigResults);

  if( trigResults.failedToGet() ) {

    //    cout << endl << "-----------------------------" << endl;
    //    cout << "--- NO TRIGGER RESULTS !! ---" << endl;
    //    cout << "-----------------------------" << endl << endl;

  }

  if( !trigResults.failedToGet() ) {

    const edm::TriggerNames & trigName = iEvent.triggerNames(*trigResults);

    for( unsigned int i_Trig = 0; i_Trig < trigResults->size(); ++i_Trig ) {

      if(trigResults->accept(i_Trig)) {

        // Check for all trigger paths

        for( int i = 0; i < N_TriggerPaths; i++ ) {

          if( trigName.triggerName(i_Trig) == hltPaths[i] )  Trigs->Fill(i);

        }

        // Check for signal & control trigger paths

        for( int j = 0; j < N_SignalPaths; ++j ) {

          if( trigName.triggerName(i_Trig) == hltPaths_sig_[j]  )  Fired_Signal_Trigger[j]  = true;

        }

        for( int k = 0; k < N_ControlPaths; ++k ) {

          if( trigName.triggerName(i_Trig) == hltPaths_trig_[k] )  Fired_Control_Trigger[k] = true;

        }

      }

    }

  }

  // -----------------------
  //  Analyze Trigger Event
  // -----------------------

  edm::Handle<TriggerEvent> triggerEvent;
  iEvent.getByLabel(triggerEvent_, triggerEvent);

  if( triggerEvent.failedToGet() ) {

    //    cout << endl << "---------------------------" << endl;
    //    cout << "--- NO TRIGGER EVENT !! ---" << endl;
    //    cout << "---------------------------" << endl << endl;

  }

  if( !triggerEvent.failedToGet() ) {

    size_t filterIndex = triggerEvent->filterIndex( triggerFilter_ );
    TriggerObjectCollection triggerObjects = triggerEvent->getObjects();

    if( filterIndex >= triggerEvent->sizeFilters() ) {

      //      cout << endl << "------------------------------" << endl;
      //      cout << "--- NO FILTERED OBJECTS !! ---" << endl;
      //      cout << "------------------------------" << endl << endl;

    }

    if( filterIndex < triggerEvent->sizeFilters() ) {

      const Keys & keys = triggerEvent->filterKeys( filterIndex );

      int N_mu = 0;

      for( size_t j = 0; j < keys.size(); j++ ) {

        TriggerObject foundObject = triggerObjects[keys[j]];

        if(      foundObject.pt()   < muon_pT_cut_  )  continue;
        if( abs( foundObject.eta()) > muon_eta_cut_ )  continue;
        if( abs( foundObject.particle().pdgId() ) != 13 )  continue;

        ++N_mu;
        Triggered_muons.push_back( foundObject.particle() );

      }

      if( Triggered_muons.size() == 2 ) {

        reco::Particle mu1 = Triggered_muons.at(0);
        reco::Particle mu2 = Triggered_muons.at(1);

        deltaR_Trig = deltaR( mu1.eta(), mu1.phi(), mu2.eta(), mu2.phi() );
        DeltaR_Trig->Fill(deltaR_Trig);

      }

    }

  }

  // ------------------------
  //  Analyze Primary Vertex
  // ------------------------

  edm::Handle<reco::VertexCollection> vertexs;
  iEvent.getByLabel(vertex_, vertexs);

  if( vertexs.failedToGet() ) {

    //      cout << endl << "----------------------------" << endl;
    //      cout << "--- NO PRIMARY VERTEX !! ---" << endl;
    //      cout << "----------------------------" << endl << endl;

  }

  if( !vertexs.failedToGet() ) {

    reco::Vertex primaryVertex = vertexs->front();

    int numberTracks = primaryVertex.tracksSize();
    //      double ndof      = primaryVertex.ndof();
    bool fake        = primaryVertex.isFake();

    NTracks->Fill(numberTracks);

    if( !fake && numberTracks > 3 ) {

      vertex_X = primaryVertex.x();
      vertex_Y = primaryVertex.y();
      vertex_Z = primaryVertex.z();

    }

  }

  // --------------------
  //  Analyze RECO Muons
  // --------------------

  edm::Handle<reco::MuonCollection> muons;
  iEvent.getByLabel(muons_, muons);

  reco::MuonCollection::const_iterator muon;

  if( muons.failedToGet() ) {

    //    cout << endl << "------------------------" << endl;
    //    cout << "--- NO RECO MUONS !! ---" << endl;
    //    cout << "------------------------" << endl << endl;

  }

  if( !muons.failedToGet() ) {

    NMuons->Fill( muons->size() );

    //    cout << "N_muons : " << muons->size() << endl;

    for(muon = muons->begin(); muon!= muons->end(); ++muon) {

      float N_muons = muons->size();
      float Q_muon  = muon->charge();

      NMuons_charge->Fill(N_muons*Q_muon);

      double track_X = 100.;
      double track_Y = 100.;
      double track_Z = 100.;

      double N_PixelHits   = 0.;
      double N_TrackerHits = 0.;

      if( muon->isGlobalMuon() && muon->isTrackerMuon() ) {

        reco::TrackRef track = muon->globalTrack();

        track_X       = track->vx();
        track_Y       = track->vy();
        track_Z       = track->vz();
        N_PixelHits   = track->hitPattern().numberOfValidPixelHits();
        N_TrackerHits = track->hitPattern().numberOfValidTrackerHits();

        VxVy_muons->Fill(track_X, track_Y);
        Vz_muons->Fill(track_Z);
        PixelHits_muons->Fill(N_PixelHits);
        TrackerHits_muons->Fill(N_TrackerHits);

      }

      // Vertex and kinematic cuts

      if(          track_X > vertex_X_cut_ )  continue;
      if(          track_Y > vertex_Y_cut_ )  continue;
      if(          track_Z > vertex_Z_cut_ )  continue;
      if(      N_PixelHits <  1.           )  continue;
      if(    N_TrackerHits < 11.           )  continue;
      if(     muon->pt()   < muon_pT_cut_  )  continue;
      if( abs(muon->eta()) > muon_eta_cut_ )  continue;

      reco::MuonIsolation muIso03 = muon->isolationR03();

      double muonCombRelIso = 1.;

      if ( muon->pt() != 0. )
        muonCombRelIso = ( muIso03.emEt + muIso03.hadEt + muIso03.hoEt + muIso03.sumPt ) / muon->pt();

      CombRelIso03->Fill( muonCombRelIso );

      if( muonCombRelIso < muon_iso_cut_ ) {

        ++N_iso_mu;
        Isolated_muons.push_back(*muon);

      }

    }

    NMuons_iso->Fill(N_iso_mu);

    //    if( Isolated_muons.size() > 1 && Fired_Control_Trigger[0] ) {
    if( Isolated_muons.size() > 1 ) {

      // Vertex cut

      if( vertex_X < vertex_X_cut_ && vertex_Y < vertex_Y_cut_ && vertex_Z < vertex_Z_cut_ ) {

        for( int i = 0; i < (static_cast<int>(Isolated_muons.size()) - 1); ++i ) {

          for( int j = i+1; j < static_cast<int>(Isolated_muons.size()); ++j ) {

            reco::MuonCollection::const_reference mu1 = Isolated_muons.at(i);
            reco::MuonCollection::const_reference mu2 = Isolated_muons.at(j);

            DilepMass = sqrt( (mu1.energy()+mu2.energy())*(mu1.energy()+mu2.energy())
                              - (mu1.px()+mu2.px())*(mu1.px()+mu2.px())
                              - (mu1.py()+mu2.py())*(mu1.py()+mu2.py())
                              - (mu1.pz()+mu2.pz())*(mu1.pz()+mu2.pz())
                              );

            if( DilepMass < 1. ) {

              if( i > 0 ) {

                Isolated_muons.erase(Isolated_muons.begin()+i);
                --i;

              }

              continue;

            }

            // Opposite muon charges -> Right Charge (RC)

            if( mu1.charge()*mu2.charge() < 0. ) {

              DiMuonMassRC->Fill( DilepMass );
              DiMuonMassRC_LOGX->Fill( DilepMass );

              if( DilepMass > MassWindow_down_ && DilepMass < MassWindow_up_ ) {

                PtMuons->Fill(  mu1.pt()  );
                PtMuons->Fill(  mu2.pt()  );
                PtMuons_LOGX->Fill( mu1.pt() );
                PtMuons_LOGX->Fill( mu2.pt() );
                EtaMuons->Fill( mu1.eta() );
                EtaMuons->Fill( mu2.eta() );
                PhiMuons->Fill( mu1.phi() );
                PhiMuons->Fill( mu2.phi() );

                DeltaEtaMuonsRC->Fill(mu1.eta()-mu2.eta());
                DeltaPhiMuonsRC->Fill( deltaPhi(mu1.phi(),mu2.phi()) );

                // Determinating relative trigger efficiencies

                for( int k = 0; k < N_SignalPaths; ++k ) {

                  if( Fired_Signal_Trigger[k] && Fired_Control_Trigger[k] )  TriggerEfficiencies_sig->Fill(k);

                  if( Fired_Control_Trigger[k] )  TriggerEfficiencies_trig->Fill(k);

                }

                // Trigger object matching

                int N_Match = 0;
                double   DR = 0.1;

                if( Isolated_muons.size() == 2 && Triggered_muons.size() > 0 ) {

                  deltaR_Reco = deltaR( mu1.eta(), mu1.phi(), mu2.eta(), mu2.phi() );
                  DeltaR_Reco->Fill(deltaR_Reco);

                  if( deltaR_Reco > 2.*DR && deltaR_Trig > 2.*DR ) {

                    for( int i = 0; i < static_cast<int>(Isolated_muons.size()); ++i ) {

                      for( int j = 0; j < static_cast<int>(Triggered_muons.size()); ++j ) {

                        reco::Particle & Trigger_mu = Triggered_muons.at(j);
                        reco::Muon     &    Reco_mu =  Isolated_muons.at(i);

                        deltaR_Match = deltaR( Trigger_mu.eta(), Trigger_mu.phi(), Reco_mu.eta(), Reco_mu.phi() );
                        DeltaR_Match->Fill(deltaR_Match);

                        if( deltaR_Match < DR) {

                          ++N_Match;
                          Matched_muons.push_back(Reco_mu);

                        }

                      }

                    }

                    Trigger_Match->Fill(N_Match);

                  }

                }


                // Muon Tag & Probe Efficiency

                if( Matched_muons.size() == 1 ) {

                  reco::MuonCollection::const_reference matched_mu1 = Matched_muons.at(0);

                  MuonEfficiency_pT_trig->Fill( matched_mu1.pt() );
                  MuonEfficiency_pT_LOGX_trig->Fill( matched_mu1.pt() );
                  MuonEfficiency_eta_trig->Fill(matched_mu1.eta());
                  MuonEfficiency_phi_trig->Fill(matched_mu1.phi());

                }

                if( Matched_muons.size() == 2 ) {

                  reco::MuonCollection::const_reference matched_mu1 = Matched_muons.at(0);
                  reco::MuonCollection::const_reference matched_mu2 = Matched_muons.at(1);

                  MuonEfficiency_pT_trig->Fill( matched_mu1.pt() );
                  MuonEfficiency_pT_trig->Fill( matched_mu2.pt() );
                  MuonEfficiency_pT_LOGX_trig->Fill( matched_mu1.pt() );
                  MuonEfficiency_pT_LOGX_trig->Fill( matched_mu2.pt() );
                  MuonEfficiency_eta_trig->Fill(matched_mu1.eta());
                  MuonEfficiency_eta_trig->Fill(matched_mu2.eta());
                  MuonEfficiency_phi_trig->Fill(matched_mu1.phi());
                  MuonEfficiency_phi_trig->Fill(matched_mu2.phi());

                  MuonEfficiency_pT_sig->Fill( matched_mu1.pt() );
                  MuonEfficiency_pT_sig->Fill( matched_mu2.pt() );
                  MuonEfficiency_pT_LOGX_sig->Fill( matched_mu1.pt() );
                  MuonEfficiency_pT_LOGX_sig->Fill( matched_mu2.pt() );
                  MuonEfficiency_eta_sig->Fill(matched_mu1.eta());
                  MuonEfficiency_eta_sig->Fill(matched_mu2.eta());
                  MuonEfficiency_phi_sig->Fill(matched_mu1.phi());
                  MuonEfficiency_phi_sig->Fill(matched_mu2.phi());

                }

              }

            }

            // Same muon charges -> Wrong Charge (WC)

            if( mu1.charge()*mu2.charge() > 0. ) {

              DiMuonMassWC->Fill( DilepMass );
              DiMuonMassWC_LOGX->Fill( DilepMass );

              if( DilepMass > MassWindow_down_ && DilepMass < MassWindow_up_ ) {

                DeltaEtaMuonsWC->Fill( mu1.eta()-mu2.eta() );
                DeltaPhiMuonsWC->Fill( deltaPhi(mu1.phi(),mu2.phi()) );

              }

            }

          }

        }

      }

    }

  }

}
void TopHLTDiMuonDQM::beginJob ( void  ) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 49 of file TopHLTDiMuonDQM.cc.

References DQMStore::book1D(), DQMStore::book2D(), dbe_, i, j, label, funct::log(), cmsCodeRules::cppFunctionSkipper::operator, funct::pow(), DQMStore::setCurrentFolder(), and MonitorElement::setTitle().

                               {

  dbe_ = Service<DQMStore>().operator->();

  if( dbe_ ) {

    dbe_->setCurrentFolder(monitorName_);

    Trigs = dbe_->book1D("Trigs", "Fired triggers", 15, 0., 15.);

    TriggerEfficiencies = dbe_->book1D("TriggerEfficiencies", "HL Trigger Efficiencies", 10, 0., 10.);

    TriggerEfficiencies->setTitle("HL Trigger Efficiencies #epsilon_{signal} = #frac{[signal] && [control]}{[control]}");

    TriggerEfficiencies_sig  = dbe_->book1D("TriggerEfficiencies_sig",  "HL Trigger Signal && Control Counts", 10, 0., 10.);
    TriggerEfficiencies_trig = dbe_->book1D("TriggerEfficiencies_trig", "HL Trigger Control Counts",           10, 0., 10.);

    const int nbins_Pt = 5;

    float bins_Pt[nbins_Pt+1] = { 0., 5., 10., 20., 50., 100. };

    MuonEfficiency_pT      = dbe_->book1D("MuonEfficiency_pT",      "Muon Efficiency P_{T}",           nbins_Pt, &bins_Pt[0]);
    MuonEfficiency_pT_sig  = dbe_->book1D("MuonEfficiency_pT_sig",  "P^{#mu}_{T} (signal triggered)",  nbins_Pt, &bins_Pt[0]);
    MuonEfficiency_pT_trig = dbe_->book1D("MuonEfficiency_pT_trig", "P^{#mu}_{T} (control triggered)", nbins_Pt, &bins_Pt[0]);

    const int nbins_eta = 7;

    float bins_eta[nbins_eta+1] = { -3.5, -2.5, -1.5, -0.5, 0.5, 1.5, 2.5, 3.5 };

    MuonEfficiency_eta      = dbe_->book1D("MuonEfficiency_eta",      "Muon Efficiency  #eta",           nbins_eta, &bins_eta[0]);
    MuonEfficiency_eta_sig  = dbe_->book1D("MuonEfficiency_eta_sig",  "#eta_{muon} (signal triggered)",  nbins_eta, &bins_eta[0]);
    MuonEfficiency_eta_trig = dbe_->book1D("MuonEfficiency_eta_trig", "#eta_{muon} (control triggered)", nbins_eta, &bins_eta[0]);

    const int nbins_phi = 9;

    float bins_phi[nbins_phi+1] = { -3.5, -3.2, -2.6, -1.56, -0.52, 0.52, 1.56, 2.6, 3.2, 3.5 };

    MuonEfficiency_phi      = dbe_->book1D("MuonEfficiency_phi",      "Muon Efficiency  #phi",           nbins_phi, &bins_phi[0]);
    MuonEfficiency_phi_sig  = dbe_->book1D("MuonEfficiency_phi_sig",  "#phi_{muon} (signal triggered)",  nbins_phi, &bins_phi[0]);
    MuonEfficiency_phi_trig = dbe_->book1D("MuonEfficiency_phi_trig", "#phi_{muon} (control triggered)", nbins_phi, &bins_phi[0]);

    const int N_TriggerPaths = hltPaths_L3_.size();
    const int N_SignalPaths  = hltPaths_sig_.size();

    for( int i = 0; i < N_TriggerPaths; i++ ) {

      const string &label = hltPaths_L3_[i];
      Trigs->setBinLabel( i+1,label.c_str() );

    }

    for( int j = 0; j < N_SignalPaths; ++j ) {

      const string &label_eff = "#frac{["+hltPaths_sig_[j]+"]}{vs. ["+hltPaths_trig_[j]+"]}";
      const string &label_sig = hltPaths_sig_[j]+"\n && "+hltPaths_trig_[j];
      TriggerEfficiencies->setBinLabel(     j+1, label_eff.c_str() );
      TriggerEfficiencies_sig->setBinLabel( j+1, label_sig.c_str() );

    }

    NMuons        = dbe_->book1D("Nmuons",        "Number of muons",             20,   0.,  10.);
    NMuons_iso    = dbe_->book1D("Nmuons_Iso",    "Number of isolated muons",    20,   0.,  10.);
    NMuons_charge = dbe_->book1D("Nmuons_Charge", "N_{muons} * Q(#mu)",          19, -10.,  10.);
    NTracks       = dbe_->book1D("Ntracks",       "Number of tracks",            50,   0.,  50.);
    VxVy_muons    = dbe_->book2D("VxVy_muons",    "Vertex x-y-positon (global)", 40,  -1.,   1., 40 , -1., 1.);
    Vz_muons      = dbe_->book1D("Vz_muons",      "Vertex z-positon (global)",   40, -20.,  20.);
    PtMuons       = dbe_->book1D("PtMuon",        "P^{#mu}_{T}",                 50,   0., 100.);
    EtaMuons      = dbe_->book1D("EtaMuon",       "#eta_{muon}",                 50,  -5.,   5.);
    PhiMuons      = dbe_->book1D("PhiMuon",       "#phi_{muon}",                 40,  -4.,   4.);

    DeltaEtaMuonsRC   = dbe_->book1D("DeltaEtaMuonsRC",    "#Delta #eta of muon pair (RC)", 50,  -5.,   5.);
    DeltaPhiMuonsRC   = dbe_->book1D("DeltaPhiMuonsRC",    "#Delta #phi of muon pair (RC)", 50,  -5.,   5.);
    DeltaEtaMuonsWC   = dbe_->book1D("DeltaEtaMuonsWC",    "#Delta #eta of muon pair (WC)", 50,  -5.,   5.);
    DeltaPhiMuonsWC   = dbe_->book1D("DeltaPhiMuonsWC",    "#Delta #phi of muon pair (WC)", 50,  -5.,   5.);
    CombRelIso03      = dbe_->book1D("MuIso_CombRelIso03", "Muon CombRelIso dR=03",         50,   0.,   1.);
    PixelHits_muons   = dbe_->book1D("NPixelHits_muons",   "Number of pixel hits",          50,   0.,  50.);
    TrackerHits_muons = dbe_->book1D("NTrackerHits_muons", "Number of hits in the tracker", 50,   0.,  50.);

    DeltaR_Trig   = dbe_->book1D("DeltaRTrigger", "#Delta R of trigger muon pair",       50, 0., 5.);
    DeltaR_Reco   = dbe_->book1D("DeltaRReco",    "#Delta R of RECO muon pair",          50, 0., 5.);
    DeltaR_Match  = dbe_->book1D("DeltaRMatch",   "#Delta R of matched muon pair",       50, 0., 5.);
    Trigger_Match = dbe_->book1D("TriggerMatch",  "Number of Trigger-RECO assignements",  6, 0., 6.);

    const int nbins_Pt_Log = 15;

    double logmin = 0.;
    double logmax = 3.;  // 10^(3.)=1000

    float bins_Pt_Log[nbins_Pt_Log+1];

    for (int i = 0; i <= nbins_Pt_Log; i++) {
      double log = logmin + (logmax-logmin)*i/nbins_Pt_Log;
      bins_Pt_Log[i] = std::pow(10.0, log);
    }

    PtMuons_LOGX  = dbe_->book1D("Pt_muon_LOGX", "P^{#mu}_{T}", nbins_Pt_Log, &bins_Pt_Log[0]);

    MuonEfficiency_pT_LOGX      = dbe_->book1D("MuonEfficiency_pT_LOGX",      "Muon Efficiency P_{T}",           nbins_Pt_Log, &bins_Pt_Log[0]);
    MuonEfficiency_pT_LOGX_sig  = dbe_->book1D("MuonEfficiency_pT_LOGX_sig",  "P^{#mu}_{T} (signal triggered)",  nbins_Pt_Log, &bins_Pt_Log[0]);
    MuonEfficiency_pT_LOGX_trig = dbe_->book1D("MuonEfficiency_pT_LOGX_trig", "P^{#mu}_{T} (control triggered)", nbins_Pt_Log, &bins_Pt_Log[0]);

    const int nbins_mass = 200;

    float bins_mass[nbins_mass+1];

    for (int i = 0; i <= nbins_mass; i++) {
      double log = logmin + (logmax-logmin)*i/nbins_mass;
      bins_mass[i] = std::pow(10.0, log);
    }

    DiMuonMassRC       = dbe_->book1D("DiMuonMassRC",      "Invariant Dimuon Mass (Right Charge)", 50, 0., 200.);
    DiMuonMassWC       = dbe_->book1D("DiMuonMassWC",      "Invariant Dimuon Mass (Wrong Charge)", 50, 0., 200.);

    DiMuonMassRC_LOGX  = dbe_->book1D("DiMuonMassRC_LOGX", "Invariant Dimuon Mass (Right Charge)", nbins_mass, &bins_mass[0]);
    DiMuonMassWC_LOGX  = dbe_->book1D("DiMuonMassWC_LOGX", "Invariant Dimuon Mass (Wrong Charge)", nbins_mass, &bins_mass[0]);

  }

}
void TopHLTDiMuonDQM::beginLuminosityBlock ( const edm::LuminosityBlock lumiSeg,
const edm::EventSetup context 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 175 of file TopHLTDiMuonDQM.cc.

                                                                                                          {

}
void TopHLTDiMuonDQM::beginRun ( const edm::Run r,
const edm::EventSetup context 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 170 of file TopHLTDiMuonDQM.cc.

                                                                            {

}
void TopHLTDiMuonDQM::endJob ( void  ) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 627 of file TopHLTDiMuonDQM.cc.

                             {

}
void TopHLTDiMuonDQM::endLuminosityBlock ( const edm::LuminosityBlock lumiSeg,
const edm::EventSetup context 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 617 of file TopHLTDiMuonDQM.cc.

                                                                                                        {

}
void TopHLTDiMuonDQM::endRun ( const edm::Run r,
const edm::EventSetup context 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 622 of file TopHLTDiMuonDQM.cc.

                                                                          {

}

Member Data Documentation

Definition at line 107 of file TopHLTDiMuonDQM.h.

Definition at line 72 of file TopHLTDiMuonDQM.h.

Definition at line 138 of file TopHLTDiMuonDQM.h.

Definition at line 140 of file TopHLTDiMuonDQM.h.

Definition at line 139 of file TopHLTDiMuonDQM.h.

Definition at line 141 of file TopHLTDiMuonDQM.h.

Definition at line 145 of file TopHLTDiMuonDQM.h.

Definition at line 144 of file TopHLTDiMuonDQM.h.

Definition at line 143 of file TopHLTDiMuonDQM.h.

Definition at line 133 of file TopHLTDiMuonDQM.h.

Definition at line 135 of file TopHLTDiMuonDQM.h.

Definition at line 134 of file TopHLTDiMuonDQM.h.

Definition at line 136 of file TopHLTDiMuonDQM.h.

Definition at line 105 of file TopHLTDiMuonDQM.h.

std::vector<std::string> TopHLTDiMuonDQM::hltPaths_L1_ [private]

Definition at line 82 of file TopHLTDiMuonDQM.h.

std::vector<std::string> TopHLTDiMuonDQM::hltPaths_L3_ [private]

Definition at line 83 of file TopHLTDiMuonDQM.h.

std::vector<std::string> TopHLTDiMuonDQM::hltPaths_sig_ [private]

Definition at line 84 of file TopHLTDiMuonDQM.h.

std::vector<std::string> TopHLTDiMuonDQM::hltPaths_trig_ [private]

Definition at line 85 of file TopHLTDiMuonDQM.h.

Definition at line 96 of file TopHLTDiMuonDQM.h.

Definition at line 95 of file TopHLTDiMuonDQM.h.

std::string TopHLTDiMuonDQM::monitorName_ [private]

Definition at line 73 of file TopHLTDiMuonDQM.h.

Definition at line 92 of file TopHLTDiMuonDQM.h.

Definition at line 93 of file TopHLTDiMuonDQM.h.

Definition at line 91 of file TopHLTDiMuonDQM.h.

Definition at line 125 of file TopHLTDiMuonDQM.h.

Definition at line 126 of file TopHLTDiMuonDQM.h.

Definition at line 127 of file TopHLTDiMuonDQM.h.

Definition at line 129 of file TopHLTDiMuonDQM.h.

Definition at line 130 of file TopHLTDiMuonDQM.h.

Definition at line 131 of file TopHLTDiMuonDQM.h.

Definition at line 117 of file TopHLTDiMuonDQM.h.

Definition at line 121 of file TopHLTDiMuonDQM.h.

Definition at line 122 of file TopHLTDiMuonDQM.h.

Definition at line 123 of file TopHLTDiMuonDQM.h.

Definition at line 118 of file TopHLTDiMuonDQM.h.

Definition at line 119 of file TopHLTDiMuonDQM.h.

Definition at line 80 of file TopHLTDiMuonDQM.h.

Definition at line 100 of file TopHLTDiMuonDQM.h.

Definition at line 101 of file TopHLTDiMuonDQM.h.

Definition at line 102 of file TopHLTDiMuonDQM.h.

Definition at line 99 of file TopHLTDiMuonDQM.h.

Definition at line 106 of file TopHLTDiMuonDQM.h.

Definition at line 110 of file TopHLTDiMuonDQM.h.

Definition at line 103 of file TopHLTDiMuonDQM.h.

Definition at line 104 of file TopHLTDiMuonDQM.h.

Definition at line 111 of file TopHLTDiMuonDQM.h.

Definition at line 146 of file TopHLTDiMuonDQM.h.

Definition at line 113 of file TopHLTDiMuonDQM.h.

Definition at line 114 of file TopHLTDiMuonDQM.h.

Definition at line 115 of file TopHLTDiMuonDQM.h.

Definition at line 76 of file TopHLTDiMuonDQM.h.

Definition at line 77 of file TopHLTDiMuonDQM.h.

Definition at line 75 of file TopHLTDiMuonDQM.h.

Definition at line 98 of file TopHLTDiMuonDQM.h.

Definition at line 79 of file TopHLTDiMuonDQM.h.

Definition at line 87 of file TopHLTDiMuonDQM.h.

Definition at line 88 of file TopHLTDiMuonDQM.h.

Definition at line 89 of file TopHLTDiMuonDQM.h.

Definition at line 108 of file TopHLTDiMuonDQM.h.

Definition at line 109 of file TopHLTDiMuonDQM.h.