HLTDoublet< T1, T2 > Class Template Reference

#include <HLTDoublet.h>

Inheritance diagram for HLTDoublet< T1, T2 >:

Public Member Functions
	HLTDoublet (const edm::ParameterSet &)
virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
	~HLTDoublet ()
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 Types
typedef std::vector< T1 >	T1Collection
typedef edm::Ref< T1Collection >	T1Ref
typedef std::vector< T2 >	T2Collection
typedef edm::Ref< T2Collection >	T2Ref

Private Attributes
const bool	cutdelr_
const bool	cutdeta_
const bool	cutdphi_
const bool	cutminv_
const bool	cutpt_
const edm::InputTag	inputTag1_
const edm::InputTag	inputTag2_
const edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	inputToken1_
const edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	inputToken2_
const double	max_DelR_
const double	max_Deta_
const double	max_Dphi_
const double	max_Minv_
const double	max_Pt_
const double	min_DelR_
const double	min_Deta_
const double	min_Dphi_
const double	min_Minv_
const int	min_N_
const double	min_Pt_
const std::vector< edm::InputTag >	originTag1_
const std::vector< edm::InputTag >	originTag2_
const bool	same_
const int	triggerType1_
const int	triggerType2_

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 T1, typename T2>
class HLTDoublet< T1, T2 >

This class is an HLTFilter (-> EDFilter) implementing a basic HLT trigger for pairs of object, evaluating all pairs with the first object from collection 1, and the second object from collection 2, cutting on variables relating to their 4-momentum representations. The object collections are assumed to be outputs of HLTSinglet single-object-type filters so that the access is thorugh RefToBases and polymorphic.

Author

Martin Grunewald

See header file for documentation

Author

Martin Grunewald

Definition at line 34 of file HLTDoublet.h.

Member Typedef Documentation

template<typename T1 , typename T2 >

typedef std::vector<T1> HLTDoublet< T1, T2 >::T1Collection

private

Definition at line 65 of file HLTDoublet.h.

template<typename T1 , typename T2 >

typedef edm::Ref<T1Collection> HLTDoublet< T1, T2 >::T1Ref

private

Definition at line 66 of file HLTDoublet.h.

template<typename T1 , typename T2 >

typedef std::vector<T2> HLTDoublet< T1, T2 >::T2Collection

private

Definition at line 67 of file HLTDoublet.h.

template<typename T1 , typename T2 >

typedef edm::Ref<T2Collection> HLTDoublet< T1, T2 >::T2Ref

private

Definition at line 68 of file HLTDoublet.h.

Constructor & Destructor Documentation

template<typename T1 , typename T2 >

HLTDoublet< T1, T2 >::HLTDoublet ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 29 of file HLTDoublet.cc.

References HLTDoublet< T1, T2 >::cutdelr_, HLTDoublet< T1, T2 >::cutdeta_, HLTDoublet< T1, T2 >::cutdphi_, HLTDoublet< T1, T2 >::cutminv_, HLTDoublet< T1, T2 >::cutpt_, edm::InputTag::encode(), HLTDoublet< T1, T2 >::inputTag1_, HLTDoublet< T1, T2 >::inputTag2_, LogDebug, HLTDoublet< T1, T2 >::max_DelR_, HLTDoublet< T1, T2 >::max_Deta_, HLTDoublet< T1, T2 >::max_Dphi_, HLTDoublet< T1, T2 >::max_Minv_, HLTDoublet< T1, T2 >::max_Pt_, HLTDoublet< T1, T2 >::min_DelR_, HLTDoublet< T1, T2 >::min_Deta_, HLTDoublet< T1, T2 >::min_Dphi_, HLTDoublet< T1, T2 >::min_Minv_, HLTDoublet< T1, T2 >::min_N_, HLTDoublet< T1, T2 >::min_Pt_, HLTDoublet< T1, T2 >::same_, HLTDoublet< T1, T2 >::triggerType1_, and HLTDoublet< T1, T2 >::triggerType2_.

                                                            : HLTFilter(iConfig),
  originTag1_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag1")),
  originTag2_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag2")),
  inputTag1_(iConfig.template getParameter<edm::InputTag>("inputTag1")),
  inputTag2_(iConfig.template getParameter<edm::InputTag>("inputTag2")),
  inputToken1_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag1_)),
  inputToken2_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag2_)),
  triggerType1_(iConfig.template getParameter<int>("triggerType1")),
  triggerType2_(iConfig.template getParameter<int>("triggerType2")),
  min_Dphi_ (iConfig.template getParameter<double>("MinDphi")),
  max_Dphi_ (iConfig.template getParameter<double>("MaxDphi")),
  min_Deta_ (iConfig.template getParameter<double>("MinDeta")),
  max_Deta_ (iConfig.template getParameter<double>("MaxDeta")),
  min_Minv_ (iConfig.template getParameter<double>("MinMinv")),
  max_Minv_ (iConfig.template getParameter<double>("MaxMinv")),
  min_DelR_ (iConfig.template getParameter<double>("MinDelR")),
  max_DelR_ (iConfig.template getParameter<double>("MaxDelR")),
  min_Pt_   (iConfig.template getParameter<double>("MinPt")),
  max_Pt_   (iConfig.template getParameter<double>("MaxPt")),
  min_N_    (iConfig.template getParameter<int>("MinN")),
  same_     (inputTag1_.encode() == inputTag2_.encode()),    // same collections to be compared?
  cutdphi_  (min_Dphi_ <= max_Dphi_), // cut active?
  cutdeta_  (min_Deta_ <= max_Deta_), // cut active?
  cutminv_  (min_Minv_ <= max_Minv_), // cut active?
  cutdelr_  (min_DelR_ <= max_DelR_), // cut active?
  cutpt_    (min_Pt_   <= max_Pt_  )  // cut active?
{
   LogDebug("") << "InputTags and cuts : "
        << inputTag1_.encode() << " " << inputTag2_.encode()
        << triggerType1_ << " " << triggerType2_
        << " Dphi [" << min_Dphi_ << " " << max_Dphi_ << "]"
                << " Deta [" << min_Deta_ << " " << max_Deta_ << "]"
                << " Minv [" << min_Minv_ << " " << max_Minv_ << "]"
                << " DelR [" << min_DelR_ << " " << max_DelR_ << "]"
                << " Pt   [" << min_Pt_   << " " << max_Pt_   << "]"
                << " MinN =" << min_N_
        << " same/dphi/deta/minv/delr/pt "
        << same_
        << cutdphi_ << cutdeta_ << cutminv_ << cutdelr_ << cutpt_;
}

template<typename T1 , typename T2 >

HLTDoublet< T1, T2 >::~HLTDoublet

(

)

Definition at line 71 of file HLTDoublet.cc.

{
}

Member Function Documentation

template<typename T1 , typename T2 >

void HLTDoublet< T1, T2 >::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 76 of file HLTDoublet.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), defaultModuleLabel(), HLT_25ns10e33_v2_cff::InputTag, HLT_25ns10e33_v2_cff::originTag1, and HLT_25ns10e33_v2_cff::originTag2.

                                                                              {
  edm::ParameterSetDescription desc;
  makeHLTFilterDescription(desc);
  std::vector<edm::InputTag> originTag1(1,edm::InputTag("hltOriginal1"));
  std::vector<edm::InputTag> originTag2(1,edm::InputTag("hltOriginal2"));
  desc.add<std::vector<edm::InputTag> >("originTag1",originTag1);
  desc.add<std::vector<edm::InputTag> >("originTag2",originTag2);
  desc.add<edm::InputTag>("inputTag1",edm::InputTag("hltFiltered1"));
  desc.add<edm::InputTag>("inputTag2",edm::InputTag("hltFiltered22"));
  desc.add<int>("triggerType1",0);
  desc.add<int>("triggerType2",0);
  desc.add<double>("MinDphi",+1.0);
  desc.add<double>("MaxDphi",-1.0);
  desc.add<double>("MinDeta",+1.0);
  desc.add<double>("MaxDeta",-1.0);
  desc.add<double>("MinMinv",+1.0);
  desc.add<double>("MaxMinv",-1.0);
  desc.add<double>("MinDelR",+1.0);
  desc.add<double>("MaxDelR",-1.0);
  desc.add<double>("MinPt"  ,+1.0);
  desc.add<double>("MaxPt"  ,-1.0);
  desc.add<int>("MinN",1);
  descriptions.add(defaultModuleLabel<HLTDoublet<T1,T2>>(), desc);
}

template<typename T1 , typename T2 >

bool HLTDoublet< T1, T2 >::hltFilter ( edm::Event &  iEvent, const edm::EventSetup &  iSetup, trigger::TriggerFilterObjectWithRefs &  filterproduct  ) const

overridevirtual

Implements HLTFilter.

Definition at line 108 of file HLTDoublet.cc.

References funct::abs(), accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), edm::InputTag::encode(), edm::Event::getByToken(), edm::Event::getProvenance(), Exhume::I, i, HLT_25ns10e33_v2_cff::InputTag, instance, diffTwoXMLs::label, M_PI, edm::Provenance::moduleLabel(), gen::n, AlCaHLTBitMon_ParallelJobs::p, p1, p2, sysUtil::pid, LaserDQM_cfg::process, edm::Provenance::processName(), edm::Provenance::productInstanceName(), reco::tau::disc::Pt(), diffTwoXMLs::r1, diffTwoXMLs::r2, dt_dqm_sourceclient_common_cff::reco, CommPDSkim_cfg::saveTags, and mathSSE::sqrt().

{
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;

   // 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.

   bool accept(false);

   LogVerbatim("HLTDoublet") << " XXX " << moduleLabel() << " 0 " << std::endl;

   std::vector<T1Ref> coll1;
   std::vector<T2Ref> coll2;

   // get hold of pre-filtered object collections
   Handle<TriggerFilterObjectWithRefs> handle1, handle2;
   if (iEvent.getByToken(inputToken1_,handle1) && iEvent.getByToken(inputToken2_,handle2)) {
     handle1->getObjects(triggerType1_, coll1);
     handle2->getObjects(triggerType2_, coll2);
     const size_type n1(coll1.size());
     const size_type n2(coll2.size());

     if (saveTags()) {
       InputTag tagOld;
       for (unsigned int i=0; i<originTag1_.size(); ++i) {
     filterproduct.addCollectionTag(originTag1_[i]);
     LogVerbatim("HLTDoublet") << " XXX " << moduleLabel() << " 1a/" << i << " " << originTag1_[i].encode() << std::endl;
       }
       tagOld=InputTag();
       for (size_type i1=0; i1!=n1; ++i1) {
     const ProductID pid(coll1[i1].id());
     const string&    label(iEvent.getProvenance(pid).moduleLabel());
     const string& instance(iEvent.getProvenance(pid).productInstanceName());
     const string&  process(iEvent.getProvenance(pid).processName());
     InputTag tagNew(InputTag(label,instance,process));
     if (tagOld.encode()!=tagNew.encode()) {
       filterproduct.addCollectionTag(tagNew);
       tagOld=tagNew;
       LogVerbatim("HLTDoublet") << " XXX " << moduleLabel() << " 1b " << tagNew.encode() << std::endl;
     }
       }
       for (unsigned int i=0; i<originTag2_.size(); ++i) {
     filterproduct.addCollectionTag(originTag2_[i]);
     LogVerbatim("HLTDoublet") << " XXX " << moduleLabel() << " 2a/" << i << " " << originTag2_[i].encode() << std::endl;
       }
       tagOld=InputTag();
       for (size_type i2=0; i2!=n2; ++i2) {
     const ProductID pid(coll2[i2].id());
     const string&    label(iEvent.getProvenance(pid).moduleLabel());
     const string& instance(iEvent.getProvenance(pid).productInstanceName());
     const string&  process(iEvent.getProvenance(pid).processName());
     InputTag tagNew(InputTag(label,instance,process));
     if (tagOld.encode()!=tagNew.encode()) {
       filterproduct.addCollectionTag(tagNew);
       tagOld=tagNew;
       LogVerbatim("HLTDoublet") << " XXX " << moduleLabel() << " 2b " << tagNew.encode() << std::endl;
     }
       }
     }

     int n(0);
     T1Ref r1;
     T2Ref r2;
     Particle::LorentzVector p1,p2,p;
     for (unsigned int i1=0; i1!=n1; i1++) {
       r1=coll1[i1];
       p1=r1->p4();
       unsigned int I(0);
       if (same_) {I=i1+1;}
       for (unsigned int i2=I; i2!=n2; i2++) {
     r2=coll2[i2];
     p2=r2->p4();

     double Dphi(std::abs(p1.phi()-p2.phi()));
     if (Dphi>M_PI) Dphi=2.0*M_PI-Dphi;
    
     double Deta(std::abs(p1.eta()-p2.eta()));
    
     p=p1+p2;
     double Minv(std::abs(p.mass()));
     double Pt(p.pt());

     double DelR(sqrt(Dphi*Dphi+Deta*Deta));
    
     if ( ( (!cutdphi_) || ((min_Dphi_<=Dphi) && (Dphi<=max_Dphi_)) ) &&
          ( (!cutdeta_) || ((min_Deta_<=Deta) && (Deta<=max_Deta_)) ) &&
          ( (!cutminv_) || ((min_Minv_<=Minv) && (Minv<=max_Minv_)) ) &&
              ( (!cutdelr_) || ((min_DelR_<=DelR) && (DelR<=max_DelR_)) ) &&
              ( (!cutpt_  ) || ((min_Pt_  <=Pt  ) && (Pt  <=max_Pt_  )) ) ) {
       n++;
       filterproduct.addObject(triggerType1_,r1);
       filterproduct.addObject(triggerType2_,r2);
     }
    
       }
     }
     // filter decision
     accept = (n>=min_N_);
   }

   return accept;
}

Member Data Documentation

template<typename T1 , typename T2 >

const bool HLTDoublet< T1, T2 >::cutdelr_

private

Definition at line 62 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const bool HLTDoublet< T1, T2 >::cutdeta_

private

Definition at line 62 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const bool HLTDoublet< T1, T2 >::cutdphi_

private

Definition at line 62 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const bool HLTDoublet< T1, T2 >::cutminv_

private

Definition at line 62 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const bool HLTDoublet< T1, T2 >::cutpt_

private

Definition at line 62 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const edm::InputTag HLTDoublet< T1, T2 >::inputTag1_

private

Definition at line 47 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const edm::InputTag HLTDoublet< T1, T2 >::inputTag2_

private

Definition at line 48 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTDoublet< T1, T2 >::inputToken1_

private

Definition at line 49 of file HLTDoublet.h.

template<typename T1 , typename T2 >

const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTDoublet< T1, T2 >::inputToken2_

private

Definition at line 50 of file HLTDoublet.h.

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::max_DelR_

private

Definition at line 56 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::max_Deta_

private

Definition at line 54 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::max_Dphi_

private

Definition at line 53 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::max_Minv_

private

Definition at line 55 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::max_Pt_

private

Definition at line 57 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::min_DelR_

private

Definition at line 56 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::min_Deta_

private

Definition at line 54 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::min_Dphi_

private

Definition at line 53 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::min_Minv_

private

Definition at line 55 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const int HLTDoublet< T1, T2 >::min_N_

private

Definition at line 58 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const double HLTDoublet< T1, T2 >::min_Pt_

private

Definition at line 57 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const std::vector<edm::InputTag> HLTDoublet< T1, T2 >::originTag1_

private

Definition at line 45 of file HLTDoublet.h.

template<typename T1 , typename T2 >

const std::vector<edm::InputTag> HLTDoublet< T1, T2 >::originTag2_

private

Definition at line 46 of file HLTDoublet.h.

template<typename T1 , typename T2 >

const bool HLTDoublet< T1, T2 >::same_

private

Definition at line 61 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const int HLTDoublet< T1, T2 >::triggerType1_

private

Definition at line 51 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().

template<typename T1 , typename T2 >

const int HLTDoublet< T1, T2 >::triggerType2_

private

Definition at line 52 of file HLTDoublet.h.

Referenced by HLTDoublet< T1, T2 >::HLTDoublet().