#include <InterestingDetIdCollectionProducer.h>

Inheritance diagram for InterestingDetIdCollectionProducer:

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

	InterestingDetIdCollectionProducer (const edm::ParameterSet &)
	ctor More...

void	produce (edm::Event &, const edm::EventSetup &) override
	producer More...

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
edm::EDGetTokenT< reco::BasicClusterCollection >	basicClustersToken_

const CaloTopology *	caloTopology_

std::string	interestingDetIdCollection_

bool	keepNextToBoundary_

bool	keepNextToDead_

int	minimalEtaSize_

int	minimalPhiSize_

edm::EDGetTokenT< EcalRecHitCollection >	recHitsToken_

const EcalSeverityLevelAlgo *	severity_

int	severityLevel_

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

Original author: Paolo Meridiani PH/CMG

Make a collection of detids to be kept tipically in a AOD rechit collection

The following classes of "interesting id" are considered

1.in a region around  the seed of the cluster collection specified
  by paramter basicClusters. The size of the region is specified by
  minimalEtaSize_, minimalPhiSize_

2. if the severity of the hit is >= severityLevel_
   If severityLevel=0 this class is ignored

3. Channels next to dead ones,  keepNextToDead_ is true
4. Channels next to the EB/EE transition if keepNextToBoundary_ is true

Definition at line 47 of file InterestingDetIdCollectionProducer.h.

Constructor & Destructor Documentation

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

explicit

ctor

Definition at line 23 of file InterestingDetIdCollectionProducer.cc.

References basicClustersToken_, edm::ParameterSet::getParameter(), interestingDetIdCollection_, keepNextToBoundary_, keepNextToDead_, minimalEtaSize_, minimalPhiSize_, recHitsToken_, severityLevel_, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
 
   recHitsToken_ = 
       consumes<EcalRecHitCollection>(iConfig.getParameter< edm::InputTag > ("recHitsLabel"));
   basicClustersToken_ = 
       consumes<reco::BasicClusterCollection>(iConfig.getParameter< edm::InputTag >("basicClustersLabel"));
 
   interestingDetIdCollection_ = iConfig.getParameter<std::string>("interestingDetIdCollection");
   
   minimalEtaSize_ = iConfig.getParameter<int> ("etaSize");
   minimalPhiSize_ = iConfig.getParameter<int> ("phiSize");
   if ( minimalPhiSize_ % 2 == 0 ||  minimalEtaSize_ % 2 == 0)
     edm::LogError("InterestingDetIdCollectionProducerError") << "Size of eta/phi should be odd numbers";
  
    //register your products
   produces< DetIdCollection > (interestingDetIdCollection_) ;
 
   severityLevel_  = iConfig.getParameter<int>("severityLevel");
   keepNextToDead_ = iConfig.getParameter<bool>("keepNextToDead");
   keepNextToBoundary_ = iConfig.getParameter<bool>("keepNextToBoundary");
 }

Member Function Documentation

void InterestingDetIdCollectionProducer::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

final

Definition at line 47 of file InterestingDetIdCollectionProducer.cc.

References caloTopology_, edm::EventSetup::get(), edm::ESHandle< T >::product(), and severity_.

 {
   edm::ESHandle<CaloTopology> theCaloTopology;
   iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
   caloTopology_ = &(*theCaloTopology); 
 
   edm::ESHandle<EcalSeverityLevelAlgo> sevLv;
   iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevLv);
   severity_ = sevLv.product();
 }

void InterestingDetIdCollectionProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

producer

Definition at line 60 of file InterestingDetIdCollectionProducer.cc.

References funct::abs(), basicClustersToken_, edm::SortedCollection< T, SORT >::begin(), caloTopology_, DetId::det(), EcalBarrel, EcalEndcap, photonPostprocessing_cfi::eMax, edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), edm::Event::getByToken(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), interestingDetIdCollection_, EcalTools::isNextToDead(), EEDetId::isNextToRingBoundary(), keepNextToBoundary_, keepNextToDead_, EcalRecHit::kTPSaturated, minimalEtaSize_, minimalPhiSize_, DetId::null(), edm::Event::put(), recHitsToken_, severity_, EcalSeverityLevelAlgo::severityLevel(), severityLevel_, DetId::subdetId(), ecaldqm::topology(), and tier0::unique().

 {
    using namespace edm;
    using namespace std;
 
    // take BasicClusters
   Handle<reco::BasicClusterCollection> pClusters;
   iEvent.getByToken(basicClustersToken_, pClusters);
   
   // take EcalRecHits
   Handle<EcalRecHitCollection> recHitsHandle;
   iEvent.getByToken(recHitsToken_,recHitsHandle);
 
   //Create empty output collections
   std::vector<DetId> indexToStore;
   indexToStore.reserve(1000);
 
   reco::BasicClusterCollection::const_iterator clusIt;
 
   std::vector<DetId> xtalsToStore;
   xtalsToStore.reserve(50);
   for (clusIt=pClusters->begin(); clusIt!=pClusters->end(); clusIt++) {
     const std::vector<std::pair<DetId,float> > &clusterDetIds = clusIt->hitsAndFractions();
     for (const auto &detidpair : clusterDetIds) {
       indexToStore.push_back(detidpair.first);
     }
     
     //below checks and additional code are only relevant for EB/EE, not for ES
     if (clusterDetIds.front().first.subdetId()==EcalBarrel || clusterDetIds.front().first.subdetId()==EcalEndcap) {
       std::vector<std::pair<DetId,float> >::const_iterator posCurrent;
       
       float eMax=0.;
       DetId eMaxId(0);        
 
       EcalRecHit testEcalRecHit;
       
       for(posCurrent = clusterDetIds.begin(); posCurrent != clusterDetIds.end(); posCurrent++)
         {
           EcalRecHitCollection::const_iterator itt = recHitsHandle->find((*posCurrent).first);
           if ((!((*posCurrent).first.null())) && (itt != recHitsHandle->end()) && ((*itt).energy() > eMax) )
             {
               eMax = (*itt).energy();
               eMaxId = (*itt).id();
             }
         }
       
       if (eMaxId.null())
       continue;
       
       const CaloSubdetectorTopology* topology  = caloTopology_->getSubdetectorTopology(eMaxId.det(),eMaxId.subdetId());
 
       xtalsToStore=topology->getWindow(eMaxId,minimalEtaSize_,minimalPhiSize_);
       
       for (const auto &detid : xtalsToStore) {
         indexToStore.push_back(detid);
       }
       
     }
   }
 
 
   if (severityLevel_>=0 || keepNextToDead_ || keepNextToBoundary_) {
     for (EcalRecHitCollection::const_iterator it = recHitsHandle->begin(); it != recHitsHandle->end(); ++it) {
       // also add recHits of dead TT if the corresponding TP is saturated
       if ( it->checkFlag(EcalRecHit::kTPSaturated) ) {
         indexToStore.push_back(it->id());
       }
       // add hits for severities above a threshold
       if ( severityLevel_>=0 && 
           severity_->severityLevel(*it) >=severityLevel_){
         
         indexToStore.push_back(it->id());
       } 
       if (keepNextToDead_) {
         // also keep channels next to dead ones
         if (EcalTools::isNextToDead(it->id(), iSetup)) {
           indexToStore.push_back(it->id());
         }
       } 
 
       if (keepNextToBoundary_){
         // keep channels around EB/EE boundary
         if (it->id().subdetId() == EcalBarrel){
           EBDetId ebid(it->id());
           if (abs(ebid.ieta())== 85)
             indexToStore.push_back(it->id());
         } else {
       
           if (EEDetId::isNextToRingBoundary(it->id()))
             indexToStore.push_back(it->id());
         }
 
       }
       
     }
   }
 
   //unify the vector
   std::sort(indexToStore.begin(),indexToStore.end());
   std::unique(indexToStore.begin(),indexToStore.end());
   
  
   iEvent.put(std::make_unique<DetIdCollection>(indexToStore), interestingDetIdCollection_);
 
 }