#include <InterestingDetIdCollectionProducer.h>

Inheritance diagram for InterestingDetIdCollectionProducer:

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

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

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

	~InterestingDetIdCollectionProducer ()

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
edm::InputTag	basicClusters_

const CaloTopology *	caloTopology_

std::string	interestingDetIdCollection_

int	minimalEtaSize_

int	minimalPhiSize_

edm::InputTag	recHitsLabel_

const EcalSeverityLevelAlgo *	severity_

int	severityLevel_

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

typedef WorkerT< EDProducer >	WorkerType

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)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Original author: Paolo Meridiani PH/CMG

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

Definition at line 34 of file InterestingDetIdCollectionProducer.h.

Constructor & Destructor Documentation

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

explicit

ctor

Definition at line 24 of file InterestingDetIdCollectionProducer.cc.

References basicClusters_, edm::ParameterSet::getParameter(), interestingDetIdCollection_, minimalEtaSize_, minimalPhiSize_, recHitsLabel_, and severityLevel_.

 {
 
   recHitsLabel_ = iConfig.getParameter< edm::InputTag > ("recHitsLabel");
   basicClusters_ = 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");
 }

InterestingDetIdCollectionProducer::~InterestingDetIdCollectionProducer ( )

Definition at line 44 of file InterestingDetIdCollectionProducer.cc.

45 {}

Member Function Documentation

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

virtual

Reimplemented from edm::EDProducer.

Definition at line 47 of file InterestingDetIdCollectionProducer.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 InterestingDetIdCollectionProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

producer

Implements edm::EDProducer.

Definition at line 60 of file InterestingDetIdCollectionProducer.cc.

References basicClusters_, caloTopology_, DetId::det(), jptDQMConfig_cff::eMax, CaloRecHit::energy(), spr::find(), edm::Event::getByLabel(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), interestingDetIdCollection_, EcalRecHit::kTPSaturated, minimalEtaSize_, minimalPhiSize_, DetId::null(), edm::Event::put(), recHitsLabel_, severity_, EcalSeverityLevelAlgo::severityLevel(), severityLevel_, and DetId::subdetId().

 {
    using namespace edm;
    using namespace std;
 
    // take BasicClusters
   Handle<reco::BasicClusterCollection> pClusters;
   iEvent.getByLabel(basicClusters_, pClusters);
   
   // take EcalRecHits
   Handle<EcalRecHitCollection> recHitsHandle;
   iEvent.getByLabel(recHitsLabel_,recHitsHandle);
 
   //Create empty output collections
   std::auto_ptr< DetIdCollection > detIdCollection (new DetIdCollection() ) ;
   
   reco::BasicClusterCollection::const_iterator clusIt;
   
   for (clusIt=pClusters->begin(); clusIt!=pClusters->end(); 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());
     std::vector<DetId> xtalsToStore=topology->getWindow(eMaxId,minimalEtaSize_,minimalPhiSize_);
     std::vector<std::pair<DetId,float > > xtalsInClus=(*clusIt).hitsAndFractions();
     
     for (unsigned int ii=0;ii<xtalsInClus.size();ii++)
       {
     if (std::find(xtalsToStore.begin(),xtalsToStore.end(),xtalsInClus[ii].first) == xtalsToStore.end())
       xtalsToStore.push_back(xtalsInClus[ii].first);
       }
     
     for (unsigned int iCry=0;iCry<xtalsToStore.size();iCry++)
       {
     if ( 
         std::find(detIdCollection->begin(),detIdCollection->end(),xtalsToStore[iCry]) == detIdCollection->end()
         )
       detIdCollection->push_back(xtalsToStore[iCry]);
       }     
   }
 
   // also add recHits of dead TT if the corresponding TP is saturated
   for (EcalRecHitCollection::const_iterator it = recHitsHandle->begin(); it != recHitsHandle->end(); ++it) {
           if ( it->checkFlag(EcalRecHit::kTPSaturated) ) {
                   if ( std::find( detIdCollection->begin(), detIdCollection->end(), it->id() ) == detIdCollection->end()
                      ) {
                           detIdCollection->push_back( it->id() );
                   }
           }else if ( severityLevel_>=0 && severity_->severityLevel(*it) >=severityLevel_){
         detIdCollection->push_back( it->id() );
       }
 
   }
   
   //  std::cout << "Interesting DetId Collection size is " << detIdCollection->size() <<  " BCs are " << pClusters->size() << std::endl;
   iEvent.put( detIdCollection, interestingDetIdCollection_ );
 
 }