SingleTopTChannelLepton_miniAOD::MonitorEnsemble Class Reference

#include <SingleTopTChannelLeptonDQM_miniAOD.h>

Public Types
enum	Level { STANDARD, VERBOSE, DEBUG }
	different verbosity levels More...

Public Member Functions
void	book (DQMStore::IBooker &ibooker)
	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)
	default contructor More...
	~MonitorEnsemble ()
	default destructor More...

Private Member Functions
bool	booked (const std::string histName) const
	check if histogram was booked 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...
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
edm::EDGetTokenT< reco::JetTagCollection >	btagCSV_
double	btagCSVWP_
edm::EDGetTokenT< reco::JetTagCollection >	btagEff_
	btag discriminator labels More...
double	btagEffWP_
	btag working points More...
edm::EDGetTokenT< reco::JetTagCollection >	btagPur_
double	btagPurWP_
edm::EDGetTokenT< reco::JetTagCollection >	btagVtx_
double	btagVtxWP_
std::string	directory_
double	eidCutValue_
std::unique_ptr< StringCutObjectSelector< pat::Electron > >	elecIso_
	extra isolation criterion on electron More...
edm::EDGetTokenT< edm::View< pat::Electron > >	elecs_
std::unique_ptr< StringCutObjectSelector< pat::Electron > >	elecSelect_
	extra selection on electrons More...
edm::EDGetTokenT< edm::ValueMap< float > >	electronId_
	electronId label More...
std::map< std::string, MonitorElement * >	hists_
	histogram container More...
bool	includeBTag_
std::string	jetCorrector_
	jetCorrector More...
edm::EDGetTokenT< reco::JetIDValueMap >	jetIDLabel_
	jetID as an extra selection type More...
std::unique_ptr< StringCutObjectSelector< reco::JetID > >	jetIDSelect_
	extra jetID selection on calo jets More...
edm::EDGetTokenT< edm::View< pat::Jet > >	jets_
	input sources for monitoring More...
std::unique_ptr< StringCutObjectSelector< pat::Jet > >	jetSelect
std::string	jetSelect_
std::string	label_
	instance label More...
int	logged_
	number of logged interesting events More...
double	lowerEdge_
	mass window upper and lower edge More...
std::vector< edm::EDGetTokenT< edm::View< pat::MET > > >	mets_
	considers a vector of METs More...
std::unique_ptr< StringCutObjectSelector< pat::Muon > >	muonIso_
	extra isolation criterion on muon More...
edm::EDGetTokenT< edm::View< pat::Muon > >	muons_
std::unique_ptr< StringCutObjectSelector< pat::Muon > >	muonSelect_
	extra selection on muons More...
edm::EDGetTokenT< edm::View< reco::Vertex > >	pvs_
std::unique_ptr< StringCutObjectSelector< reco::Vertex > >	pvSelect_
edm::InputTag	rhoTag
edm::EDConsumerBase	tmpConsumerBase
std::vector< std::string >	triggerPaths_
edm::EDGetTokenT< edm::TriggerResults >	triggerTable_
	trigger table More...
double	upperEdge_
Level	verbosity_
	verbosity level for booking More...

Detailed Description

Definition at line 32 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Member Enumeration Documentation

enum SingleTopTChannelLepton_miniAOD::MonitorEnsemble::Level

different verbosity levels

Enumerator
STANDARD
VERBOSE
DEBUG

Definition at line 35 of file SingleTopTChannelLeptonDQM_miniAOD.h.

             {
    STANDARD,
    VERBOSE,
    DEBUG
  };

Constructor & Destructor Documentation

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

default contructor

Definition at line 29 of file SingleTopTChannelLeptonDQM_miniAOD.cc.

References DEBUG, directory_, eidCutValue_, elecIso_, elecs_, elecSelect_, electronId_, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), jetCorrector_, jetIDLabel_, jetIDSelect_, jets_, jetSelect, jetSelect_, lowerEdge_, mets_, TtFullHadSignalSelMVATrainTreeSaver_cff::monitoring, muonIso_, muons_, muonSelect_, pvs_, pvSelect_, rhoTag, STANDARD, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, triggerPaths_, triggerTable_, upperEdge_, VERBOSE, and verbosity_.

    : label_(label),
      elecIso_(nullptr),
      elecSelect_(nullptr),
      pvSelect_(nullptr),
      muonIso_(nullptr),
      muonSelect_(nullptr),
      jetIDSelect_(nullptr),
      jetSelect(nullptr),
      includeBTag_(false),
      lowerEdge_(-1.),
      upperEdge_(-1.),
      logged_(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"));


  muons_ = iC.consumes<edm::View<pat::Muon> >(
    sources.getParameter<edm::InputTag>("muons"));


  elecs_ = iC.consumes<edm::View<pat::Electron> >(
      sources.getParameter<edm::InputTag>("elecs"));
  pvs_ = iC.consumes<edm::View<reco::Vertex> >(
      sources.getParameter<edm::InputTag>("pvs"));
  jets_ = iC.consumes<edm::View<pat::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<pat::MET> >(tag));
  // electronExtras 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_.reset( new StringCutObjectSelector<pat::Electron>(
          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_.reset( new StringCutObjectSelector<pat::Electron>(
          elecExtras.getParameter<std::string>("isolation")));
    }
      
    if (elecExtras.existsAs<std::string>("rho")) {
        rhoTag = elecExtras.getParameter<edm::InputTag>("rho");
    }
    // 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");
    }
  }
  // pvExtras are opetional; they may be omitted or empty
  if (cfg.existsAs<edm::ParameterSet>("pvExtras")) {
    edm::ParameterSet pvExtras =
        cfg.getParameter<edm::ParameterSet>("pvExtras");
    // select is optional; in case it's not found no
    // selection will be applied
    if (pvExtras.existsAs<std::string>("select")) {
      pvSelect_.reset(new StringCutObjectSelector<reco::Vertex>(
          pvExtras.getParameter<std::string>("select")));
    }
  }
  // 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_.reset(new StringCutObjectSelector<pat::Muon>(
          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_.reset(new StringCutObjectSelector<pat::Muon>(
          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_.reset(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");
      jetSelect.reset(new StringCutObjectSelector<pat::Jet> (jetSelect_));
    }
  }

  // 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"));
    triggerPaths_ =
        triggerExtras.getParameter<std::vector<std::string> >("paths");
  }

  // 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
  directory_ = cfg.getParameter<std::string>("directory");
  // book(ibooker);
}

SingleTopTChannelLepton_miniAOD::MonitorEnsemble::~MonitorEnsemble ( )

inline

default destructor

Definition at line 46 of file SingleTopTChannelLeptonDQM_miniAOD.h.

References book(), fill(), and GeneralSetup::setup().

{};

Member Function Documentation

void MonitorEnsemble::book

(

DQMStore::IBooker &

ibooker

)

book histograms in subdirectory directory

Definition at line 194 of file SingleTopTChannelLeptonDQM_miniAOD.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), directory_, hists_, label_, DQMStore::IBooker::setCurrentFolder(), STANDARD, AlCaHLTBitMon_QueryRunRegistry::string, triggerBinLabels(), triggerPaths_, VERBOSE, and verbosity_.

Referenced by ~MonitorEnsemble().

                                                    {
  // set up the current directory path
  std::string current(directory_);
  current += label_;
  ibooker.setCurrentFolder(current);

  // determine number of bins for trigger monitoring
  //unsigned int nPaths = triggerPaths_.size();

  // --- [STANDARD] --- //
  // Run Number
  //hists_["RunNumb_"] = ibooker.book1D("RunNumber", "Run Nr.", 1.e4, 1.5e5, 3.e5);
  // instantaneous luminosity
  //hists_["InstLumi_"] = ibooker.book1D("InstLumi", "Inst. Lumi.", 100, 0., 1.e3);
  // number of selected primary vertices
  hists_["pvMult_"] = ibooker.book1D("PvMult", "N_{good pvs}", 50, 0., 50.);
  // pt of the leading muon
  hists_["muonPt_"] = ibooker.book1D("MuonPt", "pt(#mu TightId, TightIso)", 40, 0., 200.);
  // muon multiplicity before std isolation
  hists_["muonMult_"] = ibooker.book1D("MuonMult", "N_{loose}(#mu)", 10, 0., 10.);
  // muon multiplicity after  std isolation
  //hists_["muonMultIso_"] = ibooker.book1D("MuonMultIso",
  //    "N_{TightIso}(#mu)", 10, 0., 10.);
 
  hists_["muonMultTight_"] = ibooker.book1D("MuonMultTight",
      "N_{TightIso,TightId}(#mu)", 10, 0., 10.);

 // pt of the leading electron
  hists_["elecPt_"] = ibooker.book1D("ElecPt", "pt(e TightId, TightIso)", 40, 0., 200.);
  // electron multiplicity before std isolation
  //hists_["elecMult_"] = ibooker.book1D("ElecMult", "N_{looseId}(e)", 10, 0., 10.);
  // electron multiplicity after  std isolation
  //hists_["elecMultIso_"] = ibooker.book1D("ElecMultIso", "N_{Iso}(e)", 10, 0., 10.);
  // multiplicity of jets with pt>20 (corrected to L2+L3)
  hists_["jetMult_"] = ibooker.book1D("JetMult", "N_{30}(jet)", 10, 0., 10.);
  hists_["jetLooseMult_"] = ibooker.book1D("JetLooseMult", "N_{30,loose}(jet)", 10, 0., 10.);

  // trigger efficiency estimates for single lepton triggers
  //hists_["triggerEff_"] = ibooker.book1D("TriggerEff",
  //    "Eff(trigger)", nPaths, 0., nPaths);
  // monitored trigger occupancy for single lepton triggers
  //hists_["triggerMon_"] = ibooker.book1D("TriggerMon",
  //    "Mon(trigger)", nPaths, 0., nPaths);
  // MET (calo)
  hists_["slimmedMETs_"] = ibooker.book1D("slimmedMETs", "MET_{slimmed}", 40, 0., 200.);
  // W mass estimate
  hists_["massW_"] = ibooker.book1D("MassW", "M(W)", 60, 0., 300.);
  // Top mass estimate
  hists_["massTop_"] = ibooker.book1D("MassTop", "M(Top)", 50, 0., 500.);
  // b-tagged Top mass
  //hists_["massBTop_"] = ibooker.book1D("MassBTop", "M(Top, 1 b-tag)", 50, 0., 500.);
  
    // W mass transverse estimate mu
    hists_["MTWm_"] = ibooker.book1D("MTWm", "M_{T}^{W}(#mu)", 60, 0., 300.);
    // Top mass transverse estimate mu
    hists_["mMTT_"] = ibooker.book1D("mMTT", "M_{T}^{t}(#mu)", 50, 0., 500.);
    
    // W mass transverse estimate e
    hists_["MTWe_"] = ibooker.book1D("MTWe", "M_{T}^{W}(e)", 60, 0., 300.);
    // Top mass transverse estimate e
    hists_["eMTT_"] = ibooker.book1D("eMTT", "M_{T}^{t}(e)", 50, 0., 500.);
    
    
    
  // set bin labels for trigger monitoring
  triggerBinLabels(std::string("trigger"), triggerPaths_);

  if (verbosity_ == STANDARD) return;

  // --- [VERBOSE] --- //
  // eta of the leading muon
  hists_["muonEta_"] = ibooker.book1D("MuonEta", "#eta(#mu TightId,TightIso)", 30, -3., 3.);
  // relative isolation of the candidate muon (depending on the decay channel)
  hists_["muonPhi_"] = ibooker.book1D("MuonPhi", "#phi(#mu TightId,TightIso)", 40, -4., 4.);
  hists_["muonRelIso_"] = ibooker.book1D(
      "MuonRelIso", "Iso_{Rel}(#mu TightId) (#Delta#beta Corrected)", 50, 0., 1.);
    
  // eta of the leading electron
  hists_["elecEta_"] = ibooker.book1D("ElecEta", "#eta(e TightId, TightIso)", 30, -3., 3.);
  hists_["elecPhi_"] = ibooker.book1D("ElecPhi", "#phi(e TightId, TightIso)", 40, -4., 4.);
  // std isolation variable of the leading electron
  hists_["elecRelIso_"] = ibooker.book1D("ElecRelIso", "Iso_{Rel}(e TightId)", 50, 0., 1.);
    
  hists_["elecMultTight_"] = ibooker.book1D("ElecMultTight",
                                              "N_{TightIso,TightId}(e)", 10, 0., 10.);
  
    
  // multiplicity of btagged jets (for track counting high efficiency) with
  // pt(L2L3)>30
  //hists_["jetMultBEff_"] = ibooker.book1D("JetMultBEff",
  //    "N_{30}(TCHE)", 10, 0., 10.);
  // btag discriminator for track counting high efficiency for jets with
  // pt(L2L3)>30
  //hists_["jetBDiscEff_"] = ibooker.book1D("JetBDiscEff",
  //    "Disc_{TCHE}(jet)", 100, 0., 10.);
  // eta of the 1. leading jet (corrected to L2+L3)
  hists_["jet1Eta_"] = ibooker.book1D("Jet1Eta", "#eta_{30,loose}(jet1)", 60, -3., 3.);
  // pt of the 1. leading jet (corrected to L2+L3)
  hists_["jet1Pt_"] = ibooker.book1D("Jet1Pt", "pt_{30,loose}(jet1)", 60, 0., 300.);
  // eta of the 2. leading jet (corrected to L2+L3)
  hists_["jet2Eta_"] = ibooker.book1D("Jet2Eta", "#eta_{30,loose}(jet2)", 60, -3., 3.);
  // pt of the 2. leading jet (corrected to L2+L3)
  hists_["jet2Pt_"] = ibooker.book1D("Jet2Pt", "pt_{30,loose}(jet2)", 60, 0., 300.);
  // eta of the 3. leading jet (corrected to L2+L3)
  //hists_["jet3Eta_"] = ibooker.book1D("Jet3Eta", "#eta_{30,loose}(jet3)", 60, -3., 3.);
  // pt of the 3. leading jet (corrected to L2+L3)
  //hists_["jet3Pt_"] = ibooker.book1D("Jet3Pt", "pt_{30,loose}(jet3)", 60, 0., 300.);
  // eta of the 4. leading jet (corrected to L2+L3)
  //hists_["jet4Eta_"] = ibooker.book1D("Jet4Eta", "#eta_{30,loose}(jet4)", 60, -3., 3.);
  // pt of the 4. leading jet (corrected to L2+L3)
  //hists_["jet4Pt_"] = ibooker.book1D("Jet4Pt", "pt_{30,loose}(jet4)", 60, 0., 300.);
  // MET (tc)
  hists_["slimmedMETsNoHF_"] = ibooker.book1D("slimmedMETsNoHF", "MET_{slimmedNoHF}", 40, 0., 200.);
  // MET (pflow)
  hists_["slimmedMETsPuppi_"] = ibooker.book1D("slimmedMETsPuppi", "MET_{slimmedPuppi}", 40, 0., 200.);
  // dz for muons (to suppress cosmis)
  hists_["muonDelZ_"] = ibooker.book1D("MuonDelZ", "d_{z}(#mu)", 50, -25., 25.);
  // dxy for muons (to suppress cosmics)
  hists_["muonDelXY_"] = ibooker.book2D("MuonDelXY",
      "d_{xy}(#mu)", 50, -0.1, 0.1, 50, -0.1, 0.1);

  // set axes titles for dxy for muons
  hists_["muonDelXY_"]->setAxisTitle("x [cm]", 1);
  hists_["muonDelXY_"]->setAxisTitle("y [cm]", 2);

  if (verbosity_ == VERBOSE) return;

  // --- [DEBUG] --- //
  // charged hadron isolation component of the candidate muon (depending on the
  // decay channel)
  hists_["muonChHadIso_"] = ibooker.book1D("MuonChHadIsoComp",
      "ChHad_{IsoComponent}(#mu TightId)", 50, 0., 5.);
  // neutral hadron isolation component of the candidate muon (depending on the
  // decay channel)
  hists_["muonNeHadIso_"] = ibooker.book1D("MuonNeHadIsoComp",
      "NeHad_{IsoComponent}(#mu TightId)", 50, 0., 5.);
  // photon isolation component of the candidate muon (depending on the decay
  // channel)
  hists_["muonPhIso_"] = ibooker.book1D("MuonPhIsoComp",
      "Photon_{IsoComponent}(#mu TightId)", 50, 0., 5.);
  // charged hadron isolation component of the candidate electron (depending on
  // the decay channel)
  hists_["elecChHadIso_"] = ibooker.book1D("ElectronChHadIsoComp",
      "ChHad_{IsoComponent}(e TightId)", 50, 0., 5.);
  // neutral hadron isolation component of the candidate electron (depending on
  // the decay channel)
  hists_["elecNeHadIso_"] = ibooker.book1D("ElectronNeHadIsoComp",
      "NeHad_{IsoComponent}(e TightId)", 50, 0., 5.);
  // photon isolation component of the candidate electron (depending on the
  // decay channel)
  hists_["elecPhIso_"] = ibooker.book1D("ElectronPhIsoComp",
      "Photon_{IsoComponent}(e TightId)", 50, 0., 5.);
  // multiplicity of btagged jets (for track counting high purity) with
  // pt(L2L3)>30
  //hists_["jetMultBPur_"] = ibooker.book1D("JetMultBPur",
  //    "N_{30}(TCHP)", 10, 0., 10.);
  // btag discriminator for track counting high purity
  //hists_["jetBDiscPur_"] = ibooker.book1D("JetBDiscPur",
  //    "Disc_{TCHP}(Jet)", 100, 0., 10.);
  // multiplicity of btagged jets (for simple secondary vertex) with pt(L2L3)>30
  //hists_["jetMultBVtx_"] = ibooker.book1D("JetMultBVtx",
  //    "N_{30}(SSVHE)", 10, 0., 10.);
  // btag discriminator for simple secondary vertex
  //hists_["jetBDiscVtx_"] = ibooker.book1D("JetBDiscVtx",
  //    "Disc_{SSVHE}(Jet)", 35, -1., 6.);
  // multiplicity for combined secondary vertex
  hists_["jetMultBCSVM_"] = ibooker.book1D("JetMultBCSVM", "N_{30}(CSVM)", 10, 0., 10.);
  // btag discriminator for combined secondary vertex
  hists_["jetBCSV_"] = ibooker.book1D("JetDiscCSV","BJet Disc_{CSV}(JET)", 100, -1., 2.);
  // pt of the 1. leading jet (uncorrected)
  //hists_["jet1PtRaw_"] = ibooker.book1D("Jet1PtRaw", "pt_{Raw}(jet1)", 60, 0., 300.);
  // pt of the 2. leading jet (uncorrected)
  //hists_["jet2PtRaw_"] = ibooker.book1D("Jet2PtRaw", "pt_{Raw}(jet2)", 60, 0., 300.);
  // pt of the 3. leading jet (uncorrected)
  //hists_["jet3PtRaw_"] = ibooker.book1D("Jet3PtRaw", "pt_{Raw}(jet3)", 60, 0., 300.);
  // pt of the 4. leading jet (uncorrected)
  //hists_["jet4PtRaw_"] = ibooker.book1D("Jet4PtRaw", "pt_{Raw}(jet4)", 60, 0., 300.);
  // selected events
  hists_["eventLogger_"] = ibooker.book2D("EventLogger",
      "Logged Events", 9, 0., 9., 10, 0., 10.);

  // set axes titles for selected events
  hists_["eventLogger_"]->getTH1()->SetOption("TEXT");
  hists_["eventLogger_"]->setBinLabel(1, "Run", 1);
  hists_["eventLogger_"]->setBinLabel(2, "Block", 1);
  hists_["eventLogger_"]->setBinLabel(3, "Event", 1);
  hists_["eventLogger_"]->setBinLabel(4, "pt_{L2L3}(jet1)", 1);
  hists_["eventLogger_"]->setBinLabel(5, "pt_{L2L3}(jet2)", 1);
  hists_["eventLogger_"]->setBinLabel(6, "pt_{L2L3}(jet3)", 1);
  hists_["eventLogger_"]->setBinLabel(7, "pt_{L2L3}(jet4)", 1);
  hists_["eventLogger_"]->setBinLabel(8, "M_{W}", 1);
  hists_["eventLogger_"]->setBinLabel(9, "M_{Top}", 1);
  hists_["eventLogger_"]->setAxisTitle("logged evts", 2);
  return;
}

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

inlineprivate

check if histogram was booked

Definition at line 73 of file SingleTopTChannelLeptonDQM_miniAOD.h.

References hists_.

Referenced by fill().

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

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

fill monitor histograms with electronId and jetCorrections

Definition at line 390 of file SingleTopTChannelLeptonDQM_miniAOD.cc.

References funct::abs(), pat::Jet::bDiscriminator(), pat::Jet::chargedEmEnergyFraction(), pat::Jet::chargedHadronEnergyFraction(), pat::Jet::chargedMultiplicity(), MillePedeFileConverter_cfg::e, PatTopSelectionAnalyzer_cfi::elecs, elecs_, elecSelect_, electronId_, reco::LeafCandidate::eta(), edm::EventAuxiliary::event(), edm::Event::eventAuxiliary(), objects.autophobj::float, edm::Event::getByToken(), hists_, training_settings::idx, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), metsig::jet, fwrapper::jets, jets_, jetSelect, jetSelect_, logged_, lowerEdge_, edm::EventAuxiliary::luminosityBlock(), SiStripPI::max, RazorAnalyzer::met, mets_, RPCpg::mu, metsig::muon, extraflags_cff::muons, muons_, muonSelect_, pat::Jet::neutralEmEnergyFraction(), pat::Jet::neutralHadronEnergyFraction(), pat::Jet::neutralMultiplicity(), reco::Vertex::position(), EnergyCorrector::pt, reco::LeafCandidate::pt(), MetAnalyzer::pv(), pvs_, rhoTag, edm::EventAuxiliary::run(), triggerObjects_cff::sel, triggerPaths_, triggerTable_, upperEdge_, and pseudoTop_cfi::wMass.

Referenced by fill(), selectionPath(), and ~MonitorEnsemble().

                                                       {
  // fetch trigger event if configured such
  edm::Handle<edm::TriggerResults> triggerTable;

  if (!triggerTable_.isUninitialized()) {
    if (!event.getByToken(triggerTable_, triggerTable)) return;
  }
    
  /*
  ------------------------------------------------------------

  Primary Vertex Monitoring

  ------------------------------------------------------------
  */
  // fill monitoring plots for primary verices
  edm::Handle<edm::View<reco::Vertex> > pvs;
  if (!event.getByToken(pvs_, pvs)) return;
  const reco::Vertex& pver= pvs->front();

  unsigned int pvMult = 0; 
  if(pvs.isValid()){
    for (edm::View<reco::Vertex>::const_iterator pv = pvs->begin();
       pv != pvs->end(); ++pv) {

      bool isGood = ( !(pv->isFake()) &&
              (pv->ndof() > 4.0) &&
              (abs(pv->z()) < 24.0) &&
              (abs(pv->position().Rho()) < 2.0) 
              );
      if( !isGood ) continue;
      pvMult++;
    }
    //std::cout<<" npv  "<<testn<<endl;
  }    
  
  fill("pvMult_", pvMult);

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

  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);
    
  /*
  ------------------------------------------------------------

  Electron Monitoring

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

  // fill monitoring plots for electrons
  edm::Handle<edm::View<pat::Electron> > elecs;
  if (!event.getByToken(elecs_, elecs)) return;
    
  edm::Handle< double > _rhoHandle;
  event.getByLabel(rhoTag,_rhoHandle);
  //if (!event.getByToken(elecs_, elecs)) return;
    

  // check availability of electron id
  edm::Handle<edm::ValueMap<float> > electronId;
  if (!electronId_.isUninitialized()) {
    if (!event.getByToken(electronId_, electronId)) return;
  }

  // loop electron collection
  unsigned int eMultIso = 0, eMult = 0;
  std::vector<const pat::Electron*> isoElecs;

  pat::Electron e;
    
  for (edm::View<pat::Electron>::const_iterator elec = elecs->begin();
       elec != elecs->end(); ++elec) {

      if(true){//loose id
      if (!elecSelect_ || (*elecSelect_)(*elec)) {
          
        double el_ChHadIso = elec->pfIsolationVariables().sumChargedHadronPt;
        double el_NeHadIso = elec->pfIsolationVariables().sumNeutralHadronEt;
        double el_PhIso = elec->pfIsolationVariables().sumPhotonEt;
          
        double rho = _rhoHandle.isValid() ? (float)(*_rhoHandle) : 0;
        double absEta = abs(elec->superCluster()->eta());
        double eA = 0;
        if      (absEta < 1.000) eA = 0.1703;
        else if (absEta < 1.479) eA = 0.1715;
        else if (absEta < 2.000) eA = 0.1213;
        else if (absEta < 2.200) eA = 0.1230;
        else if (absEta < 2.300) eA = 0.1635;
        else if (absEta < 2.400) eA = 0.1937;
        else if (absEta < 5.000) eA = 0.2393;
          
          
        double el_pfRelIso = (el_ChHadIso + max(0., el_NeHadIso + el_PhIso - rho * eA)) /elec->pt();
          
        ++eMult;
          
        if(eMult==1){
          fill("elecRelIso_", el_pfRelIso);
          fill("elecChHadIso_", el_ChHadIso);
          fill("elecNeHadIso_", el_NeHadIso);
          fill("elecPhIso_", el_PhIso);
        }
        //loose Iso
        //if(!((el_pfRelIso<0.0994 && absEta<1.479)||(el_pfRelIso<0.107 && absEta>1.479)))continue;
        
        //tight Iso
        if(!((el_pfRelIso<0.0588 && absEta<1.479)||(el_pfRelIso<0.0571 && absEta>1.479)))continue;
        ++eMultIso;
        
          
        if (eMultIso == 1) {
          // restrict to the leading electron
          e = *elec;
            
          fill("elecPt_", elec->pt());
          fill("elecEta_", elec->eta());
          fill("elecPhi_", elec->phi());
        }
      }
    }
  }
  //fill("elecMult_", eMult);
  fill("elecMultTight_", eMultIso);
    
    
    
    

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

  Muon Monitoring

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

  // fill monitoring plots for muons
  unsigned int mMult = 0,  mTight=0;

  edm::Handle<edm::View<pat::Muon> > muons;
    
  pat::Muon mu;
    
  edm::View<pat::Muon>::const_iterator muonit;

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

  for (edm::View<pat::Muon>::const_iterator muon = muons->begin();
       muon != muons->end(); ++muon) {


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



      // apply preselection loose muon
      if (!muonSelect_ || (*muonSelect_)(*muon)) {
          
        //loose muon count
        ++mMult;

        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!

        if(!(muon->isGlobalMuon() && muon->isPFMuon() && muon->globalTrack()->normalizedChi2() < 10. && muon->globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && muon->numberOfMatchedStations() > 1 && fabs(muon->muonBestTrack()->dxy(pver.position())) < 0.2  && fabs(muon->muonBestTrack()->dz(pver.position())) < 0.5 && muon->innerTrack()->hitPattern().numberOfValidPixelHits() > 0 && muon->innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5   ) )continue;
          
        if (mMult == 1) {
              // restrict to leading muon
          fill("muonRelIso_", pfRelIso);
          fill("muonChHadIso_", chHadPt);
          fill("muonNeHadIso_", neHadEt);
          fill("muonPhIso_", phoEt);
          fill("muonRelIso_",pfRelIso);
            
        }
          
        if(!(pfRelIso<0.15))continue;
          
        ++mTight;
          
          
        //tight id
        if (mTight == 1) {
          // restrict to leading muon
          mu = *muon;
          fill("muonPt_", muon->pt());
          fill("muonEta_", muon->eta());
          fill("muonPhi_", muon->phi());
            
        }
      }
    }
  }
  fill("muonMult_", mMult); //loose
  fill("muonMultTight_", mTight); //tight id & iso
    
  /*
  ------------------------------------------------------------

  Jet Monitoring

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

  

  // loop jet collection
  std::vector<pat::Jet> correctedJets;
  std::vector<double> JetTagValues;
  pat::Jet TaggedJetCand;
  vector<double> bJetDiscVal;
  
  unsigned int mult = 0, loosemult=0, multBCSVM = 0;

  edm::Handle<edm::View<pat::Jet> > jets;
  if (!event.getByToken(jets_, jets)) {
    return;
  }

  for (edm::View<pat::Jet>::const_iterator jet = jets->begin();
       jet != jets->end(); ++jet) {
    // check jetID for calo jets
    //unsigned int idx = jet - jets->begin();

      pat::Jet sel = *jet;

      if ( jetSelect==nullptr)
    jetSelect.reset(new StringCutObjectSelector<pat::Jet>(jetSelect_));
   
      if(!(*jetSelect)(sel))continue;
//      if (!jetSelect(sel)) continue;

    // prepare jet to fill monitor histograms
    pat::Jet monitorJet = *jet;


    ++mult;
      
    if(monitorJet.chargedHadronEnergyFraction()>0 && monitorJet.chargedMultiplicity()>0 && monitorJet.chargedEmEnergyFraction()<0.99 && monitorJet.neutralHadronEnergyFraction()<0.99 && monitorJet.neutralEmEnergyFraction()<0.99 && (monitorJet.chargedMultiplicity()+monitorJet.neutralMultiplicity())>1 ){

      correctedJets.push_back(monitorJet);
      ++loosemult;  // determine jet multiplicity
    
      fill("jetBCSV_", monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags")); //hard coded discriminator and value right now.
      if (monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags") > 0.89){
            
        if (multBCSVM == 0){
           TaggedJetCand = monitorJet;
       bJetDiscVal.push_back(monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags"));
        }else if(multBCSVM == 1){
       bJetDiscVal.push_back(monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags"));
           if (bJetDiscVal[1] > bJetDiscVal[0]) TaggedJetCand = monitorJet;
        }
          
        ++multBCSVM;
      }

        

      // Fill a vector with Jet b-tag WP for later M3+1tag calculation: CSV
      // tagger
      JetTagValues.push_back(monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags"));
      //    }
      // fill pt (raw or L2L3) for the leading four jets
      if (loosemult == 1) {
      //cout<<" jet id= "<<monitorJet.chargedHadronEnergyFraction()<<endl;

        fill("jet1Pt_", monitorJet.pt());
        //fill("jet1PtRaw_", jet->pt());
        fill("jet1Eta_", monitorJet.eta());
      };
      if (loosemult == 2) {
        fill("jet2Pt_", monitorJet.pt());
        //fill("jet2PtRaw_", jet->pt());
        fill("jet2Eta_", monitorJet.eta());
      }
      
    }
  }
  fill("jetMult_", mult);
  fill("jetLooseMult_", loosemult);
  fill("jetMultBCSVM_", multBCSVM);

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

  MET Monitoring

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

  // fill monitoring histograms for met
    
  pat::MET mET;
    
  for (std::vector<edm::EDGetTokenT<edm::View<pat::MET> > >::const_iterator
           met_ = mets_.begin();
       met_ != mets_.end(); ++met_) {
    edm::Handle<edm::View<pat::MET> > met;
    if (!event.getByToken(*met_, met)) continue;
    if (met->begin() != met->end()) {
      unsigned int idx = met_ - mets_.begin();
      if (idx == 0) fill("slimmedMETs_", met->begin()->et());
      if (idx == 1) fill("slimmedMETsNoHF_", met->begin()->et());
      if (idx == 2) fill("slimmedMETsPuppi_", met->begin()->et());
    }
  }

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

  Event Monitoring

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

  // fill W boson and top mass estimates
    
  Calculate_miniAOD eventKinematics(MAXJETS, WMASS);
  double wMass = eventKinematics.massWBoson(correctedJets);
  double topMass = eventKinematics.massTopQuark(correctedJets);
  if (wMass >= 0 && topMass >= 0) {
    fill("massW_", wMass);
    fill("massTop_", topMass);
  }

  // Fill M3 with Btag (CSV Tight) requirement

 // if (!includeBTag_) return;
  //if (correctedJets.size() != JetTagValues.size()) return;
  //double btopMass =
  //    eventKinematics.massBTopQuark(correctedJets, JetTagValues, 0.89); //hard coded CSVv2 value
  
  //if (btopMass >= 0) fill("massBTop_", btopMass);
    
  // fill plots for trigger monitoring
  if ((lowerEdge_ == -1. && upperEdge_ == -1.) ||
    (lowerEdge_ < wMass && wMass < upperEdge_)) {
    if (!triggerTable_.isUninitialized())
      fill(event, *triggerTable, "trigger", triggerPaths_);
    if (logged_ <= hists_.find("eventLogger_")->second->getNbinsY()) {
      // log runnumber, lumi block, event number & some
      // more pysics infomation for interesting events
      fill("eventLogger_", 0.5, logged_ + 0.5, event.eventAuxiliary().run());
      fill("eventLogger_", 1.5, logged_ + 0.5,event.eventAuxiliary().luminosityBlock());
      fill("eventLogger_", 2.5, logged_ + 0.5, event.eventAuxiliary().event());
      if (correctedJets.size() > 0)
        fill("eventLogger_", 3.5, logged_ + 0.5, correctedJets[0].pt());
      if (correctedJets.size() > 1)
        fill("eventLogger_", 4.5, logged_ + 0.5, correctedJets[1].pt());
      if (correctedJets.size() > 2)
        fill("eventLogger_", 5.5, logged_ + 0.5, correctedJets[2].pt());
      if (correctedJets.size() > 3)
        fill("eventLogger_", 6.5, logged_ + 0.5, correctedJets[3].pt());

      fill("eventLogger_", 7.5, logged_ + 0.5, wMass);
      fill("eventLogger_", 8.5, logged_ + 0.5, topMass);
      ++logged_;
    }
  }
    
    if (multBCSVM != 0 && mTight == 1) {
        
        double mtW = eventKinematics.tmassWBoson(&mu, mET, TaggedJetCand);
        fill("MTWm_", mtW);
        double MTT = eventKinematics.tmassTopQuark(&mu, mET, TaggedJetCand);
        fill("mMTT_", MTT);
    }
    
    if (multBCSVM != 0 && eMultIso == 1) {
        double mtW = eventKinematics.tmassWBoson(&e, mET, TaggedJetCand);
        fill("MTWe_", mtW);
        double MTT = eventKinematics.tmassTopQuark(&e, mET, TaggedJetCand);
        fill("eMTT_", MTT);
    }
  
    
    

}

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 174 of file SingleTopTChannelLeptonDQM_miniAOD.h.

References accept(), fill(), hists_, training_settings::idx, monitorPath(), and selectionPath().

                                                                           {
  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 SingleTopTChannelLepton_miniAOD::MonitorEnsemble::fill ( const std::string  histName, double  value  ) const

inlineprivate

fill histogram if it had been booked before

Definition at line 77 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

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

void SingleTopTChannelLepton_miniAOD::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 82 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

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

void SingleTopTChannelLepton_miniAOD::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 87 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

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

std::string SingleTopTChannelLepton_miniAOD::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 56 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and triggerBinLabels().

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

std::string SingleTopTChannelLepton_miniAOD::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 61 of file SingleTopTChannelLeptonDQM_miniAOD.h.

References fill(), tablePrinter::labels, AlCaHLTBitMon_QueryRunRegistry::string, and triggerBinLabels().

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 162 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

Referenced by book(), and selectionPath().

                                                          {
  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

edm::EDGetTokenT<reco::JetTagCollection> SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagCSV_

private

Definition at line 146 of file SingleTopTChannelLeptonDQM_miniAOD.h.

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagCSVWP_

private

Definition at line 149 of file SingleTopTChannelLeptonDQM_miniAOD.h.

edm::EDGetTokenT<reco::JetTagCollection> SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagEff_

private

btag discriminator labels

Definition at line 146 of file SingleTopTChannelLeptonDQM_miniAOD.h.

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagEffWP_

private

btag working points

Definition at line 149 of file SingleTopTChannelLeptonDQM_miniAOD.h.

edm::EDGetTokenT<reco::JetTagCollection> SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagPur_

private

Definition at line 146 of file SingleTopTChannelLeptonDQM_miniAOD.h.

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagPurWP_

private

Definition at line 149 of file SingleTopTChannelLeptonDQM_miniAOD.h.

edm::EDGetTokenT<reco::JetTagCollection> SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagVtx_

private

Definition at line 146 of file SingleTopTChannelLeptonDQM_miniAOD.h.

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagVtxWP_

private

Definition at line 149 of file SingleTopTChannelLeptonDQM_miniAOD.h.

std::string SingleTopTChannelLepton_miniAOD::MonitorEnsemble::directory_

private

Definition at line 159 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by book(), and MonitorEnsemble().

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::eidCutValue_

private

Definition at line 116 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<pat::Electron> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecIso_

private

extra isolation criterion on electron

Definition at line 118 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

edm::EDGetTokenT<edm::View<pat::Electron> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecs_

private

Definition at line 103 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<pat::Electron> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecSelect_

private

extra selection on electrons

Definition at line 120 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

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

private

electronId label

Definition at line 114 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

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

private

histogram container

Definition at line 157 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

bool SingleTopTChannelLepton_miniAOD::MonitorEnsemble::includeBTag_

private

include btag information or not to be determined from the cfg

Definition at line 144 of file SingleTopTChannelLeptonDQM_miniAOD.h.

std::string SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetCorrector_

private

jetCorrector

Definition at line 133 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

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

private

jetID as an extra selection type

Definition at line 135 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<reco::JetID> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetIDSelect_

private

extra jetID selection on calo jets

Definition at line 137 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

edm::EDGetTokenT<edm::View<pat::Jet> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jets_

private

input sources for monitoring

Definition at line 101 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<pat::Jet> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetSelect

private

Definition at line 141 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::string SingleTopTChannelLepton_miniAOD::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 140 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::string SingleTopTChannelLepton_miniAOD::MonitorEnsemble::label_

private

instance label

Definition at line 97 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), book(), and Mixins._Labelable::setLabel().

int SingleTopTChannelLepton_miniAOD::MonitorEnsemble::logged_

private

number of logged interesting events

Definition at line 154 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill().

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::lowerEdge_

private

mass window upper and lower edge

Definition at line 151 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::vector<edm::EDGetTokenT<edm::View<pat::MET> > > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::mets_

private

considers a vector of METs

Definition at line 99 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<pat::Muon> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muonIso_

private

extra isolation criterion on muon

Definition at line 127 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

edm::EDGetTokenT<edm::View<pat::Muon> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muons_

private

Definition at line 102 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<pat::Muon> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muonSelect_

private

extra selection on muons

Definition at line 130 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDGetTokenT<edm::View<reco::Vertex> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::pvs_

private

Definition at line 104 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<reco::Vertex> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::pvSelect_

private

extra selection on primary vertices; meant to investigate the pile-up effect

Definition at line 124 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

edm::InputTag SingleTopTChannelLepton_miniAOD::MonitorEnsemble::rhoTag

private

Definition at line 111 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDConsumerBase SingleTopTChannelLepton_miniAOD::MonitorEnsemble::tmpConsumerBase

private

Definition at line 158 of file SingleTopTChannelLeptonDQM_miniAOD.h.

std::vector<std::string> SingleTopTChannelLepton_miniAOD::MonitorEnsemble::triggerPaths_

private

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

Definition at line 109 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

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

private

trigger table

Definition at line 106 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::upperEdge_

private

Definition at line 151 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

Level SingleTopTChannelLepton_miniAOD::MonitorEnsemble::verbosity_

private

verbosity level for booking

Definition at line 91 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by book(), and MonitorEnsemble().