HeavyChHiggsToTauNuSkim Class Reference

#include <HeavyChHiggsToTauNuSkim.h>

Inheritance diagram for HeavyChHiggsToTauNuSkim:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	HeavyChHiggsToTauNuSkim (const edm::ParameterSet &)
	~HeavyChHiggsToTauNuSkim ()
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
ModuleDescription const &	moduleDescription () const
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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(), 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

(

)

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  )

virtual

Implements edm::EDFilter.

Definition at line 54 of file HeavyChHiggsToTauNuSkim.cc.

References deltaR(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), edm::Event::getByToken(), i, edm::HandleBase::isValid(), 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(), and rrapi.RRApi::get().

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.