#include <TrigResRateMon.h>
Module to monitor rates from TriggerResults
Implementation: <Notes on="" implementation>="">
Definition at line 86 of file TrigResRateMon.h.
TrigResRateMon::TrigResRateMon | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 17 of file TrigResRateMon.cc.
References DQMStore::book1D(), gather_cfg::cout, custompathnamepairs_, dbe_, dirname_, doCombineRuns_, doVBTFMuon_, dRMax_, dRMaxElectronMuon_, dxyCut_, fCustomBXPath, fIsSetup, fLumiFlag, found_zbIndex, foundReferenceTrigger_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltPaths_, isAlsoTrackerMuon_, jmsDebug, jmsFakeZBCounts, LogDebug, LSsize_, maskedPaths_, ME_HLT_BX, ME_HLT_CUSTOM_BX, ME_HLTAll_LS, muonHitsCut_, nBins2D_, nBins_, nBinsDR_, nBinsOneOverEt_, Nbx_, nLS_, nMatchesCut_, normalizedChi2Cut_, NULL, cppFunctionSkipper::operator, EgammaValidation_cff::paths, pathsIndividualHLTPathsPerLSFolder_, pathsSummaryFilterCountsFolder_, pathsSummaryFolder_, pathsSummaryHLTCorrelationsFolder_, pathsSummaryHLTPathsPerBXFolder_, pathsSummaryHLTPathsPerLSFolder_, pathsSummaryStreamsFolder_, pixelHitsCut_, plotAll_, processname_, ptMax_, ptMin_, recHitsEBTag_, recHitsEETag_, referenceBX_, referenceTrigInput_, scalersSelect, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), specialPaths_, AlCaHLTBitMon_QueryRunRegistry::string, sumEtMin_, testPaths_, testPathsFolder_, thresholdFactor_, trackerHitsCut_, triggerResultsLabel_, and triggerSummaryLabel_.
: currentRun_(-99) { LogDebug("TrigResRateMon") << "constructor...." ; fIsSetup = false; dbe_ = Service < DQMStore > ().operator->(); if ( ! dbe_ ) { LogInfo("TrigResRateMon") << "unabel to get DQMStore service?"; } if (iConfig.getUntrackedParameter < bool > ("DQMStore", false)) { dbe_->setVerbose(0); } dirname_ = iConfig.getUntrackedParameter("dirname", std::string("HLT/TrigResults/")); //dirname_ += iConfig.getParameter<std::string>("@module_label"); if (dbe_ != 0 ) { dbe_->setCurrentFolder(dirname_); } doCombineRuns_ = iConfig.getUntrackedParameter<bool>("doCombineRuns", false); doVBTFMuon_ = iConfig.getUntrackedParameter<bool>("doVBTFMuon", true); processname_ = iConfig.getParameter<std::string>("processname"); fCustomBXPath = iConfig.getUntrackedParameter<std::string>("customBXPath", std::string("HLT_MinBiasBSC")); referenceBX_ = iConfig.getUntrackedParameter<unsigned int>("referenceBX",51); Nbx_ = iConfig.getUntrackedParameter<unsigned int>("Nbx",3564); // plotting paramters ptMin_ = iConfig.getUntrackedParameter<double>("ptMin",0.); ptMax_ = iConfig.getUntrackedParameter<double>("ptMax",1000.); nBins_ = iConfig.getUntrackedParameter<unsigned int>("Nbins",20); nBins2D_ = iConfig.getUntrackedParameter<unsigned int>("Nbins2D",40); dRMax_ = iConfig.getUntrackedParameter<double>("dRMax",1.0); dRMaxElectronMuon_ = iConfig.getUntrackedParameter<double>("dRMaxElectronMuon",0.3); nBinsDR_ = iConfig.getUntrackedParameter<unsigned int>("NbinsDR",10); nBinsOneOverEt_ = iConfig.getUntrackedParameter<unsigned int>("NbinsOneOverEt",10000); nLS_ = iConfig.getUntrackedParameter<unsigned int>("NLuminositySegments",10); LSsize_ = iConfig.getUntrackedParameter<double>("LuminositySegmentSize",23); thresholdFactor_ = iConfig.getUntrackedParameter<double>("thresholdFactor",1.0); plotAll_ = iConfig.getUntrackedParameter<bool>("plotAll", 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++) { custompathnamepairs_.push_back( make_pair( pathconf->getParameter<std::string>("pathname"), pathconf->getParameter<std::string>("denompathname") ) ); } if (hltPaths_.size() > 0) { // book a histogram of scalers scalersSelect = dbe_->book1D("selectedScalers","Selected Scalers", hltPaths_.size(), 0.0, (double)hltPaths_.size()); } triggerSummaryLabel_ = iConfig.getParameter<edm::InputTag>("triggerSummaryLabel"); triggerResultsLabel_ = iConfig.getParameter<edm::InputTag>("triggerResultsLabel"); /* muonRecoCollectionName_ = iConfig.getUntrackedParameter("muonRecoCollectionName", std::string("muons")); electronEtaMax_ = iConfig.getUntrackedParameter<double>("electronEtaMax",2.5); electronEtMin_ = iConfig.getUntrackedParameter<double>("electronEtMin",3.0); electronDRMatch_ =iConfig.getUntrackedParameter<double>("electronDRMatch",0.3); electronL1DRMatch_ =iConfig.getUntrackedParameter<double>("electronL1DRMatch",0.3); muonEtaMax_ = iConfig.getUntrackedParameter<double>("muonEtaMax",2.1); muonEtMin_ = iConfig.getUntrackedParameter<double>("muonEtMin",0.0); muonDRMatch_ =iConfig.getUntrackedParameter<double>("muonDRMatch",0.3); muonL1DRMatch_ =iConfig.getUntrackedParameter<double>("muonL1DRMatch",0.3); tauEtaMax_ = iConfig.getUntrackedParameter<double>("tauEtaMax",2.5); tauEtMin_ = iConfig.getUntrackedParameter<double>("tauEtMin",3.0); tauDRMatch_ =iConfig.getUntrackedParameter<double>("tauDRMatch",0.3); tauL1DRMatch_ =iConfig.getUntrackedParameter<double>("tauL1DRMatch",0.5); jetEtaMax_ = iConfig.getUntrackedParameter<double>("jetEtaMax",5.0); jetEtMin_ = iConfig.getUntrackedParameter<double>("jetEtMin",10.0); jetDRMatch_ =iConfig.getUntrackedParameter<double>("jetDRMatch",0.3); jetL1DRMatch_ =iConfig.getUntrackedParameter<double>("jetL1DRMatch",0.5); bjetEtaMax_ = iConfig.getUntrackedParameter<double>("bjetEtaMax",2.5); bjetEtMin_ = iConfig.getUntrackedParameter<double>("bjetEtMin",10.0); bjetDRMatch_ =iConfig.getUntrackedParameter<double>("bjetDRMatch",0.3); bjetL1DRMatch_ =iConfig.getUntrackedParameter<double>("bjetL1DRMatch",0.3); photonEtaMax_ = iConfig.getUntrackedParameter<double>("photonEtaMax",2.5); photonEtMin_ = iConfig.getUntrackedParameter<double>("photonEtMin",3.0); photonDRMatch_ =iConfig.getUntrackedParameter<double>("photonDRMatch",0.3); photonL1DRMatch_ =iConfig.getUntrackedParameter<double>("photonL1DRMatch",0.3); trackEtaMax_ = iConfig.getUntrackedParameter<double>("trackEtaMax",2.5); trackEtMin_ = iConfig.getUntrackedParameter<double>("trackEtMin",3.0); trackDRMatch_ =iConfig.getUntrackedParameter<double>("trackDRMatch",0.3); trackL1DRMatch_ =iConfig.getUntrackedParameter<double>("trackL1DRMatch",0.3); metEtaMax_ = iConfig.getUntrackedParameter<double>("metEtaMax",5); metMin_ = iConfig.getUntrackedParameter<double>("metMin",10.0); metDRMatch_ =iConfig.getUntrackedParameter<double>("metDRMatch",0.5); metL1DRMatch_ =iConfig.getUntrackedParameter<double>("metL1DRMatch",0.5); htEtaMax_ = iConfig.getUntrackedParameter<double>("htEtaMax",5); htMin_ = iConfig.getUntrackedParameter<double>("htMin",10.0); htDRMatch_ =iConfig.getUntrackedParameter<double>("htDRMatch",0.5); htL1DRMatch_ =iConfig.getUntrackedParameter<double>("htL1DRMatch",0.5); */ sumEtMin_ = iConfig.getUntrackedParameter<double>("sumEtMin",10.0); // Muon quality cuts dxyCut_ = iConfig.getUntrackedParameter<double>("DxyCut", 0.2); // dxy < 0.2 cm normalizedChi2Cut_ = iConfig.getUntrackedParameter<double>("NormalizedChi2Cut", 10.); // chi2/ndof (of global fit) <10.0 trackerHitsCut_ = iConfig.getUntrackedParameter<int>("TrackerHitsCut", 11); // Tracker Hits >10 pixelHitsCut_ = iConfig.getUntrackedParameter<int>("PixelHitsCut", 1); // Pixel Hits >0 muonHitsCut_ = iConfig.getUntrackedParameter<int>("MuonHitsCut", 1); // Valid Muon Hits >0 isAlsoTrackerMuon_ = iConfig.getUntrackedParameter<bool>("IsAlsoTrackerMuon", true); nMatchesCut_ = iConfig.getUntrackedParameter<int>("NMatchesCut", 2); // At least 2 Chambers with matches specialPaths_ = iConfig.getParameter<std::vector<std::string > >("SpecialPaths"); testPathsFolder_ = iConfig.getUntrackedParameter ("testPathsFolder",std::string("HLT/TrigResults/testPaths/")); pathsSummaryFolder_ = iConfig.getUntrackedParameter ("pathsSummaryFolder",std::string("HLT/TrigResults/PathsSummary/HLT Counts/")); pathsSummaryStreamsFolder_ = iConfig.getUntrackedParameter ("pathsSummaryFolder",std::string("HLT/TrigResults/PathsSummary/")); //pathsSummaryStreamsFolder_ = iConfig.getUntrackedParameter ("pathsSummaryFolder",std::string("HLT/TrigResults/PathsSummary/Streams/")); pathsSummaryHLTCorrelationsFolder_ = iConfig.getUntrackedParameter ("hltCorrelationsFolder",std::string("HLT/TrigResults/PathsSummary/HLT Correlations/")); pathsSummaryFilterCountsFolder_ = iConfig.getUntrackedParameter ("filterCountsFolder",std::string("HLT/TrigResults/PathsSummary/Filters Counts/")); pathsSummaryHLTPathsPerLSFolder_ = iConfig.getUntrackedParameter ("individualPathsPerLSFolder",std::string("HLT/TrigResults/PathsSummary/HLT LS/")); pathsIndividualHLTPathsPerLSFolder_ = iConfig.getUntrackedParameter ("individualPathsPerLSFolder",std::string("HLT/TrigResults/PathsSummary/HLT LS/Paths/")); pathsSummaryHLTPathsPerBXFolder_ = iConfig.getUntrackedParameter ("individualPathsPerBXFolder",std::string("HLT/TrigResults/PathsSummary/HLT BX/")); // mask off some of the paths so that they don't appear in the plots maskedPaths_ = iConfig.getParameter<std::vector<std::string > >("MaskedPaths"); referenceTrigInput_ = iConfig.getParameter<std::string> ("ReferenceTrigger"); foundReferenceTrigger_ = false; //Robin testPaths_ = iConfig.getParameter<std::vector<std::string > >("testPaths"); fLumiFlag = true; ME_HLTAll_LS = NULL; ME_HLT_BX = NULL; ME_HLT_CUSTOM_BX = NULL; //jetID = new reco::helper::JetIDHelper(iConfig.getParameter<ParameterSet>("JetIDParams")); recHitsEBTag_ = iConfig.getUntrackedParameter<edm::InputTag>("RecHitsEBTag",edm::InputTag("reducedEcalRecHitsEB")); recHitsEETag_ = iConfig.getUntrackedParameter<edm::InputTag>("RecHitsEETag",edm::InputTag("reducedEcalRecHitsEE")); jmsDebug = false; jmsFakeZBCounts = false; found_zbIndex = false; if (jmsDebug ) std::cout << "Printing extra info " << std::endl; }
TrigResRateMon::~TrigResRateMon | ( | ) |
Definition at line 192 of file TrigResRateMon.cc.
{ // do anything here that needs to be done at desctruction time // (e.g. close files, deallocate resources etc.) }
void TrigResRateMon::addLumiToAverage | ( | double | lumi | ) | [private] |
Definition at line 1635 of file TrigResRateMon.cc.
References averageInstLumi, and fjr2json::lumi.
Referenced by analyze().
{ if (averageInstLumi == 0) { averageInstLumi = lumi; } else { averageInstLumi = (averageInstLumi + lumi) / 2; } }
void TrigResRateMon::analyze | ( | const edm::Event & | iEvent, |
const edm::EventSetup & | iSetup | ||
) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 207 of file TrigResRateMon.cc.
References addLumiToAverage(), gather_cfg::cout, funct::false, fillCountsPerPath(), fillHltMatrix(), edm::Event::getByLabel(), edm::InputTag::instance(), jmsDebug, edm::InputTag::label(), LogDebug, nev_, nPass_, nStream_, passAny, edm::Event::triggerNames(), patRefSel_triggerSelection_cff::triggerResults, triggerResults_, and triggerResultsLabel_.
{ //if(! fLumiFlag ) return; using namespace edm; using namespace trigger; ++nev_; LogDebug("TrigResRateMon")<< " analyze...." ; 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()) { edm::LogInfo("TrigResRateMon") << "TriggerResults not found, " "skipping event"; return; } } triggerResults_ = triggerResults; const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults); //Robin--- nStream_++; passAny = false ; // Find which index is zero bias for this run // // if(!found_zbIndex){ // // // set default to something that is out-of-bounds // // zbIndex = triggerNames.size() +2; // // for (unsigned int i = 0; i < triggerNames.size(); i ++){ // // std::string thisName = triggerNames.triggerName(i); // // TString checkName(thisName.c_str()); // // if (checkName.Contains("HLT_ZeroBias_v")){ // // zbIndex = i; // // found_zbIndex = true; // // if(jmsDebug) std::cout << "Found the ZeroBias index!!!!!!!! It is " << zbIndex <<std::endl; // // } // // } // // } // int bx = iEvent.bunchCrossing(); /* // Fill HLTPassed_Correlation Matrix bin (i,j) = (Any,Any) // -------------------------------------------------------- int anyBinNumber = ME_HLTPassPass->getTH2F()->GetXaxis()->FindBin("HLT_Any"); // any triger accepted if(triggerResults->accept()){ ME_HLTPassPass->Fill(anyBinNumber-1,anyBinNumber-1);//binNumber1 = 0 = first filter } */ fillHltMatrix(triggerNames, iEvent, iSetup); // // Experimental Testing Area // Try to get the lumi information for these events // edm::Handle<LumiScalersCollection> lumiScalers; bool lumiHandleOK = iEvent.getByLabel(InputTag("hltScalersRawToDigi","",""), lumiScalers); if (jmsDebug) std::cout << "Tried to get lumi handle result = " << lumiHandleOK << std::endl; if (lumiHandleOK) { if (jmsDebug) std::cout << "LumiScalers size is: " << lumiScalers->size() << std::endl; if (lumiScalers->size()) { LumiScalersCollection::const_iterator it3 = lumiScalers->begin(); unsigned int lumisection = it3->sectionNumber(); if(lumisection){ if (jmsDebug) std::cout << "Instanteous Lumi is " << it3->instantLumi() << std::endl; if (jmsDebug) std::cout << "Instanteous Lumi Error is " <<it3->instantLumiErr() << std::endl; if (jmsDebug) std::cout << "Lumi Fill is " <<it3->lumiFill() << std::endl; if (jmsDebug) std::cout << "Lumi Fill is " <<it3->lumiRun() << std::endl; if (jmsDebug) std::cout << "Live Lumi Fill is " <<it3->liveLumiFill() << std::endl; if (jmsDebug) std::cout << "Live Lumi Run is " <<it3->liveLumiRun() << std::endl; addLumiToAverage(it3->instantLumi()); } // end }// end if lumi scalers exist }// end if lumi handle ok fillCountsPerPath(iEvent, iSetup); if (passAny) nPass_ ++ ; return; }
void TrigResRateMon::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 320 of file TrigResRateMon.cc.
References DQMStore::book1D(), dirname_, MonitorElement::Fill(), nev_, cppFunctionSkipper::operator, and DQMStore::setCurrentFolder().
{ nev_ = 0; DQMStore *dbe = 0; dbe = Service<DQMStore>().operator->(); if (dbe) { dbe->setCurrentFolder(dirname_); } MonitorElement* reportSummaryME = dbe->book1D("reportSummaryMap","report Summary Map",2,0,2); if(reportSummaryME) reportSummaryME->Fill(1); }
void TrigResRateMon::beginLuminosityBlock | ( | const edm::LuminosityBlock & | lumiSeg, |
const edm::EventSetup & | c | ||
) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 1866 of file TrigResRateMon.cc.
References clearCountsPerPath(), clearLumiAverage(), gather_cfg::cout, jmsDebug, nPass_, and nStream_.
{ //int lumi = int(lumiSeg.id().luminosityBlock()); //if(lumi < 74 || lumi > 77) fLumiFlag = false; //else fLumiFlag = true; if (jmsDebug) std::cout << "Inside begin lumi block" << std::endl; clearCountsPerPath(); clearLumiAverage(); nStream_ = 0 ; nPass_ = 0 ; }
void TrigResRateMon::beginRun | ( | const edm::Run & | run, |
const edm::EventSetup & | c | ||
) | [private, virtual] |
add dataset name and thier triggers to the list
Reimplemented from edm::EDAnalyzer.
Definition at line 347 of file TrigResRateMon.cc.
References averageInstLumi3LS, DQMStore::book1D(), bookCountsPerPath(), bookTestHisto(), clearLumiAverage(), TrigResRateMon::DatasetInfo::countsPerPathME_Name, gather_cfg::cout, HLTConfigProvider::datasetContent(), TrigResRateMon::DatasetInfo::datasetName, makeLayoutFileForGui::datasetNames, dbe_, dirname_, doCombineRuns_, fGroupL1TempCountPair, fGroupName, fGroupNamePathsPair, fGroupTempCountPair, findReferenceTriggerIndex(), first, fIsSetup, fPathBxTempCountPair, fPathTempCountPair, g, hltConfig_, hltPathsDiagonal_, i, edm::RunBase::id(), HLTConfigProvider::init(), jmsDebug, LogDebug, maskedPaths_, meAverageLumiPerLS, meCountsDroppedPerLS, meCountsPassPerLS, meCountsStreamPerLS, meDiagnostic, meXsecStreamPerLS, nLS_, cppFunctionSkipper::operator, EgammaValidation_cff::pathName, pathsSummaryFolder_, primaryDataSetInformation, TrigResRateMon::DatasetInfo::printMaskedPaths(), processname_, TrigResRateMon::DatasetInfo::ratePerLSME_Name, TrigResRateMon::DatasetInfo::rawCountsPerPathME_Name, rawCountsPerPD, TrigResRateMon::DatasetInfo::scaledXsecPerPathME_Name, edm::second(), DQMStore::setCurrentFolder(), TrigResRateMon::DatasetInfo::setMaskedPaths(), TrigResRateMon::DatasetInfo::setPaths(), setupHltMatrix(), HLTConfigProvider::streamContent(), AlCaHLTBitMon_QueryRunRegistry::string, TotalDroppedCounts, findQualityFiles::v, and TrigResRateMon::DatasetInfo::xsecPerPathME_Name.
{ LogDebug("TrigResRateMon") << "beginRun, run " << run.id(); if(fIsSetup) return; // HLT config does not change within runs! bool changed=false; if (!hltConfig_.init(run, c, processname_, changed)) { processname_ = "FU"; if (!hltConfig_.init(run, c, processname_, changed)){ if (jmsDebug) std::cout << "HLTConfigProvider failed to initialize."; } else { if (jmsDebug) std::cout << "Initialized HLTConfigProvider with name FU " << std::endl; } // check if trigger name in (new) config // cout << "Available TriggerNames are: " << endl; // hltConfig_.dump("Triggers"); } else { if (jmsDebug) std::cout << "Initialized HLTConfigProvider with name HLT " << std::endl; } if (1) { DQMStore *dbe = 0; dbe = Service<DQMStore>().operator->(); if (dbe) { dbe->setCurrentFolder(dirname_); } // const unsigned int n(hltConfig_.size()); TotalDroppedCounts = 0; //Robin-------Diagnostic plots-------- meDiagnostic = dbe->book1D("DroppedCounts Diagnose", "LSs vs Status;Status;LSs", 3, -0.5,2.5 ); meCountsDroppedPerLS = dbe->book1D("CountsDroppedVsLS", "Counts vs LumiSec;LS;Dropped Stream Counts", nLS_, 0.5, nLS_+0.5); meCountsPassPerLS = dbe->book1D("CountsPassVsLS", "Counts vs LumiSec;LS;Passed Stream Counts", nLS_, 0.5, nLS_+0.5); meCountsStreamPerLS = dbe->book1D("CountsStreamVsLS", "Counts vs LumiSec;LS;Stream Counts", nLS_, 0.5, nLS_+0.5); meXsecStreamPerLS = dbe->book1D("XsecStreamVsLS", "Xsec vs LumiSec;LS;Stream Xsec", nLS_, 0.5, nLS_+0.5); // meXsecPerLS = dbe->book1D("XsecVsLS", "Xsec vs LumiSec;LS;Xsec", nLS_, 0.5, nLS_+0.5); // meXsec = dbe->book1D("Xsec", "histo for Xsec ", 20, 0.01, 0.06); // // meXsecPerIL = dbe->book2D("XsecVsIL", "Xsec vs Inst Lumi;#mub^{-1}*s^{-1}; Xsec", 200, 700, 900, 100, 0.01, 0.1); // TProfile tempProfile("XsecVsIL", "Xsec vs Inst Lumi;#mub^{-1}*s^{-1}; Xsec", 40, 600, 3000); // meXsecPerIL = dbe_->bookProfile("XsecVsIL", &tempProfile); bookTestHisto(); //Robin dbe->setCurrentFolder(dirname_); //Robin // JMS fill the counts per path bookCountsPerPath(); clearLumiAverage(); averageInstLumi3LS = 0 ; findReferenceTriggerIndex(); meAverageLumiPerLS = dbe->book1D("InstLumiVsLS", "Instantaneous Luminosity vs LumiSec;LS;#mub^{-1}*s^{-1}", nLS_, 0.5, nLS_+0.5); // // if (plotAll_){ // // for (unsigned int j=0; j!=n; ++j) { // // std::string pathname = hltConfig_.triggerName(j); // // string l1pathname = getL1ConditionModuleName(pathname); // // int l1ModuleIndex = hltConfig_.moduleIndex(pathname, l1pathname); // // int objectType = getTriggerTypeParsePathName(pathname); // // for (unsigned int i=0; i!=n; ++i) { // // std::string denompathname = hltConfig_.triggerName(i); // // int denomobjectType = getTriggerTypeParsePathName(denompathname); // // std::string filtername("dummy"); // // float ptMin = 0.0; // // float ptMax = 100.0; // // if (plotAll_ && denomobjectType == objectType && objectType != 0) { // // int hltThreshold = getThresholdFromName(pathname); // // int l1Threshold = getThresholdFromName(l1pathname); // // hltPaths_.push_back(PathInfo(denompathname, pathname, l1pathname, l1ModuleIndex, filtername, processname_, objectType, ptMin, ptMax, hltThreshold, l1Threshold)); // // } // // } // // } // // } // end if plotAll // // else { // plot all diagonal combinations plus any other specified pairs // // for (unsigned int i=0; i!=n; ++i) { // // std::string denompathname = ""; // // std::string pathname = hltConfig_.triggerName(i); // // //parse pathname to guess object type // // int objectType = getTriggerTypeParsePathName(pathname); // // string l1pathname = getL1ConditionModuleName(pathname); // // int l1ModuleIndex = hltConfig_.moduleIndex(pathname, l1pathname); // // std::string filtername("dummy"); // // float ptMin = 0.0; // // float ptMax = 100.0; // // if (objectType == trigger::TriggerPhoton) ptMax = 100.0; // // if (objectType == trigger::TriggerElectron) ptMax = 100.0; // // if (objectType == trigger::TriggerMuon) ptMax = 150.0; // // if (objectType == trigger::TriggerTau) ptMax = 100.0; // // if (objectType == trigger::TriggerJet) ptMax = 300.0; // // if (objectType == trigger::TriggerBJet) ptMax = 300.0; // // if (objectType == trigger::TriggerMET) ptMax = 300.0; // // if (objectType == trigger::TriggerTET) ptMax = 300.0; // // if (objectType == trigger::TriggerTrack) ptMax = 100.0; // // // keep track of all paths, except for FinalPath // // if (objectType != -1 && pathname.find("FinalPath") == std::string::npos){ // // int hltThreshold = getThresholdFromName(pathname); // // int l1Threshold = getThresholdFromName(l1pathname); // // hltPaths_.push_back(PathInfo(denompathname, pathname, l1pathname, l1ModuleIndex, filtername, processname_, objectType, ptMin, ptMax, hltThreshold, l1Threshold)); // // hltPathsDiagonal_.push_back(PathInfo(denompathname, pathname, l1pathname, l1ModuleIndex, filtername, processname_, objectType, ptMin, ptMax, hltThreshold, l1Threshold)); // // } // // } // end for i // now loop over denom/num path pairs specified in cfg, // recording the off-diagonal ones // // for (std::vector<std::pair<std::string, std::string> >::iterator custompathnamepair = custompathnamepairs_.begin(); custompathnamepair != custompathnamepairs_.end(); ++custompathnamepair) { // // std::string numpathname = custompathnamepair->first; // // std::string denompathname = custompathnamepair->second; // // if (numpathname != denompathname) { // // // check that denominator exists // // bool founddenominator = false; // // for (unsigned int k=0; k!=n; ++k) { // // string n_pathname = hltConfig_.triggerName(k); // // if (n_pathname.find(denompathname) != std::string::npos) { // // LogDebug("TrigResRateMon") << "denompathname is selected to be = " << n_pathname << endl;; // // founddenominator = true; // // break; // // } // // } // // if (!founddenominator) { // // edm::LogInfo("TrigResRateMon") << "denompathname not found, go to the next pair numearator-denominator" << endl; // // // go to the next pair // // continue; // // } // // // check that numerator exists // // bool foundnumerator = false; // // for (unsigned int j=0; j!=n; ++j) { // // string pathname = hltConfig_.triggerName(j); // // LogDebug("TrigResRateMon") << "check if path " << pathname << " is numpathname = " << numpathname << endl; // // if (hltConfig_.triggerName(j).find(numpathname)!= std::string::npos) { // // LogDebug("TrigResRateMon") << "pathname is selected to be = " << denompathname << endl;; // // foundnumerator = true; // // } // // if (!foundnumerator) { // // edm::LogInfo("TrigResRateMon") << "pathname not found, ignoring " << pathname; // // continue; // // } // // string l1pathname = getL1ConditionModuleName(pathname); // // int l1ModuleIndex = hltConfig_.moduleIndex(pathname, l1pathname); // // int objectType = getTriggerTypeParsePathName(pathname); // // std::string filtername("dummy"); // // float ptMin = 0.0; // // float ptMax = 100.0; // // if (objectType == trigger::TriggerPhoton) ptMax = 100.0; // // if (objectType == trigger::TriggerElectron) ptMax = 100.0; // // if (objectType == trigger::TriggerMuon) ptMax = 150.0; // // if (objectType == trigger::TriggerTau) ptMax = 100.0; // // if (objectType == trigger::TriggerJet) ptMax = 300.0; // // if (objectType == trigger::TriggerBJet) ptMax = 300.0; // // if (objectType == trigger::TriggerMET) ptMax = 300.0; // // if (objectType == trigger::TriggerTET) ptMax = 300.0; // // if (objectType == trigger::TriggerTrack) ptMax = 100.0; // // // monitor regardless of the objectType of the path // // if (objectType != 0) { // // int hltThreshold = getThresholdFromName(pathname); // // int l1Threshold = getThresholdFromName(l1pathname); // // hltPaths_.push_back(PathInfo(denompathname, pathname, l1pathname, l1ModuleIndex, filtername, processname_, objectType, ptMin, ptMax, hltThreshold, l1Threshold)); // // } // // } // end for j, loop over paths // // } // end if not same num and denominator // // } // // end for pair // } // end else /* vector<string> muonPaths; vector<string> egammaPaths; vector<string> tauPaths; vector<string> jetmetPaths; vector<string> restPaths; */ vector<string> allPaths; // fill vectors of Muon, Egamma, JetMET, Rest, and Special paths for(PathInfoCollection::iterator v = hltPathsDiagonal_.begin(); v!= hltPathsDiagonal_.end(); ++v ) { std::string pathName = v->getPath(); //int objectType = v->getObjectType(); vector<int> tempCount(5,0); fPathTempCountPair.push_back(make_pair(pathName,0)); fPathBxTempCountPair.push_back(make_pair(pathName,tempCount)); allPaths.push_back(pathName); /* switch (objectType) { case trigger::TriggerMuon : muonPaths.push_back(pathName); break; case trigger::TriggerElectron : case trigger::TriggerPhoton : egammaPaths.push_back(pathName); break; case trigger::TriggerTau : tauPaths.push_back(pathName); break; case trigger::TriggerJet : case trigger::TriggerMET : jetmetPaths.push_back(pathName); break; default: restPaths.push_back(pathName); } */ } fPathTempCountPair.push_back(make_pair("HLT_Any",0)); fGroupName.push_back("AllSelectedPaths"); /* fGroupName.push_back("Muon"); fGroupName.push_back("Egamma"); fGroupName.push_back("Tau"); fGroupName.push_back("JetMET"); fGroupName.push_back("Rest"); fGroupName.push_back("Special"); */ for(unsigned int g=0; g<fGroupName.size(); g++) { //fGroupTempCountPair.push_back(make_pair(fGroupName[g],0)); //fGroupL1TempCountPair.push_back(make_pair(fGroupName[g],0)); } dbe_->setCurrentFolder(pathsSummaryFolder_.c_str()); // fGroupNamePathsPair.push_back(make_pair("AllSelectedPaths",allPaths)); /* fGroupNamePathsPair.push_back(make_pair("Muon",muonPaths)); fGroupNamePathsPair.push_back(make_pair("Egamma",egammaPaths)); fGroupNamePathsPair.push_back(make_pair("Tau",tauPaths)); fGroupNamePathsPair.push_back(make_pair("JetMET",jetmetPaths)); fGroupNamePathsPair.push_back(make_pair("Rest",restPaths)); fGroupNamePathsPair.push_back(make_pair("Special",specialPaths_)); */ //vector<string> datasetNames = hltConfig_.datasetNames() ; vector<string> datasetNames = hltConfig_.streamContent("A") ; for (unsigned int i=0;i<datasetNames.size();i++) { vector<string> datasetPaths = hltConfig_.datasetContent(datasetNames[i]); fGroupNamePathsPair.push_back(make_pair(datasetNames[i],datasetPaths)); DatasetInfo tempDS; tempDS.datasetName = datasetNames[i]; tempDS.setPaths(datasetPaths); tempDS.countsPerPathME_Name = pathsSummaryFolder_ + "HLT_" + datasetNames[i] + "_Pass_Any"; tempDS.xsecPerPathME_Name = pathsSummaryFolder_ + "HLT_" + datasetNames[i] + "_Xsec"; tempDS.rawCountsPerPathME_Name = pathsSummaryFolder_ + "HLT_" + datasetNames[i] + "_RawCounts"; tempDS.scaledXsecPerPathME_Name = pathsSummaryFolder_ + "HLT_" + datasetNames[i] + "_XsecScaled"; tempDS.ratePerLSME_Name = pathsSummaryFolder_ + "HLT_" + datasetNames[i] + "_Rate"; //Robin tempDS.setMaskedPaths(maskedPaths_); if (jmsDebug) tempDS.printMaskedPaths(); primaryDataSetInformation.push_back(tempDS); rawCountsPerPD.push_back(0); //Robin } // push stream A and its PDs fGroupNamePathsPair.push_back(make_pair("A",datasetNames)); for (unsigned int g=0;g<fGroupNamePathsPair.size();g++) { fGroupTempCountPair.push_back(make_pair(fGroupNamePathsPair[g].first,0)); fGroupL1TempCountPair.push_back(make_pair(fGroupNamePathsPair[g].first,0)); setupHltMatrix(fGroupNamePathsPair[g].first,fGroupNamePathsPair[g].second); } /* // HLT matrices from Streams const std::vector<std::string> streamNames = hltConfig_.streamNames(); for (unsigned int s=0;s<streamNames.size();s++) { vector<string> hltConfig = streamDatasetNames_.streamContent(streamNames[s]) ; if(streamNames[s] == "A") setupStreamMatrix(streamNames[s],streamDatasetNames); } */ // // setupHltLsPlots(); // // setupHltBxPlots(); // for(PathInfoCollection::iterator v = hltPathsDiagonal_.begin(); v!= hltPathsDiagonal_.end(); ++v ) { // // ------------------------- // // // // Filters for each path // // // // ------------------------- // // get all modules in this HLT path // vector<string> moduleNames = hltConfig_.moduleLabels( v->getPath() ); // int numModule = 0; // string moduleName, moduleType; // unsigned int moduleIndex; // //print module name // vector<string>::const_iterator iDumpModName; // for (iDumpModName = moduleNames.begin();iDumpModName != moduleNames.end();iDumpModName++) { // moduleName = *iDumpModName; // moduleType = hltConfig_.moduleType(moduleName); // moduleIndex = hltConfig_.moduleIndex(v->getPath(), moduleName); // LogTrace ("TrigResRateMon") << "Module " << numModule // << " is called " << moduleName // << " , type = " << moduleType // << " , index = " << moduleIndex // << endl; // numModule++; // if((moduleType.find("Filter") != string::npos && moduleType.find("HLTTriggerTypeFilter") == string::npos ) || // (moduleType.find("Associator") != string::npos) || // (moduleType.find("HLTLevel1GTSeed") != string::npos) || // (moduleType.find("HLTGlobalSumsCaloMET") != string::npos) || // (moduleType.find("HLTPrescaler") != string::npos) ) { // //std::pair<std::string, int> filterIndexPair; // //filterIndexPair.first = moduleName; // //filterIndexPair.second = moduleIndex; // //v->filtersAndIndices.push_back(filterIndexPair); // v->filtersAndIndices.push_back(make_pair(moduleName,moduleIndex)); // } // }//end for modulesName // // dbe_->setCurrentFolder(pathsSummaryFilterCountsFolder_.c_str()); // // //int nbin_sub = 5; // // int nbin_sub = v->filtersAndIndices.size()+2; // // // count plots for subfilter // // MonitorElement* filters = dbe_->book1D("Filters_" + v->getPath(), // // "Filters_" + v->getPath(), // // nbin_sub+1, -0.5, 0.5+(double)nbin_sub); // // for(unsigned int filt = 0; filt < v->filtersAndIndices.size(); filt++){ // // filters->setBinLabel(filt+1, (v->filtersAndIndices[filt]).first); // // } // // book Count vs LS // dbe_->setCurrentFolder(pathsIndividualHLTPathsPerLSFolder_.c_str()); // MonitorElement* tempME = dbe_->book1D(v->getPath() + "_count_per_LS", // v->getPath() + " rate [Hz]", // nLS_, 0,nLS_); // tempME->setAxisTitle("Luminosity Section"); // //v->setFilterHistos(filters); // } // end for paths // now set up all of the histos for each path-denom // for(PathInfoCollection::iterator v = hltPaths_.begin(); v!= hltPaths_.end(); ++v ) { /* MonitorElement *NOn, *onEtOn, *onOneOverEtOn, *onEtavsonPhiOn=0; MonitorElement *NOff, *offEtOff, *offEtavsoffPhiOff=0; MonitorElement *NL1, *l1EtL1, *l1Etavsl1PhiL1=0; MonitorElement *NL1On, *l1EtL1On, *l1Etavsl1PhiL1On=0; MonitorElement *NL1Off, *offEtL1Off, *offEtavsoffPhiL1Off=0; MonitorElement *NOnOff, *offEtOnOff, *offEtavsoffPhiOnOff=0; MonitorElement *NL1OnUM, *l1EtL1OnUM, *l1Etavsl1PhiL1OnUM=0; MonitorElement *NL1OffUM, *offEtL1OffUM, *offEtavsoffPhiL1OffUM=0; MonitorElement *NOnOffUM, *offEtOnOffUM, *offEtavsoffPhiOnOffUM=0; MonitorElement *offDRL1Off, *offDROnOff, *l1DRL1On=0; */ // std::string labelname("dummy"); // labelname = v->getPath() + "_wrt_" + v->getDenomPath(); // std::string histoname(labelname+"_NOn"); // std::string title(labelname+" N online"); // double histEtaMax = 2.5; // if (v->getObjectType() == trigger::TriggerMuon || v->getObjectType() == trigger::TriggerL1Mu) { // histEtaMax = muonEtaMax_; // } // else if (v->getObjectType() == trigger::TriggerElectron || v->getObjectType() == trigger::TriggerL1NoIsoEG || v->getObjectType() == trigger::TriggerL1IsoEG ) // { // histEtaMax = electronEtaMax_; // } // else if (v->getObjectType() == trigger::TriggerTau || v->getObjectType() == trigger::TriggerL1TauJet ) // { // histEtaMax = tauEtaMax_; // } // else if (v->getObjectType() == trigger::TriggerJet || v->getObjectType() == trigger::TriggerL1CenJet || v->getObjectType() == trigger::TriggerL1ForJet ) // { // histEtaMax = jetEtaMax_; // } // else if (v->getObjectType() == trigger::TriggerBJet) // { // histEtaMax = bjetEtaMax_; // } // else if (v->getObjectType() == trigger::TriggerMET || v->getObjectType() == trigger::TriggerL1ETM ) // { // histEtaMax = metEtaMax_; // } // else if (v->getObjectType() == trigger::TriggerPhoton) // { // histEtaMax = photonEtaMax_; // } // else if (v->getObjectType() == trigger::TriggerTrack) // { // histEtaMax = trackEtaMax_; // } // TString pathfolder = dirname_ + TString("/FourVector/") + v->getPath(); /* dbe_->setCurrentFolder(pathfolder.Data()); NOn = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NOff"; title = labelname+" N Off"; NOff = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NL1"; title = labelname+" N L1"; NL1 = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NL1On"; title = labelname+" N L1On"; NL1On = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NL1Off"; title = labelname+" N L1Off"; NL1Off = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NOnOff"; title = labelname+" N OnOff"; NOnOff = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NL1OnUM"; title = labelname+" N L1OnUM"; NL1OnUM = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NL1OffUM"; title = labelname+" N L1OffUM"; NL1OffUM = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_NOnOffUM"; title = labelname+" N OnOffUM"; NOnOffUM = dbe->book1D(histoname.c_str(), title.c_str(),10, 0.5, 10.5); histoname = labelname+"_onEtOn"; title = labelname+" onE_t online"; onEtOn = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_onOneOverEtOn"; title = labelname+" 1 / onE_t online"; onOneOverEtOn = dbe->book1D(histoname.c_str(), title.c_str(),nBinsOneOverEt_, 0, 1); onOneOverEtOn->setAxisTitle("HLT 1/Et [1/GeV]"); histoname = labelname+"_offEtOff"; title = labelname+" offE_t offline"; offEtOff = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_l1EtL1"; title = labelname+" l1E_t L1"; l1EtL1 = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_onEtaonPhiOn"; title = labelname+" on#eta vs on#phi online"; onEtavsonPhiOn = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_offEtaoffPhiOff"; title = labelname+" off#eta vs off#phi offline"; offEtavsoffPhiOff = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_l1Etal1PhiL1"; title = labelname+" l1#eta vs l1#phi L1"; l1Etavsl1PhiL1 = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_l1EtL1On"; title = labelname+" l1E_t L1+online"; l1EtL1On = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_offEtL1Off"; title = labelname+" offE_t L1+offline"; offEtL1Off = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_offEtOnOff"; title = labelname+" offE_t online+offline"; offEtOnOff = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_l1Etal1PhiL1On"; title = labelname+" l1#eta vs l1#phi L1+online"; l1Etavsl1PhiL1On = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_offEtaoffPhiL1Off"; title = labelname+" off#eta vs off#phi L1+offline"; offEtavsoffPhiL1Off = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_offEtaoffPhiOnOff"; title = labelname+" off#eta vs off#phi online+offline"; offEtavsoffPhiOnOff = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_l1EtL1OnUM"; title = labelname+" l1E_t L1+onlineUM"; l1EtL1OnUM = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_offEtL1OffUM"; title = labelname+" offE_t L1+offlineUM"; offEtL1OffUM = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_offEtOnOffUM"; title = labelname+" offE_t online+offlineUM"; offEtOnOffUM = dbe->book1D(histoname.c_str(), title.c_str(),nBins_, v->getPtMin(), v->getPtMax()); histoname = labelname+"_l1Etal1PhiL1OnUM"; title = labelname+" l1#eta vs l1#phi L1+onlineUM"; l1Etavsl1PhiL1OnUM = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_offEtaoffPhiL1OffUM"; title = labelname+" off#eta vs off#phi L1+offlineUM"; offEtavsoffPhiL1OffUM = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_offEtaoffPhiOnOffUM"; title = labelname+" off#eta vs off#phi online+offlineUM"; offEtavsoffPhiOnOffUM = dbe->book2D(histoname.c_str(), title.c_str(), nBins2D_,-histEtaMax,histEtaMax, nBins2D_,-TMath::Pi(), TMath::Pi()); histoname = labelname+"_l1DRL1On"; title = labelname+" l1DR L1+online"; l1DRL1On = dbe->book1D(histoname.c_str(), title.c_str(),nBinsDR_, 0, dRMax_); histoname = labelname+"_offDRL1Off"; title = labelname+" offDR L1+offline"; offDRL1Off = dbe->book1D(histoname.c_str(), title.c_str(),nBinsDR_, 0, dRMax_); histoname = labelname+"_offDROnOff"; title = labelname+" offDR online+offline"; offDROnOff = dbe->book1D(histoname.c_str(), title.c_str(),nBinsDR_, 0, dRMax_); */ //v->setHistos( NOn, onEtOn, onOneOverEtOn, onEtavsonPhiOn, NOff, offEtOff, offEtavsoffPhiOff, NL1, l1EtL1, l1Etavsl1PhiL1, NL1On, l1EtL1On, l1Etavsl1PhiL1On, NL1Off, offEtL1Off, offEtavsoffPhiL1Off, NOnOff, offEtOnOff, offEtavsoffPhiOnOff, NL1OnUM, l1EtL1OnUM, l1Etavsl1PhiL1OnUM, NL1OffUM, offEtL1OffUM, offEtavsoffPhiL1OffUM, NOnOffUM, offEtOnOffUM, offEtavsoffPhiOnOffUM, offDRL1Off, offDROnOff, l1DRL1On ); // } // end for hltPath // HLT_Any // book Count vs LS // // dbe_->setCurrentFolder(pathsIndividualHLTPathsPerLSFolder_.c_str()); // // MonitorElement* tempME = dbe_->book1D("HLT_Any_count_per_LS", // // "HLT_Any rate [Hz]", // // nLS_, 0,nLS_); // // tempME->setAxisTitle("Luminosity Section"); } // end if(1) dummy if(doCombineRuns_) fIsSetup = true; return; }
void TrigResRateMon::bookCountsPerPath | ( | ) | [private] |
Definition at line 1544 of file TrigResRateMon.cc.
References finalCountsPerPath, hltConfig_, rawCountsPerPath, and HLTConfigProvider::size().
Referenced by beginRun().
{ for (unsigned iName = 0; iName < hltConfig_.size(); iName++) { rawCountsPerPath.push_back(0); finalCountsPerPath.push_back(0); //Robin } }
void TrigResRateMon::bookTestHisto | ( | ) | [private] |
Definition at line 1796 of file TrigResRateMon.cc.
References DQMStore::book1D(), DQMStore::bookProfile(), dbe_, MonitorElement::getTProfile(), i, meXsecPerTestPath, mergeVDriftHistosByStation::name, nLS_, timingPdfMaker::pathname, MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), testPaths_, testPathsFolder_, indexGen::title, v_ME_CountsPerLS, v_ME_Xsec, and v_ME_XsecPerLS.
Referenced by beginRun().
{ unsigned int npaths = testPaths_.size(); //pathsSummaryHLTPathsPerLSFolder_ = TString("HLT/TrigResults/PathsSummary/HLT LS/"); //dbe_->setCurrentFolder(pathsSummaryHLTPathsPerLSFolder_.c_str()); //dbe_->setCurrentFolder(pathsSummaryHLTPathsPerLSFolder_); dbe_->setCurrentFolder(testPathsFolder_); TProfile tempProfile("XsecVsTestPath", "Xsec vs Test Path;;Xsec #mub", npaths, 0.5, npaths+0.5); meXsecPerTestPath = dbe_->bookProfile("XsecVsTestPath", &tempProfile); for(unsigned int i = 0; i < npaths; i++){ const char* pname = testPaths_[i].c_str() ; TString pathname = testPaths_[i].c_str() ; meXsecPerTestPath->getTProfile()->GetXaxis()->SetBinLabel(i+1, pathname ); char name[200]; char title[200]; sprintf(name, "path_%s_XsecVsLS", pname); sprintf(title, "path_%s_XsecVsLS;LS;Xsec #mub", pname); MonitorElement* tempME = dbe_->book1D(name,title, nLS_, 0, nLS_ ); tempME->setAxisTitle("LS"); v_ME_XsecPerLS.push_back(tempME); char name2[200]; char title2[200]; sprintf(name2, "path_%s_countsVsLS", pname); sprintf(title2, "path_%s_countsVsLS;LS;Counts", pname); MonitorElement* tempME2 = dbe_->book1D(name2, title2, nLS_, 0, nLS_ ); tempME2->setAxisTitle("LS"); v_ME_CountsPerLS.push_back(tempME2); } MonitorElement* meXsec1 = dbe_->book1D("Xsec_HLT_IsoMu30_eta2p1", "HLT_IsoMu30_eta2p1 Xsec;Xsec #mub;number of LS", 10, 0.008, 0.012); MonitorElement* meXsec2 = dbe_->book1D("Xsec_HLT_Ele65_CaloIdVT_TrkIdT", "HLT_Ele65_CaloIdVT_TrkIdT Xsec;Xsec #mub;number of LS", 10, 0.002, 0.0025); MonitorElement* meXsec3 = dbe_->book1D("Xsec_HLT_MET200", "HLT_MET200 Xsec;Xsec #mub;number of LS", 10, 0.0004, 0.0008); MonitorElement* meXsec4 = dbe_->book1D("Xsec_HLT_Jet370", "HLT_Jet370 Xsec;Xsec #mub;number of LS", 10, 0.0006, 0.0008); MonitorElement* meXsec5 = dbe_->book1D("Xsec_HLT_HT600", "HLT_HT600 Xsec;Xsec #mub;number of LS", 10, 0.004, 0.005); MonitorElement* meXsec6 = dbe_->book1D("Xsec_HLT_Photon26_R9Id_Photon18_R9Id", "HLT_Photon26_R9Id_Photon18_R9Id Xsec;Xsec #mub;number of LS", 10, 0.002, 0.004); MonitorElement* meXsec7 = dbe_->book1D("Xsec_HLT_IsoMu15_eta2p1_LooseIsoPFTau20", "HLT_IsoMu15_eta2p1_LooseIsoPFTau20 Xsec;Xsec #mub;number of LS", 10, 0.0022, 0.003); MonitorElement* meXsec8 = dbe_->book1D("Xsec_HLT_PFMHT150", "HLT_PFMHT150 Xsec;Xsec #mub;number of LS", 10, 0.0005, 0.001); MonitorElement* meXsec9 = dbe_->book1D("Xsec_HLT_Photon90_CaloIdVL_IsoL", "HLT_Photon90_CaloIdVL_IsoL Xsec;Xsec #mub;number of LS", 10, 0.0015, 0.0025); v_ME_Xsec.push_back(meXsec1); v_ME_Xsec.push_back(meXsec2); v_ME_Xsec.push_back(meXsec3); v_ME_Xsec.push_back(meXsec4); v_ME_Xsec.push_back(meXsec5); v_ME_Xsec.push_back(meXsec6); v_ME_Xsec.push_back(meXsec7); v_ME_Xsec.push_back(meXsec8); v_ME_Xsec.push_back(meXsec9); }
void TrigResRateMon::clearCountsPerPath | ( | ) | [private] |
Definition at line 1607 of file TrigResRateMon.cc.
References primaryDataSetInformation, and referenceTrigCountsPS_.
Referenced by beginLuminosityBlock().
{ // for (unsigned iName = 0; iName < hltConfig_.size() ; iName++) { // rawCountsPerPath[iName] = 0; // finalCountsPerPath[iName] = 0; //Robin // } referenceTrigCountsPS_ = 0 ; // testTrigCountsPS_ = 0 ; //Robin for (std::vector<DatasetInfo>::iterator iDS = primaryDataSetInformation.begin(); iDS != primaryDataSetInformation.end(); iDS++) { iDS->clearCountsPerPath(); } }
void TrigResRateMon::clearLumiAverage | ( | ) | [private] |
Definition at line 1629 of file TrigResRateMon.cc.
References averageInstLumi.
Referenced by beginLuminosityBlock(), and beginRun().
{ averageInstLumi = 0; }
void TrigResRateMon::countHLTGroupBXHitsEndLumiBlock | ( | const int & | lumi | ) | [private] |
Definition at line 2034 of file TrigResRateMon.cc.
References b, plotTTrigFromHistos::binNumber(), fGroupNamePathsPair, fPathBxTempCountPair, MonitorElement::getTH2F(), j, gen::k, LogDebug, LogTrace, ME_HLT_BX, Nbx_, timingPdfMaker::pathname, referenceBX_, edm::second(), findQualityFiles::size, AlCaHLTBitMon_QueryRunRegistry::string, v_ME_Total_BX, and v_ME_Total_BX_Norm.
{ LogTrace("TrigResRateMon") << " countHLTGroupBXHitsEndLumiBlock() lumiSection number " << lumi << endl; if(! ME_HLT_BX) return; TH2F * hist_2d_bx = ME_HLT_BX->getTH2F(); for (std::vector<std::pair<std::string, vector<int> > >::iterator ip = fPathBxTempCountPair.begin(); ip != fPathBxTempCountPair.end(); ++ip) { // get the path and its previous count std::string pathname = ip->first; vector<int> prevCount = ip->second; // vector of 5 zeros vector<int> currCount (5,0); vector<int> diffCount (5,0); // get the current count of path up to now int pathBin = hist_2d_bx->GetYaxis()->FindBin(pathname.c_str()); if(pathBin > hist_2d_bx->GetNbinsY()) { LogTrace("TrigResRateMon") << " Cannot find the bin for path " << pathname << endl; continue; } for (unsigned int b =0;b<currCount.size();b++) { int bxOffset = b-2; int bunch = referenceBX_+bxOffset; if(bunch < 1) bunch += Nbx_ ; int bxBin = bunch +1; // add one to get the right bin currCount[b] = int(hist_2d_bx->GetBinContent(bxBin, pathBin)); // add one to get the right bin LogTrace("TrigResRateMon") << "currCount = " << currCount[b] << endl; // count due to prev lumi sec is a difference bw current and previous diffCount[b] = currCount[b] - prevCount[b]; LogTrace("TrigResRateMon") << " lumi = " << lumi << " path " << pathname << "bxOffset = " << bxOffset << " count = " << diffCount[b] << endl; } // end for bx b // set the counter in the pair to current count ip->second = currCount; // fill the 2D Group paths' BX count per LS, using currCount LogTrace("TrigResRateMon") << "Find " << pathname << endl; //check if the path is in this group //for (unsigned int j=0;j<fGroupNamePathsPair.size();j++) for (unsigned int j=0;j<v_ME_Total_BX.size();j++) { bool isMember = false; LogTrace("TrigResRateMon") << " ---- Group " << fGroupNamePathsPair[j].first << endl; // decide if pathname is member of this group for (unsigned int k = 0; k<(fGroupNamePathsPair[j].second).size();k++) { LogTrace("TrigResRateMon") << " comparing to " << fGroupNamePathsPair[j].second[k] << endl; if(fGroupNamePathsPair[j].second[k] == pathname) { isMember = true; break; } } // end for k if(!isMember) { LogTrace("TrigResRateMon") << "Could not find a group to which the path belongs, path = " << pathname << " group = " << fGroupNamePathsPair[j].first << endl; continue; } MonitorElement* ME_2d = v_ME_Total_BX[j]; if (! ME_2d) { LogDebug("TrigResRateMon") << " cannot find ME_2d for group " << fGroupNamePathsPair[j].first << endl; continue; } vector<int> updatedLumiCount(5,0); float entireBXWindowUpdatedLumiCount = 0; TH2F* hist_All = ME_2d->getTH2F(); for (unsigned int b = 0; b<diffCount.size();b++) { // find the bin int binNumber = b+1; // add one to get right bin // update the bin content (must do that since events don't ncessarily come in the order int currentLumiCount = int(hist_All->GetBinContent(lumi+1,binNumber)); updatedLumiCount[b] = currentLumiCount + diffCount[b]; hist_All->SetBinContent(lumi+1,binNumber,updatedLumiCount[b]); entireBXWindowUpdatedLumiCount += updatedLumiCount[b]; } // end for bx b MonitorElement* ME_2d_Norm = v_ME_Total_BX_Norm[j]; if (! ME_2d_Norm) { LogDebug("TrigResRateMon") << " cannot find ME_2d_Norm for group " << fGroupNamePathsPair[j].first << endl; continue; } TH2F* hist_All_Norm = ME_2d_Norm->getTH2F(); for (unsigned int b = 0; b<diffCount.size();b++) { // find the bin int binNumber = b+1; // add one to get right bin // update the bin content but normalized to the whole columb (BX windw +/- 2) hist_All_Norm->SetBinContent(lumi+1,binNumber,float(updatedLumiCount[b])/entireBXWindowUpdatedLumiCount); } // end for bx b } // end for group j } // end for ip }
void TrigResRateMon::countHLTGroupHitsEndLumiBlock | ( | const int & | lumi | ) | [private] |
Definition at line 2256 of file TrigResRateMon.cc.
References dbe_, fGroupNamePathsPair, fGroupTempCountPair, DQMStore::get(), MonitorElement::getTH1(), MonitorElement::getTH1F(), MonitorElement::getTH2F(), i, LogDebug, LogTrace, LSsize_, timingPdfMaker::pathname, pathsSummaryFolder_, and pathsSummaryHLTPathsPerLSFolder_.
{ LogTrace("TrigResRateMon") << " countHLTGroupHitsEndLumiBlock() lumiSection number " << lumi << endl; for(unsigned int i=0; i<fGroupNamePathsPair.size(); i++){ // get the count of path up to now string fullPathToME = pathsSummaryFolder_ + "HLT_" + fGroupNamePathsPair[i].first + "_Pass_Any"; MonitorElement* ME_1d = dbe_->get(fullPathToME); if(! ME_1d) { LogTrace("TrigResRateMon") << " could not find 1d matrix " << fullPathToME << endl; continue; } LogTrace("TrigResRateMon") << " Looking in histogram " << fullPathToME << endl; TH1F * hist_1d = ME_1d->getTH1F(); for (std::vector<std::pair<std::string, float> >::iterator ip = fGroupTempCountPair.begin(); ip != fGroupTempCountPair.end(); ++ip) { // get the path and its previous count string pathname = ip->first; float prevCount = ip->second; string binLabel = "Total "+pathname; LogTrace("TrigResRateMon") << " Looking for binLabel = " << binLabel << endl; // get the current count of path up to now int pathBin = hist_1d->GetXaxis()->FindBin(binLabel.c_str()); LogTrace("TrigResRateMon") << " pathBin = " << pathBin << " out of histogram total number of bins " << hist_1d->GetNbinsX() << endl; if(pathBin == -1) { binLabel = pathname; int alternativePathBin = hist_1d->GetXaxis()->FindBin(binLabel.c_str()); if(alternativePathBin == -1) { LogTrace("TrigResRateMon") << " Cannot find the bin for path " << pathname << endl; continue; } else { pathBin = alternativePathBin; } } float currCount = hist_1d->GetBinContent(pathBin)/LSsize_; // count due to prev lumi sec is a difference bw current and previous float diffCount = currCount - prevCount; LogTrace("TrigResRateMon") << " lumi = " << lumi << " path " << pathname << " diffCount " << diffCount << endl; // set the counter in the pair to current count ip->second = currCount; // fill the 1D and 2D gruop and 2d_Stream_A count per LS string fullPathToME_count = pathsSummaryHLTPathsPerLSFolder_ +"HLT_" + pathname + "_Total_LS"; MonitorElement* ME_1d = dbe_->get(fullPathToME_count); string fullPathToME_2D_count = pathsSummaryHLTPathsPerLSFolder_ +"HLT_" + pathname + "_LS"; MonitorElement* ME_2d = dbe_->get(fullPathToME_2D_count); string fullPathToME_Stream_A_2D_count = pathsSummaryHLTPathsPerLSFolder_ +"HLT_A_LS"; MonitorElement* ME_Stream_A_2d = dbe_->get(fullPathToME_Stream_A_2D_count); if ( ME_1d && ME_2d && ME_Stream_A_2d) { // update the bin content (must do that since events don't ncessarily come in the order float currentLumiCount = ME_1d->getTH1()->GetBinContent(lumi+1); float updatedLumiCount = currentLumiCount + diffCount; ME_1d->getTH1()->SetBinContent(lumi+1,updatedLumiCount); string groupBinLabel = "Total " + fGroupNamePathsPair[i].first; int groupBin = ME_2d->getTH2F()->GetYaxis()->FindBin(groupBinLabel.c_str()); if(groupBin != -1) ME_2d->getTH2F()->SetBinContent(lumi+1,groupBin,updatedLumiCount); // this is to deal with Stream A and bins with names of PDs groupBinLabel = fGroupNamePathsPair[i].first; groupBin = ME_Stream_A_2d->getTH2F()->GetYaxis()->FindBin(groupBinLabel.c_str()); if(groupBin != -1) ME_Stream_A_2d->getTH2F()->SetBinContent(lumi+1,groupBin,updatedLumiCount); } else { LogDebug("TrigResRateMon") << " cannot find ME " << fullPathToME_count << endl; } } // end for ip } // end for i }
void TrigResRateMon::countHLTGroupL1HitsEndLumiBlock | ( | const int & | lumi | ) | [private] |
Definition at line 2176 of file TrigResRateMon.cc.
References dbe_, fGroupL1TempCountPair, fGroupNamePathsPair, DQMStore::get(), MonitorElement::getTH1(), MonitorElement::getTH1F(), i, LogDebug, LogTrace, LSsize_, timingPdfMaker::pathname, pathsSummaryFolder_, and pathsSummaryHLTPathsPerLSFolder_.
{ LogTrace("TrigResRateMon") << " countHLTGroupL1HitsEndLumiBlock() lumiSection number " << lumi << endl; for(unsigned int i=0; i<fGroupNamePathsPair.size(); i++){ // get the count of path up to now string fullPathToME = pathsSummaryFolder_ + "HLT_" + fGroupNamePathsPair[i].first+ "_Pass_Any"; MonitorElement* ME_1d = dbe_->get(fullPathToME); if(! ME_1d) { LogTrace("TrigResRateMon") << " could not find 1d matrix " << fullPathToME << endl; continue; } LogTrace("TrigResRateMon") << " Looking in histogram " << fullPathToME << endl; TH1F * hist_1d = ME_1d->getTH1F(); for (std::vector<std::pair<std::string, float> >::iterator ip = fGroupL1TempCountPair.begin(); ip != fGroupL1TempCountPair.end(); ++ip) { // get the path and its previous count string pathname = ip->first; float prevCount = ip->second; string binLabel = "HLT_"+pathname+"_L1_Any"; LogTrace("TrigResRateMon") << " Looking for binLabel = " << binLabel << endl; // get the current count of path up to now int pathBin = hist_1d->GetXaxis()->FindBin(binLabel.c_str()); LogTrace("TrigResRateMon") << " pathBin = " << pathBin << " out of histogram total number of bins " << hist_1d->GetNbinsX() << endl; if(pathBin == -1) { LogTrace("TrigResRateMon") << " Cannot find the bin for path " << pathname << endl; continue; } float currCount = hist_1d->GetBinContent(pathBin)/LSsize_; // count due to prev lumi sec is a difference bw current and previous float diffCount = currCount - prevCount; LogTrace("TrigResRateMon") << " lumi = " << lumi << " path " << pathname << " count " << diffCount << endl; // set the counter in the pair to current count ip->second = currCount; // fill the 1D individual path count per LS string fullPathToME_count = pathsSummaryHLTPathsPerLSFolder_ +"HLT_" + pathname + "_L1_Total_LS"; MonitorElement* ME_1d = dbe_->get(fullPathToME_count); if ( ME_1d) { // update the bin content (must do that since events don't ncessarily come in the order float currentLumiCount = ME_1d->getTH1()->GetBinContent(lumi+1); float updatedLumiCount = currentLumiCount + diffCount; ME_1d->getTH1()->SetBinContent(lumi+1,updatedLumiCount); } else { LogDebug("TrigResRateMon") << " cannot find ME " << fullPathToME_count << endl; } } // end for ip } // end for i }
void TrigResRateMon::countHLTPathHitsEndLumiBlock | ( | const int & | lumi | ) | [private] |
Definition at line 2365 of file TrigResRateMon.cc.
References dbe_, fPathTempCountPair, DQMStore::get(), MonitorElement::getTH1(), MonitorElement::getTH2F(), i, LogDebug, LogTrace, LSsize_, ME_HLTAll_LS, timingPdfMaker::pathname, pathsIndividualHLTPathsPerLSFolder_, pathsSummaryFolder_, AlCaHLTBitMon_QueryRunRegistry::string, and v_ME_HLTAll_LS.
{ LogTrace("TrigResRateMon") << " countHLTPathHitsEndLumiBlock() lumiSection number " << lumi << endl; // get the count of path up to now string fullPathToME = pathsSummaryFolder_ + "HLT_AllSelectedPaths_PassPass"; MonitorElement* ME_2d = dbe_->get(fullPathToME); if(! ME_2d) { LogTrace("TrigResRateMon") << " could not fine 2d matrix " << fullPathToME << endl; return; } TH2F * hist_2d = ME_2d->getTH2F(); for (std::vector<std::pair<std::string, float> >::iterator ip = fPathTempCountPair.begin(); ip != fPathTempCountPair.end(); ++ip) { // get the path and its previous count std::string pathname = ip->first; float prevCount = ip->second; // get the current count of path up to now float pathBin = hist_2d->GetXaxis()->FindBin(pathname.c_str()); if(pathBin > hist_2d->GetNbinsX()) { LogTrace("TrigResRateMon") << " Cannot find the bin for path " << pathname << endl; continue; } float currCount = hist_2d->GetBinContent(pathBin, pathBin)/LSsize_; // count due to prev lumi sec is a difference bw current and previous float diffCount = currCount - prevCount; LogTrace("TrigResRateMon") << " lumi = " << lumi << " path " << pathname << " count " << diffCount << endl; // set the counter in the pair to current count ip->second = currCount; // fill the 2D All paths' count per LS if ( ME_HLTAll_LS) { TH2F* hist_All = ME_HLTAll_LS->getTH2F(); // find the bin int pathBinNumber = hist_All->GetYaxis()->FindBin(pathname.c_str()); // update the bin content (must do that since events don't ncessarily come in the order float currentLumiCount = hist_All->GetBinContent(lumi+1,pathBinNumber); float updatedLumiCount = currentLumiCount + diffCount; hist_All->SetBinContent(lumi+1,pathBinNumber,updatedLumiCount); } else { LogDebug("TrigResRateMon") << " cannot find ME_HLTAll_LS" << endl; } for (unsigned int i=0 ; i< v_ME_HLTAll_LS.size(); i++) { MonitorElement* tempME = v_ME_HLTAll_LS[i]; if ( tempME ) { TH2F* hist_All = tempME->getTH2F(); // find the bin int pathBinNumber = hist_All->GetYaxis()->FindBin(pathname.c_str()); // update the bin content (must do that since events don't ncessarily come in the order float currentLumiCount = hist_All->GetBinContent(lumi+1,pathBinNumber); float updatedLumiCount = currentLumiCount + diffCount; hist_All->SetBinContent(lumi+1,pathBinNumber,updatedLumiCount); } else { LogDebug("TrigResRateMon") << " cannot find tempME " << endl; } } // fill the 1D individual path count per LS string fullPathToME_count = pathsIndividualHLTPathsPerLSFolder_ + pathname + "_count_per_LS"; MonitorElement* ME_1d = dbe_->get(fullPathToME_count); if ( ME_1d) { // update the bin content (must do that since events don't ncessarily come in the order float currentLumiCount = ME_1d->getTH1()->GetBinContent(lumi+1); float updatedLumiCount = currentLumiCount + diffCount; ME_1d->getTH1()->SetBinContent(lumi+1,updatedLumiCount); } else { LogDebug("TrigResRateMon") << " cannot find ME " << fullPathToME_count << endl; } } // end for ip }
void TrigResRateMon::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 339 of file TrigResRateMon.cc.
References nev_.
void TrigResRateMon::endLuminosityBlock | ( | const edm::LuminosityBlock & | lumiSeg, |
const edm::EventSetup & | c | ||
) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 1881 of file TrigResRateMon.cc.
References averageInstLumi, gather_cfg::cout, dbe_, MonitorElement::DQM_KIND_INT, filltestHisto(), fillXsecPerDataset(), DQMStore::get(), MonitorElement::getIntValue(), MonitorElement::getTH1F(), edm::LuminosityBlockBase::id(), jmsDebug, MonitorElement::kind(), LogTrace, LSsize_, fjr2json::lumi, edm::LuminosityBlockID::luminosityBlock(), meAverageLumiPerLS, meCountsDroppedPerLS, meCountsPassPerLS, meCountsStreamPerLS, meDiagnostic, meXsecStreamPerLS, nPass_, nStream_, printCountsPerPathThisLumi(), and TotalDroppedCounts.
{ int lumi = int(lumiSeg.id().luminosityBlock()); LogTrace("TrigResRateMon") << " end lumiSection number " << lumi << endl; // countHLTPathHitsEndLumiBlock(lumi); // countHLTGroupHitsEndLumiBlock(lumi); //countHLTGroupL1HitsEndLumiBlock(lumi); //countHLTGroupBXHitsEndLumiBlock(lumi); // normalizeHLTMatrix(); //if (jmsDebug) printCountsPerPathThisLumi(); if (jmsDebug) printCountsPerPathThisLumi(); if (jmsDebug) std::cout << "Average lumi is " << averageInstLumi << std::endl; if (averageInstLumi > 500) { // MonitorElement* reportSumME = dbe_->get("Info/EventInfo/reportSummaryMap" ); // TH2F * reportSum = reportSumME->getTH2F(); // float physDecl = reportSum->GetBinContent(lumi,26); fillXsecPerDataset(lumi); filltestHisto(lumi); //Robin } //Robin-------Diagnostic plots-------- // TH1F* tempXsecPerLS = meXsecPerLS->getTH1F(); // double xsec = 1.0; // // std::cout << "counts for HLT_IsoMu15* is " << testTrigCountsPS_ << std::endl; // if (averageInstLumi > 0) { // xsec = testTrigCountsPS_ / (averageInstLumi*LSsize_); // } // // std::cout << "LS is " << lumi << " ; xsec is "<< xsec << std::endl; // tempXsecPerLS->SetBinContent(lumi, xsec); // TH1F* tempXsec = meXsec->getTH1F(); // tempXsec->Fill(xsec); // TProfile* tempXsecPerIL = meXsecPerIL->getTProfile(); // tempXsecPerIL->Fill(averageInstLumi,xsec); //--------- stream counts and xsec TH1F* tempCountsStreamPerLS = meCountsStreamPerLS->getTH1F(); // std::cout << "number of stream counts is " << nStream_ << std::endl; tempCountsStreamPerLS->SetBinContent(lumi, nStream_); // dropped events MonitorElement* tempDroppedEvents = dbe_->get("SM_SMPS_Stats/droppedEventsCount_HLTTrigerResults DQM Consumer" ); if (tempDroppedEvents) { TH1F* tempDiagnostic = meDiagnostic->getTH1F(); if (tempDroppedEvents->kind() == MonitorElement::DQM_KIND_INT){ tempDiagnostic->Fill(2); int64_t tempDroppedCounts = tempDroppedEvents->getIntValue(); int64_t currentDroppedCounts = tempDroppedCounts - TotalDroppedCounts; TotalDroppedCounts = tempDroppedCounts ; TH1F* tempCountsDroppedPerLS = meCountsDroppedPerLS->getTH1F(); tempCountsDroppedPerLS->SetBinContent(lumi, currentDroppedCounts); } else tempDiagnostic->Fill(1); } else { TH1F* tempDiagnostic = meDiagnostic->getTH1F(); tempDiagnostic->Fill(0); } TH1F* tempXsecStreamPerLS = meXsecStreamPerLS->getTH1F(); double xsecStream = 1.0 ; if (averageInstLumi > 0){ xsecStream = nStream_ / (averageInstLumi*LSsize_); tempXsecStreamPerLS->SetBinContent(lumi, xsecStream); } TH1F* tempCountsPassPerLS = meCountsPassPerLS->getTH1F(); // std::cout << "number of passed stream counts is " << nPass_ << std::endl; tempCountsPassPerLS->SetBinContent(lumi, nPass_); //----------- // keep track of what you thought the lumi was // assign an error of 6% TH1F* tempLumiPerLS = meAverageLumiPerLS->getTH1F(); tempLumiPerLS->SetBinContent(lumi, averageInstLumi); tempLumiPerLS->SetBinError(lumi, averageInstLumi*0.06); }
void TrigResRateMon::endRun | ( | const edm::Run & | run, |
const edm::EventSetup & | c | ||
) | [private, virtual] |
EndRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 1017 of file TrigResRateMon.cc.
References edm::RunBase::id(), and LogDebug.
void TrigResRateMon::fillCountsPerPath | ( | const edm::Event & | iEvent, |
const edm::EventSetup & | iSetup | ||
) | [private] |
Definition at line 1443 of file TrigResRateMon.cc.
References accept(), gather_cfg::cout, dbe_, finalCountsPerPath, foundReferenceTrigger_, DQMStore::get(), hltConfig_, edm::HandleBase::isValid(), jmsDebug, HLTConfigProvider::prescaleValues(), primaryDataSetInformation, rawCountsPerPath, referenceTrigCountsPS_, referenceTrigIndex_, referenceTrigName_, HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, HLTConfigProvider::triggerIndex(), HLTConfigProvider::triggerName(), and triggerResults_.
Referenced by analyze().
{ if (jmsDebug) std::cout << "Filling counts per path" << std::endl; if (!triggerResults_.isValid()) { if (jmsDebug) std::cout << "Trigger Results not valid, sorry" << std::endl; return; } for (unsigned iName = 0; iName < hltConfig_.size(); iName++) { if ( triggerResults_ -> accept ( iName ) ){ rawCountsPerPath[iName]++; //---Robin std::string thisName = hltConfig_.triggerName(iName); std::pair<int,int> psValueCombo = hltConfig_.prescaleValues(iEvent, iSetup, thisName); // if ps OK, if ( (psValueCombo.first > 0) && (psValueCombo.second > 0) ){ finalCountsPerPath[iName] += psValueCombo.first * psValueCombo.second; } else { finalCountsPerPath[iName]++; } //----------- if ( (iName == referenceTrigIndex_) && (foundReferenceTrigger_) ) { // the get the prescales, and increment the PS*counts std::pair<int,int> psValueCombo = hltConfig_.prescaleValues(iEvent, iSetup, referenceTrigName_); // if ps OK, if ( (psValueCombo.first > 0) && (psValueCombo.second > 0) ){ referenceTrigCountsPS_ += psValueCombo.first * psValueCombo.second; } else { referenceTrigCountsPS_++; } }// end if this is reference //Robin // std::string thisName = hltConfig_.triggerName(iName); // TString checkName(thisName.c_str()); // if (checkName.Contains("HLT_IsoMu24_v")){ // std::pair<int,int> psValueCombo = hltConfig_.prescaleValues(iEvent, iSetup, thisName); // // if ps OK, // if ( (psValueCombo.first > 0) && (psValueCombo.second > 0) ){ // testTrigCountsPS_ += psValueCombo.first * psValueCombo.second; // } else { // testTrigCountsPS_++; // } // } } // end if trig fired }// end loop over paths // loop over all pds for (std::vector<DatasetInfo>::iterator iDS = primaryDataSetInformation.begin(); iDS != primaryDataSetInformation.end(); iDS++) { // now loop over all paths in the PD for (std::vector<std::string>::iterator iPath = iDS->pathNames.begin(); iPath != iDS->pathNames.end(); iPath++){ unsigned trigIndex = hltConfig_.triggerIndex(*iPath); // did you pass the trigger? if ( triggerResults_ -> accept ( trigIndex ) ){ // ok, you passed, increment the raw counts plot MonitorElement * thisRawCountsPlot = dbe_->get(iDS->rawCountsPerPathME_Name); if (thisRawCountsPlot){ iDS->fillRawCountsForPath(thisRawCountsPlot, *iPath); } else { if (jmsDebug) std::cout << "sorry but couldn't find this raw counts plot"<< iDS->datasetName << std::endl; } // the get the prescales, and increment the PS*counts std::pair<int,int> psValueCombo = hltConfig_.prescaleValues(iEvent, iSetup, *iPath); // if ps OK, if ( (psValueCombo.first > 0) && (psValueCombo.second > 0) ){ iDS->incrementCountsForPath(*iPath, psValueCombo.first*psValueCombo.second); } else { iDS->incrementCountsForPath(*iPath); } } }// end for each path }// end for each pd }
void TrigResRateMon::fillHltMatrix | ( | const edm::TriggerNames & | triggerNames, |
const edm::Event & | iEvent, | ||
const edm::EventSetup & | iSetup | ||
) | [private] |
Definition at line 1198 of file TrigResRateMon.cc.
References gather_cfg::cout, dbe_, fGroupNamePathsPair, first, g, DQMStore::get(), MonitorElement::getTH1F(), hltConfig_, i, jmsDebug, LogTrace, passAny, pathsSummaryFolder_, HLTConfigProvider::prescaleValues(), rawCountsPerPD, edm::second(), HLTConfigProvider::streamContent(), edm::TriggerNames::triggerIndex(), triggerResults_, and funct::true.
Referenced by analyze().
{ string fullPathToME; std::vector <std::pair<std::string, bool> > groupAcceptPair; // This will store the prescale values std::pair<int,int> psValueCombo; // For each dataset for (unsigned int mi=0;mi<fGroupNamePathsPair.size()-1;mi++) { // // fullPathToME = pathsSummaryFolder_ + "HLT_"+fGroupNamePathsPair[mi].first+"_PassPass"; // // MonitorElement* ME_2d = dbe_->get(fullPathToME); fullPathToME = pathsSummaryFolder_ + "HLT_"+fGroupNamePathsPair[mi].first+"_Pass_Any"; MonitorElement* ME_1d = dbe_->get(fullPathToME); // // if(!ME_2d || !ME_1d) { if(!ME_1d) { LogTrace("TrigResRateMon") << " ME not valid although I gave full path" << endl; continue; } // // TH2F * hist_2d = ME_2d->getTH2F(); TH1F * hist_1d = ME_1d->getTH1F(); // Fill HLTPassed Matrix bin (i,j) = (Any,Any) // -------------------------------------------------------- // // int anyBinNumber = hist_2d->GetXaxis()->FindBin("HLT_Any"); //string groupBinLabel = "HLT_"+fGroupNamePathsPair[mi].first+"_Any"; // // string groupBinLabel = "Total "+fGroupNamePathsPair[mi].first; // // int groupBinNumber = hist_2d->GetXaxis()->FindBin(groupBinLabel.c_str()); // any triger accepted // if(triggerResults_->accept()){ // hist_2d->Fill(anyBinNumber-1,anyBinNumber-1);//binNumber1 = 0 = first filter // hist_1d->Fill(anyBinNumber-1);//binNumber1 = 0 = first filter // } //if (jmsDebug) std::cout << "Setting histograms to HLT_ZeroBias index " << zbIndex << std::endl; // // if (zbIndex < triggerResults_->size() ) { // // if (triggerResults_->accept(zbIndex) || jmsFakeZBCounts) { // // if (jmsDebug) std::cout << "Filling bin " << (groupBinNumber-1) // // << " (out of " << hist_1d->GetNbinsX() // // << ") with ZB counts in histo " // // << hist_1d->GetName() << std::endl; // // hist_1d->Fill(groupBinNumber-1, 50000); // // hist_2d->Fill(groupBinNumber-1,groupBinNumber-1, 10000);//binNumber1 = 0 = first filter // // } // // } bool groupPassed = false; //bool groupL1Passed = false; // Main loop over paths // -------------------- //for (int i=1; i< hist_2d->GetNbinsX();i++) for (unsigned int i=0; i< fGroupNamePathsPair[mi].second.size(); i++) { //string hltPathName = hist_2d->GetXaxis()->GetBinLabel(i); string hltPathName = fGroupNamePathsPair[mi].second[i]; // check if this is hlt path name //unsigned int pathByIndex = triggerNames.triggerIndex(hltPathName); unsigned int pathByIndex = triggerNames.triggerIndex(fGroupNamePathsPair[mi].second[i]); if(pathByIndex >= triggerResults_->size() ) continue; // check if its L1 passed // comment out below but set groupL1Passed to true always //if(hasL1Passed(hltPathName,triggerNames)) groupL1Passed = true; //groupL1Passed = true; // Fill HLTPassed Matrix and HLTPassFail Matrix // -------------------------------------------------------- if(triggerResults_->accept(pathByIndex)){ groupPassed = true; //groupL1Passed = true; // // hist_2d->Fill(i,anyBinNumber-1);//binNumber1 = 0 = first filter // // hist_2d->Fill(anyBinNumber-1,i);//binNumber1 = 0 = first filter // // hist_2d->Fill(i,groupBinNumber-1);//binNumber1 = 0 = first filter // // hist_2d->Fill(groupBinNumber-1,i);//binNumber1 = 0 = first filter if (jmsDebug) std::cout << "Trying to get prescales... " << std::endl; psValueCombo = hltConfig_.prescaleValues(iEvent, iSetup, hltPathName); if (jmsDebug) std::cout << "Path " << hltPathName << " L1 PS " << psValueCombo.first << " and hlt ps " << psValueCombo.second << std::endl; if ( (psValueCombo.first > 0) && (psValueCombo.second > 0) ){ hist_1d->Fill(i, psValueCombo.first * psValueCombo.second );//binNumber1 = 0 = first filter } else { hist_1d->Fill(i); } //for (int j=1; j< hist_2d->GetNbinsY();j++) //* for (unsigned int j=0; j< fGroupNamePathsPair[mi].second.size(); j++) //* { //* string crossHltPathName = fGroupNamePathsPair[mi].second[j]; //* //unsigned int crosspathByIndex = triggerNames.triggerIndex(hist_2d->GetXaxis()->GetBinLabel(j)); //* //unsigned int crosspathByIndex = triggerNames.triggerIndex(crossHltPathName); //* unsigned int crosspathByIndex = triggerNames.triggerIndex(fGroupNamePathsPair[mi].second[j]); //* if(crosspathByIndex >= triggerResults_->size() ) continue; //* if(triggerResults_->accept(crosspathByIndex)){ //* hist_2d->Fill(i,j);//binNumber1 = 0 = first filter //* } // end if j path passed //* } // end for j } // end if i passed } // end for i if(groupPassed) { rawCountsPerPD[mi]++ ; //hist_1d->Fill(groupBinNumber-1);//binNumber1 = 0 = first filter //hist_2d->Fill(groupBinNumber-1,groupBinNumber-1);//binNumber1 = 0 = first filter //hist_2d->Fill(anyBinNumber-1,groupBinNumber-1);//binNumber1 = 0 = first filter //hist_2d->Fill(groupBinNumber-1,anyBinNumber-1);//binNumber1 = 0 = first filter } // if the group belongs to stream A // store groupName and Bool if it has passed bool isGroupFromStreamA = false; vector<string> streamDatasetNames = hltConfig_.streamContent("A") ; for (unsigned int g=0;g<streamDatasetNames.size();g++) { if(streamDatasetNames[g] == fGroupNamePathsPair[mi].first) { isGroupFromStreamA = true; break; } } if(isGroupFromStreamA) groupAcceptPair.push_back(make_pair(fGroupNamePathsPair[mi].first,groupPassed)); // L1 groups - not used anymore // string groupL1BinLabel = "HLT_"+fGroupNamePathsPair[mi].first+"_L1_Any"; // // int groupL1BinNumber = hist_2d->GetXaxis()->FindBin(groupL1BinLabel.c_str()); // // if(groupL1Passed) hist_1d->Fill(groupL1BinNumber-1);//binNumber1 = 0 = first filter } // end for mi // // fullPathToME = pathsSummaryFolder_ + "HLT_A_PassPass"; // // MonitorElement* ME_2d_Stream = dbe_->get(fullPathToME); fullPathToME = pathsSummaryFolder_ + "HLT_A_Pass_Any"; MonitorElement* ME_1d_Stream = dbe_->get(fullPathToME); // // if(!ME_2d_Stream || !ME_1d_Stream) { if(!ME_1d_Stream) { LogTrace("TrigResRateMon") << " ME not valid although I gave full path" << endl; return; } else { // // TH2F * hist_2d_Stream = ME_2d_Stream->getTH2F(); TH1F * hist_1d_Stream = ME_1d_Stream->getTH1F(); int streamBinNumber = hist_1d_Stream->GetXaxis()->GetLast(); bool acceptedStreamA = false; // loop over groups for (unsigned int i=0;i<groupAcceptPair.size();i++) { if(groupAcceptPair[i].second) { acceptedStreamA = true; // // int groupBinNumber_i = hist_2d_Stream->GetXaxis()->FindBin(groupAcceptPair[i].first.c_str()); int groupBinNumber_i = hist_1d_Stream->GetXaxis()->FindBin(groupAcceptPair[i].first.c_str()); //LogTrace("TrigResRateMon") << "Accepted group X " << groupAcceptPair[i].first.c_str() << " bin number " << groupBinNumber_i << endl; // Robin what about prescale for this one? hist_1d_Stream->Fill(groupBinNumber_i-1);//binNumber1 = 0 = first filter // // hist_2d_Stream->Fill(groupBinNumber_i-1,streamBinNumber-1);//binNumber1 = 0 = first filter // // hist_2d_Stream->Fill(streamBinNumber-1,groupBinNumber_i-1);//binNumber1 = 0 = first filter // // for (unsigned int j=0;j<groupAcceptPair.size();j++) { // // if(groupAcceptPair[j].second) { // // int groupBinNumber_j = hist_2d_Stream->GetXaxis()->FindBin(groupAcceptPair[j].first.c_str()); // // //LogTrace("TrigResRateMon") << "Accepted group Y " << groupAcceptPair[j].first.c_str() << " bin number " << groupBinNumber_j << endl; // // // fill StreamMatrix(i,j) // // hist_2d_Stream->Fill(groupBinNumber_i-1,groupBinNumber_j-1);//binNumber1 = 0 = first filter // // } // end if j-th group accepted // // } // end for j } // end if i-th group accepted } // end for i if(acceptedStreamA) { // hist_2d_Stream->Fill(streamBinNumber-1,streamBinNumber-1);//binNumber1 = 0 = first filter hist_1d_Stream->Fill(streamBinNumber-1);//binNumber1 = 0 = first filter passAny = true ; //Robin } } // end else }
void TrigResRateMon::filltestHisto | ( | const int & | lumi | ) | [private] |
Definition at line 1973 of file TrigResRateMon.cc.
References averageInstLumi, averageInstLumi3LS, finalCountsPerPath, MonitorElement::getTH1F(), MonitorElement::getTProfile(), hltConfig_, i, getHLTprescales::index, LSsize_, meXsecPerTestPath, timingPdfMaker::pathname, rawCountsPerPath, HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, testPaths_, HLTConfigProvider::triggerName(), v_ME_CountsPerLS, v_ME_Xsec, and v_ME_XsecPerLS.
Referenced by endLuminosityBlock().
{ averageInstLumi3LS += averageInstLumi; if (lumi%3 == 0){ unsigned int npaths = testPaths_.size(); for(unsigned int i = 0; i < npaths; i++){ TString pathname = testPaths_[i].c_str() ; pathname += "_v" ; int index = 0 ; int rawCount = 0 ; int finalCount = 0; double xsec = 0 ; //find the index for this test trigger path for (unsigned iName = 0; iName < hltConfig_.size(); iName++) { std::string thisName = hltConfig_.triggerName(iName); TString checkName(thisName.c_str()); if (checkName.Contains(pathname)){ index = iName ; // std::cout << "==>test path name is " << checkName << std::endl; break ; } } MonitorElement* testME_XsecPerLS = v_ME_XsecPerLS[i]; MonitorElement* testME_rawCountsPerLS = v_ME_CountsPerLS[i]; MonitorElement* testME_Xsec = v_ME_Xsec[i]; // TProfile tempProfile("XsecVsTestPath", "Xsec vs Test Path", npaths, 0.5, npaths+0.5); rawCount = rawCountsPerPath[index]; finalCount = finalCountsPerPath[index]; testME_rawCountsPerLS->getTH1F()->SetBinContent(lumi, rawCount); if (averageInstLumi > 0 ) { xsec = finalCount/ (averageInstLumi3LS*LSsize_); //averageInstLumi ??? testME_XsecPerLS->getTH1F()->SetBinContent(lumi, xsec); // testME_rawCountsPerLS->getTH1F()->SetBinContent(lumi, rawCount); testME_Xsec->getTH1F()->Fill(xsec); meXsecPerTestPath->getTProfile()->Fill(i+1,xsec); } } for (unsigned iName = 0; iName < hltConfig_.size() ; iName++) { rawCountsPerPath[iName] = 0; finalCountsPerPath[iName] = 0; //Robin } averageInstLumi3LS = 0 ; } // end if 3xLS }
void TrigResRateMon::fillXsecPerDataset | ( | const int & | lumi | ) | [private] |
PD rate plot
Definition at line 1645 of file TrigResRateMon.cc.
References averageInstLumi, gather_cfg::cout, dbe_, DQMStore::get(), MonitorElement::getTH1F(), jmsDebug, LSsize_, primaryDataSetInformation, RPCpg::rate(), rawCountsPerPD, and referenceTrigCountsPS_.
Referenced by endLuminosityBlock().
{ // calculate the reference cross section double refTrigXSec = referenceTrigCountsPS_ / ( averageInstLumi * LSsize_); // string fullpath = pathsSummaryFolder_ + "HLT_A_Pass_Any"; // MonitorElement * meStreamA = dbe_->get(fullpath); // if (!meStreamA ) std::cout << "sorry but couldn't get the stream A ME" << std::endl; int iPD = 0; for (std::vector<DatasetInfo>::iterator iDS = primaryDataSetInformation.begin(); iDS != primaryDataSetInformation.end(); iDS++) { MonitorElement * thisXsecPlot = dbe_->get(iDS->xsecPerPathME_Name); MonitorElement * scaledXsecPlot = dbe_->get(iDS->scaledXsecPerPathME_Name); if (thisXsecPlot){ iDS->fillXsecPlot(thisXsecPlot, averageInstLumi, LSsize_); } else { if (jmsDebug) std::cout << "sorry but couldn't find this xsec plot"<< iDS->datasetName << std::endl; } if (scaledXsecPlot){ iDS->fillXsecPlot(scaledXsecPlot, averageInstLumi, LSsize_, refTrigXSec); } else { if (jmsDebug) std::cout << "sorry but couldn't find this scaled xsec plot"<< iDS->datasetName << std::endl; } MonitorElement * thisRatePlot = dbe_->get(iDS->ratePerLSME_Name); if (thisRatePlot) { double rate = rawCountsPerPD[iPD] / LSsize_ ; TH1F* rateHist = thisRatePlot->getTH1F(); rateHist->SetBinContent(lumi, rate); } else { if (jmsDebug) std::cout << "sorry but couldn't find this rate plot"<< iDS->datasetName << std::endl; } rawCountsPerPD[iPD] = 0 ; iPD++; } }
void TrigResRateMon::findReferenceTriggerIndex | ( | ) | [private] |
Definition at line 1558 of file TrigResRateMon.cc.
References gather_cfg::cout, foundReferenceTrigger_, hltConfig_, jmsDebug, referenceTrigCountsPS_, referenceTrigIndex_, referenceTrigInput_, referenceTrigName_, HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerName().
Referenced by beginRun().
{ if (jmsDebug) std::cout << "Looking for reference trigger " << referenceTrigInput_ << std::endl; for (unsigned iName = 0; iName < hltConfig_.size(); iName++) { std::string thisName = hltConfig_.triggerName(iName); TString tempThisName(thisName.c_str()); if (tempThisName.Contains(referenceTrigInput_)){ referenceTrigName_ = thisName; if (jmsDebug) std::cout << "Using Reference trigger " << referenceTrigName_ << std::endl; referenceTrigIndex_ = iName; foundReferenceTrigger_ = true; referenceTrigCountsPS_ = 0; break; }// end if name contains substring } if (!foundReferenceTrigger_) { std::cout << "Sorry, we couldn't find a trigger like " << referenceTrigInput_ << std::endl; } }
const std::string TrigResRateMon::getL1ConditionModuleName | ( | const std::string & | pathname | ) | [private] |
int TrigResRateMon::getThresholdFromName | ( | const std::string & | pathname | ) | [private] |
Definition at line 2577 of file TrigResRateMon.cc.
References i, mergeVDriftHistosByStation::name, timingPdfMaker::pathname, and AlCaHLTBitMon_QueryRunRegistry::string.
{ std::string pathname = name; //cout << "----------------------------------------------" << endl; //cout << pathname << endl; //remove "L1" substr if(pathname.find("L1") != std::string::npos) pathname.replace(pathname.find("L1"),2,""); //remove "L2" substr if(pathname.find("L2") != std::string::npos) pathname.replace(pathname.find("L2"),2,""); //remove "8E29" substr if(pathname.find("8E29") != std::string::npos) pathname.replace(pathname.find("8E29"),4,""); int digitLocation=0; for (unsigned int i=0; i < pathname.length(); i++) { if (isdigit(pathname.at(i))) { digitLocation = i; break; } } // get the string from the location of the first digit to the end string hltThresholdString = pathname.substr(digitLocation); int hltThreshold = 0; // get intiger at the begining of the string sscanf (hltThresholdString.c_str(),"%d%*s",&hltThreshold); //printf ("%s -> %s -> %d\n",pathname.c_str(), hltThresholdString.c_str(), hltThreshold); return hltThreshold; }
int TrigResRateMon::getTriggerTypeParsePathName | ( | const std::string & | pathname | ) | [private] |
bool TrigResRateMon::hasHLTPassed | ( | const std::string & | pathname, |
const edm::TriggerNames & | triggerNames | ||
) | [private] |
bool TrigResRateMon::hasL1Passed | ( | const std::string & | pathname, |
const edm::TriggerNames & | triggerNames | ||
) | [private] |
void TrigResRateMon::normalizeHLTMatrix | ( | ) | [private] |
Definition at line 2615 of file TrigResRateMon.cc.
References makeLayoutFileForGui::datasetNames, dbe_, DQMStore::get(), MonitorElement::getTH1F(), MonitorElement::getTH2F(), hltConfig_, i, pathsSummaryFolder_, pathsSummaryHLTCorrelationsFolder_, HLTConfigProvider::streamContent(), v_ME_HLTPass_Normalized_Any, v_ME_HLTPassPass, and v_ME_HLTPassPass_Normalized.
{ string fullPathToME; // again, get hold of dataset names //vector<string> datasetNames = hltConfig_.datasetNames() ; vector<string> datasetNames = hltConfig_.streamContent("A") ; // fill vectors of MEs needed in normalization for (unsigned int i=0;i<datasetNames.size();i++) { fullPathToME = pathsSummaryFolder_ +"HLT_"+datasetNames[i]+"_PassPass"; v_ME_HLTPassPass.push_back( dbe_->get(fullPathToME)); fullPathToME = pathsSummaryHLTCorrelationsFolder_+"HLT_"+datasetNames[i]+"_PassPass_Normalized"; v_ME_HLTPassPass_Normalized.push_back( dbe_->get(fullPathToME)); fullPathToME = pathsSummaryHLTCorrelationsFolder_+"HLT_"+datasetNames[i]+"_Pass_Normalized_Any"; v_ME_HLTPass_Normalized_Any.push_back( dbe_->get(fullPathToME)); } // add stream MEs fullPathToME = pathsSummaryFolder_ +"HLT_A_PassPass"; v_ME_HLTPassPass.push_back( dbe_->get(fullPathToME)); fullPathToME = pathsSummaryHLTCorrelationsFolder_+"HLT_A_PassPass_Normalized"; v_ME_HLTPassPass_Normalized.push_back( dbe_->get(fullPathToME)); fullPathToME = pathsSummaryHLTCorrelationsFolder_+"HLT_A_Pass_Normalized_Any"; v_ME_HLTPass_Normalized_Any.push_back( dbe_->get(fullPathToME)); for (unsigned int i =0;i<v_ME_HLTPassPass.size();i++) { MonitorElement* ME_HLTPassPass = v_ME_HLTPassPass[i]; MonitorElement* ME_HLTPassPass_Normalized = v_ME_HLTPassPass_Normalized[i]; MonitorElement* ME_HLTPass_Normalized_Any = v_ME_HLTPass_Normalized_Any[i]; if(!ME_HLTPassPass || !ME_HLTPassPass_Normalized || !ME_HLTPass_Normalized_Any) return; float passCount = 0; unsigned int nBinsX = ME_HLTPassPass->getTH2F()->GetNbinsX(); unsigned int nBinsY = ME_HLTPassPass->getTH2F()->GetNbinsY(); for(unsigned int binX = 0; binX < nBinsX+1; binX++) { passCount = ME_HLTPassPass->getTH2F()->GetBinContent(binX,binX); for(unsigned int binY = 0; binY < nBinsY+1; binY++) { if(passCount != 0) { // normalize each bin to number of passCount float normalizedBinContentPassPass = (ME_HLTPassPass->getTH2F()->GetBinContent(binX,binY))/passCount; //float normalizedBinContentPassFail = (ME_HLTPassFail_->getTH2F()->GetBinContent(binX,binY))/passCount; ME_HLTPassPass_Normalized->getTH2F()->SetBinContent(binX,binY,normalizedBinContentPassPass); //ME_HLTPassFail_Normalized_->getTH2F()->SetBinContent(binX,binY,normalizedBinContentPassFail); if(binX == nBinsX) { ME_HLTPass_Normalized_Any->getTH1F()->SetBinContent(binY,normalizedBinContentPassPass); } } else { ME_HLTPassPass_Normalized->getTH2F()->SetBinContent(binX,binY,0); //ME_HLTPassFail_Normalized_->getTH2F()->SetBinContent(binX,binY,0); } // end if else } // end for binY } // end for binX } // end for i }
void TrigResRateMon::printCountsPerPathThisLumi | ( | ) | [private] |
Definition at line 1582 of file TrigResRateMon.cc.
References gather_cfg::cout, finalCountsPerPath, hltConfig_, primaryDataSetInformation, rawCountsPerPath, referenceTrigCountsPS_, referenceTrigIndex_, referenceTrigName_, HLTConfigProvider::size(), and HLTConfigProvider::triggerName().
Referenced by endLuminosityBlock().
{ std::cout << "===> COUNTS THIS LUMI <===" << std::endl; for (unsigned iName = 0; iName < hltConfig_.size() ; iName++) { std::cout << hltConfig_.triggerName(iName) << " = " << rawCountsPerPath[iName] << "finalCounts = " << finalCountsPerPath[iName] //Robin << std::endl; } std::cout << "+++ Reference trigger " << referenceTrigName_ << " index " << referenceTrigIndex_ << " counts " << referenceTrigCountsPS_ << std::endl; // loop over all pds for (std::vector<DatasetInfo>::const_iterator iDS = primaryDataSetInformation.begin(); iDS != primaryDataSetInformation.end(); iDS++) { iDS->printCountsPerPath(); } }
void TrigResRateMon::setupHltBxPlots | ( | ) | [private] |
Definition at line 1694 of file TrigResRateMon.cc.
References DQMStore::book2D(), dbe_, MonitorElement::getTH2F(), hltPathsDiagonal_, i, ME_HLT_BX, ME_HLT_CUSTOM_BX, Nbx_, pathsSummaryFolder_, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().
{ //pathsSummaryFolder_ = TString("HLT/TrigResults/PathsSummary/"); //dbe_->setCurrentFolder(pathsSummaryFolder_.c_str()); dbe_->setCurrentFolder(pathsSummaryFolder_); // setup HLT bx plot unsigned int npaths = hltPathsDiagonal_.size(); ME_HLT_BX = dbe_->book2D("HLT_bx", "HLT counts vs Event bx", Nbx_+1, -0.5, Nbx_+1-0.5, npaths, -0.5, npaths-0.5); ME_HLT_CUSTOM_BX = dbe_->book2D("HLT_Custom_bx", "HLT counts vs Event bx", Nbx_+1, -0.5, Nbx_+1-0.5, npaths, -0.5, npaths-0.5); ME_HLT_BX->setAxisTitle("Bunch Crossing"); ME_HLT_CUSTOM_BX->setAxisTitle("Bunch Crossing"); // Set up bin labels on Y axis continuing to cover all npaths for(unsigned int i = 0; i < npaths; i++){ ME_HLT_BX->getTH2F()->GetYaxis()->SetBinLabel(i+1, (hltPathsDiagonal_[i]).getPath().c_str()); ME_HLT_CUSTOM_BX->getTH2F()->GetYaxis()->SetBinLabel(i+1, (hltPathsDiagonal_[i]).getPath().c_str()); } }
void TrigResRateMon::setupHltLsPlots | ( | ) | [private] |
Definition at line 1725 of file TrigResRateMon.cc.
References newFWLiteAna::bin, DQMStore::book2D(), dbe_, MonitorElement::getTH2F(), hltPathsDiagonal_, i, ME_HLTAll_LS, mergeVDriftHistosByStation::name, nLS_, pathsSummaryHLTPathsPerLSFolder_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), indexGen::title, and v_ME_HLTAll_LS.
{ unsigned int npaths = hltPathsDiagonal_.size(); //pathsSummaryHLTPathsPerLSFolder_ = TString("HLT/TrigResults/PathsSummary/HLT LS/"); //dbe_->setCurrentFolder(pathsSummaryHLTPathsPerLSFolder_.c_str()); dbe_->setCurrentFolder(pathsSummaryHLTPathsPerLSFolder_); ME_HLTAll_LS = dbe_->book2D("AllSelectedPaths_count_LS", "AllSelectedPaths paths rate [Hz]", nLS_, 0, nLS_, npaths+1, -0.5, npaths+1-0.5); ME_HLTAll_LS->setAxisTitle("Luminosity Section"); // Set up bin labels on Y axis continuing to cover all npaths for(unsigned int i = 0; i < npaths; i++){ ME_HLTAll_LS->getTH2F()->GetYaxis()->SetBinLabel(i+1, (hltPathsDiagonal_[i]).getPath().c_str()); } unsigned int i = npaths; ME_HLTAll_LS->getTH2F()->GetYaxis()->SetBinLabel(i+1, "HLT_Any"); int nBinsPerLSHisto = 20; int nLSHistos = npaths/nBinsPerLSHisto; for (int nh=0;nh<nLSHistos+1;nh++) { char name[200]; char title[200]; sprintf(name, "Group_%d_paths_count_LS",nLSHistos-nh); sprintf(title, "Group %d, paths rate [Hz]",nLSHistos-nh); MonitorElement* tempME = dbe_->book2D(name,title, nLS_, 0, nLS_, nBinsPerLSHisto+3, -0.5, nBinsPerLSHisto+3-0.5); tempME->setAxisTitle("LS"); // Set up bin labels on Y axis continuing to cover all npaths for(int i = nh*nBinsPerLSHisto; i < (nh+1)*nBinsPerLSHisto; i++){ if (i == int(npaths)) break; int bin; if(nh == 0){ bin = i; } else { bin = i % nBinsPerLSHisto; } tempME->setBinLabel(bin+1, hltPathsDiagonal_[i].getPath().c_str(), 2); } tempME->setBinLabel(nBinsPerLSHisto+3, "HLT_Any", 2); tempME->setBinLabel(nBinsPerLSHisto+2, "HLT_PhysicsDeclared", 2); v_ME_HLTAll_LS.push_back(tempME); } }
void TrigResRateMon::setupHltMatrix | ( | const std::string & | label, |
std::vector< std::string > & | paths | ||
) | [private] |
Referenced by beginRun().
void TrigResRateMon::setupStreamMatrix | ( | const std::string & | label, |
std::vector< std::string > & | paths | ||
) | [private] |
Definition at line 1152 of file TrigResRateMon.cc.
References DQMStore::book1D(), DQMStore::book2D(), dbe_, MonitorElement::getTH1F(), MonitorElement::getTH2F(), i, pathsSummaryFolder_, pathsSummaryHLTCorrelationsFolder_, and DQMStore::setCurrentFolder().
{ paths.push_back(""); paths.push_back("HLT_"+label+"_Any"); string h_name; string h_title; dbe_->setCurrentFolder(pathsSummaryFolder_.c_str()); h_name= "HLT_"+label+"_PassPass"; h_title = "HLT_"+label+"_PassPass (x=Pass, y=Pass)"; MonitorElement* ME = dbe_->book2D(h_name.c_str(), h_title.c_str(), paths.size(), -0.5, paths.size()-0.5, paths.size(), -0.5, paths.size()-0.5); h_name= "HLT_"+label+"_Pass_Any"; h_title = "HLT_"+label+"_Pass (x=Pass, Any=Pass) normalized to HLT_Any Pass"; MonitorElement* ME_Any = dbe_->book1D(h_name.c_str(), h_title.c_str(), paths.size(), -0.5, paths.size()-0.5); dbe_->setCurrentFolder(pathsSummaryHLTCorrelationsFolder_.c_str()); h_name= "HLT_"+label+"_PassPass_Normalized"; h_title = "HLT_"+label+"_PassPass (x=Pass, y=Pass) normalized to xBin=Pass"; MonitorElement* ME_Normalized = dbe_->book2D(h_name.c_str(), h_title.c_str(), paths.size(), -0.5, paths.size()-0.5, paths.size(), -0.5, paths.size()-0.5); h_name= "HLT_"+label+"_Pass_Normalized_Any"; h_title = "HLT_"+label+"_Pass (x=Pass, Any=Pass) normalized to HLT_Any Pass"; MonitorElement* ME_Normalized_Any = dbe_->book1D(h_name.c_str(), h_title.c_str(), paths.size(), -0.5, paths.size()-0.5); for(unsigned int i = 0; i < paths.size(); i++){ ME->getTH2F()->GetXaxis()->SetBinLabel(i+1, (paths[i]).c_str()); ME->getTH2F()->GetYaxis()->SetBinLabel(i+1, (paths[i]).c_str()); ME_Normalized->getTH2F()->GetXaxis()->SetBinLabel(i+1, (paths[i]).c_str()); ME_Normalized->getTH2F()->GetYaxis()->SetBinLabel(i+1, (paths[i]).c_str()); ME_Any->getTH1F()->GetXaxis()->SetBinLabel(i+1, (paths[i]).c_str()); ME_Normalized_Any->getTH1F()->GetXaxis()->SetBinLabel(i+1, (paths[i]).c_str()); } }
double TrigResRateMon::averageInstLumi [private] |
Definition at line 192 of file TrigResRateMon.h.
Referenced by addLumiToAverage(), clearLumiAverage(), endLuminosityBlock(), filltestHisto(), and fillXsecPerDataset().
double TrigResRateMon::averageInstLumi3LS [private] |
Definition at line 193 of file TrigResRateMon.h.
Referenced by beginRun(), and filltestHisto().
double TrigResRateMon::bjetDRMatch_ [private] |
Definition at line 267 of file TrigResRateMon.h.
double TrigResRateMon::bjetEtaMax_ [private] |
Definition at line 265 of file TrigResRateMon.h.
double TrigResRateMon::bjetEtMin_ [private] |
Definition at line 266 of file TrigResRateMon.h.
double TrigResRateMon::bjetL1DRMatch_ [private] |
Definition at line 268 of file TrigResRateMon.h.
int TrigResRateMon::currentRun_ [private] |
Definition at line 238 of file TrigResRateMon.h.
std::vector<std::pair<std::string, std::string> > TrigResRateMon::custompathnamepairs_ [private] |
Definition at line 296 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
DQMStore* TrigResRateMon::dbe_ [private] |
Definition at line 151 of file TrigResRateMon.h.
Referenced by beginRun(), bookTestHisto(), countHLTGroupHitsEndLumiBlock(), countHLTGroupL1HitsEndLumiBlock(), countHLTPathHitsEndLumiBlock(), endLuminosityBlock(), fillCountsPerPath(), fillHltMatrix(), fillXsecPerDataset(), normalizeHLTMatrix(), setupHltBxPlots(), setupHltLsPlots(), GeneralHLTOffline::setupHltMatrix(), setupStreamMatrix(), and TrigResRateMon().
std::string TrigResRateMon::dirname_ [private] |
Definition at line 312 of file TrigResRateMon.h.
Referenced by beginJob(), beginRun(), and TrigResRateMon().
bool TrigResRateMon::doCombineRuns_ [private] |
Definition at line 236 of file TrigResRateMon.h.
Referenced by beginRun(), and TrigResRateMon().
bool TrigResRateMon::doVBTFMuon_ [private] |
Definition at line 237 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::dRMax_ [private] |
Definition at line 246 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::dRMaxElectronMuon_ [private] |
Definition at line 247 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::dxyCut_ [private] |
Definition at line 288 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::electronDRMatch_ [private] |
Definition at line 251 of file TrigResRateMon.h.
double TrigResRateMon::electronEtaMax_ [private] |
Definition at line 249 of file TrigResRateMon.h.
double TrigResRateMon::electronEtMin_ [private] |
Definition at line 250 of file TrigResRateMon.h.
double TrigResRateMon::electronL1DRMatch_ [private] |
Definition at line 252 of file TrigResRateMon.h.
std::string TrigResRateMon::fCustomBXPath [private] |
Definition at line 183 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
std::vector<std::pair<std::string, float> > TrigResRateMon::fGroupL1TempCountPair [private] |
Definition at line 303 of file TrigResRateMon.h.
Referenced by beginRun(), and countHLTGroupL1HitsEndLumiBlock().
std::vector<std::string> TrigResRateMon::fGroupName [private] |
Definition at line 185 of file TrigResRateMon.h.
Referenced by beginRun().
std::vector<std::pair<std::string, std::vector<std::string> > > TrigResRateMon::fGroupNamePathsPair [private] |
Definition at line 308 of file TrigResRateMon.h.
Referenced by beginRun(), countHLTGroupBXHitsEndLumiBlock(), countHLTGroupHitsEndLumiBlock(), countHLTGroupL1HitsEndLumiBlock(), and fillHltMatrix().
std::vector<std::pair<std::string, float> > TrigResRateMon::fGroupTempCountPair [private] |
Definition at line 302 of file TrigResRateMon.h.
Referenced by beginRun(), and countHLTGroupHitsEndLumiBlock().
std::vector<unsigned> TrigResRateMon::finalCountsPerPath [private] |
Definition at line 191 of file TrigResRateMon.h.
Referenced by bookCountsPerPath(), fillCountsPerPath(), filltestHisto(), and printCountsPerPathThisLumi().
bool TrigResRateMon::fIsSetup [private] |
Definition at line 153 of file TrigResRateMon.h.
Referenced by beginRun(), and TrigResRateMon().
bool TrigResRateMon::fLumiFlag [private] |
Definition at line 152 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
bool TrigResRateMon::found_zbIndex [private] |
Definition at line 161 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
bool TrigResRateMon::foundReferenceTrigger_ [private] |
Definition at line 222 of file TrigResRateMon.h.
Referenced by fillCountsPerPath(), findReferenceTriggerIndex(), and TrigResRateMon().
std::vector<std::pair<std::string, std::vector<int> > > TrigResRateMon::fPathBxTempCountPair [private] |
Definition at line 301 of file TrigResRateMon.h.
Referenced by beginRun(), and countHLTGroupBXHitsEndLumiBlock().
std::vector<std::pair<std::string, float> > TrigResRateMon::fPathTempCountPair [private] |
Definition at line 300 of file TrigResRateMon.h.
Referenced by beginRun(), and countHLTPathHitsEndLumiBlock().
Definition at line 94 of file TrigResRateMon.h.
HLTConfigProvider TrigResRateMon::hltConfig_ [private] |
Definition at line 320 of file TrigResRateMon.h.
Referenced by beginRun(), bookCountsPerPath(), fillCountsPerPath(), fillHltMatrix(), filltestHisto(), findReferenceTriggerIndex(), JetMETHLTOfflineSource::getL1ConditionModuleName(), normalizeHLTMatrix(), and printCountsPerPathThisLumi().
Definition at line 857 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
Definition at line 859 of file TrigResRateMon.h.
Referenced by beginRun(), FourVectorHLTOffline::hasL1Passed(), setupHltBxPlots(), and setupHltLsPlots().
double TrigResRateMon::htDRMatch_ [private] |
Definition at line 283 of file TrigResRateMon.h.
double TrigResRateMon::htEtaMax_ [private] |
Definition at line 281 of file TrigResRateMon.h.
double TrigResRateMon::htL1DRMatch_ [private] |
Definition at line 284 of file TrigResRateMon.h.
double TrigResRateMon::htMin_ [private] |
Definition at line 282 of file TrigResRateMon.h.
bool TrigResRateMon::isAlsoTrackerMuon_ [private] |
Definition at line 293 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::jetDRMatch_ [private] |
Definition at line 263 of file TrigResRateMon.h.
double TrigResRateMon::jetEtaMax_ [private] |
Definition at line 261 of file TrigResRateMon.h.
double TrigResRateMon::jetEtMin_ [private] |
Definition at line 262 of file TrigResRateMon.h.
double TrigResRateMon::jetL1DRMatch_ [private] |
Definition at line 264 of file TrigResRateMon.h.
bool TrigResRateMon::jmsDebug [private] |
Definition at line 158 of file TrigResRateMon.h.
Referenced by analyze(), beginLuminosityBlock(), beginRun(), endLuminosityBlock(), fillCountsPerPath(), fillHltMatrix(), fillXsecPerDataset(), findReferenceTriggerIndex(), and TrigResRateMon().
bool TrigResRateMon::jmsFakeZBCounts [private] |
Definition at line 159 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::LSsize_ [private] |
Definition at line 230 of file TrigResRateMon.h.
Referenced by countHLTGroupHitsEndLumiBlock(), countHLTGroupL1HitsEndLumiBlock(), countHLTPathHitsEndLumiBlock(), endLuminosityBlock(), filltestHisto(), fillXsecPerDataset(), and TrigResRateMon().
std::vector< std::string > TrigResRateMon::maskedPaths_ [private] |
Definition at line 214 of file TrigResRateMon.h.
Referenced by beginRun(), and TrigResRateMon().
MonitorElement* TrigResRateMon::ME_HLT_BX [private] |
Definition at line 164 of file TrigResRateMon.h.
Referenced by countHLTGroupBXHitsEndLumiBlock(), setupHltBxPlots(), and TrigResRateMon().
MonitorElement* TrigResRateMon::ME_HLT_CUSTOM_BX [private] |
Definition at line 165 of file TrigResRateMon.h.
Referenced by setupHltBxPlots(), and TrigResRateMon().
MonitorElement* TrigResRateMon::ME_HLTAll_LS [private] |
Definition at line 163 of file TrigResRateMon.h.
Referenced by countHLTPathHitsEndLumiBlock(), setupHltLsPlots(), and TrigResRateMon().
MonitorElement* TrigResRateMon::meAverageLumiPerLS [private] |
Definition at line 194 of file TrigResRateMon.h.
Referenced by beginRun(), and endLuminosityBlock().
Definition at line 199 of file TrigResRateMon.h.
Referenced by beginRun(), and endLuminosityBlock().
MonitorElement* TrigResRateMon::meCountsPassPerLS [private] |
Definition at line 200 of file TrigResRateMon.h.
Referenced by beginRun(), and endLuminosityBlock().
Definition at line 201 of file TrigResRateMon.h.
Referenced by beginRun(), and endLuminosityBlock().
MonitorElement* TrigResRateMon::meDiagnostic [private] |
Definition at line 198 of file TrigResRateMon.h.
Referenced by beginRun(), and endLuminosityBlock().
double TrigResRateMon::metDRMatch_ [private] |
Definition at line 279 of file TrigResRateMon.h.
double TrigResRateMon::metEtaMax_ [private] |
Definition at line 277 of file TrigResRateMon.h.
double TrigResRateMon::metL1DRMatch_ [private] |
Definition at line 280 of file TrigResRateMon.h.
double TrigResRateMon::metMin_ [private] |
Definition at line 278 of file TrigResRateMon.h.
MonitorElement* TrigResRateMon::meXsecPerTestPath [private] |
Definition at line 207 of file TrigResRateMon.h.
Referenced by bookTestHisto(), and filltestHisto().
MonitorElement* TrigResRateMon::meXsecStreamPerLS [private] |
Definition at line 202 of file TrigResRateMon.h.
Referenced by beginRun(), and endLuminosityBlock().
bool TrigResRateMon::monitorDaemon_ [private] |
Definition at line 315 of file TrigResRateMon.h.
double TrigResRateMon::muonDRMatch_ [private] |
Definition at line 255 of file TrigResRateMon.h.
double TrigResRateMon::muonEtaMax_ [private] |
Definition at line 253 of file TrigResRateMon.h.
double TrigResRateMon::muonEtMin_ [private] |
Definition at line 254 of file TrigResRateMon.h.
int TrigResRateMon::muonHitsCut_ [private] |
Definition at line 292 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::muonL1DRMatch_ [private] |
Definition at line 256 of file TrigResRateMon.h.
std::string TrigResRateMon::muonRecoCollectionName_ [private] |
Definition at line 314 of file TrigResRateMon.h.
unsigned int TrigResRateMon::nBins2D_ [private] |
Definition at line 241 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
unsigned int TrigResRateMon::nBins_ [private] |
Definition at line 240 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
unsigned int TrigResRateMon::nBinsDR_ [private] |
Definition at line 242 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
unsigned int TrigResRateMon::nBinsOneOverEt_ [private] |
Definition at line 243 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
unsigned int TrigResRateMon::Nbx_ [private] |
Definition at line 233 of file TrigResRateMon.h.
Referenced by countHLTGroupBXHitsEndLumiBlock(), setupHltBxPlots(), and TrigResRateMon().
int TrigResRateMon::nev_ [private] |
Definition at line 147 of file TrigResRateMon.h.
Referenced by analyze(), beginJob(), and endJob().
unsigned int TrigResRateMon::nLS_ [private] |
Definition at line 229 of file TrigResRateMon.h.
Referenced by beginRun(), bookTestHisto(), setupHltLsPlots(), and TrigResRateMon().
int TrigResRateMon::nMatchesCut_ [private] |
Definition at line 294 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::normalizedChi2Cut_ [private] |
Definition at line 289 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
int TrigResRateMon::nPass_ [private] |
Definition at line 149 of file TrigResRateMon.h.
Referenced by analyze(), beginLuminosityBlock(), and endLuminosityBlock().
int TrigResRateMon::nStream_ [private] |
Definition at line 148 of file TrigResRateMon.h.
Referenced by analyze(), beginLuminosityBlock(), and endLuminosityBlock().
bool TrigResRateMon::passAny [private] |
Definition at line 150 of file TrigResRateMon.h.
Referenced by analyze(), and fillHltMatrix().
std::string TrigResRateMon::pathsIndividualHLTPathsPerLSFolder_ [private] |
Definition at line 181 of file TrigResRateMon.h.
Referenced by countHLTPathHitsEndLumiBlock(), and TrigResRateMon().
std::string TrigResRateMon::pathsSummaryFilterCountsFolder_ [private] |
Definition at line 179 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
std::string TrigResRateMon::pathsSummaryFilterEfficiencyFolder_ [private] |
Definition at line 178 of file TrigResRateMon.h.
std::string TrigResRateMon::pathsSummaryFolder_ [private] |
Definition at line 175 of file TrigResRateMon.h.
Referenced by beginRun(), countHLTGroupHitsEndLumiBlock(), countHLTGroupL1HitsEndLumiBlock(), countHLTPathHitsEndLumiBlock(), fillHltMatrix(), normalizeHLTMatrix(), setupHltBxPlots(), GeneralHLTOffline::setupHltMatrix(), setupStreamMatrix(), and TrigResRateMon().
std::string TrigResRateMon::pathsSummaryHLTCorrelationsFolder_ [private] |
Definition at line 177 of file TrigResRateMon.h.
Referenced by normalizeHLTMatrix(), setupStreamMatrix(), and TrigResRateMon().
std::string TrigResRateMon::pathsSummaryHLTPathsPerBXFolder_ [private] |
Definition at line 182 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
std::string TrigResRateMon::pathsSummaryHLTPathsPerLSFolder_ [private] |
Definition at line 180 of file TrigResRateMon.h.
Referenced by countHLTGroupHitsEndLumiBlock(), countHLTGroupL1HitsEndLumiBlock(), setupHltLsPlots(), and TrigResRateMon().
std::string TrigResRateMon::pathsSummaryStreamsFolder_ [private] |
Definition at line 176 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::photonDRMatch_ [private] |
Definition at line 271 of file TrigResRateMon.h.
double TrigResRateMon::photonEtaMax_ [private] |
Definition at line 269 of file TrigResRateMon.h.
double TrigResRateMon::photonEtMin_ [private] |
Definition at line 270 of file TrigResRateMon.h.
double TrigResRateMon::photonL1DRMatch_ [private] |
Definition at line 272 of file TrigResRateMon.h.
int TrigResRateMon::pixelHitsCut_ [private] |
Definition at line 291 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
bool TrigResRateMon::plotAll_ [private] |
Definition at line 235 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
std::vector<DatasetInfo> TrigResRateMon::primaryDataSetInformation [private] |
Definition at line 495 of file TrigResRateMon.h.
Referenced by beginRun(), clearCountsPerPath(), fillCountsPerPath(), fillXsecPerDataset(), and printCountsPerPathThisLumi().
std::string TrigResRateMon::processname_ [private] |
Definition at line 313 of file TrigResRateMon.h.
Referenced by beginRun(), and TrigResRateMon().
double TrigResRateMon::ptMax_ [private] |
Definition at line 245 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::ptMin_ [private] |
Definition at line 244 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
std::vector<unsigned> TrigResRateMon::rawCountsPerPath [private] |
Definition at line 190 of file TrigResRateMon.h.
Referenced by bookCountsPerPath(), fillCountsPerPath(), filltestHisto(), and printCountsPerPathThisLumi().
std::vector<unsigned> TrigResRateMon::rawCountsPerPD [private] |
Definition at line 189 of file TrigResRateMon.h.
Referenced by beginRun(), fillHltMatrix(), and fillXsecPerDataset().
edm::InputTag TrigResRateMon::recHitsEBTag_ [private] |
Definition at line 319 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
edm::InputTag TrigResRateMon::recHitsEETag_ [private] |
Definition at line 319 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
unsigned int TrigResRateMon::referenceBX_ [private] |
Definition at line 232 of file TrigResRateMon.h.
Referenced by countHLTGroupBXHitsEndLumiBlock(), and TrigResRateMon().
unsigned TrigResRateMon::referenceTrigCountsPS_ [private] |
Definition at line 223 of file TrigResRateMon.h.
Referenced by clearCountsPerPath(), fillCountsPerPath(), fillXsecPerDataset(), findReferenceTriggerIndex(), and printCountsPerPathThisLumi().
unsigned TrigResRateMon::referenceTrigIndex_ [private] |
Definition at line 221 of file TrigResRateMon.h.
Referenced by fillCountsPerPath(), findReferenceTriggerIndex(), and printCountsPerPathThisLumi().
std::string TrigResRateMon::referenceTrigInput_ [private] |
Definition at line 219 of file TrigResRateMon.h.
Referenced by findReferenceTriggerIndex(), and TrigResRateMon().
std::string TrigResRateMon::referenceTrigName_ [private] |
Definition at line 220 of file TrigResRateMon.h.
Referenced by fillCountsPerPath(), findReferenceTriggerIndex(), and printCountsPerPathThisLumi().
MonitorElement* TrigResRateMon::scalersSelect [private] |
Definition at line 322 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
std::vector<std::string> TrigResRateMon::specialPaths_ [private] |
Definition at line 310 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::sumEtMin_ [private] |
Definition at line 285 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::tauDRMatch_ [private] |
Definition at line 259 of file TrigResRateMon.h.
double TrigResRateMon::tauEtaMax_ [private] |
Definition at line 257 of file TrigResRateMon.h.
double TrigResRateMon::tauEtMin_ [private] |
Definition at line 258 of file TrigResRateMon.h.
double TrigResRateMon::tauL1DRMatch_ [private] |
Definition at line 260 of file TrigResRateMon.h.
std::vector< std::string > TrigResRateMon::testPaths_ [private] |
Definition at line 215 of file TrigResRateMon.h.
Referenced by bookTestHisto(), filltestHisto(), and TrigResRateMon().
std::string TrigResRateMon::testPathsFolder_ [private] |
Definition at line 174 of file TrigResRateMon.h.
Referenced by bookTestHisto(), and TrigResRateMon().
int TrigResRateMon::theHLTOutputType [private] |
Definition at line 316 of file TrigResRateMon.h.
double TrigResRateMon::thresholdFactor_ [private] |
Definition at line 231 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
int64_t TrigResRateMon::TotalDroppedCounts [private] |
Definition at line 197 of file TrigResRateMon.h.
Referenced by beginRun(), and endLuminosityBlock().
double TrigResRateMon::trackDRMatch_ [private] |
Definition at line 275 of file TrigResRateMon.h.
int TrigResRateMon::trackerHitsCut_ [private] |
Definition at line 290 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
double TrigResRateMon::trackEtaMax_ [private] |
Definition at line 273 of file TrigResRateMon.h.
double TrigResRateMon::trackEtMin_ [private] |
Definition at line 274 of file TrigResRateMon.h.
double TrigResRateMon::trackL1DRMatch_ [private] |
Definition at line 276 of file TrigResRateMon.h.
std::vector<std::vector <uint> > TrigResRateMon::triggerFilterIndices_ [private] |
Definition at line 299 of file TrigResRateMon.h.
std::vector<std::vector <std::string> > TrigResRateMon::triggerFilters_ [private] |
Definition at line 298 of file TrigResRateMon.h.
Definition at line 325 of file TrigResRateMon.h.
Referenced by analyze(), fillCountsPerPath(), fillHltMatrix(), FourVectorHLTOffline::hasHLTPassed(), and FourVectorHLTOffline::hasL1Passed().
Definition at line 318 of file TrigResRateMon.h.
Referenced by analyze(), and TrigResRateMon().
Definition at line 317 of file TrigResRateMon.h.
Referenced by TrigResRateMon().
std::vector<MonitorElement*> TrigResRateMon::v_ME_CountsPerLS [private] |
Definition at line 209 of file TrigResRateMon.h.
Referenced by bookTestHisto(), and filltestHisto().
std::vector<MonitorElement*> TrigResRateMon::v_ME_HLTAll_LS [private] |
Definition at line 166 of file TrigResRateMon.h.
Referenced by countHLTPathHitsEndLumiBlock(), and setupHltLsPlots().
std::vector<MonitorElement*> TrigResRateMon::v_ME_HLTPass_Normalized_Any [private] |
Definition at line 172 of file TrigResRateMon.h.
Referenced by normalizeHLTMatrix().
std::vector<MonitorElement*> TrigResRateMon::v_ME_HLTPassPass [private] |
Definition at line 170 of file TrigResRateMon.h.
Referenced by normalizeHLTMatrix().
std::vector<MonitorElement*> TrigResRateMon::v_ME_HLTPassPass_Normalized [private] |
Definition at line 171 of file TrigResRateMon.h.
Referenced by normalizeHLTMatrix().
std::vector<MonitorElement*> TrigResRateMon::v_ME_Total_BX [private] |
Definition at line 167 of file TrigResRateMon.h.
Referenced by countHLTGroupBXHitsEndLumiBlock().
std::vector<MonitorElement*> TrigResRateMon::v_ME_Total_BX_Norm [private] |
Definition at line 168 of file TrigResRateMon.h.
Referenced by countHLTGroupBXHitsEndLumiBlock().
std::vector<MonitorElement*> TrigResRateMon::v_ME_Xsec [private] |
Definition at line 210 of file TrigResRateMon.h.
Referenced by bookTestHisto(), and filltestHisto().
std::vector<MonitorElement*> TrigResRateMon::v_ME_XsecPerLS [private] |
Definition at line 208 of file TrigResRateMon.h.
Referenced by bookTestHisto(), and filltestHisto().
unsigned int TrigResRateMon::zbIndex [private] |
Definition at line 160 of file TrigResRateMon.h.