|
|
|
#include <HLTJetMETDQMSource.h>
Definition at line 38 of file HLTJetMETDQMSource.h.
| HLTJetMETDQMSource::HLTJetMETDQMSource | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 33 of file HLTJetMETDQMSource.cc.
References custompathnamemu_, custompathnamepairs_, dbe_, dirname_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), htMin_, jetDRMatch_, jetEtaMax_, jetEtMin_, LogDebug, metMin_, muonDRMatch_, muonEtaMax_, muonEtMin_, nBins_, cmsCodeRules::cppFunctionSkipper::operator, EgammaValidation_cff::paths, plotAll_, plotEff_, plotwrtMu_, prescUsed_, processname_, ptMax_, ptMin_, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), sumEtMin_, triggerResultsLabel_, triggerSummaryLabel_, and verbose_.
: resetMe_(true), currentRun_(-99) { LogDebug("HLTJetMETDQMSource") << "constructor...."; dbe_ = Service < DQMStore > ().operator->(); if ( ! dbe_ ) { LogDebug("HLTJetMETDQMSource") << "unabel to get DQMStore service?"; } if (iConfig.getUntrackedParameter < bool > ("DQMStore", false)) { dbe_->setVerbose(0); } dirname_ = iConfig.getUntrackedParameter("dirname", std::string("HLT/JetMET/")); //dirname_ += iConfig.getParameter<std::string>("@module_label"); if (dbe_ != 0 ) { dbe_->setCurrentFolder(dirname_); } //hltconfigchanged=false; processname_ = iConfig.getParameter<std::string>("processname"); //HLTPathNamesConfigPreVal_= iConfig.getParameter<std::vector<std::string > >("HLTPaths"), verbose_ = iConfig.getUntrackedParameter< bool >("verbose", false); // plotting paramters ptMin_ = iConfig.getUntrackedParameter<double>("ptMin",0.); ptMax_ = iConfig.getUntrackedParameter<double>("ptMax",1000.); nBins_ = iConfig.getUntrackedParameter<unsigned int>("Nbins",50); plotAll_ = iConfig.getUntrackedParameter<bool>("plotAll", false); plotwrtMu_ = iConfig.getUntrackedParameter<bool>("plotwrtMu", false); plotEff_ = iConfig.getUntrackedParameter<bool>("plotEff", false); // this is the list of paths to look at. std::vector<edm::ParameterSet> paths = iConfig.getParameter<std::vector<edm::ParameterSet> >("paths"); for(std::vector<edm::ParameterSet>::iterator pathconf = paths.begin() ; pathconf != paths.end(); pathconf++) { std::pair<std::string, std::string> custompathnamepair; int prescaleused; custompathnamepair.first =pathconf->getParameter<std::string>("pathname"); custompathnamepair.second = pathconf->getParameter<std::string>("denompathname"); custompathnamepairs_.push_back(custompathnamepair); //prescaleEvt_ = parameters_.getUntrackedParameter<int>("prescaleEvt", -1); prescaleused = pathconf->getParameter<int>("prescaleused"); prescUsed_.push_back(prescaleused); // customdenompathnames_.push_back(pathconf->getParameter<std::string>("denompathname")); // custompathnames_.push_back(pathconf->getParameter<std::string>("pathname")); } custompathnamemu_ = iConfig.getUntrackedParameter("pathnameMuon", std::string("HLT_L1Mu")); triggerSummaryLabel_ = iConfig.getParameter<edm::InputTag>("triggerSummaryLabel"); triggerResultsLabel_ = iConfig.getParameter<edm::InputTag>("triggerResultsLabel"); muonEtaMax_ = iConfig.getUntrackedParameter<double>("muonEtaMax",2.5); muonEtMin_ = iConfig.getUntrackedParameter<double>("muonEtMin",0.0); muonDRMatch_ =iConfig.getUntrackedParameter<double>("muonDRMatch",0.3); jetEtaMax_ = iConfig.getUntrackedParameter<double>("jetEtaMax",5.0); jetEtMin_ = iConfig.getUntrackedParameter<double>("jetEtMin",0.0); jetDRMatch_ =iConfig.getUntrackedParameter<double>("jetDRMatch",0.3); metMin_ = iConfig.getUntrackedParameter<double>("metMin",0.0); htMin_ = iConfig.getUntrackedParameter<double>("htMin",0.0); sumEtMin_ = iConfig.getUntrackedParameter<double>("sumEtMin",0.0); }
| HLTJetMETDQMSource::~HLTJetMETDQMSource | ( | ) |
Definition at line 106 of file HLTJetMETDQMSource.cc.
{
//
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
}
| void HLTJetMETDQMSource::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 117 of file HLTJetMETDQMSource.cc.
References gather_cfg::cout, reco::deltaR(), eta(), edm::Event::getByLabel(), getTH1F(), MonitorElement::getTH1F(), hltConfig_, hltPathsAll_, hltPathsEff_, hltPathswrtMu_, i, getHLTprescales::index, edm::InputTag::instance(), j, gen::k, edm::InputTag::label(), HLTConfigProvider::moduleLabels(), N, nev_, NULL, phi, plotAll_, plotEff_, plotwrtMu_, edm::InputTag::process(), processname_, rate_Denom, rate_Eff, rate_Num, MonitorElement::setBinContent(), MonitorElement::setBinError(), mathSSE::sqrt(), trigger::TriggerJet, trigger::TriggerL1CenJet, trigger::TriggerL1ETM, trigger::TriggerL1ETT, trigger::TriggerL1ForJet, trigger::TriggerL1Mu, trigger::TriggerL1TauJet, trigger::TriggerMET, trigger::TriggerMuon, edm::Event::triggerNames(), patRefSel_triggerSelection_cff::triggerResults, triggerResultsLabel_, triggerSummaryLabel_, trigger::TriggerTET, v, and verbose_.
{
using namespace edm;
using namespace trigger;
++nev_;
if (verbose_) std::cout << " --------------- N Event ------------" << nev_ << std::endl; // LogDebug("HLTJetMETDQMSource")<< "HLTJetMETDQMSource: analyze...." ;
//currun=iEvent.id().run();
//if ( nev_==1) runno = iEvent.id().run();
/*
bool changed(true);
if (hltConfig_.init(iEvent,processname_,changed)) {
if (changed) {
hltconfigchanged=true;
//if(verbose_)
//hltConfig_.dump("Triggers");
LogWarning("HLTJetMETDQMSource") << "HLTJetMETDQMSource:analyze: The number of valid triggers has changed since beginning of job." << std::endl;
// << "Processing of events halted for summary histograms" << std::endl;
//<< "Summary histograms do not support changing configurations." << std::endl
}
}
if ( hltConfig_.size() <= 0 ) {
LogError("HLTJetMETDQMSource") << "HLT config size error" << std::endl;
return;
}
// if ( currun != runno) runchanged=true;
*/
edm::Handle<TriggerResults> triggerResults;
iEvent.getByLabel(triggerResultsLabel_,triggerResults);
if(!triggerResults.isValid()) {
edm::InputTag triggerResultsLabelFU(triggerResultsLabel_.label(),triggerResultsLabel_.instance(), "FU");
iEvent.getByLabel(triggerResultsLabelFU,triggerResults);
if(!triggerResults.isValid()) {
LogWarning("HLTJetMETDQMSource") << "TriggerResults not found, skipping event" << std::endl;
return;
}
}
const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
unsigned int npath = triggerResults->size();
edm::Handle<TriggerEvent> triggerObj;
iEvent.getByLabel(triggerSummaryLabel_,triggerObj);
if(!triggerObj.isValid()) {
edm::InputTag triggerSummaryLabelFU(triggerSummaryLabel_.label(),triggerSummaryLabel_.instance(), "FU");
iEvent.getByLabel(triggerSummaryLabelFU,triggerObj);
if(!triggerObj.isValid()) {
LogWarning("HLTJetMETDQMSource") << "TriggerEvent not found, skipping event" << std::endl;
return;
}
}
//----plot all the ditributions for JET and MET -----------------------
const trigger::TriggerObjectCollection & toc(triggerObj->getObjects());
int testint=0;
if (plotAll_) testint=1;
if (verbose_) std::cout << " plots all " << testint << std::endl;
if (plotAll_) {
if (verbose_) std::cout << " Look at basic distributions " << std::endl;
int N =-1;
int NL1=-1;
for(PathInfoCollection::iterator v = hltPathsAll_.begin(); v!= hltPathsAll_.end(); ++v ){
NL1++;
N++;
int triggertype = 0;
triggertype = v->getObjectType();
bool l1accept = false;
edm::InputTag l1testTag(v->getl1Path(),"",processname_);
const int l1index = triggerObj->filterIndex(l1testTag);
if ( l1index >= triggerObj->sizeFilters() ) {
if (verbose_) std::cout<< "no index "<< l1index << " of that name "
<< v->getl1Path() << "\t" << "\t" << l1testTag << std::endl;
continue;
}
const trigger::Vids & idtype = triggerObj->filterIds(l1index);
const trigger::Keys & l1k = triggerObj->filterKeys(l1index);
//if (verbose_) std::cout << "filterID "<< idtype << " keys " << l1k << std::endl;
//if (verbose_) std::cout << " keys " << l1k << std::endl;
l1accept = l1k.size() > 0;
//if (l1accept)
//l1accept = true ;
//
//-----------------------------------------------------
bool passed = false;
for(unsigned int i = 0; i < npath; ++i) {
if ( triggerNames.triggerName(i).find(v->getPath()) != std::string::npos
&& triggerResults->accept(i)){
passed = true;
if (verbose_) std::cout << " i " << i << " trigger name " << v->getPath() << std::endl;
break;
}
}
if(passed){
//if(!hltconfigchanged)
//fill later//rate_All->Fill(N);
if (verbose_) std::cout << " N " << N << " trigger name " << v->getPath() << std::endl;
if (!l1accept) {
edm::LogInfo("HLTJetMETDQMSource") << "l1 seed path not accepted for hlt path "
<< v->getPath() << "\t" << v->getl1Path();
}
edm::InputTag filterTag = v->getTag();
//const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getPath());
if (v->getLabel() == "dummy"){
const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getPath());
//loop over labels
for (std::vector<std::string>::const_iterator labelIter= filterLabels.begin();
labelIter!=filterLabels.end(); labelIter++) {
edm::InputTag testTag(*labelIter,"",processname_);
if (verbose_) std::cout << v->getPath() << "\t" << testTag.label() << "\t"
<< testTag.process() << std::endl;
int testindex = triggerObj->filterIndex(testTag);
if ( !(testindex >= triggerObj->sizeFilters()) ) {
if (verbose_) std::cout << "found one! " << v->getPath() << "\t"
<< testTag.label() << std::endl;
filterTag = testTag;
v->setLabel(*labelIter);
}
}
}
const int index = triggerObj->filterIndex(filterTag);
if (verbose_) std::cout << "filter index "<< index << " of that name "
<< filterTag << std::endl;
if ( index >= triggerObj->sizeFilters() ) {
if (verbose_) std::cout << "WTF no index "<< index << " of that name "
<< filterTag << std::endl;
continue; //
}
const trigger::Keys & k = triggerObj->filterKeys(index);
for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
if (verbose_) cout << " filling HLT " << v->getPath() << "\t" << toc[*ki].pt()
<< "\t" << toc[*ki].eta() << "\t" << toc[*ki].phi() << std::endl;
if ( triggertype == trigger::TriggerJet
|| triggertype == trigger::TriggerL1TauJet
|| triggertype == trigger::TriggerL1CenJet
|| triggertype == trigger::TriggerL1ForJet)
{
if (verbose_) cout << " filling HLT " << v->getPath() << "\t"
<< " *ki " << *ki << " Nki "<< N<< " "
<< toc[*ki].pt() << "\t" << toc[*ki].eta() << "\t"
<< toc[*ki].phi() << std::endl;
v->getEtHisto()->Fill(toc[*ki].pt());
v->getEtaHisto()->Fill(toc[*ki].eta());
v->getPhiHisto()->Fill(toc[*ki].phi());
v->getEtaPhiHisto()->Fill(toc[*ki].eta(), toc[*ki].phi());
//---L1 Histograms
trigger::Vids::const_iterator idtypeiter = idtype.begin();
for (trigger::Keys::const_iterator l1ki = l1k.begin();
l1ki !=l1k.end(); ++l1ki ) {
if (*idtypeiter == trigger::TriggerL1TauJet
|| *idtypeiter == trigger::TriggerL1ForJet
|| *idtypeiter == trigger::TriggerL1CenJet)
{
if (verbose_) cout << " filling L1 HLT " << v->getPath()
<< "\t" << " *l1ki " << *l1ki << " NL1 "
<< NL1<< " " << toc[*l1ki].pt() << "\t"
<< toc[*l1ki].eta() << "\t" << toc[*l1ki].phi()
<< std::endl;
if (reco::deltaR(toc[*l1ki].eta(),toc[*l1ki].phi(),toc[*ki].eta(),toc[*ki].phi()) < 0.3){
//if(!hltconfigchanged)
//fill later//rate_All_L1->Fill(NL1);
v->getL1EtHisto()->Fill(toc[*l1ki].pt());
v->getL1PhiHisto()->Fill(toc[*l1ki].phi());
v->getL1EtaHisto()->Fill(toc[*l1ki].eta());
v->getL1EtaPhiHisto()->Fill(toc[*l1ki].eta(),toc[*l1ki].phi());
}
}
++idtypeiter;
}//end L1
}
//----------------------------
else if (triggertype == trigger::TriggerMuon
||triggertype == trigger::TriggerL1Mu )
{
if (verbose_) cout << " filling HLT " << v->getPath() << "\t" << " *ki "
<< *ki << " Nki "<< N << " " << toc[*ki].pt()
<< "\t" << toc[*ki].eta() << "\t" << toc[*ki].phi() << std::endl;
v->getEtHisto()->Fill(toc[*ki].pt());
//v->getEtaHisto()->Fill(toc[*ki].eta());
//v->getPhiHisto()->Fill(toc[*ki].phi());
v->getEtaPhiHisto()->Fill(toc[*ki].eta(), toc[*ki].phi());
//---L1 Histograms
trigger::Vids::const_iterator idtypeiter = idtype.begin();
for (trigger::Keys::const_iterator l1ki = l1k.begin(); l1ki !=l1k.end(); ++l1ki ) {
if (*idtypeiter == trigger::TriggerL1Mu)
{
if (reco::deltaR(toc[*l1ki].eta(),toc[*l1ki].phi(),toc[*ki].eta(),toc[*ki].phi()) < 0.3) {
//if(!hltconfigchanged)
//fill later//rate_All_L1->Fill(NL1);
v->getL1EtHisto()->Fill(toc[*l1ki].pt());
//v->getL1PhiHisto()->Fill(toc[*l1ki].phi()());
//v->getL1EtaHisto()->Fill(toc[*l1ki].eta()());
v->getL1EtaPhiHisto()->Fill(toc[*l1ki].eta(),toc[*l1ki].phi());
}
}
++idtypeiter;
}//end L1
}
//--------------------------------------
else if ( triggertype == trigger::TriggerMET
|| triggertype == trigger::TriggerL1ETM
|| triggertype == trigger::TriggerTET){
v->getEtHisto()->Fill(toc[*ki].pt());
v->getPhiHisto()->Fill(toc[*ki].phi());
//---L1 Histograms
trigger::Vids::const_iterator idtypeiter = idtype.begin();
for (trigger::Keys::const_iterator l1ki = l1k.begin(); l1ki !=l1k.end(); ++l1ki ) {
if (*idtypeiter == trigger::TriggerL1ETM)
{
if (reco::deltaR(toc[*l1ki].eta(),toc[*l1ki].phi(),toc[*ki].eta(),toc[*ki].phi()) < 0.3) {
//NL1[npth]++;
//if(!hltconfigchanged)
//fill later//rate_All_L1->Fill(NL1);
v->getL1EtHisto()->Fill(toc[*l1ki].pt());
v->getL1PhiHisto()->Fill(toc[*l1ki].phi());
}
}
++idtypeiter;
}//end L1
}
//-------------------------------------------
else if ( triggertype == trigger::TriggerTET
|| triggertype == trigger::TriggerL1ETT){
v->getEtHisto()->Fill(toc[*ki].pt());
v->getPhiHisto()->Fill(toc[*ki].phi());
//---L1 Histograms
trigger::Vids::const_iterator idtypeiter = idtype.begin();
for (trigger::Keys::const_iterator l1ki = l1k.begin(); l1ki !=l1k.end(); ++l1ki ) {
if (*idtypeiter == TriggerL1ETT)
{
if (reco::deltaR(toc[*l1ki].eta(),toc[*l1ki].phi(),toc[*ki].eta(),toc[*ki].phi()) < 0.3) {
//NL1[npth]++;
//if(!hltconfigchanged)
//fill later//rate_All_L1->Fill(NL1);
v->getL1EtHisto()->Fill(toc[*l1ki].pt());
v->getL1PhiHisto()->Fill(toc[*l1ki].phi());
}
}
++idtypeiter;
}//end L1
}
}
}//ifpassed
}
}//----------------- end if plotAll
//----plot all the ditributions for JET and MET with respect to a muon trigger --------------
if (plotwrtMu_) {
// if (verbose_) std::cout << " Look at basic distributions with respect to muon " << std::endl;
int N =-1;
int NL1=-1;
for(PathInfoCollection::iterator v = hltPathswrtMu_.begin(); v!= hltPathswrtMu_.end(); ++v ){
NL1++;
N++;
int triggertype = 0;
triggertype = v->getObjectType();
bool l1accept = false;
edm::InputTag l1testTag(v->getl1Path(),"",processname_);
const int l1index = triggerObj->filterIndex(l1testTag);
if ( l1index >= triggerObj->sizeFilters() ) {
if (verbose_) std::cout<< "no index "<< l1index << " of that name (wrtMu)" << v->getl1Path() << "\t" << "\t" << l1testTag << std::endl;
continue;
}
//const trigger::Vids & idtype = triggerObj->filterIds(l1index);
const trigger::Keys & l1k = triggerObj->filterKeys(l1index);
//if (verbose_) std::cout << "filterID "<< idtype << " keys " << l1k << std::endl;
//if (verbose_) std::cout << " keys " << l1k << std::endl;
l1accept = l1k.size() > 0;
//if (l1accept) rate_All_L1->Fill(NL1);
//l1accept = true ;
//
//-----------------------------------------------------
bool denompassed = false;
bool passed = false;
for(unsigned int i = 0; i < npath; ++i) {
if ( triggerNames.triggerName(i).find(v->getPath()) != std::string::npos && triggerResults->accept(i)){
passed = true; break ;
}
}
for(unsigned int i = 0; i < npath; ++i) {
if ( triggerNames.triggerName(i).find(v->getDenomPath()) != std::string::npos && triggerResults->accept(i)){
denompassed = true; break ;
}
}
if(denompassed){
if (verbose_) std::cout << " N " << N << " trigger wrt mu demom name "
<< v->getDenomPath() << std::endl;
if(passed){
//if(!hltconfigchanged)
//fill later//rate_wrtMu->Fill(N);
if (verbose_) std::cout << " N " << N << " trigger name "
<< v->getPath() << std::endl;
if (!l1accept) {
edm::LogInfo("HLTJetMETDQMSource") << "l1 seed path not accepted for hlt path "
<< v->getPath() << "\t" << v->getl1Path();
}
edm::InputTag filterTag = v->getTag();
//const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getPath());
if (v->getLabel() == "dummy"){
const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getPath());
//loop over labels
for (std::vector<std::string>::const_iterator labelIter= filterLabels.begin(); labelIter!=filterLabels.end(); labelIter++) {
edm::InputTag testTag(*labelIter,"",processname_);
if (verbose_) std::cout << v->getPath() << "\t" << testTag.label() << "\t" << testTag.process() << std::endl;
int testindex = triggerObj->filterIndex(testTag);
if ( !(testindex >= triggerObj->sizeFilters()) ) {
if (verbose_) std::cout << "found one! " << v->getPath() << "\t" << testTag.label() << std::endl;
filterTag = testTag;
v->setLabel(*labelIter);
}
}
}
const int index = triggerObj->filterIndex(filterTag);
if (verbose_) std::cout << "filter index "<< index << " of that name (wrtMu)"
<< filterTag << std::endl;
if ( index >= triggerObj->sizeFilters() ) {
if (verbose_) std::cout << "WTF no index "<< index << " of that name (wrtMu)"
<< filterTag << std::endl;
continue; //
}
const trigger::Keys & k = triggerObj->filterKeys(index);
for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
if (verbose_) cout << " filling HLT (wrtMu)" << v->getPath()
<< "\t" << toc[*ki].pt() << "\t" << toc[*ki].eta()
<< "\t" << toc[*ki].phi() << std::endl;
if ( triggertype == trigger::TriggerJet
|| triggertype == trigger::TriggerL1TauJet
|| triggertype == trigger::TriggerL1CenJet
|| triggertype == trigger::TriggerL1ForJet)
{
if (verbose_) cout << " filling HLT (wrtMu)" << v->getPath()
<< "\t" << " *ki " << *ki << " Nki "<< N
<< " " << toc[*ki].pt() << "\t" << toc[*ki].eta()
<< "\t" << toc[*ki].phi() << std::endl;
v->getEtwrtMuHisto()->Fill(toc[*ki].pt());
//v->getEtawrtMuHisto()->Fill(toc[*ki].eta());
//v->getPhiwrtMuHisto()->Fill(toc[*ki].phi());
v->getEtaPhiwrtMuHisto()->Fill(toc[*ki].eta(), toc[*ki].phi());
}
//----------------------------
else if (triggertype == trigger::TriggerMuon
||triggertype == trigger::TriggerL1Mu )
{
if (verbose_) cout << " filling HLT (wrtMu)" << v->getPath() << "\t" << " *ki " << *ki << " Nki "<< N << " " << toc[*ki].pt() << "\t" << toc[*ki].eta() << "\t" << toc[*ki].phi() << std::endl;
v->getEtwrtMuHisto()->Fill(toc[*ki].pt());
//v->getEtaHisto()->Fill(toc[*ki].eta());
//v->getPhiHisto()->Fill(toc[*ki].phi());
v->getEtaPhiwrtMuHisto()->Fill(toc[*ki].eta(), toc[*ki].phi());
}
//--------------------------------------
else if ( triggertype == trigger::TriggerMET
|| triggertype == trigger::TriggerL1ETM)
{
v->getEtwrtMuHisto()->Fill(toc[*ki].pt());
v->getPhiwrtMuHisto()->Fill(toc[*ki].phi());
}
//-------------------------------------------
else if ( triggertype == trigger::TriggerTET
|| triggertype == trigger::TriggerL1ETT)
{
v->getEtwrtMuHisto()->Fill(toc[*ki].pt());
}
}
}//ifpassed
}//if denom passed
}
}//----------------- end if plowrtMu
//----------plot efficiency---------------
if (plotEff_) {
if (verbose_) std::cout << " Look at basic distributions for Eff " << std::endl;
int N =-1;
int Ndenom =-1;
int Nnum =-1;
for(PathInfoCollection::iterator v = hltPathsEff_.begin(); v!= hltPathsEff_.end(); ++v ){
Ndenom++;
Nnum++;
N++;
int triggertype = 0;
triggertype = v->getObjectType();
bool l1accept = false;
edm::InputTag l1testTag(v->getl1Path(),"",processname_);
const int l1index = triggerObj->filterIndex(l1testTag);
if ( l1index >= triggerObj->sizeFilters() ) {
if (verbose_) std::cout<< "no index "<< l1index << " of that name "
<< v->getl1Path() << "\t" << "\t" << l1testTag << std::endl;
continue;
}
//const trigger::Vids & idtype = triggerObj->filterIds(l1index);
const trigger::Keys & l1k = triggerObj->filterKeys(l1index);
//if (verbose_) std::cout << "filterID "<< idtype << " keys " << l1k << std::endl;
//if (verbose_) std::cout << " keys " << l1k << std::endl;
l1accept = l1k.size() > 0;
bool passed = false;
for(unsigned int i = 0; i < npath; ++i) {
if ( triggerNames.triggerName(i).find(v->getPath()) != std::string::npos
&& triggerResults->accept(i)){
passed = true;
if (verbose_) std::cout << " i " << i << " trigger name " << v->getPath() << std::endl;
break;
}
}
bool denompassed = false;
for(unsigned int i = 0; i < npath; ++i) {
if ( triggerNames.triggerName(i).find(v->getDenomPath()) != std::string::npos
&& triggerResults->accept(i)){
denompassed = true;
if (verbose_) std::cout << " i " << i << " trigger name " << v->getDenomPath() << std::endl;
break;
}
}
if (denompassed){
//if(!hltconfigchanged)
//fill later//rate_Denom->Fill(Ndenom);
if (verbose_) std::cout << " N " << N << " trigger name " << v->getDenomPath() << std::endl;
if (!l1accept) {
edm::LogInfo("HLTJetMETDQMSource") << "l1 seed path not accepted for hlt path "
<< v->getPath() << "\t" << v->getl1Path();
}
edm::InputTag filterTag = v->getDenomTag();
//const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getPath());
if (verbose_) std::cout << v->getDenomPath() << "\t" << v->getDenomLabel() << "\t" << std::endl;
if (v->getDenomLabel() == "denomdummy"){
const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getDenomPath());
//loop over labels
for (std::vector<std::string>::const_iterator labelIter= filterLabels.begin(); labelIter!=filterLabels.end(); labelIter++) {
edm::InputTag testTag(*labelIter,"",processname_);
if (verbose_) std::cout << v->getDenomPath() << "\t" << testTag.label() << "\t"
<< testTag.process() << std::endl;
int testindex = triggerObj->filterIndex(testTag);
if ( !(testindex >= triggerObj->sizeFilters()) ) {
if (verbose_) std::cout << "found one! " << v->getDenomPath() << "\t"
<< testTag.label() << std::endl;
filterTag = testTag;
v->setDenomLabel(*labelIter);
}
}
}
const int index = triggerObj->filterIndex(filterTag);
if (verbose_) std::cout << "filter index "<< index << " of that name "
<< filterTag << std::endl;
if ( index >= triggerObj->sizeFilters() ) {
if (verbose_) std::cout << "WTF no index "<< index << " of that name "
<< filterTag << std::endl;
continue; //
}
const trigger::Keys & k = triggerObj->filterKeys(index);
int numobj=-1;
for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
if (verbose_) cout << " filling HLT Denom " << v->getDenomPath() << "\t"
<< toc[*ki].pt() << "\t" << toc[*ki].eta() << "\t"
<< toc[*ki].phi() << std::endl;
if ( triggertype == trigger::TriggerJet
|| triggertype == trigger::TriggerL1TauJet
|| triggertype == trigger::TriggerL1CenJet
|| triggertype == trigger::TriggerL1ForJet)
{
if (verbose_) cout << " filling HLT Denom" << v->getDenomPath()
<< "\t" << " *kiDeno " << *ki << " NkiDenom "
<< Ndenom<< " " << toc[*ki].pt() << "\t"
<< toc[*ki].eta() << "\t" << toc[*ki].phi() << std::endl;
numobj=*ki;
v->getEtDenomHisto()->Fill(toc[*ki].pt());
v->getEtaDenomHisto()->Fill(toc[*ki].eta());
v->getPhiDenomHisto()->Fill(toc[*ki].phi());
}
//--------------------------------------
else if ( triggertype == trigger::TriggerMET
|| triggertype == trigger::TriggerL1ETM
|| triggertype == trigger::TriggerTET){
v->getEtDenomHisto()->Fill(toc[*ki].pt());
v->getPhiDenomHisto()->Fill(toc[*ki].phi());
}
//-------------------------------------------
else if ( triggertype == trigger::TriggerTET
|| triggertype == trigger::TriggerL1ETT){
v->getEtDenomHisto()->Fill(toc[*ki].pt());
v->getPhiDenomHisto()->Fill(toc[*ki].phi());
}
if ( numobj != -1) break;
}
}//if denom passed
if(denompassed){
if (passed){
//if(!hltconfigchanged)
//fill later//rate_Num->Fill(Nnum);
if (verbose_) std::cout << " N " << N << " trigger name " << v->getPath() << std::endl;
if (!l1accept) {
edm::LogInfo("HLTJetMETDQMSource") << "l1 seed path not accepted for hlt path "<< v->getPath() << "\t" << v->getl1Path();
}
edm::InputTag filterTag = v->getTag();
//const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getPath());
if (verbose_) std::cout << v->getPath() << "\t" << v->getLabel() << std::endl;
if (v->getLabel() == "dummy"){
const std::vector<std::string> filterLabels = hltConfig_.moduleLabels(v->getPath());
//loop over labels
for (std::vector<std::string>::const_iterator labelIter= filterLabels.begin(); labelIter!=filterLabels.end(); labelIter++) {
edm::InputTag testTag(*labelIter,"",processname_);
if (verbose_) std::cout << v->getPath() << "\t" << testTag.label() << "\t"
<< testTag.process() << std::endl;
int testindex = triggerObj->filterIndex(testTag);
if ( !(testindex >= triggerObj->sizeFilters()) ) {
if (verbose_) std::cout << "found one! " << v->getPath() << "\t"
<< testTag.label() << std::endl;
filterTag = testTag;
v->setLabel(*labelIter);
}
}
}
const int index = triggerObj->filterIndex(filterTag);
if (verbose_) std::cout << "filter index "<< index << " of that name "
<< filterTag << std::endl;
if ( index >= triggerObj->sizeFilters() ) {
if (verbose_) std::cout << "WTF no index "<< index << " of that name "
<< filterTag << std::endl;
continue; //
}
const trigger::Keys & k = triggerObj->filterKeys(index);
int numobj = -1;
for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
if (verbose_) cout << " filling HLT " << v->getPath() << "\t"
<< toc[*ki].pt() << "\t" << toc[*ki].eta()
<< "\t" << toc[*ki].phi() << std::endl;
if ( triggertype == trigger::TriggerJet
|| triggertype == trigger::TriggerL1TauJet
|| triggertype == trigger::TriggerL1CenJet
|| triggertype == trigger::TriggerL1ForJet)
{
if (verbose_) cout << " filling HLT " << v->getPath() << "\t"
<< " *ki " << *ki << " Nki "<< N<< " "
<< toc[*ki].pt() << "\t" << toc[*ki].eta()
<< "\t" << toc[*ki].phi() << std::endl;
numobj=*ki;
v->getEtNumHisto()->Fill(toc[*ki].pt());
v->getEtaNumHisto()->Fill(toc[*ki].eta());
v->getPhiNumHisto()->Fill(toc[*ki].phi());
}
//--------------------------------------
else if ( triggertype == trigger::TriggerMET
|| triggertype == trigger::TriggerL1ETM
|| triggertype == trigger::TriggerTET){
v->getEtNumHisto()->Fill(toc[*ki].pt());
v->getPhiNumHisto()->Fill(toc[*ki].phi());
}
//-------------------------------------------
else if ( triggertype == trigger::TriggerTET
|| triggertype == trigger::TriggerL1ETT){
v->getEtNumHisto()->Fill(toc[*ki].pt());
v->getPhiNumHisto()->Fill(toc[*ki].phi());
}
if ( numobj != -1) break;
}
}//end if passed
} //if denompassed
}
//----------plot efficiency---------------
if ( (nev_ % 1000) == 0 ){
if (verbose_) std::cout << " Calculating Eff.... " << nev_ << std::endl;
//if (nev_ % 1000)
//if (verbose_) std::cout << "starting the luminosity block -----" << std::endl;
TH1F *rate_denom=NULL;
TH1F *rate_num=NULL;
rate_denom = rate_Denom->getTH1F();
rate_num = rate_Num->getTH1F();
//rate_denom->Sumw2();rate_num->Sumw2();
if ((rate_denom->Integral() != 0.) && (rate_num->Integral() != 0.) ) {
if (verbose_) std::cout << " Nonzero rate summary -----" << std::endl;
for(int j=1; j <= rate_denom->GetXaxis()->GetNbins();j++ ){
double y1 = rate_num->GetBinContent(j);
double y2 = rate_denom->GetBinContent(j);
double eff = y2 > 0. ? y1/y2 : 0.;
rate_Eff->setBinContent(j, eff);
double y1err = rate_num->GetBinError(j);
double y2err = rate_denom->GetBinError(j);
double efferr = 0.0;
if (y2 && y1 > 0.) efferr = (y1/y2)* sqrt ((y1err/y1)*(y1err/y1)
+ (y2err/y2)*(y2err/y2)) ;
rate_Eff->setBinError(j, efferr);
}
}
for(PathInfoCollection::iterator v = hltPathsEff_.begin(); v!= hltPathsEff_.end(); ++v ) {
int triggertype = 0;
triggertype = v->getObjectType();
if (verbose_) std::cout << " triggertype ----- " << triggertype << std::endl;
if ( triggertype == trigger::TriggerJet
|| triggertype == trigger::TriggerL1TauJet
|| triggertype == trigger::TriggerL1CenJet
|| triggertype == trigger::TriggerL1ForJet
|| triggertype == trigger::TriggerL1Mu
|| triggertype == trigger::TriggerMuon ){
TH1F *EtNum=NULL;
TH1F *EtaNum=NULL;
TH1F *PhiNum=NULL;
TH1F *EtDenom=NULL;
TH1F *EtaDenom=NULL;
TH1F *PhiDenom=NULL;
EtNum= ( v->getEtNumHisto())->getTH1F() ;
EtaNum= ( v->getEtaNumHisto())->getTH1F() ;
PhiNum= ( v->getPhiNumHisto())->getTH1F() ;
EtDenom= ( v->getEtDenomHisto())->getTH1F() ;
EtaDenom= ( v->getEtaDenomHisto())->getTH1F() ;
PhiDenom= ( v->getPhiDenomHisto())->getTH1F() ;
if ((EtNum->Integral() != 0.) && (EtDenom->Integral() != 0.) ) {
if (verbose_) std::cout << " Nonzero Jet Et -----" << std::endl;
for(int j=1; j <= EtNum->GetXaxis()->GetNbins();j++ ){
double y1 = EtNum->GetBinContent(j);
double y2 = EtDenom->GetBinContent(j);
double eff = y2 > 0. ? y1/y2 : 0.;
v->getEtEffHisto()->setBinContent(j, eff);
double y1err = EtNum->GetBinError(j);
double y2err = EtDenom->GetBinError(j);
double efferr = 0.0;
if (y2 && y1 > 0.) efferr = (y1/y2)* sqrt ((y1err/y1)*(y1err/y1)
+ (y2err/y2)*(y2err/y2)) ;
v->getEtEffHisto()->setBinError(j, efferr);
if(verbose_) std::cout << eff << " "<< efferr << " "
<< y1 << " " << y2 << " "<< y1err << " " << y2err << std::endl;
}
}
if (EtaNum->Integral() != 0. && EtaDenom->Integral() != 0. ) {
for(int j=1; j <= EtaNum->GetXaxis()->GetNbins();j++ ){
double y1 = EtaNum->GetBinContent(j);
double y2 = EtaDenom->GetBinContent(j);
double eff = y2 > 0. ? y1/y2 : 0.;
v->getEtaEffHisto()->setBinContent(j, eff);
double y1err = EtaNum->GetBinError(j);
double y2err = EtaDenom->GetBinError(j);
double efferr = 0.0;
if (y2 && y1 > 0.) efferr = (y1/y2)* sqrt ((y1err/y1)*(y1err/y1)
+ (y2err/y2)*(y2err/y2)) ;
v->getEtaEffHisto()->setBinError(j, efferr);
//
}
}
if (PhiNum->Integral() != 0. && PhiDenom->Integral() != 0. ) {
for(int j=1; j <= PhiNum->GetXaxis()->GetNbins();j++ ){
double y1 = PhiNum->GetBinContent(j);
double y2 = PhiDenom->GetBinContent(j);
double eff = y2 > 0. ? y1/y2 : 0.;
v->getPhiEffHisto()->setBinContent(j, eff);
double y1err = PhiNum->GetBinError(j);
double y2err = PhiDenom->GetBinError(j);
double efferr = 0.0;
if (y2 && y1 > 0.) efferr = (y1/y2)* sqrt ((y1err/y1)*(y1err/y1)
+ (y2err/y2)*(y2err/y2)) ;
v->getPhiEffHisto()->setBinError(j, efferr);
//
}
}
}
if ( triggertype == trigger::TriggerMET
|| triggertype == trigger::TriggerTET ){
TH1F *EtNum=NULL;
TH1F *PhiNum=NULL;
TH1F *EtDenom=NULL;
TH1F *PhiDenom=NULL;
EtNum= ( v->getEtNumHisto())->getTH1F() ;
PhiNum= ( v->getPhiNumHisto())->getTH1F() ;
EtDenom= ( v->getEtDenomHisto())->getTH1F() ;
PhiDenom= ( v->getPhiDenomHisto())->getTH1F() ;
if (EtNum->Integral() != 0. && EtDenom->Integral() != 0. ) {
if (verbose_) std::cout << " Nonzero Met Et -----" << std::endl;
for(int j=1; j <= EtNum->GetXaxis()->GetNbins();j++ ){
double y1 = EtNum->GetBinContent(j);
double y2 = EtDenom->GetBinContent(j);
double eff = y2 > 0. ? y1/y2 : 0.;
v->getEtEffHisto()->setBinContent(j, eff);
double y1err = EtNum->GetBinError(j);
double y2err = EtDenom->GetBinError(j);
double efferr = 0.0;
if (y2 && y1 > 0.) efferr = (y1/y2)* sqrt ((y1err/y1)*(y1err/y1)
+ (y2err/y2)*(y2err/y2)) ;
v->getEtEffHisto()->setBinError(j, efferr);
}
}
if (PhiNum->Integral() != 0. && PhiDenom->Integral() != 0. ) {
for(int j=1; j <= PhiNum->GetXaxis()->GetNbins();j++ ){
double y1 = PhiNum->GetBinContent(j);
double y2 = PhiDenom->GetBinContent(j);
double eff = y2 > 0. ? y1/y2 : 0.;
v->getPhiEffHisto()->setBinContent(j, eff);
double y1err = PhiNum->GetBinError(j);
double y2err = PhiDenom->GetBinError(j);
double efferr = 0.0;
if (y2 && y1 > 0.) efferr = (y1/y2)* sqrt ((y1err/y1)*(y1err/y1)
+ (y2err/y2)*(y2err/y2)) ;
v->getPhiEffHisto()->setBinError(j, efferr);
}
}
}// met triggers
}
}
}
}
| void HLTJetMETDQMSource::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 892 of file HLTJetMETDQMSource.cc.
References dirname_, nev_, cmsCodeRules::cppFunctionSkipper::operator, DQMStore::rmdir(), and DQMStore::setCurrentFolder().
{
nev_ = 0;
DQMStore *dbe = 0;
dbe = Service<DQMStore>().operator->();
if (dbe) {
dbe->setCurrentFolder(dirname_);
dbe->rmdir(dirname_);
}
if (dbe) {
dbe->setCurrentFolder(dirname_);
}
/*
hltconfigchanged=false;
//if (!hltConfig_.init(processname_)) {
if (!hltConfig_.init(run, c, processname_, hltconfigchanged)) {
processname_ = "FU";
// if (!hltConfig_.init(processname_)){
if (!hltConfig_.init(run, c, processname_, hltconfigchanged))
LogError("HLTJetMETDQMSource") << "HLTConfigProvider failed to initialize." << std::endl;
}
// hltConfig_.dump("Triggers");
}
const unsigned int n(hltConfig_.size());
if (verbose_) std::cout << " hltConfig_.size() " << n << std::endl;
//-------------if plot all jetmet trigger pt,eta, phi-------------------
if (plotAll_){
if (verbose_) std::cout << " booking histos All " << std::endl;
std::string foldernm = "/All/";
if (dbe) {
dbe->setCurrentFolder(dirname_ + foldernm);
}
for (unsigned int j=0; j!=n; ++j) {
std::string pathname = hltConfig_.triggerName(j);
std::string l1pathname = "dummy";
//for (unsigned int i=0; i!=n; ++i) {
// std::cout << hltConfig_.triggerName(i) << std::endl;
//std::string denompathname = hltConfig_.triggerName(i);
std::string denompathname = "";
unsigned int usedPresscale = 1;
unsigned int objectType = 0;
//int denomobjectType = 0;
//parse pathname to guess object type
if ( pathname.find("Mu") && (pathname=="HLT_L1MuOpen"))
objectType = trigger::TriggerMuon;
else if (pathname.find("MET") != std::string::npos) objectType = trigger::TriggerMET;
else if (pathname.find("L1MET") != std::string::npos) objectType = trigger::TriggerL1ETM;
else if (pathname.find("SumET") != std::string::npos) objectType = trigger::TriggerTET;
else if (pathname.find("Jet") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet30") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet50") != std::string::npos) objectType = trigger::TriggerJet;
//if ((pathname.find("HLT_Mu3") != std::string::npos) || (pathname.find("HLT_L2Mu9") != std::string::npos) ) objectType = trigger::TriggerMuon;
else continue;
//std::cout << "objecttye " << objectType << std::endl;
// find L1 condition for numpath with numpath objecttype
// find PSet for L1 global seed for numpath,
// list module labels for numpath
std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
numpathmodule!= numpathmodules.end();
++numpathmodule ) {
// std::cout << pathname << "\t" << *numpathmodule << "\t"
//<< hltConfig_.moduleType(*numpathmodule) << std::endl;
if (hltConfig_.moduleType(*numpathmodule) == "HLTLevel1GTSeed") {
edm::ParameterSet l1GTPSet = hltConfig_.modulePSet(*numpathmodule);
//std::cout << l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression") << std::endl;
// l1pathname = l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression");
l1pathname = *numpathmodule;
break;
}
}
std::string filtername("dummy");
std::string Denomfiltername("denomdummy");
float ptMin = 0.0;
float ptMax = 300.0;
if (objectType != 0 )
hltPathsAll_.push_back(PathInfo(usedPresscale, denompathname,
pathname, l1pathname, filtername,
Denomfiltername, processname_,
objectType, ptMin, ptMax));
}
std::string histonm="JetMET_rate_All";
std::string histonmL1="JetMET_rate_All_L1";
std::string histot="JetMET Rate Summary";
std::string histoL1t="JetMET L1 Rate Summary";
//fill later//rate_All = dbe->book1D(histonm.c_str(),histot.c_str(),
hltPathsAll_.size()+1,-0.5,hltPathsAll_.size()+1-0.5);
//fill later//rate_All_L1 = dbe->book1D(histonmL1.c_str(),histoL1t.c_str(),
hltPathsAll_.size()+1,-0.5,hltPathsAll_.size()+1-0.5);
//rate_All->setBinLabel(hltPathsAll_.size()+1,"Rate",1);
unsigned int nname=0;
unsigned int nnameL1=0;
for(PathInfoCollection::iterator v = hltPathsAll_.begin(); v!= hltPathsAll_.end(); ++v ){
std::string labelnm("dummy");
labelnm = v->getPath();
rate_All->setBinLabel(nname+1,labelnm);
nname++;
std::string labelnml1("dummyl1");
labelnml1 = v->getl1Path();
rate_All_L1->setBinLabel(nnameL1+1,labelnml1);
nnameL1++;
}
// now set up all of the histos for each path
for(PathInfoCollection::iterator v = hltPathsAll_.begin(); v!= hltPathsAll_.end(); ++v ) {
MonitorElement *N = 0;
MonitorElement *Et = 0;
MonitorElement *EtaPhi = 0;
MonitorElement *Eta = 0;
MonitorElement *Phi = 0;
MonitorElement *NL1 = 0;
MonitorElement *l1Et = 0;
MonitorElement *l1EtaPhi = 0;
MonitorElement *l1Eta = 0;
MonitorElement *l1Phi = 0;
std::string labelname("dummy");
labelname = v->getPath();
std::string histoname(labelname+"");
std::string title(labelname+"");
double histEtaMax = 2.5;
if (v->getObjectType() == trigger::TriggerMuon
|| v->getObjectType() == trigger::TriggerL1Mu) {
histEtaMax = muonEtaMax_;
nBins_ = 25 ;
}
else if (v->getObjectType() == trigger::TriggerJet
|| v->getObjectType() == trigger::TriggerL1CenJet
|| v->getObjectType() == trigger::TriggerL1ForJet ) {
histEtaMax = jetEtaMax_;
nBins_ = 60 ;
}
else if (v->getObjectType() == trigger::TriggerMET
|| v->getObjectType() == trigger::TriggerL1ETM ) {
histEtaMax = 5.0;
nBins_ = 60 ;
}
TString pathfolder = dirname_ + foldernm + v->getPath();
dbe_->setCurrentFolder(pathfolder.Data());
if (verbose_) std::cout << "Booking Histos in Directory " << pathfolder.Data() << std::endl;
int nBins2D = 10;
histoname = labelname+"_Et";
title = labelname+" E_t";
Et = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
histoname = labelname+"_l1Et";
title = labelname+" L1 E_t";
l1Et = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
if (labelname.find("Jet") != std::string::npos
|| labelname.find("Mu") != std::string::npos) {
histoname = labelname+"_EtaPhi";
title = labelname+" #eta vs #phi";
EtaPhi = dbe->book2D(histoname.c_str(),
title.c_str(),
nBins2D,-histEtaMax,histEtaMax,
nBins2D,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_l1EtaPhi";
title = labelname+"L1 #eta vs L1 #phi";
l1EtaPhi = dbe->book2D(histoname.c_str(),
title.c_str(),
nBins2D,-histEtaMax,histEtaMax,
nBins2D,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Phi";
title = labelname+" #phi";
Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_l1Phi";
title = labelname+"L1 #phi";
l1Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Eta";
title = labelname+" #eta";
Eta = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax
);
histoname = labelname+"_l1Eta";
title = labelname+"L1 #eta";
l1Eta = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
}
else if( (labelname.find("MET") != std::string::npos)
|| (labelname.find("SumET") != std::string::npos) ){
histoname = labelname+"_phi";
title = labelname+" #phi";
Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_l1Phi";
title = labelname+"L1 #phi";
l1Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
v->setHistos( N, Et, EtaPhi, Eta, Phi, NL1, l1Et, l1EtaPhi,l1Eta, l1Phi);
}
if (verbose_) std::cout << "Done booking histos All " << std::endl;
}
//----------------plot all jetmet trigger wrt some muon trigger-----
if (plotwrtMu_){
if (verbose_) std::cout << " booking histos wrt Muon " << std::endl;
std::string foldernm = "/wrtMuon/";
if (dbe) {
dbe->setCurrentFolder(dirname_ + foldernm);
}
for (unsigned int j=0; j!=n; ++j) {
std::string pathname = hltConfig_.triggerName(j);
std::string l1pathname = "dummy";
//for (unsigned int i=0; i!=n; ++i)
//{
if (verbose_) std::cout << hltConfig_.triggerName(j) << std::endl;
std::string denompathname = custompathnamemu_ ;
int objectType = 0;
int usedPresscale = 1;
//int denomobjectType = 0;
//parse pathname to guess object type
if ( pathname.find("Mu") && (pathname=="HLT_L1MuOpen"))
objectType = trigger::TriggerMuon;
else if (pathname.find("MET") != std::string::npos) objectType = trigger::TriggerMET;
else if (pathname.find("L1MET") != std::string::npos) objectType = trigger::TriggerL1ETM;
else if (pathname.find("SumET") != std::string::npos) objectType = trigger::TriggerTET;
else if (pathname.find("Jet") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet30") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet50") != std::string::npos) objectType = trigger::TriggerJet;
//if ((pathname.find("HLT_Mu3") != std::string::npos) || (pathname.find("HLT_L2Mu9") != std::string::npos) ) objectType = trigger::TriggerMuon;
else continue;
// find L1 condition for numpath with numpath objecttype
// find PSet for L1 global seed for numpath,
// list module labels for numpath
std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
// std::cout << pathname << "\t" << *numpathmodule << "\t" << hltConfig_.moduleType(*numpathmodule) << std::endl;
if (hltConfig_.moduleType(*numpathmodule) == "HLTLevel1GTSeed") {
edm::ParameterSet l1GTPSet = hltConfig_.modulePSet(*numpathmodule);
//std::cout << l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression") << std::endl;
// l1pathname = l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression");
l1pathname = *numpathmodule;
break;
}
}
std::string filtername("dummy");
std::string Denomfiltername("denomdummy");
float ptMin = 0.0;
float ptMax = 300.0;
if ( objectType != 0){
// if (verbose_) std::cout << " wrt muon PathInfo(denompathname, pathname, l1pathname, filtername, denomfiltername, processname_, objectType, ptMin, ptMax " << denompathname << " "<< pathname << " "<< l1pathname << " " << filtername << " " << Denomfiltername << " " << processname_ << " " << objectType << " " << ptMin << " " << ptMax<< std::endl;
hltPathswrtMu_.push_back(PathInfo(usedPresscale, denompathname, pathname, l1pathname,
filtername, Denomfiltername, processname_,
objectType, ptMin, ptMax));
}
}
std::string histonm="JetMET_rate_wrt_" + custompathnamemu_ + "_Summary";
std::string histt="JetMET Rate wrt " + custompathnamemu_ + "Summary";
rate_wrtMu = dbe->book1D(histonm.c_str(),histt.c_str(),
hltPathswrtMu_.size()+1,-0.5,hltPathswrtMu_.size()+1.0-0.5);
int nname=0;
for(PathInfoCollection::iterator v = hltPathswrtMu_.begin(); v!= hltPathswrtMu_.end(); ++v ){
std::string labelnm("dummy");
labelnm = v->getPath();
rate_wrtMu->setBinLabel(nname+1,labelnm);
nname++;
}
// now set up all of the histos for each path
for(PathInfoCollection::iterator v = hltPathswrtMu_.begin(); v!= hltPathswrtMu_.end(); ++v ) {
MonitorElement *NwrtMu = 0;
MonitorElement *EtwrtMu = 0;
MonitorElement *EtaPhiwrtMu = 0;
MonitorElement *PhiwrtMu = 0;
std::string labelname("dummy");
labelname = v->getPath() + "_wrt_" + v->getDenomPath();
std::string histoname(labelname+"");
std::string title(labelname+"");
double histEtaMax = 2.5;
if (v->getObjectType() == trigger::TriggerMuon
|| v->getObjectType() == trigger::TriggerL1Mu) {
histEtaMax = muonEtaMax_;nBins_ = 20 ;
}
else if (v->getObjectType() == trigger::TriggerJet
|| v->getObjectType() == trigger::TriggerL1CenJet
|| v->getObjectType() == trigger::TriggerL1ForJet ){
histEtaMax = jetEtaMax_; nBins_ = 60 ;
}
else if (v->getObjectType() == trigger::TriggerMET
|| v->getObjectType() == trigger::TriggerL1ETM ) {
histEtaMax = 5.0; nBins_ = 60 ;
}
TString pathfolder = dirname_ + foldernm + v->getPath();
dbe_->setCurrentFolder(pathfolder.Data());
if (verbose_) std::cout << "Booking Histos in Directory " << pathfolder.Data() << std::endl;
int nBins2D = 10;
//pathfolder = dirname_ + TString("/wrtMuon/") + v->getPath();
histoname = labelname+"_Et_wrtMu";
title = labelname+" E_t";
EtwrtMu = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
if ((v->getPath()).find("Jet") != std::string::npos
|| (v->getPath()).find("Mu") != std::string::npos) {
histoname = labelname+"_EtaPhi_wrtMu";
title = labelname+" #eta vs #phi";
EtaPhiwrtMu = dbe->book2D(histoname.c_str(),
title.c_str(),
nBins2D,-histEtaMax,histEtaMax,
nBins2D,-TMath::Pi(), TMath::Pi());
}
else if( ((v->getPath()).find("MET") != std::string::npos)
|| ((v->getPath()).find("SumET") != std::string::npos) ){
histoname = labelname+"_phi_wrtMu";
title = labelname+" #phi";
PhiwrtMu = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
v->setHistoswrtMu( NwrtMu, EtwrtMu, EtaPhiwrtMu, PhiwrtMu);
}
if (verbose_) std::cout << "Done booking histos wrt Muon " << std::endl;
}
if (plotEff_) {
// plot efficiency for specified HLT path pairs
if (verbose_) std::cout << " booking histos for Efficiency " << std::endl;
std::string foldernm = "/Efficiency/";
if (dbe) {
dbe->setCurrentFolder(dirname_ + foldernm);
}
// now loop over denom/num path pairs specified in cfg,
int countN = 0;
for (std::vector<std::pair<std::string, std::string> >::iterator
custompathnamepair = custompathnamepairs_.begin();
custompathnamepair != custompathnamepairs_.end();
++custompathnamepair) {
std::string denompathname = custompathnamepair->second;
std::string pathname = custompathnamepair->first;
int usedPrescale = prescUsed_[countN];
if (verbose_) std::cout << " ------prescale used -----" << usedPrescale << std::endl;
// check that these exist
bool foundfirst = false;
bool foundsecond = false;
for (unsigned int i=0; i!=n; ++i) {
if (hltConfig_.triggerName(i) == denompathname) foundsecond = true;
if (hltConfig_.triggerName(i) == pathname) foundfirst = true;
}
if (!foundfirst) {
edm::LogInfo("HLTJetMETDQMSource")
<< "pathname not found, ignoring "
<< pathname;
continue;
}
if (!foundsecond) {
edm::LogInfo("HLTJetMETDQMSource")
<< "denompathname not found, ignoring "
<< pathname;
continue;
}
//if (verbose_) std::cout << pathname << "\t" << denompathname << std::endl;
std::string l1pathname = "dummy";
int objectType = 0;
//int denomobjectType = 0;
//parse pathname to guess object type
if ( pathname.find("Mu") && (pathname=="HLT_L1MuOpen"))
objectType = trigger::TriggerMuon;
else if (pathname.find("MET") != std::string::npos) objectType = trigger::TriggerMET;
else if (pathname.find("L1MET") != std::string::npos) objectType = trigger::TriggerL1ETM;
else if (pathname.find("SumET") != std::string::npos) objectType = trigger::TriggerTET;
else if (pathname.find("Jet") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet30") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet50") != std::string::npos) objectType = trigger::TriggerJet;
//if ((pathname.find("HLT_Mu3") != std::string::npos) || (pathname.find("HLT_L2Mu9") != std::string::npos) ) objectType = trigger::TriggerMuon;
else continue;
// find L1 condition for numpath with numpath objecttype
// find PSet for L1 global seed for numpath,
// list module labels for numpath
std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
numpathmodule!= numpathmodules.end();
++numpathmodule ) {
// if (verbose_) std::cout << pathname << "\t" << *numpathmodule << "\t" << hltConfig_.moduleType(*numpathmodule) << std::endl;
if (hltConfig_.moduleType(*numpathmodule) == "HLTLevel1GTSeed"){
edm::ParameterSet l1GTPSet = hltConfig_.modulePSet(*numpathmodule);
// if (verbose_) std::cout << l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression") << std::endl;
// l1pathname = l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression");
l1pathname = *numpathmodule;
//if (verbose_) std::cout << *numpathmodule << std::endl;
break;
}
}
std::string filtername("dummy");
std::string Denomfiltername("denomdummy");
float ptMin = 0.0;
float ptMax = 300.0;
if (objectType == trigger::TriggerMuon) ptMax = 300.0;
if (objectType == trigger::TriggerJet) ptMax = 300.0;
if (objectType == trigger::TriggerMET) ptMax = 300.0;
if (objectType == trigger::TriggerTET) ptMax = 300.0;
if (objectType != 0){
if (verbose_) std::cout << " PathInfo(denompathname, pathname, l1pathname, filtername, Denomfiltername, processname_, objectType, ptMin, ptMax " << denompathname << " "<< pathname << " "<< l1pathname << " " << filtername << " " << Denomfiltername << " " << processname_ << " " << objectType << " " << ptMin << " " << ptMax<< std::endl;
hltPathsEff_.push_back(PathInfo(usedPrescale, denompathname, pathname, l1pathname, filtername, Denomfiltername, processname_, objectType, ptMin, ptMax));
}
countN++;
}
std::string histonm="JetMET_Efficiency_Summary";
std::string histonmDenom="Denom_passed_Summary";
std::string histonmNum="Num_passed_Summary";
rate_Denom = dbe->book1D(histonmDenom.c_str(),histonmDenom.c_str(),
hltPathsEff_.size(),0,hltPathsEff_.size());
rate_Num = dbe->book1D(histonmNum.c_str(),histonmNum.c_str(),
hltPathsEff_.size(),-0.5,hltPathsEff_.size()+1-0.5);
rate_Eff = dbe->book1D(histonm.c_str(),histonm.c_str(),
hltPathsEff_.size(),-0.5,hltPathsEff_.size()+1-0.5);
//rate_Eff = dbe_->bookProfile("Efficiency_Summary","Efficiency_Summary", hltPathsEff_.size(), -0.5, hltPathsEff_.size()-0.5, 1000, 0.0, 1.0);
int nname=0;
for(PathInfoCollection::iterator v = hltPathsEff_.begin(); v!= hltPathsEff_.end(); ++v ){
std::string labelnm("dummy");
std::string labeldenom("dummy");
labelnm = v->getPath();
labeldenom = v->getDenomPath();
//rate_Eff->getTProfile()->GetXaxis()->SetBinLabel(nname+1,labelnm.c_str());
rate_Eff->setBinLabel(nname+1,labelnm);
rate_Denom->setBinLabel(nname+1,labeldenom);
rate_Num->setBinLabel(nname+1,labelnm);
nname++;
}
for(PathInfoCollection::iterator v = hltPathsEff_.begin(); v!= hltPathsEff_.end(); ++v ) {
MonitorElement *NEff=0;
MonitorElement *EtEff=0;
MonitorElement *EtaEff=0;
MonitorElement *PhiEff=0;
MonitorElement *NNum =0;
MonitorElement *EtNum =0;
MonitorElement *EtaNum =0;
MonitorElement *PhiNum =0;
MonitorElement *NDenom=0;
MonitorElement *EtDenom=0;
MonitorElement *EtaDenom=0;
MonitorElement *PhiDenom=0;
std::string labelname("dummy");
labelname = "Eff_" + v->getPath() + "_wrt_" + v->getDenomPath();
std::string histoname(labelname+"");
std::string title(labelname+"");
double histEtaMax = 5.0;
if (v->getObjectType() == trigger::TriggerMuon
|| v->getObjectType() == trigger::TriggerL1Mu) {
histEtaMax = muonEtaMax_; nBins_ = 20 ;
}
else if (v->getObjectType() == trigger::TriggerJet
|| v->getObjectType() == trigger::TriggerL1CenJet
|| v->getObjectType() == trigger::TriggerL1ForJet ){
histEtaMax = jetEtaMax_; nBins_ = 60 ;
}
else if (v->getObjectType() == trigger::TriggerMET
|| v->getObjectType() == trigger::TriggerL1ETM ) {
histEtaMax = 5.0; nBins_ = 60 ;
}
TString pathfolder = dirname_ + foldernm + v->getPath();
dbe_->setCurrentFolder(pathfolder.Data());
if (verbose_) std::cout << "Booking Histos in Directory " << pathfolder.Data() << std::endl;
//pathfolder = dirname_ + TString("/Eff/") + v->getPath();
//int nBins2D = 10;
histoname = labelname+"_Et_Eff";
title = labelname+" E_t Eff";
EtEff = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
histoname = labelname+"_Et_Num";
title = labelname+" E_t Num";
EtNum = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
histoname = labelname+"_Et_Denom";
title = labelname+" E_t Denom";
EtDenom = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
if ((v->getPath()).find("Jet") != std::string::npos
|| (v->getPath()).find("Mu") != std::string::npos) {
histoname = labelname+"_Eta_Eff";
title = labelname+" #eta Eff";
EtaEff = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
histoname = labelname+"_Phi_Eff";
title = labelname+" #phi Eff";
PhiEff = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Eta_Num";
title = labelname+" #eta Num";
EtaNum = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
histoname = labelname+"_Phi_Num";
title = labelname+" #phi Num";
PhiNum = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Eta_Denom";
title = labelname+" #eta Denom";
EtaDenom = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
histoname = labelname+"_Phi_Denom";
title = labelname+" #phi Denom";
PhiDenom = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
else if( ((v->getPath()).find("MET") != std::string::npos)
|| ((v->getPath()).find("SumET") != std::string::npos) ){
histoname = labelname+"_Phi_Eff";
title = labelname+" #phi Eff";
PhiEff = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Phi_Num";
title = labelname+" #phi Num";
PhiNum = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Phi_Denom";
title = labelname+" #phi Denom";
PhiDenom = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
v->setHistosEff( NEff, EtEff, EtaEff, PhiEff, NNum, EtNum, EtaNum, PhiNum, NDenom, EtDenom, EtaDenom, PhiDenom);
}
if (verbose_) std::cout << "Done booking histos for Efficiency " << std::endl;
}
*/
}
| void HLTJetMETDQMSource::beginLuminosityBlock | ( | const edm::LuminosityBlock & | lumiSeg, |
| const edm::EventSetup & | c | ||
| ) | [private, virtual] |
| void HLTJetMETDQMSource::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | c | ||
| ) | [private, virtual] |
from begin job
Reimplemented from edm::EDAnalyzer.
Definition at line 1521 of file HLTJetMETDQMSource.cc.
References DQMStore::book1D(), DQMStore::book2D(), gather_cfg::cout, custompathnamemu_, custompathnamepairs_, dbe_, dirname_, reco::tau::disc::Eta(), spr::find(), SingleJetValidation_cfi::foldernm, hltConfig_, hltconfigchanged, hltPathsAll_, hltPathsEff_, hltPathswrtMu_, i, edm::RunBase::id(), HLTConfigProvider::init(), j, jetEtaMax_, HLTConfigProvider::moduleLabels(), HLTConfigProvider::modulePSet(), HLTConfigProvider::moduleType(), muonEtaMax_, N, n, nBins_, lumiNorm::nname, cmsCodeRules::cppFunctionSkipper::operator, timingPdfMaker::pathname, colinearityKinematic::Phi, Pi, plotAll_, plotEff_, plotwrtMu_, prescUsed_, processname_, jptDQMConfig_cff::ptMax, PtMinSelector_cfg::ptMin, rate_Denom, rate_Eff, rate_Num, MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), HLTConfigProvider::size(), indexGen::title, trigger::TriggerJet, trigger::TriggerL1CenJet, trigger::TriggerL1ETM, trigger::TriggerL1ForJet, trigger::TriggerL1Mu, trigger::TriggerMET, trigger::TriggerMuon, HLTConfigProvider::triggerName(), trigger::TriggerTET, v, and verbose_.
{
if (verbose_) std::cout << "beginRun, run " << run.id() << std::endl;
/*
// HLT config does not change within runs!
bool changed(true);
//if (hltConfig_.init(iEvent,processname_,changed)) {
if (hltConfig_.init(run,c,processname_,changed)) {
if (changed) {
hltconfigchanged=true;
//if(verbose_)
//hltConfig_.dump("Triggers");
LogWarning("HLTJetMETDQMSource") << "HLTJetMETDQMSource:analyze: The number of valid triggers has changed since beginning of job." << std::endl;
// << "Processing of events halted for summary histograms" << std::endl;
//<< "Summary histograms do not support changing configurations." << std::endl
}
}
else {
LogError("HLTJetMETDQMSource") << "HLT config init error" << std::endl;
return;
}
if ( hltConfig_.size() <= 0 ) {
LogError("HLTJetMETDQMSource") << "HLT config size error" << std::endl;
return;
}
// if ( currun != runno) runchanged=true;
*/
hltconfigchanged=false;
//if (!hltConfig_.init(processname_)) {
if (!hltConfig_.init(run, c, processname_, hltconfigchanged)) {
processname_ = "FU";
// if (!hltConfig_.init(processname_)){
if (!hltConfig_.init(run, c, processname_, hltconfigchanged)) {
LogError("HLTJetMETDQMSource") << "HLTConfigProvider failed to initialize." << std::endl;
}
// hltConfig_.dump("Triggers");
}
DQMStore *dbe = 0;
dbe = Service<DQMStore>().operator->();
if (dbe) {
dbe->setCurrentFolder(dirname_);
}
const unsigned int n(hltConfig_.size());
if (verbose_) std::cout << " hltConfig_.size() " << n << std::endl;
//-------------if plot all jetmet trigger pt,eta, phi-------------------
if (plotAll_){
if (verbose_) std::cout << " booking histos All " << std::endl;
std::string foldernm = "/All/";
if (dbe) {
dbe->setCurrentFolder(dirname_ + foldernm);
}
for (unsigned int j=0; j!=n; ++j) {
std::string pathname = hltConfig_.triggerName(j);
std::string l1pathname = "dummy";
//for (unsigned int i=0; i!=n; ++i) {
// std::cout << hltConfig_.triggerName(i) << std::endl;
//std::string denompathname = hltConfig_.triggerName(i);
std::string denompathname = "";
std::string denompathnamemu = custompathnamemu_;
unsigned int usedPresscale = 1;
unsigned int objectType = 0;
//int denomobjectType = 0;
//parse pathname to guess object type
if ( pathname.find("Mu") && (pathname==denompathnamemu))
objectType = trigger::TriggerMuon;
else if (pathname.find("MET") != std::string::npos) objectType = trigger::TriggerMET;
else if (pathname.find("L1MET") != std::string::npos) objectType = trigger::TriggerL1ETM;
else if (pathname.find("SumET") != std::string::npos) objectType = trigger::TriggerTET;
else if (pathname.find("Jet") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet30") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet50") != std::string::npos) objectType = trigger::TriggerJet;
//if ((pathname.find("HLT_Mu3") != std::string::npos) || (pathname.find("HLT_L2Mu9") != std::string::npos) ) objectType = trigger::TriggerMuon;
else continue;
//std::cout << "objecttye " << objectType << std::endl;
// find L1 condition for numpath with numpath objecttype
// find PSet for L1 global seed for numpath,
// list module labels for numpath
std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
numpathmodule!= numpathmodules.end();
++numpathmodule ) {
// std::cout << pathname << "\t" << *numpathmodule << "\t"
//<< hltConfig_.moduleType(*numpathmodule) << std::endl;
if (hltConfig_.moduleType(*numpathmodule) == "HLTLevel1GTSeed") {
edm::ParameterSet l1GTPSet = hltConfig_.modulePSet(*numpathmodule);
//std::cout << l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression") << std::endl;
// l1pathname = l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression");
l1pathname = *numpathmodule;
break;
}
}
std::string filtername("dummy");
std::string Denomfiltername("denomdummy");
float ptMin = 0.0;
float ptMax = 300.0;
if (objectType != 0 )
hltPathsAll_.push_back(PathInfo(usedPresscale, denompathname,
pathname, l1pathname, filtername,
Denomfiltername, processname_,
objectType, ptMin, ptMax));
}
std::string histonm="JetMET_rate_All";
std::string histonmL1="JetMET_rate_All_L1";
std::string histot="JetMET Rate Summary";
std::string histoL1t="JetMET L1 Rate Summary";
//fill later//rate_All = dbe->book1D(histonm.c_str(),histot.c_str(),
//fill later// hltPathsAll_.size()+1,-0.5,hltPathsAll_.size()+1-0.5);
//fill later//rate_All_L1 = dbe->book1D(histonmL1.c_str(),histoL1t.c_str(),
//fill later// hltPathsAll_.size()+1,-0.5,hltPathsAll_.size()+1-0.5);
//rate_All->setBinLabel(hltPathsAll_.size()+1,"Rate",1);
unsigned int nname=0;
unsigned int nnameL1=0;
for(PathInfoCollection::iterator v = hltPathsAll_.begin(); v!= hltPathsAll_.end(); ++v ){
std::string labelnm("dummy");
labelnm = v->getPath();
//fill later//rate_All->setBinLabel(nname+1,labelnm);
nname++;
std::string labelnml1("dummyl1");
labelnml1 = v->getl1Path();
//fill later//rate_All_L1->setBinLabel(nnameL1+1,labelnml1);
nnameL1++;
}
// now set up all of the histos for each path
for(PathInfoCollection::iterator v = hltPathsAll_.begin(); v!= hltPathsAll_.end(); ++v ) {
MonitorElement *N = 0;
MonitorElement *Et = 0;
MonitorElement *EtaPhi = 0;
MonitorElement *Eta = 0;
MonitorElement *Phi = 0;
MonitorElement *NL1 = 0;
MonitorElement *l1Et = 0;
MonitorElement *l1EtaPhi = 0;
MonitorElement *l1Eta = 0;
MonitorElement *l1Phi = 0;
std::string labelname("dummy");
labelname = v->getPath();
std::string histoname(labelname+"");
std::string title(labelname+"");
double histEtaMax = 2.5;
if (v->getObjectType() == trigger::TriggerMuon
|| v->getObjectType() == trigger::TriggerL1Mu) {
histEtaMax = muonEtaMax_;
nBins_ = 25 ;
}
else if (v->getObjectType() == trigger::TriggerJet
|| v->getObjectType() == trigger::TriggerL1CenJet
|| v->getObjectType() == trigger::TriggerL1ForJet ) {
histEtaMax = jetEtaMax_;
nBins_ = 60 ;
}
else if (v->getObjectType() == trigger::TriggerMET
|| v->getObjectType() == trigger::TriggerL1ETM ) {
histEtaMax = 5.0;
nBins_ = 60 ;
}
TString pathfolder = dirname_ + foldernm + v->getPath();
dbe_->setCurrentFolder(pathfolder.Data());
if (verbose_) std::cout << "Booking Histos in Directory " << pathfolder.Data() << std::endl;
int nBins2D = 10;
histoname = labelname+"_Et";
title = labelname+" E_t";
Et = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
histoname = labelname+"_l1Et";
title = labelname+" L1 E_t";
l1Et = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
if (labelname.find("Jet") != std::string::npos
|| labelname.find("Mu") != std::string::npos) {
histoname = labelname+"_EtaPhi";
title = labelname+" #eta vs #phi";
EtaPhi = dbe->book2D(histoname.c_str(),
title.c_str(),
nBins2D,-histEtaMax,histEtaMax,
nBins2D,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_l1EtaPhi";
title = labelname+"L1 #eta vs L1 #phi";
l1EtaPhi = dbe->book2D(histoname.c_str(),
title.c_str(),
nBins2D,-histEtaMax,histEtaMax,
nBins2D,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Phi";
title = labelname+" #phi";
Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_l1Phi";
title = labelname+"L1 #phi";
l1Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Eta";
title = labelname+" #eta";
Eta = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax
);
histoname = labelname+"_l1Eta";
title = labelname+"L1 #eta";
l1Eta = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
}
else if( (labelname.find("MET") != std::string::npos)
|| (labelname.find("SumET") != std::string::npos) ){
histoname = labelname+"_phi";
title = labelname+" #phi";
Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_l1Phi";
title = labelname+"L1 #phi";
l1Phi = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
v->setHistos( N, Et, EtaPhi, Eta, Phi, NL1, l1Et, l1EtaPhi,l1Eta, l1Phi);
}
if (verbose_) std::cout << "Done booking histos All " << std::endl;
}
//----------------plot all jetmet trigger wrt some muon trigger-----
if (plotwrtMu_){
if (verbose_) std::cout << " booking histos wrt Muon " << std::endl;
std::string foldernm = "/wrtMuon/";
if (dbe) {
dbe->setCurrentFolder(dirname_ + foldernm);
}
for (unsigned int j=0; j!=n; ++j) {
std::string pathname = hltConfig_.triggerName(j);
std::string l1pathname = "dummy";
//for (unsigned int i=0; i!=n; ++i)
//{
if (verbose_) std::cout << hltConfig_.triggerName(j) << std::endl;
std::string denompathname = custompathnamemu_ ;
int objectType = 0;
int usedPresscale = 1;
//int denomobjectType = 0;
//parse pathname to guess object type
if ( pathname.find("Mu") && (pathname==denompathname))
objectType = trigger::TriggerMuon;
else if (pathname.find("MET") != std::string::npos) objectType = trigger::TriggerMET;
else if (pathname.find("L1MET") != std::string::npos) objectType = trigger::TriggerL1ETM;
else if (pathname.find("SumET") != std::string::npos) objectType = trigger::TriggerTET;
else if (pathname.find("Jet") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet30") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet50") != std::string::npos) objectType = trigger::TriggerJet;
//if ((pathname.find("HLT_Mu3") != std::string::npos) || (pathname.find("HLT_L2Mu9") != std::string::npos) ) objectType = trigger::TriggerMuon;
else continue;
// find L1 condition for numpath with numpath objecttype
// find PSet for L1 global seed for numpath,
// list module labels for numpath
std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
// std::cout << pathname << "\t" << *numpathmodule << "\t" << hltConfig_.moduleType(*numpathmodule) << std::endl;
if (hltConfig_.moduleType(*numpathmodule) == "HLTLevel1GTSeed") {
edm::ParameterSet l1GTPSet = hltConfig_.modulePSet(*numpathmodule);
//std::cout << l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression") << std::endl;
// l1pathname = l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression");
l1pathname = *numpathmodule;
break;
}
}
std::string filtername("dummy");
std::string Denomfiltername("denomdummy");
float ptMin = 0.0;
float ptMax = 300.0;
if ( objectType != 0){
// if (verbose_) std::cout << " wrt muon PathInfo(denompathname, pathname, l1pathname, filtername, denomfiltername, processname_, objectType, ptMin, ptMax " << denompathname << " "<< pathname << " "<< l1pathname << " " << filtername << " " << Denomfiltername << " " << processname_ << " " << objectType << " " << ptMin << " " << ptMax<< std::endl;
hltPathswrtMu_.push_back(PathInfo(usedPresscale, denompathname, pathname, l1pathname,
filtername, Denomfiltername, processname_,
objectType, ptMin, ptMax));
}
}
std::string histonm="JetMET_rate_wrt_" + custompathnamemu_ + "_Summary";
std::string histt="JetMET Rate wrt " + custompathnamemu_ + "Summary";
//fill later//rate_wrtMu = dbe->book1D(histonm.c_str(),histt.c_str(),
//fill later// hltPathswrtMu_.size()+1,-0.5,hltPathswrtMu_.size()+1.0-0.5);
int nname=0;
for(PathInfoCollection::iterator v = hltPathswrtMu_.begin(); v!= hltPathswrtMu_.end(); ++v ){
std::string labelnm("dummy");
labelnm = v->getPath();
//fill later//rate_wrtMu->setBinLabel(nname+1,labelnm);
nname++;
}
// now set up all of the histos for each path
for(PathInfoCollection::iterator v = hltPathswrtMu_.begin(); v!= hltPathswrtMu_.end(); ++v ) {
MonitorElement *NwrtMu = 0;
MonitorElement *EtwrtMu = 0;
MonitorElement *EtaPhiwrtMu = 0;
MonitorElement *PhiwrtMu = 0;
std::string labelname("dummy");
labelname = v->getPath() + "_wrt_" + v->getDenomPath();
std::string histoname(labelname+"");
std::string title(labelname+"");
double histEtaMax = 2.5;
if (v->getObjectType() == trigger::TriggerMuon
|| v->getObjectType() == trigger::TriggerL1Mu) {
histEtaMax = muonEtaMax_;nBins_ = 20 ;
}
else if (v->getObjectType() == trigger::TriggerJet
|| v->getObjectType() == trigger::TriggerL1CenJet
|| v->getObjectType() == trigger::TriggerL1ForJet ){
histEtaMax = jetEtaMax_; nBins_ = 60 ;
}
else if (v->getObjectType() == trigger::TriggerMET
|| v->getObjectType() == trigger::TriggerL1ETM ) {
histEtaMax = 5.0; nBins_ = 60 ;
}
TString pathfolder = dirname_ + foldernm + v->getPath();
dbe_->setCurrentFolder(pathfolder.Data());
if (verbose_) std::cout << "Booking Histos in Directory " << pathfolder.Data() << std::endl;
int nBins2D = 10;
//pathfolder = dirname_ + TString("/wrtMuon/") + v->getPath();
histoname = labelname+"_Et_wrtMu";
title = labelname+" E_t";
EtwrtMu = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
if ((v->getPath()).find("Jet") != std::string::npos
|| (v->getPath()).find("Mu") != std::string::npos) {
histoname = labelname+"_EtaPhi_wrtMu";
title = labelname+" #eta vs #phi";
EtaPhiwrtMu = dbe->book2D(histoname.c_str(),
title.c_str(),
nBins2D,-histEtaMax,histEtaMax,
nBins2D,-TMath::Pi(), TMath::Pi());
}
else if( ((v->getPath()).find("MET") != std::string::npos)
|| ((v->getPath()).find("SumET") != std::string::npos) ){
histoname = labelname+"_phi_wrtMu";
title = labelname+" #phi";
PhiwrtMu = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
v->setHistoswrtMu( NwrtMu, EtwrtMu, EtaPhiwrtMu, PhiwrtMu);
}
if (verbose_) std::cout << "Done booking histos wrt Muon " << std::endl;
}
if (plotEff_) {
// plot efficiency for specified HLT path pairs
if (verbose_) std::cout << " booking histos for Efficiency " << std::endl;
std::string foldernm = "/Efficiency/";
if (dbe) {
dbe->setCurrentFolder(dirname_ + foldernm);
}
// now loop over denom/num path pairs specified in cfg,
int countN = 0;
for (std::vector<std::pair<std::string, std::string> >::iterator
custompathnamepair = custompathnamepairs_.begin();
custompathnamepair != custompathnamepairs_.end();
++custompathnamepair) {
std::string denompathname = custompathnamepair->second;
std::string denompathnamemu = custompathnamemu_;
std::string pathname = custompathnamepair->first;
int usedPrescale = prescUsed_[countN];
if (verbose_) std::cout << " ------prescale used -----" << usedPrescale << std::endl;
// check that these exist
bool foundfirst = false;
bool foundsecond = false;
for (unsigned int i=0; i!=n; ++i) {
if (hltConfig_.triggerName(i) == denompathname) foundsecond = true;
if (hltConfig_.triggerName(i) == pathname) foundfirst = true;
}
if (!foundfirst) {
edm::LogInfo("HLTJetMETDQMSource")
<< "pathname not found, ignoring "
<< pathname;
continue;
}
if (!foundsecond) {
edm::LogInfo("HLTJetMETDQMSource")
<< "denompathname not found, ignoring "
<< pathname;
continue;
}
//if (verbose_) std::cout << pathname << "\t" << denompathname << std::endl;
std::string l1pathname = "dummy";
int objectType = 0;
//int denomobjectType = 0;
//parse pathname to guess object type
if ( pathname.find("Mu") && (pathname==denompathnamemu))
objectType = trigger::TriggerMuon;
else if (pathname.find("MET") != std::string::npos) objectType = trigger::TriggerMET;
else if (pathname.find("L1MET") != std::string::npos) objectType = trigger::TriggerL1ETM;
else if (pathname.find("SumET") != std::string::npos) objectType = trigger::TriggerTET;
else if (pathname.find("Jet") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet30") != std::string::npos) objectType = trigger::TriggerJet;
//if (pathname.find("HLT_Jet50") != std::string::npos) objectType = trigger::TriggerJet;
//if ((pathname.find("HLT_Mu3") != std::string::npos) || (pathname.find("HLT_L2Mu9") != std::string::npos) ) objectType = trigger::TriggerMuon;
else continue;
// find L1 condition for numpath with numpath objecttype
// find PSet for L1 global seed for numpath,
// list module labels for numpath
std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
numpathmodule!= numpathmodules.end();
++numpathmodule ) {
// if (verbose_) std::cout << pathname << "\t" << *numpathmodule << "\t" << hltConfig_.moduleType(*numpathmodule) << std::endl;
if (hltConfig_.moduleType(*numpathmodule) == "HLTLevel1GTSeed"){
edm::ParameterSet l1GTPSet = hltConfig_.modulePSet(*numpathmodule);
// if (verbose_) std::cout << l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression") << std::endl;
// l1pathname = l1GTPSet.getParameter<std::string>("L1SeedsLogicalExpression");
l1pathname = *numpathmodule;
//if (verbose_) std::cout << *numpathmodule << std::endl;
break;
}
}
std::string filtername("dummy");
std::string Denomfiltername("denomdummy");
float ptMin = 0.0;
float ptMax = 300.0;
if (objectType == trigger::TriggerMuon) ptMax = 300.0;
if (objectType == trigger::TriggerJet) ptMax = 300.0;
if (objectType == trigger::TriggerMET) ptMax = 300.0;
if (objectType == trigger::TriggerTET) ptMax = 300.0;
if (objectType != 0){
if (verbose_) std::cout << " PathInfo(denompathname, pathname, l1pathname, filtername, Denomfiltername, processname_, objectType, ptMin, ptMax " << denompathname << " "<< pathname << " "<< l1pathname << " " << filtername << " " << Denomfiltername << " " << processname_ << " " << objectType << " " << ptMin << " " << ptMax<< std::endl;
hltPathsEff_.push_back(PathInfo(usedPrescale, denompathname, pathname, l1pathname, filtername, Denomfiltername, processname_, objectType, ptMin, ptMax));
}
countN++;
}
std::string histonm="JetMET_Efficiency_Summary";
std::string histonmDenom="Denom_passed_Summary";
std::string histonmNum="Num_passed_Summary";
rate_Denom = dbe->book1D(histonmDenom.c_str(),histonmDenom.c_str(),
hltPathsEff_.size(),0,hltPathsEff_.size());
rate_Num = dbe->book1D(histonmNum.c_str(),histonmNum.c_str(),
hltPathsEff_.size(),-0.5,hltPathsEff_.size()+1-0.5);
rate_Eff = dbe->book1D(histonm.c_str(),histonm.c_str(),
hltPathsEff_.size(),-0.5,hltPathsEff_.size()+1-0.5);
//rate_Eff = dbe_->bookProfile("Efficiency_Summary","Efficiency_Summary", hltPathsEff_.size(), -0.5, hltPathsEff_.size()-0.5, 1000, 0.0, 1.0);
int nname=0;
for(PathInfoCollection::iterator v = hltPathsEff_.begin(); v!= hltPathsEff_.end(); ++v ){
std::string labelnm("dummy");
std::string labeldenom("dummy");
labelnm = v->getPath();
labeldenom = v->getDenomPath();
//rate_Eff->getTProfile()->GetXaxis()->SetBinLabel(nname+1,labelnm.c_str());
rate_Eff->setBinLabel(nname+1,labelnm);
rate_Denom->setBinLabel(nname+1,labeldenom);
rate_Num->setBinLabel(nname+1,labelnm);
nname++;
}
for(PathInfoCollection::iterator v = hltPathsEff_.begin(); v!= hltPathsEff_.end(); ++v ) {
MonitorElement *NEff=0;
MonitorElement *EtEff=0;
MonitorElement *EtaEff=0;
MonitorElement *PhiEff=0;
MonitorElement *NNum =0;
MonitorElement *EtNum =0;
MonitorElement *EtaNum =0;
MonitorElement *PhiNum =0;
MonitorElement *NDenom=0;
MonitorElement *EtDenom=0;
MonitorElement *EtaDenom=0;
MonitorElement *PhiDenom=0;
std::string labelname("dummy");
labelname = "Eff_" + v->getPath() + "_wrt_" + v->getDenomPath();
std::string histoname(labelname+"");
std::string title(labelname+"");
double histEtaMax = 5.0;
if (v->getObjectType() == trigger::TriggerMuon
|| v->getObjectType() == trigger::TriggerL1Mu) {
histEtaMax = muonEtaMax_; nBins_ = 20 ;
}
else if (v->getObjectType() == trigger::TriggerJet
|| v->getObjectType() == trigger::TriggerL1CenJet
|| v->getObjectType() == trigger::TriggerL1ForJet ){
histEtaMax = jetEtaMax_; nBins_ = 60 ;
}
else if (v->getObjectType() == trigger::TriggerMET
|| v->getObjectType() == trigger::TriggerL1ETM ) {
histEtaMax = 5.0; nBins_ = 60 ;
}
TString pathfolder = dirname_ + foldernm + v->getPath();
dbe_->setCurrentFolder(pathfolder.Data());
if (verbose_) std::cout << "Booking Histos in Directory " << pathfolder.Data() << std::endl;
//pathfolder = dirname_ + TString("/Eff/") + v->getPath();
//int nBins2D = 10;
histoname = labelname+"_Et_Eff";
title = labelname+" E_t Eff";
EtEff = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
histoname = labelname+"_Et_Num";
title = labelname+" E_t Num";
EtNum = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
histoname = labelname+"_Et_Denom";
title = labelname+" E_t Denom";
EtDenom = dbe->book1D(histoname.c_str(),
title.c_str(),nBins_,
v->getPtMin(),
v->getPtMax());
if ((v->getPath()).find("Jet") != std::string::npos
|| (v->getPath()).find("Mu") != std::string::npos) {
histoname = labelname+"_Eta_Eff";
title = labelname+" #eta Eff";
EtaEff = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
histoname = labelname+"_Phi_Eff";
title = labelname+" #phi Eff";
PhiEff = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Eta_Num";
title = labelname+" #eta Num";
EtaNum = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
histoname = labelname+"_Phi_Num";
title = labelname+" #phi Num";
PhiNum = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Eta_Denom";
title = labelname+" #eta Denom";
EtaDenom = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-histEtaMax,histEtaMax);
histoname = labelname+"_Phi_Denom";
title = labelname+" #phi Denom";
PhiDenom = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
else if( ((v->getPath()).find("MET") != std::string::npos)
|| ((v->getPath()).find("SumET") != std::string::npos) ){
histoname = labelname+"_Phi_Eff";
title = labelname+" #phi Eff";
PhiEff = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Phi_Num";
title = labelname+" #phi Num";
PhiNum = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
histoname = labelname+"_Phi_Denom";
title = labelname+" #phi Denom";
PhiDenom = dbe->book1D(histoname.c_str(),
title.c_str(),
nBins_,-TMath::Pi(), TMath::Pi());
}
v->setHistosEff( NEff, EtEff, EtaEff, PhiEff, NNum, EtNum, EtaNum, PhiNum, NDenom, EtDenom, EtaDenom, PhiDenom);
}
if (verbose_) std::cout << "Done booking histos for Efficiency " << std::endl;
}
//end from begin job
if (verbose_) std::cout << "End BeginRun ---------------- " << std::endl;
}
| void HLTJetMETDQMSource::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 2195 of file HLTJetMETDQMSource.cc.
References nev_.
| void HLTJetMETDQMSource::endLuminosityBlock | ( | const edm::LuminosityBlock & | lumiSeg, |
| const edm::EventSetup & | c | ||
| ) | [private, virtual] |
DQM Client Diagnostic.
Reimplemented from edm::EDAnalyzer.
Definition at line 2190 of file HLTJetMETDQMSource.cc.
{
}
| void HLTJetMETDQMSource::endRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | c | ||
| ) | [private, virtual] |
EndRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 2201 of file HLTJetMETDQMSource.cc.
References gather_cfg::cout, edm::RunBase::id(), and verbose_.
| void HLTJetMETDQMSource::histobooking | ( | const edm::EventSetup & | c | ) | [private] |
int HLTJetMETDQMSource::currentRun_ [private] |
Definition at line 76 of file HLTJetMETDQMSource.h.
std::string HLTJetMETDQMSource::custompathnamemu_ [private] |
Definition at line 96 of file HLTJetMETDQMSource.h.
Referenced by beginRun(), and HLTJetMETDQMSource().
std::vector<std::pair<std::string, std::string> > HLTJetMETDQMSource::custompathnamepairs_ [private] |
Definition at line 97 of file HLTJetMETDQMSource.h.
Referenced by beginRun(), and HLTJetMETDQMSource().
DQMStore* HLTJetMETDQMSource::dbe_ [private] |
Definition at line 68 of file HLTJetMETDQMSource.h.
Referenced by beginRun(), and HLTJetMETDQMSource().
std::string HLTJetMETDQMSource::dirname_ [private] |
Definition at line 101 of file HLTJetMETDQMSource.h.
Referenced by beginJob(), beginRun(), and HLTJetMETDQMSource().
Definition at line 109 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and beginRun().
bool HLTJetMETDQMSource::hltconfigchanged [private] |
Definition at line 78 of file HLTJetMETDQMSource.h.
Referenced by beginRun().
Definition at line 401 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and beginRun().
Definition at line 402 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and beginRun().
Definition at line 403 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and beginRun().
double HLTJetMETDQMSource::htMin_ [private] |
Definition at line 93 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
double HLTJetMETDQMSource::jetDRMatch_ [private] |
Definition at line 90 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
double HLTJetMETDQMSource::jetEtaMax_ [private] |
Definition at line 88 of file HLTJetMETDQMSource.h.
Referenced by beginRun(), and HLTJetMETDQMSource().
double HLTJetMETDQMSource::jetEtMin_ [private] |
Definition at line 89 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
double HLTJetMETDQMSource::metMin_ [private] |
Definition at line 92 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
bool HLTJetMETDQMSource::monitorDaemon_ [private] |
Definition at line 104 of file HLTJetMETDQMSource.h.
double HLTJetMETDQMSource::muonDRMatch_ [private] |
Definition at line 86 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
double HLTJetMETDQMSource::muonEtaMax_ [private] |
Definition at line 84 of file HLTJetMETDQMSource.h.
Referenced by beginRun(), and HLTJetMETDQMSource().
double HLTJetMETDQMSource::muonEtMin_ [private] |
Definition at line 85 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
unsigned int HLTJetMETDQMSource::nBins_ [private] |
Definition at line 79 of file HLTJetMETDQMSource.h.
Referenced by beginRun(), and HLTJetMETDQMSource().
int HLTJetMETDQMSource::nev_ [private] |
Definition at line 67 of file HLTJetMETDQMSource.h.
Referenced by analyze(), beginJob(), and endJob().
int HLTJetMETDQMSource::nhltcfg [private] |
Definition at line 77 of file HLTJetMETDQMSource.h.
bool HLTJetMETDQMSource::plotAll_ [private] |
Definition at line 72 of file HLTJetMETDQMSource.h.
Referenced by analyze(), beginRun(), and HLTJetMETDQMSource().
bool HLTJetMETDQMSource::plotEff_ [private] |
Definition at line 73 of file HLTJetMETDQMSource.h.
Referenced by analyze(), beginRun(), and HLTJetMETDQMSource().
bool HLTJetMETDQMSource::plotwrtMu_ [private] |
Definition at line 74 of file HLTJetMETDQMSource.h.
Referenced by analyze(), beginRun(), and HLTJetMETDQMSource().
std::vector<int> HLTJetMETDQMSource::prescUsed_ [private] |
Definition at line 98 of file HLTJetMETDQMSource.h.
Referenced by beginRun(), and HLTJetMETDQMSource().
std::string HLTJetMETDQMSource::processname_ [private] |
Definition at line 102 of file HLTJetMETDQMSource.h.
Referenced by analyze(), beginRun(), and HLTJetMETDQMSource().
double HLTJetMETDQMSource::ptMax_ [private] |
Definition at line 81 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
double HLTJetMETDQMSource::ptMin_ [private] |
Definition at line 80 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
MonitorElement* HLTJetMETDQMSource::rate_All [private] |
Definition at line 404 of file HLTJetMETDQMSource.h.
MonitorElement* HLTJetMETDQMSource::rate_All_L1 [private] |
Definition at line 405 of file HLTJetMETDQMSource.h.
MonitorElement* HLTJetMETDQMSource::rate_Denom [private] |
Definition at line 407 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and beginRun().
MonitorElement* HLTJetMETDQMSource::rate_Eff [private] |
Definition at line 406 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and beginRun().
MonitorElement* HLTJetMETDQMSource::rate_Num [private] |
Definition at line 408 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and beginRun().
MonitorElement* HLTJetMETDQMSource::rate_wrtMu [private] |
Definition at line 409 of file HLTJetMETDQMSource.h.
bool HLTJetMETDQMSource::resetMe_ [private] |
Definition at line 75 of file HLTJetMETDQMSource.h.
MonitorElement* HLTJetMETDQMSource::scalersSelect [private] |
Definition at line 111 of file HLTJetMETDQMSource.h.
double HLTJetMETDQMSource::sumEtMin_ [private] |
Definition at line 94 of file HLTJetMETDQMSource.h.
Referenced by HLTJetMETDQMSource().
int HLTJetMETDQMSource::theHLTOutputType [private] |
Definition at line 105 of file HLTJetMETDQMSource.h.
MonitorElement* HLTJetMETDQMSource::total_ [private] |
Definition at line 70 of file HLTJetMETDQMSource.h.
Definition at line 108 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and HLTJetMETDQMSource().
Definition at line 107 of file HLTJetMETDQMSource.h.
Referenced by analyze(), and HLTJetMETDQMSource().
bool HLTJetMETDQMSource::verbose_ [private] |
Definition at line 103 of file HLTJetMETDQMSource.h.
Referenced by analyze(), beginRun(), endRun(), and HLTJetMETDQMSource().
1.7.3