DQMExample_Step1 Class Reference

#include <DQMExample_Step1.h>

Inheritance diagram for DQMExample_Step1:

Public Member Functions
	DQMExample_Step1 (const edm::ParameterSet &ps)
	~DQMExample_Step1 () override
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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	bookHistos (DQMStore::IBooker &)
double	calcDeltaPhi (double phi1, double phi2)
double	Distance (const reco::Candidate &c1, const reco::Candidate &c2)
double	DistancePhi (const reco::Candidate &c1, const reco::Candidate &c2)
bool	MediumEle (const edm::Event &iEvent, const edm::EventSetup &iESetup, const reco::GsfElectron &electron)

Private Attributes
MonitorElement *	h_eEta_leading
MonitorElement *	h_eEta_leading_HLT
MonitorElement *	h_eEta_leading_HLT_matched
MonitorElement *	h_eEta_leading_matched
MonitorElement *	h_eMultiplicity
MonitorElement *	h_eMultiplicity_HLT
MonitorElement *	h_ePhi_leading
MonitorElement *	h_ePhi_leading_HLT
MonitorElement *	h_ePhi_leading_HLT_matched
MonitorElement *	h_ePhi_leading_matched
MonitorElement *	h_ePt_diff
MonitorElement *	h_ePt_leading
MonitorElement *	h_ePt_leading_HLT
MonitorElement *	h_ePt_leading_HLT_matched
MonitorElement *	h_ePt_leading_matched
MonitorElement *	h_jEta_leading
MonitorElement *	h_jMultiplicity
MonitorElement *	h_jPhi_leading
MonitorElement *	h_jPt_leading
MonitorElement *	h_pfMet
MonitorElement *	h_vertex_number
int	nBJets
int	nElectrons
double	ptThrJet_
double	ptThrL1_
double	ptThrL2_
double	ptThrMet_
math::XYZPoint	PVPoint_
edm::EDGetTokenT< reco::BeamSpot >	theBSCollection_
edm::EDGetTokenT< reco::CaloJetCollection >	theCaloJetCollection_
edm::EDGetTokenT< reco::ConversionCollection >	theConversionCollection_
edm::EDGetTokenT< reco::GsfElectronCollection >	theElectronCollection_
edm::EDGetTokenT< reco::PFMETCollection >	thePfMETCollection_
edm::EDGetTokenT< reco::VertexCollection >	thePVCollection_
edm::EDGetTokenT< trigger::TriggerEvent >	triggerEvent_
edm::InputTag	triggerFilter_
std::string	triggerPath_
edm::EDGetTokenT< edm::TriggerResults >	triggerResults_

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 >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 46 of file DQMExample_Step1.h.

Constructor & Destructor Documentation

DQMExample_Step1()

DQMExample_Step1::DQMExample_Step1

(

const edm::ParameterSet &

ps

)

Definition at line 21 of file DQMExample_Step1.cc.

                                                            {
  edm::LogInfo("DQMExample_Step1") << "Constructor  DQMExample_Step1::DQMExample_Step1 " << std::endl;
 
  // Get parameters from configuration file
  theElectronCollection_ = consumes<reco::GsfElectronCollection>(ps.getParameter<edm::InputTag>("electronCollection"));
  theCaloJetCollection_ = consumes<reco::CaloJetCollection>(ps.getParameter<edm::InputTag>("caloJetCollection"));
  thePfMETCollection_ = consumes<reco::PFMETCollection>(ps.getParameter<edm::InputTag>("pfMETCollection"));
  theConversionCollection_ =
      consumes<reco::ConversionCollection>(ps.getParameter<edm::InputTag>("conversionsCollection"));
  thePVCollection_ = consumes<reco::VertexCollection>(ps.getParameter<edm::InputTag>("PVCollection"));
  theBSCollection_ = consumes<reco::BeamSpot>(ps.getParameter<edm::InputTag>("beamSpotCollection"));
  triggerEvent_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("TriggerEvent"));
  triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
  triggerFilter_ = ps.getParameter<edm::InputTag>("TriggerFilter");
  triggerPath_ = ps.getParameter<std::string>("TriggerPath");
 
  // cuts:
  ptThrL1_ = ps.getUntrackedParameter<double>("PtThrL1");
  ptThrL2_ = ps.getUntrackedParameter<double>("PtThrL2");
  ptThrJet_ = ps.getUntrackedParameter<double>("PtThrJet");
  ptThrMet_ = ps.getUntrackedParameter<double>("PtThrMet");
}

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ptThrJet_, ptThrL1_, ptThrL2_, ptThrMet_, AlCaHLTBitMon_QueryRunRegistry::string, theBSCollection_, theCaloJetCollection_, theConversionCollection_, theElectronCollection_, thePfMETCollection_, thePVCollection_, triggerEvent_, triggerFilter_, triggerPath_, and triggerResults_.

~DQMExample_Step1()

DQMExample_Step1::~DQMExample_Step1 ( )

override

Definition at line 47 of file DQMExample_Step1.cc.

                                    {
  edm::LogInfo("DQMExample_Step1") << "Destructor DQMExample_Step1::~DQMExample_Step1 " << std::endl;
}

Member Function Documentation

analyze()

void DQMExample_Step1::analyze ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 72 of file DQMExample_Step1.cc.

                                                                             {
  edm::LogInfo("DQMExample_Step1") << "DQMExample_Step1::analyze" << std::endl;
 
  //-------------------------------
  //--- Vertex Info
  //-------------------------------
  edm::Handle<reco::VertexCollection> vertexHandle;
  e.getByToken(thePVCollection_, vertexHandle);
  if (!vertexHandle.isValid()) {
    edm::LogError("DQMClientExample") << "invalid collection: vertex"
                                      << "\n";
    return;
  }
 
  int vertex_number = vertexHandle->size();
  reco::VertexCollection::const_iterator v = vertexHandle->begin();
 
  math::XYZPoint PVPoint(-999, -999, -999);
  if (vertex_number != 0)
    PVPoint = math::XYZPoint(v->position().x(), v->position().y(), v->position().z());
 
  PVPoint_ = PVPoint;
 
  //-------------------------------
  //--- MET
  //-------------------------------
  edm::Handle<reco::PFMETCollection> pfMETCollection;
  e.getByToken(thePfMETCollection_, pfMETCollection);
  if (!pfMETCollection.isValid()) {
    edm::LogError("DQMClientExample") << "invalid collection: MET"
                                      << "\n";
    return;
  }
  //-------------------------------
  //--- Electrons
  //-------------------------------
  edm::Handle<reco::GsfElectronCollection> electronCollection;
  e.getByToken(theElectronCollection_, electronCollection);
  if (!electronCollection.isValid()) {
    edm::LogError("DQMClientExample") << "invalid collection: electrons"
                                      << "\n";
    return;
  }
 
  float nEle = 0;
  int posEle = 0, negEle = 0;
  const reco::GsfElectron *ele1 = nullptr;
  const reco::GsfElectron *ele2 = nullptr;
  for (reco::GsfElectronCollection::const_iterator recoElectron = electronCollection->begin();
       recoElectron != electronCollection->end();
       ++recoElectron) {
    // decreasing pT
    if (MediumEle(e, eSetup, *recoElectron)) {
      if (!ele1 && recoElectron->pt() > ptThrL1_)
        ele1 = &(*recoElectron);
 
      else if (!ele2 && recoElectron->pt() > ptThrL2_)
        ele2 = &(*recoElectron);
    }
 
    if (recoElectron->charge() == 1)
      posEle++;
    else if (recoElectron->charge() == -1)
      negEle++;
 
  }  // end of loop over electrons
 
  nEle = posEle + negEle;
 
  //-------------------------------
  //--- Jets
  //-------------------------------
  edm::Handle<reco::CaloJetCollection> caloJetCollection;
  e.getByToken(theCaloJetCollection_, caloJetCollection);
  if (!caloJetCollection.isValid()) {
    edm::LogError("DQMClientExample") << "invalid collection: jets"
                                      << "\n";
    return;
  }
 
  int nJet = 0;
  const reco::CaloJet *jet1 = nullptr;
  const reco::CaloJet *jet2 = nullptr;
 
  for (reco::CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin();
       i_calojet != caloJetCollection->end();
       ++i_calojet) {
    // remove jet-ele matching
    if (ele1)
      if (Distance(*i_calojet, *ele1) < 0.3)
        continue;
 
    if (ele2)
      if (Distance(*i_calojet, *ele2) < 0.3)
        continue;
 
    if (i_calojet->pt() < ptThrJet_)
      continue;
 
    nJet++;
 
    if (!jet1)
      jet1 = &(*i_calojet);
 
    else if (!jet2)
      jet2 = &(*i_calojet);
  }
 
  // ---------------------------
  // ---- Analyze Trigger Event
  // ---------------------------
 
  // check what is in the menu
  edm::Handle<edm::TriggerResults> hltresults;
  e.getByToken(triggerResults_, hltresults);
 
  if (!hltresults.isValid()) {
    edm::LogError("DQMClientExample") << "invalid collection: TriggerResults"
                                      << "\n";
    return;
  }
 
  bool hasFired = false;
  const edm::TriggerNames &trigNames = e.triggerNames(*hltresults);
  unsigned int numTriggers = trigNames.size();
 
  for (unsigned int hltIndex = 0; hltIndex < numTriggers; ++hltIndex) {
    if (trigNames.triggerName(hltIndex) == triggerPath_ && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex))
      hasFired = true;
  }
 
  // access the trigger event
  edm::Handle<trigger::TriggerEvent> triggerEvent;
  e.getByToken(triggerEvent_, triggerEvent);
  if (triggerEvent.failedToGet()) {
    edm::LogError("DQMClientExample") << "invalid collection: TriggerEvent"
                                      << "\n";
    return;
  }
 
  reco::Particle *ele1_HLT = nullptr;
  int nEle_HLT = 0;
 
  size_t filterIndex = triggerEvent->filterIndex(triggerFilter_);
  trigger::TriggerObjectCollection triggerObjects = triggerEvent->getObjects();
  if (!(filterIndex >= triggerEvent->sizeFilters())) {
    const trigger::Keys &keys = triggerEvent->filterKeys(filterIndex);
    std::vector<reco::Particle> triggeredEle;
 
    for (size_t j = 0; j < keys.size(); ++j) {
      trigger::TriggerObject foundObject = triggerObjects[keys[j]];
      if (abs(foundObject.particle().pdgId()) != 11)
        continue;  // make sure that it is an electron
 
      triggeredEle.push_back(foundObject.particle());
      ++nEle_HLT;
    }
 
    if (!triggeredEle.empty())
      ele1_HLT = &(triggeredEle.at(0));
  }
 
  //-------------------------------
  //--- Fill the histos
  //-------------------------------
 
  // vertex
  h_vertex_number->Fill(vertex_number);
 
  // met
  h_pfMet->Fill(pfMETCollection->begin()->et());
 
  // multiplicities
  h_eMultiplicity->Fill(nEle);
  h_jMultiplicity->Fill(nJet);
  h_eMultiplicity_HLT->Fill(nEle_HLT);
 
  // leading not matched
  if (ele1) {
    h_ePt_leading->Fill(ele1->pt());
    h_eEta_leading->Fill(ele1->eta());
    h_ePhi_leading->Fill(ele1->phi());
  }
  if (ele1_HLT) {
    h_ePt_leading_HLT->Fill(ele1_HLT->pt());
    h_eEta_leading_HLT->Fill(ele1_HLT->eta());
    h_ePhi_leading_HLT->Fill(ele1_HLT->phi());
  }
  // leading Jet
  if (jet1) {
    h_jPt_leading->Fill(jet1->pt());
    h_jEta_leading->Fill(jet1->eta());
    h_jPhi_leading->Fill(jet1->phi());
  }
 
  // fill only when the trigger candidate mathes with the reco one
  if (ele1 && ele1_HLT && deltaR(*ele1_HLT, *ele1) < 0.3 && hasFired == true) {
    h_ePt_leading_matched->Fill(ele1->pt());
    h_eEta_leading_matched->Fill(ele1->eta());
    h_ePhi_leading_matched->Fill(ele1->phi());
 
    h_ePt_leading_HLT_matched->Fill(ele1_HLT->pt());
    h_eEta_leading_HLT_matched->Fill(ele1_HLT->eta());
    h_ePhi_leading_HLT_matched->Fill(ele1_HLT->phi());
 
    h_ePt_diff->Fill(ele1->pt() - ele1_HLT->pt());
  }
}

References funct::abs(), pdwgLeptonRecoSkim_cfi::caloJetCollection, PbPb_ZMuSkimMuonDPG_cff::deltaR, Distance(), MillePedeFileConverter_cfg::e, pdwgLeptonRecoSkim_cfi::electronCollection, reco::Particle::eta(), reco::LeafCandidate::eta(), dqm::impl::MonitorElement::Fill(), h_eEta_leading, h_eEta_leading_HLT, h_eEta_leading_HLT_matched, h_eEta_leading_matched, h_eMultiplicity, h_eMultiplicity_HLT, h_ePhi_leading, h_ePhi_leading_HLT, h_ePhi_leading_HLT_matched, h_ePhi_leading_matched, h_ePt_diff, h_ePt_leading, h_ePt_leading_HLT, h_ePt_leading_HLT_matched, h_ePt_leading_matched, h_jEta_leading, h_jMultiplicity, h_jPhi_leading, h_jPt_leading, h_pfMet, h_vertex_number, HLTBitAnalyser_cfi::hltresults, edm::HandleBase::isValid(), dqmiolumiharvest::j, relativeConstraints::keys, MediumEle(), trigger::TriggerObject::particle(), reco::Particle::pdgId(), B2GDQM_cfi::pfMETCollection, reco::Particle::phi(), reco::LeafCandidate::phi(), reco::Particle::pt(), reco::LeafCandidate::pt(), ptThrJet_, ptThrL1_, ptThrL2_, PVPoint_, theCaloJetCollection_, theElectronCollection_, thePfMETCollection_, thePVCollection_, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEvent_, triggerFilter_, triggerMatchMonitor_cfi::triggerObjects, triggerPath_, triggerResults_, trigNames, and findQualityFiles::v.

bookHistograms()

void DQMExample_Step1::bookHistograms ( DQMStore::IBooker &  ibooker_, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 62 of file DQMExample_Step1.cc.

                                                                                                      {
  edm::LogInfo("DQMExample_Step1") << "DQMExample_Step1::bookHistograms" << std::endl;
 
  // book at beginRun
  bookHistos(ibooker_);
}

References bookHistos().

bookHistos()

void DQMExample_Step1::bookHistos ( DQMStore::IBooker &  ibooker_ )

private

Definition at line 289 of file DQMExample_Step1.cc.

                                                           {
  ibooker_.cd();
  ibooker_.setCurrentFolder("Physics/TopTest");
 
  h_vertex_number = ibooker_.book1D("Vertex_number", "Number of event vertices in collection", 40, -0.5, 39.5);
 
  h_pfMet = ibooker_.book1D("pfMet", "Pf Missing E_{T}; GeV", 20, 0.0, 100);
 
  h_eMultiplicity = ibooker_.book1D("NElectrons", "# of electrons per event", 10, 0., 10.);
  h_ePt_leading_matched = ibooker_.book1D("ElePt_leading_matched", "Pt of leading electron", 50, 0., 100.);
  h_eEta_leading_matched = ibooker_.book1D("EleEta_leading_matched", "Eta of leading electron", 50, -5., 5.);
  h_ePhi_leading_matched = ibooker_.book1D("ElePhi_leading_matched", "Phi of leading electron", 50, -3.5, 3.5);
 
  h_ePt_leading = ibooker_.book1D("ElePt_leading", "Pt of leading electron", 50, 0., 100.);
  h_eEta_leading = ibooker_.book1D("EleEta_leading", "Eta of leading electron", 50, -5., 5.);
  h_ePhi_leading = ibooker_.book1D("ElePhi_leading", "Phi of leading electron", 50, -3.5, 3.5);
 
  h_jMultiplicity = ibooker_.book1D("NJets", "# of electrons per event", 10, 0., 10.);
  h_jPt_leading = ibooker_.book1D("JetPt_leading", "Pt of leading Jet", 150, 0., 300.);
  h_jEta_leading = ibooker_.book1D("JetEta_leading", "Eta of leading Jet", 50, -5., 5.);
  h_jPhi_leading = ibooker_.book1D("JetPhi_leading", "Phi of leading Jet", 50, -3.5, 3.5);
 
  h_eMultiplicity_HLT = ibooker_.book1D("NElectrons_HLT", "# of electrons per event @HLT", 10, 0., 10.);
  h_ePt_leading_HLT = ibooker_.book1D("ElePt_leading_HLT", "Pt of leading electron @HLT", 50, 0., 100.);
  h_eEta_leading_HLT = ibooker_.book1D("EleEta_leading_HLT", "Eta of leading electron @HLT", 50, -5., 5.);
  h_ePhi_leading_HLT = ibooker_.book1D("ElePhi_leading_HLT", "Phi of leading electron @HLT", 50, -3.5, 3.5);
 
  h_ePt_leading_HLT_matched = ibooker_.book1D("ElePt_leading_HLT_matched", "Pt of leading electron @HLT", 50, 0., 100.);
  h_eEta_leading_HLT_matched =
      ibooker_.book1D("EleEta_leading_HLT_matched", "Eta of leading electron @HLT", 50, -5., 5.);
  h_ePhi_leading_HLT_matched =
      ibooker_.book1D("ElePhi_leading_HLT_matched", "Phi of leading electron @HLT", 50, -3.5, 3.5);
 
  h_ePt_diff = ibooker_.book1D("ElePt_diff_matched", "pT(RECO) - pT(HLT) for mathed candidates", 100, -10, 10.);
 
  ibooker_.cd();
}

References dqm::implementation::IBooker::book1D(), dqm::implementation::NavigatorBase::cd(), h_eEta_leading, h_eEta_leading_HLT, h_eEta_leading_HLT_matched, h_eEta_leading_matched, h_eMultiplicity, h_eMultiplicity_HLT, h_ePhi_leading, h_ePhi_leading_HLT, h_ePhi_leading_HLT_matched, h_ePhi_leading_matched, h_ePt_diff, h_ePt_leading, h_ePt_leading_HLT, h_ePt_leading_HLT_matched, h_ePt_leading_matched, h_jEta_leading, h_jMultiplicity, h_jPhi_leading, h_jPt_leading, h_pfMet, h_vertex_number, and dqm::implementation::NavigatorBase::setCurrentFolder().

Referenced by bookHistograms().

calcDeltaPhi()

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

private

Definition at line 338 of file DQMExample_Step1.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.

Distance()

double DQMExample_Step1::Distance ( const reco::Candidate &  c1, const reco::Candidate &  c2  )

private

Definition at line 331 of file DQMExample_Step1.cc.

{ return deltaR(c1, c2); }

References alignmentValidation::c1, and PbPb_ZMuSkimMuonDPG_cff::deltaR.

Referenced by analyze().

DistancePhi()

double DQMExample_Step1::DistancePhi ( const reco::Candidate &  c1, const reco::Candidate &  c2  )

private

Definition at line 333 of file DQMExample_Step1.cc.

                                                                                     {
  return deltaPhi(c1.p4().phi(), c2.p4().phi());
}

References alignmentValidation::c1, SiPixelRawToDigiRegional_cfi::deltaPhi, and reco::Candidate::p4().

dqmBeginRun()

void DQMExample_Step1::dqmBeginRun ( edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 55 of file DQMExample_Step1.cc.

                                                                        {
  edm::LogInfo("DQMExample_Step1") << "DQMExample_Step1::beginRun" << std::endl;
}

MediumEle()

bool DQMExample_Step1::MediumEle ( const edm::Event &  iEvent, const edm::EventSetup &  iESetup, const reco::GsfElectron &  electron  )

private

Definition at line 352 of file DQMExample_Step1.cc.

                                                                  {
  //********* CONVERSION TOOLS
  edm::Handle<reco::ConversionCollection> conversions_h;
  iEvent.getByToken(theConversionCollection_, conversions_h);
 
  bool isMediumEle = false;
 
  float pt = electron.pt();
  float eta = electron.eta();
 
  int isEB = electron.isEB();
  float sigmaIetaIeta = electron.sigmaIetaIeta();
  float DetaIn = electron.deltaEtaSuperClusterTrackAtVtx();
  float DphiIn = electron.deltaPhiSuperClusterTrackAtVtx();
  float HOverE = electron.hadronicOverEm();
  float ooemoop = (1.0 / electron.ecalEnergy() - electron.eSuperClusterOverP() / electron.ecalEnergy());
 
  int mishits = electron.gsfTrack()->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
  int nAmbiguousGsfTracks = electron.ambiguousGsfTracksSize();
 
  reco::GsfTrackRef eleTrack = electron.gsfTrack();
  float dxy = eleTrack->dxy(PVPoint_);
  float dz = eleTrack->dz(PVPoint_);
 
  edm::Handle<reco::BeamSpot> BSHandle;
  iEvent.getByToken(theBSCollection_, BSHandle);
  const reco::BeamSpot BS = *BSHandle;
 
  bool isConverted = ConversionTools::hasMatchedConversion(electron, *conversions_h, BS.position());
 
  // default
  if ((pt > 12.) && (fabs(eta) < 2.5) &&
      (((isEB == 1) && (fabs(DetaIn) < 0.004)) || ((isEB == 0) && (fabs(DetaIn) < 0.007))) &&
      (((isEB == 1) && (fabs(DphiIn) < 0.060)) || ((isEB == 0) && (fabs(DphiIn) < 0.030))) &&
      (((isEB == 1) && (sigmaIetaIeta < 0.010)) || ((isEB == 0) && (sigmaIetaIeta < 0.030))) &&
      (((isEB == 1) && (HOverE < 0.120)) || ((isEB == 0) && (HOverE < 0.100))) &&
      (((isEB == 1) && (fabs(ooemoop) < 0.050)) || ((isEB == 0) && (fabs(ooemoop) < 0.050))) &&
      (((isEB == 1) && (fabs(dxy) < 0.020)) || ((isEB == 0) && (fabs(dxy) < 0.020))) &&
      (((isEB == 1) && (fabs(dz) < 0.100)) || ((isEB == 0) && (fabs(dz) < 0.100))) &&
      (((isEB == 1) && (!isConverted)) || ((isEB == 0) && (!isConverted))) && (mishits == 0) &&
      (nAmbiguousGsfTracks == 0))
    isMediumEle = true;
 
  return isMediumEle;
}

References PVValHelper::dxy, PVValHelper::dz, HPSPFTauProducerPuppi_cfi::electron, PVValHelper::eta, ConversionTools::hasMatchedConversion(), iEvent, reco::HitPattern::MISSING_INNER_HITS, reco::BeamSpot::position(), DiDispStaMuonMonitor_cfi::pt, PVPoint_, theBSCollection_, and theConversionCollection_.

Referenced by analyze().

Member Data Documentation

h_eEta_leading

MonitorElement* DQMExample_Step1::h_eEta_leading

private

Definition at line 96 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_eEta_leading_HLT

MonitorElement* DQMExample_Step1::h_eEta_leading_HLT

private

Definition at line 104 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_eEta_leading_HLT_matched

MonitorElement* DQMExample_Step1::h_eEta_leading_HLT_matched

private

Definition at line 107 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_eEta_leading_matched

MonitorElement* DQMExample_Step1::h_eEta_leading_matched

private

Definition at line 99 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_eMultiplicity

MonitorElement* DQMExample_Step1::h_eMultiplicity

private

Definition at line 94 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_eMultiplicity_HLT

MonitorElement* DQMExample_Step1::h_eMultiplicity_HLT

private

Definition at line 102 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePhi_leading

MonitorElement* DQMExample_Step1::h_ePhi_leading

private

Definition at line 97 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePhi_leading_HLT

MonitorElement* DQMExample_Step1::h_ePhi_leading_HLT

private

Definition at line 105 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePhi_leading_HLT_matched

MonitorElement* DQMExample_Step1::h_ePhi_leading_HLT_matched

private

Definition at line 108 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePhi_leading_matched

MonitorElement* DQMExample_Step1::h_ePhi_leading_matched

private

Definition at line 100 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePt_diff

MonitorElement* DQMExample_Step1::h_ePt_diff

private

Definition at line 115 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePt_leading

MonitorElement* DQMExample_Step1::h_ePt_leading

private

Definition at line 95 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePt_leading_HLT

MonitorElement* DQMExample_Step1::h_ePt_leading_HLT

private

Definition at line 103 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePt_leading_HLT_matched

MonitorElement* DQMExample_Step1::h_ePt_leading_HLT_matched

private

Definition at line 106 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_ePt_leading_matched

MonitorElement* DQMExample_Step1::h_ePt_leading_matched

private

Definition at line 98 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_jEta_leading

MonitorElement* DQMExample_Step1::h_jEta_leading

private

Definition at line 112 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_jMultiplicity

MonitorElement* DQMExample_Step1::h_jMultiplicity

private

Definition at line 110 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_jPhi_leading

MonitorElement* DQMExample_Step1::h_jPhi_leading

private

Definition at line 113 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_jPt_leading

MonitorElement* DQMExample_Step1::h_jPt_leading

private

Definition at line 111 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_pfMet

MonitorElement* DQMExample_Step1::h_pfMet

private

Definition at line 92 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

h_vertex_number

MonitorElement* DQMExample_Step1::h_vertex_number

private

Definition at line 90 of file DQMExample_Step1.h.

Referenced by analyze(), and bookHistos().

nBJets

int DQMExample_Step1::nBJets

private

Definition at line 87 of file DQMExample_Step1.h.

nElectrons

int DQMExample_Step1::nElectrons

private

Definition at line 86 of file DQMExample_Step1.h.

ptThrJet_

double DQMExample_Step1::ptThrJet_

private

Definition at line 83 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

ptThrL1_

double DQMExample_Step1::ptThrL1_

private

Definition at line 81 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

ptThrL2_

double DQMExample_Step1::ptThrL2_

private

Definition at line 82 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

ptThrMet_

double DQMExample_Step1::ptThrMet_

private

Definition at line 84 of file DQMExample_Step1.h.

Referenced by DQMExample_Step1().

PVPoint_

math::XYZPoint DQMExample_Step1::PVPoint_

private

Definition at line 67 of file DQMExample_Step1.h.

Referenced by analyze(), and MediumEle().

theBSCollection_

edm::EDGetTokenT<reco::BeamSpot> DQMExample_Step1::theBSCollection_

private

Definition at line 75 of file DQMExample_Step1.h.

Referenced by DQMExample_Step1(), and MediumEle().

theCaloJetCollection_

edm::EDGetTokenT<reco::CaloJetCollection> DQMExample_Step1::theCaloJetCollection_

private

Definition at line 72 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

theConversionCollection_

edm::EDGetTokenT<reco::ConversionCollection> DQMExample_Step1::theConversionCollection_

private

Definition at line 71 of file DQMExample_Step1.h.

Referenced by DQMExample_Step1(), and MediumEle().

theElectronCollection_

edm::EDGetTokenT<reco::GsfElectronCollection> DQMExample_Step1::theElectronCollection_

private

Definition at line 70 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

thePfMETCollection_

edm::EDGetTokenT<reco::PFMETCollection> DQMExample_Step1::thePfMETCollection_

private

Definition at line 73 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

thePVCollection_

edm::EDGetTokenT<reco::VertexCollection> DQMExample_Step1::thePVCollection_

private

Definition at line 74 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

triggerEvent_

edm::EDGetTokenT<trigger::TriggerEvent> DQMExample_Step1::triggerEvent_

private

Definition at line 76 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

triggerFilter_

edm::InputTag DQMExample_Step1::triggerFilter_

private

Definition at line 78 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

triggerPath_

std::string DQMExample_Step1::triggerPath_

private

Definition at line 79 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().

triggerResults_

edm::EDGetTokenT<edm::TriggerResults> DQMExample_Step1::triggerResults_

private

Definition at line 77 of file DQMExample_Step1.h.

Referenced by analyze(), and DQMExample_Step1().