#include <TopSingleLeptonHLTOfflineDQM.h>

Public Member Functions
void	book (DQMStore::IBooker &store_)
	book histograms in subdirectory directory More...

void	fill (const edm::Event &event, const edm::EventSetup &setup, const HLTConfigProvider &hltConfig, const std::vector< std::string > triggerPaths)
	fill monitor histograms with electronId and jetCorrections More...

	MonitorSingleLepton (const char *label, const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
	default contructor More...

	~MonitorSingleLepton ()
	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 >	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_

int	eidPattern_

StringCutObjectSelector < reco::GsfElectron > *	elecIso_
	extra isolation criterion on electron More...

edm::EDGetTokenT< edm::View < reco::GsfElectron > >	elecs_
	input sources for monitoring More...

StringCutObjectSelector < reco::GsfElectron > *	elecSelect_
	extra selection on electrons More...

edm::EDGetTokenT < edm::ValueMap< float > >	electronId_
	electronId label More...

trigger::Vids	electronIds_

trigger::VRelectron	electronRefs_

std::string	folder_

bool	hasRawTriggerSummary

std::map< std::string, MonitorElement * >	hists_
	histogram container More...

bool	includeBTag_

std::string	jetCorrector_
	jetCorrector More...

edm::EDGetTokenT < reco::JetIDValueMap >	jetIDLabel_
	jetID as an extra selection type More...

StringCutObjectSelector < reco::JetID > *	jetIDSelect_
	extra jetID selection on calo jets More...

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

std::string	jetSelect_

std::string	label_
	instance label More...

int	logged_
	number of logged interesting events More...

double	lowerEdge_
	mass window upper and lower edge More...

std::vector< edm::EDGetTokenT < edm::View< reco::MET > > >	mets_
	considers a vector of METs More...

trigger::Vids	muonIds_

StringCutObjectSelector < reco::Muon > *	muonIso_
	extra isolation criterion on muon More...

trigger::VRmuon	muonRefs_

edm::EDGetTokenT< edm::View < reco::Muon > >	muons_

StringCutObjectSelector < reco::Muon > *	muonSelect_
	extra selection on muons More...

trigger::Vids	pfjetIds_

trigger::VRpfjet	pfjetRefs_

std::string	processName_
	hlt objects More...

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

StringCutObjectSelector < reco::Vertex > *	pvSelect_
	extra selection on primary vertices; meant to investigate the pile-up effect More...

std::vector< std::string >	triggerPaths_

edm::EDGetTokenT < trigger::TriggerEventWithRefs >	triggerSummaryTokenAOD

edm::EDGetTokenT < trigger::TriggerEventWithRefs >	triggerSummaryTokenRAW

edm::EDGetTokenT < edm::TriggerResults >	triggerTable_
	trigger table More...

double	upperEdge_

Detailed Description

Definition at line 46 of file TopSingleLeptonHLTOfflineDQM.h.

Constructor & Destructor Documentation

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

default contructor

Definition at line 37 of file TopSingleLeptonHLTOfflineDQM.cc.

References btagEff_, btagEffWP_, btagVtx_, btagVtxWP_, eidPattern_, elecIso_, elecs_, elecSelect_, electronId_, edm::ParameterSet::existsAs(), folder_, edm::ParameterSet::getParameter(), includeBTag_, HLT_25ns10e33_v2_cff::InputTag, jetCorrector_, jetIDLabel_, jetIDSelect_, jets_, jetSelect_, lowerEdge_, objects.METAnalyzer::met, mets_, muonIso_, muons_, muonSelect_, LaserDQM_cfg::process, processName_, pvs_, pvSelect_, AlCaHLTBitMon_QueryRunRegistry::string, triggerPaths_, triggerSummaryTokenAOD, triggerSummaryTokenRAW, triggerTable_, and upperEdge_.

                                                                                                                    : 
     label_(label), elecIso_(0), elecSelect_(0), pvSelect_(0), muonIso_(0), muonSelect_(0), jetIDSelect_(0), 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::Muon> >(sources.getParameter<edm::InputTag>("muons"));
     elecs_= iC.consumes< edm::View<reco::GsfElectron> >(sources.getParameter<edm::InputTag>("elecs"));
     jets_ = iC.consumes< edm::View<reco::Jet> >(sources.getParameter<edm::InputTag>("jets" ));
 
     const auto& mets = sources.getParameter<std::vector<edm::InputTag>>("mets");
     for (const auto& met: mets) {
       mets_.push_back(iC.consumes<edm::View<reco::MET>>(met));
     }
 
     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") ){
       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_= iC.consumes< edm::ValueMap<float> >(elecId.getParameter<edm::InputTag>("src"));
         eidPattern_= elecId.getParameter<int>("pattern");
       }
     }
     // 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_= new StringCutObjectSelector<reco::Vertex>(pvExtras.getParameter<std::string>("select"));
       }
     }
     // muonExtras are optional; they may be omitted or empty
     if( cfg.existsAs<edm::ParameterSet>("muonExtras") ){
       edm::ParameterSet muonExtras=cfg.getParameter<edm::ParameterSet>("muonExtras");
       // select is optional; in case it's not found no
       // selection will be applied
       if( muonExtras.existsAs<std::string>("select") ){
         muonSelect_= 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_ = iC.consumes< reco::JetIDValueMap >(jetID.getParameter<edm::InputTag>("label"));
         jetIDSelect_= new StringCutObjectSelector<reco::JetID>(jetID.getParameter<std::string>("select"));
       }
       // select is optional; in case it's not found no
       // selection will be applied (only implemented for 
       // CaloJets at the moment)
       if( jetExtras.existsAs<std::string>("select") ){
         jetSelect_= jetExtras.getParameter<std::string>("select");
       }
       // 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 btagEff=jetExtras.getParameter<edm::ParameterSet>("jetBTaggers").getParameter<edm::ParameterSet>("trackCountingEff");
         btagEff_= iC.consumes< reco::JetTagCollection >(btagEff.getParameter<edm::InputTag>("label")); btagEffWP_= btagEff.getParameter<double>("workingPoint");
 //        edm::ParameterSet btagPur=jetExtras.getParameter<edm::ParameterSet>("jetBTaggers").getParameter<edm::ParameterSet>("trackCountingPur");
 //        btagPur_= iC.consumes< reco::JetTagCollection >(btagPur.getParameter<edm::InputTag>("label")); btagPurWP_= btagPur.getParameter<double>("workingPoint");
         edm::ParameterSet btagVtx=jetExtras.getParameter<edm::ParameterSet>("jetBTaggers").getParameter<edm::ParameterSet>("secondaryVertex" );
         btagVtx_= iC.consumes< reco::JetTagCollection >(btagVtx.getParameter<edm::InputTag>("label")); btagVtxWP_= btagVtx.getParameter<double>("workingPoint");
       }
     }
 
     // triggerExtras are optional; they may be omitted or empty
     processName_ = "HLT";
     if( cfg.existsAs<edm::ParameterSet>("triggerExtras") ){
       edm::ParameterSet triggerExtras=cfg.getParameter<edm::ParameterSet>("triggerExtras");
       triggerTable_= iC.consumes< edm::TriggerResults >(triggerExtras.getParameter<edm::InputTag>("src"));
       processName_ = triggerExtras.getParameter<edm::InputTag>("src").process();
       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");
     }
 
     // and don't forget to do the histogram booking
     folder_=cfg.getParameter<std::string>("directory");
 
 //    triggerEventWithRefsTag_ = iC.consumes< trigger::TriggerEventWithRefs >(edm::InputTag("hltTriggerSummaryRAW","",processName_));
     triggerSummaryTokenRAW = iC.consumes <trigger::TriggerEventWithRefs>(edm::InputTag("hltTriggerSummaryRAW","",processName_));
     triggerSummaryTokenAOD = iC.consumes <trigger::TriggerEventWithRefs>(edm::InputTag("hltTriggerSummaryAOD","",processName_));
 
   }

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::~MonitorSingleLepton ( )

inline

default destructor

Definition at line 53 of file TopSingleLeptonHLTOfflineDQM.h.

53 {};

Member Function Documentation

void MonitorSingleLepton::book ( DQMStore::IBooker & store_ )

book histograms in subdirectory directory

Definition at line 163 of file TopSingleLeptonHLTOfflineDQM.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), folder_, hists_, label_, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, triggerBinLabels(), and triggerPaths_.

     {
       //set up the current directory path
       std::string current(folder_); current+=label_;
       store_.setCurrentFolder(current);
 
       // determine number of bins for trigger monitoring
       unsigned int nPaths=triggerPaths_.size();
 
       // number of selected primary vertices
       hists_["pvMult_"     ] = store_.book1D("PvMult"     , "N_{pvs}"          ,     100,     0.,    100.);  
       // multiplicity of jets with pt>20 (corrected to L2+L3)
       hists_["jetMult_"    ] = store_.book1D("JetMult"    , "N_{20}(jet)"      ,     10,     0.,     10.);   
       // // trigger efficiency estimates for single lepton triggers
       // hists_["triggerEff_" ] = store_.book1D("TriggerEff" , "Eff(trigger)"     , nPaths,     0.,  nPaths);
       // monitored trigger occupancy for single lepton triggers
       hists_["triggerMon_" ] = store_.book1D("TriggerMon" , "Mon(trigger)"     , nPaths,     0.,  nPaths);
       // W mass estimate
       hists_["massW_"      ] = store_.book1D("MassW"      , "M(W)"             ,     60,     0.,    300.);   
       // Top mass estimate
       hists_["massTop_"    ] = store_.book1D("MassTop"    , "M(Top)"           ,     50,     0.,    500.);   
       // Mlb mu 
       hists_["mMub_"       ] = store_.book1D("mMub"       , "m_{#mub}"         ,     50,     0.,    500.);
       // W mass transverse estimate mu
       hists_["MTWm_"       ] = store_.book1D("MTWm"       , "M_{T}^{W}(#mu)"   ,     60,     0.,    300.);
       // Top mass transverse estimate mu
       hists_["mMTT_"       ] = store_.book1D("mMTT"       , "M_{T}^{t}(#mu)"   ,     50,     0.,    500.);
       // Mlb e 
       hists_["mEb_"        ] = store_.book1D("mEb"        , "m_{eb}"           ,     50,     0.,    500.);
       // W mass transverse estimate e
       hists_["MTWe_"       ] = store_.book1D("MTWe"       , "M_{T}^{W}(e)"     ,     60,     0.,    300.);
       // Top mass transverse estimate e
       hists_["eMTT_"       ] = store_.book1D("eMTT"       , "M_{T}^{t}(e)"     ,     50,     0.,    500.);
       // set bin labels for trigger monitoring
       triggerBinLabels(std::string("trigger"), triggerPaths_);
       // multiplicity of btagged jets (for track counting high efficiency) with pt(L2L3)>20
       hists_["jetMultBEff_"] = store_.book1D("JetMultBProb", "N_{20}(b/prob)"    ,     10,     0.,     10.);   
       // btag discriminator for track counting high efficiency for jets with pt(L2L3)>20
       hists_["jetBDiscEff_"] = store_.book1D("JetBDiscProb", "Disc_{b/prob}(jet)",     25,     0.,     2.5);   
       // multiplicity of btagged jets (for track counting high purity) with pt(L2L3)>20
 //      hists_["jetMultBPur_"] = store_.book1D("JetMultBPur", "N_{20}(b/pur)"    ,     10,     0.,     10.);   
       // btag discriminator for track counting high purity
 //      hists_["jetBDiscPur_"] = store_.book1D("JetBDiscPur", "Disc_{b/pur}(Jet)",     100,     0.,     10.);   
       // multiplicity of btagged jets (for simple secondary vertex) with pt(L2L3)>20
       hists_["jetMultBVtx_"] = store_.book1D("JetMultBVtx", "N_{20}(b/vtx)"    ,     10,     0.,     10.);   
       // btag discriminator for simple secondary vertex
       hists_["jetBDiscVtx_"] = store_.book1D("JetBDiscVtx", "Disc_{b/vtx}(Jet)",     35,    -1.,      6.);   
       // selected events
       hists_["eventLogger_"] = store_.book2D("EventLogger", "Logged Events"    ,      3,     0.,      3.,   4,   0.,   4.);
       // 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_"]->setAxisTitle("logged evts"         , 2);
 
       // deltaR min between hlt iso lepton and reco iso lepton wrt eta
       hists_["leptDeltaREta_"] = store_.book2D("DeltaRMinEtaLepton", "#Delta R_{min}(leptons) wrt #eta", 30, -3, 3, 10, 0., 0.1);   
       // deltaR min between hlt jets and reco jets wrt eta
       hists_["jetDeltaREta_"] = store_.book2D("DeltaRMinEtaJet", "#Delta R_{min}(jets) wrt #eta", 30, -3, 3, 10, 0., 0.1);   
       // resolution in pT for matched isolated leptons
       hists_["leptResolution_"] = store_.book1D("ResIsoLeptons", "#Delta p_{T}/p_{T}(matched leptons)", 20, 0., 0.1);   
       // resolution in pT for matched jets
       hists_["jetResolution_"] = store_.book1D("ResIsoJets", "#Delta p_{T}/p_{T}(matched jets)", 20, 0., 0.1);   
       // matching monitoring
       hists_["matchingMon_"] = store_.book1D("MatchingMon", "Mon(matching)", 5, 0., 5.);   
       // set axes titles for matching monitoring
       hists_["matchingMon_"]->setBinLabel( 1 , "iso lepton" );
       hists_["matchingMon_"]->setBinLabel( 2 , "1st jet" );
       hists_["matchingMon_"]->setBinLabel( 3 , "2nd jet" );
       hists_["matchingMon_"]->setBinLabel( 4 , "3rd jet" );
       hists_["matchingMon_"]->setBinLabel( 5 , "all " );
 
       return;
     }

bool HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::booked ( const std::string histName ) const

inlineprivate

check if histogram was booked

Definition at line 74 of file TopSingleLeptonHLTOfflineDQM.h.

References hists_.

Referenced by fill().

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

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::hists_

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

histogram container

Definition: TopSingleLeptonHLTOfflineDQM.h:156

void MonitorSingleLepton::fill	(	const edm::Event &	event,
		const edm::EventSetup &	setup,
		const HLTConfigProvider &	hltConfig,
		const std::vector< std::string >	triggerPaths
	)

fill monitor histograms with electronId and jetCorrections

Definition at line 240 of file TopSingleLeptonHLTOfflineDQM.cc.

References btagEff_, btagEffWP_, btagVtx_, btagVtxWP_, JetCorrector::correction(), mvaPFMET_cff::corrector, deltaR(), HLTOfflineDQMTopSingleLepton::DRMIN, alignCSCRings::e, eidPattern_, elecIso_, elecs_, elecSelect_, electronId_, electronIds_, electronRefs_, eta, edm::EventSetup::find(), edm::Event::getByToken(), JetCorrector::getJetCorrector(), hasRawTriggerSummary, hists_, i, customizeTrackingMonitorSeedNumber::idx, includeBTag_, trackerHitRTTI::isMatched(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), j, metsig::jet, jetCorrector_, jetIDLabel_, jetIDSelect_, fwrapper::jets, jets_, jetSelect_, relval_2017::k, WDecay::kElec, WDecay::kMuon, logged_, lowerEdge_, CalculateHLT::masslb(), CalculateHLT::massTopQuark(), CalculateHLT::massWBoson(), HLTOfflineDQMTopSingleLepton::MAXJETS, objects.METAnalyzer::met, mets_, HLTConfigProvider::moduleLabels(), HLTConfigProvider::moduleType(), RPCpg::mu, VarParsing::mult, metsig::muon, muonIds_, muonIso_, muonRefs_, patZpeak::muons, muons_, muonSelect_, mergeVDriftHistosByStation::name, pfjetIds_, pfjetRefs_, processName_, EnergyCorrector::pt, MetAnalyzer::pv(), pvs_, pvSelect_, reco::Jet::scaleEnergy(), EgammaValidation_Wenu_cff::sel, findQualityFiles::size, CalculateHLT::tmassTopQuark(), CalculateHLT::tmassWBoson(), edm::TriggerNames::triggerIndex(), edm::TriggerNames::triggerNames(), triggerPaths_, triggerSummaryTokenAOD, triggerSummaryTokenRAW, triggerTable_, upperEdge_, and HLTOfflineDQMTopSingleLepton::WMASS.

Referenced by fill().

     {
       // fetch trigger event if configured such 
       edm::Handle<edm::TriggerResults> triggerTable;
       if(!triggerTable_.isUninitialized()) {
         if( !event.getByToken(triggerTable_, triggerTable) ) return;
       }
 
       /*
          ------------------------------------------------------------
 
          Primary Vertex Monitoring
 
          ------------------------------------------------------------
          */
       // fill monitoring plots for primary verices
       edm::Handle<edm::View<reco::Vertex> > pvs;
       if( !event.getByToken(pvs_, pvs) ) {
         edm::LogWarning( "TopSingleLeptonHLTOfflineDQM" ) 
           << "Vertex collection not found \n";
         return;
       }
       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::GsfElectron> > elecs;
       if( !event.getByToken(elecs_, elecs) ) {
         edm::LogWarning( "TopSingleLeptonHLTOfflineDQM" ) 
           << "Electron collection not found \n";
         return;
       }
 
 
       // check availability of electron id
       edm::Handle<edm::ValueMap<float> > electronId; 
       if(!electronId_.isUninitialized()) {
         if( !event.getByToken(electronId_, electronId) ) return;
       }
 
       // loop electron collection
       unsigned int eMultIso=0;
       std::vector<const reco::GsfElectron*> isoElecs;
       reco::GsfElectron e;
       for(edm::View<reco::GsfElectron>::const_iterator elec=elecs->begin(); elec!=elecs->end(); ++elec){
         unsigned int idx = elec-elecs->begin();
         // restrict to electrons with good electronId
         if( electronId_.isUninitialized() ? true : ((int)(*electronId)[elecs->refAt(idx)] & eidPattern_) ){
           if(!elecSelect_ || (*elecSelect_)(*elec)){
             if(!elecIso_ || (*elecIso_)(*elec)){ if(eMultIso == 0) e = *elec;  isoElecs.push_back(&(*elec)); ++eMultIso;}
           }
         }
       }
 
       /* 
          ------------------------------------------------------------
 
          Muon Monitoring
 
          ------------------------------------------------------------
          */
 
       // fill monitoring plots for muons
       unsigned int mMultIso=0;
       edm::Handle<edm::View<reco::Muon> > muons;
       std::vector<const reco::Muon*> isoMuons;
       if( !event.getByToken(muons_, muons) ) {
         edm::LogWarning( "TopSingleLeptonHLTOfflineDQM" ) 
           << "Muon collection not found \n";
         return;
       }
       reco::Muon mu;
       for(edm::View<reco::Muon>::const_iterator muon=muons->begin(); muon!=muons->end(); ++muon){
         // restrict to globalMuons
         if( muon->isGlobalMuon() ){ 
           // apply preselection
           if(!muonSelect_ || (*muonSelect_)(*muon)){
             if(!muonIso_ || (*muonIso_)(*muon)) {if(mMultIso == 0) mu = *muon; isoMuons.push_back(&(*muon)); ++mMultIso;}
           }
         }
       }
 
       /* 
          ------------------------------------------------------------
 
          Jet Monitoring
 
          ------------------------------------------------------------
          */
 
       // check availability of the btaggers
       edm::Handle<reco::JetTagCollection> btagEff, btagPur, btagVtx;
       if( includeBTag_ ){ 
         if( !event.getByToken(btagEff_, btagEff) ) return;
 //        if( !event.getByToken(btagPur_, btagPur) ) return;
         if( !event.getByToken(btagVtx_, btagVtx) ) return;
       }
       // load jet corrector if configured such
       const JetCorrector* corrector=0;
       if(!jetCorrector_.empty()){
         // check whether a jet corrector is in the event setup or not
         if(setup.find( edm::eventsetup::EventSetupRecordKey::makeKey<JetCorrectionsRecord>() )){
           corrector = JetCorrector::getJetCorrector(jetCorrector_, setup);
         }
         else{ 
           edm::LogVerbatim( "TopSingleLeptonHLTOfflineDQM" ) 
             << "\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;
       unsigned int mult=0, multBEff=0, multBVtx=0; //multBPur=0;
 
       edm::Handle<edm::View<reco::Jet> > jets; 
       if( !event.getByToken(jets_, jets) ) {
         edm::LogWarning( "TopSingleLeptonHLTOfflineDQM" ) 
           << "Jet collection not found \n";
         return;
       }
 
       edm::Handle<reco::JetIDValueMap> jetID; 
       if(jetIDSelect_){ 
         if( !event.getByToken(jetIDLabel_, jetID) ) return;
       }
       reco::Jet bJetCand;   
       for(edm::View<reco::Jet>::const_iterator jet=jets->begin(); jet!=jets->end(); ++jet){
         // check jetID for calo jets
         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;
         }
 
         // prepare jet to fill monitor histograms
         reco::Jet monitorJet = *jet; monitorJet.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
         correctedJets.push_back(monitorJet);
         ++mult; // determine jet multiplicity
         if( includeBTag_ ){
           // fill b-discriminators
           edm::RefToBase<reco::Jet> jetRef = jets->refAt(idx);  
           // for the btagEff collection
           double btagEffDisc = (*btagEff)[jetRef];
           fill("jetBDiscEff_", btagEffDisc);
           if( (*btagEff)[jetRef]>btagEffWP_ ) ++multBEff; 
           // for the btagPur collection
 //          double btagPurDisc = (*btagPur)[jetRef];
 //          fill("jetBDiscPur_", btagPurDisc);
 //          if( (*btagPur)[jetRef]>btagPurWP_ ) {if(multBPur == 0) bJetCand = *jet; ++multBPur;} 
 
           // for the btagVtx collection
           double btagVtxDisc = (*btagVtx)[jetRef];
           fill("jetBDiscVtx_", btagVtxDisc);
           if( (*btagVtx)[jetRef]>btagVtxWP_ ) {if(multBVtx == 0) bJetCand = *jet; ++multBVtx;}
 
         }   
       }
 
       fill("jetMult_"    , mult    );
       fill("jetMultBEff_", multBEff);
 //      fill("jetMultBPur_", multBPur);
       fill("jetMultBVtx_", multBVtx);
         
       /* 
          ------------------------------------------------------------
 
          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()){
           mET = *(met->begin());
         }
       }
 
       /* 
          ------------------------------------------------------------
 
          Event Monitoring
 
          ------------------------------------------------------------
          */
 
       // fill W boson and top mass estimates
       CalculateHLT 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
           // We're doing a static_cast here to denote the explicity of the cast
           double runID = static_cast<double>(event.eventAuxiliary().run());
           double luminosityBlockID = static_cast<double>(event.eventAuxiliary().luminosityBlock());
           double eventID = static_cast<double>(event.eventAuxiliary().event());
           fill("eventLogger_", 0.5, logged_+0.5, runID); 
           fill("eventLogger_", 1.5, logged_+0.5, luminosityBlockID); 
           fill("eventLogger_", 2.5, logged_+0.5, eventID); 
           ++logged_;
         }
       }
       if(multBVtx != 0 && mMultIso == 1 ){
         double mtW = eventKinematics.tmassWBoson(&mu,mET,bJetCand); if (mtW == mtW) fill("MTWm_",mtW);
         double Mlb = eventKinematics.masslb(&mu,mET,bJetCand); if (Mlb == Mlb) fill("mMub_", Mlb);
         double MTT = eventKinematics.tmassTopQuark(&mu,mET,bJetCand); if (MTT == MTT) fill("mMTT_", MTT);
       }
 
       if(multBVtx != 0 && eMultIso == 1 ){
         double mtW = eventKinematics.tmassWBoson(&mu,mET,bJetCand); if (mtW == mtW)fill("MTWe_",mtW);
         double Mlb = eventKinematics.masslb(&mu,mET,bJetCand); if (Mlb == Mlb) fill("mEb_", Mlb);
         double MTT = eventKinematics.tmassTopQuark(&mu,mET,bJetCand); if (MTT == MTT) fill("eMTT_", MTT);
       }
 
 
       /* 
          ------------------------------------------------------------
 
          HLT Objects Monitoring
 
          ------------------------------------------------------------
          */
 
 //      edm::Handle<trigger::TriggerEventWithRefs> triggerEventWithRefsHandle;
 //      if(event.getByToken(triggerEventWithRefsTag_,triggerEventWithRefsHandle)) {
 
         edm::Handle<trigger::TriggerEventWithRefs> rawTriggerEvent;
         event.getByToken(triggerSummaryTokenRAW,rawTriggerEvent);
      
         edm::Handle<trigger::TriggerEventWithRefs> aodTriggerEvent;
         event.getByToken(triggerSummaryTokenAOD,aodTriggerEvent);
       
         hasRawTriggerSummary=true;
         if(!rawTriggerEvent.isValid()){ 
             hasRawTriggerSummary=false;
 //          edm::LogWarning( "TopSingleLeptonHLTOfflineDQM" ) 
 //          << "No RAW trigger summary found! Returning... \n";
        }
 
        if(!hasRawTriggerSummary && !aodTriggerEvent.isValid()){
 //      edm::LogWarning( "TopSingleLeptonHLTOfflineDQM" ) 
 //              << "No AOD trigger summary found! Returning... \n";
                 return;   
         }       
         
         const edm::TriggerNames& triggerNames = event.triggerNames(*triggerTable);
         // loop over trigger paths 
         for(unsigned int i=0; i<triggerNames.triggerNames().size(); ++i){
           // consider only path from triggerPaths
           string name = triggerNames.triggerNames()[i];
           bool isInteresting = false;
           for (unsigned int j=0; j<triggerPaths.size(); j++) {
             if (TString(name.c_str()).Contains(TString(triggerPaths[j]), TString::kIgnoreCase)) isInteresting = true; 
           }
           if (!isInteresting) continue;
           // dump infos on the considered trigger path 
           const unsigned int triggerIndex = triggerNames.triggerIndex(name);
           // get modules for the considered trigger path
           const vector<string>& moduleLabels(hltConfig.moduleLabels(triggerIndex));
           const unsigned int moduleIndex(triggerTable->index(triggerIndex));
           // Results from TriggerEventWithRefs product
           electronIds_.clear(); electronRefs_.clear();
           muonIds_.clear();     muonRefs_.clear();
           pfjetIds_.clear();    pfjetRefs_.clear();
           // look only for modules actually run in this path
           unsigned int kElec=0;
           unsigned int kMuon=0;
           unsigned int kJet=0;
           for (unsigned int k=0; k<=moduleIndex; ++k) {
             const string& moduleLabel(moduleLabels[k]);
             const string  moduleType(hltConfig.moduleType(moduleLabel));
             // check whether the module is packed up in TriggerEventWithRef product
 //            const unsigned int filterIndex(triggerEventWithRefsHandle->filterIndex(edm::InputTag(moduleLabel,"",processName_)));
             if(hasRawTriggerSummary){
                 const unsigned int filterIndex(rawTriggerEvent->filterIndex(edm::InputTag(moduleLabel,"",processName_)));
                 if (filterIndex<rawTriggerEvent->size()) {
                   rawTriggerEvent->getObjects(filterIndex,electronIds_,electronRefs_);
                   const unsigned int nElectrons(electronIds_.size());
                   if (nElectrons>0) kElec = k;
 
                   rawTriggerEvent->getObjects(filterIndex,muonIds_,muonRefs_);
                   const unsigned int nMuons(muonIds_.size());
                   if (nMuons>0) kMuon = k;
 
                   rawTriggerEvent->getObjects(filterIndex,pfjetIds_,pfjetRefs_);
                   const unsigned int nPFJets(pfjetIds_.size());
                   if (nPFJets>0) kJet = k;
                 }
             }
             
             else{
                 const unsigned int filterIndex(aodTriggerEvent->filterIndex(edm::InputTag(moduleLabel,"",processName_)));
                 if (filterIndex<aodTriggerEvent->size()) {
                   aodTriggerEvent->getObjects(filterIndex,electronIds_,electronRefs_);
                   const unsigned int nElectrons(electronIds_.size());
                   if (nElectrons>0) kElec = k;
 
                   aodTriggerEvent->getObjects(filterIndex,muonIds_,muonRefs_);
                   const unsigned int nMuons(muonIds_.size());
                   if (nMuons>0) kMuon = k;
 
                   aodTriggerEvent->getObjects(filterIndex,pfjetIds_,pfjetRefs_);
                   const unsigned int nPFJets(pfjetIds_.size());
                   if (nPFJets>0) kJet = k;
                 }
             }
                     
           }
           bool isMatched = true;
           bool lMatched = false;
           bool j1Matched = false;
           bool j2Matched = false;
           bool j3Matched = false;
 
 // access to hlt elecs
           double eDeltaRMin = 500.;
           unsigned int eIndMatched = 500;
           electronIds_.clear(); electronRefs_.clear();
           if (kElec > 0) {
             const string& moduleLabelElec(moduleLabels[kElec]);
             const string  moduleTypeElec(hltConfig.moduleType(moduleLabelElec));
 //            const unsigned int filterIndexElec(triggerEventWithRefsHandle->filterIndex(edm::InputTag(moduleLabelElec,"",processName_)));
 //            triggerEventWithRefsHandle->getObjects(filterIndexElec,electronIds_,electronRefs_);
 
             if(hasRawTriggerSummary){
                 const unsigned int filterIndexElec(rawTriggerEvent->filterIndex(edm::InputTag(moduleLabelElec,"",processName_)));
                 rawTriggerEvent->getObjects(filterIndexElec,electronIds_,electronRefs_);
             }
             else{
                 const unsigned int filterIndexElec(aodTriggerEvent->filterIndex(edm::InputTag(moduleLabelElec,"",processName_)));
                 aodTriggerEvent->getObjects(filterIndexElec,electronIds_,electronRefs_);
             }               
             
             for (unsigned int inde = 0; inde < isoElecs.size(); inde++) {
               double deltar = deltaR(*electronRefs_[0],*isoElecs[inde]); 
               if (deltar < eDeltaRMin) {
                 eDeltaRMin = deltar;
                 eIndMatched = inde;
               }
             }
             if (eDeltaRMin < DRMIN) lMatched = true;
           }
 
 // access to hlt muons
           muonIds_.clear(); muonRefs_.clear();
           double mDeltaRMin = 500.;
           unsigned int mIndMatched = 500;
           if (kMuon > 0) {
             const string& moduleLabelMuon(moduleLabels[kMuon]);
             const string  moduleTypeMuon(hltConfig.moduleType(moduleLabelMuon));
 //            const unsigned int filterIndexMuon(triggerEventWithRefsHandle->filterIndex(edm::InputTag(moduleLabelMuon,"",processName_)));
 //            triggerEventWithRefsHandle->getObjects(filterIndexMuon,muonIds_,muonRefs_);
             
             if(hasRawTriggerSummary){
                 const unsigned int filterIndexMuon(rawTriggerEvent->filterIndex(edm::InputTag(moduleLabelMuon,"",processName_)));
                 rawTriggerEvent->getObjects(filterIndexMuon,muonIds_,muonRefs_);
             }
             else{
                 const unsigned int filterIndexMuon(aodTriggerEvent->filterIndex(edm::InputTag(moduleLabelMuon,"",processName_)));
                 aodTriggerEvent->getObjects(filterIndexMuon,muonIds_,muonRefs_);
             }       
             
             for (unsigned int indm = 0; indm < isoMuons.size(); indm++) {
               double deltar = deltaR(*muonRefs_[0],*isoMuons[indm]); 
               if (deltar < mDeltaRMin) {
                 mDeltaRMin = deltar;
                 mIndMatched = indm;
               }
             }
             if (mDeltaRMin < DRMIN) lMatched = true;
           }          
 
 // access to hlt pf jets
           const unsigned int nPFJets(pfjetIds_.size());
           pfjetIds_.clear();    pfjetRefs_.clear();
           double j1DeltaRMin = 500.;   
           double j2DeltaRMin = 500.;  
           double j3DeltaRMin = 500.; 
           unsigned int j1IndMatched = 500;
           unsigned int j2IndMatched = 500;
           unsigned int j3IndMatched = 500;
           if (kJet > 0) {
             const string& moduleLabelJet(moduleLabels[kJet]);
             const string  moduleTypeJet(hltConfig.moduleType(moduleLabelJet));
 //            const unsigned int filterIndexJet(triggerEventWithRefsHandle->filterIndex(edm::InputTag(moduleLabelJet,"",processName_)));
 //            triggerEventWithRefsHandle->getObjects(filterIndexJet,pfjetIds_,pfjetRefs_);
 
             if(hasRawTriggerSummary){
                 const unsigned int filterIndexJet(rawTriggerEvent->filterIndex(edm::InputTag(moduleLabelJet,"",processName_)));
                 rawTriggerEvent->getObjects(filterIndexJet,pfjetIds_,pfjetRefs_);
             }
             else{
                 const unsigned int filterIndexJet(aodTriggerEvent->filterIndex(edm::InputTag(moduleLabelJet,"",processName_)));
                 aodTriggerEvent->getObjects(filterIndexJet,pfjetIds_,pfjetRefs_);               
             }
                 
             for (unsigned int indj = 0; indj < correctedJets.size(); indj++) {
               double deltar1 = deltaR(*pfjetRefs_[0],correctedJets[indj]); 
               if (deltar1 < j1DeltaRMin) {j1DeltaRMin = deltar1; j1IndMatched = indj;}
               if (nPFJets > 1) { 
                 double deltar2 = deltaR(*pfjetRefs_[1],correctedJets[indj]); 
                 if (deltar2 < j2DeltaRMin) {j2DeltaRMin = deltar2; j2IndMatched = indj;}
                 if (nPFJets > 2) {
                   double deltar3 = deltaR(*pfjetRefs_[2],correctedJets[indj]); 
                   if (deltar3 < j3DeltaRMin) {j3DeltaRMin = deltar3; j3IndMatched = indj;}
                 }
               }
             }
             if (nPFJets > 0 && j1DeltaRMin < DRMIN) j1Matched = true;
             if (nPFJets > 1 && j2DeltaRMin < DRMIN) j2Matched = true;
             if (nPFJets > 2 && j3DeltaRMin < DRMIN) j3Matched = true;
           }
 
           if (eIndMatched < 500) {
             fill("leptDeltaREta_", isoElecs[eIndMatched]->eta(), eDeltaRMin);   
             if (lMatched) fill("leptResolution_", fabs(isoElecs[eIndMatched]->pt()-electronRefs_[0]->pt())/isoElecs[eIndMatched]->pt() );   
           }
           
           if (mIndMatched < 500) {
             fill("leptDeltaREta_", isoMuons[mIndMatched]->eta(), mDeltaRMin);   
             if (lMatched) fill("leptResolution_", fabs(isoMuons[mIndMatched]->pt()-muonRefs_[0]->pt())/isoMuons[mIndMatched]->pt() );   
           }
           
           if (lMatched) fill("matchingMon_", 0.5 );
           else isMatched = false;
           
           if (j1IndMatched < 500) {
             fill("jetDeltaREta_", correctedJets[j1IndMatched].eta(), j1DeltaRMin);   
             if (j1Matched) {
               fill("jetResolution_", fabs(correctedJets[j1IndMatched].pt()-pfjetRefs_[0]->pt())/correctedJets[j1IndMatched].pt() );   
               fill("matchingMon_", 1.5 );
             }
             else isMatched = false;
             if (j2IndMatched < 500) {
               fill("jetDeltaREta_", correctedJets[j2IndMatched].eta(), j2DeltaRMin);   
               if (j2Matched) {
                 fill("jetResolution_", fabs(correctedJets[j2IndMatched].pt()-pfjetRefs_[1]->pt())/correctedJets[j2IndMatched].pt() );   
                 fill("matchingMon_", 2.5 );
               }
               else isMatched = false;
               if (j3IndMatched < 500) {
                 fill("jetDeltaREta_", correctedJets[j3IndMatched].eta(), j3DeltaRMin);   
                 if (j3Matched) {
                   fill("jetResolution_", fabs(correctedJets[j3IndMatched].pt()-pfjetRefs_[2]->pt())/correctedJets[j3IndMatched].pt() );   
                   fill("matchingMon_", 3.5 );
                 }
                 else isMatched = false;
               }
             }
           }
           if (isMatched) fill("matchingMon_", 4.5 );
 
         }
 
 //      }
     }

void MonitorSingleLepton::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 177 of file TopSingleLeptonHLTOfflineDQM.h.

References acceptHLT(), fill(), customizeTrackingMonitorSeedNumber::idx, and monitorPath().

     {
       for(unsigned int idx=0; idx<labels.size(); ++idx){
         if( acceptHLT(event, triggerTable, monitorPath(labels[idx])) ){
           fill((channel+"Mon_").c_str(), idx+0.5 );
         }
       }
     }

void HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::fill	(	const std::string	histName,
		double	value
	)		const

inlineprivate

fill histogram if it had been booked before

Definition at line 76 of file TopSingleLeptonHLTOfflineDQM.h.

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

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

relativeConstraints.value

tuple value

Definition: relativeConstraints.py:54

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::hists_

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

histogram container

Definition: TopSingleLeptonHLTOfflineDQM.h:156

edm::second

U second(std::pair< T, U > const &p)

Definition: ParameterSet.cc:250

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::booked

bool booked(const std::string histName) const

check if histogram was booked

Definition: TopSingleLeptonHLTOfflineDQM.h:74

void HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::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 78 of file TopSingleLeptonHLTOfflineDQM.h.

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

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

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::hists_

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

histogram container

Definition: TopSingleLeptonHLTOfflineDQM.h:156

edm::second

U second(std::pair< T, U > const &p)

Definition: ParameterSet.cc:250

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::booked

bool booked(const std::string histName) const

check if histogram was booked

Definition: TopSingleLeptonHLTOfflineDQM.h:74

void HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::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 80 of file TopSingleLeptonHLTOfflineDQM.h.

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

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::hists_

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

histogram container

Definition: TopSingleLeptonHLTOfflineDQM.h:156

edm::second

U second(std::pair< T, U > const &p)

Definition: ParameterSet.cc:250

HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::booked

bool booked(const std::string histName) const

check if histogram was booked

Definition: TopSingleLeptonHLTOfflineDQM.h:74

std::string HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::monitorPath ( const std::string & label ) const

inlineprivate

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

Definition at line 63 of file TopSingleLeptonHLTOfflineDQM.h.

Referenced by fill(), and triggerBinLabels().

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

diffTwoXMLs.label

list label

Definition: diffTwoXMLs.py:42

std::string HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::selectionPath ( const std::string & label ) const

inlineprivate

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

Definition at line 66 of file TopSingleLeptonHLTOfflineDQM.h.

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

diffTwoXMLs.label

list label

Definition: diffTwoXMLs.py:42

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

inlineprivate

set configurable labels for trigger monitoring histograms

Definition at line 169 of file TopSingleLeptonHLTOfflineDQM.h.

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

Referenced by book().

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

Member Data Documentation

edm::EDGetTokenT< reco::JetTagCollection > HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::btagEff_

private

btag discriminator labels

Definition at line 145 of file TopSingleLeptonHLTOfflineDQM.h.

Referenced by fill(), and MonitorSingleLepton().

double HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::btagEffWP_

private

btag working points

Definition at line 149 of file TopSingleLeptonHLTOfflineDQM.h.

Referenced by fill(), and MonitorSingleLepton().

edm::EDGetTokenT< reco::JetTagCollection > HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::btagPur_

private

Definition at line 146 of file TopSingleLeptonHLTOfflineDQM.h.

double HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::btagPurWP_

private

Definition at line 149 of file TopSingleLeptonHLTOfflineDQM.h.

edm::EDGetTokenT< reco::JetTagCollection > HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::btagVtx_

private

Definition at line 147 of file TopSingleLeptonHLTOfflineDQM.h.

Referenced by fill(), and MonitorSingleLepton().

double HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::btagVtxWP_

private

Definition at line 149 of file TopSingleLeptonHLTOfflineDQM.h.

Referenced by fill(), and MonitorSingleLepton().

int HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::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 118 of file TopSingleLeptonHLTOfflineDQM.h.

Referenced by fill(), and MonitorSingleLepton().

StringCutObjectSelector<reco::GsfElectron>* HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::elecIso_

private

extra isolation criterion on electron

Definition at line 120 of file TopSingleLeptonHLTOfflineDQM.h.

Referenced by fill(), and MonitorSingleLepton().