#include <EwkElecDQM.h>

Inheritance diagram for EwkElecDQM:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

virtual void	beginRun (const edm::Run &, const edm::EventSetup &)

double	calcDeltaPhi (double phi1, double phi2)

virtual void	endJob ()

virtual void	endRun (const edm::Run &, const edm::EventSetup &)

	EwkElecDQM (const edm::ParameterSet &)

void	init_histograms ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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_

HLTConfigProvider	hltConfigProvider_

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_

DQMStore *	theDbe

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 edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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 ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

DQM offline for EWK Electrons

Definition at line 24 of file EwkElecDQM.h.

Constructor & Destructor Documentation

EwkElecDQM::EwkElecDQM ( const edm::ParameterSet & cfg )

Definition at line 39 of file EwkElecDQM.cc.

References init_histograms(), isValidHltConfig_, cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), and theDbe.

                                                  :
       // 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))
 
 //       caloJetCollection_(cfg.getUntrackedParameter<edm:InputTag>("CaloJetCollection","sisCone5CaloJets"))
 
 
 
 {
   isValidHltConfig_ = false;
 
   // access to dbe
   theDbe = Service<DQMStore>().operator->();
   theDbe->setCurrentFolder("Physics/EwkElecDQM");
   init_histograms();
 
 }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 376 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;
 
       /* 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=hltConfigProvider_.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 = 11; // 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[8] = true;
       mt_before_->Fill(massT);
       if (met_et>metMin_ && met_et<metMax_) electron_sel[9] = 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[10] = 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[8] || flags_passed==NFLAGS)
         {
           mt_after_->Fill(massT);
         }
           if (!electron_sel[9] || 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[10] || 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;
 
 }

void EwkElecDQM::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 140 of file EwkElecDQM.cc.

                           {
 
 }

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 121 of file EwkElecDQM.cc.

References hltConfigProvider_, HLTConfigProvider::init(), isValidHltConfig_, LogTrace, nall, neid, niso, nrec, and nsel.

                                                                  {
       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_ = hltConfigProvider_.init( iRun, iSet, "HLT", isConfigChanged );
 
      LogTrace("") << "isValidHltConfig_=" << isValidHltConfig_ << "\n";
 }

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

Definition at line 1046 of file EwkElecDQM.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi.

Referenced by analyze().

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

void EwkElecDQM::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 285 of file EwkElecDQM.cc.

285 {

286 }

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 288 of file EwkElecDQM.cc.

References cond::ecalcond::all, nall, neid, niso, nrec, nsel, pileupDistInMC::num, and mathSSE::sqrt().

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

void EwkElecDQM::init_histograms ( )

Definition at line 144 of file EwkElecDQM.cc.

References DQMStore::book1D(), DQMStore::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_, sieiebarrel_after_, sieiebarrel_before_, sieieendcap_after_, sieieendcap_before_, theDbe, trig_after_, trig_before_, trkisobarrel_after_, trkisobarrel_before_, trkisoendcap_after_, and trkisoendcap_before_.

Referenced by EwkElecDQM().

                                  {
 
     char chtitle[256] = "";
     //             pt_before_ = theDbe->book1D("PT_BEFORECUTS","Muon transverse momentum (global muon) [GeV],100,0.,100.);
     //             pt_after_ = theDbe->book1D("PT_LASTCUT","Muon transverse momentum (global muon) [GeV],100,0.,100.);
 
     pt_before_ = theDbe->book1D("PT_BEFORECUTS","Electron transverse momentum [GeV]",100,0.,100.);
     pt_after_ = theDbe->book1D("PT_LASTCUT","Electron transverse momentum [GeV]",100,0.,100.);
 
     eta_before_ = theDbe->book1D("ETA_BEFORECUTS","Electron pseudo-rapidity",50,-2.5,2.5);
     eta_after_ = theDbe->book1D("ETA_LASTCUT","Electron pseudo-rapidity",50,-2.5,2.5);
 
     sieiebarrel_before_ = theDbe->book1D("SIEIEBARREL_BEFORECUTS","Electron #sigma_{i#etai#eta} (barrel)",70,0.,0.07);
     sieiebarrel_after_ = theDbe->book1D("SIEIEBARREL_LASTCUT","Electron #sigma_{i#etai#eta} (barrel)",70,0.,0.07);
 
     sieieendcap_before_ = theDbe->book1D("SIEIEENDCAP_BEFORECUTS","Electron #sigma_{i#etai#eta} (endcap)",70,0.,0.07);
     sieieendcap_after_ = theDbe->book1D("SIEIEENDCAP_LASTCUT","Electron #sigma_{i#etai#eta} (endcap)",70,0.,0.07);
 
     detainbarrel_before_ = theDbe->book1D("DETAINBARREL_BEFORECUTS","Electron #Delta#eta_{in} (barrel)",40,-0.02,0.02);
     detainbarrel_after_ = theDbe->book1D("DETAINBARREL_LASTCUT","Electron #Delta#eta_{in} (barrel)",40,-0.02,0.02);
 
     detainendcap_before_ = theDbe->book1D("DETAINENDCAP_BEFORECUTS","Electron #Delta#eta_{in} (endcap)",40,-0.02,0.02);
     detainendcap_after_ = theDbe->book1D("DETAINENDCAP_LASTCUT","Electron #Delta#eta_{in} (endcap)",40,-0.02,0.02);
 
     //             dxy_before_ = theDbe->book1D("DXY_BEFORECUTS","Muon transverse distance to beam spot [cm]",100,-0.5,0.5);
     //             dxy_after_ = theDbe->book1D("DXY_LASTCUT","Muon transverse distance to beam spot [cm]",100,-0.5,0.5);
 
     //             chi2_before_ = theDbe->book1D("CHI2_BEFORECUTS","Normalized Chi2, inner track fit",100,0.,100.);
     //             chi2_after_ = theDbe->book1D("CHI2_LASTCUT","Normalized Chi2, inner track fit",100,0.,100.);
 
     //             nhits_before_ = theDbe->book1D("NHITS_BEFORECUTS","Number of hits, inner track",40,-0.5,39.5);
     //             nhits_after_ = theDbe->book1D("NHITS_LASTCUT","Number of hits, inner track",40,-0.5,39.5);
 
     //             tkmu_before_ = theDbe->book1D("TKMU_BEFORECUTS","Tracker-muon flag (for global muons)",2,-0.5,1.5);
     //             tkmu_after_ = theDbe->book1D("TKMU_LASTCUT","Tracker-muon flag (for global muons)",2,-0.5,1.5);
 
     //             if (isRelativeIso_) {
     //               if (isCombinedIso_) {
     //                 snprintf(chtitle, 255, "Relative (combined) isolation variable");
     //                 iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Relative (combined) isolation variable",100, 0., 1.);
     //                 iso_after_ = theDbe->book1D("ISO_LASTCUT","Relative (combined) isolation variable",100, 0., 1.);
     //               } else {
     //                 snprintf(chtitle, 255, "Relative (tracker) isolation variable");
     //                 iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Relative (tracker) isolation variable",100, 0., 1.);
     //                 iso_after_ = theDbe->book1D("ISO_LASTCUT","Relative (tracker) isolation variable",100, 0., 1.);
     //               }
     //             } else {
     //               if (isCombinedIso_) {
     //                 iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Absolute (combined) isolation variable [GeV]",100, 0., 20.);
     //                 iso_after_ = theDbe->book1D("ISO_LASTCUT","Absolute (combined) isolation variable [GeV]",100, 0., 20.);
     //               } else {
     //                 iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Absolute (tracker) isolation variable [GeV]",100, 0., 20.);
     //                 iso_after_ = theDbe->book1D("ISO_LASTCUT","Absolute (tracker) isolation variable [GeV]",100, 0., 20.);
     //               }
     //             }
 
             ecalisobarrel_before_ = theDbe->book1D("ECALISOBARREL_BEFORECUTS","Absolute electron ECAL isolation variable (barrel) [GeV]",50,0.,50.);
             ecalisobarrel_after_ = theDbe->book1D("ECALISOBARREL_LASTCUT","Absolute electron ECAL isolation variable (barrel) [GeV]",50,0.,50.);
 
             ecalisoendcap_before_ = theDbe->book1D("ECALISOENDCAP_BEFORECUTS","Absolute electron ECAL isolation variable (endcap) [GeV]",50,0.,50.);
             ecalisoendcap_after_ = theDbe->book1D("ECALISOENDCAP_LASTCUT","Absolute electron ECAL isolation variable (endcap) [GeV]",50,0.,50.);
 
             hcalisobarrel_before_ = theDbe->book1D("HCALISOBARREL_BEFORECUTS","Absolute electron HCAL isolation variable (barrel) [GeV]",50,0.,50.);
             hcalisobarrel_after_ = theDbe->book1D("HCALISOBARREL_LASTCUT","Absolute electron HCAL isolation variable (barrel) [GeV]",50,0.,50.);
 
             hcalisoendcap_before_ = theDbe->book1D("HCALISOENDCAP_BEFORECUTS","Absolute electron HCAL isolation variable (endcap) [GeV]",50,0.,50.);
             hcalisoendcap_after_ = theDbe->book1D("HCALISOENDCAP_LASTCUT","Absolute electron HCAL isolation variable (endcap) [GeV]",50,0.,50.);
 
             trkisobarrel_before_ = theDbe->book1D("TRKISOBARREL_BEFORECUTS","Absolute electron track isolation variable (barrel) [GeV]",50,0.,50.);
             trkisobarrel_after_ = theDbe->book1D("TRKISOBARREL_LASTCUT","Absolute electron track isolation variable (barrel) [GeV]",50,0.,50.);
 
             trkisoendcap_before_ = theDbe->book1D("TRKISOENDCAP_BEFORECUTS","Absolute electron track isolation variable (endcap) [GeV]",50,0.,50.);
             trkisoendcap_after_ = theDbe->book1D("TRKISOENDCAP_LASTCUT","Absolute electron track isolation variable (endcap) [GeV]",50,0.,50.);
 
 //             snprintf(chtitle, 255, "Trigger response (bit %s)", muonTrig_.data());
 //             trig_before_ = theDbe->book1D("TRIG_BEFORECUTS",chtitle,2,-0.5,1.5);
 //             trig_after_ = theDbe->book1D("TRIG_LASTCUT",chtitle,2,-0.5,1.5);
 
             //snprintf(chtitle, 255, "Trigger response (bit %s)", elecTrig_.data());
             trig_before_ = theDbe->book1D("TRIG_BEFORECUTS","Trigger response",2,-0.5,1.5); // elecTrig_ is now a vector of strings!
             trig_after_ = theDbe->book1D("TRIG_LASTCUT","Trigger response",2,-0.5,1.5);
 
             invmass_before_ = theDbe->book1D("INVMASS_BEFORECUTS","Di-electron invariant mass [GeV]",100,0.,200.);
             invmass_after_ = theDbe->book1D("INVMASS_AFTERCUTS","Di-electron invariant mass [GeV]",100,0.,200.);
 
             invmassPU_before_ = theDbe->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_ = theDbe->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_ = theDbe->book1D("NPVs_BEFORECUTS","Number of Valid Primary Vertices; nGoodPVs",PUMax_+1,-0.5,PUMax_+0.5);
         //npvs_afterW_ = theDbe->book1D("NPVs_AFTERWCUTS","Number of Valid Primary Vertices",PUMax_+1,-0.5,PUMax_+0.5);
         npvs_afterZ_ = theDbe->book1D("NPVs_AFTERZCUTS","Number of Valid Primary Vertices; nGoodPVs",PUMax_+1,-0.5,PUMax_+0.5);
 
 
         nelectrons_before_ = theDbe->book1D("NELECTRONS_BEFORECUTS","Number of electrons in event",10,-0.5,9.5);
         nelectrons_after_ = theDbe->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_ = theDbe->book1D("MT_BEFORECUTS",chtitle,150,0.,300.);
             mt_after_ = theDbe->book1D("MT_LASTCUT",chtitle,150,0.,300.);
 
 //             snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
 //             mt_before_ = theDbe->book1D("MT_BEFORECUTS",chtitle,150,0.,300.);
 //             mt_after_ = theDbe->book1D("MT_LASTCUT",chtitle,150,0.,300.);
 
         snprintf(chtitle, 255, "Missing transverse energy (%s) [GeV]", metTag_.label().data());
         met_before_ = theDbe->book1D("MET_BEFORECUTS",chtitle,100,0.,200.);
         met_after_ = theDbe->book1D("MET_LASTCUT",chtitle,100,0.,200.);
 
 //             snprintf(chtitle, 255, "MU-MET (%s) acoplanarity", metTag_.label().data());
 //             acop_before_ = theDbe->book1D("ACOP_BEFORECUTS",chtitle,50,0.,M_PI);
 //             acop_after_ = theDbe->book1D("ACOP_LASTCUT",chtitle,50,0.,M_PI);
 
 //             snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ1_);
 //             nz1_before_ = theDbe->book1D("NZ1_BEFORECUTS",chtitle,10,-0.5,9.5);
 //             nz1_after_ = theDbe->book1D("NZ1_LASTCUT",chtitle,10,-0.5,9.5);
 
 //             snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ2_);
 //             nz2_before_ = theDbe->book1D("NZ2_BEFORECUTS",chtitle,10,-0.5,9.5);
 //             nz2_after_ = theDbe->book1D("NZ2_LASTCUT",chtitle,10,-0.5,9.5);
 
              snprintf(chtitle, 255, "Number of jets (%s) above %.2f GeV", jetTag_.label().data(), eJetMin_);
              njets_before_ = theDbe->book1D("NJETS_BEFORECUTS",chtitle,10,-0.5,9.5);
              njets_after_  = theDbe->book1D("NJETS_LASTCUT",chtitle,10,-0.5,9.5);
 
          snprintf(chtitle, 255, "Jet with highest E_{T} (%s)", jetTag_.label().data());
          jet_et_before_       = theDbe->book1D("JETET1_BEFORECUTS",chtitle, 20, 0., 200.0);
          jet_et_after_       = theDbe->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_       = theDbe->book1D("JETETA1_BEFORECUTS",chtitle, 20, -5, 5);
          jet_eta_after_       = theDbe->book1D("JETETA1_AFTERCUTS",chtitle, 20, -5, 5);
 
 //       snprintf(chtitle, 255, "Jet with 2nd highest E_{T} (%s)", jetTag_.label().data());
 //       jet2_et_before      = theDbe->book1D("JETET2_BEFORECUTS",chtitle, 20, 0., 200.0);
 //       jet2_et_after       = theDbe->book1D("JETET2_AFTERCUTS",chtitle, 20, 0., 200.0);
 
 
 
 }