#include <TopDiLeptonOfflineDQM.h>

Public Types
enum	DecayChannel { NONE, DIMUON, DIELEC, ELECMU, NONE, DIMUON, DIELEC, ELECMU }
	different decay channels More...
enum	DecayChannel { NONE, DIMUON, DIELEC, ELECMU, NONE, DIMUON, DIELEC, ELECMU }
	different decay channels More...
enum	Level { STANDARD, VERBOSE, DEBUG, STANDARD, VERBOSE, DEBUG }
	different verbosity levels More...
enum	Level { STANDARD, VERBOSE, DEBUG, STANDARD, VERBOSE, DEBUG }
	different verbosity levels More...
typedef reco::LeafCandidate::LorentzVector	LorentzVector
typedef reco::LeafCandidate::LorentzVector	LorentzVector
Public Member Functions
void	book (std::string directory)
	book histograms in subdirectory _directory_
void	book (std::string directory)
	book histograms in subdirectory _directory_
void	fill (const edm::Event &event, const edm::EventSetup &setup)
	fill monitor histograms with electronId and jetCorrections
void	fill (const edm::Event &event, const edm::EventSetup &setup)
	fill monitor histograms with electronId and jetCorrections
	MonitorEnsemble (const char *label, const edm::ParameterSet &cfg)
	default contructor
	MonitorEnsemble (const char *label, const edm::ParameterSet &cfg)
	default contructor
	~MonitorEnsemble ()
	default destructor
	~MonitorEnsemble ()
	default destructor
Private Member Functions
bool	booked (const std::string histName) const
	check if histogram was booked
bool	booked (const std::string histName) const
	check if histogram was booked
DecayChannel	decayChannel (const std::vector< const reco::Muon * > &muons, const std::vector< const reco::GsfElectron * > &elecs) const
	determine dileptonic decay channel
DecayChannel	decayChannel (const std::vector< const reco::Muon * > &muons, const std::vector< const reco::GsfElectron * > &elecs) const
	determine dileptonic decay channel
void	fill (const std::string histName, double value) const
	fill histogram if it had been booked before
void	fill (const std::string histName, double xValue, double yValue) const
	fill histogram if it had been booked before (2-dim version)
void	fill (const edm::Event &event, const edm::TriggerResults &triggerTable, std::string channel, const std::vector< std::string > labels) const
	fill trigger monitoring histograms
void	fill (const std::string histName, double value) const
	fill histogram if it had been booked before
void	fill (const std::string histName, double xValue, double yValue, double zValue) const
	fill histogram if it had been booked before (2-dim version)
void	fill (const std::string histName, double xValue, double yValue, double zValue) const
	fill histogram if it had been booked before (2-dim version)
void	fill (const std::string histName, double xValue, double yValue) const
	fill histogram if it had been booked before (2-dim version)
void	fill (const edm::Event &event, const edm::TriggerResults &triggerTable, std::string channel, const std::vector< std::string > labels) const
	fill trigger monitoring histograms
void	loggerBinLabels (std::string hist)
	set labels for event logging histograms
void	loggerBinLabels (std::string hist)
	set labels for event logging histograms
std::string	monitorPath (const std::string &label) const
std::string	monitorPath (const std::string &label) const
std::string	selectionPath (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
void	triggerBinLabels (std::string channel, const std::vector< std::string > labels)
	set configurable labels for trigger monitoring histograms
Private Attributes
int	diElecLogged_
int	diMuonLogged_
std::vector< std::string >	diMuonPaths_
	trigger paths for di muon channel
int	eidPattern_
StringCutObjectSelector < reco::GsfElectron > *	elecIso_
	extra isolation criterion on electron
int	elecMuLogged_
	number of logged interesting events
std::vector< std::string >	elecMuPaths_
edm::InputTag	elecs_
	input sources for monitoring
StringCutObjectSelector < reco::GsfElectron > *	elecSelect_
	extra selection on electrons
edm::InputTag	electronId_
	electronId label
std::map< std::string, MonitorElement * >	hists_
	histogram container
std::string	jetCorrector_
	jetCorrector
edm::InputTag	jetIDLabel_
	jetID as an extra selection type
StringCutObjectSelector < reco::JetID > *	jetIDSelect_
	extra jetID selection on calo jets
edm::InputTag	jets_
std::string	jetSelect_
std::string	label_
	instance label
double	lowerEdge_
	mass window upper and lower edge
std::vector< edm::InputTag >	mets_
	considers a vector of METs
StringCutObjectSelector < reco::Muon > *	muonIso_
	extra isolation criterion on muon
edm::InputTag	muons_
StringCutObjectSelector < reco::Muon > *	muonSelect_
	extra selection on muons
DQMStore *	store_
	storage manager
edm::InputTag	triggerTable_
	trigger table
double	upperEdge_
Level	verbosity_
	verbosity level for booking

Detailed Description

Definition at line 39 of file TopDiLeptonOfflineDQM.h.

Member Typedef Documentation

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

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

Definition at line 45 of file TopDiLeptonOfflineDQM.h.

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

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

Definition at line 46 of file TopHLTDiLeptonOfflineDQM.h.

Member Enumeration Documentation

enum TopDiLeptonOffline::MonitorEnsemble::DecayChannel

different decay channels

Enumerator:

NONE
DIMUON
DIELEC
ELECMU
NONE
DIMUON
DIELEC
ELECMU

Definition at line 47 of file TopDiLeptonOfflineDQM.h.

{ NONE, DIMUON, DIELEC, ELECMU };

enum TopDiLeptonOffline::MonitorEnsemble::DecayChannel

different decay channels

Enumerator:

NONE
DIMUON
DIELEC
ELECMU
NONE
DIMUON
DIELEC
ELECMU

Definition at line 48 of file TopHLTDiLeptonOfflineDQM.h.

{ NONE, DIMUON, DIELEC, ELECMU };

enum TopDiLeptonOffline::MonitorEnsemble::Level

different verbosity levels

Enumerator:

STANDARD
VERBOSE
DEBUG
STANDARD
VERBOSE
DEBUG

Definition at line 43 of file TopHLTDiLeptonOfflineDQM.h.

{ STANDARD, VERBOSE, DEBUG };

enum TopDiLeptonOffline::MonitorEnsemble::Level

different verbosity levels

Enumerator:

STANDARD
VERBOSE
DEBUG
STANDARD
VERBOSE
DEBUG

Definition at line 42 of file TopDiLeptonOfflineDQM.h.

{ STANDARD, VERBOSE, DEBUG };

Constructor & Destructor Documentation

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

default contructor

Definition at line 11 of file TopDiLeptonOfflineDQM.cc.

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

                                                                                : 
   label_(label), eidPattern_(0), elecIso_(0), elecSelect_(0), muonIso_(0), muonSelect_(0), jetIDSelect_(0), 
   lowerEdge_(-1.), upperEdge_(-1.), elecMuLogged_(0), diMuonLogged_(0), diElecLogged_(0)
  {
    // sources have to be given; this PSet is not optional
    edm::ParameterSet sources=cfg.getParameter<edm::ParameterSet>("sources");
    muons_= sources.getParameter<edm::InputTag>("muons");
    elecs_= sources.getParameter<edm::InputTag>("elecs");
    jets_ = sources.getParameter<edm::InputTag>("jets" );
    mets_ = sources.getParameter<std::vector<edm::InputTag> >("mets" );

    // elecExtras are optional; they may be omitted or empty
    if( cfg.existsAs<edm::ParameterSet>("elecExtras") ){
      edm::ParameterSet elecExtras=cfg.getParameter<edm::ParameterSet>("elecExtras");
      // select is optional; in case it's not found no
      // selection will be applied
      if( elecExtras.existsAs<std::string>("select") ){
        elecSelect_= new StringCutObjectSelector<reco::GsfElectron>(elecExtras.getParameter<std::string>("select"));
      }
      // isolation is optional; in case it's not found no
      // isolation will be applied
      if( elecExtras.existsAs<std::string>("isolation") ){
        elecIso_= new StringCutObjectSelector<reco::GsfElectron>(elecExtras.getParameter<std::string>("isolation"));
      }
      // electronId is optional; in case it's not found the 
      // InputTag will remain empty
      if( elecExtras.existsAs<edm::ParameterSet>("electronId") ){
        edm::ParameterSet elecId=elecExtras.getParameter<edm::ParameterSet>("electronId");
        electronId_= elecId.getParameter<edm::InputTag>("src");
        eidPattern_= elecId.getParameter<int>("pattern");
      }
    }
    // muonExtras are optional; they may be omitted or empty
    if( cfg.existsAs<edm::ParameterSet>("muonExtras") ){
      edm::ParameterSet muonExtras=cfg.getParameter<edm::ParameterSet>("muonExtras");
      // select is optional; in case it's not found no
      // selection will be applied
      if( muonExtras.existsAs<std::string>("select") ){
        muonSelect_= new StringCutObjectSelector<reco::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_= new StringCutObjectSelector<reco::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_ =jetID.getParameter<edm::InputTag>("label");
        jetIDSelect_= new StringCutObjectSelector<reco::JetID>(jetID.getParameter<std::string>("select"));
      }
      // select is optional; in case it's not found no
      // selection will be applied (only implemented for 
      // CaloJets at the moment)
      if( jetExtras.existsAs<std::string>("select") ){
        jetSelect_= jetExtras.getParameter<std::string>("select");
      }
    }
    // triggerExtras are optional; they may be omitted or empty
    if( cfg.existsAs<edm::ParameterSet>("triggerExtras") ){
      edm::ParameterSet triggerExtras=cfg.getParameter<edm::ParameterSet>("triggerExtras");
      triggerTable_=triggerExtras.getParameter<edm::InputTag>("src");
      elecMuPaths_ =triggerExtras.getParameter<std::vector<std::string> >("pathsELECMU");
      diMuonPaths_ =triggerExtras.getParameter<std::vector<std::string> >("pathsDIMUON");
    }
    // massExtras is optional; in case it's not found no mass
    // window cuts are applied for the same flavor monitor
    // histograms
    if( cfg.existsAs<edm::ParameterSet>("massExtras") ){
      edm::ParameterSet massExtras=cfg.getParameter<edm::ParameterSet>("massExtras");
      lowerEdge_= massExtras.getParameter<double>("lowerEdge");
      upperEdge_= massExtras.getParameter<double>("upperEdge");
    }

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

TopDiLeptonOffline::MonitorEnsemble::~MonitorEnsemble ( ) [inline]

default destructor

Definition at line 53 of file TopDiLeptonOfflineDQM.h.

{};

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

default contructor

TopDiLeptonOffline::MonitorEnsemble::~MonitorEnsemble ( ) [inline]

default destructor

Definition at line 54 of file TopHLTDiLeptonOfflineDQM.h.

{};

Member Function Documentation

void MonitorEnsemble::book ( std::string directory )

book histograms in subdirectory _directory_

Definition at line 113 of file TopDiLeptonOfflineDQM.cc.

References DQMStore::book1D(), DQMStore::book2D(), cond::rpcobimon::current, diMuonPaths_, elecMuPaths_, hists_, label_, loggerBinLabels(), cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), STANDARD, store_, triggerBinLabels(), VERBOSE, and verbosity_.

Referenced by MonitorEnsemble().

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

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

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

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

    if( verbosity_==STANDARD) return;

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

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

    if( verbosity_==VERBOSE) return;

    // --- [DEBUG] --- //
    // electron multiplicity after std isolation
    hists_["elecMultIso_" ] = store_->book1D("ElecMultIso" , "N_{Iso}(e)"              ,       10,   0.,      10.);
    // muon multiplicity after std isolation
    hists_["muonMultIso_" ] = store_->book1D("MuonMultIso" , "N_{Iso}(#mu)"            ,       10,   0.,      10.);
    // calo isolation of the candidate muon (depending on the decay channel)
    hists_["muonCalIso_"  ] = store_->book1D("MuonCalIso"  , "Iso_{Cal}(#mu)"          ,       50,   0.,       1.);
    // track isolation of the candidate muon (depending on the decay channel)
    hists_["muonTrkIso_"  ] = store_->book1D("MuonTrkIso"  , "Iso_{Trk}(#mu)"          ,       50,   0.,       1.);
    // calo isolation of the candidate electron (depending on the decay channel)
    hists_["elecCalIso_"  ] = store_->book1D("ElecCalIso"  , "Iso_{Cal}(e)"            ,       50,   0.,       1.);
    // track isolation of the candidate electron (depending on the decay channel)
    hists_["elecTrkIso_"  ] = store_->book1D("ElecTrkIso"  , "Iso_{Trk}(e)"            ,       50,   0.,       1.);
    // eta of the leading jet
    hists_["jet1Eta_"     ] = store_->book1D("Jet1Eta"     , "#eta(jet1)"              ,       30,  -5.,       5.); 
    // eta of the 2. leading jet
    hists_["jet2Eta_"     ] = store_->book1D("Jet2Eta"     , "#eta(jet2)"              ,       30,  -5.,       5.);
    // pt of the 1. leading jet (not corrected)
    hists_["jet1PtRaw_"   ] = store_->book1D("Jet1PtRaw"   , "pt_{Raw}(jet1)"          ,       60,   0.,     300.);   
    // pt of the 2. leading jet (not corrected)     
    hists_["jet2PtRaw_"   ] = store_->book1D("Jet2PtRaw"   , "pt_{Raw}(jet2)"          ,       60,   0.,     300.);
    // deltaEta between the 2 leading jets
    hists_["dEtaJet1Jet2_"] = store_->book1D("DEtaJet1Jet2", "#Delta#eta(jet1,jet2)"   ,       30,   0.,       3.);
    // deltaEta between the lepton and the leading jet
    hists_["dEtaJet1Lep1_"] = store_->book1D("DEtaJet1Lep1", "#Delta#eta(jet1,lep1)"   ,       30,   0.,       3.);
    // deltaEta between the lepton and MET
    hists_["dEtaLep1MET_" ] = store_->book1D("DEtaLep1MET" , "#Delta#eta(lep1,MET)"    ,       30,   0.,       3.);
    // deltaEta between leading jet and MET
    hists_["dEtaJet1MET_" ] = store_->book1D("DEtaJet1MET" , "#Delta#eta(jet1,MET)"    ,       30,   0.,       3.);
    // deltaPhi of 2 leading jets
    hists_["dPhiJet1Jet2_"] = store_->book1D("DPhiJet1Jet2", "#Delta#phi(jet1,jet2)"   ,       32,   0.,      3.2);
    // deltaPhi of 1. lepton and 1. jet
    hists_["dPhiJet1Lep1_"] = store_->book1D("DPhiJet1Lep1", "#Delta#phi(jet1,lep1)"   ,       32,   0.,      3.2);
    // deltaPhi of 1. lepton and MET
    hists_["dPhiLep1MET_" ] = store_->book1D("DPhiLep1MET" , "#Delta#phi(lep1,MET)"    ,       32,   0.,      3.2);
    // deltaPhi of 1. jet and MET
    hists_["dPhiJet1MET_" ] = store_->book1D("DPhiJet1MET" , "#Delta#phi(jet1,MET)"    ,       32,   0.,      3.2);
    // selected dimuon events
    hists_["diMuonLogger_"] = store_->book2D("DiMuonLogger", "Logged DiMuon Events"    ,        8,   0.,       8.,   10,   0.,   10.);
    // selected dielec events
    hists_["diElecLogger_"] = store_->book2D("DiElecLogger", "Logged DiElec Events"    ,        8,   0.,       8.,   10,   0.,   10.);
    // selected elemu events
    hists_["elecMuLogger_"] = store_->book2D("ElecMuLogger", "Logged ElecMu Events"    ,        8,   0.,       8.,   10,   0.,   10.);

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

void TopDiLeptonOffline::MonitorEnsemble::book ( std::string directory )

book histograms in subdirectory _directory_

bool TopDiLeptonOffline::MonitorEnsemble::booked ( const std::string histName ) const [inline, private]

check if histogram was booked

Definition at line 78 of file TopDiLeptonOfflineDQM.h.

References hists_.

Referenced by fill().

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

bool TopDiLeptonOffline::MonitorEnsemble::booked ( const std::string histName ) const [inline, private]

check if histogram was booked

Definition at line 79 of file TopHLTDiLeptonOfflineDQM.h.

References hists_.

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

MonitorEnsemble::DecayChannel MonitorEnsemble::decayChannel	(	const std::vector< const reco::Muon * > &	muons,
		const std::vector< const reco::GsfElectron * > &	elecs
	)		const `[inline, private]`

determine dileptonic decay channel

Definition at line 198 of file TopDiLeptonOfflineDQM.h.

References DIELEC, DIMUON, ELECMU, and NONE.

Referenced by fill().

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

DecayChannel TopDiLeptonOffline::MonitorEnsemble::decayChannel	(	const std::vector< const reco::Muon * > &	muons,
		const std::vector< const reco::GsfElectron * > &	elecs
	)		const `[private]`

determine dileptonic decay channel

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

fill trigger monitoring histograms

void TopDiLeptonOffline::MonitorEnsemble::fill	(	const std::string	histName,
		double	value
	)		const `[inline, private]`

fill histogram if it had been booked before

Definition at line 80 of file TopDiLeptonOfflineDQM.h.

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

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

void TopDiLeptonOffline::MonitorEnsemble::fill	(	const std::string	histName,
		double	value
	)		const `[inline, private]`

fill histogram if it had been booked before

Definition at line 81 of file TopHLTDiLeptonOfflineDQM.h.

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

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

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

fill trigger monitoring histograms

Definition at line 184 of file TopDiLeptonOfflineDQM.h.

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

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

void TopDiLeptonOffline::MonitorEnsemble::fill	(	const std::string	histName,
		double	xValue,
		double	yValue
	)		const `[inline, private]`

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

Definition at line 83 of file TopHLTDiLeptonOfflineDQM.h.

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

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

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

fill monitor histograms with electronId and jetCorrections

Definition at line 251 of file TopDiLeptonOfflineDQM.cc.

References DeDxDiscriminatorTools::charge(), decayChannel(), reco::deltaPhi(), reco::deltaR(), DIELEC, diElecLogged_, DIMUON, diMuonLogged_, diMuonPaths_, eidPattern_, elecIso_, ELECMU, elecMuLogged_, elecMuPaths_, elecs_, elecSelect_, electronId_, reco::LeafCandidate::et(), eta(), edm::Event::eventAuxiliary(), edm::EventSetup::find(), edm::Event::getByLabel(), JetCorrector::getJetCorrector(), hists_, UserOptions_cff::idx, metsig::jet, jetCorrector_, jetIDLabel_, jetIDSelect_, analyzePatCleaning_cfg::jets, jets_, jetSelect_, edm::InputTag::label(), lowerEdge_, scaleCards::mass, CaloMET_cfi::met, mets_, VarParsing::mult, metsig::muon, muonIso_, patZpeak::muons, muons_, muonSelect_, muon::overlap(), phi, reco::LeafCandidate::pt(), reco::Jet::scaleEnergy(), EgammaValidation_Wenu_cff::sel, triggerTable_, and upperEdge_.

Referenced by fill().

  {
    // fetch trigger event if configured such 
    edm::Handle<edm::TriggerResults> triggerTable;
    if(!triggerTable_.label().empty()) {
      if( !event.getByLabel(triggerTable_, triggerTable) ) return;
    }

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

    Muon Selection

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

    // buffer isolated muons
    std::vector<const reco::Muon*> isoMuons;

    edm::Handle<edm::View<reco::Muon> > muons;
    if( !event.getByLabel(muons_, muons) ) return;

    for(edm::View<reco::Muon>::const_iterator muon=muons->begin(); muon!=muons->end(); ++muon){
      // restrict to globalMuons
      if( muon->isGlobalMuon() ){ 
        fill("muonDelZ_" , muon->globalTrack()->vz());
        fill("muonDelXY_", muon->globalTrack()->vx(), muon->globalTrack()->vy());
        // apply preselection
        if(!muonSelect_ || (*muonSelect_)(*muon)){
          double isolationTrk = muon->pt()/(muon->pt()+muon->isolationR03().sumPt);
          double isolationCal = muon->pt()/(muon->pt()+muon->isolationR03().emEt+muon->isolationR03().hadEt);
          double isolationRel = (muon->isolationR03().sumPt+muon->isolationR03().emEt+muon->isolationR03().hadEt)/muon->pt();
          fill("muonTrkIso_" , isolationTrk); fill("muonCalIso_" , isolationCal); fill("muonRelIso_" , isolationRel);
          if(!muonIso_ || (*muonIso_)(*muon)) isoMuons.push_back(&(*muon));
        }
      }
    }
    fill("muonMultIso_", isoMuons.size());

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

    Electron Selection

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

    // buffer isolated electronss
    std::vector<const reco::GsfElectron*> isoElecs;
    edm::Handle<edm::ValueMap<float> > electronId; 
    if(!electronId_.label().empty()) {
      if( !event.getByLabel(electronId_, electronId) ) return;
    }

    edm::Handle<edm::View<reco::GsfElectron> > elecs;
    if( !event.getByLabel(elecs_, elecs) ) return;

    for(edm::View<reco::GsfElectron>::const_iterator elec=elecs->begin(); elec!=elecs->end(); ++elec){
      // restrict to electrons with good electronId
      int idx = elec-elecs->begin();
      if( electronId_.label().empty() ? true : ((int)(*electronId)[elecs->refAt(idx)] & eidPattern_) ){
        // apply preselection
        if(!elecSelect_ || (*elecSelect_)(*elec)){
          double isolationTrk = elec->pt()/(elec->pt()+elec->dr03TkSumPt());
          double isolationCal = elec->pt()/(elec->pt()+elec->dr03EcalRecHitSumEt()+elec->dr03HcalTowerSumEt());
          double isolationRel = (elec->dr03TkSumPt()+elec->dr03EcalRecHitSumEt()+elec->dr03HcalTowerSumEt())/elec->pt();
          fill("elecTrkIso_" , isolationTrk); fill("elecCalIso_" , isolationCal); fill("elecRelIso_" , isolationRel);
          if(!elecIso_ || (*elecIso_)(*elec)) isoElecs.push_back(&(*elec));
        }
      }
    }
    fill("elecMultIso_", isoElecs.size());

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

    Jet Selection

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

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

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

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

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

    MET Selection

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

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

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

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


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

    Event Monitoring

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

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

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

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

void TopDiLeptonOffline::MonitorEnsemble::fill	(	const std::string	histName,
		double	xValue,
		double	yValue
	)		const `[inline, private]`

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

Definition at line 82 of file TopDiLeptonOfflineDQM.h.

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

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

void TopDiLeptonOffline::MonitorEnsemble::fill	(	const std::string	histName,
		double	xValue,
		double	yValue,
		double	zValue
	)		const `[inline, private]`

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

Definition at line 84 of file TopDiLeptonOfflineDQM.h.

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

void TopDiLeptonOffline::MonitorEnsemble::fill	(	const std::string	histName,
		double	xValue,
		double	yValue,
		double	zValue
	)		const `[inline, private]`

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

Definition at line 85 of file TopHLTDiLeptonOfflineDQM.h.

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

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

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

fill monitor histograms with electronId and jetCorrections

void MonitorEnsemble::loggerBinLabels ( std::string hist ) [inline, private]

set labels for event logging histograms

Definition at line 148 of file TopDiLeptonOfflineDQM.h.

References hists_.

Referenced by book().

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

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

void TopDiLeptonOffline::MonitorEnsemble::loggerBinLabels ( std::string hist ) [private]

set labels for event logging histograms

std::string TopDiLeptonOffline::MonitorEnsemble::monitorPath ( const std::string & label ) const [inline, private]

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

Definition at line 64 of file TopHLTDiLeptonOfflineDQM.h.

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

std::string TopDiLeptonOffline::MonitorEnsemble::monitorPath ( const std::string & label ) const [inline, private]

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

Definition at line 63 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and triggerBinLabels().

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

std::string TopDiLeptonOffline::MonitorEnsemble::selectionPath ( const std::string & label ) const [inline, private]

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

Definition at line 67 of file TopHLTDiLeptonOfflineDQM.h.

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

std::string TopDiLeptonOffline::MonitorEnsemble::selectionPath ( const std::string & label ) const [inline, private]

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

Definition at line 66 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and triggerBinLabels().

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

void TopDiLeptonOffline::MonitorEnsemble::triggerBinLabels	(	std::string	channel,
		const std::vector< std::string >	labels
	)		`[private]`

set configurable labels for trigger monitoring histograms

void MonitorEnsemble::triggerBinLabels	(	std::string	channel,
		const std::vector< std::string >	labels
	)		`[inline, private]`

set configurable labels for trigger monitoring histograms

Definition at line 175 of file TopDiLeptonOfflineDQM.h.

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

Referenced by book().

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

Member Data Documentation

int MonitorEnsemble::diElecLogged_ [private]

Definition at line 140 of file TopDiLeptonOfflineDQM.h.

Referenced by fill().

int MonitorEnsemble::diMuonLogged_ [private]

Definition at line 140 of file TopDiLeptonOfflineDQM.h.

Referenced by fill().

std::vector< std::string > MonitorEnsemble::diMuonPaths_ [private]

trigger paths for di muon channel

Definition at line 102 of file TopDiLeptonOfflineDQM.h.

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

int MonitorEnsemble::eidPattern_ [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 116 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

StringCutObjectSelector< reco::GsfElectron > * MonitorEnsemble::elecIso_ [private]

extra isolation criterion on electron

Definition at line 118 of file TopDiLeptonOfflineDQM.h.

Referenced by fill(), and MonitorEnsemble().

int MonitorEnsemble::elecMuLogged_ [private]

number of logged interesting events

Definition at line 140 of file TopDiLeptonOfflineDQM.h.

Referenced by fill().

std::vector< std::string > MonitorEnsemble::elecMuPaths_ [private]

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