EwkDQM Class Reference

#include <EwkDQM.h>

Inheritance diagram for EwkDQM:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Get the analysis. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
double	calcDeltaPhi (double phi1, double phi2)
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	EwkDQM (const edm::ParameterSet &)
	Constructor. More...
	~EwkDQM () override
	Destructor. More...
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
double	eJetMin_
MonitorElement *	h_e1_et
MonitorElement *	h_e1_eta
MonitorElement *	h_e1_phi
MonitorElement *	h_e2_et
MonitorElement *	h_e2_eta
MonitorElement *	h_e2_phi
MonitorElement *	h_e_invWMass
MonitorElement *	h_ee_invMass
MonitorElement *	h_jet2_et
MonitorElement *	h_jet2_eta
MonitorElement *	h_jet2_phi
MonitorElement *	h_jet_count
MonitorElement *	h_jet_et
MonitorElement *	h_jet_eta
MonitorElement *	h_jet_phi
MonitorElement *	h_jet_pt
MonitorElement *	h_m1_eta
MonitorElement *	h_m1_phi
MonitorElement *	h_m1_pt
MonitorElement *	h_m2_eta
MonitorElement *	h_m2_phi
MonitorElement *	h_m2_pt
MonitorElement *	h_m_invWMass
MonitorElement *	h_met
MonitorElement *	h_met_phi
MonitorElement *	h_mumu_invMass
MonitorElement *	h_vertex_chi2
MonitorElement *	h_vertex_d0
MonitorElement *	h_vertex_number
MonitorElement *	h_vertex_numTrks
MonitorElement *	h_vertex_sumTrks
HLTConfigProvider	hltConfigProvider_
bool	isValidHltConfig_
std::string	logTraceName
edm::InputTag	theCaloMETCollectionLabel_
edm::EDGetTokenT< edm::View< reco::MET > >	theCaloMETCollectionToken_
std::vector< std::string >	theElecTriggerPathToPass_
edm::EDGetTokenT< reco::GsfElectronCollection >	theElectronCollectionLabel_
edm::EDGetTokenT< reco::MuonCollection >	theMuonCollectionLabel_
std::vector< std::string >	theMuonTriggerPathToPass_
edm::InputTag	thePFJetCollectionLabel_
edm::EDGetTokenT< edm::View< reco::Jet > >	thePFJetCollectionToken_
edm::InputTag	theTriggerResultsCollection_
edm::EDGetTokenT< edm::TriggerResults >	theTriggerResultsToken_
edm::EDGetTokenT< reco::VertexCollection >	theVertexToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
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
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

DQM offline for SMP V+Jets

Author

Valentina Gori, University of Firenze

Definition at line 29 of file EwkDQM.h.

Constructor & Destructor Documentation

EwkDQM()

EwkDQM::EwkDQM

(

const edm::ParameterSet &

parameters

)

Constructor.

Definition at line 43 of file EwkDQM.cc.

                                             {
  eJetMin_ = parameters.getUntrackedParameter<double>("EJetMin", 999999.);
 
  // riguardare questa sintassi
  // Get parameters from configuration file
  thePFJetCollectionLabel_ = parameters.getParameter<InputTag>("PFJetCollection");
  theCaloMETCollectionLabel_ = parameters.getParameter<InputTag>("caloMETCollection");
  theTriggerResultsCollection_ = parameters.getParameter<InputTag>("triggerResultsCollection");
 
  theElecTriggerPathToPass_ = parameters.getParameter<std::vector<string> >("elecTriggerPathToPass");
  theMuonTriggerPathToPass_ = parameters.getParameter<std::vector<string> >("muonTriggerPathToPass");
  theTriggerResultsToken_ =
      consumes<edm::TriggerResults>(parameters.getParameter<InputTag>("triggerResultsCollection"));
  theMuonCollectionLabel_ = consumes<reco::MuonCollection>(parameters.getParameter<InputTag>("muonCollection"));
  theElectronCollectionLabel_ =
      consumes<reco::GsfElectronCollection>(parameters.getParameter<InputTag>("electronCollection"));
  thePFJetCollectionToken_ = consumes<edm::View<reco::Jet> >(parameters.getParameter<InputTag>("PFJetCollection"));
  theCaloMETCollectionToken_ = consumes<edm::View<reco::MET> >(parameters.getParameter<InputTag>("caloMETCollection"));
  theVertexToken_ = consumes<reco::VertexCollection>(parameters.getParameter<InputTag>("vertexCollection"));
 
  // just to initialize
  isValidHltConfig_ = false;
 
  h_vertex_number = nullptr;
  h_vertex_chi2 = nullptr;
  h_vertex_numTrks = nullptr;
  h_vertex_sumTrks = nullptr;
  h_vertex_d0 = nullptr;
 
  h_jet_count = nullptr;
  h_jet_et = nullptr;
  h_jet_pt = nullptr;
  h_jet_eta = nullptr;
  h_jet_phi = nullptr;
  h_jet2_et = nullptr;
  // h_jet2_pt = 0;
  h_jet2_eta = nullptr;
  h_jet2_phi = nullptr;
 
  h_e1_et = nullptr;
  h_e2_et = nullptr;
  h_e1_eta = nullptr;
  h_e2_eta = nullptr;
  h_e1_phi = nullptr;
  h_e2_phi = nullptr;
 
  h_m1_pt = nullptr;
  h_m2_pt = nullptr;
  h_m1_eta = nullptr;
  h_m2_eta = nullptr;
  h_m1_phi = nullptr;
  h_m2_phi = nullptr;
 
  // h_t1_et = 0;
  // h_t1_eta = 0;
  // h_t1_phi = 0;
 
  h_met = nullptr;
  h_met_phi = nullptr;
 
  h_e_invWMass = nullptr;
  h_m_invWMass = nullptr;
  h_mumu_invMass = nullptr;
  h_ee_invMass = nullptr;
}

~EwkDQM()

EwkDQM::~EwkDQM ( )

override

Destructor.

Definition at line 109 of file EwkDQM.cc.

{}

Member Function Documentation

analyze()

void EwkDQM::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overridevirtual

Get the analysis.

Reimplemented from DQMEDAnalyzer.

Definition at line 203 of file EwkDQM.cc.

                                                                  {
  // short-circuit if hlt problems
  if (!isValidHltConfig_)
    return;
 
  LogTrace(logTraceName) << "Analysis of event # ";
  // Did it pass certain HLT path?
  Handle<TriggerResults> HLTresults;
  iEvent.getByToken(theTriggerResultsToken_, HLTresults);
  if (!HLTresults.isValid())
    return;
 
  const edm::TriggerNames& trigNames = iEvent.triggerNames(*HLTresults);
 
  // a temporary, until we have a list of triggers of interest
  std::vector<std::string> eleTrigPathNames;
  std::vector<std::string> muTrigPathNames;
 
  // eleTrigPathNames.push_back(theElecTriggerPathToPass_);
  // muTrigPathNames.push_back(theMuonTriggerPathToPass_);
  // end of temporary
 
  bool passed_electron_HLT = false;
  bool passed_muon_HLT = false;
  for (unsigned int i = 0; i < HLTresults->size(); i++) {
    const std::string& trigName = trigNames.triggerName(i);
    // check if triggerName matches electronPath
    for (unsigned int index = 0; index < theElecTriggerPathToPass_.size() && !passed_electron_HLT; index++) {
      // 0 if found, pos if not
      size_t trigPath = trigName.find(theElecTriggerPathToPass_[index]);
      if (trigPath == 0) {
        //      cout << "MuonTrigger passed (=trigName): " << trigName <<endl;
        passed_electron_HLT = HLTresults->accept(i);
      }
    }
    // check if triggerName matches muonPath
    for (unsigned int index = 0; index < theMuonTriggerPathToPass_.size() && !passed_muon_HLT; index++) {
      // 0 if found, pos if not
      size_t trigPath = trigName.find(theMuonTriggerPathToPass_[index]);
      if (trigPath == 0) {
        //      cout << "MuonTrigger passed (=trigName): " << trigName <<endl;
        passed_muon_HLT = HLTresults->accept(i);
      }
    }
  }
 
  // we are interested in events with a valid electron or muon
  if (!(passed_electron_HLT || passed_muon_HLT))
    return;
 
  // Vertex information
  Handle<VertexCollection> vertexHandle;
  iEvent.getByToken(theVertexToken_, vertexHandle);
  if (!vertexHandle.isValid())
    return;
  VertexCollection vertexCollection = *(vertexHandle.product());
  VertexCollection::const_iterator v = vertexCollection.begin();
  int vertex_number = vertexCollection.size();
  double vertex_chi2 = v->normalizedChi2();  // v->chi2();
  double vertex_d0 = sqrt(v->x() * v->x() + v->y() * v->y());
  double vertex_numTrks = v->tracksSize();
  double vertex_sumTrks = 0.0;
  // std::cout << "vertex_d0=" << vertex_d0 << "\n";
  // double vertex_ndof    = v->ndof();cout << "ndof="<<vertex_ndof<<endl;
  for (Vertex::trackRef_iterator vertex_curTrack = v->tracks_begin(); vertex_curTrack != v->tracks_end();
       vertex_curTrack++)
    vertex_sumTrks += (*vertex_curTrack)->pt();
 
  // Missing ET
  Handle<View<MET> > caloMETCollection;
  iEvent.getByToken(theCaloMETCollectionToken_, caloMETCollection);
  if (!caloMETCollection.isValid())
    return;
  float missing_et = caloMETCollection->begin()->et();
  float met_phi = caloMETCollection->begin()->phi();
 
  // grab "gaussian sum fitting" electrons
  Handle<GsfElectronCollection> electronCollection;
  iEvent.getByToken(theElectronCollectionLabel_, electronCollection);
  if (!electronCollection.isValid())
    return;
 
  // Find the highest and 2nd highest electron
  float electron_et = -8.0;
  float electron_eta = -8.0;
  float electron_phi = -8.0;
  float electron2_et = -9.0;
  float electron2_eta = -9.0;
  float electron2_phi = -9.0;
  float ee_invMass = -9.0;
  TLorentzVector e1, e2;
 
  // If it passed electron HLT and the collection was found, find electrons near
  // Z mass
  if (passed_electron_HLT) {
    for (reco::GsfElectronCollection::const_iterator recoElectron = electronCollection->begin();
         recoElectron != electronCollection->end();
         recoElectron++) {
      // Require electron to pass some basic cuts
      if (recoElectron->et() < 20 || fabs(recoElectron->eta()) > 2.5)
        continue;
 
      // Tighter electron cuts
      if (recoElectron->deltaPhiSuperClusterTrackAtVtx() > 0.58 ||
          recoElectron->deltaEtaSuperClusterTrackAtVtx() > 0.01 || recoElectron->sigmaIetaIeta() > 0.027)
        continue;
 
      if (recoElectron->et() > electron_et) {
        electron2_et = electron_et;  // 2nd highest gets values from current highest
        electron2_eta = electron_eta;
        electron2_phi = electron_phi;
        electron_et = recoElectron->et();  // 1st highest gets values from new highest
        electron_eta = recoElectron->eta();
        electron_phi = recoElectron->phi();
        e1 = TLorentzVector(recoElectron->momentum().x(),
                            recoElectron->momentum().y(),
                            recoElectron->momentum().z(),
                            recoElectron->p());
      } else if (recoElectron->et() > electron2_et) {
        electron2_et = recoElectron->et();
        electron2_eta = recoElectron->eta();
        electron2_phi = recoElectron->phi();
        e2 = TLorentzVector(recoElectron->momentum().x(),
                            recoElectron->momentum().y(),
                            recoElectron->momentum().z(),
                            recoElectron->p());
      }
    }  // end of loop over electrons
    if (electron2_et > 0.0) {
      TLorentzVector pair = e1 + e2;
      ee_invMass = pair.M();
    }
  }  // end of "are electrons valid"
 
  // Take the STA muon container
  Handle<MuonCollection> muonCollection;
  iEvent.getByToken(theMuonCollectionLabel_, muonCollection);
  if (!muonCollection.isValid())
    return;
 
  // Find the highest pt muons
  float mm_invMass = -9.0;
  float muon_pt = -9.0;
  float muon_eta = -9.0;
  float muon_phi = -9.0;
  float muon2_pt = -9.0;
  float muon2_eta = -9.0;
  float muon2_phi = -9.0;
  TLorentzVector m1, m2;
 
  if (passed_muon_HLT) {
    for (reco::MuonCollection::const_iterator recoMuon = muonCollection->begin(); recoMuon != muonCollection->end();
         recoMuon++) {
      // Require muon to pass some basic cuts
      if (recoMuon->pt() < 20 || !recoMuon->isGlobalMuon())
        continue;
      // Some tighter muon cuts
      if (recoMuon->globalTrack()->normalizedChi2() > 10)
        continue;
 
      if (recoMuon->pt() > muon_pt) {
        muon2_pt = muon_pt;  // 2nd highest gets values from current highest
        muon2_eta = muon_eta;
        muon2_phi = muon_phi;
        muon_pt = recoMuon->pt();  // 1st highest gets values from new highest
        muon_eta = recoMuon->eta();
        muon_phi = recoMuon->phi();
        m1 =
            TLorentzVector(recoMuon->momentum().x(), recoMuon->momentum().y(), recoMuon->momentum().z(), recoMuon->p());
      } else if (recoMuon->pt() > muon2_pt) {
        muon2_pt = recoMuon->pt();
        muon2_eta = recoMuon->eta();
        muon2_phi = recoMuon->phi();
        m2 =
            TLorentzVector(recoMuon->momentum().x(), recoMuon->momentum().y(), recoMuon->momentum().z(), recoMuon->p());
      }
    }
  }
  if (muon2_pt > 0.0) {
    TLorentzVector pair = m1 + m2;
    mm_invMass = pair.M();
  }
 
  // Find the highest et jet
 
  //  Handle<CaloJetCollection> caloJetCollection;
  Handle<View<Jet> > PFJetCollection;
  iEvent.getByToken(thePFJetCollectionToken_, PFJetCollection);
  if (!PFJetCollection.isValid())
    return;
 
  unsigned int muonCollectionSize = muonCollection->size();
  // unsigned int jetCollectionSize = jetCollection->size();
  unsigned int PFJetCollectionSize = PFJetCollection->size();
  int jet_count = 0;
  // int LEADJET=-1;  double max_pt=0;
 
  float jet_et = -80.0;
  float jet_pt = -80.0;   // prova
  float jet_eta = -80.0;  // now USED
  float jet_phi = -80.0;  // now USED
  float jet2_et = -90.0;
  float jet2_eta = -90.0;  // now USED
  float jet2_phi = -90.0;  // now USED
  //  for (CaloJetCollection::const_iterator i_calojet =
  // caloJetCollection->begin();
  //       i_calojet != caloJetCollection->end(); i_calojet++) {
  //  for (PFJetCollection::const_iterator i_pfjet = PFJetCollection->begin();
  //       i_pfjet != PFJetCollection->end(); i_pfjet++) {
  //  float jet_current_et = i_calojet->et();
  //  float jet_current_et = i_pfjet->et();            // e` identico a jet.et()
  //    jet_count++;
 
  // cleaning: va messo prima del riempimento dell'istogramma // This is in
  // order to use PFJets
  for (unsigned int i = 0; i < PFJetCollectionSize; i++) {
    const Jet& jet = PFJetCollection->at(i);
    // la classe "jet" viene definita qui!!!
    double minDistance = 99999;
    for (unsigned int j = 0; j < muonCollectionSize; j++) {
      const Muon& mu = muonCollection->at(j);
      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
    // se la distanza muone-cono del jet e` minore di 0.3, passo avanti e non
    // conteggio il mio jet
 
    // If it overlaps with ELECTRON, it is not a jet
    if (electron_et > 0.0 && fabs(jet.eta() - electron_eta) < 0.2 && calcDeltaPhi(jet.phi(), electron_phi) < 0.2)
      continue;
    if (electron2_et > 0.0 && fabs(jet.eta() - electron2_eta) < 0.2 && calcDeltaPhi(jet.phi(), electron2_phi) < 0.2)
      continue;
 
    // provo a cambiare la parte degli elettroni in modo simmetrico alla parte
    // per i muoni
 
    // ...
    // ...
 
    // if it has too low Et, throw away
    if (jet.et() < eJetMin_)
      continue;
    jet_count++;
 
    // ovvero: incrementa jet_count se:
    //   - non c'e un muone entro 0.3 di distanza dal cono del jet;
    //   - se il jet non si sovrappone ad un elettrone;
    //   - se l'energia trasversa e` maggiore della soglia impostata (15?)
 
    // if(jet.et()>max_pt) { LEADJET=i; max_pt=jet.et();}
    // se l'energia del jet e` maggiore di max_pt, diventa "i"
    // l'indice del jet piu` energetico e max_pt la sua energia
 
    // riguardare questo!!!
    // fino ad ora, jet_et era inizializzato a -8.0
    if (jet.et() > jet_et) {
      jet2_et = jet_et;  // 2nd highest jet gets et from current highest
      // perche` prende l'energia del primo jet??
      jet2_eta = jet_eta;  // now USED
      jet2_phi = jet_phi;  // now USED
      // jet_et = i_calojet->et(); // current highest jet gets
      // et from the new highest
      jet_et = jet.et();  // current highest jet gets et from the new highest
      // ah, ok! lo riaggiorna solo dopo!
      jet_pt = jet.pt();                          // e` il pT del leading jet
      jet_eta = jet.eta();                        // now USED
      jet_phi = jet.phi() * (Geom::pi() / 180.);  // now USED
    } else if (jet.et() > jet2_et) {
      //      jet2_et  = i_calojet->et();
      jet2_et = jet.et();
      //      jet2_eta = i_calojet->eta();  // UNUSED
      //      jet2_phi = i_calojet->phi();  // UNUSED
      jet2_eta = jet.eta();  // now USED
      jet2_phi = jet.phi();  // now USED
    }
    // questo elseif funziona
  }
 
  //                 Fill Histograms //
 
  bool fill_e1 = false;
  bool fill_e2 = false;
  bool fill_m1 = false;
  bool fill_m2 = false;
  bool fill_met = false;
 
  // Was Z->ee found?
  if (ee_invMass > 0.0) {
    h_ee_invMass->Fill(ee_invMass);
    fill_e1 = true;
    fill_e2 = true;
  }
 
  // Was Z->mu mu found?
  if (mm_invMass > 0.0) {
    h_mumu_invMass->Fill(mm_invMass);
    fill_m1 = true;
    fill_m2 = true;
    h_jet2_et->Fill(jet2_et);
  }
 
  // Was W->e nu found?
  if (electron_et > 0.0 && missing_et > 20.0) {
    float dphiW = fabs(met_phi - electron_phi);
    float W_mt_e = sqrt(2 * missing_et * electron_et * (1 - cos(dphiW)));
    h_e_invWMass->Fill(W_mt_e);
    fill_e1 = true;
    fill_met = true;
  }
 
  // Was W->mu nu found?
  if (muon_pt > 0.0 && missing_et > 20.0) {
    float dphiW = fabs(met_phi - muon_phi);
    float W_mt_m = sqrt(2 * missing_et * muon_pt * (1 - cos(dphiW)));
    h_m_invWMass->Fill(W_mt_m);
    fill_m1 = true;
    fill_met = true;
  }
 
  if (jet_et > -10.0) {
    h_jet_et->Fill(jet_et);
    h_jet_count->Fill(jet_count);
  }
 
  if (jet_pt > 0.) {
    h_jet_pt->Fill(jet_pt);
  }
 
  if (jet_eta > -50.) {
    h_jet_eta->Fill(jet_eta);
  }
 
  if (jet_phi > -10.) {
    h_jet_phi->Fill(jet_phi);
  }
 
  if (jet2_et > -10.0) {
    h_jet2_et->Fill(jet2_et);
  }
 
  // if (jet2_pt>0.) {
  //  h_jet2_pt   ->Fill(jet2_pt);
  // }
 
  if (jet2_eta > -50.) {
    h_jet2_eta->Fill(jet2_eta);
  }
 
  if (jet2_phi > -10.) {
    h_jet2_phi->Fill(jet2_phi);
  }
 
  if (fill_e1 || fill_m1) {
    h_vertex_number->Fill(vertex_number);
    h_vertex_chi2->Fill(vertex_chi2);
    h_vertex_d0->Fill(vertex_d0);
    h_vertex_numTrks->Fill(vertex_numTrks);
    h_vertex_sumTrks->Fill(vertex_sumTrks);
  }
 
  if (fill_e1) {
    h_e1_et->Fill(electron_et);
    h_e1_eta->Fill(electron_eta);
    h_e1_phi->Fill(electron_phi);
  }
  if (fill_e2) {
    h_e2_et->Fill(electron2_et);
    h_e2_eta->Fill(electron2_eta);
    h_e2_phi->Fill(electron2_phi);
  }
  if (fill_m1) {
    h_m1_pt->Fill(muon_pt);
    h_m1_eta->Fill(muon_eta);
    h_m1_phi->Fill(muon_phi);
  }
  if (fill_m2) {
    h_m2_pt->Fill(muon2_pt);
    h_m2_eta->Fill(muon2_eta);
    h_m2_phi->Fill(muon2_phi);
  }
  if (fill_met) {
    h_met->Fill(missing_et);
    h_met_phi->Fill(met_phi);
  }
}

References edm::HLTGlobalStatus::accept(), calcDeltaPhi(), ewkDQM_cfi::caloMETCollection, funct::cos(), HLT_2018_cff::distance, StorageManager_cfg::e1, pdwgLeptonRecoSkim_cfi::electronCollection, mps_fire::i, iEvent, edm::HandleBase::isValid(), dqmiolumiharvest::j, metsig::jet, LogTrace, amptDefaultParameters_cff::mu, pdwgLeptonRecoSkim_cfi::muonCollection, Geom::pi(), edm::Handle< T >::product(), edm::HLTGlobalStatus::size(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, EgHLTOffTrigSelection_cfi::trigName, trigNames, findQualityFiles::v, and spclusmultinvestigator_cfi::vertexCollection.

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 111 of file EwkDQM.cc.

                                                                                         {
  ibooker.setCurrentFolder("Physics/EwkDQM");
 
  char chtitle[256] = "";
  const size_t title_s = sizeof(chtitle);
 
  logTraceName = "EwkAnalyzer";
 
  LogTrace(logTraceName) << "Parameters initialization";
 
  const float pi = 4 * atan(1);
 
  // Keep the number of plots and number of bins to a minimum!
 
  h_vertex_number = ibooker.book1D("vertex_number", "Number of event vertices in collection", 10, -0.5, 9.5);
  h_vertex_chi2 = ibooker.book1D("vertex_chi2", "Event Vertex #chi^{2}/n.d.o.f.", 20, 0.0, 2.0);
  h_vertex_numTrks = ibooker.book1D("vertex_numTrks", "Event Vertex, number of tracks", 20, -0.5, 59.5);
  h_vertex_sumTrks = ibooker.book1D("vertex_sumTrks", "Event Vertex, sum of track pt", 20, 0.0, 100.0);
  h_vertex_d0 = ibooker.book1D("vertex_d0", "Event Vertex d0", 20, 0.0, 0.05);
 
  h_jet_count = ibooker.book1D("jet_count", chtitle, 8, -0.5, 7.5);
 
  snprintf(
      chtitle, title_s, "Leading jet E_{T} (from %s);E_{T}(1^{st} jet) (GeV)", thePFJetCollectionLabel_.label().data());
  h_jet_et = ibooker.book1D("jet_et", chtitle, 20, 0., 200.0);
 
  snprintf(chtitle,
           title_s,
           "Leading jet p_{T} (from %s);p_{T}(1^{st} jet) (GeV/c)",
           thePFJetCollectionLabel_.label().data());
  h_jet_pt = ibooker.book1D("jet_pt", chtitle, 20, 0., 200.0);
 
  snprintf(chtitle, title_s, "Leading jet #eta (from %s); #eta (1^{st} jet)", thePFJetCollectionLabel_.label().data());
  h_jet_eta = ibooker.book1D("jet_eta", chtitle, 20, -10., 10.0);
 
  snprintf(chtitle, title_s, "Leading jet #phi (from %s); #phi(1^{st} jet)", thePFJetCollectionLabel_.label().data());
  h_jet_phi = ibooker.book1D("jet_phi", chtitle, 22, -1.1 * pi, 1.1 * pi);
 
  snprintf(chtitle,
           title_s,
           "2^{nd} leading jet E_{T} (from %s);E_{T}(2^{nd} jet) (GeV)",
           thePFJetCollectionLabel_.label().data());
  h_jet2_et = ibooker.book1D("jet2_et", chtitle, 20, 0., 200.0);
 
  snprintf(chtitle,
           title_s,
           "2^{nd} leading jet #eta (from %s); #eta (2^{nd} jet)",
           thePFJetCollectionLabel_.label().data());
  h_jet2_eta = ibooker.book1D("jet2_eta", chtitle, 20, -10., 10.0);
 
  snprintf(
      chtitle, title_s, "2^{nd} leading jet #phi (from %s); #phi(2^{nd} jet)", thePFJetCollectionLabel_.label().data());
  h_jet2_phi = ibooker.book1D("jet2_phi", chtitle, 22, -1.1 * pi, 1.1 * pi);
 
  h_e1_et = ibooker.book1D("e1_et", "E_{T} of Leading Electron;E_{T} (GeV)", 20, 0.0, 100.0);
  h_e2_et = ibooker.book1D("e2_et", "E_{T} of Second Electron;E_{T} (GeV)", 20, 0.0, 100.0);
  h_e1_eta = ibooker.book1D("e1_eta", "#eta of Leading Electron;#eta", 20, -4.0, 4.0);
 
  h_e2_eta = ibooker.book1D("e2_eta", "#eta of Second Electron;#eta", 20, -4.0, 4.0);
 
  h_e1_phi = ibooker.book1D("e1_phi", "#phi of Leading Electron;#phi", 22, -1.1 * pi, 1.1 * pi);
  h_e2_phi = ibooker.book1D("e2_phi", "#phi of Second Electron;#phi", 22, -1.1 * pi, 1.1 * pi);
  h_m1_pt = ibooker.book1D("m1_pt", "p_{T} of Leading Muon;p_{T}(1^{st} #mu) (GeV)", 20, 0.0, 100.0);
  h_m2_pt = ibooker.book1D("m2_pt", "p_{T} of Second Muon;p_{T}(2^{nd} #mu) (GeV)", 20, 0.0, 100.0);
  h_m1_eta = ibooker.book1D("m1_eta", "#eta of Leading Muon;#eta(1^{st} #mu)", 20, -4.0, 4.0);
  h_m2_eta = ibooker.book1D("m2_eta", "#eta of Second Muon;#eta(2^{nd} #mu)", 20, -4.0, 4.0);
  h_m1_phi = ibooker.book1D(
      "m1_phi", "#phi of Leading Muon;#phi(1^{st} #mu)", 20, (-1. - 1. / 10.) * pi, (1. + 1. / 10.) * pi);
  h_m2_phi =
      ibooker.book1D("m2_phi", "#phi of Second Muon;#phi(2^{nd} #mu)", 20, (-1. - 1. / 10.) * pi, (1. + 1. / 10.) * pi);
 
  snprintf(chtitle, title_s, "Missing E_{T} (%s); GeV", theCaloMETCollectionLabel_.label().data());
  h_met = ibooker.book1D("met", chtitle, 20, 0.0, 100);
  h_met_phi =
      ibooker.book1D("met_phi", "Missing E_{T} #phi;#phi(MET)", 22, (-1. - 1. / 10.) * pi, (1. + 1. / 10.) * pi);
 
  h_e_invWMass = ibooker.book1D("we_invWMass", "W-> e #nu Transverse Mass;M_{T} (GeV)", 20, 0.0, 140.0);
  h_m_invWMass = ibooker.book1D("wm_invWMass", "W-> #mu #nu Transverse Mass;M_{T} (GeV)", 20, 0.0, 140.0);
  h_mumu_invMass = ibooker.book1D("z_mm_invMass", "#mu#mu Invariant Mass;InvMass (GeV)", 20, 40.0, 140.0);
  h_ee_invMass = ibooker.book1D("z_ee_invMass", "ee Invariant Mass;InvMass (Gev)", 20, 40.0, 140.0);
}

References dqm::implementation::IBooker::book1D(), LogTrace, pi, and dqm::implementation::NavigatorBase::setCurrentFolder().

calcDeltaPhi()

double EwkDQM::calcDeltaPhi	(	double	phi1,
		double	phi2
	)

Definition at line 605 of file EwkDQM.cc.

                                                    {
  double deltaPhi = phi1 - phi2;
 
  if (deltaPhi < 0)
    deltaPhi = -deltaPhi;
 
  if (deltaPhi > 3.1415926)
    deltaPhi = 2 * 3.1415926 - deltaPhi;
 
  return deltaPhi;
}

References SiPixelRawToDigiRegional_cfi::deltaPhi.

dqmBeginRun()

void EwkDQM::dqmBeginRun ( const edm::Run &  theRun, const edm::EventSetup &  theSetup  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 194 of file EwkDQM.cc.

                                                                            {
  // passed as parameter to HLTConfigProvider::init(), not yet used
  bool isConfigChanged = false;
 
  // isValidHltConfig_ used to short-circuit analyze() in case of problems
  const std::string hltProcessName(theTriggerResultsCollection_.process());
  isValidHltConfig_ = hltConfigProvider_.init(theRun, theSetup, hltProcessName, isConfigChanged);
}

References HigPhotonJetHLTOfflineSource_cfi::hltProcessName, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

eJetMin_

double EwkDQM::eJetMin_

private

Definition at line 69 of file EwkDQM.h.

h_e1_et

MonitorElement* EwkDQM::h_e1_et

private

Definition at line 89 of file EwkDQM.h.

h_e1_eta

MonitorElement* EwkDQM::h_e1_eta

private

Definition at line 91 of file EwkDQM.h.

h_e1_phi

MonitorElement* EwkDQM::h_e1_phi

private

Definition at line 93 of file EwkDQM.h.

h_e2_et

MonitorElement* EwkDQM::h_e2_et

private

Definition at line 90 of file EwkDQM.h.

h_e2_eta

MonitorElement* EwkDQM::h_e2_eta

private

Definition at line 92 of file EwkDQM.h.

h_e2_phi

MonitorElement* EwkDQM::h_e2_phi

private

Definition at line 94 of file EwkDQM.h.

h_e_invWMass

MonitorElement* EwkDQM::h_e_invWMass

private

Definition at line 110 of file EwkDQM.h.

h_ee_invMass

MonitorElement* EwkDQM::h_ee_invMass

private

Definition at line 113 of file EwkDQM.h.

h_jet2_et

MonitorElement* EwkDQM::h_jet2_et

private

Definition at line 84 of file EwkDQM.h.

h_jet2_eta

MonitorElement* EwkDQM::h_jet2_eta

private

Definition at line 86 of file EwkDQM.h.

h_jet2_phi

MonitorElement* EwkDQM::h_jet2_phi

private

Definition at line 87 of file EwkDQM.h.

h_jet_count

MonitorElement* EwkDQM::h_jet_count

private

Definition at line 78 of file EwkDQM.h.

h_jet_et

MonitorElement* EwkDQM::h_jet_et

private

Definition at line 79 of file EwkDQM.h.

h_jet_eta

MonitorElement* EwkDQM::h_jet_eta

private

Definition at line 81 of file EwkDQM.h.

h_jet_phi

MonitorElement* EwkDQM::h_jet_phi

private

Definition at line 82 of file EwkDQM.h.

h_jet_pt

MonitorElement* EwkDQM::h_jet_pt

private

Definition at line 80 of file EwkDQM.h.

h_m1_eta

MonitorElement* EwkDQM::h_m1_eta

private

Definition at line 98 of file EwkDQM.h.

h_m1_phi

MonitorElement* EwkDQM::h_m1_phi

private

Definition at line 100 of file EwkDQM.h.

h_m1_pt

MonitorElement* EwkDQM::h_m1_pt

private

Definition at line 96 of file EwkDQM.h.

h_m2_eta

MonitorElement* EwkDQM::h_m2_eta

private

Definition at line 99 of file EwkDQM.h.

h_m2_phi

MonitorElement* EwkDQM::h_m2_phi

private

Definition at line 101 of file EwkDQM.h.

h_m2_pt

MonitorElement* EwkDQM::h_m2_pt

private

Definition at line 97 of file EwkDQM.h.

h_m_invWMass

MonitorElement* EwkDQM::h_m_invWMass

private

Definition at line 111 of file EwkDQM.h.

h_met

MonitorElement* EwkDQM::h_met

private

Definition at line 107 of file EwkDQM.h.

h_met_phi

MonitorElement* EwkDQM::h_met_phi

private

Definition at line 108 of file EwkDQM.h.

h_mumu_invMass

MonitorElement* EwkDQM::h_mumu_invMass

private

Definition at line 112 of file EwkDQM.h.

h_vertex_chi2

MonitorElement* EwkDQM::h_vertex_chi2

private

Definition at line 73 of file EwkDQM.h.

h_vertex_d0

MonitorElement* EwkDQM::h_vertex_d0

private

Definition at line 76 of file EwkDQM.h.

h_vertex_number

MonitorElement* EwkDQM::h_vertex_number

private

Definition at line 72 of file EwkDQM.h.

h_vertex_numTrks

MonitorElement* EwkDQM::h_vertex_numTrks

private

Definition at line 74 of file EwkDQM.h.

h_vertex_sumTrks

MonitorElement* EwkDQM::h_vertex_sumTrks

private

Definition at line 75 of file EwkDQM.h.

hltConfigProvider_

HLTConfigProvider EwkDQM::hltConfigProvider_

private

Definition at line 53 of file EwkDQM.h.

isValidHltConfig_

bool EwkDQM::isValidHltConfig_

private

Definition at line 54 of file EwkDQM.h.

logTraceName

std::string EwkDQM::logTraceName

private

Definition at line 51 of file EwkDQM.h.

theCaloMETCollectionLabel_

edm::InputTag EwkDQM::theCaloMETCollectionLabel_

private

Definition at line 60 of file EwkDQM.h.

theCaloMETCollectionToken_

edm::EDGetTokenT<edm::View<reco::MET> > EwkDQM::theCaloMETCollectionToken_

private

Definition at line 66 of file EwkDQM.h.

theElecTriggerPathToPass_

std::vector<std::string> EwkDQM::theElecTriggerPathToPass_

private

Definition at line 57 of file EwkDQM.h.

theElectronCollectionLabel_

edm::EDGetTokenT<reco::GsfElectronCollection> EwkDQM::theElectronCollectionLabel_

private

Definition at line 64 of file EwkDQM.h.

theMuonCollectionLabel_

edm::EDGetTokenT<reco::MuonCollection> EwkDQM::theMuonCollectionLabel_

private

Definition at line 63 of file EwkDQM.h.

theMuonTriggerPathToPass_

std::vector<std::string> EwkDQM::theMuonTriggerPathToPass_

private

Definition at line 58 of file EwkDQM.h.

thePFJetCollectionLabel_

edm::InputTag EwkDQM::thePFJetCollectionLabel_

private

Definition at line 59 of file EwkDQM.h.

thePFJetCollectionToken_

edm::EDGetTokenT<edm::View<reco::Jet> > EwkDQM::thePFJetCollectionToken_

private

Definition at line 65 of file EwkDQM.h.

theTriggerResultsCollection_

edm::InputTag EwkDQM::theTriggerResultsCollection_

private

Definition at line 61 of file EwkDQM.h.

theTriggerResultsToken_

edm::EDGetTokenT<edm::TriggerResults> EwkDQM::theTriggerResultsToken_

private

Definition at line 62 of file EwkDQM.h.

theVertexToken_

edm::EDGetTokenT<reco::VertexCollection> EwkDQM::theVertexToken_

private

Definition at line 67 of file EwkDQM.h.