#include <ElectronSeedProducer.h>

Inheritance diagram for ElectronSeedProducer:

Public Member Functions
void	beginRun (edm::Run const &, edm::EventSetup const &) final

	ElectronSeedProducer (const edm::ParameterSet &)

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	produce (edm::Event &, const edm::EventSetup &) final

	~ElectronSeedProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

Private Member Functions
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)

void	filterSeeds (edm::Event &e, const edm::EventSetup &setup, reco::SuperClusterRefVector &sclRefs)

Private Attributes
bool	allowHGCal_

bool	applyHOverECut_

bool	applySigmaIEtaIEtaCut_

edm::EDGetTokenT< reco::BeamSpot >	beamSpotTag_

edm::ESHandle< CaloGeometry >	caloGeom_

unsigned long long	caloGeomCacheId_

edm::ESHandle< CaloTopology >	caloTopo_

unsigned long long	caloTopoCacheId_

edm::EDGetTokenT< EcalRecHitCollection >	ebRecHitCollection_

edm::EDGetTokenT< EcalRecHitCollection >	eeRecHitCollection_

bool	fromTrackerSeeds_

ElectronHcalHelper *	hcalHelper_

std::unique_ptr< hgcal::ClusterTools >	hgcClusterTools_

std::vector< edm::EDGetTokenT< TrajectorySeedCollection > >	initialSeeds_

ElectronSeedGenerator *	matcher_

double	maxHBarrel_

double	maxHEndcaps_

double	maxHOverEBarrel_

double	maxHOverEEndcaps_

double	maxSigmaIEtaIEtaBarrel_

double	maxSigmaIEtaIEtaEndcaps_

bool	prefilteredSeeds_

double	SCEtCut_

std::unique_ptr< SeedFilter >	seedFilter_

edm::EDGetTokenT< reco::SuperClusterCollection >	superClusters_ [2]

TrajectorySeedCollection *	theInitialSeedColl

std::vector< const TrajectorySeedCollection * >	theInitialSeedColls

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: MeasurementEstimator for Pixel Barrel, ported from ORCA Class defining the search area in the barrel in the pixel match Implementation: <Notes on="" implementation>="">

Description: MeasurementEstimator for Pixel Barrel, ported from ORCA

Implementation: <Notes on="" implementation>="">

Description: EDProducer of ElectronSeed objects

Implementation: <Notes on="" implementation>="">

Definition at line 42 of file ElectronSeedProducer.h.

Constructor & Destructor Documentation

ElectronSeedProducer::ElectronSeedProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 61 of file ElectronSeedProducer.cc.

References allowHGCal_, applyHOverECut_, applySigmaIEtaIEtaCut_, beamSpotTag_, ebRecHitCollection_, eeRecHitCollection_, edm::ParameterSet::empty(), fromTrackerSeeds_, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), hcalHelper_, ElectronHcalHelper::Configuration::hcalTowers, hgcClusterTools_, ElectronHcalHelper::Configuration::hOverEConeSize, ElectronHcalHelper::Configuration::hOverEPtMin, initialSeeds_, matcher_, maxHBarrel_, maxHEndcaps_, maxHOverEBarrel_, maxHOverEEndcaps_, maxSigmaIEtaIEtaBarrel_, maxSigmaIEtaIEtaEndcaps_, prefilteredSeeds_, SCEtCut_, seedFilter_, superClusters_, GlobalPosition_Frontier_DevDB_cff::tag, SeedFilter::Tokens::token_bs, ElectronSeedGenerator::Tokens::token_bs, SeedFilter::Tokens::token_vtx, funct::true, ElectronHcalHelper::Configuration::useTowers, and edm::vector_transform().

  : applyHOverECut_(true), hcalHelper_(nullptr),
    caloGeom_(nullptr), caloGeomCacheId_(0), caloTopo_(nullptr), caloTopoCacheId_(0)
  {
   auto const& conf = iConfig.getParameter<edm::ParameterSet>("SeedConfiguration") ;
 
   auto legacyConfSeeds = conf.getParameter<edm::InputTag>("initialSeeds");
   if (legacyConfSeeds.label().empty())
     {//new format
       initialSeeds_ = edm::vector_transform(conf.getParameter<std::vector<edm::InputTag> >( "initialSeedsVector" ), 
                                             [this](edm::InputTag const & tag){return consumes<TrajectorySeedCollection >(tag);});
     } 
   else
     {
       initialSeeds_ = {consumes<TrajectorySeedCollection>(conf.getParameter<edm::InputTag>("initialSeeds"))};
     }
 
   SCEtCut_ = conf.getParameter<double>("SCEtCut");
   fromTrackerSeeds_ = conf.getParameter<bool>("fromTrackerSeeds") ;
   prefilteredSeeds_ = conf.getParameter<bool>("preFilteredSeeds") ;
 
   auto theconsumes = consumesCollector();
 
   // new beamSpot tag
   beamSpotTag_ = consumes<reco::BeamSpot>(conf.getParameter<edm::InputTag>("beamSpot")); 
 
   // for H/E
   applyHOverECut_ = conf.getParameter<bool>("applyHOverECut") ;
   if (applyHOverECut_)
    {
      ElectronHcalHelper::Configuration hcalCfg ;
      hcalCfg.hOverEConeSize = conf.getParameter<double>("hOverEConeSize") ;
      if (hcalCfg.hOverEConeSize>0)
        {
      hcalCfg.useTowers = true ;
      hcalCfg.hcalTowers = 
        consumes<CaloTowerCollection>(conf.getParameter<edm::InputTag>("hcalTowers")) ;
      hcalCfg.hOverEPtMin = conf.getParameter<double>("hOverEPtMin") ;
        }
      hcalHelper_ = new ElectronHcalHelper(hcalCfg) ;
      
      allowHGCal_ = conf.getParameter<bool>("allowHGCal");
      if( allowHGCal_ ) {
        const edm::ParameterSet& hgcCfg = conf.getParameterSet("HGCalConfig");
        hgcClusterTools_.reset( new hgcal::ClusterTools(hgcCfg, theconsumes) );
      }
      
      maxHOverEBarrel_=conf.getParameter<double>("maxHOverEBarrel") ;
      maxHOverEEndcaps_=conf.getParameter<double>("maxHOverEEndcaps") ;
      maxHBarrel_=conf.getParameter<double>("maxHBarrel") ;
      maxHEndcaps_=conf.getParameter<double>("maxHEndcaps") ;
    }
 
   applySigmaIEtaIEtaCut_ = conf.getParameter<bool>("applySigmaIEtaIEtaCut");
 
   // apply sigma_ieta_ieta cut
   if (applySigmaIEtaIEtaCut_ == true)
     {
       maxSigmaIEtaIEtaBarrel_ = conf.getParameter<double>("maxSigmaIEtaIEtaBarrel");
       maxSigmaIEtaIEtaEndcaps_ = conf.getParameter<double>("maxSigmaIEtaIEtaEndcaps");
     }
 
   edm::ParameterSet rpset = conf.getParameter<edm::ParameterSet>("RegionPSet");
 
   ElectronSeedGenerator::Tokens esg_tokens;
   esg_tokens.token_bs = beamSpotTag_;
   esg_tokens.token_vtx = mayConsume<reco::VertexCollection>(conf.getParameter<edm::InputTag>("vertices"));
   esg_tokens.token_measTrkEvt= consumes<MeasurementTrackerEvent>(conf.getParameter<edm::InputTag>("measurementTrackerEvent"));
 
   matcher_ = new ElectronSeedGenerator(conf,esg_tokens) ;
 
   //  get collections from config
   if (applySigmaIEtaIEtaCut_ == true) {
     ebRecHitCollection_ = consumes<EcalRecHitCollection> (iConfig.getParameter<edm::InputTag>("ebRecHitCollection"));
     eeRecHitCollection_ = consumes<EcalRecHitCollection> (iConfig.getParameter<edm::InputTag>("eeRecHitCollection"));
   }
 
   superClusters_[0]=
     consumes<reco::SuperClusterCollection>(iConfig.getParameter<edm::InputTag>("barrelSuperClusters")) ;
   superClusters_[1]=
     consumes<reco::SuperClusterCollection>(iConfig.getParameter<edm::InputTag>("endcapSuperClusters")) ;
 
   // Construction of SeedFilter was in beginRun() with the comment
   // below, but it has to be done here because of ConsumesCollector
   //
   // FIXME: because of a bug presumably in tracker seeding,
   // perhaps in CombinedHitPairGenerator, badly caching some EventSetup product,
   // we must redo the SeedFilter for each run.
   if (prefilteredSeeds_) {
     SeedFilter::Tokens sf_tokens;
     sf_tokens.token_bs  = beamSpotTag_;
     sf_tokens.token_vtx = consumes<std::vector<reco::Vertex> >(rpset.getParameter<edm::InputTag> ("VertexProducer"));
 
     edm::ConsumesCollector iC = consumesCollector();
     seedFilter_.reset(new SeedFilter(conf, sf_tokens, iC));
   }
 
   //register your products
   produces<ElectronSeedCollection>() ;
 }

ElectronSeedProducer::~ElectronSeedProducer ( )

override

Definition at line 168 of file ElectronSeedProducer.cc.

References hcalHelper_, and matcher_.

  {
   delete hcalHelper_ ;
   delete matcher_ ;
  }

Member Function Documentation

void ElectronSeedProducer::beginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

final

Definition at line 162 of file ElectronSeedProducer.cc.

163 {}

void ElectronSeedProducer::endRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

final

Definition at line 165 of file ElectronSeedProducer.cc.

166 {}

void ElectronSeedProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 348 of file ElectronSeedProducer.cc.

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

                                                                                  {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("endcapSuperClusters",edm::InputTag("particleFlowSuperClusterECAL","particleFlowSuperClusterECALEndcapWithPreshower"));
   {
     edm::ParameterSetDescription psd0, psd1, psd2, psd3, psd4;
     psd1.add<unsigned int>("maxElement", 0);
     psd1.add<std::string>("ComponentName", std::string("StandardHitPairGenerator"));
     psd1.addUntracked<int>("useOnDemandTracker", 0);
     psd1.add<edm::InputTag>("SeedingLayers", edm::InputTag("hltMixedLayerPairs"));
     psd0.add<edm::ParameterSetDescription>("OrderedHitsFactoryPSet", psd1);
 
     psd2.add<double>("deltaPhiRegion", 0.4);
     psd2.add<double>("originHalfLength", 15.0);
     psd2.add<bool>("useZInVertex", true);
     psd2.add<double>("deltaEtaRegion", 0.1);
     psd2.add<double>("ptMin", 1.5 );
     psd2.add<double>("originRadius", 0.2);
     psd2.add<edm::InputTag>("VertexProducer", edm::InputTag("dummyVertices"));
     psd0.add<edm::ParameterSetDescription>("RegionPSet", psd2);
      
     psd0.add<double>("PhiMax2B",0.002);
     psd0.add<double>("hOverEPtMin",0.0);
     psd0.add<double>("PhiMax2F",0.003);
     psd0.add<bool>("searchInTIDTEC",true);
     psd0.add<double>("pPhiMax1",0.125);
     psd0.add<double>("HighPtThreshold",35.0);
     psd0.add<double>("r2MinF",-0.15);
     psd0.add<double>("maxHBarrel",0.0);
     psd0.add<double>("DeltaPhi1Low",0.23);
     psd0.add<double>("DeltaPhi1High",0.08);
     psd0.add<double>("ePhiMin1",-0.125);
     psd0.add<edm::InputTag>("hcalTowers",edm::InputTag("towerMaker"));
     psd0.add<double>("LowPtThreshold",5.0);
     psd0.add<double>("maxHOverEBarrel",0.15);
     psd0.add<double>("maxSigmaIEtaIEtaBarrel", 0.5);
     psd0.add<double>("maxSigmaIEtaIEtaEndcaps", 0.5);
     psd0.add<bool>("dynamicPhiRoad",true);
     psd0.add<double>("ePhiMax1",0.075);
     psd0.add<std::string>("measurementTrackerName","");
     psd0.add<double>("SizeWindowENeg",0.675);
     psd0.add<double>("nSigmasDeltaZ1",5.0);
     psd0.add<double>("rMaxI",0.2);
     psd0.add<double>("maxHEndcaps",0.0);
     psd0.add<bool>("preFilteredSeeds",false);
     psd0.add<double>("r2MaxF",0.15);
     psd0.add<double>("hOverEConeSize",0.15);
     psd0.add<double>("pPhiMin1",-0.075);
     psd0.add<edm::InputTag>("initialSeeds",edm::InputTag(""));//keep for be compatibility
     psd0.add<std::vector<edm::InputTag>>("initialSeedsVector",{{"newCombinedSeeds"}});
     psd0.add<double>("deltaZ1WithVertex",25.0);
     psd0.add<double>("SCEtCut",0.0);
     psd0.add<double>("z2MaxB",0.09);
     psd0.add<bool>("fromTrackerSeeds",true);
     psd0.add<edm::InputTag>("hcalRecHits",edm::InputTag("hbhereco"));
     psd0.add<double>("z2MinB",-0.09);
     psd0.add<double>("rMinI",-0.2);
     psd0.add<double>("maxHOverEEndcaps",0.15);
     psd0.add<double>("hOverEHBMinE",0.7);
     psd0.add<bool>("useRecoVertex",false);
     psd0.add<edm::InputTag>("beamSpot",edm::InputTag("offlineBeamSpot"));
     psd0.add<edm::InputTag>("measurementTrackerEvent",edm::InputTag("MeasurementTrackerEvent"));
     psd0.add<edm::InputTag>("vertices",edm::InputTag("offlinePrimaryVerticesWithBS"));
     psd0.add<bool>("applyHOverECut",true);
     psd0.add<edm::InputTag>("ebRecHitCollection", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
     psd0.add<edm::InputTag>("eeRecHitCollection", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
     psd0.add<bool>("applySigmaIEtaIEtaCut", false);
     psd0.add<double>("DeltaPhi2F",0.012);
     psd0.add<double>("PhiMin2F",-0.003);
     psd0.add<double>("hOverEHFMinE",0.8);
     psd0.add<double>("DeltaPhi2B",0.008);
     psd0.add<double>("PhiMin2B",-0.002);
     psd0.add<bool>("allowHGCal",false);
     
     psd3.add<std::string>("ComponentName",std::string("SeedFromConsecutiveHitsCreator"));
     psd3.add<std::string>("propagator",std::string("PropagatorWithMaterial"));
     psd3.add<double>("SeedMomentumForBOFF",5.0);
     psd3.add<double>("OriginTransverseErrorMultiplier",1.0);
     psd3.add<double>("MinOneOverPtError",1.0);
     psd3.add<std::string>("magneticField",std::string(""));
     psd3.add<std::string>("TTRHBuilder",std::string("WithTrackAngle"));
     psd3.add<bool>("forceKinematicWithRegionDirection",false);
 
     psd4.add<edm::InputTag>("HGCEEInput",edm::InputTag("HGCalRecHit","HGCEERecHits"));
     psd4.add<edm::InputTag>("HGCFHInput",edm::InputTag("HGCalRecHit","HGCHEFRecHits"));
     psd4.add<edm::InputTag>("HGCBHInput",edm::InputTag("HGCalRecHit","HGCHEBRecHits"));
     
     psd0.add<edm::ParameterSetDescription>("HGCalConfig",psd4);
 
     psd0.add<edm::ParameterSetDescription>("SeedCreatorPSet",psd3);
 
     desc.add<edm::ParameterSetDescription>("SeedConfiguration",psd0);
   }
   desc.add<edm::InputTag>("barrelSuperClusters",edm::InputTag("particleFlowSuperClusterECAL","particleFlowSuperClusterECALBarrel"));
   descriptions.add("ecalDrivenElectronSeeds",desc);
 }

void ElectronSeedProducer::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
	)

private

Definition at line 270 of file ElectronSeedProducer.cc.

References allowHGCal_, applyHOverECut_, applySigmaIEtaIEtaCut_, gamEcalExtractorBlocks_cff::detector, ebRecHitCollection_, EcalBarrel, EcalEndcap, eeRecHitCollection_, reco::CaloCluster::energy(), PVValHelper::eta, funct::false, filterSeeds(), ElectronHcalHelper::hcalESumDepth1(), ElectronHcalHelper::hcalESumDepth2(), hcalHelper_, hgcClusterTools_, reco::CaloCluster::hitsAndFractions(), mps_fire::i, infinity, EcalTools::isHGCalDet(), edm::isNotFinite(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::localCovariances(), LogDebug, maxHBarrel_, maxHEndcaps_, maxHOverEBarrel_, maxHOverEEndcaps_, reco::BeamSpot::position(), reco::CaloCluster::position(), edm::RefVector< C, T, F >::push_back(), SCEtCut_, reco::SuperCluster::seed(), edm::RefVector< C, T, F >::size(), and mathSSE::sqrt().

Referenced by produce().

  {
 
    std::vector<float> sigmaIEtaIEtaEB_;
    std::vector<float> sigmaIEtaIEtaEE_;
 
   for (unsigned int i=0;i<superClusters->size();++i)
    {
     const SuperCluster & scl = (*superClusters)[i] ;
     double sclEta = EleRelPoint(scl.position(),bs.position()).eta() ;
     if (scl.energy()/cosh(sclEta)>SCEtCut_)
      {
 //      if ((applyHOverECut_==true)&&((hcalHelper_->hcalESum(scl)/scl.energy()) > maxHOverE_))
 //       { continue ; }
 //      sclRefs.push_back(edm::Ref<reco::SuperClusterCollection>(superClusters,i)) ;
        double had1, had2, had, scle ;
 
        bool HoeVeto = false ;
        if (applyHOverECut_==true)
         {
          had1 = hcalHelper_->hcalESumDepth1(scl);
          had2 = hcalHelper_->hcalESumDepth2(scl);
          had = had1+had2 ;
          scle = scl.energy() ;
          int det_group = scl.seed()->hitsAndFractions()[0].first.det() ;
          int detector = scl.seed()->hitsAndFractions()[0].first.subdetId() ;
          if ( detector==EcalBarrel && (had<maxHBarrel_ || had/scle<maxHOverEBarrel_)) HoeVeto=true;
          else if( !allowHGCal_ && detector==EcalEndcap && (had<maxHEndcaps_ || had/scle<maxHOverEEndcaps_) ) HoeVeto=true;
          else if( allowHGCal_ && EcalTools::isHGCalDet((DetId::Detector)det_group) ) {
            float had_fraction = hgcClusterTools_->getClusterHadronFraction(*(scl.seed()));
            had1 = had_fraction*scl.seed()->energy();
            had2 = 0.;
            HoeVeto= ( had_fraction >= 0.f && had_fraction < maxHOverEEndcaps_ );
          }
          if (HoeVeto)
           {
            sclRefs.push_back(edm::Ref<reco::SuperClusterCollection>(superClusters,i)) ;
            hoe1s.push_back(had1/scle) ;
            hoe2s.push_back(had2/scle) ;
           }
         }
        else
         {
          sclRefs.push_back(edm::Ref<reco::SuperClusterCollection>(superClusters,i)) ;
          hoe1s.push_back(std::numeric_limits<float>::infinity()) ;
          hoe2s.push_back(std::numeric_limits<float>::infinity()) ;
         }
      }
 
     if (applySigmaIEtaIEtaCut_ == true)
       {
     noZS::EcalClusterLazyTools lazyTool_noZS(event, setup, ebRecHitCollection_, eeRecHitCollection_);
     std::vector<float> vCov = lazyTool_noZS.localCovariances(*(scl.seed()));
     int detector = scl.seed()->hitsAndFractions()[0].first.subdetId() ;
     if (detector==EcalBarrel) sigmaIEtaIEtaEB_ .push_back(edm::isNotFinite(vCov[0]) ? 0. : sqrt(vCov[0]));
     if (detector==EcalEndcap) sigmaIEtaIEtaEE_ .push_back(edm::isNotFinite(vCov[0]) ? 0. : sqrt(vCov[0]));
       }
    }
   LogDebug("ElectronSeedProducer")<<"Filtered out "<<sclRefs.size()<<" superclusters from "<<superClusters->size() ;
  }

void ElectronSeedProducer::filterSeeds	(	edm::Event &	e,
		const edm::EventSetup &	setup,
		reco::SuperClusterRefVector &	sclRefs
	)

private

Definition at line 336 of file ElectronSeedProducer.cc.

References mps_fire::i, LogDebug, seedFilter_, edm::RefVector< C, T, F >::size(), theInitialSeedColl, and theInitialSeedColls.

Referenced by filterClusters(), and produce().

  {
   for ( unsigned int i=0 ; i<sclRefs.size() ; ++i )
    {
     seedFilter_->seeds(event,setup,sclRefs[i],theInitialSeedColl) ;
     theInitialSeedColls.push_back(theInitialSeedColl);
     LogDebug("ElectronSeedProducer")<<"Number of Seeds: "<<theInitialSeedColl->size() ;
    }
  }

void ElectronSeedProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	iSetup
	)

final

Definition at line 174 of file ElectronSeedProducer.cc.

References allowHGCal_, beamSpotTag_, edm::eventsetup::EventSetupRecord::cacheIdentifier(), caloGeom_, caloGeomCacheId_, caloTopo_, caloTopoCacheId_, edm::RefToBase< T >::castTo(), ElectronHcalHelper::checkSetup(), fastPrimaryVertexProducer_cfi::clusters, filterClusters(), filterSeeds(), fromTrackerSeeds_, edm::EventSetup::get(), edm::Event::getByToken(), hcalHelper_, hgcClusterTools_, mps_fire::i, edm::Ref< C, T, F >::id(), initialSeeds_, LogDebug, matcher_, eostools::move(), nullptr, prefilteredSeeds_, edm::Handle< T >::product(), edm::Event::put(), ElectronHcalHelper::readEvent(), ElectronSeedGenerator::run(), ElectronSeedGenerator::setupES(), superClusters_, theInitialSeedColl, and theInitialSeedColls.

  {
   LogDebug("ElectronSeedProducer") <<"[ElectronSeedProducer::produce] entering " ;
 
   edm::Handle<reco::BeamSpot> theBeamSpot ;
   e.getByToken(beamSpotTag_,theBeamSpot) ;
 
   if (hcalHelper_)
    {
     hcalHelper_->checkSetup(iSetup) ;
     hcalHelper_->readEvent(e) ;
     if( allowHGCal_ ) {
       hgcClusterTools_->getEventSetup(iSetup);
       hgcClusterTools_->getEvent(e);
     }
    }
 
   // get calo geometry
   if (caloGeomCacheId_!=iSetup.get<CaloGeometryRecord>().cacheIdentifier()) {
     iSetup.get<CaloGeometryRecord>().get(caloGeom_);
     caloGeomCacheId_=iSetup.get<CaloGeometryRecord>().cacheIdentifier();
   }
   if (caloTopoCacheId_!=iSetup.get<CaloTopologyRecord>().cacheIdentifier()){
     caloTopoCacheId_=iSetup.get<CaloTopologyRecord>().cacheIdentifier();
     iSetup.get<CaloTopologyRecord>().get(caloTopo_);
   }
 
   matcher_->setupES(iSetup);
 
   // get initial TrajectorySeeds if necessary
   if (fromTrackerSeeds_)
    {
     if (!prefilteredSeeds_)
      {
        theInitialSeedColls.clear();
        for (auto const& seeds : initialSeeds_)
          {
            edm::Handle<TrajectorySeedCollection> hSeeds;
            e.getByToken(seeds, hSeeds);
            theInitialSeedColls.push_back(hSeeds.product());
          }
        theInitialSeedColl = nullptr;// not needed in this case
 
      }
     else
      { 
        theInitialSeedColls.clear();//reset later
        theInitialSeedColl = new TrajectorySeedCollection ; 
      }
    }
   else
    {
      // not needed in this case
      theInitialSeedColls.clear();
      theInitialSeedColl = nullptr ; 
    }
 
   ElectronSeedCollection * seeds = new ElectronSeedCollection ;
 
   // loop over barrel + endcap
   for (unsigned int i=0; i<2; i++) {
     edm::Handle<SuperClusterCollection> clusters ;
     e.getByToken(superClusters_[i],clusters);
     SuperClusterRefVector clusterRefs ;
     std::vector<float> hoe1s, hoe2s ;
     filterClusters(*theBeamSpot,clusters,/*mhbhe_,*/clusterRefs,hoe1s,hoe2s,e, iSetup);
     if ((fromTrackerSeeds_) && (prefilteredSeeds_))
       { filterSeeds(e,iSetup,clusterRefs) ; }
     matcher_->run(e,iSetup,clusterRefs,hoe1s,hoe2s,theInitialSeedColls,*seeds);
   }
 
   // store the accumulated result
   std::unique_ptr<ElectronSeedCollection> pSeeds(seeds);
   ElectronSeedCollection::iterator is ;
   for ( is=pSeeds->begin() ; is!=pSeeds->end() ; is++ ) {
     edm::RefToBase<CaloCluster> caloCluster = is->caloCluster() ;
     SuperClusterRef superCluster = caloCluster.castTo<SuperClusterRef>() ;
     LogDebug("ElectronSeedProducer")
       << "new seed with "
       << (*is).nHits() << " hits"
       << ", charge " << (*is).getCharge()
       << " and cluster energy " << superCluster->energy()
       << " PID "<<superCluster.id() ;
   }
   e.put(std::move(pSeeds));
   if (fromTrackerSeeds_ && prefilteredSeeds_) delete theInitialSeedColl;
  }