SingleTopTChannelLepton::MonitorEnsemble Class Reference

#include <SingleTopTChannelLeptonDQM.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, const edm::VParameterSet &vcfg, 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 >	btagCombVtx_
double	btagCombVtxWP_
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 < reco::PFCandidate, true > >	elecIso
std::string	elecIso_
	extra isolation criterion on electron More...
edm::EDGetTokenT< edm::View < reco::PFCandidate > >	elecs_
edm::EDGetTokenT< edm::View < reco::GsfElectron > >	elecs_gsf_
std::unique_ptr < StringCutObjectSelector < reco::PFCandidate, true > >	elecSelect
std::unique_ptr < StringCutObjectSelector < reco::PFCandidate > >	elecSelect_
	extra selection on electrons More...
edm::EDGetTokenT < edm::ValueMap< float > >	electronId_
	electronId label More...
std::map< std::string, MonitorElement * >	hists_
bool	includeBTag_
edm::ESGetToken< JetCorrector, JetCorrectionsRecord >	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...
std::unique_ptr < StringCutObjectSelector < reco::PFJet > >	jetlooseSelection_
edm::EDGetTokenT< edm::View < reco::Jet > >	jets_
	input sources for monitoring More...
std::string	jetSelect_
std::unique_ptr < StringCutObjectSelector < reco::PFJet > >	jetSelection_
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< reco::MET > > >	mets_
	considers a vector of METs More...
std::unique_ptr < StringCutObjectSelector < reco::PFCandidate, true > >	muonIso
std::unique_ptr < StringCutObjectSelector < reco::PFCandidate > >	muonIso_
	extra isolation criterion on muon More...
edm::EDGetTokenT< edm::View < reco::PFCandidate > >	muons_
std::unique_ptr < StringCutObjectSelector < reco::PFCandidate, true > >	muonSelect
std::unique_ptr < StringCutObjectSelector < reco::PFCandidate > >	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 55 of file SingleTopTChannelLeptonDQM.h.

Member Enumeration Documentation

enum SingleTopTChannelLepton::MonitorEnsemble::Level

different verbosity levels

Enumerator
STANDARD
VERBOSE
DEBUG

Definition at line 58 of file SingleTopTChannelLeptonDQM.h.

{ STANDARD, VERBOSE, DEBUG };

Constructor & Destructor Documentation

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

default contructor

Definition at line 26 of file SingleTopTChannelLeptonDQM.cc.

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

      : label_(label),
        elecSelect_(nullptr),
        pvSelect_(nullptr),
        muonIso_(nullptr),
        muonSelect_(nullptr),
        jetIDSelect_(nullptr),
        jetlooseSelection_(nullptr),
        jetSelection_(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"));
    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));
    pvs_ = iC.consumes<edm::View<reco::Vertex>>(sources.getParameter<edm::InputTag>("pvs"));
    // electronExtras are optional; they may be omitted or
    // empty
    if (cfg.existsAs<edm::ParameterSet>("elecExtras")) {
      // rho for PF isolation with EA corrections
      // eventrhoToken_ =
      // iC.consumes<double>(edm::InputTag("fixedGridRhoFastjetAll"));

      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_ = std::make_unique<StringCutObjectSelector<reco::PFCandidate>>(
            elecExtras.getParameter<std::string>("select"));
      }

      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 optional; 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_ =
            std::make_unique<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_ = std::make_unique<StringCutObjectSelector<reco::PFCandidate>>(
            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_ = std::make_unique<StringCutObjectSelector<reco::PFCandidate>>(
            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_ = iC.esConsumes(edm::ESInputTag("", 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_ =
            std::make_unique<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");
        jetSelection_ = std::make_unique<StringCutObjectSelector<reco::PFJet>>(jetSelect_);
        jetlooseSelection_ = std::make_unique<StringCutObjectSelector<reco::PFJet>>(
            "chargedHadronEnergyFraction()>0 && chargedMultiplicity()>0 && chargedEmEnergyFraction()<0.99 && "
            "neutralHadronEnergyFraction()<0.99 && neutralEmEnergyFraction()<0.99 && "
            "(chargedMultiplicity()+neutralMultiplicity())>1");
      }

      // jetBDiscriminators are optional; in case they are
      // not found the InputTag will remain empty; they
      // consist of pairs of edm::JetFlavorAssociation's &
      // corresponding working points
      includeBTag_ = jetExtras.existsAs<edm::ParameterSet>("jetBTaggers");
      if (includeBTag_) {
        edm::ParameterSet btagCSV =
            jetExtras.getParameter<edm::ParameterSet>("jetBTaggers").getParameter<edm::ParameterSet>("cvsVertex");
        btagCSV_ = iC.consumes<reco::JetTagCollection>(btagCSV.getParameter<edm::InputTag>("label"));
        btagCSVWP_ = btagCSV.getParameter<double>("workingPoint");
      }
    }

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

SingleTopTChannelLepton::MonitorEnsemble::~MonitorEnsemble ( )

inline

default destructor

Definition at line 67 of file SingleTopTChannelLeptonDQM.h.

{}

Member Function Documentation

void MonitorEnsemble::book

(

DQMStore::IBooker &

ibooker

)

book histograms in subdirectory directory

Definition at line 178 of file SingleTopTChannelLeptonDQM.cc.

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

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

    // --- [STANDARD] --- //
    // 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 tight Id and tight Iso
    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.);
    // multiplicity of jets with pt>30 (corrected to L1+L2+L3)
    hists_["jetMult_"] = ibooker.book1D("JetMult", "N_{30}(jet)", 10, 0., 10.);
    // multiplicity of loose Id jets with pt>30
    hists_["jetLooseMult_"] = ibooker.book1D("JetLooseMult", "N_{30,loose}(jet)", 10, 0., 10.);

    // MET (pflow)
    hists_["metPflow_"] = ibooker.book1D("METPflow", "MET_{Pflow}", 50, 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.);
    // 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_["muonRelIso_"] = ibooker.book1D("MuonRelIso", "Iso_{Rel}(#mu TightId) (#Delta#beta Corrected)", 50, 0., 1.);
    // phi of the leading muon
    hists_["muonPhi_"] = ibooker.book1D("MuonPhi", "#phi(#mu TightId, TightIso)", 40, -4., 4.);
    // eta of the leading electron
    hists_["elecEta_"] = ibooker.book1D("ElecEta", "#eta(e tightId, TightIso)", 30, -3., 3.);
    // std isolation variable of the leading electron
    hists_["elecRelIso_"] = ibooker.book1D("ElecRelIso", "Iso_{Rel}(e TightId)", 50, 0., 1.);
    // phi of the leading electron
    hists_["elecPhi_"] = ibooker.book1D("ElecPhi", "#phi(e tightId, TightIso)", 40, -4., 4.);
    // multiplicity of tight Id, tight Iso electorns
    hists_["elecMultTight_"] = ibooker.book1D("ElecMultTight", "N_{TightIso,TightId}(e)", 10, 0., 10.);

    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.);

    // 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 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_{30,loose}(jet1)", 1);
    hists_["eventLogger_"]->setBinLabel(5, "pt_{30,loose}(jet2)", 1);
    hists_["eventLogger_"]->setBinLabel(6, "pt_{30,loose}(jet3)", 1);
    hists_["eventLogger_"]->setBinLabel(7, "pt_{30,loose}(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::MonitorEnsemble::booked ( const std::string  histName ) const

inlineprivate

check if histogram was booked

Definition at line 91 of file SingleTopTChannelLeptonDQM.h.

References hists_.

Referenced by fill().

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

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

fill monitor histograms with electronId and jetCorrections

Definition at line 309 of file SingleTopTChannelLeptonDQM.cc.

References funct::abs(), btagCSV_, btagCSVWP_, JetCorrector::correction(), pfClustersFromHGC3DClusters_cfi::corrector, alignCSCRings::e, eidCutValue_, elecs_, elecSelect_, electronId_, reco::LeafCandidate::eta(), edm::EventAuxiliary::event(), edm::Event::eventAuxiliary(), edm::EventSetup::find(), edm::Event::getByToken(), edm::EventSetup::getData(), edm::ESGetToken< T, R >::hasValidIndex(), hists_, includeBTag_, edm::ESGetToken< T, R >::isInitialized(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), metsig::jet, jetCorrector_, jetlooseSelection_, fwrapper::jets, jets_, jetSelection_, logged_, lowerEdge_, edm::EventAuxiliary::luminosityBlock(), Calculate::massTopQuark(), Calculate::massWBoson(), SiStripPI::max, SingleTopTChannelLepton::MAXJETS, objects.METAnalyzer::met, mets_, RPCpg::mu, VarParsing::mult, HLT_FULL_cff::muon, patZpeak::muons, muons_, muonSelect_, reco::Vertex::position(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), MetAnalyzer::pv(), pvs_, pvSelect_, rhoTag, edm::EventAuxiliary::run(), reco::Jet::scaleEnergy(), EgammaValidation_Wenu_cff::sel, Calculate::tmassTopQuark(), Calculate::tmassWBoson(), triggerPaths_, triggerTable_, upperEdge_, trackerHitRTTI::vector, and SingleTopTChannelLepton::WMASS.

Referenced by fill().

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

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

    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& Pvertex = pvs->front();
    unsigned int pvMult = 0;
    for (edm::View<reco::Vertex>::const_iterator pv = pvs->begin(); pv != pvs->end(); ++pv) {
      if (!pvSelect_ || (*pvSelect_)(*pv))
        pvMult++;
    }
    fill("pvMult_", pvMult);
    /*
     ------------------------------------------------------------

     Electron Monitoring

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

    // fill monitoring plots for electrons
    edm::Handle<edm::View<reco::PFCandidate>> elecs;
    reco::GsfElectron e;
    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 eMult = 0, eMultIso = 0;
    std::vector<const reco::PFCandidate*> isoElecs;

    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_)) {  //This Electron Id is not currently used, but we can keep this for future needs
        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 absEta = std::abs(gsf_el->superCluster()->eta());

          //Effective Area computation
          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 rho = _rhoHandle.isValid() ? (float)(*_rhoHandle) : 0;
          double el_pfRelIso = (el_ChHadIso + max(0., el_NeHadIso + el_PhIso - rho * eA)) / gsf_el->pt();

          //Only TightId
          if (eMult == 0) {  // Restricted to the leading tight electron
            fill("elecRelIso_", el_pfRelIso);
            fill("elecChHadIso_", el_ChHadIso);
            fill("elecNeHadIso_", el_NeHadIso);
            fill("elecPhIso_", el_PhIso);
          }
          ++eMult;

          if (!((el_pfRelIso < 0.0588 && absEta < 1.479) || (el_pfRelIso < 0.0571 && absEta > 1.479)))
            continue;  // PF Isolation with Effective Area Corrections according to https://twiki.cern.ch/twiki/bin/viewauth/CMS/CutBasedElectronIdentificationRun2

          // TightId and TightIso
          if (eMultIso == 0) {  //Only leading
            e = *gsf_el;
            fill("elecPt_", gsf_el->pt());
            fill("elecEta_", gsf_el->eta());
            fill("elecPhi_", gsf_el->phi());
          }
          ++eMultIso;
        }
      }
    }
    //fill("elecMult_", eMult);
    fill("elecMultTight_", eMultIso);

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

     Muon Monitoring

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

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

    edm::Handle<edm::View<reco::PFCandidate>> muons;
    edm::View<reco::PFCandidate>::const_iterator muonit;
    reco::MuonRef muon;
    reco::Muon mu;

    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();

      // 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
        if ((!muonSelect_ || (*muonSelect_)(*muonit))) {
          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 &&
                muon->innerTrack()->hitPattern().numberOfValidPixelHits() > 0 &&
                muon->innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
                fabs(muon->muonBestTrack()->dxy(Pvertex.position())) < 0.2 &&
                fabs(muon->muonBestTrack()->dz(Pvertex.position())) < 0.5))
            continue;  //Only tight muons

          if (mTightId == 0) {
            fill("muonRelIso_", pfRelIso);
            fill("muonChHadIso_", chHadPt);
            fill("muonNeHadIso_", neHadEt);
            fill("muonPhIso_", phoEt);
          }
          mTightId++;

          if (!(pfRelIso < 0.15))
            continue;  //Tight Iso

          if (mTight == 0) {  //Leading tightId tightIso muon
            mu = *(muon);
            fill("muonPt_", muon->pt());
            fill("muonEta_", muon->eta());
            fill("muonPhi_", muon->phi());
          }
          mTight++;
        }
      }
    }
    fill("muonMult_", mMult);
    fill("muonMultTight_", mTight);

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

  Jet Monitoring

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

    // check availability of the btaggers
    edm::Handle<reco::JetTagCollection> btagEff, btagPur, btagVtx, btagCSV;
    if (includeBTag_) {
      if (!event.getByToken(btagCSV_, btagCSV))
        return;
    }

    // load jet
    // corrector if configured such
    const JetCorrector* corrector = nullptr;
    if (!jetCorrector_.isInitialized() && jetCorrector_.hasValidIndex()) {
      // check whether a jet correcto is in the event setup or not
      if (setup.find(edm::eventsetup::EventSetupRecordKey::makeKey<JetCorrectionsRecord>())) {
        corrector = &setup.getData(jetCorrector_);
        ;
      } else {
        edm::LogVerbatim("SingleTopTChannelLeptonDQM")
            << "\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";
      }
    }
    // loop jet collection
    std::vector<reco::Jet> correctedJets;
    std::vector<double> JetTagValues;
    reco::Jet TaggedJetCand;
    unsigned int mult = 0, multLoose = 0, multCSV = 0;
    vector<double> bJetDiscVal;

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

    for (edm::View<reco::Jet>::const_iterator jet = jets->begin(); jet != jets->end(); ++jet) {
      bool isLoose = false;
      // check jetID for calo jets, keep functionality if we ever want to go back to those
      //    unsigned int idx = jet - jets->begin();
      //    if (dynamic_cast<const reco::CaloJet*>(&*jet)) {
      //      if (jetIDSelect_ &&
      //          dynamic_cast<const reco::CaloJet*>(jets->refAt(idx).get())) {
      //        if (!(*jetIDSelect_)((*jetID)[jets->refAt(idx)])) continue;
      //      }
      //    }

      // check additional jet selection for pf jets
      if (dynamic_cast<const reco::PFJet*>(&*jet)) {
        reco::PFJet sel = dynamic_cast<const reco::PFJet&>(*jet);
        if ((*jetlooseSelection_)(sel))
          isLoose = true;
        sel.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
        if (!(*jetSelection_)(sel))
          continue;
      }
      // prepare jet to fill monitor histograms
      reco::Jet monitorJet = *jet;

      ++mult;  // determine jet (no Id) multiplicity
      monitorJet.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);

      if (isLoose) {  //Loose Id
        unsigned int idx = jet - jets->begin();
        correctedJets.push_back(monitorJet);
        if (includeBTag_) {
          // fill b-discriminators
          edm::RefToBase<reco::Jet> jetRef = jets->refAt(idx);
          fill("jetBCSV_", (*btagCSV)[jetRef]);
          if ((*btagCSV)[jetRef] > btagCSVWP_) {
            if (multCSV == 0) {
              TaggedJetCand = monitorJet;
              bJetDiscVal.push_back((*btagCSV)[jetRef]);
            } else if (multCSV == 1) {
              bJetDiscVal.push_back((*btagCSV)[jetRef]);
              if (bJetDiscVal[1] > bJetDiscVal[0])
                TaggedJetCand = monitorJet;
            }
            ++multCSV;
          }
          JetTagValues.push_back((*btagCSV)[jetRef]);
        }

        // fill pt/eta for the leading four jets
        if (multLoose == 0) {
          fill("jet1Pt_", monitorJet.pt());
          fill("jet1Eta_", monitorJet.eta());
        };
        if (multLoose == 1) {
          fill("jet2Pt_", monitorJet.pt());
          fill("jet2Eta_", monitorJet.eta());
        }
        multLoose++;
      }
    }
    fill("jetMult_", mult);
    fill("jetMultLoose_", multLoose);
    fill("jetMultBCSVM_", multCSV);

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

  MET Monitoring

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

    // fill monitoring histograms for met
    reco::MET mET;
    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) {
          fill("metPflow_", met->begin()->et());
          mET = *(met->begin());
        }
      }
    }

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

     Event Monitoring

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

    // fill W boson and top mass estimates
    Calculate 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 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.empty())
          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 (multCSV != 0 && mTight == 1) {
      double mtW = eventKinematics.tmassWBoson(&mu, mET, TaggedJetCand);
      fill("MTWm_", mtW);
      double MTT = eventKinematics.tmassTopQuark(&mu, mET, TaggedJetCand);
      fill("mMTT_", MTT);
    }

    if (multCSV != 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 213 of file SingleTopTChannelLeptonDQM.h.

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

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

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

inlineprivate

fill histogram if it had been booked before

Definition at line 93 of file SingleTopTChannelLeptonDQM.h.

References booked(), and hists_.

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

void SingleTopTChannelLepton::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 98 of file SingleTopTChannelLeptonDQM.h.

References booked(), and hists_.

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

void SingleTopTChannelLepton::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 103 of file SingleTopTChannelLeptonDQM.h.

References booked(), and hists_.

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

std::string SingleTopTChannelLepton::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 77 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and triggerBinLabels().

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

std::string SingleTopTChannelLepton::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 80 of file SingleTopTChannelLeptonDQM.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 205 of file SingleTopTChannelLeptonDQM.h.

References hists_, monitorPath(), and selectionPath().

Referenced by book().

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

Member Data Documentation

edm::EDGetTokenT<reco::JetTagCollection> SingleTopTChannelLepton::MonitorEnsemble::btagCombVtx_

private

Definition at line 182 of file SingleTopTChannelLeptonDQM.h.

double SingleTopTChannelLepton::MonitorEnsemble::btagCombVtxWP_

private

Definition at line 185 of file SingleTopTChannelLeptonDQM.h.

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

private

Definition at line 182 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

double SingleTopTChannelLepton::MonitorEnsemble::btagCSVWP_

private

Definition at line 185 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

btag discriminator labels

Definition at line 182 of file SingleTopTChannelLeptonDQM.h.

double SingleTopTChannelLepton::MonitorEnsemble::btagEffWP_

private

btag working points

Definition at line 185 of file SingleTopTChannelLeptonDQM.h.

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

private

Definition at line 182 of file SingleTopTChannelLeptonDQM.h.

double SingleTopTChannelLepton::MonitorEnsemble::btagPurWP_

private

Definition at line 185 of file SingleTopTChannelLeptonDQM.h.

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

private

Definition at line 182 of file SingleTopTChannelLeptonDQM.h.

double SingleTopTChannelLepton::MonitorEnsemble::btagVtxWP_

private

Definition at line 185 of file SingleTopTChannelLeptonDQM.h.

std::string SingleTopTChannelLepton::MonitorEnsemble::directory_

private

Definition at line 196 of file SingleTopTChannelLeptonDQM.h.

Referenced by book(), and MonitorEnsemble().

double SingleTopTChannelLepton::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 149 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate, true> > SingleTopTChannelLepton::MonitorEnsemble::elecIso

private

Definition at line 202 of file SingleTopTChannelLeptonDQM.h.

std::string SingleTopTChannelLepton::MonitorEnsemble::elecIso_

private

extra isolation criterion on electron

Definition at line 151 of file SingleTopTChannelLeptonDQM.h.

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

private

Definition at line 120 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDGetTokenT<edm::View<reco::GsfElectron> > SingleTopTChannelLepton::MonitorEnsemble::elecs_gsf_

private

Definition at line 119 of file SingleTopTChannelLeptonDQM.h.

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate, true> > SingleTopTChannelLepton::MonitorEnsemble::elecSelect

private

Definition at line 201 of file SingleTopTChannelLeptonDQM.h.

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate> > SingleTopTChannelLepton::MonitorEnsemble::elecSelect_

private

extra selection on electrons

Definition at line 153 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

electronId label

Definition at line 135 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

storage manager histogram container

Definition at line 193 of file SingleTopTChannelLeptonDQM.h.

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

bool SingleTopTChannelLepton::MonitorEnsemble::includeBTag_

private

include btag information or not to be determined from the cfg

Definition at line 179 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::ESGetToken<JetCorrector, JetCorrectionsRecord> SingleTopTChannelLepton::MonitorEnsemble::jetCorrector_

private

jetCorrector

Definition at line 166 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

jetID as an extra selection type

Definition at line 168 of file SingleTopTChannelLeptonDQM.h.

Referenced by MonitorEnsemble().

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

private

extra jetID selection on calo jets

Definition at line 171 of file SingleTopTChannelLeptonDQM.h.

Referenced by MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<reco::PFJet> > SingleTopTChannelLepton::MonitorEnsemble::jetlooseSelection_

private

Definition at line 172 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

input sources for monitoring

Definition at line 117 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::string SingleTopTChannelLepton::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 176 of file SingleTopTChannelLeptonDQM.h.

Referenced by MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<reco::PFJet> > SingleTopTChannelLepton::MonitorEnsemble::jetSelection_

private

Definition at line 173 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::string SingleTopTChannelLepton::MonitorEnsemble::label_

private

instance label

Definition at line 112 of file SingleTopTChannelLeptonDQM.h.

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

int SingleTopTChannelLepton::MonitorEnsemble::logged_

private

number of logged interesting events

Definition at line 190 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill().

double SingleTopTChannelLepton::MonitorEnsemble::lowerEdge_

private

mass window upper and lower edge

Definition at line 187 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

considers a vector of METs

Definition at line 114 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate, true> > SingleTopTChannelLepton::MonitorEnsemble::muonIso

private

Definition at line 199 of file SingleTopTChannelLeptonDQM.h.

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate> > SingleTopTChannelLepton::MonitorEnsemble::muonIso_

private

extra isolation criterion on muon

Definition at line 161 of file SingleTopTChannelLeptonDQM.h.

Referenced by MonitorEnsemble().

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

private

Definition at line 118 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate, true> > SingleTopTChannelLepton::MonitorEnsemble::muonSelect

private

Definition at line 198 of file SingleTopTChannelLeptonDQM.h.

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate> > SingleTopTChannelLepton::MonitorEnsemble::muonSelect_

private

extra selection on muons

Definition at line 164 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

Definition at line 121 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

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

private

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

Definition at line 158 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::InputTag SingleTopTChannelLepton::MonitorEnsemble::rhoTag

private

Definition at line 154 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

edm::EDConsumerBase SingleTopTChannelLepton::MonitorEnsemble::tmpConsumerBase

private

Definition at line 194 of file SingleTopTChannelLeptonDQM.h.

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

private

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

Definition at line 131 of file SingleTopTChannelLeptonDQM.h.

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

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

private

trigger table

Definition at line 127 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

double SingleTopTChannelLepton::MonitorEnsemble::upperEdge_

private

Definition at line 187 of file SingleTopTChannelLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

Level SingleTopTChannelLepton::MonitorEnsemble::verbosity_

private

verbosity level for booking

Definition at line 106 of file SingleTopTChannelLeptonDQM.h.

Referenced by book(), and MonitorEnsemble().