CMS 3D CMS Logo

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

#include <EgammaHLTIslandClusterProducer.h>

Inheritance diagram for EgammaHLTIslandClusterProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 EgammaHLTIslandClusterProducer (const edm::ParameterSet &ps)
 
void produce (edm::Event &, const edm::EventSetup &) override
 
 ~EgammaHLTIslandClusterProducer () override
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducer () 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
 
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)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 

Private Member Functions

void clusterizeECALPart (edm::Event &evt, const edm::EventSetup &es, edm::EDGetTokenT< EcalRecHitCollection > &hitToken, const std::string &clusterCollection, const std::vector< RectangularEtaPhiRegion > &regions, const IslandClusterAlgo::EcalPart &ecalPart)
 
bool counterExceeded () const
 
const EcalRecHitCollectiongetCollection (edm::Event &evt, edm::EDGetTokenT< EcalRecHitCollection > &hitToken)
 

Private Attributes

std::string barrelClusterCollection_
 
edm::InputTag barrelHitCollection_
 
edm::EDGetTokenT< EcalRecHitCollectionbarrelHitToken_
 
bool doBarrel_
 
bool doEndcaps_
 
bool doIsolated_
 
std::string endcapClusterCollection_
 
edm::InputTag endcapHitCollection_
 
edm::EDGetTokenT< EcalRecHitCollectionendcapHitToken_
 
IslandClusterAlgoisland_p
 
double l1LowerThr_
 
double l1LowerThrIgnoreIsolation_
 
edm::EDGetTokenT< l1extra::L1EmParticleCollectionl1TagIsolated_
 
edm::EDGetTokenT< l1extra::L1EmParticleCollectionl1TagNonIsolated_
 
double l1UpperThr_
 
int nEvt_
 
int nMaxPrintout_
 
PositionCalc posCalculator_
 
double regionEtaMargin_
 
double regionPhiMargin_
 
IslandClusterAlgo::VerbosityLevel verbosity
 

Additional Inherited Members

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

Definition at line 25 of file EgammaHLTIslandClusterProducer.h.

Constructor & Destructor Documentation

EgammaHLTIslandClusterProducer::EgammaHLTIslandClusterProducer ( const edm::ParameterSet ps)

Definition at line 36 of file EgammaHLTIslandClusterProducer.cc.

References barrelClusterCollection_, barrelHitCollection_, barrelHitToken_, doBarrel_, doEndcaps_, doIsolated_, endcapClusterCollection_, endcapHitCollection_, endcapHitToken_, edm::ParameterSet::getParameter(), island_p, l1LowerThr_, l1LowerThrIgnoreIsolation_, l1TagIsolated_, l1TagNonIsolated_, l1UpperThr_, nEvt_, IslandClusterAlgo::pDEBUG, IslandClusterAlgo::pERROR, IslandClusterAlgo::pINFO, posCalculator_, IslandClusterAlgo::pWARNING, regionEtaMargin_, regionPhiMargin_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

37 {
38  // The verbosity level
39  std::string verbosityString = ps.getParameter<std::string>("VerbosityLevel");
40  if (verbosityString == "DEBUG") verbosity = IslandClusterAlgo::pDEBUG;
41  else if (verbosityString == "WARNING") verbosity = IslandClusterAlgo::pWARNING;
42  else if (verbosityString == "INFO") verbosity = IslandClusterAlgo::pINFO;
44 
45  doBarrel_ = ps.getParameter<bool>("doBarrel");
46  doEndcaps_ = ps.getParameter<bool>("doEndcaps");
47  doIsolated_ = ps.getParameter<bool>("doIsolated");
48 
49  // Parameters to identify the hit collections
50  barrelHitCollection_ = ps.getParameter<edm::InputTag>("barrelHitProducer");
51  endcapHitCollection_ = ps.getParameter<edm::InputTag>("endcapHitProducer");
52  barrelHitToken_ = consumes<EcalRecHitCollection>(barrelHitCollection_);
53  endcapHitToken_ = consumes<EcalRecHitCollection>(endcapHitCollection_);
54 
55  // The names of the produced cluster collections
56  barrelClusterCollection_ = ps.getParameter<std::string>("barrelClusterCollection");
57  endcapClusterCollection_ = ps.getParameter<std::string>("endcapClusterCollection");
58 
59  // Island algorithm parameters
60  double barrelSeedThreshold = ps.getParameter<double>("IslandBarrelSeedThr");
61  double endcapSeedThreshold = ps.getParameter<double>("IslandEndcapSeedThr");
62 
63  // L1 matching parameters
64  l1TagIsolated_ = consumes<l1extra::L1EmParticleCollection>(ps.getParameter< edm::InputTag > ("l1TagIsolated"));
65  l1TagNonIsolated_ = consumes<l1extra::L1EmParticleCollection>(ps.getParameter< edm::InputTag > ("l1TagNonIsolated"));
66  l1LowerThr_ = ps.getParameter<double> ("l1LowerThr");
67  l1UpperThr_ = ps.getParameter<double> ("l1UpperThr");
68  l1LowerThrIgnoreIsolation_ = ps.getParameter<double> ("l1LowerThrIgnoreIsolation");
69 
70  regionEtaMargin_ = ps.getParameter<double>("regionEtaMargin");
71  regionPhiMargin_ = ps.getParameter<double>("regionPhiMargin");
72 
73  // Parameters for the position calculation:
74  posCalculator_ = PositionCalc( ps.getParameter<edm::ParameterSet>("posCalcParameters") );
75 
76  // Produces a collection of barrel and a collection of endcap clusters
77 
78  produces< reco::BasicClusterCollection >(endcapClusterCollection_);
79  produces< reco::BasicClusterCollection >(barrelClusterCollection_);
80 
81  island_p = new IslandClusterAlgo(barrelSeedThreshold, endcapSeedThreshold, posCalculator_,verbosity);
82 
83  nEvt_ = 0;
84 }
T getParameter(std::string const &) const
edm::EDGetTokenT< EcalRecHitCollection > endcapHitToken_
edm::EDGetTokenT< l1extra::L1EmParticleCollection > l1TagNonIsolated_
IslandClusterAlgo::VerbosityLevel verbosity
edm::EDGetTokenT< l1extra::L1EmParticleCollection > l1TagIsolated_
edm::EDGetTokenT< EcalRecHitCollection > barrelHitToken_
EgammaHLTIslandClusterProducer::~EgammaHLTIslandClusterProducer ( )
override

Definition at line 86 of file EgammaHLTIslandClusterProducer.cc.

References island_p.

86  {
87  delete island_p;
88 }

Member Function Documentation

void EgammaHLTIslandClusterProducer::clusterizeECALPart ( edm::Event evt,
const edm::EventSetup es,
edm::EDGetTokenT< EcalRecHitCollection > &  hitToken,
const std::string &  clusterCollection,
const std::vector< RectangularEtaPhiRegion > &  regions,
const IslandClusterAlgo::EcalPart ecalPart 
)
private

Definition at line 261 of file EgammaHLTIslandClusterProducer.cc.

References IslandClusterAlgo::barrel, barrelClusterCollection_, fastPrimaryVertexProducer_cfi::clusters, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, endcapClusterCollection_, edm::EventSetup::get(), getCollection(), CaloGeometry::getSubdetectorGeometry(), island_p, IslandClusterAlgo::makeClusters(), eostools::move(), and edm::Event::put().

Referenced by produce().

266 {
267  // get the hit collection from the event:
268  const EcalRecHitCollection *hitCollection_p = getCollection(evt, hitToken);
269 
270  // get the geometry and topology from the event setup:
271  edm::ESHandle<CaloGeometry> geoHandle;
272  es.get<CaloGeometryRecord>().get(geoHandle);
273 
274  const CaloSubdetectorGeometry *geometry_p;
275  CaloSubdetectorTopology *topology_p;
276 
277  if (ecalPart == IslandClusterAlgo::barrel)
278  {
279  geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
280  topology_p = new EcalBarrelTopology(geoHandle);
281  }
282  else
283  {
284  geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
285  topology_p = new EcalEndcapTopology(geoHandle);
286  }
287 
288  const CaloSubdetectorGeometry *geometryES_p;
289  geometryES_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
290 
291  // Run the clusterization algorithm:
293  clusters = island_p->makeClusters(hitCollection_p, geometry_p, topology_p, geometryES_p, ecalPart, true, regions);
294 
295  // create an unique_ptr to a BasicClusterCollection, copy the barrel clusters into it and put in the Event:
296  auto clusters_p = std::make_unique<reco::BasicClusterCollection>();
297  clusters_p->assign(clusters.begin(), clusters.end());
299  if (ecalPart == IslandClusterAlgo::barrel)
300  bccHandle = evt.put(std::move(clusters_p), barrelClusterCollection_);
301  else
302  bccHandle = evt.put(std::move(clusters_p), endcapClusterCollection_);
303 
304  delete topology_p;
305 }
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:49
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:137
std::vector< reco::BasicCluster > makeClusters(const EcalRecHitCollection *hits, const CaloSubdetectorGeometry *geometry, const CaloSubdetectorTopology *topology_p, const CaloSubdetectorGeometry *geometryES_p, EcalPart ecalPart, bool regional=false, const std::vector< RectangularEtaPhiRegion > &regions=std::vector< RectangularEtaPhiRegion >())
const EcalRecHitCollection * getCollection(edm::Event &evt, edm::EDGetTokenT< EcalRecHitCollection > &hitToken)
std::vector< BasicCluster > BasicClusterCollection
collection of BasicCluster objects
T get() const
Definition: EventSetup.h:63
def move(src, dest)
Definition: eostools.py:510
bool EgammaHLTIslandClusterProducer::counterExceeded ( ) const
inlineprivate
void EgammaHLTIslandClusterProducer::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 90 of file EgammaHLTIslandClusterProducer.cc.

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

90  {
91 
93  desc.add<std::string>("VerbosityLevel", "ERROR");
94  desc.add<bool>("doBarrel", true);
95  desc.add<bool>("doEndcaps", true);
96  desc.add<bool>("doIsolated", true);
97  desc.add<edm::InputTag>("barrelHitProducer", edm::InputTag("islandEndcapBasicClusters", "EcalRecHitsEB"));
98  desc.add<edm::InputTag>("endcapHitProducer", edm::InputTag("islandEndcapBasicClusters", "EcalRecHitsEB"));
99  desc.add<std::string>("barrelClusterCollection", "islandBarrelBasicClusters");
100  desc.add<std::string>("endcapClusterCollection", "islandEndcapBasicClusters");
101  desc.add<double>("IslandBarrelSeedThr", 0.5);
102  desc.add<double>("IslandEndcapSeedThr", 0.18);
103  desc.add<edm::InputTag>("l1TagIsolated", edm::InputTag("l1extraParticles","Isolated"));
104  desc.add<edm::InputTag>("l1TagNonIsolated", edm::InputTag("l1extraParticles","NonIsolated"));
105  desc.add<double>("l1LowerThr", 0.0);
106  desc.add<double>("l1UpperThr", 9999.0);
107  desc.add<double>("l1LowerThrIgnoreIsolation", 9999.0);
108  desc.add<double>("regionEtaMargin", 0.3);
109  desc.add<double>("regionPhiMargin", 0.4);
110  //desc.add<edm::ParameterSet>("posCalcParameters"), edm::ParameterSet());
111  descriptions.add("hltEgammaHLTIslandClusterProducer", desc);
112 }
ParameterDescriptionBase * add(U const &iLabel, T const &value)
void add(std::string const &label, ParameterSetDescription const &psetDescription)
const EcalRecHitCollection * EgammaHLTIslandClusterProducer::getCollection ( edm::Event evt,
edm::EDGetTokenT< EcalRecHitCollection > &  hitToken 
)
private

Definition at line 246 of file EgammaHLTIslandClusterProducer.cc.

References gather_cfg::cout, edm::Event::getByToken(), edm::HandleBase::isValid(), and edm::Handle< T >::product().

Referenced by clusterizeECALPart().

248 {
250  evt.getByToken(hitToken_, rhcHandle);
251  if (!(rhcHandle.isValid()))
252  {
253  std::cout << "could not get a handle on the EcalRecHitCollection!" << std::endl;
254  edm::LogError("EgammaHLTIslandClusterProducerError") << "Error! can't get the product ";
255  return nullptr;
256  }
257  return rhcHandle.product();
258 }
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:579
bool isValid() const
Definition: HandleBase.h:74
T const * product() const
Definition: Handle.h:81
void EgammaHLTIslandClusterProducer::produce ( edm::Event evt,
const edm::EventSetup es 
)
override

Definition at line 114 of file EgammaHLTIslandClusterProducer.cc.

References IslandClusterAlgo::barrel, barrelClusterCollection_, barrelHitToken_, clusterizeECALPart(), doBarrel_, doEndcaps_, doIsolated_, L1CaloGeometry::emJetPhiBinHighEdge(), L1CaloGeometry::emJetPhiBinLowEdge(), IslandClusterAlgo::endcap, endcapClusterCollection_, endcapHitToken_, L1CaloGeometry::etaBinHighEdge(), L1CaloGeometry::etaBinLowEdge(), edm::EventSetup::get(), edm::Event::getByToken(), l1LowerThr_, l1LowerThrIgnoreIsolation_, l1TagIsolated_, l1TagNonIsolated_, l1UpperThr_, nEvt_, regionEtaMargin_, and regionPhiMargin_.

114  {
115 
116  //Get the L1 EM Particle Collection
118  if(doIsolated_)
119  evt.getByToken(l1TagIsolated_, emIsolColl);
120 
121  //Get the L1 EM Particle Collection
123  evt.getByToken(l1TagNonIsolated_, emNonIsolColl);
124 
125  // Get the CaloGeometry
126  edm::ESHandle<L1CaloGeometry> l1CaloGeom ;
127  es.get<L1CaloGeometryRecord>().get(l1CaloGeom) ;
128 
129  std::vector<RectangularEtaPhiRegion> barrelRegions;
130  std::vector<RectangularEtaPhiRegion> endcapRegions;
131 
132  if(doIsolated_) {
133  for( l1extra::L1EmParticleCollection::const_iterator emItr = emIsolColl->begin(); emItr != emIsolColl->end() ;++emItr ){
134 
135  if (emItr->et() > l1LowerThr_ && emItr->et() < l1UpperThr_) {
136 
137  // Access the GCT hardware object corresponding to the L1Extra EM object.
138  int etaIndex = emItr->gctEmCand()->etaIndex() ;
139 
140  int phiIndex = emItr->gctEmCand()->phiIndex() ;
141  // Use the L1CaloGeometry to find the eta, phi bin boundaries.
142  double etaLow = l1CaloGeom->etaBinLowEdge( etaIndex ) ;
143  double etaHigh = l1CaloGeom->etaBinHighEdge( etaIndex ) ;
144  double phiLow = l1CaloGeom->emJetPhiBinLowEdge( phiIndex ) ;
145  double phiHigh = l1CaloGeom->emJetPhiBinHighEdge( phiIndex ) ;
146 
147  //Attention isForward does not work
148  int isforw=0;
149  int isbarl=0;
150  if((float)(etaHigh)>1.479 || (float)(etaLow)<-1.479) isforw=1;
151  if(((float)(etaLow)>-1.479 && (float)(etaLow)<1.479) ||
152  ((float)(etaHigh)>-1.479 && (float)(etaHigh)<1.479)) isbarl=1;
153 
154  //std::cout<<"Island etaindex "<<etaIndex<<" low hig : "<<etaLow<<" "<<etaHigh<<" phi low hig" <<phiLow<<" " << phiHigh<<" isforw "<<emItr->gctEmCand()->regionId().isForward()<<" isforwnew" <<isforw<< std::endl;
155 
156  etaLow -= regionEtaMargin_;
157  etaHigh += regionEtaMargin_;
158  phiLow -= regionPhiMargin_;
159  phiHigh += regionPhiMargin_;
160 
161  //if (emItr->gctEmCand()->regionId().isForward()) {
162  if (isforw) {
163  if (etaHigh>-1.479 && etaHigh<1.479) etaHigh=-1.479;
164  if ( etaLow>-1.479 && etaLow<1.479) etaLow=1.479;
165  RectangularEtaPhiRegion region(etaLow,etaHigh,phiLow,phiHigh);
166  endcapRegions.push_back(region);
167  }
168  if (isbarl) {
169  if (etaHigh>1.479) etaHigh=1.479;
170  if (etaLow<-1.479) etaLow=-1.479;
171  RectangularEtaPhiRegion region(etaLow,etaHigh,phiLow,phiHigh);
172  barrelRegions.push_back(region);
173  }
174  RectangularEtaPhiRegion region(etaLow,etaHigh,phiLow,phiHigh);
175 
176  }
177  }
178  }
179 
180 
182  for( l1extra::L1EmParticleCollection::const_iterator emItr = emNonIsolColl->begin(); emItr != emNonIsolColl->end() ;++emItr ){
183 
184  if(doIsolated_&&emItr->et()<l1LowerThrIgnoreIsolation_) continue;
185 
186  if (emItr->et() > l1LowerThr_ && emItr->et() < l1UpperThr_) {
187 
188  // Access the GCT hardware object corresponding to the L1Extra EM object.
189  int etaIndex = emItr->gctEmCand()->etaIndex() ;
190 
191 
192  int phiIndex = emItr->gctEmCand()->phiIndex() ;
193  // Use the L1CaloGeometry to find the eta, phi bin boundaries.
194  double etaLow = l1CaloGeom->etaBinLowEdge( etaIndex ) ;
195  double etaHigh = l1CaloGeom->etaBinHighEdge( etaIndex ) ;
196  double phiLow = l1CaloGeom->emJetPhiBinLowEdge( phiIndex ) ;
197  double phiHigh = l1CaloGeom->emJetPhiBinHighEdge( phiIndex ) ;
198 
199 
200  int isforw=0;
201  int isbarl=0;
202  if((float)(etaHigh)>1.479 || (float)(etaLow)<-1.479) isforw=1;
203  if(((float)(etaLow)>-1.479 && (float)(etaLow)<1.479) ||
204  ((float)(etaHigh)>-1.479 && (float)(etaHigh)<1.479)) isbarl=1;
205 
206  //std::cout<<"Island etaindex "<<etaIndex<<" low hig : "<<etaLow<<" "<<etaHigh<<" phi low hig" <<phiLow<<" " << phiHigh<<" isforw "<<emItr->gctEmCand()->regionId().isForward()<<" isforwnew" <<isforw<< std::endl;
207 
208  etaLow -= regionEtaMargin_;
209  etaHigh += regionEtaMargin_;
210  phiLow -= regionPhiMargin_;
211  phiHigh += regionPhiMargin_;
212 
213  //if (emItr->gctEmCand()->regionId().isForward()) {
214  if (isforw) {
215  if (etaHigh>-1.479 && etaHigh<1.479) etaHigh=-1.479;
216  if ( etaLow>-1.479 && etaLow<1.479) etaLow=1.479;
217  RectangularEtaPhiRegion region(etaLow,etaHigh,phiLow,phiHigh);
218  endcapRegions.push_back(region);
219  }
220  if (isbarl) {
221  if (etaHigh>1.479) etaHigh=1.479;
222  if (etaLow<-1.479) etaLow=-1.479;
223  RectangularEtaPhiRegion region(etaLow,etaHigh,phiLow,phiHigh);
224  barrelRegions.push_back(region);
225  }
226 
227  }
228  }
229  }
230 
231  if (doEndcaps_
232  //&&endcapRegions.size()!=0
233  ) {
234 
236  }
237  if (doBarrel_
238  //&& barrelRegions.size()!=0
239  ) {
241  }
242  nEvt_++;
243 }
edm::EDGetTokenT< EcalRecHitCollection > endcapHitToken_
double etaBinHighEdge(unsigned int etaIndex, bool central=true) const
double etaBinLowEdge(unsigned int etaIndex, bool central=true) const
edm::EDGetTokenT< l1extra::L1EmParticleCollection > l1TagNonIsolated_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:579
edm::EDGetTokenT< l1extra::L1EmParticleCollection > l1TagIsolated_
void clusterizeECALPart(edm::Event &evt, const edm::EventSetup &es, edm::EDGetTokenT< EcalRecHitCollection > &hitToken, const std::string &clusterCollection, const std::vector< RectangularEtaPhiRegion > &regions, const IslandClusterAlgo::EcalPart &ecalPart)
T get() const
Definition: EventSetup.h:63
double emJetPhiBinLowEdge(unsigned int phiIndex) const
edm::EDGetTokenT< EcalRecHitCollection > barrelHitToken_
double emJetPhiBinHighEdge(unsigned int phiIndex) const

Member Data Documentation

std::string EgammaHLTIslandClusterProducer::barrelClusterCollection_
private
edm::InputTag EgammaHLTIslandClusterProducer::barrelHitCollection_
private

Definition at line 43 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer().

edm::EDGetTokenT<EcalRecHitCollection> EgammaHLTIslandClusterProducer::barrelHitToken_
private

Definition at line 45 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

bool EgammaHLTIslandClusterProducer::doBarrel_
private

Definition at line 39 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

bool EgammaHLTIslandClusterProducer::doEndcaps_
private

Definition at line 40 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

bool EgammaHLTIslandClusterProducer::doIsolated_
private

Definition at line 41 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

std::string EgammaHLTIslandClusterProducer::endcapClusterCollection_
private
edm::InputTag EgammaHLTIslandClusterProducer::endcapHitCollection_
private

Definition at line 44 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer().

edm::EDGetTokenT<EcalRecHitCollection> EgammaHLTIslandClusterProducer::endcapHitToken_
private

Definition at line 46 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

IslandClusterAlgo* EgammaHLTIslandClusterProducer::island_p
private
double EgammaHLTIslandClusterProducer::l1LowerThr_
private

Definition at line 53 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

double EgammaHLTIslandClusterProducer::l1LowerThrIgnoreIsolation_
private

Definition at line 55 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

edm::EDGetTokenT<l1extra::L1EmParticleCollection> EgammaHLTIslandClusterProducer::l1TagIsolated_
private

Definition at line 51 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

edm::EDGetTokenT<l1extra::L1EmParticleCollection> EgammaHLTIslandClusterProducer::l1TagNonIsolated_
private

Definition at line 52 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

double EgammaHLTIslandClusterProducer::l1UpperThr_
private

Definition at line 54 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

int EgammaHLTIslandClusterProducer::nEvt_
private

Definition at line 35 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

int EgammaHLTIslandClusterProducer::nMaxPrintout_
private

Definition at line 34 of file EgammaHLTIslandClusterProducer.h.

PositionCalc EgammaHLTIslandClusterProducer::posCalculator_
private

Definition at line 60 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer().

double EgammaHLTIslandClusterProducer::regionEtaMargin_
private

Definition at line 57 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

double EgammaHLTIslandClusterProducer::regionPhiMargin_
private

Definition at line 58 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer(), and produce().

IslandClusterAlgo::VerbosityLevel EgammaHLTIslandClusterProducer::verbosity
private

Definition at line 37 of file EgammaHLTIslandClusterProducer.h.

Referenced by EgammaHLTIslandClusterProducer().