#include <L1Trigger/L1TNtuples/src/L1MetFilterRecoTreeProducer.cc>

Inheritance diagram for L1MetFilterRecoTreeProducer:

Public Member Functions
	L1MetFilterRecoTreeProducer (const edm::ParameterSet &)

	~L1MetFilterRecoTreeProducer () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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)

Public Attributes
L1Analysis::L1AnalysisRecoMetFilterDataFormat *	metFilter_data

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob (void) override

void	doMetFilters (edm::Handle< edm::TriggerResults > trigRes, edm::TriggerNames trigNames)

void	endJob () override

Private Attributes
bool	badChCandFilterResultMissing_

edm::EDGetTokenT< bool >	badChCandFilterToken_

bool	badChCandRes = false

bool	badPFMuonFilterResultMissing_

edm::EDGetTokenT< bool >	badPFMuonFilterToken_

bool	badPFMuonRes = false

edm::Service< TFileService >	fs_

TTree *	tree_

bool	triggerResultsMissing_

edm::EDGetTokenT< edm::TriggerResults >	triggerResultsToken_

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

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

Description: Produces tree containing reco quantities

Definition at line 50 of file L1MetFilterRecoTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 94 of file L1MetFilterRecoTreeProducer.cc.

References badChCandFilterToken_, badPFMuonFilterToken_, fs_, edm::ParameterSet::getUntrackedParameter(), TFileService::make(), metFilter_data, tree_, and triggerResultsToken_.

                                                                                       :
   triggerResultsMissing_(false)
   //hbheNoiseFilterResultMissing_(false)
 {
 
   triggerResultsToken_ = consumes<edm::TriggerResults>(iConfig.getUntrackedParameter("triggerResultsToken",edm::InputTag("TriggerResults")));
   
   //hbheNoiseFilterResultToken_ = consumes<bool>(iConfig.getUntrackedParameter("hbheNoiseFilterResultToken",edm::InputTag("HBHENoiseFilterResultProducer:HBHENoiseFilterResult")));
 
   badChCandFilterToken_ = consumes<bool>(iConfig.getUntrackedParameter<edm::InputTag>("badChargedCandidateFilterToken"));
 
   badPFMuonFilterToken_ = consumes<bool>(iConfig.getUntrackedParameter<edm::InputTag>("badPFMuonFilterToken"));
 
   metFilter_data = new L1Analysis::L1AnalysisRecoMetFilterDataFormat();
 
   // set up output
   tree_=fs_->make<TTree>("MetFilterRecoTree", "MetFilterRecoTree");
   tree_->Branch("MetFilters", "L1Analysis::L1AnalysisRecoMetFilterDataFormat", &metFilter_data, 32000, 3);
 
 }

L1MetFilterRecoTreeProducer::~L1MetFilterRecoTreeProducer ( )

override

Definition at line 116 of file L1MetFilterRecoTreeProducer.cc.

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

Member Function Documentation

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

overrideprivate

Definition at line 130 of file L1MetFilterRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoMetFilterDataFormat::badChCandFilter, badChCandFilterResultMissing_, badChCandFilterToken_, badChCandRes, L1Analysis::L1AnalysisRecoMetFilterDataFormat::badPFMuonFilter, badPFMuonFilterResultMissing_, badPFMuonFilterToken_, badPFMuonRes, doMetFilters(), edm::Event::getByToken(), edm::HandleBase::isValid(), metFilter_data, L1Analysis::L1AnalysisRecoMetFilterDataFormat::Reset(), tree_, edm::Event::triggerNames(), triggerResultsMissing_, triggerResultsToken_, and trigNames.

 {
 
   metFilter_data->Reset();
 
   //hbheNFRes = false;
   badPFMuonRes = false;
   badChCandRes = false;
 
 
   // get trigger results
   edm::Handle<edm::TriggerResults> trigRes;
   iEvent.getByToken(triggerResultsToken_, trigRes);
   
   if (trigRes.isValid()) {
     // get trigger names
     const edm::TriggerNames& trigNames = iEvent.triggerNames(*trigRes);
     doMetFilters(trigRes, trigNames);
   }
   else {
     if (!triggerResultsMissing_) {edm::LogWarning("MissingProduct") << "Trigger results not found.  Some MET filters will not be filled" << std::endl;}
     triggerResultsMissing_ = true;
   }
   
   // // get hbhe noise filter result
   // edm::Handle<bool> hbheNoiseFilterResult;
   // iEvent.getByToken(hbheNoiseFilterResultToken_, hbheNoiseFilterResult);
   
   // if(hbheNoiseFilterResult.isValid()){
   //   hbheNFRes = *hbheNoiseFilterResult;
   // }
   // else {
   //   if(!hbheNoiseFilterResultMissing_){edm::LogWarning("MissingProduct") << "HBHE Noise Filter Result not found.  Branch will not be filled" << std::endl;}
   //   hbheNoiseFilterResultMissing_=true;
   // }
 
   edm::Handle<bool> badChCandFilterResult;
   iEvent.getByToken(badChCandFilterToken_, badChCandFilterResult);
   
   if(badChCandFilterResult.isValid()){
     badChCandRes = *badChCandFilterResult;
   }
   else {
     if(!badChCandFilterResultMissing_){edm::LogWarning("MissingProduct") << "Bad Charged Hadron Candidate Filter Result not found.  Branch will not be filled" << std::endl;}
     badChCandFilterResultMissing_=true;
   }
 
   edm::Handle<bool> badPFMuonFilterResult;
   iEvent.getByToken(badPFMuonFilterToken_, badPFMuonFilterResult);
 
   if(badPFMuonFilterResult.isValid()){
     badPFMuonRes = *badPFMuonFilterResult;
   }
   else {
     if(!badPFMuonFilterResultMissing_){edm::LogWarning("MissingProduct") << "Bad PF Muon Filter Result not found.  Branch will not be filled" << std::endl;}
     badPFMuonFilterResultMissing_=true;
   }
   
   //metFilter_data->hbheNoiseFilter                 = hbheNFRes;
   metFilter_data->badChCandFilter                 = badChCandRes;
   metFilter_data->badPFMuonFilter                 = badPFMuonRes;
 
   tree_->Fill();
 }

void L1MetFilterRecoTreeProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 231 of file L1MetFilterRecoTreeProducer.cc.

232 {

233 }

void L1MetFilterRecoTreeProducer::doMetFilters	(	edm::Handle< edm::TriggerResults >	trigRes,
		edm::TriggerNames	trigNames
	)

private

Definition at line 199 of file L1MetFilterRecoTreeProducer.cc.

Referenced by analyze().

                                                                                                               {
 
   //get array size
   unsigned int numTrigs = trigNames.triggerNames().size();
 
   //get indices of flags from event parameter set
   unsigned int hbheNoiseFilterIndex                = trigNames.triggerIndex("Flag_HBHENoiseFilter");
   unsigned int hbheNoiseIsoFilterIndex             = trigNames.triggerIndex("Flag_HBHENoiseIsoFilter");
   unsigned int cscTightHalo2015FilterIndex         = trigNames.triggerIndex("Flag_CSCTightHalo2015Filter");
   unsigned int globalSuperTightHalo2016FilterIndex = trigNames.triggerIndex("Flag_globalSuperTightHalo2016Filter");
   unsigned int ecalDeadCellTPFilterIndex           = trigNames.triggerIndex("Flag_EcalDeadCellTriggerPrimitiveFilter");
   unsigned int goodVerticesFilterIndex             = trigNames.triggerIndex("Flag_goodVertices");
   unsigned int eeBadScFilterIndex                  = trigNames.triggerIndex("Flag_eeBadScFilter");
   unsigned int chHadTrackResFilterIndex            = trigNames.triggerIndex("Flag_chargedHadronTrackResolutionFilter");
   unsigned int muonBadTrackFilterIndex             = trigNames.triggerIndex("Flag_muonBadTrackFilter");                 
 
   //set flag
   metFilter_data->hbheNoiseFilter                 = hbheNoiseFilterIndex                  <  numTrigs ? trigRes->accept(hbheNoiseFilterIndex)                : false ;     
   metFilter_data->hbheNoiseIsoFilter              = hbheNoiseIsoFilterIndex               <  numTrigs ? trigRes->accept(hbheNoiseIsoFilterIndex)             : false ;     
   metFilter_data->cscTightHalo2015Filter          = cscTightHalo2015FilterIndex           <  numTrigs ? trigRes->accept(cscTightHalo2015FilterIndex)         : false ; 
   metFilter_data->globalSuperTightHalo2016Filter  = globalSuperTightHalo2016FilterIndex   <  numTrigs ? trigRes->accept(globalSuperTightHalo2016FilterIndex) : false ; 
   metFilter_data->ecalDeadCellTPFilter            = ecalDeadCellTPFilterIndex             <  numTrigs ? trigRes->accept(ecalDeadCellTPFilterIndex)           : false ;  
   metFilter_data->goodVerticesFilter              = goodVerticesFilterIndex               <  numTrigs ? trigRes->accept(goodVerticesFilterIndex)             : false ;  
   metFilter_data->eeBadScFilter                   = eeBadScFilterIndex                <  numTrigs ? trigRes->accept(eeBadScFilterIndex)                  : false ;  
   metFilter_data->chHadTrackResFilter             = chHadTrackResFilterIndex              <  numTrigs ? trigRes->accept(chHadTrackResFilterIndex)            : false ;  
   metFilter_data->muonBadTrackFilter              = muonBadTrackFilterIndex               <  numTrigs ? trigRes->accept(muonBadTrackFilterIndex)             : false ;  
   
 }