TopDiLeptonOffline::MonitorEnsemble Class Reference

#include <TopDiLeptonOfflineDQM.h>

Public Types
enum	DecayChannel { NONE, DIMUON, DIELEC, ELECMU }
	different decay channels More...
enum	Level { STANDARD, VERBOSE, DEBUG }
	different verbosity levels More...
typedef reco::LeafCandidate::LorentzVector	LorentzVector

Public Member Functions
void	book (std::string directory)
	book histograms in subdirectory directory More...
void	fill (const edm::Event &event, const edm::EventSetup &setup)
	fill monitor histograms with electronId and jetCorrections More...
	MonitorEnsemble (const char *label, const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
	~MonitorEnsemble ()
	default destructor More...

Private Member Functions
bool	booked (const std::string histName) const
	check if histogram was booked More...
DecayChannel	decayChannel (const std::vector< const reco::PFCandidate * > &muons, const std::vector< const reco::PFCandidate * > &elecs) const
	determine dileptonic decay channel More...
void	fill (const edm::Event &event, const edm::TriggerResults &triggerTable, std::string channel, const std::vector< std::string > labels) const
	fill trigger monitoring histograms More...
void	fill (const std::string histName, double value) const
	fill histogram if it had been booked before More...
void	fill (const std::string histName, double xValue, double yValue) const
	fill histogram if it had been booked before (2-dim version) More...
void	fill (const std::string histName, double xValue, double yValue, double zValue) const
	fill histogram if it had been booked before (2-dim version) More...
void	loggerBinLabels (std::string hist)
	set labels for event logging histograms More...
std::string	monitorPath (const std::string &label) const
std::string	selectionPath (const std::string &label) const
void	triggerBinLabels (std::string channel, const std::vector< std::string > labels)
	set configurable labels for trigger monitoring histograms More...

Private Attributes
int	diElecLogged_
int	diMuonLogged_
std::vector< std::string >	diMuonPaths_
	trigger paths for di muon channel More...
double	eidCutValue_
StringCutObjectSelector < reco::PFCandidate > *	elecIso_
	extra isolation criterion on electron More...
int	elecMuLogged_
	number of logged interesting events More...
std::vector< std::string >	elecMuPaths_
edm::EDGetTokenT< edm::View < reco::PFCandidate > >	elecs_
StringCutObjectSelector < reco::PFCandidate > *	elecSelect_
	extra selection on electrons More...
edm::EDGetTokenT < edm::ValueMap< float > >	electronId_
	electronId label More...
std::map< std::string, MonitorElement * >	hists_
	histogram container More...
std::string	jetCorrector_
	jetCorrector More...
edm::EDGetTokenT < reco::JetIDValueMap >	jetIDLabel_
	jetID as an extra selection type More...
StringCutObjectSelector < reco::JetID > *	jetIDSelect_
	extra jetID selection on calo jets More...
edm::EDGetTokenT< edm::View < reco::Jet > >	jets_
	input sources for monitoring More...
std::string	jetSelect_
std::string	label_
	instance label More...
double	lowerEdge_
	mass window upper and lower edge More...
std::vector< edm::EDGetTokenT < edm::View< reco::MET > > >	mets_
	considers a vector of METs More...
StringCutObjectSelector < reco::PFCandidate, true > *	muonIso_
	extra isolation criterion on muon More...
edm::EDGetTokenT< edm::View < reco::PFCandidate > >	muons_
StringCutObjectSelector < reco::PFCandidate, true > *	muonSelect_
	extra selection on muons More...
DQMStore *	store_
	storage manager More...
edm::EDGetTokenT < edm::TriggerResults >	triggerTable_
	trigger table More...
double	upperEdge_
Level	verbosity_
	verbosity level for booking More...

Detailed Description

Definition at line 52 of file TopDiLeptonOfflineDQM.h.

Member Typedef Documentation

typedef reco::LeafCandidate::LorentzVector TopDiLeptonOffline::MonitorEnsemble::LorentzVector

make clear which LorentzVector to use for jet, electrons and muon buffering

Definition at line 62 of file TopDiLeptonOfflineDQM.h.

Member Enumeration Documentation

enum TopDiLeptonOffline::MonitorEnsemble::DecayChannel

different decay channels

Enumerator
NONE
DIMUON
DIELEC
ELECMU

Definition at line 64 of file TopDiLeptonOfflineDQM.h.

                    {
    NONE,
    DIMUON,
    DIELEC,
    ELECMU
  };

enum TopDiLeptonOffline::MonitorEnsemble::Level

different verbosity levels

Enumerator
STANDARD
VERBOSE
DEBUG

Definition at line 55 of file TopDiLeptonOfflineDQM.h.

             {
    STANDARD,
    VERBOSE,
    DEBUG
  };

Constructor & Destructor Documentation

MonitorEnsemble::MonitorEnsemble	(	const char *	label,
		const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	iC
	)

default contructor MonitorEnsemble(const char* label, const edm::ParameterSet& cfg);

Definition at line 14 of file TopDiLeptonOfflineDQM.cc.

References book(), DEBUG, diMuonPaths_, eidCutValue_, elecIso_, elecMuPaths_, elecs_, elecSelect_, electronId_, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), jetCorrector_, bTagSequences_cff::jetID, jetIDLabel_, jetIDSelect_, jets_, jetSelect_, lowerEdge_, mets_, muonIso_, muons_, muonSelect_, STANDARD, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, triggerTable_, upperEdge_, VERBOSE, and verbosity_.

    : label_(label),
      eidCutValue_(0.),
      elecIso_(0),
      elecSelect_(0),
      muonIso_(0),
      muonSelect_(0),
      jetIDSelect_(0),
      lowerEdge_(-1.),
      upperEdge_(-1.),
      elecMuLogged_(0),
      diMuonLogged_(0),
      diElecLogged_(0) {
  // sources have to be given; this PSet is not optional
  edm::ParameterSet sources = cfg.getParameter<edm::ParameterSet>("sources");
  muons_ = iC.consumes<edm::View<reco::PFCandidate> >(
      sources.getParameter<edm::InputTag>("muons"));
  elecs_ = iC.consumes<edm::View<reco::PFCandidate> >(
      sources.getParameter<edm::InputTag>("elecs"));
  jets_ = iC.consumes<edm::View<reco::Jet> >(
      sources.getParameter<edm::InputTag>("jets"));
  for (edm::InputTag const& tag :
       sources.getParameter<std::vector<edm::InputTag> >("mets"))
    mets_.push_back(iC.consumes<edm::View<reco::MET> >(tag));

  // elecExtras are optional; they may be omitted or empty
  if (cfg.existsAs<edm::ParameterSet>("elecExtras")) {
    edm::ParameterSet elecExtras =
        cfg.getParameter<edm::ParameterSet>("elecExtras");
    // select is optional; in case it's not found no
    // selection will be applied
    if (elecExtras.existsAs<std::string>("select")) {
      elecSelect_ = new StringCutObjectSelector<reco::PFCandidate>(
          elecExtras.getParameter<std::string>("select"));
    }
    // isolation is optional; in case it's not found no
    // isolation will be applied
    if (elecExtras.existsAs<std::string>("isolation")) {
      elecIso_ = new StringCutObjectSelector<reco::PFCandidate>(
          elecExtras.getParameter<std::string>("isolation"));
    }
    // electronId is optional; in case it's not found the
    // InputTag will remain empty
    if (elecExtras.existsAs<edm::ParameterSet>("electronId")) {
      edm::ParameterSet elecId =
          elecExtras.getParameter<edm::ParameterSet>("electronId");
      electronId_ = iC.consumes<edm::ValueMap<float> >(
          elecId.getParameter<edm::InputTag>("src"));
      eidCutValue_ = elecId.getParameter<double>("cutValue");
      //    eidPattern_= elecId.getParameter<int>("pattern");
    }
  }
  // muonExtras are optional; they may be omitted or empty
  if (cfg.existsAs<edm::ParameterSet>("muonExtras")) {
    edm::ParameterSet muonExtras =
        cfg.getParameter<edm::ParameterSet>("muonExtras");
    // select is optional; in case it's not found no
    // selection will be applied
    if (muonExtras.existsAs<std::string>("select")) {
      muonSelect_ = new StringCutObjectSelector<reco::PFCandidate, true>(
          muonExtras.getParameter<std::string>("select"));
    }
    // isolation is optional; in case it's not found no
    // isolation will be applied
    if (muonExtras.existsAs<std::string>("isolation")) {
      muonIso_ = new StringCutObjectSelector<reco::PFCandidate, true>(
          muonExtras.getParameter<std::string>("isolation"));
    }
  }
  // jetExtras are optional; they may be omitted or empty
  if (cfg.existsAs<edm::ParameterSet>("jetExtras")) {
    edm::ParameterSet jetExtras =
        cfg.getParameter<edm::ParameterSet>("jetExtras");
    // jetCorrector is optional; in case it's not found
    // the InputTag will remain empty
    if (jetExtras.existsAs<std::string>("jetCorrector")) {
      jetCorrector_ = jetExtras.getParameter<std::string>("jetCorrector");
    }
    // read jetID information if it exists
    if (jetExtras.existsAs<edm::ParameterSet>("jetID")) {
      edm::ParameterSet jetID =
          jetExtras.getParameter<edm::ParameterSet>("jetID");
      jetIDLabel_ = iC.consumes<reco::JetIDValueMap>(
          jetID.getParameter<edm::InputTag>("label"));
      jetIDSelect_ = new StringCutObjectSelector<reco::JetID>(
          jetID.getParameter<std::string>("select"));
    }
    // select is optional; in case it's not found no
    // selection will be applied (only implemented for
    // CaloJets at the moment)
    if (jetExtras.existsAs<std::string>("select")) {
      jetSelect_ = jetExtras.getParameter<std::string>("select");
    }
  }
  // triggerExtras are optional; they may be omitted or empty
  if (cfg.existsAs<edm::ParameterSet>("triggerExtras")) {
    edm::ParameterSet triggerExtras =
        cfg.getParameter<edm::ParameterSet>("triggerExtras");
    triggerTable_ = iC.consumes<edm::TriggerResults>(
        triggerExtras.getParameter<edm::InputTag>("src"));
    elecMuPaths_ =
        triggerExtras.getParameter<std::vector<std::string> >("pathsELECMU");
    diMuonPaths_ =
        triggerExtras.getParameter<std::vector<std::string> >("pathsDIMUON");
  }
  // massExtras is optional; in case it's not found no mass
  // window cuts are applied for the same flavor monitor
  // histograms
  if (cfg.existsAs<edm::ParameterSet>("massExtras")) {
    edm::ParameterSet massExtras =
        cfg.getParameter<edm::ParameterSet>("massExtras");
    lowerEdge_ = massExtras.getParameter<double>("lowerEdge");
    upperEdge_ = massExtras.getParameter<double>("upperEdge");
  }

  // setup the verbosity level for booking histograms;
  // per default the verbosity level will be set to
  // STANDARD. This will also be the chosen level in
  // the case when the monitoring PSet is not found
  verbosity_ = STANDARD;
  if (cfg.existsAs<edm::ParameterSet>("monitoring")) {
    edm::ParameterSet monitoring =
        cfg.getParameter<edm::ParameterSet>("monitoring");
    if (monitoring.getParameter<std::string>("verbosity") == "DEBUG")
      verbosity_ = DEBUG;
    if (monitoring.getParameter<std::string>("verbosity") == "VERBOSE")
      verbosity_ = VERBOSE;
    if (monitoring.getParameter<std::string>("verbosity") == "STANDARD")
      verbosity_ = STANDARD;
  }
  // and don't forget to do the histogram booking
  book(cfg.getParameter<std::string>("directory"));
}

TopDiLeptonOffline::MonitorEnsemble::~MonitorEnsemble ( )

inline

default destructor

Definition at line 77 of file TopDiLeptonOfflineDQM.h.

{};

Member Function Documentation

void MonitorEnsemble::book

(

std::string

directory

)

book histograms in subdirectory directory

Definition at line 150 of file TopDiLeptonOfflineDQM.cc.

References cond::rpcobimon::current, diMuonPaths_, elecMuPaths_, hists_, label_, loggerBinLabels(), cppFunctionSkipper::operator, STANDARD, store_, AlCaHLTBitMon_QueryRunRegistry::string, triggerBinLabels(), VERBOSE, and verbosity_.

Referenced by MonitorEnsemble().

                                              {
  // set up the current directory path
  std::string current(directory);
  current += label_;
  store_ = edm::Service<DQMStore>().operator->();
  store_->setCurrentFolder(current);

  // determine number of bins for trigger monitoring
  unsigned int nElecMu = elecMuPaths_.size();
  unsigned int nDiMuon = diMuonPaths_.size();

  // --- [STANDARD] --- //
  // Run Number
  hists_["RunNumb_"] =
      store_->book1D("RunNumber", "Run Nr.", 1.e4, 1.5e5, 3.e5);
  // invariant mass of opposite charge lepton pair (only filled for same flavor)
  hists_["invMass_"] = store_->book1D("InvMass", "M(lep1, lep2)", 80, 0., 320.);
  // invariant mass of opposite charge lepton pair (only filled for same flavor)
  hists_["invMassLog_"] =
      store_->book1D("InvMassLog", "log_{10}(M(lep1, lep2))", 80, .1, 2.5);
  // invariant mass of same charge lepton pair (log10 for low mass region, only
  // filled for same flavor)
  hists_["invMassWC_"] =
      store_->book1D("InvMassWC", "M_{WC}(L1, L2)", 80, 0., 320.);
  // invariant mass of same charge lepton pair (log10 for low mass region, only
  // filled for same flavor)
  hists_["invMassWCLog_"] =
      store_->book1D("InvMassLogWC", "log_{10}(M_{WC})", 80, .1, 2.5);
  // decay channel [1]: muon/muon, [2]:elec/elec, [3]:elec/muon
  hists_["decayChannel_"] =
      store_->book1D("DecayChannel", "Decay Channel", 3, 0, 3);
  // trigger efficiency estimates for the electron muon channel
  hists_["elecMuEff_"] =
      store_->book1D("ElecMuEff", "Eff(e/#mu paths)", nElecMu, 0., nElecMu);
  // monitored trigger occupancy for the electron muon channel
  hists_["elecMuMon_"] =
      store_->book1D("ElecMuMon", "Mon(e/#mu paths)", nElecMu, 0., nElecMu);
  // trigger efficiency estimates for the di muon channel
  hists_["diMuonEff_"] =
      store_->book1D("DiMuonEff", "Eff(#mu/#mu paths)", nDiMuon, 0., nDiMuon);
  // monitored trigger occupancy for the di muon channel
  hists_["diMuonMon_"] =
      store_->book1D("DiMuonMon", "Mon(#mu/#mu paths)", nDiMuon, 0., nDiMuon);
  // pt of the leading lepton
  hists_["lep1Pt_"] = store_->book1D("Lep1Pt", "pt(lep1)", 50, 0., 200.);
  // pt of the 2. leading lepton
  hists_["lep2Pt_"] = store_->book1D("Lep2Pt", "pt(lep2)", 50, 0., 200.);
  // multiplicity of jets with pt>30 (corrected to L2+L3)
  hists_["jetMult_"] = store_->book1D("JetMult", "N_{30}(jet)", 21, -0.5, 20.5);
  // MET (calo)
  hists_["metCalo_"] = store_->book1D("METCalo", "MET_{Calo}", 50, 0., 200.);

  // set bin labels for trigger monitoring
  triggerBinLabels(std::string("elecMu"), elecMuPaths_);
  triggerBinLabels(std::string("diMuon"), diMuonPaths_);
  // set bin labels for decayChannel_
  hists_["decayChannel_"]->setBinLabel(1, "#mu e", 1);
  hists_["decayChannel_"]->setBinLabel(2, "#mu #mu", 1);
  hists_["decayChannel_"]->setBinLabel(3, "e e", 1);

  if (verbosity_ == STANDARD) return;

  // --- [VERBOSE] --- //
  // mean eta of the candidate leptons
  hists_["sumEtaL1L2_"] =
      store_->book1D("SumEtaL1L2", "<#eta>(lep1, lep2)", 100, -5., 5.);
  // deltaEta between the 2 candidate leptons
  hists_["dEtaL1L2_"] =
      store_->book1D("DEtaL1L2", "#Delta#eta(lep1,lep2)", 80, -4., 4.);
  // deltaPhi between the 2 candidate leptons
  hists_["dPhiL1L2_"] =
      store_->book1D("DPhiL1L2", "#Delta#phi(lep1,lep2)", 64, -3.2, 3.2);
  // pt of the candidate electron (depending on the decay channel)
  hists_["elecPt_"] = store_->book1D("ElecPt", "pt(e)", 50, 0., 200.);
  // relative isolation of the candidate electron (depending on the decay
  // channel)
  hists_["elecRelIso_"] =
      store_->book1D("ElecRelIso", "Iso_{Rel}(e)", 50, 0., 1.);
  // pt of the canddiate muon (depending on the decay channel)
  hists_["muonPt_"] = store_->book1D("MuonPt", "pt(#mu)", 50, 0., 200.);
  // relative isolation of the candidate muon (depending on the decay channel)
  hists_["muonRelIso_"] = store_->book1D(
      "MuonRelIso", "Iso_{Rel}(#mu) (#Delta#beta Corrected)", 50, 0., 1.);
  // pt of the 1. leading jet (corrected to L2+L3)
  hists_["jet1Pt_"] = store_->book1D("Jet1Pt", "pt_{L2L3}(jet1)", 60, 0., 300.);
  // pt of the 2. leading jet (corrected to L2+L3)
  hists_["jet2Pt_"] = store_->book1D("Jet2Pt", "pt_{L2L3}(jet2)", 60, 0., 300.);
  // MET (PF)
  hists_["metPflow_"] = store_->book1D("METPflow", "MET_{Pflow}", 50, 0., 200.);
  // MET (TC)
  hists_["metTC_"] = store_->book1D("METTC", "MET_{TC}", 50, 0., 200.);
  // dz for muons (to suppress cosmis)
  hists_["muonDelZ_"] = store_->book1D("MuonDelZ", "d_{z}(#mu)", 50, -25., 25.);
  // dxy for muons (to suppress cosmics)
  hists_["muonDelXY_"] =
      store_->book2D("MuonDelXY", "d_{xy}(#mu)", 50, -1., 1., 50, -1., 1.);
  // lepton multiplicity after std isolation
  hists_["lepMultIso_"] = store_->book2D(
      "LepMultIso", "N_{Iso}(e) vs N_{Iso}(#mu)", 5, 0., 5., 5, 0., 5.);

  // set axes titles for dxy for muons
  hists_["muonDelXY_"]->setAxisTitle("x [cm]", 1);
  hists_["muonDelXY_"]->setAxisTitle("y [cm]", 2);
  // set axes titles for lepton multiplicity after std isolation
  hists_["lepMultIso_"]->setAxisTitle("N_{Iso}(#mu)", 1);
  hists_["lepMultIso_"]->setAxisTitle("N_{Iso}(elec)", 2);

  if (verbosity_ == VERBOSE) return;

  // --- [DEBUG] --- //
  // electron multiplicity after std isolation
  hists_["elecMultIso_"] =
      store_->book1D("ElecMultIso", "N_{Iso}(e)", 11, -0.5, 10.5);
  // muon multiplicity after std isolation
  hists_["muonMultIso_"] =
      store_->book1D("MuonMultIso", "N_{Iso}(#mu)", 11, -0.5, 10.5);
  // charged hadron isolation component of the candidate muon (depending on the
  // decay channel)
  hists_["muonChHadIso_"] = store_->book1D(
      "MuonChHadIsoComp", "ChHad_{IsoComponent}(#mu)", 50, 0., 5.);
  // neutral hadron isolation component of the candidate muon (depending on the
  // decay channel)
  hists_["muonNeHadIso_"] = store_->book1D(
      "MuonNeHadIsoComp", "NeHad_{IsoComponent}(#mu)", 50, 0., 5.);
  // photon isolation component of the candidate muon (depending on the decay
  // channel)
  hists_["muonPhIso_"] =
      store_->book1D("MuonPhIsoComp", "Photon_{IsoComponent}(#mu)", 50, 0., 5.);
  // charged hadron isolation component of the candidate electron (depending on
  // the decay channel)
  hists_["elecChHadIso_"] = store_->book1D(
      "ElectronChHadIsoComp", "ChHad_{IsoComponent}(e)", 50, 0., 5.);
  // neutral hadron isolation component of the candidate electron (depending on
  // the decay channel)
  hists_["elecNeHadIso_"] = store_->book1D(
      "ElectronNeHadIsoComp", "NeHad_{IsoComponent}(e)", 50, 0., 5.);
  // photon isolation component of the candidate electron (depending on the
  // decay channel)
  hists_["elecPhIso_"] = store_->book1D("ElectronPhIsoComp",
                                        "Photon_{IsoComponent}(e)", 50, 0., 5.);
  // eta of the leading jet
  hists_["jet1Eta_"] = store_->book1D("Jet1Eta", "#eta(jet1)", 30, -5., 5.);
  // eta of the 2. leading jet
  hists_["jet2Eta_"] = store_->book1D("Jet2Eta", "#eta(jet2)", 30, -5., 5.);
  // pt of the 1. leading jet (not corrected)
  hists_["jet1PtRaw_"] =
      store_->book1D("Jet1PtRaw", "pt_{Raw}(jet1)", 60, 0., 300.);
  // pt of the 2. leading jet (not corrected)
  hists_["jet2PtRaw_"] =
      store_->book1D("Jet2PtRaw", "pt_{Raw}(jet2)", 60, 0., 300.);
  // deltaEta between the 2 leading jets
  hists_["dEtaJet1Jet2_"] =
      store_->book1D("DEtaJet1Jet2", "#Delta#eta(jet1,jet2)", 80, -4., 4.);
  // deltaEta between the lepton and the leading jet
  hists_["dEtaJet1Lep1_"] =
      store_->book1D("DEtaJet1Lep1", "#Delta#eta(jet1,lep1)", 80, -4., 4.);
  // deltaEta between the lepton and MET
  hists_["dEtaLep1MET_"] =
      store_->book1D("DEtaLep1MET", "#Delta#eta(lep1,MET)", 80, -4., 4.);
  // deltaEta between leading jet and MET
  hists_["dEtaJet1MET_"] =
      store_->book1D("DEtaJet1MET", "#Delta#eta(jet1,MET)", 80, -4., 4.);
  // deltaPhi of 2 leading jets
  hists_["dPhiJet1Jet2_"] =
      store_->book1D("DPhiJet1Jet2", "#Delta#phi(jet1,jet2)", 64, -3.2, 3.2);
  // deltaPhi of 1. lepton and 1. jet
  hists_["dPhiJet1Lep1_"] =
      store_->book1D("DPhiJet1Lep1", "#Delta#phi(jet1,lep1)", 64, -3.2, 3.2);
  // deltaPhi of 1. lepton and MET
  hists_["dPhiLep1MET_"] =
      store_->book1D("DPhiLep1MET", "#Delta#phi(lep1,MET)", 64, -3.2, 3.2);
  // deltaPhi of 1. jet and MET
  hists_["dPhiJet1MET_"] =
      store_->book1D("DPhiJet1MET", "#Delta#phi(jet1,MET)", 64, -3.2, 3.2);
  // selected dimuon events
  hists_["diMuonLogger_"] = store_->book2D(
      "DiMuonLogger", "Logged DiMuon Events", 8, 0., 8., 10, 0., 10.);
  // selected dielec events
  hists_["diElecLogger_"] = store_->book2D(
      "DiElecLogger", "Logged DiElec Events", 8, 0., 8., 10, 0., 10.);
  // selected elemu events
  hists_["elecMuLogger_"] = store_->book2D(
      "ElecMuLogger", "Logged ElecMu Events", 8, 0., 8., 10, 0., 10.);

  // set bin labels for trigger monitoring
  loggerBinLabels(std::string("diMuonLogger_"));
  loggerBinLabels(std::string("diElecLogger_"));
  loggerBinLabels(std::string("elecMuLogger_"));
  return;
}

bool TopDiLeptonOffline::MonitorEnsemble::booked ( const std::string  histName ) const

inlineprivate

check if histogram was booked

Definition at line 110 of file TopDiLeptonOfflineDQM.h.

References hists_.

Referenced by fill().

                                              {
    return hists_.find(histName.c_str()) != hists_.end();
  };

MonitorEnsemble::DecayChannel MonitorEnsemble::decayChannel ( const std::vector< const reco::PFCandidate * > &  muons, const std::vector< const reco::PFCandidate * > &  elecs  ) const

inlineprivate

determine dileptonic decay channel

Definition at line 255 of file TopDiLeptonOfflineDQM.h.

References DIELEC, DIMUON, ELECMU, and NONE.

Referenced by fill().

                                                          {
  DecayChannel type = NONE;
  if (muons.size() > 1) {
    type = DIMUON;
  } else if (elecs.size() > 1) {
    type = DIELEC;
  } else if (!elecs.empty() && !muons.empty()) {
    type = ELECMU;
  }
  return type;
}

void MonitorEnsemble::fill	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

fill monitor histograms with electronId and jetCorrections

Definition at line 341 of file TopDiLeptonOfflineDQM.cc.

References decayChannel(), reco::deltaPhi(), reco::deltaR(), DIELEC, diElecLogged_, DIMUON, diMuonLogged_, diMuonPaths_, eidCutValue_, elecIso_, ELECMU, elecMuLogged_, elecMuPaths_, elecs_, elecSelect_, electronId_, reco::LeafCandidate::et(), eta(), edm::EventAuxiliary::event(), edm::Event::eventAuxiliary(), edm::EventSetup::find(), edm::Event::getByToken(), JetCorrector::getJetCorrector(), hists_, customizeTrackingMonitorSeedNumber::idx, edm::Ref< C, T, F >::isNull(), edm::EDGetTokenT< T >::isUninitialized(), metsig::jet, jetCorrector_, bTagSequences_cff::jetID, jetIDLabel_, jetIDSelect_, fwrapper::jets, jets_, jetSelect_, lowerEdge_, edm::EventAuxiliary::luminosityBlock(), edm::eventsetup::EventSetupRecordKey::makeKey(), max(), caloMETBenchmarkGeneric_cfi::met, mets_, VarParsing::mult, trackingTruthProducerFastSim_cfi::muon, muonIso_, patZpeak::muons, muons_, muonSelect_, muon::overlap(), phi, RecoTauCleanerPlugins::pt, reco::LeafCandidate::pt(), edm::EventAuxiliary::run(), reco::Jet::scaleEnergy(), EgammaValidation_Wenu_cff::sel, triggerTable_, and upperEdge_.

Referenced by fill().

                                                       {
  // fetch trigger event if configured such
  edm::Handle<edm::TriggerResults> triggerTable;
  if (!triggerTable_.isUninitialized()) {
    if (!event.getByToken(triggerTable_, triggerTable)) return;
  }
  /*
  ------------------------------------------------------------

  Run and Inst. Luminosity information (Inst. Lumi. filled now with a dummy
  value=5.0)

  ------------------------------------------------------------
  */

  if (!event.eventAuxiliary().run()) return;
  fill("RunNumb_", event.eventAuxiliary().run());

  double dummy = 5.;
  fill("InstLumi_", dummy);

  /*
  ------------------------------------------------------------

  Muon Selection

  ------------------------------------------------------------
  */

  std::vector<const reco::PFCandidate*> isoMuons;

  edm::Handle<edm::View<reco::PFCandidate> > muons;
  edm::View<reco::PFCandidate>::const_iterator muonit;

  if (!event.getByToken(muons_, muons)) return;

  for (edm::View<reco::PFCandidate>::const_iterator muonit = muons->begin();
       muonit != muons->end(); ++muonit) {

    if (muonit->muonRef().isNull()) continue;
    reco::MuonRef muon = muonit->muonRef();

    if (muon->innerTrack().isNull()) continue;

    if (muon->isGlobalMuon()) {
      fill("muonDelZ_", muon->innerTrack()->vz());  // CB using inner track!
      fill("muonDelXY_", muon->innerTrack()->vx(), muon->innerTrack()->vy());

      // apply selection
      if (!muonSelect_ || (*muonSelect_)(*muonit)) {

        double chHadPt = muon->pfIsolationR04().sumChargedHadronPt;
        double neHadEt = muon->pfIsolationR04().sumNeutralHadronEt;
        double phoEt = muon->pfIsolationR04().sumPhotonEt;

        double pfRelIso =
            (chHadPt +
             max(0., neHadEt + phoEt - 0.5 * muon->pfIsolationR04().sumPUPt)) /
            muon->pt();  // CB dBeta corrected iso!

        fill("muonRelIso_", pfRelIso);

        fill("muonChHadIso_", chHadPt);
        fill("muonNeHadIso_", neHadEt);
        fill("muonPhIso_", phoEt);

        if (!muonIso_ || (*muonIso_)(*muonit)) isoMuons.push_back(&(*muonit));
      }
    }
  }

  fill("muonMultIso_", isoMuons.size());

  /*
  ------------------------------------------------------------

  Electron Selection

  ------------------------------------------------------------
  */

  // buffer isolated electronss
  std::vector<const reco::PFCandidate*> isoElecs;
  edm::Handle<edm::ValueMap<float> > electronId;
  if (!electronId_.isUninitialized()) {
    if (!event.getByToken(electronId_, electronId)) return;
  }
  edm::Handle<edm::View<reco::PFCandidate> > elecs;
  if (!event.getByToken(elecs_, elecs)) return;

  for (edm::View<reco::PFCandidate>::const_iterator elec = elecs->begin();
       elec != elecs->end(); ++elec) {
    if (elec->gsfElectronRef().isNull()) {
      continue;
    }
    reco::GsfElectronRef gsf_el = elec->gsfElectronRef();
    // restrict to electrons with good electronId
    if (electronId_.isUninitialized() ? true : ((double)(*electronId)[gsf_el] >=
                                                eidCutValue_)) {
      // apply preselection
      if (!elecSelect_ || (*elecSelect_)(*elec)) {
        double el_ChHadIso = gsf_el->pfIsolationVariables().sumChargedHadronPt;
        double el_NeHadIso = gsf_el->pfIsolationVariables().sumNeutralHadronEt;
        double el_PhIso = gsf_el->pfIsolationVariables().sumPhotonEt;
        double el_pfRelIso =
            (el_ChHadIso +
             max(0., el_NeHadIso + el_PhIso -
                         0.5 * gsf_el->pfIsolationVariables().sumPUPt)) /
            gsf_el->pt();
        fill("elecRelIso_", el_pfRelIso);
        fill("elecChHadIso_", el_ChHadIso);
        fill("elecNeHadIso_", el_NeHadIso);
        fill("elecPhIso_", el_PhIso);
        if (!elecIso_ || (*elecIso_)(*elec)) isoElecs.push_back(&(*elec));
      }
    }
  }
  fill("elecMultIso_", isoElecs.size());

  /*
  ------------------------------------------------------------

  Jet Selection

  ------------------------------------------------------------
  */

  const JetCorrector* corrector = 0;
  if (!jetCorrector_.empty()) {
    // check whether a jet correcto is in the event setup or not
    if (setup.find(edm::eventsetup::EventSetupRecordKey::makeKey<
            JetCorrectionsRecord>())) {
      corrector = JetCorrector::getJetCorrector(jetCorrector_, setup);
    } else {
      edm::LogVerbatim("TopDiLeptonOfflineDQM")
          << "\n"
          << "-----------------------------------------------------------------"
             "-------------------- \n"
          << " No JetCorrectionsRecord available from EventSetup:              "
             "                     \n"
          << "  - Jets will not be corrected.                                  "
             "                     \n"
          << "  - If you want to change this add the following lines to your "
             "cfg file:              \n"
          << "                                                                 "
             "                     \n"
          << "  ## load jet corrections                                        "
             "                     \n"
          << "  "
             "process.load(\"JetMETCorrections.Configuration."
             "JetCorrectionServicesAllAlgos_cff\") \n"
          << "  process.prefer(\"ak5CaloL2L3\")                                "
             "                     \n"
          << "                                                                 "
             "                     \n"
          << "-----------------------------------------------------------------"
             "-------------------- \n";
    }
  }

  unsigned int mult = 0;
  // buffer leadingJets
  std::vector<reco::Jet> leadingJets;
  edm::Handle<edm::View<reco::Jet> > jets;
  if (!event.getByToken(jets_, jets)) return;

  edm::Handle<reco::JetIDValueMap> jetID;
  if (jetIDSelect_) {
    if (!event.getByToken(jetIDLabel_, jetID)) return;
  }

  for (edm::View<reco::Jet>::const_iterator jet = jets->begin();
       jet != jets->end(); ++jet) {
    unsigned int idx = jet - jets->begin();
    if (jetIDSelect_ &&
        dynamic_cast<const reco::CaloJet*>(jets->refAt(idx).get())) {
      if (!(*jetIDSelect_)((*jetID)[jets->refAt(idx)])) continue;
    }
    // chekc additional jet selection for calo, pf and bare reco jets
    if (dynamic_cast<const reco::CaloJet*>(&*jet)) {
      reco::CaloJet sel = dynamic_cast<const reco::CaloJet&>(*jet);
      sel.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
      StringCutObjectSelector<reco::CaloJet> jetSelect(jetSelect_);
      if (!jetSelect(sel)) {
        continue;
      }
    } else if (dynamic_cast<const reco::PFJet*>(&*jet)) {
      reco::PFJet sel = dynamic_cast<const reco::PFJet&>(*jet);
      sel.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
      StringCutObjectSelector<reco::PFJet> jetSelect(jetSelect_);
      if (!jetSelect(sel)) continue;
    } else {
      reco::Jet sel = *jet;
      sel.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
      StringCutObjectSelector<reco::Jet> jetSelect(jetSelect_);
      if (!jetSelect(sel)) continue;
    }
    // check for overlaps
    bool overlap = false;
    for (std::vector<const reco::PFCandidate*>::const_iterator elec =
             isoElecs.begin();
         elec != isoElecs.end(); ++elec) {
      if (reco::deltaR((*elec)->eta(), (*elec)->phi(), jet->eta(), jet->phi()) <
          0.4) {
        overlap = true;
        break;
      }
    }
    if (overlap) {
      continue;
    }
    // prepare jet to fill monitor histograms
    reco::Jet monitorJet = *jet;
    monitorJet.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
    ++mult;  // determine jet multiplicity
    if (idx == 0) {
      leadingJets.push_back(monitorJet);
      fill("jet1Pt_", monitorJet.pt());
      fill("jet1PtRaw_", jet->pt());
      fill("jet1Eta_", jet->eta());
    }
    if (idx == 1) {
      leadingJets.push_back(monitorJet);
      fill("jet2Pt_", monitorJet.pt());
      fill("jet2PtRaw_", jet->pt());
      fill("jet2Eta_", jet->eta());
    }
  }
  if (leadingJets.size() > 1) {
    fill("dEtaJet1Jet2_", leadingJets[0].eta() - leadingJets[1].eta());
    fill("dPhiJet1Jet2_",
         reco::deltaPhi(leadingJets[0].phi(), leadingJets[1].phi()));
    if (!isoMuons.empty()) {
      if (isoElecs.empty() || isoMuons[0]->pt() > isoElecs[0]->pt()) {
        fill("dEtaJet1Lep1_", isoMuons[0]->eta() - leadingJets[0].eta());
        fill("dPhiJet1Lep1_",
             reco::deltaPhi(isoMuons[0]->phi(), leadingJets[0].phi()));
      }
    }
    if (!isoElecs.empty()) {
      if (isoMuons.empty() || isoElecs[0]->pt() > isoMuons[0]->pt()) {
        fill("dEtaJet1Lep1_", isoElecs[0]->eta() - leadingJets[0].eta());
        fill("dPhiJet1Lep1_",
             reco::deltaPhi(isoElecs[0]->phi(), leadingJets[0].phi()));
      }
    }
  }
  fill("jetMult_", mult);

  /*
  ------------------------------------------------------------

  MET Selection

  ------------------------------------------------------------
  */

  // buffer for event logging
  reco::MET caloMET;
  for (std::vector<edm::EDGetTokenT<edm::View<reco::MET> > >::const_iterator
           met_ = mets_.begin();
       met_ != mets_.end(); ++met_) {

    edm::Handle<edm::View<reco::MET> > met;
    if (!event.getByToken(*met_, met)) continue;

    if (met->begin() != met->end()) {
      unsigned int idx = met_ - mets_.begin();
      if (idx == 0) {
        caloMET = *met->begin();
        fill("metCalo_", met->begin()->et());
        if (!leadingJets.empty()) {
          fill("dEtaJet1MET_", leadingJets[0].eta() - met->begin()->eta());
          fill("dPhiJet1MET_",
               reco::deltaPhi(leadingJets[0].phi(), met->begin()->phi()));
        }
        if (!isoMuons.empty()) {
          if (isoElecs.empty() || isoMuons[0]->pt() > isoElecs[0]->pt()) {
            fill("dEtaLep1MET_", isoMuons[0]->eta() - met->begin()->eta());
            fill("dPhiLep1MET_",
                 reco::deltaPhi(isoMuons[0]->phi(), met->begin()->phi()));
          }
        }
        if (!isoElecs.empty()) {
          if (isoMuons.empty() || isoElecs[0]->pt() > isoMuons[0]->pt()) {
            fill("dEtaLep1MET_", isoElecs[0]->eta() - met->begin()->eta());
            fill("dPhiLep1MET_",
                 reco::deltaPhi(isoElecs[0]->phi(), met->begin()->phi()));
          }
        }
      }
      if (idx == 1) {
        fill("metTC_", met->begin()->et());
      }
      if (idx == 2) {
        fill("metPflow_", met->begin()->et());
      }
    }
  }

  /*
  ------------------------------------------------------------

  Event Monitoring

  ------------------------------------------------------------
  */

  // check number of isolated leptons
  fill("lepMultIso_", isoMuons.size(), isoElecs.size());
  // ELECMU channel
  if (decayChannel(isoMuons, isoElecs) == ELECMU) {
    fill("decayChannel_", 0.5);
    double mass = (isoElecs[0]->p4() + isoMuons[0]->p4()).mass();
    if ((lowerEdge_ == -1. && upperEdge_ == -1.) ||
        (lowerEdge_ < mass && mass < upperEdge_)) {

      fill("dEtaL1L2_", isoElecs[0]->eta() - isoMuons[0]->eta());
      fill("sumEtaL1L2_", (isoElecs[0]->eta() + isoMuons[0]->eta()) / 2);
      fill("dPhiL1L2_", reco::deltaPhi(isoElecs[0]->phi(), isoMuons[0]->eta()));
      fill("elecPt_", isoElecs[0]->pt());
      fill("muonPt_", isoMuons[0]->pt());
      fill("lep1Pt_", isoElecs[0]->pt() > isoMuons[0]->pt()
                          ? isoElecs[0]->pt()
                          : isoMuons[0]->pt());
      fill("lep2Pt_", isoElecs[0]->pt() > isoMuons[0]->pt()
                          ? isoMuons[0]->pt()
                          : isoElecs[0]->pt());
      // fill plots for trigger monitoring
      if (!triggerTable_.isUninitialized())
        fill(event, *triggerTable, "elecMu", elecMuPaths_);
      if (elecMuLogged_ <= hists_.find("elecMuLogger_")->second->getNbinsY()) {
        // log runnumber, lumi block, event number & some
        // more pysics infomation for interesting events
        fill("elecMuLogger_", 0.5, elecMuLogged_ + 0.5,
             event.eventAuxiliary().run());
        fill("elecMuLogger_", 1.5, elecMuLogged_ + 0.5,
             event.eventAuxiliary().luminosityBlock());
        fill("elecMuLogger_", 2.5, elecMuLogged_ + 0.5,
             event.eventAuxiliary().event());
        fill("elecMuLogger_", 3.5, elecMuLogged_ + 0.5, isoMuons[0]->pt());
        fill("elecMuLogger_", 4.5, elecMuLogged_ + 0.5, isoElecs[0]->pt());
        if (leadingJets.size() > 0)
          fill("elecMuLogger_", 5.5, elecMuLogged_ + 0.5, leadingJets[0].pt());
        if (leadingJets.size() > 1)
          fill("elecMuLogger_", 6.5, elecMuLogged_ + 0.5, leadingJets[1].pt());
        fill("elecMuLogger_", 7.5, elecMuLogged_ + 0.5, caloMET.et());
        ++elecMuLogged_;
      }
    }
  }

  // DIMUON channel
  if (decayChannel(isoMuons, isoElecs) == DIMUON) {
    fill("decayChannel_", 1.5);
    int charge = isoMuons[0]->charge() * isoMuons[1]->charge();
    double mass = (isoMuons[0]->p4() + isoMuons[1]->p4()).mass();

    fill(charge < 0 ? "invMass_" : "invMassWC_", mass);
    fill(charge < 0 ? "invMassLog_" : "invMassWCLog_", log10(mass));
    if ((lowerEdge_ == -1. && upperEdge_ == -1.) ||
        (lowerEdge_ < mass && mass < upperEdge_)) {
      fill("dEtaL1L2_", isoMuons[0]->eta() - isoMuons[1]->eta());
      fill("sumEtaL1L2_", (isoMuons[0]->eta() + isoMuons[1]->eta()) / 2);
      fill("dPhiL1L2_", reco::deltaPhi(isoMuons[0]->phi(), isoMuons[1]->phi()));
      fill("muonPt_", isoMuons[0]->pt());
      fill("muonPt_", isoMuons[1]->pt());
      fill("lep1Pt_", isoMuons[0]->pt());
      fill("lep2Pt_", isoMuons[1]->pt());
      // fill plots for trigger monitoring
      if (!triggerTable_.isUninitialized())
        fill(event, *triggerTable, "diMuon", diMuonPaths_);
      if (diMuonLogged_ <= hists_.find("diMuonLogger_")->second->getNbinsY()) {
        // log runnumber, lumi block, event number & some
        // more pysics infomation for interesting events
        fill("diMuonLogger_", 0.5, diMuonLogged_ + 0.5,
             event.eventAuxiliary().run());
        fill("diMuonLogger_", 1.5, diMuonLogged_ + 0.5,
             event.eventAuxiliary().luminosityBlock());
        fill("diMuonLogger_", 2.5, diMuonLogged_ + 0.5,
             event.eventAuxiliary().event());
        fill("diMuonLogger_", 3.5, diMuonLogged_ + 0.5, isoMuons[0]->pt());
        fill("diMuonLogger_", 4.5, diMuonLogged_ + 0.5, isoMuons[1]->pt());
        if (leadingJets.size() > 0)
          fill("diMuonLogger_", 5.5, diMuonLogged_ + 0.5, leadingJets[0].pt());
        if (leadingJets.size() > 1)
          fill("diMuonLogger_", 6.5, diMuonLogged_ + 0.5, leadingJets[1].pt());
        fill("diMuonLogger_", 7.5, diMuonLogged_ + 0.5, caloMET.et());
        ++diMuonLogged_;
      }
    }
  }

  // DIELEC channel
  if (decayChannel(isoMuons, isoElecs) == DIELEC) {
    fill("decayChannel_", 2.5);
    int charge = isoElecs[0]->charge() * isoElecs[1]->charge();
    double mass = (isoElecs[0]->p4() + isoElecs[1]->p4()).mass();
    fill(charge < 0 ? "invMass_" : "invMassWC_", mass);
    fill(charge < 0 ? "invMassLog_" : "invMassWCLog_", log10(mass));
    if ((lowerEdge_ == -1. && upperEdge_ == -1.) ||
        (lowerEdge_ < mass && mass < upperEdge_)) {
      fill("dEtaL1L2_", isoElecs[0]->eta() - isoElecs[1]->eta());
      fill("sumEtaL1L2_", (isoElecs[0]->eta() + isoElecs[1]->eta()) / 2);
      fill("dPhiL1L2_", reco::deltaPhi(isoElecs[0]->phi(), isoElecs[1]->phi()));
      fill("elecPt_", isoElecs[0]->pt());
      fill("elecPt_", isoElecs[1]->pt());
      fill("lep1Pt_", isoElecs[0]->pt());
      fill("lep2Pt_", isoElecs[1]->pt());
      if (diElecLogged_ <= hists_.find("diElecLogger_")->second->getNbinsY()) {
        // log runnumber, lumi block, event number & some
        // more pysics infomation for interesting events
        fill("diElecLogger_", 0.5, diElecLogged_ + 0.5,
             event.eventAuxiliary().run());
        fill("diElecLogger_", 1.5, diElecLogged_ + 0.5,
             event.eventAuxiliary().luminosityBlock());
        fill("diElecLogger_", 2.5, diElecLogged_ + 0.5,
             event.eventAuxiliary().event());
        fill("diElecLogger_", 3.5, diElecLogged_ + 0.5, isoElecs[0]->pt());
        fill("diElecLogger_", 4.5, diElecLogged_ + 0.5, isoElecs[1]->pt());
        if (leadingJets.size() > 0)
          fill("diElecLogger_", 5.5, diElecLogged_ + 0.5, leadingJets[0].pt());
        if (leadingJets.size() > 1)
          fill("diElecLogger_", 6.5, diElecLogged_ + 0.5, leadingJets[1].pt());
        fill("diElecLogger_", 7.5, diElecLogged_ + 0.5, caloMET.et());
        ++diElecLogged_;
      }
    }
  }
}

void MonitorEnsemble::fill ( const edm::Event &  event, const edm::TriggerResults &  triggerTable, std::string  channel, const std::vector< std::string >  labels  ) const

inlineprivate

fill trigger monitoring histograms

Definition at line 235 of file TopDiLeptonOfflineDQM.h.

References accept(), fill(), hists_, customizeTrackingMonitorSeedNumber::idx, monitorPath(), selectionPath(), and relativeConstraints::value.

                                                                           {
  for (unsigned int idx = 0; idx < labels.size(); ++idx) {
    if (accept(event, triggerTable, monitorPath(labels[idx]))) {
      fill((channel + "Mon_").c_str(), idx + 0.5);
      // take care to fill triggerMon_ before evts is being called
      int evts = hists_.find((channel + "Mon_").c_str())
                     ->second->getBinContent(idx + 1);
      double value = hists_.find((channel + "Eff_").c_str())
                         ->second->getBinContent(idx + 1);
      fill(
          (channel + "Eff_").c_str(), idx + 0.5,
          1. / evts * (accept(event, triggerTable, selectionPath(labels[idx])) -
                       value));
    }
  }
}

void TopDiLeptonOffline::MonitorEnsemble::fill ( const std::string  histName, double  value  ) const

inlineprivate

fill histogram if it had been booked before

Definition at line 114 of file TopDiLeptonOfflineDQM.h.

References booked(), hists_, and edm::second().

                                                          {
    if (booked(histName.c_str()))
      hists_.find(histName.c_str())->second->Fill(value);
  };

void TopDiLeptonOffline::MonitorEnsemble::fill ( const std::string  histName, double  xValue, double  yValue  ) const

inlineprivate

fill histogram if it had been booked before (2-dim version)

Definition at line 119 of file TopDiLeptonOfflineDQM.h.

References booked(), hists_, and edm::second().

                                                                          {
    if (booked(histName.c_str()))
      hists_.find(histName.c_str())->second->Fill(xValue, yValue);
  };

void TopDiLeptonOffline::MonitorEnsemble::fill ( const std::string  histName, double  xValue, double  yValue, double  zValue  ) const

inlineprivate

fill histogram if it had been booked before (2-dim version)

Definition at line 124 of file TopDiLeptonOfflineDQM.h.

References booked(), hists_, and edm::second().

                                 {
    if (booked(histName.c_str()))
      hists_.find(histName.c_str())->second->Fill(xValue, yValue, zValue);
  };

void MonitorEnsemble::loggerBinLabels ( std::string  hist )

inlineprivate

set labels for event logging histograms

Definition at line 198 of file TopDiLeptonOfflineDQM.h.

References hists_.

Referenced by book().

                                                           {
  // set axes titles for selected events
  hists_[hist.c_str()]->getTH1()->SetOption("TEXT");
  hists_[hist.c_str()]->setBinLabel(1, "Run", 1);
  hists_[hist.c_str()]->setBinLabel(2, "Block", 1);
  hists_[hist.c_str()]->setBinLabel(3, "Event", 1);
  hists_[hist.c_str()]->setBinLabel(6, "pt_{L2L3}(jet1)", 1);
  hists_[hist.c_str()]->setBinLabel(7, "pt_{L2L3}(jet2)", 1);
  hists_[hist.c_str()]->setBinLabel(8, "MET_{Calo}", 1);
  hists_[hist.c_str()]->setAxisTitle("logged evts", 2);

  if (hist == "diMuonLogger_") {
    hists_[hist.c_str()]->setBinLabel(4, "pt(muon)", 1);
    hists_[hist.c_str()]->setBinLabel(5, "pt(muon)", 1);
  }
  if (hist == "diElecLogger_") {
    hists_[hist.c_str()]->setBinLabel(4, "pt(elec)", 1);
    hists_[hist.c_str()]->setBinLabel(5, "pt(elec)", 1);
  }
  if (hist == "elecMuLogger_") {
    hists_[hist.c_str()]->setBinLabel(4, "pt(elec)", 1);
    hists_[hist.c_str()]->setBinLabel(5, "pt(muon)", 1);
  }
}

std::string TopDiLeptonOffline::MonitorEnsemble::monitorPath ( const std::string &  label ) const

inlineprivate

deduce monitorPath from label, the label is expected to be of type 'selectionPath:monitorPath'

Definition at line 87 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and triggerBinLabels().

                                                      {
    return label.substr(label.find(':') + 1);
  };

std::string TopDiLeptonOffline::MonitorEnsemble::selectionPath ( const std::string &  label ) const

inlineprivate

deduce selectionPath from label, the label is expected to be of type 'selectionPath:monitorPath'

Definition at line 92 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and triggerBinLabels().

                                                        {
    return label.substr(0, label.find(':'));
  };

void MonitorEnsemble::triggerBinLabels ( std::string  channel, const std::vector< std::string >  labels  )

inlineprivate

set configurable labels for trigger monitoring histograms

Definition at line 223 of file TopDiLeptonOfflineDQM.h.

References hists_, customizeTrackingMonitorSeedNumber::idx, monitorPath(), and selectionPath().

Referenced by book().

                                                          {
  for (unsigned int idx = 0; idx < labels.size(); ++idx) {
    hists_[(channel + "Mon_").c_str()]
        ->setBinLabel(idx + 1, "[" + monitorPath(labels[idx]) + "]", 1);
    hists_[(channel + "Eff_").c_str()]
        ->setBinLabel(idx + 1, "[" + selectionPath(labels[idx]) + "]|[" +
                                   monitorPath(labels[idx]) + "]",
                      1);
  }
}

Member Data Documentation

int TopDiLeptonOffline::MonitorEnsemble::diElecLogged_

private

Definition at line 191 of file TopDiLeptonOfflineDQM.h.

Referenced by fill().

int TopDiLeptonOffline::MonitorEnsemble::diMuonLogged_

private

Definition at line 191 of file TopDiLeptonOfflineDQM.h.

Referenced by fill().

std::vector<std::string> TopDiLeptonOffline::MonitorEnsemble::diMuonPaths_

private

trigger paths for di muon channel

Definition at line 149 of file TopDiLeptonOfflineDQM.h.

Referenced by book(), fill(), and MonitorEnsemble().

double TopDiLeptonOffline::MonitorEnsemble::eidCutValue_

private

electronId pattern we expect the following pattern: 0: fails 1: passes electron ID only 2: passes electron Isolation only 3: passes electron ID and Isolation only 4: passes conversion rejection 5: passes conversion rejection and ID 6: passes conversion rejection and Isolation 7: passes the whole selection As described on https://twiki.cern.ch/twiki/bin/view/CMS/SimpleCutBasedEleID

Definition at line 166 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

StringCutObjectSelector<reco::PFCandidate>* TopDiLeptonOffline::MonitorEnsemble::elecIso_

private

extra isolation criterion on electron

Definition at line 168 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

int TopDiLeptonOffline::MonitorEnsemble::elecMuLogged_

private

number of logged interesting events

Definition at line 191 of file TopDiLeptonOfflineDQM.h.

Referenced by fill().

std::vector<std::string> TopDiLeptonOffline::MonitorEnsemble::elecMuPaths_

private

trigger paths for monitoring, expected to be of form signalPath:MonitorPath

Definition at line 147 of file TopDiLeptonOfflineDQM.h.

Referenced by book(), fill(), and MonitorEnsemble().

edm::EDGetTokenT<edm::View<reco::PFCandidate> > TopDiLeptonOffline::MonitorEnsemble::elecs_

private

Definition at line 138 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

StringCutObjectSelector<reco::PFCandidate>* TopDiLeptonOffline::MonitorEnsemble::elecSelect_

private

extra selection on electrons

Definition at line 170 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDGetTokenT<edm::ValueMap<float> > TopDiLeptonOffline::MonitorEnsemble::electronId_

private

electronId label

Definition at line 152 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::map<std::string, MonitorElement*> TopDiLeptonOffline::MonitorEnsemble::hists_

private

histogram container

Definition at line 195 of file TopDiLeptonOfflineDQM.h.

Referenced by book(), booked(), fill(), loggerBinLabels(), and triggerBinLabels().

std::string TopDiLeptonOffline::MonitorEnsemble::jetCorrector_

private

jetCorrector

Definition at line 179 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDGetTokenT<reco::JetIDValueMap> TopDiLeptonOffline::MonitorEnsemble::jetIDLabel_

private

jetID as an extra selection type

Definition at line 181 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

StringCutObjectSelector<reco::JetID>* TopDiLeptonOffline::MonitorEnsemble::jetIDSelect_

private

extra jetID selection on calo jets

Definition at line 183 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDGetTokenT<edm::View<reco::Jet> > TopDiLeptonOffline::MonitorEnsemble::jets_

private

input sources for monitoring

Definition at line 136 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::string TopDiLeptonOffline::MonitorEnsemble::jetSelect_

private

extra selection on jets (here given as std::string as it depends on the the jet type, which selections are valid and which not)

Definition at line 186 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::string TopDiLeptonOffline::MonitorEnsemble::label_

private

instance label

Definition at line 134 of file TopDiLeptonOfflineDQM.h.

Referenced by Mixins._Labelable::_findDependencies(), book(), Mixins._TypedParameterizable::dumpPython(), and Mixins._Labelable::setLabel().

double TopDiLeptonOffline::MonitorEnsemble::lowerEdge_

private

mass window upper and lower edge

Definition at line 188 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::vector<edm::EDGetTokenT<edm::View<reco::MET> > > TopDiLeptonOffline::MonitorEnsemble::mets_

private

considers a vector of METs

Definition at line 141 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

StringCutObjectSelector<reco::PFCandidate, true>* TopDiLeptonOffline::MonitorEnsemble::muonIso_

private

extra isolation criterion on muon

Definition at line 173 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDGetTokenT<edm::View<reco::PFCandidate> > TopDiLeptonOffline::MonitorEnsemble::muons_

private

Definition at line 137 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

StringCutObjectSelector<reco::PFCandidate, true>* TopDiLeptonOffline::MonitorEnsemble::muonSelect_

private

extra selection on muons

Definition at line 176 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

DQMStore* TopDiLeptonOffline::MonitorEnsemble::store_

private

storage manager

Definition at line 193 of file TopDiLeptonOfflineDQM.h.

Referenced by book().

edm::EDGetTokenT<edm::TriggerResults> TopDiLeptonOffline::MonitorEnsemble::triggerTable_

private

trigger table

Definition at line 144 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

double TopDiLeptonOffline::MonitorEnsemble::upperEdge_

private

Definition at line 188 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

Level TopDiLeptonOffline::MonitorEnsemble::verbosity_

private

verbosity level for booking

Definition at line 128 of file TopDiLeptonOfflineDQM.h.

Referenced by book(), and MonitorEnsemble().