#include "DQM/Physics/src/EwkTauDQM.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/View.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Common/interface/TriggerNames.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/TauReco/interface/PFTau.h"
#include "DataFormats/TauReco/interface/PFTauFwd.h"
#include "DataFormats/TauReco/interface/PFTauDiscriminator.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETFwd.h"
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "TMath.h"
#include <iostream>
#include <iomanip>
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"

Functions
double	calcDeltaPhi (double phi1, double phi2)

double	calcMt (double px1, double py1, double px2, double py2)

double	calcPzeta (const reco::Candidate::LorentzVector &p1, const reco::Candidate::LorentzVector &p2, double pxMEt, double pyMEt)

std::string	dqmDirectoryName (const std::string &dqmRootDirectory, const std::string &dqmSubDirectory)

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

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

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

const reco::PFTau *	getTheTauJet (const reco::PFTauCollection &tauJets, double tauJetEtaCut, double tauJetPtCut, int &theTauJetIndex)

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

bool	passesElectronId (const reco::GsfElectron &electron)

bool	passesElectronPreId (const reco::GsfElectron &electron)

Variables
const std::string	dqmSeparator = "/"

Function Documentation

◆ calcDeltaPhi()

double calcDeltaPhi	(	double	phi1,
		double	phi2
	)

Definition at line 1023 of file EwkTauDQM.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, and Pi.

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

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

◆ calcMt()

double calcMt	(	double	px1,
		double	py1,
		double	px2,
		double	py2
	)

Definition at line 1035 of file EwkTauDQM.cc.

References HLT_2022v15_cff::pt1, and HLT_2022v15_cff::pt2.

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

                                                               {
   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));
 }

◆ calcPzeta()

double calcPzeta	(	const reco::Candidate::LorentzVector &	p1,
		const reco::Candidate::LorentzVector &	p2,
		double	pxMEt,
		double	pyMEt
	)

Definition at line 1045 of file EwkTauDQM.cc.

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

                                {
   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;
 }

◆ dqmDirectoryName()

std::string dqmDirectoryName	(	const std::string &	dqmRootDirectory,
		const std::string &	dqmSubDirectory
	)

Definition at line 10 of file EwkTauDQM.cc.

References TrackerOfflineValidation_Dqm_cff::dirName, dqmRootDirectory, dqmSeparator, ewkTauDQM_cfi::dqmSubDirectory, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by TauDQMFileLoader::endRun(), TauDQMHistPlotter::endRun(), EwkTauDQM::EwkTauDQM(), and mapSubDirectoryStructure().

                                                                                                 {
   //--- concatenate names of dqmRootDirectory and dqmSubDirectory;
   //    add "/" separator inbetween if necessary
   std::string dirName = dqmRootDirectory;
   if (!dirName.empty() && dirName.find_last_of(dqmSeparator) != (dirName.length() - 1))
     dirName.append(dqmSeparator);
   dirName.append(dqmSubDirectory);
   return dirName;
 }

◆ getIsoMode()

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

Definition at line 1009 of file EwkTauDQM.cc.

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

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

                                                             {
   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;
 }

◆ getTheElectron()

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

Definition at line 1099 of file EwkTauDQM.cc.

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

Referenced by EwkElecTauHistManager::fillHistograms().

                                                               {
   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;
 }

◆ getTheMuon()

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

Definition at line 1116 of file EwkTauDQM.cc.

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

Referenced by EwkMuTauHistManager::fillHistograms().

                                                                                                  {
   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;
 }

◆ getTheTauJet()

const reco::PFTau* getTheTauJet	(	const reco::PFTauCollection &	tauJets,
		double	tauJetEtaCut,
		double	tauJetPtCut,
		int &	theTauJetIndex
	)

Definition at line 1129 of file EwkTauDQM.cc.

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

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

                                                      {
   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;
 }

◆ getVertexD0()

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

Definition at line 1151 of file EwkTauDQM.cc.

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

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

◆ passesElectronId()

bool passesElectronId ( const reco::GsfElectron & electron )

Definition at line 1082 of file EwkTauDQM.cc.

References HPSPFTauProducerPuppi_cfi::electron, and passesElectronPreId().

Referenced by EwkElecTauHistManager::fillHistograms().

                                                        {
   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;
   }
 }

◆ passesElectronPreId()

bool passesElectronPreId ( const reco::GsfElectron & electron )

Definition at line 1072 of file EwkTauDQM.cc.

References HPSPFTauProducerPuppi_cfi::electron.

Referenced by getTheElectron(), and passesElectronId().

                                                           {
   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;
   }
 }

Variable Documentation

◆ dqmSeparator

const std::string dqmSeparator = "/"

Definition at line 8 of file EwkTauDQM.cc.

Referenced by dqmDirectoryName(), TauDQMHistEffProducer::dqmEndJob(), dqmSubDirectoryName_merged(), and separateHistogramFromDirectoryName().

Functions

Variables