EwkElecDQM Class Reference

#include <EwkElecDQM.h>

Inheritance diagram for EwkElecDQM:

Public Member Functions
	EwkElecDQM (const edm::ParameterSet &)
Public Member Functions inherited from DQMOneEDAnalyzer<>
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
	DQMOneEDAnalyzer ()
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
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
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
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
void	dqmEndRun (const edm::Run &, const edm::EventSetup &) override
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotTag_
MonitorElement *	detainbarrel_after_
MonitorElement *	detainbarrel_before_
double	detainCutBarrel_
double	detainCutEndcap_
MonitorElement *	detainendcap_after_
MonitorElement *	detainendcap_before_
MonitorElement *	ecalisobarrel_after_
MonitorElement *	ecalisobarrel_before_
double	ecalIsoCutBarrel_
double	ecalIsoCutEndcap_
MonitorElement *	ecalisoendcap_after_
MonitorElement *	ecalisoendcap_before_
double	eJetMin_
edm::EDGetTokenT< edm::View< reco::GsfElectron > >	elecTag_
const std::vector< std::string >	elecTrig_
MonitorElement *	eta_after_
MonitorElement *	eta_before_
double	etaCut_
MonitorElement *	hcalisobarrel_after_
MonitorElement *	hcalisobarrel_before_
double	hcalIsoCutBarrel_
double	hcalIsoCutEndcap_
MonitorElement *	hcalisoendcap_after_
MonitorElement *	hcalisoendcap_before_
HLTPrescaleProvider	hltPrescaleProvider_
MonitorElement *	invmass_after_
MonitorElement *	invmass_before_
MonitorElement *	invmassPU_afterZ_
MonitorElement *	invmassPU_before_
bool	isValidHltConfig_
MonitorElement *	jet_et_after_
MonitorElement *	jet_et_before_
MonitorElement *	jet_eta_after_
MonitorElement *	jet_eta_before_
edm::InputTag	jetTag_
edm::EDGetTokenT< edm::View< reco::Jet > >	jetToken_
MonitorElement *	met_after_
MonitorElement *	met_before_
bool	metIncludesMuons_
double	metMax_
double	metMin_
edm::InputTag	metTag_
edm::EDGetTokenT< edm::View< reco::MET > >	metToken_
MonitorElement *	mt_after_
MonitorElement *	mt_before_
double	mtMax_
double	mtMin_
unsigned int	nall
unsigned int	neid
MonitorElement *	nelectrons_after_
MonitorElement *	nelectrons_before_
unsigned int	nGoodElectrons
unsigned int	niso
int	nJetMax_
MonitorElement *	njets_after_
MonitorElement *	njets_before_
MonitorElement *	npvs_afterZ_
MonitorElement *	npvs_before_
unsigned int	nrec
unsigned int	nsel
MonitorElement *	pt_after_
MonitorElement *	pt_before_
double	ptCut_
unsigned int	PUBinCount_
unsigned int	PUMax_
MonitorElement *	sieiebarrel_after_
MonitorElement *	sieiebarrel_before_
double	sieieCutBarrel_
double	sieieCutEndcap_
MonitorElement *	sieieendcap_after_
MonitorElement *	sieieendcap_before_
MonitorElement *	trig_after_
MonitorElement *	trig_before_
edm::EDGetTokenT< edm::TriggerResults >	trigTag_
MonitorElement *	trkisobarrel_after_
MonitorElement *	trkisobarrel_before_
double	trkIsoCutBarrel_
double	trkIsoCutEndcap_
MonitorElement *	trkisoendcap_after_
MonitorElement *	trkisoendcap_before_
edm::EDGetTokenT< edm::View< reco::Vertex > >	vertexTag_

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer<>
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMOneEDAnalyzer<>
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

DQM offline for EWK Electrons

Definition at line 27 of file EwkElecDQM.h.

Constructor & Destructor Documentation

EwkElecDQM()

EwkElecDQM::EwkElecDQM

(

const edm::ParameterSet &

cfg

)

Definition at line 36 of file EwkElecDQM.cc.

    :  // Input collections
      metTag_(cfg.getUntrackedParameter<edm::InputTag>("METTag", edm::InputTag("met"))),
      jetTag_(cfg.getUntrackedParameter<edm::InputTag>("JetTag", edm::InputTag("sisCone5CaloJets"))),
      //      trigTag_(consumes<edm::TriggerResults>(
      //        cfg.getUntrackedParameter<edm::InputTag>(
      //         "TrigTag", edm::InputTag("TriggerResults::HLT")))),
      elecTag_(consumes<edm::View<reco::GsfElectron> >(
          cfg.getUntrackedParameter<edm::InputTag>("ElecTag", edm::InputTag("gsfElectrons")))),
      metToken_(
          consumes<edm::View<reco::MET> >(cfg.getUntrackedParameter<edm::InputTag>("METTag", edm::InputTag("met")))),
      jetToken_(consumes<edm::View<reco::Jet> >(
          cfg.getUntrackedParameter<edm::InputTag>("JetTag", edm::InputTag("sisCone5CaloJets")))),
      vertexTag_(consumes<edm::View<reco::Vertex> >(
          cfg.getUntrackedParameter<edm::InputTag>("VertexTag", edm::InputTag("offlinePrimaryVertices")))),
      beamSpotTag_(
          consumes<reco::BeamSpot>(cfg.getUntrackedParameter<edm::InputTag>("BeamSpotTag", edm::InputTag("BeamSpot")))),
 
      // Main cuts
      //      muonTrig_(cfg.getUntrackedParameter<std::string> ("MuonTrig",
      // "HLT_Mu9")),
      // elecTrig_(cfg.getUntrackedParameter<std::vector< std::string >
      // >("ElecTrig", "HLT_Ele10_SW_L1R")),
      elecTrig_(cfg.getUntrackedParameter<std::vector<std::string> >("ElecTrig")),
      //      ptCut_(cfg.getUntrackedParameter<double>("PtCut", 25.)),
      ptCut_(cfg.getUntrackedParameter<double>("PtCut", 10.)),
      //      etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.1)),
      etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.4)),
      sieieCutBarrel_(cfg.getUntrackedParameter<double>("SieieBarrel", 0.01)),
      sieieCutEndcap_(cfg.getUntrackedParameter<double>("SieieEndcap", 0.028)),
      detainCutBarrel_(cfg.getUntrackedParameter<double>("DetainBarrel", 0.0071)),
      detainCutEndcap_(cfg.getUntrackedParameter<double>("DetainEndcap", 0.0066)),
      //      isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso",
      // true)),
      //      isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso",
      // false)),
      //      isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03", 0.1)),
      ecalIsoCutBarrel_(cfg.getUntrackedParameter<double>("EcalIsoCutBarrel", 5.7)),
      ecalIsoCutEndcap_(cfg.getUntrackedParameter<double>("EcalIsoCutEndcap", 5.0)),
      hcalIsoCutBarrel_(cfg.getUntrackedParameter<double>("HcalIsoCutBarrel", 8.1)),
      hcalIsoCutEndcap_(cfg.getUntrackedParameter<double>("HcalIsoCutEndcap", 3.4)),
      trkIsoCutBarrel_(cfg.getUntrackedParameter<double>("TrkIsoCutBarrel", 7.2)),
      trkIsoCutEndcap_(cfg.getUntrackedParameter<double>("TrkIsoCutEndcap", 5.1)),
      mtMin_(cfg.getUntrackedParameter<double>("MtMin", -999999)),
      mtMax_(cfg.getUntrackedParameter<double>("MtMax", 999999.)),
      metMin_(cfg.getUntrackedParameter<double>("MetMin", -999999.)),
      metMax_(cfg.getUntrackedParameter<double>("MetMax", 999999.)),
      //      acopCut_(cfg.getUntrackedParameter<double>("AcopCut", 2.)),
 
      // Muon quality cuts
      //      dxyCut_(cfg.getUntrackedParameter<double>("DxyCut", 0.2)),
      //      normalizedChi2Cut_(cfg.getUntrackedParameter<double>("NormalizedChi2Cut",
      // 10.)),
      //      trackerHitsCut_(cfg.getUntrackedParameter<int>("TrackerHitsCut",
      // 11)),
      //      isAlsoTrackerMuon_(cfg.getUntrackedParameter<bool>("IsAlsoTrackerMuon",
      // true)),
 
      // Z rejection
      //      ptThrForZ1_(cfg.getUntrackedParameter<double>("PtThrForZ1", 20.)),
      //      ptThrForZ2_(cfg.getUntrackedParameter<double>("PtThrForZ2", 10.)),
 
      // Top rejection
      eJetMin_(cfg.getUntrackedParameter<double>("EJetMin", 999999.)),
      nJetMax_(cfg.getUntrackedParameter<int>("NJetMax", 999999)),
      PUMax_(cfg.getUntrackedParameter<unsigned int>("PUMax", 60)),
      PUBinCount_(cfg.getUntrackedParameter<unsigned int>("PUBinCount", 12)),
      hltPrescaleProvider_(cfg, consumesCollector(), *this)
//       caloJetCollection_(cfg.getUntrackedParameter<edm:InputTag>("CaloJetCollection","sisCone5CaloJets"))
{
  isValidHltConfig_ = false;
}

References isValidHltConfig_.

Member Function Documentation

analyze()

void EwkElecDQM::analyze ( const edm::Event &  ev, const edm::EventSetup &  iSet  )

overrideprotectedvirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 347 of file EwkElecDQM.cc.

                                                                {
  // Reset global event selection flags
  bool rec_sel = false;
  bool eid_sel = false;
  bool iso_sel = false;
  bool all_sel = false;
 
  // Electron collection
  Handle<View<GsfElectron> > electronCollection;
  if (!ev.getByToken(elecTag_, electronCollection)) {
    // LogWarning("") << ">>> Electron collection does not exist !!!";
    return;
  }
  unsigned int electronCollectionSize = electronCollection->size();
 
  // Beam spot
  Handle<reco::BeamSpot> beamSpotHandle;
  if (!ev.getByToken(beamSpotTag_, beamSpotHandle)) {
    // LogWarning("") << ">>> No beam spot found !!!";
    return;
  }
 
  // MET
  double met_px = 0.;
  double met_py = 0.;
  Handle<View<MET> > metCollection;
  if (!ev.getByToken(metToken_, metCollection)) {
    // LogWarning("") << ">>> MET collection does not exist !!!";
    return;
  }
 
  const MET& met = metCollection->at(0);
  met_px = met.px();
  met_py = met.py();
  //       if (!metIncludesMuons_) {
  //             for (unsigned int i=0; i<muonCollectionSize; i++) {
  //                   const Muon& mu = muonCollection->at(i);
  //                   if (!mu.isGlobalMuon()) continue;
  //                   met_px -= mu.px();
  //                   met_py -= mu.py();
  //             }
  //       }
  double met_et = sqrt(met_px * met_px + met_py * met_py);
  LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met_px << ", " << met_py << " [GeV]";
  met_before_->Fill(met_et);
 
  // Vertices in the event
  int npvCount = 0;
  Handle<View<reco::Vertex> > vertexCollection;
  if (!ev.getByToken(vertexTag_, vertexCollection)) {
    LogError("") << ">>> Vertex collection does not exist !!!";
    return;
  }
 
  for (unsigned int i = 0; i < vertexCollection->size(); i++) {
    const Vertex& vertex = vertexCollection->at(i);
    if (vertex.isValid())
      npvCount++;
  }
  npvs_before_->Fill(npvCount);
 
  // Trigger
  // Handle<TriggerResults> triggerResults;
  // if (!ev.getByToken(trigTag_, triggerResults)) {
  // LogWarning("") << ">>> TRIGGER collection does not exist !!!";
  return;
  // }
  // const edm::TriggerNames& trigNames = ev.triggerNames(*triggerResults);
  // bool trigger_fired = false;
 
  // HLTConfigProvider const&  hltConfigProvider = hltPrescaleProvider_.hltConfigProvider();
 
  /* very old code
  for (unsigned int i=0; i<triggerResults->size(); i++) {
        if (triggerResults->accept(i)) {
              LogTrace("") << "Accept by: " << i << ", Trigger: " <<
  trigNames.triggerName(i);
        }
  }
  */
 
  // the following gives error on CRAFT08 data where itrig1=19 (vector index out
  // of range)
  /*
  int itrig1 = trigNames.triggerIndex(muonTrig_);
  if (triggerResults->accept(itrig1)) trigger_fired = true;
  */
  // suggested replacement: lm250909
  /*  Fix buggy trigger logic
  for (unsigned int i=0; i<triggerResults->size(); i++)
    {
      std::string trigName = trigNames.triggerName(i);
      bool found=false;
 
//    for (unsigned int j = 0; j < elecTrig_.size(); j++)
//      {
//        if ( trigName == elecTrig_.at(j) && triggerResults->accept(i))
//      {
//        trigger_fired = true;
//      }
//      }
 
 
     for(unsigned int index=0; index<elecTrig_.size() && found==false; index++)
{
       size_t trigPath = trigName.find(elecTrig_.at(index)); // 0 if found, pos
if not
               if (trigPath==0) found=true;
          }
          if(!found) continue;
 
          bool prescaled=false;
          for (unsigned int ps= 0; ps<  hltConfigProvider.prescaleSize();
ps++){
              const unsigned int prescaleValue =
hltConfigProvider.prescaleValue(ps, trigName) ;
              if (prescaleValue != 1) prescaled =true;
          }
 
          if(triggerResults->accept(i) && !prescaled) trigger_fired=true;
 
    }
  */
 
  // get the prescale set for this event
  const int prescaleSet = hltPrescaleProvider_.prescaleSet(ev, iSet);
  if (prescaleSet == -1) {
    LogTrace("") << "Failed to determine prescaleSet\n";
    // std::cout << "Failed to determine prescaleSet. Check cmsRun GlobalTag\n";
    return;
  }
 
  // for (unsigned int i = 0;
  //    (i < triggerResults->size()) && (trigger_fired == false); i++) {
  // skip trigger, if it did not fire
  //  if (!triggerResults->accept(i)) continue;
 
  // skip trigger, if it is not on our list
  // bool found = false;
  // const std::string trigName = trigNames.triggerName(i);
  // for (unsigned int index = 0; index < elecTrig_.size() && found == false;
  //    index++) {
  // if (trigName.find(elecTrig_.at(index)) == 0) found = true;
  // }
  // if (!found) continue;
 
  // skip trigger, if it is prescaled
  // if (hltConfigProvider.prescaleValue(prescaleSet, trigName) != 1) continue;
 
  // std::cout << "found unprescaled trigger that fired: " << trigName <<
  // "\n";
  // trigger_fired = true;
  // }
 
  /*  LogTrace("") << ">>> Trigger bit: " << trigger_fired << " for one of ( ";
  for (unsigned int k = 0; k < elecTrig_.size(); k++) {
    LogTrace("") << elecTrig_.at(k) << " ";
  }
  LogTrace("") << ")";
  trig_before_->Fill(trigger_fired);
*/
  // Jet collection
  Handle<View<Jet> > jetCollection;
  if (!ev.getByToken(jetToken_, jetCollection)) {
    // LogError("") << ">>> JET collection does not exist !!!";
    return;
  }
  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;
  // need to get some electron info so jets can be cleaned of them
  for (unsigned int i = 0; i < electronCollectionSize; i++) {
    const GsfElectron& elec = electronCollection->at(i);
 
    if (i < 1) {
      electron_et = elec.pt();
      electron_eta = elec.eta();
      electron_phi = elec.phi();
    }
    if (i == 2) {
      electron2_et = elec.pt();
      electron2_eta = elec.eta();
      electron2_phi = elec.phi();
    }
  }
 
  float jet_et = -8.0;
  float jet_eta = -8.0;
  int jet_count = 0;
  float jet2_et = -9.0;
  unsigned int jetCollectionSize = jetCollection->size();
  int njets = 0;
  for (unsigned int i = 0; i < jetCollectionSize; i++) {
    const Jet& jet = jetCollection->at(i);
 
    float jet_current_et = jet.et();
    //  cout << "jet_current_et " << jet_current_et << endl;
    // 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;
 
    // if it has too low Et, throw away
    //  if (jet_current_et < eJetMin_) continue; //Keep if only want to plot
    // above jet cut
 
    if (jet.et() > eJetMin_) {
      njets++;
      jet_count++;
    }
    if (jet_current_et > jet_et) {
      jet2_et = jet_et;   // 2nd highest jet get's et from current highest
      jet_et = jet.et();  // current highest jet gets et from the new highest
      jet_eta = jet.eta();
    } else if (jet_current_et > jet2_et) {
      jet2_et = jet.et();
    }
  }
 
  // Fill After all electron cuts (or both before and after)
  if (jet_et > 10)  // don't want low energy "jets"
  {
    jet_et_before_->Fill(jet_et);
    //    jet2_et_before_  ->Fill(jet2_et);
    jet_eta_before_->Fill(jet_eta);
  }
 
  LogTrace("") << ">>> Total number of jets: " << jetCollectionSize;
  LogTrace("") << ">>> Number of jets above " << eJetMin_ << " [GeV]: " << njets;
  njets_before_->Fill(njets);
 
  // Start counting
  nall++;
 
  // Histograms per event should be done only once, so keep track of them
  //bool hlt_hist_done = false;
  // bool minv_hist_done = false;
  bool met_hist_done = false;
  //       bool nz1_hist_done = false;
  //       bool nz2_hist_done = false;
  bool njets_hist_done = false;
 
  // Central selection criteria
  // const int NFLAGS = 13; // number of individual selection criteria
  const int NFLAGS = 10;  // number of individual selection criteria
  // 0: pt cut           | rec
  // 1: eta cut          | rec
  // 2: sieie            | eid
  // 3: detain           | eid
  // 4: ecal iso         | iso
  // 5: hcal iso         | iso
  // 6: trk iso          | iso
  // 7: trigger fired    | hlt/all
  bool electron_sel[NFLAGS];
 
  // for invariant mass calculation
  // keep track of highest-pt electrons for initial (RECO) electrons
  // and "good" electrons (passing all cuts)
  // first array dimension is for first or second good electron
  // second array dimension is for relevant quantities of good electron
  //    [0]: 1 for electron found or 0 for not found (if 0, no other quantities
  // filled)
  //    [1]: mSqr
  //    [2]: E
  //    [3]: px
  //    [4]: py
  //    [5]: pz
  // inv mass = sqrt(m_1^2 + m_2^2 + 2*(E_1*E_2 - (px1*px2 + py1*py2 + pz1+pz2)
  // ) )
  double electron[2][6];
  double goodElectron[2][6];
  nGoodElectrons = 0;
  for (unsigned int i = 0; i < 2; i++) {
    for (unsigned int j = 0; j < 6; j++) {
      electron[i][j] = 0.;
      goodElectron[i][j] = 0.;
    }
  }
 
  for (unsigned int i = 0; i < electronCollectionSize; i++) {
    for (int j = 0; j < NFLAGS; ++j) {
      electron_sel[j] = false;
    }
 
    const GsfElectron& elec = electronCollection->at(i);
    // if (!mu.isGlobalMuon()) continue;
    // if (mu.globalTrack().isNull()) continue;
    // if (mu.innerTrack().isNull()) continue;
 
    LogTrace("") << "> elec: processing electron number " << i << "...";
    // reco::TrackRef gm = mu.globalTrack();
    // reco::TrackRef tk = mu.innerTrack();
    // should have stuff for electron track?
 
    if (i < 2) {
      electron[i][0] = 1.;
      electron[i][1] = elec.massSqr();
      electron[i][2] = elec.energy();
      electron[i][3] = elec.px();
      electron[i][4] = elec.py();
      electron[i][5] = elec.pz();
    }
 
    // Pt,eta cuts
    double pt = elec.pt();
    double px = elec.px();
    double py = elec.py();
    double eta = elec.eta();
    LogTrace("") << "\t... pt, eta: " << pt << " [GeV], " << eta;
    ;
    if (pt > ptCut_)
      electron_sel[0] = true;
    if (fabs(eta) < etaCut_)
      electron_sel[1] = true;
 
    bool isBarrel = false;
    bool isEndcap = false;
    if (eta < 1.4442 && eta > -1.4442) {
      isBarrel = true;
    } else if ((eta > 1.56 && eta < 2.4) || (eta < -1.56 && eta > -2.4)) {
      isEndcap = true;
    }
 
    //             // d0, chi2, nhits quality cuts
    //             double dxy = tk->dxy(beamSpotHandle->position());
    //             double normalizedChi2 = gm->normalizedChi2();
    //             double trackerHits = tk->numberOfValidHits();
    //             LogTrace("") << "\t... dxy, normalizedChi2, trackerHits,
    // isTrackerMuon?: " << dxy << " [cm], " << normalizedChi2 << ", " <<
    // trackerHits << ", " << mu.isTrackerMuon();
    //             if (fabs(dxy)<dxyCut_) muon_sel[2] = true;
    //             if (normalizedChi2<normalizedChi2Cut_) muon_sel[3] = true;
    //             if (trackerHits>=trackerHitsCut_) muon_sel[4] = true;
    //             if (mu.isTrackerMuon()) muon_sel[5] = true;
 
    pt_before_->Fill(pt);
    eta_before_->Fill(eta);
    //             dxy_before_->Fill(dxy);
    //             chi2_before_->Fill(normalizedChi2);
    //             nhits_before_->Fill(trackerHits);
    //             tkmu_before_->Fill(mu.isTrackerMuon());
 
    // Electron ID cuts
    double sieie = (double)elec.sigmaIetaIeta();
    double detain = (double)elec.deltaEtaSuperClusterTrackAtVtx();  // think this is detain
    if (sieie < sieieCutBarrel_ && isBarrel)
      electron_sel[2] = true;
    if (sieie < sieieCutEndcap_ && isEndcap)
      electron_sel[2] = true;
    if (detain < detainCutBarrel_ && isBarrel)
      electron_sel[3] = true;
    if (detain < detainCutEndcap_ && isEndcap)
      electron_sel[3] = true;
    if (isBarrel) {
      LogTrace("") << "\t... sieie value " << sieie << " (barrel), pass? " << electron_sel[2];
      LogTrace("") << "\t... detain value " << detain << " (barrel), pass? " << electron_sel[3];
    } else if (isEndcap) {
      LogTrace("") << "\t... sieie value " << sieie << " (endcap), pass? " << electron_sel[2];
      LogTrace("") << "\t... detain value " << detain << " (endcap), pass? " << electron_sel[2];
    }
 
    if (isBarrel) {
      sieiebarrel_before_->Fill(sieie);
      detainbarrel_before_->Fill(detain);
    } else if (isEndcap) {
      sieieendcap_before_->Fill(sieie);
      detainendcap_before_->Fill(detain);
    }
 
    // Isolation cuts
    // double isovar = mu.isolationR03().sumPt;
    double ecalisovar = elec.dr03EcalRecHitSumEt();  // picked one set!
    double hcalisovar = elec.dr03HcalTowerSumEt();   // try others if
    double trkisovar = elec.dr04TkSumPt();           // doesn't work
    // if (isCombinedIso_) {
    // isovar += mu.isolationR03().emEt;
    // isovar += mu.isolationR03().hadEt;
    //}
    // if (isRelativeIso_) isovar /= pt;
    if (ecalisovar < ecalIsoCutBarrel_ && isBarrel)
      electron_sel[4] = true;
    if (ecalisovar < ecalIsoCutEndcap_ && isEndcap)
      electron_sel[4] = true;
    if (hcalisovar < hcalIsoCutBarrel_ && isBarrel)
      electron_sel[5] = true;
    if (hcalisovar < hcalIsoCutEndcap_ && isEndcap)
      electron_sel[5] = true;
    if (trkisovar < trkIsoCutBarrel_ && isBarrel)
      electron_sel[6] = true;
    if (trkisovar < trkIsoCutEndcap_ && isEndcap)
      electron_sel[6] = true;
    if (isBarrel) {
      LogTrace("") << "\t... ecal isolation value " << ecalisovar << " (barrel), pass? " << electron_sel[4];
      LogTrace("") << "\t... hcal isolation value " << hcalisovar << " (barrel), pass? " << electron_sel[5];
      LogTrace("") << "\t... trk isolation value " << trkisovar << " (barrel), pass? " << electron_sel[6];
    } else if (isEndcap) {
      LogTrace("") << "\t... ecal isolation value " << ecalisovar << " (endcap), pass? " << electron_sel[4];
      LogTrace("") << "\t... hcal isolation value " << hcalisovar << " (endcap), pass? " << electron_sel[5];
      LogTrace("") << "\t... trk isolation value " << trkisovar << " (endcap), pass? " << electron_sel[6];
    }
 
    // iso_before_->Fill(isovar);
    if (isBarrel) {
      ecalisobarrel_before_->Fill(ecalisovar);
      hcalisobarrel_before_->Fill(hcalisovar);
      trkisobarrel_before_->Fill(trkisovar);
    } else if (isEndcap) {
      ecalisoendcap_before_->Fill(ecalisovar);
      hcalisoendcap_before_->Fill(hcalisovar);
      trkisoendcap_before_->Fill(trkisovar);
    }
 
    // HLT
    // if (trigger_fired) electron_sel[7] = true;
 
    //             // MET/MT cuts
    double w_et = met_et + pt;
    double w_px = met_px + px;
    double w_py = met_py + py;
 
    double massT = w_et * w_et - w_px * w_px - w_py * w_py;
    massT = (massT > 0) ? sqrt(massT) : 0;
 
    LogTrace("") << "\t... W mass, W_et, W_px, W_py: " << massT << ", " << w_et << ", " << w_px << ", " << w_py
                 << " [GeV]";
    if (massT > mtMin_ && massT < mtMax_)
      electron_sel[7] = true;
    mt_before_->Fill(massT);
    if (met_et > metMin_ && met_et < metMax_)
      electron_sel[8] = true;
 
    //             // Acoplanarity cuts
    //             Geom::Phi<double> deltaphi(mu.phi()-atan2(met_py,met_px));
    //             double acop = deltaphi.value();
    //             if (acop<0) acop = - acop;
    //             acop = M_PI - acop;
    //             LogTrace("") << "\t... acoplanarity: " << acop;
    //             if (acop<acopCut_) muon_sel[10] = true;
    //             acop_before_->Fill(acop);
 
    //             // Remaining flags (from global event information)
    //             if (nmuonsForZ1<1 || nmuonsForZ2<2) muon_sel[11] = true;
    if (njets <= nJetMax_)
      electron_sel[9] = true;
 
    // Collect necessary flags "per electron"
    int flags_passed = 0;
    bool rec_sel_this = true;
    bool eid_sel_this = true;
    bool iso_sel_this = true;
    bool all_sel_this = true;
    for (int j = 0; j < NFLAGS; ++j) {
      if (electron_sel[j])
        flags_passed += 1;
      if (j < 2 && !electron_sel[j])
        rec_sel_this = false;
      if (j < 4 && !electron_sel[j])
        eid_sel_this = false;
      if (j < 7 && !electron_sel[j])
        iso_sel_this = false;
      if (!electron_sel[j])
        all_sel_this = false;
    }
 
    if (all_sel_this) {
      if (nGoodElectrons < 2) {
        goodElectron[nGoodElectrons][0] = 1.;
        goodElectron[nGoodElectrons][1] = elec.massSqr();
        goodElectron[nGoodElectrons][2] = elec.energy();
        goodElectron[nGoodElectrons][3] = elec.px();
        goodElectron[nGoodElectrons][4] = elec.py();
        goodElectron[nGoodElectrons][5] = elec.pz();
      }
      nGoodElectrons++;
    }
 
    //             // "rec" => pt,eta and quality cuts are satisfied
    //             if (rec_sel_this) rec_sel = true;
    //             // "iso" => "rec" AND "muon is isolated"
    //             if (iso_sel_this) iso_sel = true;
    //             // "hlt" => "iso" AND "event is triggered"
    //             if (hlt_sel_this) hlt_sel = true;
    //           // "all" => "met" AND "Z/top rejection cuts"
    //             if (all_sel_this) all_sel = true;
 
    // "rec" => pt,eta cuts are satisfied
    if (rec_sel_this)
      rec_sel = true;
    // "eid" => "rec" AND "electron passes ID"
    if (eid_sel_this)
      iso_sel = true;
    // "iso" => "eid" AND "electron is isolated"
    if (iso_sel_this)
      iso_sel = true;
    // "met" => "iso" AND "MET/MT"
    // "all" => "met" AND "event is triggered"
    if (all_sel_this)
      all_sel = true;
 
    // Do N-1 histograms now (and only once for global event quantities)
    if (flags_passed >= (NFLAGS - 1)) {
      if (!electron_sel[0] || flags_passed == NFLAGS) {
        pt_after_->Fill(pt);
      }
      if (!electron_sel[1] || flags_passed == NFLAGS) {
        eta_after_->Fill(eta);
      }
      if (!electron_sel[2] || flags_passed == NFLAGS) {
        if (isBarrel) {
          sieiebarrel_after_->Fill(sieie);
        } else if (isEndcap) {
          sieieendcap_after_->Fill(sieie);
        }
      }
      if (!electron_sel[3] || flags_passed == NFLAGS) {
        if (isBarrel) {
          detainbarrel_after_->Fill(detain);
        } else if (isEndcap) {
          detainendcap_after_->Fill(detain);
        }
      }
      if (!electron_sel[4] || flags_passed == NFLAGS) {
        if (isBarrel) {
          ecalisobarrel_after_->Fill(ecalisovar);
        } else if (isEndcap) {
          ecalisoendcap_after_->Fill(ecalisovar);
        }
      }
      if (!electron_sel[5] || flags_passed == NFLAGS) {
        if (isBarrel) {
          hcalisobarrel_after_->Fill(hcalisovar);
        } else if (isEndcap) {
          hcalisoendcap_after_->Fill(hcalisovar);
        }
      }
      if (!electron_sel[6] || flags_passed == NFLAGS) {
        if (isBarrel) {
          trkisobarrel_after_->Fill(trkisovar);
        } else if (isEndcap) {
          trkisoendcap_after_->Fill(trkisovar);
        }
      }
      //        if (!electron_sel[3] || flags_passed==NFLAGS)
      //          {
      //            detain_after_->Fill(detain);
      //          }
      //        if (!electron_sel[4] || flags_passed==NFLAGS)
      //          {
      //            ecaliso_after_->Fill(trackerHits);
      //          }
      //        if (!electron_sel[5] || flags_passed==NFLAGS)
      //          {
      //            tkelectr_after_->Fill(electr.isTrackerElectron());
      //          }
      //        if (!electron_sel[6] || flags_passed==NFLAGS)
      //          {
      //            iso_after_->Fill(isovar);
      //          }
      /*  if (!electron_sel[7] || flags_passed == NFLAGS) {
        if (!hlt_hist_done) {
          trig_after_->Fill(trigger_fired);
        }
      }*/
      //hlt_hist_done = true;
      if (!electron_sel[7] || flags_passed == NFLAGS) {
        mt_after_->Fill(massT);
      }
      if (!electron_sel[8] || flags_passed == NFLAGS) {
        if (!met_hist_done) {
          met_after_->Fill(met_et);
        }
      }
      met_hist_done = true;
      //                   if (!muon_sel[10] || flags_passed==NFLAGS)
      //                         acop_after_->Fill(acop);
      //                   if (!muon_sel[11] || flags_passed==NFLAGS)
      //                         if (!nz1_hist_done)
      // nz1_after_->Fill(nmuonsForZ1);
      //                         nz1_hist_done = true;
      //                   if (!muon_sel[11] || flags_passed==NFLAGS)
      //                         if (!nz2_hist_done)
      // nz2_after_->Fill(nmuonsForZ2);
      //                         nz2_hist_done = true;
      if (!electron_sel[9] || flags_passed == NFLAGS) {
        if (!njets_hist_done) {
          njets_after_->Fill(njets);
          if (jet_et > 10)  // don't want low energy "jets"
          {
            jet_et_after_->Fill(jet_et);
            jet_eta_after_->Fill(jet_eta);
          }
        }
      }
      njets_hist_done = true;
 
    }  // end N-1 histos block
 
  }  // end loop through electrons
 
  // inv mass = sqrt(m_1^2 + m_2^2 + 2*(E_1*E_2 - (px1*px2 + py1*py2 + pz1+pz2)
  // ) )
  double invMass = 0;
 
  nelectrons_before_->Fill(electronCollectionSize);
  if (electronCollectionSize > 1) {
    invMass =
        sqrt(electron[0][1] + electron[1][1] +
             2 * (electron[0][2] * electron[1][2] - (electron[0][3] * electron[1][3] + electron[0][4] * electron[1][4] +
                                                     electron[0][5] * electron[1][5])));
    invmass_before_->Fill(invMass);
    invmassPU_before_->Fill(invMass, npvCount);
  }
 
  nelectrons_after_->Fill(nGoodElectrons);
  if (nGoodElectrons > 1) {
    invMass = sqrt(goodElectron[0][1] + goodElectron[1][1] +
                   2 * (goodElectron[0][2] * goodElectron[1][2] -
                        (goodElectron[0][3] * goodElectron[1][3] + goodElectron[0][4] * goodElectron[1][4] +
                         goodElectron[0][5] * goodElectron[1][5])));
    invmass_after_->Fill(invMass);
    invmassPU_afterZ_->Fill(invMass, npvCount);
    npvs_afterZ_->Fill(npvCount);
  }
 
  // Collect final flags
  if (rec_sel)
    nrec++;
  if (eid_sel)
    neid++;
  if (iso_sel)
    niso++;
  //      if (hlt_sel) nhlt++;
  //      if (met_sel) nmet++;
 
  if (all_sel) {
    nsel++;
    LogTrace("") << ">>>> Event ACCEPTED";
  } else {
    LogTrace("") << ">>>> Event REJECTED";
  }
 
  return;
}

References beamSpotTag_, calcDeltaPhi(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), detainbarrel_after_, detainbarrel_before_, detainCutBarrel_, detainCutEndcap_, detainendcap_after_, detainendcap_before_, reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::dr03HcalTowerSumEt(), reco::GsfElectron::dr04TkSumPt(), ecalisobarrel_after_, ecalisobarrel_before_, ecalIsoCutBarrel_, ecalIsoCutEndcap_, ecalisoendcap_after_, ecalisoendcap_before_, eJetMin_, elecTag_, HPSPFTauProducerPuppi_cfi::electron, pdwgLeptonRecoSkim_cfi::electronCollection, reco::LeafCandidate::energy(), PVValHelper::eta, reco::LeafCandidate::eta(), eta_after_, eta_before_, etaCut_, ev, dqm::impl::MonitorElement::Fill(), hcalisobarrel_after_, hcalisobarrel_before_, hcalIsoCutBarrel_, hcalIsoCutEndcap_, hcalisoendcap_after_, hcalisoendcap_before_, hltPrescaleProvider_, mps_fire::i, invmass_after_, invmass_before_, invmassPU_afterZ_, invmassPU_before_, PixelPluginsPhase0_cfi::isBarrel, GeomDetEnumerators::isEndcap(), dqmiolumiharvest::j, metsig::jet, jet_et_after_, jet_et_before_, jet_eta_after_, jet_eta_before_, jetfilter_cfi::jetCollection, jetToken_, LogTrace, reco::LeafCandidate::massSqr(), BTaggingMonitor_cfi::met, met_after_, met_before_, susyDQM_cfi::metCollection, metMax_, metMin_, metToken_, mt_after_, mt_before_, mtMax_, mtMin_, nall, neid, nelectrons_after_, nelectrons_before_, nGoodElectrons, niso, nJetMax_, BTaggingMonitoring_cff::njets, njets_after_, njets_before_, npvs_afterZ_, npvs_before_, nrec, nsel, reco::LeafCandidate::phi(), HLTPrescaleProvider::prescaleSet(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), pt_after_, pt_before_, ptCut_, multPhiCorr_741_25nsDY_cfi::px, reco::LeafCandidate::px(), multPhiCorr_741_25nsDY_cfi::py, reco::LeafCandidate::py(), reco::LeafCandidate::pz(), electrons_cff::sieie, sieiebarrel_after_, sieiebarrel_before_, sieieCutBarrel_, sieieCutEndcap_, sieieendcap_after_, sieieendcap_before_, reco::GsfElectron::sigmaIetaIeta(), mathSSE::sqrt(), trkisobarrel_after_, trkisobarrel_before_, trkIsoCutBarrel_, trkIsoCutEndcap_, trkisoendcap_after_, trkisoendcap_before_, bphysicsOniaDQM_cfi::vertex, spclusmultinvestigator_cfi::vertexCollection, and vertexTag_.

bookHistograms()

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

overrideprotectedvirtual

1.5); // elecTrig_ is now a vector of strings!

Implements DQMOneEDAnalyzer<>.

Definition at line 128 of file EwkElecDQM.cc.

                                                                                             {
  ibooker.setCurrentFolder("Physics/EwkElecDQM");
 
  char chtitle[256] = "";
 
  pt_before_ = ibooker.book1D("PT_BEFORECUTS", "Electron transverse momentum [GeV]", 100, 0., 100.);
  pt_after_ = ibooker.book1D("PT_LASTCUT", "Electron transverse momentum [GeV]", 100, 0., 100.);
 
  eta_before_ = ibooker.book1D("ETA_BEFORECUTS", "Electron pseudo-rapidity", 50, -2.5, 2.5);
  eta_after_ = ibooker.book1D("ETA_LASTCUT", "Electron pseudo-rapidity", 50, -2.5, 2.5);
 
  sieiebarrel_before_ = ibooker.book1D("SIEIEBARREL_BEFORECUTS", "Electron #sigma_{i#etai#eta} (barrel)", 70, 0., 0.07);
  sieiebarrel_after_ = ibooker.book1D("SIEIEBARREL_LASTCUT", "Electron #sigma_{i#etai#eta} (barrel)", 70, 0., 0.07);
 
  sieieendcap_before_ = ibooker.book1D("SIEIEENDCAP_BEFORECUTS", "Electron #sigma_{i#etai#eta} (endcap)", 70, 0., 0.07);
  sieieendcap_after_ = ibooker.book1D("SIEIEENDCAP_LASTCUT", "Electron #sigma_{i#etai#eta} (endcap)", 70, 0., 0.07);
 
  detainbarrel_before_ =
      ibooker.book1D("DETAINBARREL_BEFORECUTS", "Electron #Delta#eta_{in} (barrel)", 40, -0.02, 0.02);
  detainbarrel_after_ = ibooker.book1D("DETAINBARREL_LASTCUT", "Electron #Delta#eta_{in} (barrel)", 40, -0.02, 0.02);
 
  detainendcap_before_ =
      ibooker.book1D("DETAINENDCAP_BEFORECUTS", "Electron #Delta#eta_{in} (endcap)", 40, -0.02, 0.02);
  detainendcap_after_ = ibooker.book1D("DETAINENDCAP_LASTCUT", "Electron #Delta#eta_{in} (endcap)", 40, -0.02, 0.02);
 
  ecalisobarrel_before_ = ibooker.book1D(
      "ECALISOBARREL_BEFORECUTS", "Absolute electron ECAL isolation variable (barrel) [GeV]", 50, 0., 50.);
  ecalisobarrel_after_ =
      ibooker.book1D("ECALISOBARREL_LASTCUT", "Absolute electron ECAL isolation variable (barrel) [GeV]", 50, 0., 50.);
 
  ecalisoendcap_before_ = ibooker.book1D(
      "ECALISOENDCAP_BEFORECUTS", "Absolute electron ECAL isolation variable (endcap) [GeV]", 50, 0., 50.);
  ecalisoendcap_after_ =
      ibooker.book1D("ECALISOENDCAP_LASTCUT", "Absolute electron ECAL isolation variable (endcap) [GeV]", 50, 0., 50.);
 
  hcalisobarrel_before_ = ibooker.book1D(
      "HCALISOBARREL_BEFORECUTS", "Absolute electron HCAL isolation variable (barrel) [GeV]", 50, 0., 50.);
  hcalisobarrel_after_ =
      ibooker.book1D("HCALISOBARREL_LASTCUT", "Absolute electron HCAL isolation variable (barrel) [GeV]", 50, 0., 50.);
 
  hcalisoendcap_before_ = ibooker.book1D(
      "HCALISOENDCAP_BEFORECUTS", "Absolute electron HCAL isolation variable (endcap) [GeV]", 50, 0., 50.);
  hcalisoendcap_after_ =
      ibooker.book1D("HCALISOENDCAP_LASTCUT", "Absolute electron HCAL isolation variable (endcap) [GeV]", 50, 0., 50.);
 
  trkisobarrel_before_ = ibooker.book1D(
      "TRKISOBARREL_BEFORECUTS", "Absolute electron track isolation variable (barrel) [GeV]", 50, 0., 50.);
  trkisobarrel_after_ =
      ibooker.book1D("TRKISOBARREL_LASTCUT", "Absolute electron track isolation variable (barrel) [GeV]", 50, 0., 50.);
 
  trkisoendcap_before_ = ibooker.book1D(
      "TRKISOENDCAP_BEFORECUTS", "Absolute electron track isolation variable (endcap) [GeV]", 50, 0., 50.);
  trkisoendcap_after_ =
      ibooker.book1D("TRKISOENDCAP_LASTCUT", "Absolute electron track isolation variable (endcap) [GeV]", 50, 0., 50.);
 
  // trig_before_ = ibooker.book1D("TRIG_BEFORECUTS", "Trigger response", 2, -0.5,
  //  trig_after_ = ibooker.book1D("TRIG_LASTCUT", "Trigger response", 2, -0.5, 1.5);
 
  invmass_before_ = ibooker.book1D("INVMASS_BEFORECUTS", "Di-electron invariant mass [GeV]", 100, 0., 200.);
  invmass_after_ = ibooker.book1D("INVMASS_AFTERCUTS", "Di-electron invariant mass [GeV]", 100, 0., 200.);
 
  invmassPU_before_ = ibooker.book2D("INVMASS_PU_BEFORECUTS",
                                     "Di-electron invariant mass [GeV] vs PU; mass [GeV]; PU count",
                                     100,
                                     0.,
                                     200.,
                                     PUBinCount_,
                                     -0.5,
                                     PUMax_ + 0.5);
  invmassPU_afterZ_ = ibooker.book2D("INVMASS_PU_AFTERZCUTS",
                                     "Di-electron invariant mass [GeV] vs PU; mass [GeV]; PU count",
                                     100,
                                     0.,
                                     200.,
                                     PUBinCount_,
                                     -0.5,
                                     PUMax_ + 0.5);
 
  npvs_before_ =
      ibooker.book1D("NPVs_BEFORECUTS", "Number of Valid Primary Vertices; nGoodPVs", PUMax_ + 1, -0.5, PUMax_ + 0.5);
 
  npvs_afterZ_ =
      ibooker.book1D("NPVs_AFTERZCUTS", "Number of Valid Primary Vertices; nGoodPVs", PUMax_ + 1, -0.5, PUMax_ + 0.5);
 
  nelectrons_before_ = ibooker.book1D("NELECTRONS_BEFORECUTS", "Number of electrons in event", 10, -0.5, 9.5);
  nelectrons_after_ = ibooker.book1D("NELECTRONS_AFTERCUTS", "Number of electrons in event", 10, -0.5, 9.5);
 
  snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
  mt_before_ = ibooker.book1D("MT_BEFORECUTS", chtitle, 150, 0., 300.);
  mt_after_ = ibooker.book1D("MT_LASTCUT", chtitle, 150, 0., 300.);
 
  snprintf(chtitle, 255, "Missing transverse energy (%s) [GeV]", metTag_.label().data());
  met_before_ = ibooker.book1D("MET_BEFORECUTS", chtitle, 100, 0., 200.);
  met_after_ = ibooker.book1D("MET_LASTCUT", chtitle, 100, 0., 200.);
 
  snprintf(chtitle, 255, "Number of jets (%s) above %.2f GeV", jetTag_.label().data(), eJetMin_);
  njets_before_ = ibooker.book1D("NJETS_BEFORECUTS", chtitle, 10, -0.5, 9.5);
  njets_after_ = ibooker.book1D("NJETS_LASTCUT", chtitle, 10, -0.5, 9.5);
 
  snprintf(chtitle, 255, "Jet with highest E_{T} (%s)", jetTag_.label().data());
  jet_et_before_ = ibooker.book1D("JETET1_BEFORECUTS", chtitle, 20, 0., 200.0);
  jet_et_after_ = ibooker.book1D("JETET1_AFTERCUTS", chtitle, 20, 0., 200.0);
 
  snprintf(chtitle, 255, "Eta of Jet with highest E_{T} (%s)", jetTag_.label().data());
  jet_eta_before_ = ibooker.book1D("JETETA1_BEFORECUTS", chtitle, 20, -5, 5);
  jet_eta_after_ = ibooker.book1D("JETETA1_AFTERCUTS", chtitle, 20, -5, 5);
}

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), detainbarrel_after_, detainbarrel_before_, detainendcap_after_, detainendcap_before_, ecalisobarrel_after_, ecalisobarrel_before_, ecalisoendcap_after_, ecalisoendcap_before_, eJetMin_, eta_after_, eta_before_, hcalisobarrel_after_, hcalisobarrel_before_, hcalisoendcap_after_, hcalisoendcap_before_, invmass_after_, invmass_before_, invmassPU_afterZ_, invmassPU_before_, jet_et_after_, jet_et_before_, jet_eta_after_, jet_eta_before_, jetTag_, edm::InputTag::label(), met_after_, met_before_, metTag_, mt_after_, mt_before_, nelectrons_after_, nelectrons_before_, njets_after_, njets_before_, npvs_afterZ_, npvs_before_, pt_after_, pt_before_, PUBinCount_, PUMax_, dqm::implementation::NavigatorBase::setCurrentFolder(), sieiebarrel_after_, sieiebarrel_before_, sieieendcap_after_, sieieendcap_before_, trkisobarrel_after_, trkisobarrel_before_, trkisoendcap_after_, and trkisoendcap_before_.

calcDeltaPhi()

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

protected

Definition at line 996 of file EwkElecDQM.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.

Referenced by analyze().

dqmBeginRun()

void EwkElecDQM::dqmBeginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSet  )

overrideprotectedvirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 109 of file EwkElecDQM.cc.

                                                                    {
  nall = 0;
  nsel = 0;
 
  nrec = 0;
  neid = 0;
  niso = 0;
  //       nhlt = 0;
  //       nmet = 0;
 
  // passed as parameter to HLTConfigProvider::init(), not yet used
  bool isConfigChanged = false;
  // isValidHltConfig_ could be used to short-circuit analyze() in case of
  // problems
  isValidHltConfig_ = hltPrescaleProvider_.init(iRun, iSet, "HLT", isConfigChanged);
 
  LogTrace("") << "isValidHltConfig_=" << isValidHltConfig_ << "\n";
}

References hltPrescaleProvider_, HLTPrescaleProvider::init(), isValidHltConfig_, LogTrace, nall, neid, niso, nrec, and nsel.

dqmEndRun()

void EwkElecDQM::dqmEndRun ( const edm::Run &  r, const edm::EventSetup &    )

overrideprotectedvirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 237 of file EwkElecDQM.cc.

                                                          {
  // overall
  double all = nall;
  double esel = nsel / all;
  LogVerbatim("") << "\n>>>>>> SELECTION SUMMARY BEGIN >>>>>>>>>>>>>>>";
  LogVerbatim("") << "Total number of events analyzed: " << nall << " [events]";
  LogVerbatim("") << "Total number of events selected: " << nsel << " [events]";
  LogVerbatim("") << "Overall efficiency:             "
                  << "(" << setprecision(4) << esel * 100. << " +/- " << setprecision(2)
                  << sqrt(esel * (1 - esel) / all) * 100. << ")%";
 
  double erec = nrec / all;
  double eeid = neid / all;
  double eiso = niso / all;
  //   double ehlt = nhlt/all;
  //   double emet = nmet/all;
 
  // general reconstruction step??
  double num = nrec;
  double eff = erec;
  double err = sqrt(eff * (1 - eff) / all);
  LogVerbatim("") << "Passing Pt/Eta/Quality cuts:    " << num << " [events], (" << setprecision(4) << eff * 100.
                  << " +/- " << setprecision(2) << err * 100. << ")%";
 
  // electron ID step
  num = neid;
  eff = eeid;
  err = sqrt(eff * (1 - eff) / all);
  double effstep = 0.;
  double errstep = 0.;
  if (nrec > 0)
    effstep = eeid / erec;
  if (nrec > 0)
    errstep = sqrt(effstep * (1 - effstep) / nrec);
  LogVerbatim("") << "Passing eID cuts:         " << num << " [events], (" << setprecision(4) << eff * 100. << " +/- "
                  << setprecision(2) << err * 100. << ")%, to previous step: (" << setprecision(4) << effstep * 100.
                  << " +/- " << setprecision(2) << errstep * 100. << ")%";
 
  // isolation step
  num = niso;
  eff = eiso;
  err = sqrt(eff * (1 - eff) / all);
  effstep = 0.;
  errstep = 0.;
  if (neid > 0)
    effstep = eiso / eeid;
  if (neid > 0)
    errstep = sqrt(effstep * (1 - effstep) / neid);
  LogVerbatim("") << "Passing isolation cuts:         " << num << " [events], (" << setprecision(4) << eff * 100.
                  << " +/- " << setprecision(2) << err * 100. << ")%, to previous step: (" << setprecision(4)
                  << effstep * 100. << " +/- " << setprecision(2) << errstep * 100. << ")%";
 
  //   // trigger step
  //   num = nhlt;
  //   eff = ehlt;
  //   err = sqrt(eff*(1-eff)/all);
  //   effstep = 0.;
  //   errstep = 0.;
  //   if (niso>0) effstep = ehlt/eiso;
  //   if (niso>0) errstep = sqrt(effstep*(1-effstep)/niso);
  //   LogVerbatim("") << "Passing HLT criteria:           " << num << "
  // [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) <<
  // err*100. << ")%, to previous step: (" <<  setprecision(4) << effstep*100.
  // << " +/- "<< setprecision(2) << errstep*100. <<")%";
 
  // trigger step
  num = nsel;
  eff = esel;
  err = sqrt(eff * (1 - eff) / all);
  effstep = 0.;
  errstep = 0.;
  if (niso > 0)
    effstep = esel / eiso;
  if (niso > 0)
    errstep = sqrt(effstep * (1 - effstep) / niso);
  LogVerbatim("") << "Passing HLT criteria:           " << num << " [events], (" << setprecision(4) << eff * 100.
                  << " +/- " << setprecision(2) << err * 100. << ")%, to previous step: (" << setprecision(4)
                  << effstep * 100. << " +/- " << setprecision(2) << errstep * 100. << ")%";
 
  //   // met/acoplanarity cuts
  //   num = nmet;
  //   eff = emet;
  //   err = sqrt(eff*(1-eff)/all);
  //   effstep = 0.;
  //   errstep = 0.;
  //   if (nhlt>0) effstep = emet/ehlt;
  //   if (nhlt>0) errstep = sqrt(effstep*(1-effstep)/nhlt);
  //   LogVerbatim("") << "Passing MET/acoplanarity cuts:  " << num << "
  // [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) <<
  // err*100. << ")%, to previous step: (" <<  setprecision(4) << effstep*100.
  // << " +/- "<< setprecision(2) << errstep*100. <<")%";
 
  //   // Z/top selection cuts ALSO LAST STEP so "sel" for "selection"
  //   num = nsel;
  //   eff = esel;
  //   err = sqrt(eff*(1-eff)/all);
  //   effstep = 0.;
  //   errstep = 0.;
  //   if (nmet>0) effstep = esel/emet;
  //   if (nmet>0) errstep = sqrt(effstep*(1-effstep)/nmet);
  //   LogVerbatim("") << "Passing Z/top rejection cuts:   " << num << "
  // [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) <<
  // err*100. << ")%, to previous step: (" <<  setprecision(4) << effstep*100.
  // << " +/- "<< setprecision(2) << errstep*100. <<")%";
 
  LogVerbatim("") << ">>>>>> SELECTION SUMMARY END   >>>>>>>>>>>>>>>\n";
}

References python.cmstools::all(), submitPVResolutionJobs::err, nall, neid, niso, nrec, nsel, EgammaValidation_cff::num, and mathSSE::sqrt().

Member Data Documentation

beamSpotTag_

edm::EDGetTokenT<reco::BeamSpot> EwkElecDQM::beamSpotTag_

private

Definition at line 49 of file EwkElecDQM.h.

Referenced by analyze().

detainbarrel_after_

MonitorElement* EwkElecDQM::detainbarrel_after_

private

Definition at line 119 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

detainbarrel_before_

MonitorElement* EwkElecDQM::detainbarrel_before_

private

Definition at line 118 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

detainCutBarrel_

double EwkElecDQM::detainCutBarrel_

private

Definition at line 59 of file EwkElecDQM.h.

Referenced by analyze().

detainCutEndcap_

double EwkElecDQM::detainCutEndcap_

private

Definition at line 60 of file EwkElecDQM.h.

Referenced by analyze().

detainendcap_after_

MonitorElement* EwkElecDQM::detainendcap_after_

private

Definition at line 122 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

detainendcap_before_

MonitorElement* EwkElecDQM::detainendcap_before_

private

Definition at line 121 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

ecalisobarrel_after_

MonitorElement* EwkElecDQM::ecalisobarrel_after_

private

Definition at line 137 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

ecalisobarrel_before_

MonitorElement* EwkElecDQM::ecalisobarrel_before_

private

Definition at line 136 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

ecalIsoCutBarrel_

double EwkElecDQM::ecalIsoCutBarrel_

private

Definition at line 66 of file EwkElecDQM.h.

Referenced by analyze().

ecalIsoCutEndcap_

double EwkElecDQM::ecalIsoCutEndcap_

private

Definition at line 67 of file EwkElecDQM.h.

Referenced by analyze().

ecalisoendcap_after_

MonitorElement* EwkElecDQM::ecalisoendcap_after_

private

Definition at line 140 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

ecalisoendcap_before_

MonitorElement* EwkElecDQM::ecalisoendcap_before_

private

Definition at line 139 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

eJetMin_

double EwkElecDQM::eJetMin_

private

Definition at line 86 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

elecTag_

edm::EDGetTokenT<edm::View<reco::GsfElectron> > EwkElecDQM::elecTag_

private

Definition at line 45 of file EwkElecDQM.h.

Referenced by analyze().

elecTrig_

const std::vector<std::string> EwkElecDQM::elecTrig_

private

Definition at line 53 of file EwkElecDQM.h.

eta_after_

MonitorElement* EwkElecDQM::eta_after_

private

Definition at line 110 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

eta_before_

MonitorElement* EwkElecDQM::eta_before_

private

Definition at line 109 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

etaCut_

double EwkElecDQM::etaCut_

private

Definition at line 55 of file EwkElecDQM.h.

Referenced by analyze().

hcalisobarrel_after_

MonitorElement* EwkElecDQM::hcalisobarrel_after_

private

Definition at line 143 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

hcalisobarrel_before_

MonitorElement* EwkElecDQM::hcalisobarrel_before_

private

Definition at line 142 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

hcalIsoCutBarrel_

double EwkElecDQM::hcalIsoCutBarrel_

private

Definition at line 68 of file EwkElecDQM.h.

Referenced by analyze().

hcalIsoCutEndcap_

double EwkElecDQM::hcalIsoCutEndcap_

private

Definition at line 69 of file EwkElecDQM.h.

Referenced by analyze().

hcalisoendcap_after_

MonitorElement* EwkElecDQM::hcalisoendcap_after_

private

Definition at line 146 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

hcalisoendcap_before_

MonitorElement* EwkElecDQM::hcalisoendcap_before_

private

Definition at line 145 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

hltPrescaleProvider_

HLTPrescaleProvider EwkElecDQM::hltPrescaleProvider_

private

Definition at line 93 of file EwkElecDQM.h.

Referenced by analyze(), and dqmBeginRun().

invmass_after_

MonitorElement* EwkElecDQM::invmass_after_

private

Definition at line 158 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

invmass_before_

MonitorElement* EwkElecDQM::invmass_before_

private

Definition at line 157 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

invmassPU_afterZ_

MonitorElement* EwkElecDQM::invmassPU_afterZ_

private

Definition at line 160 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

invmassPU_before_

MonitorElement* EwkElecDQM::invmassPU_before_

private

Definition at line 159 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

isValidHltConfig_

bool EwkElecDQM::isValidHltConfig_

private

Definition at line 92 of file EwkElecDQM.h.

Referenced by dqmBeginRun(), and EwkElecDQM().

jet_et_after_

MonitorElement* EwkElecDQM::jet_et_after_

private

Definition at line 187 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

jet_et_before_

MonitorElement* EwkElecDQM::jet_et_before_

private

Definition at line 186 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

jet_eta_after_

MonitorElement* EwkElecDQM::jet_eta_after_

private

Definition at line 189 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

jet_eta_before_

MonitorElement* EwkElecDQM::jet_eta_before_

private

Definition at line 188 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

jetTag_

edm::InputTag EwkElecDQM::jetTag_

private

Definition at line 43 of file EwkElecDQM.h.

Referenced by bookHistograms().

jetToken_

edm::EDGetTokenT<edm::View<reco::Jet> > EwkElecDQM::jetToken_

private

Definition at line 47 of file EwkElecDQM.h.

Referenced by analyze().

met_after_

MonitorElement* EwkElecDQM::met_after_

private

Definition at line 173 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

met_before_

MonitorElement* EwkElecDQM::met_before_

private

Definition at line 172 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

metIncludesMuons_

bool EwkElecDQM::metIncludesMuons_

private

Definition at line 50 of file EwkElecDQM.h.

metMax_

double EwkElecDQM::metMax_

private

Definition at line 75 of file EwkElecDQM.h.

Referenced by analyze().

metMin_

double EwkElecDQM::metMin_

private

Definition at line 74 of file EwkElecDQM.h.

Referenced by analyze().

metTag_

edm::InputTag EwkElecDQM::metTag_

private

Definition at line 42 of file EwkElecDQM.h.

Referenced by bookHistograms().

metToken_

edm::EDGetTokenT<edm::View<reco::MET> > EwkElecDQM::metToken_

private

Definition at line 46 of file EwkElecDQM.h.

Referenced by analyze().

mt_after_

MonitorElement* EwkElecDQM::mt_after_

private

Definition at line 170 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

mt_before_

MonitorElement* EwkElecDQM::mt_before_

private

Definition at line 169 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

mtMax_

double EwkElecDQM::mtMax_

private

Definition at line 73 of file EwkElecDQM.h.

Referenced by analyze().

mtMin_

double EwkElecDQM::mtMin_

private

Definition at line 72 of file EwkElecDQM.h.

Referenced by analyze().

nall

unsigned int EwkElecDQM::nall

private

Definition at line 95 of file EwkElecDQM.h.

Referenced by analyze(), dqmBeginRun(), and dqmEndRun().

neid

unsigned int EwkElecDQM::neid

private

Definition at line 97 of file EwkElecDQM.h.

Referenced by analyze(), dqmBeginRun(), and dqmEndRun().

nelectrons_after_

MonitorElement* EwkElecDQM::nelectrons_after_

private

Definition at line 167 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

nelectrons_before_

MonitorElement* EwkElecDQM::nelectrons_before_

private

Definition at line 166 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

nGoodElectrons

unsigned int EwkElecDQM::nGoodElectrons

private

Definition at line 104 of file EwkElecDQM.h.

Referenced by analyze().

niso

unsigned int EwkElecDQM::niso

private

Definition at line 98 of file EwkElecDQM.h.

Referenced by analyze(), dqmBeginRun(), and dqmEndRun().

nJetMax_

int EwkElecDQM::nJetMax_

private

Definition at line 87 of file EwkElecDQM.h.

Referenced by analyze().

njets_after_

MonitorElement* EwkElecDQM::njets_after_

private

Definition at line 185 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

njets_before_

MonitorElement* EwkElecDQM::njets_before_

private

Definition at line 184 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

npvs_afterZ_

MonitorElement* EwkElecDQM::npvs_afterZ_

private

Definition at line 164 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

npvs_before_

MonitorElement* EwkElecDQM::npvs_before_

private

Definition at line 162 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

nrec

unsigned int EwkElecDQM::nrec

private

Definition at line 96 of file EwkElecDQM.h.

Referenced by analyze(), dqmBeginRun(), and dqmEndRun().

nsel

unsigned int EwkElecDQM::nsel

private

Definition at line 101 of file EwkElecDQM.h.

Referenced by analyze(), dqmBeginRun(), and dqmEndRun().

pt_after_

MonitorElement* EwkElecDQM::pt_after_

private

Definition at line 107 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

pt_before_

MonitorElement* EwkElecDQM::pt_before_

private

Definition at line 106 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

ptCut_

double EwkElecDQM::ptCut_

private

Definition at line 54 of file EwkElecDQM.h.

Referenced by analyze(), and badGlobalMuonTaggerFWLite.BadGlobalMuonTagger::badMuons().

PUBinCount_

unsigned int EwkElecDQM::PUBinCount_

private

Definition at line 90 of file EwkElecDQM.h.

Referenced by bookHistograms().

PUMax_

unsigned int EwkElecDQM::PUMax_

private

Definition at line 90 of file EwkElecDQM.h.

Referenced by bookHistograms().

sieiebarrel_after_

MonitorElement* EwkElecDQM::sieiebarrel_after_

private

Definition at line 113 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

sieiebarrel_before_

MonitorElement* EwkElecDQM::sieiebarrel_before_

private

Definition at line 112 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

sieieCutBarrel_

double EwkElecDQM::sieieCutBarrel_

private

Definition at line 57 of file EwkElecDQM.h.

Referenced by analyze().

sieieCutEndcap_

double EwkElecDQM::sieieCutEndcap_

private

Definition at line 58 of file EwkElecDQM.h.

Referenced by analyze().

sieieendcap_after_

MonitorElement* EwkElecDQM::sieieendcap_after_

private

Definition at line 116 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

sieieendcap_before_

MonitorElement* EwkElecDQM::sieieendcap_before_

private

Definition at line 115 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

trig_after_

MonitorElement* EwkElecDQM::trig_after_

private

Definition at line 155 of file EwkElecDQM.h.

trig_before_

MonitorElement* EwkElecDQM::trig_before_

private

Definition at line 154 of file EwkElecDQM.h.

trigTag_

edm::EDGetTokenT<edm::TriggerResults> EwkElecDQM::trigTag_

private

Definition at line 44 of file EwkElecDQM.h.

trkisobarrel_after_

MonitorElement* EwkElecDQM::trkisobarrel_after_

private

Definition at line 149 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

trkisobarrel_before_

MonitorElement* EwkElecDQM::trkisobarrel_before_

private

Definition at line 148 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

trkIsoCutBarrel_

double EwkElecDQM::trkIsoCutBarrel_

private

Definition at line 70 of file EwkElecDQM.h.

Referenced by analyze().

trkIsoCutEndcap_

double EwkElecDQM::trkIsoCutEndcap_

private

Definition at line 71 of file EwkElecDQM.h.

Referenced by analyze().

trkisoendcap_after_

MonitorElement* EwkElecDQM::trkisoendcap_after_

private

Definition at line 152 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

trkisoendcap_before_

MonitorElement* EwkElecDQM::trkisoendcap_before_

private

Definition at line 151 of file EwkElecDQM.h.

Referenced by analyze(), and bookHistograms().

vertexTag_

edm::EDGetTokenT<edm::View<reco::Vertex> > EwkElecDQM::vertexTag_

private

Definition at line 48 of file EwkElecDQM.h.

Referenced by analyze().