Inheritance diagram for AlCaHEMuonFilter:

Public Member Functions
	AlCaHEMuonFilter (edm::ParameterSet const &, const AlCaHEMuons::Counters *count)

void	endStream () override

bool	filter (edm::Event &, edm::EventSetup const &) override

	~AlCaHEMuonFilter () override

Public Member Functions inherited from edm::stream::EDFilter< edm::GlobalCache< AlCaHEMuons::Counters > >
	EDFilter ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

	~EDFilterBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

static void	globalEndJob (const AlCaHEMuons::Counters *counters)

static std::unique_ptr< AlCaHEMuons::Counters >	initializeGlobalCache (edm::ParameterSet const &)

Static Public Member Functions inherited from edm::stream::EDFilterBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	endRun (edm::Run const &, edm::EventSetup const &) override

Private Attributes
double	caloIsoCut_

HLTConfigProvider	hltConfig_

std::vector< std::string >	HLTNames_

edm::InputTag	labelMuon_

double	muonetaCut_

double	muonptCut_

unsigned int	nAll_

unsigned int	nFinal_

unsigned int	nGood_

unsigned int	nRun_

bool	pfCut_

double	pfIsoCut_

int	preScale_

std::string	processName_

edm::EDGetTokenT< reco::MuonCollection >	tok_Muon_

edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt

edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes_

double	trackIsoCut_

edm::InputTag	triggerResults_

std::vector< std::string >	trigNames_

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter< edm::GlobalCache< AlCaHEMuons::Counters > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

Definition at line 49 of file AlCaHEMuonFilter.cc.

Constructor & Destructor Documentation

AlCaHEMuonFilter::AlCaHEMuonFilter	(	edm::ParameterSet const &	iConfig,
		const AlCaHEMuons::Counters *	count
	)

explicit

Definition at line 92 of file AlCaHEMuonFilter.cc.

References caloIsoCut_, edm::ParameterSet::getParameter(), gen::k, labelMuon_, muonetaCut_, muonptCut_, pfCut_, pfIsoCut_, preScale_, processName_, AlCaHLTBitMon_QueryRunRegistry::string, tok_Muon_, tok_trigRes_, trackIsoCut_, triggerResults_, and trigNames_.

     : nRun_(0), nAll_(0), nGood_(0), nFinal_(0) {
   //now do what ever initialization is needed
   trigNames_ = iConfig.getParameter<std::vector<std::string> >("triggers");
   processName_ = iConfig.getParameter<std::string>("processName");
   triggerResults_ = iConfig.getParameter<edm::InputTag>("triggerResultLabel");
   labelMuon_ = iConfig.getParameter<edm::InputTag>("muonLabel");
   muonptCut_ = iConfig.getParameter<double>("muonPtCut");
   muonetaCut_ = iConfig.getParameter<double>("muonEtaCut");
   pfCut_ = iConfig.getParameter<bool>("pfCut");
   pfIsoCut_ = iConfig.getParameter<double>("pfIsolationCut");
   trackIsoCut_ = iConfig.getParameter<double>("trackIsolationCut");
   caloIsoCut_ = iConfig.getParameter<double>("caloIsolationCut");
   preScale_ = iConfig.getParameter<int>("preScale");
   if (preScale_ < 1)
     preScale_ = 1;
 
   // define tokens for access
   tok_trigRes_ = consumes<edm::TriggerResults>(triggerResults_);
   tok_Muon_ = consumes<reco::MuonCollection>(labelMuon_);
   edm::LogInfo("HEMuon") << "Parameters read from config file \n"
                          << "Process " << processName_ << "  Prescale " << preScale_ << "  Isolation Cuts "
                          << trackIsoCut_ << ":" << caloIsoCut_ << "\n";
   for (unsigned int k = 0; k < trigNames_.size(); ++k)
     edm::LogInfo("HEMuon") << "Trigger[" << k << "] " << trigNames_[k] << "\n";
 }  // AlCaHEMuonFilter::AlCaHEMuonFilter  constructor

AlCaHEMuonFilter::~AlCaHEMuonFilter ( )

override

Definition at line 119 of file AlCaHEMuonFilter.cc.

119 {}

Member Function Documentation

void AlCaHEMuonFilter::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 239 of file AlCaHEMuonFilter.cc.

References RemoveAddSevLevel::flag, hltConfig_, HLTConfigProvider::init(), nRun_, processName_, and edm::RunBase::run().

                                                                                {
   bool changed(false);
   bool flag = hltConfig_.init(iRun, iSetup, processName_, changed);
   edm::LogInfo("HEMuon") << "Run[" << nRun_ << "] " << iRun.run() << " hltconfig.init " << flag << std::endl;
 }

void AlCaHEMuonFilter::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &
	)

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 246 of file AlCaHEMuonFilter.cc.

References nRun_, and edm::RunBase::run().

                                                                       {
   edm::LogInfo("HEMuon") << "endRun[" << nRun_ << "] " << iRun.run() << "\n";
   nRun_++;
 }

void AlCaHEMuonFilter::endStream ( )

overridevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 227 of file AlCaHEMuonFilter.cc.

References nAll_, nFinal_, and nGood_.

                                  {
   globalCache()->nAll_ += nAll_;
   globalCache()->nGood_ += nGood_;
   globalCache()->nFinal_ += nFinal_;
 }

void AlCaHEMuonFilter::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 252 of file AlCaHEMuonFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                   {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   std::vector<std::string> triggers = {"HLT_IsoMu", "HLT_Mu"};
   desc.add<std::string>("processName", "HLT");
   desc.add<edm::InputTag>("triggerResultLabel", edm::InputTag("TriggerResults", "", "HLT"));
   desc.add<edm::InputTag>("muonLabel", edm::InputTag("muons"));
   desc.add<double>("minimumMuonP", 10.0);
   desc.add<std::vector<std::string> >("triggers", triggers);
   desc.add<double>("muonPtCut", 20.0);
   desc.add<double>("muonEtaCut", 1.305);
   desc.add<bool>("pfCut", true);
   desc.add<double>("pfIsolationCut", 0.15);
   desc.add<double>("trackIsolationCut", 3.0);
   desc.add<double>("caloIsolationCut", 5.0);
   desc.add<int>("preScale", 1);
   descriptions.add("alcaHEMuonFilter", desc);
 }

bool AlCaHEMuonFilter::filter	(	edm::Event &	iEvent,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Implements edm::stream::EDFilterBase.

Definition at line 126 of file AlCaHEMuonFilter.cc.

References accept(), edm::HLTGlobalStatus::accept(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), caloIsoCut_, fftjetproducer_cfi::etaCut, edm::EventID::event(), spr::find(), edm::EventSetup::get(), edm::Event::getByToken(), ValidationMatrix::hlt, mps_fire::i, edm::EventBase::id(), edm::HandleBase::isValid(), edm::EventBase::luminosityBlock(), SiStripPI::max, electrons_cff::modules, muonetaCut_, muonptCut_, nAll_, nFinal_, nGood_, convertSQLiteXML::ok, spr::propagatedTrackID::okECAL, spr::propagatedTrackID::okHCAL, pfCut_, pfIsoCut_, preScale_, edm::ESHandle< T >::product(), spr::propagateCALO(), HiGenCleaner_cff::ptCut, edm::EventID::run(), edm::HLTGlobalStatus::size(), tok_Muon_, tok_trigRes_, trackIsoCut_, edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), triggerResults, and trigNames_.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

                                                                            {
   bool accept(false);
   ++nAll_;
 #ifdef DebugLog
   edm::LogInfo("HEMuon") << "AlCaHEMuonFilter::Run " << iEvent.id().run() << " Event " << iEvent.id().event()
                          << " Luminosity " << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing()
                          << std::endl;
 #endif
   //Step1: Find if the event passes one of the chosen triggers
   edm::Handle<edm::TriggerResults> triggerResults;
   iEvent.getByToken(tok_trigRes_, triggerResults);
   if (triggerResults.isValid()) {
     bool ok(false);
     std::vector<std::string> modules;
     const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
     const std::vector<std::string>& triggerNames_ = triggerNames.triggerNames();
     for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
       int hlt = triggerResults->accept(iHLT);
       for (unsigned int i = 0; i < trigNames_.size(); ++i) {
         if (triggerNames_[iHLT].find(trigNames_[i]) != std::string::npos) {
           if (hlt > 0) {
             ok = true;
           }
 #ifdef DebugLog
           edm::LogInfo("HEMuon") << "AlCaHEMuonFilter::Trigger " << triggerNames_[iHLT] << " Flag " << hlt << ":" << ok
                                  << std::endl;
 #endif
         }
       }
     }
     if (ok) {
       //Step2: Get geometry/B-field information
       //Get magnetic field
       edm::ESHandle<MagneticField> bFieldH;
       iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
       const MagneticField* bField = bFieldH.product();
       // get handles to calogeometry
       edm::ESHandle<CaloGeometry> pG;
       iSetup.get<CaloGeometryRecord>().get(pG);
       const CaloGeometry* geo = pG.product();
 
       // Relevant blocks from iEvent
       edm::Handle<reco::MuonCollection> _Muon;
       iEvent.getByToken(tok_Muon_, _Muon);
 #ifdef DebugLog
       edm::LogInfo("HEMuon") << "AlCaHEMuonFilter::Muon Handle " << _Muon.isValid() << std::endl;
 #endif
       if (_Muon.isValid()) {
         for (reco::MuonCollection::const_iterator RecMuon = _Muon->begin(); RecMuon != _Muon->end(); ++RecMuon) {
 #ifdef DebugLog
           edm::LogInfo("HEMuon") << "AlCaHEMuonFilter::Muon:Track " << RecMuon->track().isNonnull() << " innerTrack "
                                  << RecMuon->innerTrack().isNonnull() << " outerTrack "
                                  << RecMuon->outerTrack().isNonnull() << " globalTrack "
                                  << RecMuon->globalTrack().isNonnull() << std::endl;
 #endif
           if ((RecMuon->track().isNonnull()) && (RecMuon->innerTrack().isNonnull()) &&
               (RecMuon->outerTrack().isNonnull()) && (RecMuon->globalTrack().isNonnull())) {
             bool ptCut = RecMuon->pt() < muonptCut_;
             bool etaCut = fabs(RecMuon->eta()) < muonetaCut_;
             if (ptCut || etaCut)
               continue;
 
             const reco::Track* pTrack = (RecMuon->innerTrack()).get();
             spr::propagatedTrackID trackID = spr::propagateCALO(pTrack, geo, bField, false);
 #ifdef DebugLog
             edm::LogInfo("HEMuon") << "AlCaHEMuonFilter::Propagate: ECAL " << trackID.okECAL << " to HCAL "
                                    << trackID.okHCAL << std::endl;
 #endif
             double trackIso = RecMuon->isolationR03().sumPt;
             double caloIso = RecMuon->isolationR03().emEt + RecMuon->isolationR03().hadEt;
             double isolR04 =
                 ((RecMuon->pfIsolationR04().sumChargedHadronPt +
                   std::max(0.,
                            RecMuon->pfIsolationR04().sumNeutralHadronEt + RecMuon->pfIsolationR04().sumPhotonEt -
                                (0.5 * RecMuon->pfIsolationR04().sumPUPt))) /
                  RecMuon->pt());
             bool isoCut = (pfCut_) ? (isolR04 < pfIsoCut_) : ((trackIso < trackIsoCut_) && (caloIso < caloIsoCut_));
             if ((trackID.okECAL) && (trackID.okHCAL) && isoCut) {
               accept = true;
               break;
             }
           }
         }
       }
     }
   }
   // Step 4:  Return the acceptance flag
   if (accept) {
     ++nGood_;
     if (((nGood_ - 1) % preScale_) != 0) {
       accept = false;
     } else {
       ++nFinal_;
     }
   }
   return accept;
 
 }  // AlCaHEMuonFilter::filter

void AlCaHEMuonFilter::globalEndJob ( const AlCaHEMuons::Counters * counters )

static

Definition at line 233 of file AlCaHEMuonFilter.cc.

References AlCaHEMuons::Counters::nAll_, AlCaHEMuons::Counters::nFinal_, and AlCaHEMuons::Counters::nGood_.

                                                                     {
   edm::LogInfo("HEMuon") << "Selects " << count->nFinal_ << " out of " << count->nGood_ << " good events out of "
                          << count->nAll_ << " total # of events\n";
 }

static std::unique_ptr<AlCaHEMuons::Counters> AlCaHEMuonFilter::initializeGlobalCache ( edm::ParameterSet const & )

inlinestatic

Definition at line 54 of file AlCaHEMuonFilter.cc.

References Ecal2004TBTDCRanges_v1_cff::endRun, fillDescriptions(), and ALCARECOTkAlBeamHalo_cff::filter.

                                                                                           {
     return std::make_unique<AlCaHEMuons::Counters>();
   }