#include <InterestingDetIdFromSuperClusterProducer.h>

Inheritance diagram for InterestingDetIdFromSuperClusterProducer:

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

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

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

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
const CaloTopology *	caloTopology_

std::string	interestingDetIdCollection_

bool	keepNextToBoundary_

bool	keepNextToDead_

int	minimalEtaSize_

int	minimalPhiSize_

edm::EDGetTokenT < EcalRecHitCollection >	recHitsToken_

const EcalSeverityLevelAlgo *	severity_

int	severityLevel_

edm::EDGetTokenT < reco::SuperClusterCollection >	superClustersToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

Adapted from InterestingDetIdCollectionProducer by J.Bendavid

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

The following classes of "interesting id" are considered

1.All rechits included in all subclusters, plus in a region around  the seed of each subcluster
  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 InterestingDetIdFromSuperClusterProducer.h.

Constructor & Destructor Documentation

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

explicit

ctor

Definition at line 23 of file InterestingDetIdFromSuperClusterProducer.cc.

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

 {
 
   recHitsToken_ = 
       consumes<EcalRecHitCollection>(iConfig.getParameter< edm::InputTag > ("recHitsLabel"));
   superClustersToken_ = 
       consumes<reco::SuperClusterCollection>(iConfig.getParameter< edm::InputTag >("superClustersLabel"));
 
   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("InterestingDetIdFromSuperClusterProducerError") << "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 InterestingDetIdFromSuperClusterProducer::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

finaloverridevirtual

Reimplemented from edm::EDProducer.

Definition at line 47 of file InterestingDetIdFromSuperClusterProducer.cc.

References caloTopology_, edm::EventSetup::get(), edm::ESHandle< class >::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 InterestingDetIdFromSuperClusterProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

producer

Implements edm::EDProducer.

Definition at line 60 of file InterestingDetIdFromSuperClusterProducer.cc.

References funct::abs(), caloTopology_, DetId::det(), EcalBarrel, EcalRecHit::energy(), first, edm::Event::getByToken(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), cuy::ii, interestingDetIdCollection_, EcalTools::isNextToDead(), EEDetId::isNextToRingBoundary(), keepNextToBoundary_, keepNextToDead_, EcalRecHit::kTPSaturated, minimalEtaSize_, minimalPhiSize_, DetId::null(), edm::Event::put(), recHitsToken_, severity_, EcalSeverityLevelAlgo::severityLevel(), severityLevel_, python.multivaluedict::sort(), DetId::subdetId(), superClustersToken_, and ecaldqm::topology().

 {
    using namespace edm;
    using namespace std;
 
    // take BasicClusters
   Handle<reco::SuperClusterCollection> pClusters;
   iEvent.getByToken(superClustersToken_, pClusters);
   
   // take EcalRecHits
   Handle<EcalRecHitCollection> recHitsHandle;
   iEvent.getByToken(recHitsToken_,recHitsHandle);
 
   //Create empty output collections
   std::vector<DetId> indexToStore;
   indexToStore.reserve(1000);
 
   reco::SuperClusterCollection::const_iterator sclusIt;
 
   std::vector<DetId> xtalsToStore;
   xtalsToStore.reserve(50);
   
   //loop over superclusters
   for (sclusIt=pClusters->begin(); sclusIt!=pClusters->end(); sclusIt++) {
     //loop over subclusters
     for (reco::CaloCluster_iterator clusIt = sclusIt->clustersBegin(); clusIt!=sclusIt->clustersEnd(); ++clusIt) {   
       //PG barrel
       
       float eMax=0.;
       DetId eMaxId(0);
 
       std::vector<std::pair<DetId,float> > clusterDetIds = (*clusIt)->hitsAndFractions();
       std::vector<std::pair<DetId,float> >::iterator posCurrent;
 
       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_);
       std::vector<std::pair<DetId,float > > xtalsInClus=(*clusIt)->hitsAndFractions();
       
       for (unsigned int ii=0;ii<xtalsInClus.size();ii++)
         {
             xtalsToStore.push_back(xtalsInClus[ii].first);
         }
       
       indexToStore.insert(indexToStore.end(),xtalsToStore.begin(),xtalsToStore.end());
     }
   }
 
 
  
   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());
   
   std::auto_ptr< DetIdCollection > detIdCollection (new DetIdCollection(indexToStore) ) ;
 
  
   iEvent.put( detIdCollection, interestingDetIdCollection_ );
 
 }