#include <HLTGlobalSums.h>

Inheritance diagram for HLTGlobalSums< T >:

Public Member Functions
virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override

	HLTGlobalSums (const edm::ParameterSet &)

	~HLTGlobalSums ()

Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)

int	module (edm::Event const &) const

const std::string *	moduleLabel () const

int	path (edm::Event const &) const

const std::string *	pathName (edm::Event const &) const

std::pair< int, int >	pmid (edm::Event const &) const

bool	saveTags () const

virtual	~HLTFilter ()

Public Member Functions inherited from edm::global::EDFilter<>
	EDFilter ()=default

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

ModuleDescription const &	moduleDescription () const

virtual	~EDFilterBase ()

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

	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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 ()

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

Static Public Member Functions inherited from HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
edm::InputTag	inputTag_

edm::EDGetTokenT< std::vector < T > >	inputToken_

double	max_

double	min_

int	min_N_

std::string	observable_

int	tid_

int	triggerType_

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	ModuleType

Public Types inherited from edm::ProducerBase
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 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

template<typename T>
class HLTGlobalSums< T >

This class is an HLTFilter (-> EDFilter) implementing cuts on global sums such as the scalar sum of Et (a.k.a. H_T), available in the T=CaloMET or T=MET object.

Author: Martin Grunewald

Definition at line 26 of file HLTGlobalSums.h.

Constructor & Destructor Documentation

template<typename T >

HLTGlobalSums< T >::HLTGlobalSums ( const edm::ParameterSet & iConfig )

explicit

Definition at line 28 of file HLTGlobalSums.cc.

References edm::InputTag::encode(), HLTGlobalSums< T >::inputTag_, LogDebug, HLTGlobalSums< T >::max_, HLTGlobalSums< T >::min_, HLTGlobalSums< T >::min_N_, HLTGlobalSums< T >::observable_, HLTGlobalSums< T >::tid_, trigger::TriggerELongit, trigger::TriggerHLongit, trigger::TriggerMETSig, trigger::TriggerMHTSig, trigger::TriggerTET, trigger::TriggerTHT, and HLTGlobalSums< T >::triggerType_.

                                                               : HLTFilter(iConfig),
   inputTag_   (iConfig.template getParameter<edm::InputTag>("inputTag")),
   inputToken_ (consumes<std::vector<T> >(inputTag_)),
   triggerType_(iConfig.template getParameter<int>("triggerType")),
   observable_ (iConfig.template getParameter<std::string>("observable")),
   min_        (iConfig.template getParameter<double>("Min")),
   max_        (iConfig.template getParameter<double>("Max")),
   min_N_      (iConfig.template getParameter<int>("MinN")),
   tid_(triggerType_)
 {
    LogDebug("") << "InputTags and cuts : "
         << inputTag_.encode() << " "
         << triggerType_ << " "
         << observable_
         << " Range [" << min_ << " " << max_ << "]"
                 << " MinN =" << min_N_ ;
 
    if (observable_=="sumEt") {
      tid_=triggerType_;
    } else if (observable_=="mEtSig") {
      if (triggerType_==trigger::TriggerTET) {
        tid_=trigger::TriggerMETSig;
      } else if (triggerType_==trigger::TriggerTHT) {
        tid_=trigger::TriggerMHTSig;
      } else {
        tid_=triggerType_;
      }
    } else if (observable_=="e_longitudinal") {
      if (triggerType_==trigger::TriggerTET) {
        tid_=trigger::TriggerELongit;
      } else if (triggerType_==trigger::TriggerTHT) {
        tid_=trigger::TriggerHLongit;
      } else {
        tid_=triggerType_;
      }
    } else {
      tid_=triggerType_;
    }
 }

template<typename T >

HLTGlobalSums< T >::~HLTGlobalSums ( )

Definition at line 69 of file HLTGlobalSums.cc.

70 {

71 }

Member Function Documentation

template<typename T >

void HLTGlobalSums< T >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 75 of file HLTGlobalSums.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), defaultModuleLabel(), HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                              {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("inputTag",edm::InputTag("hltCollection"));
   desc.add<int>("triggerType",0);
   desc.add<std::string>("observable","");
   desc.add<double>("Min",-1e125);
   desc.add<double>("Max",+1e125);
   desc.add<int>("MinN",1);
   descriptions.add(defaultModuleLabel<HLTGlobalSums<T>>(), desc);
 }

template<typename T >

bool HLTGlobalSums< T >::hltFilter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)		const

overridevirtual

Implements HLTFilter.

Definition at line 94 of file HLTGlobalSums.cc.

References funct::abs(), accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), HLT_FULL_cff::distance, edm::Event::getByToken(), LogDebug, gen::n, dt_dqm_sourceclient_common_cff::reco, CommPDSkim_cfg::saveTags, trigger::TriggerELongit, trigger::TriggerHLongit, trigger::TriggerMETSig, trigger::TriggerMHTSig, trigger::TriggerTET, trigger::TriggerTHT, and relativeConstraints::value.

 {
    using namespace std;
    using namespace edm;
    using namespace reco;
    using namespace trigger;
 
    typedef vector<T> TCollection;
    typedef Ref<TCollection> TRef;
 
    // All HLT filters must create and fill an HLT filter object,
    // recording any reconstructed physics objects satisfying (or not)
    // this HLT filter, and place it in the Event.
 
    // The filter object
    if (saveTags()) filterproduct.addCollectionTag(inputTag_);
    // Ref to Candidate object to be recorded in filter object
    TRef ref;
 
 
    // get hold of MET product from Event
    Handle<TCollection>   objects;
    iEvent.getByToken(inputToken_,objects);
    if (!objects.isValid()) {
      LogDebug("") << inputTag_ << " collection not found!";
      return false;
    }
 
    LogDebug("") << "Size of MET collection: " << objects->size();
    if (objects->size()==0) {
      LogDebug("") << "MET collection does not contain a MET object!";
    } else if (objects->size()>1) {
      LogDebug("") << "MET collection contains more than one MET object!";
    }
 
    int n(0);
    double value(0.0);
    typename TCollection::const_iterator ibegin(objects->begin());
    typename TCollection::const_iterator iend(objects->end());
    typename TCollection::const_iterator iter;
    for (iter=ibegin; iter!=iend; iter++) {
 
      // get hold of value of observable to cut on
      if ( (tid_==TriggerTET) || (tid_==TriggerTHT) ) {
        value=iter->sumEt();
      } else if ( (tid_==TriggerMETSig) || (tid_==TriggerMHTSig) ) {
        value=iter->mEtSig();
      } else if ( (tid_==TriggerELongit) || (tid_==TriggerHLongit) ) {
        value=iter->e_longitudinal();
      } else {
        value=0.0;
      }
 
      value=std::abs(value);
 
      if ( ( (min_<0.0) || (min_<=value) ) &&
       ( (max_<0.0) || (value<=max_) ) ) {
        n++;
        ref=TRef(objects,distance(ibegin,iter));
        filterproduct.addObject(tid_,ref);
      }
 
    }
 
    // filter decision
    const bool accept(n>=min_N_);
 
    return accept;
 }