EwkTauDQM.h File Reference

#include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include <string>
#include <Math/VectorUtil.h>
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/TauReco/interface/PFTau.h"
#include "DataFormats/TauReco/interface/PFTauFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"

Go to the source code of this file.

Classes
class	EwkElecTauHistManager
class	EwkMuTauHistManager
class	EwkTauDQM

Enumerations
enum	{ kAbsoluteIso, kRelativeIso, kUndefinedIso }

Functions
double	calcDeltaPhi (double, double)
double	calcMt (double, double, double, double)
double	calcPzeta (const reco::Candidate::LorentzVector &, const reco::Candidate::LorentzVector &, double, double)
int	getIsoMode (const std::string &, int &)
const reco::GsfElectron *	getTheElectron (const reco::GsfElectronCollection &, double, double)
const reco::Muon *	getTheMuon (const reco::MuonCollection &, double, double)
const reco::PFTau *	getTheTauJet (const reco::PFTauCollection &, double, double, int &)
double	getVertexD0 (const reco::Vertex &, const reco::BeamSpot &)
bool	passesElectronId (const reco::GsfElectron &)
bool	passesElectronPreId (const reco::GsfElectron &)
template<typename T >
void	readEventData (const edm::Event &evt, const edm::InputTag &src, edm::Handle< T > &handle, long &numWarnings, int maxNumWarnings, bool &error, const char *errorMessage)

Enumeration Type Documentation

anonymous enum

anonymous enum

Auxiliary functions to compute quantities used by EWK Tau DQM (shared by different channels)

Author

Joshua Swanson

Enumerator
kAbsoluteIso
kRelativeIso
kUndefinedIso

Definition at line 370 of file EwkTauDQM.h.

{ kAbsoluteIso, kRelativeIso, kUndefinedIso };

Function Documentation

calcDeltaPhi()

double calcDeltaPhi	(	double	,
		double
	)

Definition at line 1026 of file EwkTauDQM.cc.

                                              {
  double deltaPhi = phi1 - phi2;
 
  if (deltaPhi < 0.)
    deltaPhi = -deltaPhi;
 
  if (deltaPhi > TMath::Pi())
    deltaPhi = 2 * TMath::Pi() - deltaPhi;
 
  return deltaPhi;
}

References SiPixelRawToDigiRegional_cfi::deltaPhi, and Pi.

Referenced by EwkDQM::analyze(), EwkElecTauHistManager::fillHistograms(), and EwkMuTauHistManager::fillHistograms().

calcMt()

double calcMt	(	double	,
		double	,
		double	,
		double
	)

Definition at line 1038 of file EwkTauDQM.cc.

                                                              {
  double pt1 = TMath::Sqrt(px1 * px1 + py1 * py1);
  double pt2 = TMath::Sqrt(px2 * px2 + py2 * py2);
 
  double p1Dotp2 = px1 * px2 + py1 * py2;
  double cosAlpha = p1Dotp2 / (pt1 * pt2);
 
  return TMath::Sqrt(2 * pt1 * pt2 * (1 - cosAlpha));
}

References HLT_2018_cff::pt1, and HLT_2018_cff::pt2.

Referenced by EwkElecTauHistManager::fillHistograms(), and EwkMuTauHistManager::fillHistograms().

calcPzeta()

double calcPzeta	(	const reco::Candidate::LorentzVector &	,
		const reco::Candidate::LorentzVector &	,
		double	,
		double
	)

Definition at line 1048 of file EwkTauDQM.cc.

                               {
  double cosPhi1 = cos(p1.phi());
  double sinPhi1 = sin(p1.phi());
  double cosPhi2 = cos(p2.phi());
  double sinPhi2 = sin(p2.phi());
  double zetaX = cosPhi1 + cosPhi2;
  double zetaY = sinPhi1 + sinPhi2;
  double zetaR = TMath::Sqrt(zetaX * zetaX + zetaY * zetaY);
  if (zetaR > 0.) {
    zetaX /= zetaR;
    zetaY /= zetaR;
  }
 
  double pxVis = p1.px() + p2.px();
  double pyVis = p1.py() + p2.py();
  double pZetaVis = pxVis * zetaX + pyVis * zetaY;
 
  double px = pxVis + pxMEt;
  double py = pyVis + pyMEt;
  double pZeta = px * zetaX + py * zetaY;
 
  return pZeta - 1.5 * pZetaVis;
}

References funct::cos(), p1, p2, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, and funct::sin().

getIsoMode()

int getIsoMode	(	const std::string &	,
		int &
	)

Definition at line 1012 of file EwkTauDQM.cc.

                                                            {
  int isoMode_int;
  if (isoMode_string == "absoluteIso") {
    isoMode_int = kAbsoluteIso;
  } else if (isoMode_string == "relativeIso") {
    isoMode_int = kRelativeIso;
  } else {
    edm::LogError("getIsoMode") << " Failed to decode isoMode string = " << isoMode_string << " !!";
    isoMode_int = kUndefinedIso;
    error = 1;
  }
  return isoMode_int;
}

References relativeConstraints::error, kAbsoluteIso, kRelativeIso, and kUndefinedIso.

Referenced by EwkElecTauHistManager::EwkElecTauHistManager(), and EwkMuTauHistManager::EwkMuTauHistManager().

getTheElectron()

const reco::GsfElectron* getTheElectron	(	const reco::GsfElectronCollection &	,
		double	,
		double
	)

Definition at line 1102 of file EwkTauDQM.cc.

                                                              {
  const reco::GsfElectron* theElectron = nullptr;
 
  for (reco::GsfElectronCollection::const_iterator electron = electrons.begin(); electron != electrons.end();
       ++electron) {
    if (TMath::Abs(electron->eta()) < electronEtaCut && electron->pt() > electronPtCut &&
        passesElectronPreId(*electron)) {
      if (theElectron == nullptr || electron->pt() > theElectron->pt())
        theElectron = &(*electron);
    }
  }
 
  return theElectron;
}

References Abs(), metsig::electron, ewkTauDQM_cfi::electronEtaCut, ewkTauDQM_cfi::electronPtCut, pwdgSkimBPark_cfi::electrons, passesElectronPreId(), and reco::LeafCandidate::pt().

Referenced by EwkElecTauHistManager::fillHistograms().

getTheMuon()

const reco::Muon* getTheMuon	(	const reco::MuonCollection &	,
		double	,
		double
	)

Definition at line 1119 of file EwkTauDQM.cc.

                                                                                                 {
  const reco::Muon* theMuon = nullptr;
 
  for (reco::MuonCollection::const_iterator muon = muons.begin(); muon != muons.end(); ++muon) {
    if (TMath::Abs(muon->eta()) < muonEtaCut && muon->pt() > muonPtCut) {
      if (theMuon == nullptr || muon->pt() > theMuon->pt())
        theMuon = &(*muon);
    }
  }
 
  return theMuon;
}

References Abs(), HLT_2018_cff::muon, ewkTauDQM_cfi::muonEtaCut, L1TObjectsTiming_cfi::muonPtCut, PDWG_BPHSkim_cff::muons, and reco::LeafCandidate::pt().

Referenced by EwkMuTauHistManager::fillHistograms().

getTheTauJet()

const reco::PFTau* getTheTauJet	(	const reco::PFTauCollection &	,
		double	,
		double	,
		int &
	)

Definition at line 1132 of file EwkTauDQM.cc.

                                                     {
  const reco::PFTau* theTauJet = nullptr;
  theTauJetIndex = -1;
 
  int numTauJets = tauJets.size();
  for (int iTauJet = 0; iTauJet < numTauJets; ++iTauJet) {
    const reco::PFTau& tauJet = tauJets.at(iTauJet);
 
    if (fabs(tauJet.eta()) < tauJetEtaCut && tauJet.pt() > tauJetPtCut) {
      if (theTauJet == nullptr || tauJet.pt() > theTauJet->pt()) {
        theTauJet = &tauJet;
        theTauJetIndex = iTauJet;
      }
    }
  }
 
  return theTauJet;
}

References reco::LeafCandidate::eta(), reco::LeafCandidate::pt(), reco::Candidate::size, ewkTauDQM_cfi::tauJetEtaCut, and ewkTauDQM_cfi::tauJetPtCut.

Referenced by EwkElecTauHistManager::fillHistograms(), and EwkMuTauHistManager::fillHistograms().

getVertexD0()

double getVertexD0	(	const reco::Vertex &	,
		const reco::BeamSpot &
	)

Definition at line 1154 of file EwkTauDQM.cc.

                                                                           {
  double dX = vertex.x() - beamSpot.x0();
  double dY = vertex.y() - beamSpot.y0();
  return TMath::Sqrt(dX * dX + dY * dY);
}

References pwdgSkimBPark_cfi::beamSpot, and bphysicsOniaDQM_cfi::vertex.

passesElectronId()

bool passesElectronId

(

const reco::GsfElectron &

)

Definition at line 1085 of file EwkTauDQM.cc.

                                                       {
  if (passesElectronPreId(electron) && ((TMath::Abs(electron.eta()) > 1.566 &&  // electron reconstructed in ECAL
                                                                                // endcap
                                         electron.sigmaEtaEta() < 0.03 && electron.hcalOverEcal() < 0.05 &&
                                         TMath::Abs(electron.deltaEtaSuperClusterTrackAtVtx()) < 0.009 &&
                                         TMath::Abs(electron.deltaPhiSuperClusterTrackAtVtx()) < 0.7) ||
                                        (TMath::Abs(electron.eta()) < 1.479 &&  // electron reconstructed in ECAL
                                                                                // barrel
                                         electron.sigmaEtaEta() < 0.01 && electron.hcalOverEcal() < 0.12 &&
                                         TMath::Abs(electron.deltaEtaSuperClusterTrackAtVtx()) < 0.007 &&
                                         TMath::Abs(electron.deltaPhiSuperClusterTrackAtVtx()) < 0.8))) {
    return true;
  } else {
    return false;
  }
}

References Abs(), metsig::electron, and passesElectronPreId().

Referenced by EwkElecTauHistManager::fillHistograms().

passesElectronPreId()

bool passesElectronPreId

(

const reco::GsfElectron &

)

Definition at line 1075 of file EwkTauDQM.cc.

                                                          {
  if ((TMath::Abs(electron.eta()) < 1.479 || TMath::Abs(electron.eta()) > 1.566) &&  // cut ECAL barrel/endcap crack
      electron.deltaPhiSuperClusterTrackAtVtx() < 0.8 && electron.deltaEtaSuperClusterTrackAtVtx() < 0.01 &&
      electron.sigmaIetaIeta() < 0.03) {
    return true;
  } else {
    return false;
  }
}

References Abs(), and metsig::electron.

Referenced by getTheElectron(), and passesElectronId().

readEventData()

template<typename T >

void readEventData	(	const edm::Event &	evt,
		const edm::InputTag &	src,
		edm::Handle< T > &	handle,
		long &	numWarnings,
		int	maxNumWarnings,
		bool &	error,
		const char *	errorMessage
	)

Definition at line 373 of file EwkTauDQM.h.

                                             {
  if (!evt.getByLabel(src, handle)) {
    if (numWarnings < maxNumWarnings || maxNumWarnings == -1)
      edm::LogWarning("readEventData") << errorMessage << " !!";
    ++numWarnings;
    error = true;
  }
}

References relativeConstraints::error, edm::Event::getByLabel(), patZpeak::handle, ewkTauDQM_cfi::maxNumWarnings, and TrackRefitter_38T_cff::src.

Referenced by EwkElecTauHistManager::fillHistograms(), and EwkMuTauHistManager::fillHistograms().