FourVectorHLT Class Reference

#include <DQM/FourVectorHLT/src/FourVectorHLT.cc>

Inheritance diagram for FourVectorHLT:

Classes
class	PathInfo
class	PathInfoCollection

Public Member Functions
	FourVectorHLT (const edm::ParameterSet &)
	~FourVectorHLT ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
void	beginRun (const edm::Run &run, const edm::EventSetup &c)
virtual void	endJob ()
void	endRun (const edm::Run &run, const edm::EventSetup &c)
	EndRun. More...

Private Attributes
DQMStore *	dbe_
std::string	dirname_
PathInfoCollection	hltPaths_
unsigned int	nBins_
int	nev_
bool	plotAll_
double	ptMax_
double	ptMin_
edm::InputTag	triggerResultLabel_
edm::InputTag	triggerSummaryLabel_
edm::EDGetTokenT< trigger::TriggerEvent >	triggerSummaryToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: This is a DQM source meant to plot high-level HLT trigger quantities as stored in the HLT results object TriggerResults

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

Definition at line 49 of file FourVectorHLT.h.

Constructor & Destructor Documentation

FourVectorHLT::FourVectorHLT ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 15 of file FourVectorHLT.cc.

References dbe_, filters, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), LogDebug, Utilities::operator, AlignmentTrackSelector_cfi::ptMax, ALCARECOTkAlBeamHalo_cff::ptMin, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  LogDebug("FourVectorHLT") << "constructor...." ;

  dbe_ = Service < DQMStore > ().operator->();
  if ( ! dbe_ ) {
    LogWarning("Status") << "unable to get DQMStore service?";
  }
  if (iConfig.getUntrackedParameter < bool > ("DQMStore", false)) {
    dbe_->setVerbose(0);
  }
  
  
  dirname_="HLT/FourVectorHLT" ;
  
  if (dbe_ != 0 ) {
    LogDebug("Status") << "Setting current directory to " << dirname_;
    dbe_->setCurrentFolder(dirname_);
  }
  
  
  // plotting paramters
  ptMin_ = iConfig.getUntrackedParameter<double>("ptMin",0.);
  ptMax_ = iConfig.getUntrackedParameter<double>("ptMax",200.);
  nBins_ = iConfig.getUntrackedParameter<unsigned int>("Nbins",50);
  
  plotAll_ = iConfig.getUntrackedParameter<bool>("plotAll", false);

  // this is the list of paths to look at.
  std::vector<edm::ParameterSet> filters = 
    iConfig.getParameter<std::vector<edm::ParameterSet> >("filters");
  for(std::vector<edm::ParameterSet>::iterator 
    filterconf = filters.begin() ; filterconf != filters.end(); 
      filterconf++) {
    std::string me  = filterconf->getParameter<std::string>("name");
    int objectType = filterconf->getParameter<unsigned int>("type");
    float ptMin = filterconf->getUntrackedParameter<double>("ptMin");
    float ptMax = filterconf->getUntrackedParameter<double>("ptMax");
    hltPaths_.push_back(PathInfo(me, objectType, ptMin, ptMax));
  }
  if ( hltPaths_.size() && plotAll_) {
    // these two ought to be mutually exclusive....
    LogWarning("Configuration") << "Using both plotAll and a list. "
      "list will be ignored." ;
    hltPaths_.clear();
  }
  triggerSummaryLabel_ = 
    iConfig.getParameter<edm::InputTag>("triggerSummaryLabel");

  //set Token(-s)
  triggerSummaryToken_ = consumes<trigger::TriggerEvent>(iConfig.getParameter<edm::InputTag>("triggerSummaryLabel"));
}

FourVectorHLT::~FourVectorHLT

(

)

Definition at line 69 of file FourVectorHLT.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void FourVectorHLT::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Definition at line 84 of file FourVectorHLT.cc.

References dbe_, MillePedeFileConverter_cfg::e, edm::InputTag::encode(), stringResolutionProvider_cfi::et, PVValHelper::eta, trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), diffTreeTool::index, edm::HandleBase::isValid(), gen::k, LogDebug, dataset::name, AlCaHLTBitMon_ParallelJobs::p, EnergyCorrector::pt, trigger::TriggerEvent::sizeFilters(), AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, and findQualityFiles::v.

{
  using namespace edm;
  using namespace trigger;
  ++nev_;
  LogDebug("Status")<< "analyze" ;
  
  edm::Handle<TriggerEvent> triggerObj;
  iEvent.getByToken(triggerSummaryToken_, triggerObj);
  if(!triggerObj.isValid()) { 
    edm::LogInfo("Status") << "Summary HLT object (TriggerEvent) not found, "
      "skipping event"; 
    return;
  }
  
  const trigger::TriggerObjectCollection & toc(triggerObj->getObjects());

  if ( plotAll_ ) {
    for ( size_t ia = 0; ia < triggerObj->sizeFilters(); ++ ia) {
      std::string fullname = triggerObj->filterTag(ia).encode();
      // the name can have in it the module label as well as the process and
      // other labels - strip 'em
      std::string name;
      size_t p = fullname.find_first_of(':');
      if ( p != std::string::npos) {
    name = fullname.substr(0, p);
      }
      else {
    name = fullname;
      }

      LogDebug("Parameter") << "filter " << ia << ", full name = " << fullname
                << ", p = " << p 
                << ", abbreviated = " << name ;
      
      PathInfoCollection::iterator pic =  hltPaths_.find(name);
      if ( pic == hltPaths_.end() ) {
    // doesn't exist - add it
    MonitorElement *et(0), *eta(0), *phi(0), *etavsphi(0);
    
    std::string histoname(name+"_et");
    LogDebug("Status") << "new histo with name "<<  histoname ;
    dbe_->setCurrentFolder(dirname_);
    std::string title(name+" E_{T}");
    et =  dbe_->book1D(histoname.c_str(),
              title.c_str(),nBins_, 0, 100);
      
    histoname = name+"_eta";
    title = name+" #eta";
    eta =  dbe_->book1D(histoname.c_str(),
               title.c_str(),nBins_,-2.7,2.7);
      
    histoname = name+"_phi";
    title = name+" #phi";
    phi =  dbe_->book1D(histoname.c_str(),
               title.c_str(),nBins_,-3.14,3.14);
      
      
    histoname = name+"_etaphi";
    title = name+" #eta vs #phi";
    etavsphi =  dbe_->book2D(histoname.c_str(),
                title.c_str(),
                nBins_,-2.7,2.7,
                nBins_,-3.14, 3.14);
      
    // no idea how to get the bin boundries in this mode
    PathInfo e(name,0, et, eta, phi, etavsphi, ptMin_,ptMax_);
    hltPaths_.push_back(e);  
    pic = hltPaths_.begin() + hltPaths_.size()-1;
      }
      const trigger::Keys & k = triggerObj->filterKeys(ia);
      for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
    LogDebug("Parameters") << "pt, eta, phi = " 
                   << toc[*ki].pt() << ", "
                   << toc[*ki].eta() << ", "
                   << toc[*ki].phi() ;
    pic->getEtHisto()->Fill(toc[*ki].pt());
    pic->getEtaHisto()->Fill(toc[*ki].eta());
    pic->getPhiHisto()->Fill(toc[*ki].phi());
    pic->getEtaVsPhiHisto()->Fill(toc[*ki].eta(), toc[*ki].phi());
      }  

    }

  }
  else { // not plotAll_
    for(PathInfoCollection::iterator v = hltPaths_.begin();
    v!= hltPaths_.end(); ++v ) {
      const int index = triggerObj->filterIndex(v->getName());
      if ( index >= triggerObj->sizeFilters() ) {
    continue; // not in this event
      }
      LogDebug("Status") << "filling ... " ;
      const trigger::Keys & k = triggerObj->filterKeys(index);
      for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
    v->getEtHisto()->Fill(toc[*ki].pt());
    v->getEtaHisto()->Fill(toc[*ki].eta());
    v->getPhiHisto()->Fill(toc[*ki].phi());
    v->getEtaVsPhiHisto()->Fill(toc[*ki].eta(), toc[*ki].phi());
      }  
    }
  }
}

void FourVectorHLT::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 191 of file FourVectorHLT.cc.

References fastjetJetProducer_validation_cfg::DQMStore, stringResolutionProvider_cfi::et, PVValHelper::eta, Utilities::operator, AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, and findQualityFiles::v.

{
  nev_ = 0;
  DQMStore *dbe = 0;
  dbe = Service<DQMStore>().operator->();
  
  if (dbe) {
    dbe->setCurrentFolder(dirname_);
    dbe->rmdir(dirname_);
  }
  
  
  if (dbe) {
    dbe->setCurrentFolder(dirname_);

    if ( ! plotAll_ ) {
      for(PathInfoCollection::iterator v = hltPaths_.begin();
      v!= hltPaths_.end(); ++v ) {
    MonitorElement *et, *eta, *phi, *etavsphi=0;
    std::string histoname(v->getName()+"_et");
    std::string title(v->getName()+" E_t");
    et =  dbe->book1D(histoname.c_str(),
              title.c_str(),nBins_,
              v->getPtMin(),
              v->getPtMax());
      
    histoname = v->getName()+"_eta";
    title = v->getName()+" #eta";
    eta =  dbe->book1D(histoname.c_str(),
               title.c_str(),nBins_,-2.7,2.7);

    histoname = v->getName()+"_phi";
    title = v->getName()+" #phi";
    phi =  dbe->book1D(histoname.c_str(),
               histoname.c_str(),nBins_,-3.14,3.14);
 

    histoname = v->getName()+"_etaphi";
    title = v->getName()+" #eta vs #phi";
    etavsphi =  dbe->book2D(histoname.c_str(),
                title.c_str(),
                nBins_,-2.7,2.7,
                nBins_,-3.14, 3.14);
      
    v->setHistos( et, eta, phi, etavsphi);
      } 
    } // ! plotAll_ - for plotAll we discover it during the event
  }
}

void FourVectorHLT::beginRun ( const edm::Run &  run, const edm::EventSetup &  c  )

private

Definition at line 251 of file FourVectorHLT.cc.

References edm::RunBase::id(), and LogDebug.

{
  LogDebug("Status") << "beginRun, run " << run.id();
}

void FourVectorHLT::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 243 of file FourVectorHLT.cc.

{
   LogInfo("Status") << "endJob: analyzed " << nev_ << " events";
   return;
}

void FourVectorHLT::endRun ( const edm::Run &  run, const edm::EventSetup &  c  )

private

EndRun.

Definition at line 257 of file FourVectorHLT.cc.

References edm::RunBase::id(), and LogDebug.

{
  LogDebug("Status") << "endRun, run " << run.id();
}

Member Data Documentation

DQMStore* FourVectorHLT::dbe_

private

Definition at line 69 of file FourVectorHLT.h.

std::string FourVectorHLT::dirname_

private

Definition at line 77 of file FourVectorHLT.h.

PathInfoCollection FourVectorHLT::hltPaths_

private

Definition at line 166 of file FourVectorHLT.h.

unsigned int FourVectorHLT::nBins_

private

Definition at line 73 of file FourVectorHLT.h.

int FourVectorHLT::nev_

private

Definition at line 68 of file FourVectorHLT.h.

bool FourVectorHLT::plotAll_

private

Definition at line 71 of file FourVectorHLT.h.

double FourVectorHLT::ptMax_

private

Definition at line 75 of file FourVectorHLT.h.

double FourVectorHLT::ptMin_

private

Definition at line 74 of file FourVectorHLT.h.

edm::InputTag FourVectorHLT::triggerResultLabel_

private

Definition at line 79 of file FourVectorHLT.h.

edm::InputTag FourVectorHLT::triggerSummaryLabel_

private

Definition at line 78 of file FourVectorHLT.h.

edm::EDGetTokenT<trigger::TriggerEvent> FourVectorHLT::triggerSummaryToken_

private

Definition at line 82 of file FourVectorHLT.h.