HeavyChHiggsToTauNuSkim Class Reference

#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
std::vector< bool > const &	recordProvenanceList () 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)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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)

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	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (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::HeavyChHiggsToTauNuSkim ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 25 of file HeavyChHiggsToTauNuSkim.cc.

                                                                               {
  // 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;
}

References debug, edm::ParameterSet::getParameter(), MssmHbbBtagTriggerMonitor_cfi::jetEtaMax, heavyChHiggsToTauNu_Filter_cfi::jetEtaMin, heavyChHiggsToTauNu_Filter_cfi::jetEtMin, unpackData-CaloStage1::jetToken, heavyChHiggsToTauNu_Filter_cfi::minDRFromTau, and nEvents.

~HeavyChHiggsToTauNuSkim()

HeavyChHiggsToTauNuSkim::~HeavyChHiggsToTauNuSkim ( )

override

Definition at line 41 of file HeavyChHiggsToTauNuSkim.cc.

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

References nEvents.

Member Function Documentation

deltaPhi()

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

inlineprivate

Definition at line 40 of file HeavyChHiggsToTauNuSkim.h.

                                            {
    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;
  }

References PI.

Referenced by deltaR().

deltaR()

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

inlineprivate

Definition at line 57 of file HeavyChHiggsToTauNuSkim.h.

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

References deltaPhi(), HLT_FULL_cff::eta1, HLT_FULL_cff::eta2, and mathSSE::sqrt().

filter()

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

overridevirtual

Implements edm::EDFilter.

Definition at line 47 of file HeavyChHiggsToTauNuSkim.cc.

                                                                                  {
  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;
}

References PbPb_ZMuSkimMuonDPG_cff::deltaR, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), mps_fire::i, iEvent, edm::HandleBase::isValid(), MssmHbbBtagTriggerMonitor_cfi::jetEtaMax, heavyChHiggsToTauNu_Filter_cfi::jetEtaMin, heavyChHiggsToTauNu_Filter_cfi::jetEtMin, singleTopDQM_cfi::jets, unpackData-CaloStage1::jetToken, heavyionUCCDQM_cfi::maxEt, heavyChHiggsToTauNu_Filter_cfi::minDRFromTau, nEvents, reco::LeafCandidate::phi(), and edm::Handle< T >::product().

Member Data Documentation

debug

bool HeavyChHiggsToTauNuSkim::debug

private

Definition at line 63 of file HeavyChHiggsToTauNuSkim.h.

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

hltTauToken

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

private

Definition at line 65 of file HeavyChHiggsToTauNuSkim.h.

jetEtaMax

double HeavyChHiggsToTauNuSkim::jetEtaMax

private

Definition at line 70 of file HeavyChHiggsToTauNuSkim.h.

jetEtaMin

double HeavyChHiggsToTauNuSkim::jetEtaMin

private

Definition at line 69 of file HeavyChHiggsToTauNuSkim.h.

jetEtMin

double HeavyChHiggsToTauNuSkim::jetEtMin

private

Definition at line 68 of file HeavyChHiggsToTauNuSkim.h.

jetToken

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

private

Definition at line 66 of file HeavyChHiggsToTauNuSkim.h.

minDRFromTau

double HeavyChHiggsToTauNuSkim::minDRFromTau

private

Definition at line 71 of file HeavyChHiggsToTauNuSkim.h.

minNumberOfjets

int HeavyChHiggsToTauNuSkim::minNumberOfjets

private

Definition at line 67 of file HeavyChHiggsToTauNuSkim.h.

nEvents

int HeavyChHiggsToTauNuSkim::nEvents

private

Definition at line 73 of file HeavyChHiggsToTauNuSkim.h.

Referenced by looper.Looper::loop().

nSelectedEvents

int HeavyChHiggsToTauNuSkim::nSelectedEvents

private

Definition at line 73 of file HeavyChHiggsToTauNuSkim.h.