#include <HeavyChHiggsToTauNuSkim.h>

Inheritance diagram for HeavyChHiggsToTauNuSkim:

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

	HeavyChHiggsToTauNuSkim (const edm::ParameterSet &)

	~HeavyChHiggsToTauNuSkim () override

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
double	deltaPhi (double phi1, double phi2)

double	deltaR (double eta1, double eta2, double phi1, double phi2)

Private Attributes
bool	debug

edm::EDGetTokenT< reco::IsolatedTauTagInfoCollection >	hltTauToken

double	jetEtaMax

double	jetEtaMin

double	jetEtMin

edm::EDGetTokenT< reco::CaloJetCollection >	jetToken

double	minDRFromTau

int	minNumberOfjets

int	nEvents

int	nSelectedEvents

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	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

Static Public Member Functions inherited from edm::EDFilter
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 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

Filter to select events passing L1 single tau HLT tau+MET 3 offline jets

Author: Sami Lehti - HIP Helsinki

This class is an EDFilter for heavy H+->taunu events

Author: Sami Lehti - HIP Helsinki

Definition at line 32 of file HeavyChHiggsToTauNuSkim.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 26 of file HeavyChHiggsToTauNuSkim.cc.

References debug, edm::ParameterSet::getParameter(), mhtProducer_cfi::jetEtaMax, and nEvents.

                                                                                {
 
     // Local Debug flag
     debug           = iConfig.getParameter<bool>("DebugHeavyChHiggsToTauNuSkim");
 
     hltTauToken = consumes<IsolatedTauTagInfoCollection>(iConfig.getParameter<InputTag>("HLTTauCollection"));
     jetToken        = consumes<CaloJetCollection>(iConfig.getParameter<InputTag>("JetTagCollection"));
     minNumberOfjets = iConfig.getParameter<int>("minNumberOfJets");
     jetEtMin        = iConfig.getParameter<double>("jetEtMin");
     jetEtaMin       = iConfig.getParameter<double>("jetEtaMin");
     jetEtaMax       = iConfig.getParameter<double>("jetEtaMax");
     minDRFromTau    = iConfig.getParameter<double>("minDRFromTau");
 
         nEvents         = 0;
         nSelectedEvents = 0;
 
 }

HeavyChHiggsToTauNuSkim::~HeavyChHiggsToTauNuSkim ( )

override

Definition at line 45 of file HeavyChHiggsToTauNuSkim.cc.

References nEvents.

                                                  {
   edm::LogVerbatim("HeavyChHiggsToTauNuSkim")
   << " Number_events_read " << nEvents
   << " Number_events_kept " << nSelectedEvents
   << " Efficiency         " << ((double)nSelectedEvents)/((double) nEvents + 0.01) << std::endl;
 
 }

Member Function Documentation

double HeavyChHiggsToTauNuSkim::deltaPhi	(	double	phi1,
		double	phi2
	)

inlineprivate

Definition at line 41 of file HeavyChHiggsToTauNuSkim.h.

References PI.

Referenced by deltaR().

                                              {
         const double PI = 3.1415926535;
         // in ORCA phi = [0,2pi], in TLorentzVector phi = [-pi,pi].
         // With the conversion below deltaPhi works ok despite the
         // 2*pi difference in phi definitions.
         if(phi1 < 0) phi1 += 2*PI;
         if(phi2 < 0) phi2 += 2*PI;
 
         double dphi = fabs(phi1-phi2);
 
         if(dphi > PI) dphi = 2*PI - dphi;
         return dphi;
     }

double HeavyChHiggsToTauNuSkim::deltaR	(	double	eta1,
		double	eta2,
		double	phi1,
		double	phi2
	)

inlineprivate

Definition at line 55 of file HeavyChHiggsToTauNuSkim.h.

References deltaPhi(), and mathSSE::sqrt().

                                                                      {
         double dphi = deltaPhi(phi1,phi2);
         double deta = fabs(eta1-eta2);
         return sqrt(dphi*dphi + deta*deta);
     }

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

override

Definition at line 54 of file HeavyChHiggsToTauNuSkim.cc.

References deltaR(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), edm::Event::getByToken(), mps_fire::i, edm::HandleBase::isValid(), mhtProducer_cfi::jetEtaMax, fwrapper::jets, nEvents, reco::LeafCandidate::phi(), and edm::Handle< T >::product().

                                                                                   {
 
   nEvents++;
 
   Handle<IsolatedTauTagInfoCollection> tauTagL3Handle;
   iEvent.getByToken(hltTauToken, tauTagL3Handle);
 
   if ( !tauTagL3Handle.isValid() ) return false;
 
 
   Jet theTau;
   double maxEt = 0;
   if (tauTagL3Handle.isValid() ) {
     const IsolatedTauTagInfoCollection & L3Taus = *(tauTagL3Handle.product());
     IsolatedTauTagInfoCollection::const_iterator i;
     for ( i = L3Taus.begin(); i != L3Taus.end(); i++ ) {
       if (i->discriminator() == 0) continue;
       Jet taujet = *(i->jet().get());
       if (taujet.et() > maxEt) {
         maxEt = taujet.et();
         theTau = taujet;
       }
     }
   }
 
   if (maxEt == 0) return false;
 
   // jets
 
   Handle<CaloJetCollection> jetHandle;
   iEvent.getByToken(jetToken,jetHandle);
 
   if ( !jetHandle.isValid() ) return false;
 
   bool accepted = false;
 
   if (jetHandle.isValid() ) {
     int nJets = 0;
     const reco::CaloJetCollection & jets = *(jetHandle.product());
     CaloJetCollection::const_iterator iJet;
     for (iJet = jets.begin(); iJet!= jets.end(); iJet++ ) {
       if (iJet->et()  > jetEtMin  &&
           iJet->eta() > jetEtaMin &&
       iJet->eta() < jetEtaMax ) {
         double DR = deltaR(theTau.eta(),iJet->eta(),theTau.phi(),iJet->phi());
         if (DR > minDRFromTau) nJets++;
       }
     }
     if (nJets >= minNumberOfjets) {
       accepted = true;
       nSelectedEvents++;
     }
   }
   return accepted;
 }

Member Data Documentation

bool HeavyChHiggsToTauNuSkim::debug

private

Definition at line 61 of file HeavyChHiggsToTauNuSkim.h.

Referenced by rrapi.RRApi::dprint(), rrapi.RRApi::get(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), and runTauIdMVA.TauIDEmbedder::runTauID().

edm::EDGetTokenT<reco::IsolatedTauTagInfoCollection> HeavyChHiggsToTauNuSkim::hltTauToken

private

Definition at line 63 of file HeavyChHiggsToTauNuSkim.h.

double HeavyChHiggsToTauNuSkim::jetEtaMax

private

Definition at line 68 of file HeavyChHiggsToTauNuSkim.h.

double HeavyChHiggsToTauNuSkim::jetEtaMin

private

Definition at line 67 of file HeavyChHiggsToTauNuSkim.h.

double HeavyChHiggsToTauNuSkim::jetEtMin

private

Definition at line 66 of file HeavyChHiggsToTauNuSkim.h.

edm::EDGetTokenT<reco::CaloJetCollection> HeavyChHiggsToTauNuSkim::jetToken

private

Definition at line 64 of file HeavyChHiggsToTauNuSkim.h.

double HeavyChHiggsToTauNuSkim::minDRFromTau

private

Definition at line 69 of file HeavyChHiggsToTauNuSkim.h.

int HeavyChHiggsToTauNuSkim::minNumberOfjets

private

Definition at line 65 of file HeavyChHiggsToTauNuSkim.h.

int HeavyChHiggsToTauNuSkim::nEvents

private

Definition at line 71 of file HeavyChHiggsToTauNuSkim.h.

Referenced by looper.Looper::loop().

int HeavyChHiggsToTauNuSkim::nSelectedEvents

private

Definition at line 71 of file HeavyChHiggsToTauNuSkim.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation