CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
HLTDisplacedEgammaFilter Class Reference

#include <HLTDisplacedEgammaFilter.h>

Inheritance diagram for HLTDisplacedEgammaFilter:
HLTFilter edm::global::EDFilter<> edm::global::EDFilterBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 HLTDisplacedEgammaFilter (const edm::ParameterSet &)
 
bool hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
 
 ~HLTDisplacedEgammaFilter () override
 
- 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
 
 ~HLTFilter () override
 
- Public Member Functions inherited from edm::global::EDFilter<>
 EDFilter ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsStreamLuminosityBlocks () const final
 
bool wantsStreamRuns () const final
 
- Public Member Functions inherited from edm::global::EDFilterBase
 EDFilterBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDFilterBase () 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< ConsumesInfoconsumesInfo () 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
 
ESProxyIndex const * esGetTokenIndices (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 const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (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)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

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

bool EBOnly
 
edm::InputTag inputTag_
 
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefsinputToken_
 
edm::InputTag inputTrk
 
edm::EDGetTokenT< reco::TrackCollectioninputTrkToken_
 
edm::InputTag l1EGTag_
 
int maxTrkCut
 
int ncandcut_
 
edm::InputTag rechitsEB
 
edm::EDGetTokenT< EcalRecHitCollectionrechitsEBToken_
 
edm::InputTag rechitsEE
 
edm::EDGetTokenT< EcalRecHitCollectionrechitsEEToken_
 
double seedTimeMax
 
double seedTimeMin
 
double sMaj_max
 
double sMaj_min
 
double sMin_max
 
double sMin_min
 
double trkdRCut
 
double trkPtCut
 

Additional Inherited Members

- Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase 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
 
- 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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
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

Author
Monica Vazquez Acosta (CERN)

Definition at line 23 of file HLTDisplacedEgammaFilter.h.

Constructor & Destructor Documentation

HLTDisplacedEgammaFilter::HLTDisplacedEgammaFilter ( const edm::ParameterSet iConfig)
explicit

Definition at line 22 of file HLTDisplacedEgammaFilter.cc.

References EBOnly, edm::ParameterSet::getParameter(), inputTag_, inputToken_, inputTrk, inputTrkToken_, l1EGTag_, maxTrkCut, ncandcut_, rechitsEB, rechitsEBToken_, rechitsEE, rechitsEEToken_, seedTimeMax, seedTimeMin, sMaj_max, sMaj_min, sMin_max, sMin_min, trkdRCut, trkPtCut, and ~HLTDisplacedEgammaFilter().

22  : HLTFilter(iConfig)
23 {
24  inputTag_ = iConfig.getParameter< edm::InputTag > ("inputTag");
25  ncandcut_ = iConfig.getParameter<int> ("ncandcut");
26  l1EGTag_ = iConfig.getParameter< edm::InputTag > ("l1EGCand");
27 
28  inputTrk = iConfig.getParameter< edm::InputTag > ("inputTrack");
29  trkPtCut = iConfig.getParameter<double> ("trackPtCut");
30  trkdRCut = iConfig.getParameter<double> ("trackdRCut");
31  maxTrkCut = iConfig.getParameter<int> ("maxTrackCut");
32 
33  rechitsEB = iConfig.getParameter< edm::InputTag > ("RecHitsEB");
34  rechitsEE = iConfig.getParameter< edm::InputTag > ("RecHitsEE");
35 
36  EBOnly = iConfig.getParameter<bool> ("EBOnly") ;
37  sMin_min = iConfig.getParameter<double> ("sMin_min");
38  sMin_max = iConfig.getParameter<double> ("sMin_max");
39  sMaj_min = iConfig.getParameter<double> ("sMaj_min");
40  sMaj_max = iConfig.getParameter<double> ("sMaj_max");
41  seedTimeMin = iConfig.getParameter<double> ("seedTimeMin");
42  seedTimeMax = iConfig.getParameter<double> ("seedTimeMax");
43 
44  inputToken_ = consumes<trigger::TriggerFilterObjectWithRefs> (inputTag_);
45  rechitsEBToken_ = consumes<EcalRecHitCollection>(rechitsEB);
46  rechitsEEToken_ = consumes<EcalRecHitCollection>(rechitsEE);
47  inputTrkToken_ = consumes<reco::TrackCollection>(inputTrk);
48 
49 }
T getParameter(std::string const &) const
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > inputToken_
edm::EDGetTokenT< EcalRecHitCollection > rechitsEBToken_
edm::EDGetTokenT< EcalRecHitCollection > rechitsEEToken_
edm::EDGetTokenT< reco::TrackCollection > inputTrkToken_
HLTFilter(const edm::ParameterSet &config)
Definition: HLTFilter.cc:20
HLTDisplacedEgammaFilter::~HLTDisplacedEgammaFilter ( )
overridedefault

Member Function Documentation

void HLTDisplacedEgammaFilter::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 54 of file HLTDisplacedEgammaFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLTFilter::makeHLTFilterDescription().

54  {
57  desc.add<edm::InputTag>("inputTag", edm::InputTag("hltEGRegionalL1SingleEG22"));
58  desc.add<edm::InputTag>("l1EGCand", edm::InputTag("hltL1IsoRecoEcalCandidate"));
59  desc.add<edm::InputTag>("RecHitsEB", edm::InputTag("hltEcalRecHit", "EcalRecHitsEB"));
60  desc.add<edm::InputTag>("RecHitsEE", edm::InputTag("hltEcalRecHit", "EcalRecHitsEE"));
61  desc.add<edm::InputTag>("inputTrack", edm::InputTag("hltL1SeededEgammaRegionalCTFFinalFitWithMaterial"));
62  desc.add<int>("ncandcut", 1);
63  desc.add<bool>("EBOnly", false);
64  desc.add<double>("sMin_min", 0.1);
65  desc.add<double>("sMin_max", 0.4);
66  desc.add<double>("sMaj_min", 0.0);
67  desc.add<double>("sMaj_max", 999.0);
68  desc.add<double>("seedTimeMin", -25.0);
69  desc.add<double>("seedTimeMax", 25.0);
70  desc.add<int>("maxTrackCut", 0);
71  desc.add<double>("trackPtCut", 3.0);
72  desc.add<double>("trackdRCut", 0.5);
73  descriptions.add("hltDisplacedEgammaFilter", desc);
74 }
ParameterDescriptionBase * add(U const &iLabel, T const &value)
static void makeHLTFilterDescription(edm::ParameterSetDescription &desc)
Definition: HLTFilter.cc:29
void add(std::string const &label, ParameterSetDescription const &psetDescription)
bool HLTDisplacedEgammaFilter::hltFilter ( edm::Event iEvent,
const edm::EventSetup iSetup,
trigger::TriggerFilterObjectWithRefs filterproduct 
) const
overridevirtual

Implements HLTFilter.

Definition at line 79 of file HLTDisplacedEgammaFilter.cc.

References funct::abs(), accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), DEFINE_FWK_MODULE, HiRegitMuonDetachedTripletStep_cff::DeltaR, EBOnly, edm::SortedCollection< T, SORT >::find(), objects.autophobj::float, edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), inputToken_, inputTrkToken_, l1EGTag_, maxTrkCut, gen::n, ncandcut_, edm::Handle< T >::product(), TrackInfoProducer_cfi::rechits, rechitsEBToken_, rechitsEEToken_, HLTFilter::saveTags(), seedTimeMax, Cluster2ndMoments::sMaj, sMaj_max, Cluster2ndMoments::sMin, sMin_max, l1t::tracks, trigger::TriggerCluster, trigger::TriggerPhoton, trkdRCut, and trkPtCut.

80 {
81  using namespace trigger;
82 
83  // The filter object
84  if (saveTags()) {
85  filterproduct.addCollectionTag(l1EGTag_);
86  }
87 
88  // get hold of filtered candidates
89  //edm::Handle<reco::HLTFilterObjectWithRefs> recoecalcands;
91  iEvent.getByToken (inputToken_,PrevFilterOutput);
92 
93  // get hold of collection of objects
95  iEvent.getByToken( inputTrkToken_ , tracks);
96 
97  // get the EcalRecHit
100  iEvent.getByToken( rechitsEBToken_, rechitsEB_ );
101  iEvent.getByToken( rechitsEEToken_, rechitsEE_ );
102 
103  std::vector<edm::Ref<reco::RecoEcalCandidateCollection> > recoecalcands;
104  PrevFilterOutput->getObjects(TriggerCluster, recoecalcands);
105  if(recoecalcands.empty()) PrevFilterOutput->getObjects(TriggerPhoton, recoecalcands);
106 
107  // look at all candidates, check cuts and add to filter object
108  int n(0);
109 
110  for (auto const & ref : recoecalcands) {
111 
112  if ( EBOnly && std::abs( ref->eta() ) >= 1.479 ) continue ;
113 
114 
115  // S_Minor Cuts from the seed cluster
116  reco::CaloClusterPtr SCseed = ref->superCluster()->seed() ;
117  const EcalRecHitCollection* rechits = ( std::abs( ref->eta() ) < 1.479 ) ? rechitsEB_.product() : rechitsEE_.product() ;
118 
119  Cluster2ndMoments moments = EcalClusterTools::cluster2ndMoments(*SCseed, *rechits);
120  float sMin = moments.sMin ;
121  float sMaj = moments.sMaj ;
122  if ( sMin < sMin_min || sMin > sMin_max ) continue ;
123  if ( sMaj < sMaj_min || sMaj > sMaj_max ) continue ;
124 
125  // seed Time
126  std::pair<DetId, float> maxRH = EcalClusterTools::getMaximum( *SCseed, rechits );
127  DetId seedCrystalId = maxRH.first;
128  auto seedRH = rechits->find(seedCrystalId);
129  float seedTime = (float)seedRH->time();
130  if ( seedTime < seedTimeMin || seedTime > seedTimeMax ) continue ;
131 
132  //Track Veto
133 
134  int nTrk = 0 ;
135  for (auto const & it : *tracks) {
136  if ( it.pt() < trkPtCut ) continue ;
137  LorentzVector trkP4( it.px(), it.py(), it.pz(), it.p() ) ;
138  double dR = ROOT::Math::VectorUtil::DeltaR( trkP4 , ref->p4() ) ;
139  if ( dR < trkdRCut ) nTrk++ ;
140  if ( nTrk > maxTrkCut ) break ;
141  }
142  if ( nTrk > maxTrkCut ) continue ;
143 
144 
145  n++;
146  // std::cout << "Passed eta: " << ref->eta() << std::endl;
147  filterproduct.addObject(TriggerCluster, ref);
148  }
149 
150 
151  // filter decision
152  bool accept(n>=ncandcut_);
153 
154  return accept;
155 }
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > inputToken_
void getObjects(Vids &ids, VRphoton &refs) const
various physics-level getters:
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
edm::EDGetTokenT< EcalRecHitCollection > rechitsEBToken_
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:30
edm::EDGetTokenT< EcalRecHitCollection > rechitsEEToken_
void addObject(int id, const reco::RecoEcalCandidateRef &ref)
setters for L3 collections: (id=physics type, and Ref<C>)
edm::EDGetTokenT< reco::TrackCollection > inputTrkToken_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
math::XYZTLorentzVector LorentzVector
Definition: DetId.h:18
T const * product() const
Definition: Handle.h:74
void addCollectionTag(const edm::InputTag &collectionTag)
collectionTags
iterator find(key_type k)
bool saveTags() const
Definition: HLTFilter.h:45

Member Data Documentation

bool HLTDisplacedEgammaFilter::EBOnly
private

Definition at line 41 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

edm::InputTag HLTDisplacedEgammaFilter::inputTag_
private

Definition at line 32 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTDisplacedEgammaFilter::inputToken_
private

Definition at line 33 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

edm::InputTag HLTDisplacedEgammaFilter::inputTrk
private

Definition at line 49 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter().

edm::EDGetTokenT<reco::TrackCollection> HLTDisplacedEgammaFilter::inputTrkToken_
private

Definition at line 50 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

edm::InputTag HLTDisplacedEgammaFilter::l1EGTag_
private

Definition at line 35 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

int HLTDisplacedEgammaFilter::maxTrkCut
private

Definition at line 53 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

int HLTDisplacedEgammaFilter::ncandcut_
private

Definition at line 34 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

edm::InputTag HLTDisplacedEgammaFilter::rechitsEB
private

Definition at line 36 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter().

edm::EDGetTokenT<EcalRecHitCollection> HLTDisplacedEgammaFilter::rechitsEBToken_
private

Definition at line 38 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

edm::InputTag HLTDisplacedEgammaFilter::rechitsEE
private

Definition at line 37 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter().

edm::EDGetTokenT<EcalRecHitCollection> HLTDisplacedEgammaFilter::rechitsEEToken_
private

Definition at line 39 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

double HLTDisplacedEgammaFilter::seedTimeMax
private

Definition at line 47 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

double HLTDisplacedEgammaFilter::seedTimeMin
private

Definition at line 46 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter().

double HLTDisplacedEgammaFilter::sMaj_max
private

Definition at line 45 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

double HLTDisplacedEgammaFilter::sMaj_min
private

Definition at line 44 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter().

double HLTDisplacedEgammaFilter::sMin_max
private

Definition at line 43 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

double HLTDisplacedEgammaFilter::sMin_min
private

Definition at line 42 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter().

double HLTDisplacedEgammaFilter::trkdRCut
private

Definition at line 52 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().

double HLTDisplacedEgammaFilter::trkPtCut
private

Definition at line 51 of file HLTDisplacedEgammaFilter.h.

Referenced by HLTDisplacedEgammaFilter(), and hltFilter().