CMS 3D CMS Logo

Classes | Public Member Functions | Public Attributes | Private Member Functions | Private Attributes

TrigResRateMon Class Reference

#include <TrigResRateMon.h>

Inheritance diagram for TrigResRateMon:
edm::EDAnalyzer edm::EDConsumerBase

List of all members.

Classes

class  DatasetInfo
class  PathInfo
class  PathInfoCollection

Public Member Functions

 TrigResRateMon (const edm::ParameterSet &)
 ~TrigResRateMon ()

Public Attributes

edm::Handle
< trigger::TriggerEvent
fTriggerObj
PathInfoCollection hltPaths_
PathInfoCollection hltPathsDiagonal_

Private Member Functions

void addLumiToAverage (double lumi)
virtual void analyze (const edm::Event &, const edm::EventSetup &)
virtual void beginJob ()
void beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
void beginRun (const edm::Run &run, const edm::EventSetup &c)
void bookCountsPerPath ()
void bookTestHisto ()
void clearCountsPerPath ()
void clearLumiAverage ()
void countHLTGroupBXHitsEndLumiBlock (const int &lumi)
void countHLTGroupHitsEndLumiBlock (const int &lumi)
void countHLTGroupL1HitsEndLumiBlock (const int &lumi)
void countHLTPathHitsEndLumiBlock (const int &lumi)
virtual void endJob ()
void endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
void endRun (const edm::Run &run, const edm::EventSetup &c)
 EndRun.
void fillCountsPerPath (const edm::Event &iEvent, const edm::EventSetup &iSetup)
void fillHltMatrix (const edm::TriggerNames &triggerNames, const edm::Event &iEvent, const edm::EventSetup &iSetup)
void filltestHisto (const int &lumi)
void fillXsecPerDataset (const int &lumi)
void findReferenceTriggerIndex ()
const std::string getL1ConditionModuleName (const std::string &pathname)
int getThresholdFromName (const std::string &pathname)
int getTriggerTypeParsePathName (const std::string &pathname)
bool hasHLTPassed (const std::string &pathname, const edm::TriggerNames &triggerNames)
bool hasL1Passed (const std::string &pathname, const edm::TriggerNames &triggerNames)
void normalizeHLTMatrix ()
void printCountsPerPathThisLumi ()
void setupHltBxPlots ()
void setupHltLsPlots ()
void setupHltMatrix (const std::string &label, std::vector< std::string > &paths)
void setupStreamMatrix (const std::string &label, std::vector< std::string > &paths)

Private Attributes

double averageInstLumi
double averageInstLumi3LS
double bjetDRMatch_
double bjetEtaMax_
double bjetEtMin_
double bjetL1DRMatch_
int currentRun_
std::vector< std::pair
< std::string, std::string > > 
custompathnamepairs_
DQMStoredbe_
std::string dirname_
bool doCombineRuns_
bool doVBTFMuon_
double dRMax_
double dRMaxElectronMuon_
double dxyCut_
double electronDRMatch_
double electronEtaMax_
double electronEtMin_
double electronL1DRMatch_
std::string fCustomBXPath
std::vector< std::pair
< std::string, float > > 
fGroupL1TempCountPair
std::vector< std::string > fGroupName
std::vector< std::pair
< std::string, std::vector
< std::string > > > 
fGroupNamePathsPair
std::vector< std::pair
< std::string, float > > 
fGroupTempCountPair
std::vector< unsigned > finalCountsPerPath
bool fIsSetup
bool fLumiFlag
bool found_zbIndex
bool foundReferenceTrigger_
std::vector< std::pair
< std::string, std::vector
< int > > > 
fPathBxTempCountPair
std::vector< std::pair
< std::string, float > > 
fPathTempCountPair
HLTConfigProvider hltConfig_
double htDRMatch_
double htEtaMax_
double htL1DRMatch_
double htMin_
bool isAlsoTrackerMuon_
double jetDRMatch_
double jetEtaMax_
double jetEtMin_
double jetL1DRMatch_
bool jmsDebug
bool jmsFakeZBCounts
double LSsize_
std::vector< std::string > maskedPaths_
MonitorElementME_HLT_BX
MonitorElementME_HLT_CUSTOM_BX
MonitorElementME_HLTAll_LS
MonitorElementmeAverageLumiPerLS
MonitorElementmeCountsDroppedPerLS
MonitorElementmeCountsPassPerLS
MonitorElementmeCountsStreamPerLS
MonitorElementmeDiagnostic
double metDRMatch_
double metEtaMax_
double metL1DRMatch_
double metMin_
MonitorElementmeXsecPerTestPath
MonitorElementmeXsecStreamPerLS
bool monitorDaemon_
double muonDRMatch_
double muonEtaMax_
double muonEtMin_
int muonHitsCut_
double muonL1DRMatch_
std::string muonRecoCollectionName_
unsigned int nBins2D_
unsigned int nBins_
unsigned int nBinsDR_
unsigned int nBinsOneOverEt_
unsigned int Nbx_
int nev_
unsigned int nLS_
int nMatchesCut_
double normalizedChi2Cut_
int nPass_
int nStream_
bool passAny
std::string pathsIndividualHLTPathsPerLSFolder_
std::string pathsSummaryFilterCountsFolder_
std::string pathsSummaryFilterEfficiencyFolder_
std::string pathsSummaryFolder_
std::string pathsSummaryHLTCorrelationsFolder_
std::string pathsSummaryHLTPathsPerBXFolder_
std::string pathsSummaryHLTPathsPerLSFolder_
std::string pathsSummaryStreamsFolder_
double photonDRMatch_
double photonEtaMax_
double photonEtMin_
double photonL1DRMatch_
int pixelHitsCut_
bool plotAll_
std::vector< DatasetInfoprimaryDataSetInformation
std::string processname_
double ptMax_
double ptMin_
std::vector< unsigned > rawCountsPerPath
std::vector< unsigned > rawCountsPerPD
edm::InputTag recHitsEBTag_
edm::InputTag recHitsEETag_
unsigned int referenceBX_
unsigned referenceTrigCountsPS_
unsigned referenceTrigIndex_
std::string referenceTrigInput_
std::string referenceTrigName_
MonitorElementscalersSelect
std::vector< std::string > specialPaths_
double sumEtMin_
double tauDRMatch_
double tauEtaMax_
double tauEtMin_
double tauL1DRMatch_
std::vector< std::string > testPaths_
std::string testPathsFolder_
int theHLTOutputType
double thresholdFactor_
int64_t TotalDroppedCounts
double trackDRMatch_
int trackerHitsCut_
double trackEtaMax_
double trackEtMin_
double trackL1DRMatch_
std::vector< std::vector< uint > > triggerFilterIndices_
std::vector< std::vector
< std::string > > 
triggerFilters_
edm::Handle< edm::TriggerResultstriggerResults_
edm::InputTag triggerResultsLabel_
edm::InputTag triggerSummaryLabel_
std::vector< MonitorElement * > v_ME_CountsPerLS
std::vector< MonitorElement * > v_ME_HLTAll_LS
std::vector< MonitorElement * > v_ME_HLTPass_Normalized_Any
std::vector< MonitorElement * > v_ME_HLTPassPass
std::vector< MonitorElement * > v_ME_HLTPassPass_Normalized
std::vector< MonitorElement * > v_ME_Total_BX
std::vector< MonitorElement * > v_ME_Total_BX_Norm
std::vector< MonitorElement * > v_ME_Xsec
std::vector< MonitorElement * > v_ME_XsecPerLS
unsigned int zbIndex

Detailed Description

Module to monitor rates from TriggerResults

Implementation: <Notes on="" implementation>="">

Definition at line 86 of file TrigResRateMon.h.


Constructor & Destructor Documentation

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

}

Member Function Documentation

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_.

{
   LogInfo("TrigResRateMon") << "analyzed " << nev_ << " events";
   return;
}
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.

{

  LogDebug("TrigResRateMon") << "endRun, run " << run.id();

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

    }

}

Member Data Documentation

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.

Definition at line 268 of file TrigResRateMon.h.

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

std::string TrigResRateMon::dirname_ [private]

Definition at line 312 of file TrigResRateMon.h.

Referenced by beginJob(), beginRun(), and TrigResRateMon().

Definition at line 236 of file TrigResRateMon.h.

Referenced by beginRun(), and TrigResRateMon().

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

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

Definition at line 251 of file TrigResRateMon.h.

Definition at line 249 of file TrigResRateMon.h.

Definition at line 250 of file TrigResRateMon.h.

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]
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]
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().

Definition at line 161 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

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.

Definition at line 857 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

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.

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.

Definition at line 264 of file TrigResRateMon.h.

bool TrigResRateMon::jmsDebug [private]

Definition at line 159 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

double TrigResRateMon::LSsize_ [private]
std::vector< std::string > TrigResRateMon::maskedPaths_ [private]

Definition at line 214 of file TrigResRateMon.h.

Referenced by beginRun(), and TrigResRateMon().

Definition at line 165 of file TrigResRateMon.h.

Referenced by setupHltBxPlots(), and TrigResRateMon().

Definition at line 163 of file TrigResRateMon.h.

Referenced by countHLTPathHitsEndLumiBlock(), setupHltLsPlots(), and TrigResRateMon().

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

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

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.

Definition at line 280 of file TrigResRateMon.h.

double TrigResRateMon::metMin_ [private]

Definition at line 278 of file TrigResRateMon.h.

Definition at line 207 of file TrigResRateMon.h.

Referenced by bookTestHisto(), and filltestHisto().

Definition at line 202 of file TrigResRateMon.h.

Referenced by beginRun(), and endLuminosityBlock().

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.

Definition at line 292 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

Definition at line 256 of file TrigResRateMon.h.

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]
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().

Definition at line 294 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

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

Definition at line 181 of file TrigResRateMon.h.

Referenced by countHLTPathHitsEndLumiBlock(), and TrigResRateMon().

Definition at line 179 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

Definition at line 178 of file TrigResRateMon.h.

std::string TrigResRateMon::pathsSummaryFolder_ [private]

Definition at line 177 of file TrigResRateMon.h.

Referenced by normalizeHLTMatrix(), setupStreamMatrix(), and TrigResRateMon().

Definition at line 182 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

Definition at line 176 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

Definition at line 271 of file TrigResRateMon.h.

Definition at line 269 of file TrigResRateMon.h.

double TrigResRateMon::photonEtMin_ [private]

Definition at line 270 of file TrigResRateMon.h.

Definition at line 272 of file TrigResRateMon.h.

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::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]
std::vector<unsigned> TrigResRateMon::rawCountsPerPD [private]

Definition at line 189 of file TrigResRateMon.h.

Referenced by beginRun(), fillHltMatrix(), and fillXsecPerDataset().

Definition at line 319 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

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

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

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

Definition at line 316 of file TrigResRateMon.h.

Definition at line 231 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

Definition at line 197 of file TrigResRateMon.h.

Referenced by beginRun(), and endLuminosityBlock().

Definition at line 275 of file TrigResRateMon.h.

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.

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 318 of file TrigResRateMon.h.

Referenced by analyze(), and TrigResRateMon().

Definition at line 317 of file TrigResRateMon.h.

Referenced by TrigResRateMon().

Definition at line 209 of file TrigResRateMon.h.

Referenced by bookTestHisto(), and filltestHisto().

Definition at line 166 of file TrigResRateMon.h.

Referenced by countHLTPathHitsEndLumiBlock(), and setupHltLsPlots().

Definition at line 172 of file TrigResRateMon.h.

Referenced by normalizeHLTMatrix().

Definition at line 170 of file TrigResRateMon.h.

Referenced by normalizeHLTMatrix().

Definition at line 171 of file TrigResRateMon.h.

Referenced by normalizeHLTMatrix().

Definition at line 167 of file TrigResRateMon.h.

Referenced by countHLTGroupBXHitsEndLumiBlock().

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

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.