CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
ElectronSeedProducer.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: ElectronProducers
4 // Class: ElectronSeedProducer
5 //
13 //
14 // Original Author: Ursula Berthon, Claude Charlot
15 // Created: Mon Mar 27 13:22:06 CEST 2006
16 //
17 //
18 
19 #include <vector>
20 
21 #include "ElectronSeedProducer.h"
22 
23 #include "RecoEgamma/EgammaIsolationAlgos/interface/EgammaHcalIsolation.h"
24 
25 #include "RecoEgamma/EgammaElectronAlgos/interface/ElectronSeedGenerator.h"
26 #include "RecoEgamma/EgammaElectronAlgos/interface/ElectronHcalHelper.h"
27 #include "RecoEgamma/EgammaElectronAlgos/interface/SeedFilter.h"
28 #include "RecoEgamma/EgammaElectronAlgos/interface/ElectronUtilities.h"
29 #include "RecoTracker/MeasurementDet/interface/MeasurementTrackerEvent.h"
30 
31 #include "Geometry/Records/interface/CaloGeometryRecord.h"
32 #include "Geometry/Records/interface/CaloTopologyRecord.h"
33 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
34 #include "RecoCaloTools/Selectors/interface/CaloConeSelector.h"
35 
36 #include "DataFormats/EgammaReco/interface/ElectronSeed.h"
37 #include "DataFormats/EgammaReco/interface/ElectronSeedFwd.h"
38 #include "DataFormats/EcalDetId/interface/EBDetId.h"
39 #include "DataFormats/EcalDetId/interface/EEDetId.h"
40 #include "DataFormats/ForwardDetId/interface/ForwardSubdetector.h"
41 #include "DataFormats/ForwardDetId/interface/HGCalDetId.h"
42 
43 #include "FWCore/Framework/interface/Frameworkfwd.h"
44 #include "FWCore/Framework/interface/EDProducer.h"
45 #include "FWCore/Framework/interface/Event.h"
46 #include "FWCore/Framework/interface/MakerMacros.h"
47 #include "FWCore/Framework/interface/ConsumesCollector.h"
48 #include "FWCore/MessageLogger/interface/MessageLogger.h"
49 #include "FWCore/ParameterSet/interface/ParameterSet.h"
50 #include "FWCore/Utilities/interface/isFinite.h"
51 
52 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
53 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
54 
55 #include <string>
56 
57 using namespace reco ;
58 
59 ElectronSeedProducer::ElectronSeedProducer( const edm::ParameterSet& iConfig )
60  : //conf_(iConfig),
61  applyHOverECut_(true), hcalHelper_(0),
62  caloGeom_(0), caloGeomCacheId_(0), caloTopo_(0), caloTopoCacheId_(0)
63  {
64  conf_ = iConfig.getParameter<edm::ParameterSet>("SeedConfiguration") ;
65 
66  initialSeeds_ = consumes<TrajectorySeedCollection>(conf_.getParameter<edm::InputTag>("initialSeeds")) ;
67  SCEtCut_ = conf_.getParameter<double>("SCEtCut");
68  fromTrackerSeeds_ = conf_.getParameter<bool>("fromTrackerSeeds") ;
69  prefilteredSeeds_ = conf_.getParameter<bool>("preFilteredSeeds") ;
70 
71  auto theconsumes = consumesCollector();
72 
73  // new beamSpot tag
74  beamSpotTag_ = consumes<reco::BeamSpot>(conf_.getParameter<edm::InputTag>("beamSpot"));
75 
76  // for H/E
77  applyHOverECut_ = conf_.getParameter<bool>("applyHOverECut") ;
78  if (applyHOverECut_)
79  {
80  ElectronHcalHelper::Configuration hcalCfg ;
81  hcalCfg.hOverEConeSize = conf_.getParameter<double>("hOverEConeSize") ;
82  if (hcalCfg.hOverEConeSize>0)
83  {
84  hcalCfg.useTowers = true ;
85  hcalCfg.hcalTowers =
86  consumes<CaloTowerCollection>(conf_.getParameter<edm::InputTag>("hcalTowers")) ;
87  hcalCfg.hOverEPtMin = conf_.getParameter<double>("hOverEPtMin") ;
88  }
89  hcalHelper_ = new ElectronHcalHelper(hcalCfg) ;
90 
91  allowHGCal_ = conf_.getParameter<bool>("allowHGCal");
92  if( allowHGCal_ ) {
93  const edm::ParameterSet& hgcCfg = conf_.getParameterSet("HGCalConfig");
94  hgcClusterTools_.reset( new hgcal::ClusterTools(hgcCfg, theconsumes) );
95  }
96 
97  maxHOverEBarrel_=conf_.getParameter<double>("maxHOverEBarrel") ;
98  maxHOverEEndcaps_=conf_.getParameter<double>("maxHOverEEndcaps") ;
99  maxHBarrel_=conf_.getParameter<double>("maxHBarrel") ;
100  maxHEndcaps_=conf_.getParameter<double>("maxHEndcaps") ;
101  }
102 
103  applySigmaIEtaIEtaCut_ = conf_.getParameter<bool>("applySigmaIEtaIEtaCut");
104 
105  // apply sigma_ieta_ieta cut
106  if (applySigmaIEtaIEtaCut_ == true)
107  {
108  maxSigmaIEtaIEtaBarrel_ = conf_.getParameter<double>("maxSigmaIEtaIEtaBarrel");
109  maxSigmaIEtaIEtaEndcaps_ = conf_.getParameter<double>("maxSigmaIEtaIEtaEndcaps");
110  }
111 
112  edm::ParameterSet rpset = conf_.getParameter<edm::ParameterSet>("RegionPSet");
113  filterVtxTag_ = consumes<std::vector<reco::Vertex> >(rpset.getParameter<edm::InputTag> ("VertexProducer"));
114 
115  ElectronSeedGenerator::Tokens esg_tokens;
116  esg_tokens.token_bs = beamSpotTag_;
117  esg_tokens.token_vtx = mayConsume<reco::VertexCollection>(conf_.getParameter<edm::InputTag>("vertices"));
118  esg_tokens.token_measTrkEvt= consumes<MeasurementTrackerEvent>(conf_.getParameter<edm::InputTag>("measurementTrackerEvent"));
119 
120  matcher_ = new ElectronSeedGenerator(conf_,esg_tokens) ;
121 
122  // get collections from config
123  if (applySigmaIEtaIEtaCut_ == true) {
124  ebRecHitCollection_ = consumes<EcalRecHitCollection> (iConfig.getParameter<edm::InputTag>("ebRecHitCollection"));
125  eeRecHitCollection_ = consumes<EcalRecHitCollection> (iConfig.getParameter<edm::InputTag>("eeRecHitCollection"));
126  }
127 
128  superClusters_[0]=
129  consumes<reco::SuperClusterCollection>(iConfig.getParameter<edm::InputTag>("barrelSuperClusters")) ;
130  superClusters_[1]=
131  consumes<reco::SuperClusterCollection>(iConfig.getParameter<edm::InputTag>("endcapSuperClusters")) ;
132 
133  // Construction of SeedFilter was in beginRun() with the comment
134  // below, but it has to be done here because of ConsumesCollector
135  //
136  // FIXME: because of a bug presumably in tracker seeding,
137  // perhaps in CombinedHitPairGenerator, badly caching some EventSetup product,
138  // we must redo the SeedFilter for each run.
139  if (prefilteredSeeds_) {
140  SeedFilter::Tokens sf_tokens;
141  sf_tokens.token_bs = beamSpotTag_;
142  sf_tokens.token_vtx = filterVtxTag_;
143  edm::ConsumesCollector iC = consumesCollector();
144  seedFilter_.reset(new SeedFilter(conf_, sf_tokens, iC));
145  }
146 
147  //register your products
148  produces<ElectronSeedCollection>() ;
149 }
150 
151 void ElectronSeedProducer::beginRun(edm::Run const&, edm::EventSetup const&)
152 {}
153 
154 void ElectronSeedProducer::endRun(edm::Run const&, edm::EventSetup const&)
155 {}
156 
157 ElectronSeedProducer::~ElectronSeedProducer()
158  {
159  delete hcalHelper_ ;
160  delete matcher_ ;
161  }
162 
163 void ElectronSeedProducer::produce(edm::Event& e, const edm::EventSetup& iSetup)
164  {
165  LogDebug("ElectronSeedProducer") <<"[ElectronSeedProducer::produce] entering " ;
166 
167  edm::Handle<reco::BeamSpot> theBeamSpot ;
168  e.getByToken(beamSpotTag_,theBeamSpot) ;
169 
170  if (hcalHelper_)
171  {
172  hcalHelper_->checkSetup(iSetup) ;
173  hcalHelper_->readEvent(e) ;
174  if( allowHGCal_ ) {
175  hgcClusterTools_->getEventSetup(iSetup);
176  hgcClusterTools_->getEvent(e);
177  }
178  }
179 
180  // get calo geometry
181  if (caloGeomCacheId_!=iSetup.get<CaloGeometryRecord>().cacheIdentifier()) {
182  iSetup.get<CaloGeometryRecord>().get(caloGeom_);
183  caloGeomCacheId_=iSetup.get<CaloGeometryRecord>().cacheIdentifier();
184  }
185  if (caloTopoCacheId_!=iSetup.get<CaloTopologyRecord>().cacheIdentifier()){
186  caloTopoCacheId_=iSetup.get<CaloTopologyRecord>().cacheIdentifier();
187  iSetup.get<CaloTopologyRecord>().get(caloTopo_);
188  }
189 
190  matcher_->setupES(iSetup);
191 
192  // get initial TrajectorySeeds if necessary
193  if (fromTrackerSeeds_)
194  {
195  if (!prefilteredSeeds_)
196  {
197  edm::Handle<TrajectorySeedCollection> hSeeds;
198  e.getByToken(initialSeeds_, hSeeds);
199  theInitialSeedColl = const_cast<TrajectorySeedCollection *> (hSeeds.product());
200  }
201  else
202  { theInitialSeedColl = new TrajectorySeedCollection ; }
203  }
204  else
205  { theInitialSeedColl = 0 ; } // not needed in this case
206 
207  ElectronSeedCollection * seeds = new ElectronSeedCollection ;
208 
209  // loop over barrel + endcap
210  for (unsigned int i=0; i<2; i++) {
211  edm::Handle<SuperClusterCollection> clusters ;
212  e.getByToken(superClusters_[i],clusters);
213  SuperClusterRefVector clusterRefs ;
214  std::vector<float> hoe1s, hoe2s ;
215  filterClusters(*theBeamSpot,clusters,/*mhbhe_,*/clusterRefs,hoe1s,hoe2s,e, iSetup);
216  if ((fromTrackerSeeds_) && (prefilteredSeeds_))
217  { filterSeeds(e,iSetup,clusterRefs) ; }
218  matcher_->run(e,iSetup,clusterRefs,hoe1s,hoe2s,theInitialSeedColl,*seeds);
219  }
220 
221  // store the accumulated result
222  std::unique_ptr<ElectronSeedCollection> pSeeds(seeds);
223  ElectronSeedCollection::iterator is ;
224  for ( is=pSeeds->begin() ; is!=pSeeds->end() ; is++ ) {
225  edm::RefToBase<CaloCluster> caloCluster = is->caloCluster() ;
226  SuperClusterRef superCluster = caloCluster.castTo<SuperClusterRef>() ;
227  LogDebug("ElectronSeedProducer")
228  << "new seed with "
229  << (*is).nHits() << " hits"
230  << ", charge " << (*is).getCharge()
231  << " and cluster energy " << superCluster->energy()
232  << " PID "<<superCluster.id() ;
233  }
234  e.put(std::move(pSeeds));
235  if (fromTrackerSeeds_ && prefilteredSeeds_) delete theInitialSeedColl;
236  }
237 
238 
239 //===============================
240 // Filter the superclusters
241 // - with EtCut
242 // - with HoE using calo cone
243 //===============================
244 
245 void ElectronSeedProducer::filterClusters
246  ( const reco::BeamSpot & bs,
247  const edm::Handle<reco::SuperClusterCollection> & superClusters,
248  SuperClusterRefVector & sclRefs,
249  std::vector<float> & hoe1s, std::vector<float> & hoe2s,
250  edm::Event & event, const edm::EventSetup & setup)
251  {
252 
253  std::vector<float> sigmaIEtaIEtaEB_;
254  std::vector<float> sigmaIEtaIEtaEE_;
255 
256  for (unsigned int i=0;i<superClusters->size();++i)
257  {
258  const SuperCluster & scl = (*superClusters)[i] ;
259  double sclEta = EleRelPoint(scl.position(),bs.position()).eta() ;
260  if (scl.energy()/cosh(sclEta)>SCEtCut_)
261  {
262 // if ((applyHOverECut_==true)&&((hcalHelper_->hcalESum(scl)/scl.energy()) > maxHOverE_))
263 // { continue ; }
264 // sclRefs.push_back(edm::Ref<reco::SuperClusterCollection>(superClusters,i)) ;
265  double had1, had2, had, scle ;
266 
267  bool HoeVeto = false ;
268  if (applyHOverECut_==true)
269  {
270  had1 = hcalHelper_->hcalESumDepth1(scl);
271  had2 = hcalHelper_->hcalESumDepth2(scl);
272  had = had1+had2 ;
273  scle = scl.energy() ;
274  int det_group = scl.seed()->hitsAndFractions()[0].first.det() ;
275  int detector = scl.seed()->hitsAndFractions()[0].first.subdetId() ;
276  if (detector==EcalBarrel && (had<maxHBarrel_ || had/scle<maxHOverEBarrel_)) HoeVeto=true;
277  else if( detector==EcalEndcap && (had<maxHEndcaps_ || had/scle<maxHOverEEndcaps_) ) HoeVeto=true;
278  else if( allowHGCal_ && (detector==HcalEndcap || det_group == DetId::Forward) ) {
279  float had_fraction = hgcClusterTools_->getClusterHadronFraction(*(scl.seed()));
280  had1 = had_fraction*scl.seed()->energy();
281  had2 = 0.;
282  HoeVeto= ( had_fraction >= 0.f && had_fraction < maxHOverEEndcaps_ );
283  }
284  if (HoeVeto)
285  {
286  sclRefs.push_back(edm::Ref<reco::SuperClusterCollection>(superClusters,i)) ;
287  hoe1s.push_back(had1/scle) ;
288  hoe2s.push_back(had2/scle) ;
289  }
290  }
291  else
292  {
293  sclRefs.push_back(edm::Ref<reco::SuperClusterCollection>(superClusters,i)) ;
294  hoe1s.push_back(std::numeric_limits<float>::infinity()) ;
295  hoe2s.push_back(std::numeric_limits<float>::infinity()) ;
296  }
297  }
298 
299  if (applySigmaIEtaIEtaCut_ == true)
300  {
301  noZS::EcalClusterLazyTools lazyTool_noZS(event, setup, ebRecHitCollection_, eeRecHitCollection_);
302  std::vector<float> vCov = lazyTool_noZS.localCovariances(*(scl.seed()));
303  int detector = scl.seed()->hitsAndFractions()[0].first.subdetId() ;
304  if (detector==EcalBarrel) sigmaIEtaIEtaEB_ .push_back(edm::isNotFinite(vCov[0]) ? 0. : sqrt(vCov[0]));
305  if (detector==EcalEndcap) sigmaIEtaIEtaEE_ .push_back(edm::isNotFinite(vCov[0]) ? 0. : sqrt(vCov[0]));
306  }
307  }
308  LogDebug("ElectronSeedProducer")<<"Filtered out "<<sclRefs.size()<<" superclusters from "<<superClusters->size() ;
309  }
310 
311 void ElectronSeedProducer::filterSeeds
312  ( edm::Event & event, const edm::EventSetup & setup,
313  reco::SuperClusterRefVector & sclRefs )
314  {
315  for ( unsigned int i=0 ; i<sclRefs.size() ; ++i )
316  {
317  seedFilter_->seeds(event,setup,sclRefs[i],theInitialSeedColl) ;
318  LogDebug("ElectronSeedProducer")<<"Number of Seeds: "<<theInitialSeedColl->size() ;
319  }
320  }
321 
322 void
323 ElectronSeedProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
324  edm::ParameterSetDescription desc;
325  desc.add<edm::InputTag>("endcapSuperClusters",edm::InputTag("particleFlowSuperClusterECAL","particleFlowSuperClusterECALEndcapWithPreshower"));
326  {
327  edm::ParameterSetDescription psd0, psd1, psd2, psd3, psd4;
328  psd1.add<unsigned int>("maxElement", 0);
329  psd1.add<std::string>("ComponentName", std::string("StandardHitPairGenerator"));
330  psd1.addUntracked<int>("useOnDemandTracker", 0);
331  psd1.add<edm::InputTag>("SeedingLayers", edm::InputTag("hltMixedLayerPairs"));
332  psd0.add<edm::ParameterSetDescription>("OrderedHitsFactoryPSet", psd1);
333 
334  psd2.add<double>("deltaPhiRegion", 0.4);
335  psd2.add<double>("originHalfLength", 15.0);
336  psd2.add<bool>("useZInVertex", true);
337  psd2.add<double>("deltaEtaRegion", 0.1);
338  psd2.add<double>("ptMin", 1.5 );
339  psd2.add<double>("originRadius", 0.2);
340  psd2.add<edm::InputTag>("VertexProducer", edm::InputTag("dummyVertices"));
341  psd0.add<edm::ParameterSetDescription>("RegionPSet", psd2);
342 
343  psd0.add<double>("PhiMax2B",0.002);
344  psd0.add<double>("hOverEPtMin",0.0);
345  psd0.add<double>("PhiMax2F",0.003);
346  psd0.add<bool>("searchInTIDTEC",true);
347  psd0.add<double>("pPhiMax1",0.125);
348  psd0.add<double>("HighPtThreshold",35.0);
349  psd0.add<double>("r2MinF",-0.15);
350  psd0.add<double>("maxHBarrel",0.0);
351  psd0.add<double>("DeltaPhi1Low",0.23);
352  psd0.add<double>("DeltaPhi1High",0.08);
353  psd0.add<double>("ePhiMin1",-0.125);
354  psd0.add<edm::InputTag>("hcalTowers",edm::InputTag("towerMaker"));
355  psd0.add<double>("LowPtThreshold",5.0);
356  psd0.add<double>("maxHOverEBarrel",0.15);
357  psd0.add<double>("maxSigmaIEtaIEtaBarrel", 0.5);
358  psd0.add<double>("maxSigmaIEtaIEtaEndcaps", 0.5);
359  psd0.add<bool>("dynamicPhiRoad",true);
360  psd0.add<double>("ePhiMax1",0.075);
361  psd0.add<std::string>("measurementTrackerName","");
362  psd0.add<double>("SizeWindowENeg",0.675);
363  psd0.add<double>("nSigmasDeltaZ1",5.0);
364  psd0.add<double>("rMaxI",0.2);
365  psd0.add<double>("maxHEndcaps",0.0);
366  psd0.add<bool>("preFilteredSeeds",false);
367  psd0.add<double>("r2MaxF",0.15);
368  psd0.add<double>("hOverEConeSize",0.15);
369  psd0.add<double>("pPhiMin1",-0.075);
370  psd0.add<edm::InputTag>("initialSeeds",edm::InputTag("newCombinedSeeds"));
371  psd0.add<double>("deltaZ1WithVertex",25.0);
372  psd0.add<double>("SCEtCut",0.0);
373  psd0.add<double>("z2MaxB",0.09);
374  psd0.add<bool>("fromTrackerSeeds",true);
375  psd0.add<edm::InputTag>("hcalRecHits",edm::InputTag("hbhereco"));
376  psd0.add<double>("z2MinB",-0.09);
377  psd0.add<double>("rMinI",-0.2);
378  psd0.add<double>("maxHOverEEndcaps",0.15);
379  psd0.add<double>("hOverEHBMinE",0.7);
380  psd0.add<bool>("useRecoVertex",false);
381  psd0.add<edm::InputTag>("beamSpot",edm::InputTag("offlineBeamSpot"));
382  psd0.add<edm::InputTag>("measurementTrackerEvent",edm::InputTag("MeasurementTrackerEvent"));
383  psd0.add<edm::InputTag>("vertices",edm::InputTag("offlinePrimaryVerticesWithBS"));
384  psd0.add<bool>("applyHOverECut",true);
385  psd0.add<edm::InputTag>("ebRecHitCollection", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
386  psd0.add<edm::InputTag>("eeRecHitCollection", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
387  psd0.add<bool>("applySigmaIEtaIEtaCut", false);
388  psd0.add<double>("DeltaPhi2F",0.012);
389  psd0.add<double>("PhiMin2F",-0.003);
390  psd0.add<double>("hOverEHFMinE",0.8);
391  psd0.add<double>("DeltaPhi2B",0.008);
392  psd0.add<double>("PhiMin2B",-0.002);
393  psd0.add<bool>("allowHGCal",false);
394 
395  psd3.add<std::string>("ComponentName",std::string("SeedFromConsecutiveHitsCreator"));
396  psd3.add<std::string>("propagator",std::string("PropagatorWithMaterial"));
397  psd3.add<double>("SeedMomentumForBOFF",5.0);
398  psd3.add<double>("OriginTransverseErrorMultiplier",1.0);
399  psd3.add<double>("MinOneOverPtError",1.0);
400  psd3.add<std::string>("magneticField",std::string(""));
401  psd3.add<std::string>("TTRHBuilder",std::string("WithTrackAngle"));
402  psd3.add<bool>("forceKinematicWithRegionDirection",false);
403 
404  psd4.add<edm::InputTag>("HGCEEInput",edm::InputTag("HGCalRecHit","HGCEERecHits"));
405  psd4.add<edm::InputTag>("HGCFHInput",edm::InputTag("HGCalRecHit","HGCHEFRecHits"));
406  psd4.add<edm::InputTag>("HGCBHInput",edm::InputTag("HGCalRecHit","HGCHEBRecHits"));
407 
408  psd0.add<edm::ParameterSetDescription>("HGCalConfig",psd4);
409 
410  psd0.add<edm::ParameterSetDescription>("SeedCreatorPSet",psd3);
411 
412  desc.add<edm::ParameterSetDescription>("SeedConfiguration",psd0);
413  }
414  desc.add<edm::InputTag>("barrelSuperClusters",edm::InputTag("particleFlowSuperClusterECAL","particleFlowSuperClusterECALBarrel"));
415  descriptions.add("ecalDrivenElectronSeeds",desc);
416 }
LogDebug
#define LogDebug(id)
Definition: PFTrackAlgoTools.cc:501
ElectronSeedProducer::filterClusters
void filterClusters(const reco::BeamSpot &bs, const edm::Handle< reco::SuperClusterCollection > &superClusters, reco::SuperClusterRefVector &sclRefs, std::vector< float > &hoe1s, std::vector< float > &hoe2s, edm::Event &e, const edm::EventSetup &setup)
Definition: ElectronSeedProducer.cc:246
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::eventsetup::EventSetupRecord::cacheIdentifier
unsigned long long cacheIdentifier() const
Definition: EventSetupRecord.h:174
reco::CaloCluster::position
const math::XYZPoint & position() const
cluster centroid position
Definition: CaloCluster.h:126
i
int i
Definition: DBlmapReader.cc:9
ElectronSeedProducer::endRun
virtual void endRun(edm::Run const &, edm::EventSetup const &) overridefinal
Definition: ElectronSeedProducer.cc:154
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
ElectronSeedGenerator::Tokens::token_bs
edm::EDGetTokenT< reco::BeamSpot > token_bs
Definition: ElectronSeedGenerator.h:55
ElectronHcalHelper::readEvent
void readEvent(const edm::Event &)
Definition: ElectronHcalHelper.cc:38
edm::ParameterSetDescription::addUntracked
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:94
ElectronSeedGenerator::setupES
void setupES(const edm::EventSetup &setup)
Definition: ElectronSeedGenerator.cc:146
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
ElectronSeedProducer::beginRun
virtual void beginRun(edm::Run const &, edm::EventSetup const &) overridefinal
Definition: ElectronSeedProducer.cc:151
ElectronSeedProducer::caloGeomCacheId_
unsigned long long caloGeomCacheId_
Definition: ElectronSeedProducer.h:86
ElectronSeedProducer::filterVtxTag_
edm::EDGetTokenT< std::vector< reco::Vertex > > filterVtxTag_
Definition: ElectronSeedProducer.h:67
ElectronHcalHelper::checkSetup
void checkSetup(const edm::EventSetup &)
Definition: ElectronHcalHelper.cc:15
ElectronSeedProducer::initialSeeds_
edm::EDGetTokenT< TrajectorySeedCollection > initialSeeds_
Definition: ElectronSeedProducer.h:66
ElectronSeedGenerator::run
void run(edm::Event &, const edm::EventSetup &setup, const reco::SuperClusterRefVector &, const std::vector< float > &hoe1s, const std::vector< float > &hoe2s, TrajectorySeedCollection *seeds, reco::ElectronSeedCollection &)
Definition: ElectronSeedGenerator.cc:221
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::Ref::id
ProductID id() const
Accessor for product ID.
Definition: Ref.h:258
ElectronSeedGenerator::Tokens::token_measTrkEvt
edm::EDGetTokenT< MeasurementTrackerEvent > token_measTrkEvt
Definition: ElectronSeedGenerator.h:56
ElectronSeedProducer::eeRecHitCollection_
edm::EDGetTokenT< EcalRecHitCollection > eeRecHitCollection_
Definition: ElectronSeedProducer.h:70
HLT_FULL_cff.clusters
tuple clusters
Definition: HLT_FULL_cff.py:10237
edm::isNotFinite
bool isNotFinite(T x)
Definition: isFinite.h:10
ElectronSeedProducer::caloTopo_
edm::ESHandle< CaloTopology > caloTopo_
Definition: ElectronSeedProducer.h:87
TrajectorySeedCollection
std::vector< TrajectorySeed > TrajectorySeedCollection
Definition: TrajectorySeedCollection.h:6
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:72
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
ElectronSeedProducer::ebRecHitCollection_
edm::EDGetTokenT< EcalRecHitCollection > ebRecHitCollection_
Definition: ElectronSeedProducer.h:69
eostools.move
def move
Definition: eostools.py:510
infinity
const double infinity
Definition: CSCChamberFitter.cc:10
ElectronSeedProducer::ElectronSeedProducer
ElectronSeedProducer(const edm::ParameterSet &)
Definition: ElectronSeedProducer.cc:59
reco::CaloCluster::energy
double energy() const
cluster energy
Definition: CaloCluster.h:121
reco::ElectronSeedCollection
std::vector< ElectronSeed > ElectronSeedCollection
collection of ElectronSeed objects
Definition: ElectronSeedFwd.h:10
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:89
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
ElectronSeedProducer::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: ElectronSeedProducer.cc:323
ElectronSeedProducer::seedFilter_
std::unique_ptr< SeedFilter > seedFilter_
Definition: ElectronSeedProducer.h:74
ElectronSeedProducer::conf_
edm::ParameterSet conf_
Definition: ElectronSeedProducer.h:72
ElectronSeedProducer::caloGeom_
edm::ESHandle< CaloGeometry > caloGeom_
Definition: ElectronSeedProducer.h:85
edm::Handle::product
T const * product() const
Definition: Handle.h:81
ElectronSeedProducer::matcher_
ElectronSeedGenerator * matcher_
Definition: ElectronSeedProducer.h:73
edm::ParameterSet::getParameterSet
ParameterSet const & getParameterSet(std::string const &) const
Definition: ParameterSet.cc:2364
ElectronSeedProducer::theInitialSeedColl
TrajectorySeedCollection * theInitialSeedColl
Definition: ElectronSeedProducer.h:76
EcalClusterLazyToolsT::localCovariances
std::vector< float > localCovariances(const reco::BasicCluster &cluster, float w0=4.7)
Definition: EcalClusterLazyTools.h:363
edm::RefToBase::castTo
REF castTo() const
Definition: RefToBase.h:286
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:65
ElectronSeedProducer::superClusters_
edm::EDGetTokenT< reco::SuperClusterCollection > superClusters_[2]
Definition: ElectronSeedProducer.h:65
ElectronSeedProducer::caloTopoCacheId_
unsigned long long caloTopoCacheId_
Definition: ElectronSeedProducer.h:88
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:83208
SeedFilter::Tokens::token_bs
edm::EDGetTokenT< reco::BeamSpot > token_bs
Definition: SeedFilter.h:37
ElectronSeedProducer::filterSeeds
void filterSeeds(edm::Event &e, const edm::EventSetup &setup, reco::SuperClusterRefVector &sclRefs)
Definition: ElectronSeedProducer.cc:312
edm::RefVector::push_back
void push_back(value_type const &ref)
Add a Ref&lt;C, T&gt; to the RefVector.
Definition: RefVector.h:69
edm::RefVector::size
size_type size() const
Size of the RefVector.
Definition: RefVector.h:107
ElectronHcalHelper::Configuration::hcalTowers
edm::EDGetTokenT< CaloTowerCollection > hcalTowers
Definition: ElectronHcalHelper.h:32
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
reco::SuperCluster::seed
const CaloClusterPtr & seed() const
seed BasicCluster
Definition: SuperCluster.h:66
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
ElectronSeedProducer::beamSpotTag_
edm::EDGetTokenT< reco::BeamSpot > beamSpotTag_
Definition: ElectronSeedProducer.h:68
ElectronSeedGenerator::Tokens::token_vtx
edm::EDGetTokenT< std::vector< reco::Vertex > > token_vtx
Definition: ElectronSeedGenerator.h:54
ElectronSeedProducer::hcalHelper_
ElectronHcalHelper * hcalHelper_
Definition: ElectronSeedProducer.h:82
SeedFilter::Tokens::token_vtx
edm::EDGetTokenT< std::vector< reco::Vertex > > token_vtx
Definition: SeedFilter.h:36
ElectronSeedProducer::hgcClusterTools_
std::unique_ptr< hgcal::ClusterTools > hgcClusterTools_
Definition: ElectronSeedProducer.h:110
edm::Run
Definition: Run.h:42
ElectronSeedProducer::produce
virtual void produce(edm::Event &, const edm::EventSetup &) overridefinal
Definition: ElectronSeedProducer.cc:163
Generated for CMSSW Reference Manual by doxygen 1.8.5