#include <ReducedESRecHitCollectionProducer.h>

Inheritance diagram for ReducedESRecHitCollectionProducer:

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

void	collectIds (const ESDetId strip1, const ESDetId strip2, const int &row=0)

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

	ReducedESRecHitCollectionProducer (const edm::ParameterSet &pset)

virtual	~ReducedESRecHitCollectionProducer ()

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

	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
std::vector< ConsumesInfo >	consumesInfo () const

	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
std::set< DetId >	collectedIds_

const EcalPreshowerGeometry *	geometry_p

edm::EDGetTokenT < ESRecHitCollection >	InputRecHitES_

edm::EDGetTokenT < reco::SuperClusterCollection >	InputSuperClusterEE_

std::vector< edm::EDGetTokenT < DetIdCollection > >	interestingDetIdCollections_

std::vector< edm::EDGetTokenT < DetIdCollection > >	interestingDetIdCollectionsNotToClean_

std::string	OutputLabelES_

double	scEtThresh_

CaloSubdetectorTopology *	topology_p

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 24 of file ReducedESRecHitCollectionProducer.h.

Constructor & Destructor Documentation

ReducedESRecHitCollectionProducer::ReducedESRecHitCollectionProducer ( const edm::ParameterSet & pset )

Definition at line 20 of file ReducedESRecHitCollectionProducer.cc.

References edm::ParameterSet::getParameter(), InputRecHitES_, InputSuperClusterEE_, interestingDetIdCollections_, OutputLabelES_, scEtThresh_, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and edm::vector_transform().

                                                                                              :
   geometry_p(0),
   topology_p(0)
 {
 
  scEtThresh_          = ps.getParameter<double>("scEtThreshold");
 
  InputRecHitES_       = 
      consumes<ESRecHitCollection>(ps.getParameter<edm::InputTag>("EcalRecHitCollectionES"));
  InputSuperClusterEE_ = 
      consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("EndcapSuperClusterCollection")); 
 
  OutputLabelES_       = ps.getParameter<std::string>("OutputLabel_ES");
  
  interestingDetIdCollections_  = 
      edm::vector_transform(
            ps.getParameter<std::vector<edm::InputTag>>("interestingDetIds"),
            [this](edm::InputTag const & tag) { 
                return consumes<DetIdCollection>(tag); 
            }
            );
 
  interestingDetIdCollectionsNotToClean_ = edm::vector_transform(ps.getParameter<std::vector<edm::InputTag>>("interestingDetIdsNotToClean"),
                                 [this](edm::InputTag const & tag) 
                                 { return consumes<DetIdCollection>(tag); }
                                 );
 
  produces< EcalRecHitCollection > (OutputLabelES_);
  
 }

ReducedESRecHitCollectionProducer::~ReducedESRecHitCollectionProducer ( )

virtual

Definition at line 51 of file ReducedESRecHitCollectionProducer.cc.

References topology_p.

                                                                       {
   if (topology_p) delete topology_p;
 }

Member Function Documentation

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

finaloverridevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 55 of file ReducedESRecHitCollectionProducer.cc.

References DetId::Ecal, EcalPreshower, geometry, geometry_p, edm::EventSetup::get(), and topology_p.

                                                                                           {
   ESHandle<CaloGeometry> geoHandle;
   iSetup.get<CaloGeometryRecord>().get(geoHandle);
   const CaloSubdetectorGeometry *geometry = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
   geometry_p = dynamic_cast<const EcalPreshowerGeometry *>(geometry);
   if (!geometry_p){
     edm::LogError("WrongGeometry")<<
       "could not cast the subdet geometry to preshower geometry";
   }
   
   if (geometry) topology_p = new EcalPreshowerTopology(geoHandle);
   
 }

void ReducedESRecHitCollectionProducer::collectIds	(	const ESDetId	strip1,
		const ESDetId	strip2,
		const int &	row = `0`
	)

Definition at line 169 of file ReducedESRecHitCollectionProducer.cc.

References collectedIds_, CaloNavigator< T, TOPO >::east(), CaloNavigator< T, TOPO >::home(), i, GetRecoTauVFromDQM_MC_cff::next, CaloNavigator< T, TOPO >::north(), CaloNavigator< T, TOPO >::setHome(), CaloNavigator< T, TOPO >::south(), topology_p, and CaloNavigator< T, TOPO >::west().

Referenced by produce().

                                                                                                                   {
 
   //cout<<row<<endl;
   
   map<DetId,const EcalRecHit*>::iterator it;
   map<DetId, int>::iterator itu;
   ESDetId next;
   ESDetId strip1;
   ESDetId strip2;
 
   strip1 = esDetId1;
   strip2 = esDetId2;
 
   EcalPreshowerNavigator theESNav1(strip1, topology_p);
   theESNav1.setHome(strip1);
   
   EcalPreshowerNavigator theESNav2(strip2, topology_p);
   theESNav2.setHome(strip2);
 
   if (row == 1) {
     if (strip1 != ESDetId(0)) strip1 = theESNav1.north();
     if (strip2 != ESDetId(0)) strip2 = theESNav2.east();
   } else if (row == -1) {
     if (strip1 != ESDetId(0)) strip1 = theESNav1.south();
     if (strip2 != ESDetId(0)) strip2 = theESNav2.west();
   }
 
   // Plane 1 
   if (strip1 == ESDetId(0)) {
   } else {
     collectedIds_.insert(strip1);
     //cout<<"center : "<<strip1<<endl;
     // east road 
     for (int i=0; i<15; ++i) {
       next = theESNav1.east();
       //cout<<"east : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       }
     }
 
     // west road 
     theESNav1.setHome(strip1);
     theESNav1.home();
     for (int i=0; i<15; ++i) {
       next = theESNav1.west();
       //cout<<"west : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       }
     }
 
   }
 
   if (strip2 == ESDetId(0)) {
   } else {
     collectedIds_.insert(strip2);
     //cout<<"center : "<<strip2<<endl;
     // north road 
     for (int i=0; i<15; ++i) {
       next = theESNav2.north();
       //cout<<"north : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       } 
     }
 
     // south road 
     theESNav2.setHome(strip2);
     theESNav2.home();
     for (int i=0; i<15; ++i) {
       next = theESNav2.south();
       //cout<<"south : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       }
     }
   }
 }

void ReducedESRecHitCollectionProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

Implements edm::stream::EDProducerBase.

Definition at line 69 of file ReducedESRecHitCollectionProducer.cc.

References collectedIds_, collectIds(), geometry_p, edm::Event::getByToken(), EcalPreshowerGeometry::getClosestCellInPlane(), hit::id, InputRecHitES_, InputSuperClusterEE_, interestingDetIdCollections_, interestingDetIdCollectionsNotToClean_, edm::HandleBase::isValid(), HLT_25ns10e33_v2_cff::labels, edm::EDConsumerBase::labelsForToken(), edm::ProductLabels::module, convertSQLitetoXML_cfg::output, OutputLabelES_, point, edm::ProductLabels::process, edm::Handle< T >::product(), edm::ProductLabels::productInstance, edm::Event::put(), scEtThresh_, lumiQTWidget::t, and unpackBuffers-CaloStage2::token.

                                                                                         {
 
 
   edm::Handle<ESRecHitCollection> ESRecHits_;
   e.getByToken(InputRecHitES_, ESRecHits_);
   
   std::auto_ptr<EcalRecHitCollection> output(new EcalRecHitCollection);
 
   edm::Handle<reco::SuperClusterCollection> pEndcapSuperClusters;
   e.getByToken(InputSuperClusterEE_, pEndcapSuperClusters);
   {
     const reco::SuperClusterCollection* eeSuperClusters = pEndcapSuperClusters.product();
     
     for (reco::SuperClusterCollection::const_iterator isc = eeSuperClusters->begin(); isc != eeSuperClusters->end(); ++isc) {
 
       if (isc->energy() < scEtThresh_) continue;
       if (fabs(isc->eta()) < 1.65 || fabs(isc->eta()) > 2.6) continue;
       //cout<<"SC energy : "<<isc->energy()<<" "<<isc->eta()<<endl;
 
       //Int_t nBC = 0;
       reco::CaloCluster_iterator ibc = isc->clustersBegin();
       for ( ; ibc != isc->clustersEnd(); ++ibc ) {  
 
     //cout<<"BC : "<<nBC<<endl;
 
     const GlobalPoint point((*ibc)->x(),(*ibc)->y(),(*ibc)->z());
 
     ESDetId esId1 = geometry_p->getClosestCellInPlane(point, 1);
     ESDetId esId2 = geometry_p->getClosestCellInPlane(point, 2);
     
     collectIds(esId1, esId2, 0);
     collectIds(esId1, esId2, 1);
     collectIds(esId1, esId2, -1);
 
     //nBC++;
       }
       
     }
     
   }
 
 
   edm::Handle<DetIdCollection > detId;
   for( unsigned int t = 0; t < interestingDetIdCollections_.size(); ++t )
     {
       e.getByToken(interestingDetIdCollections_[t],detId);    
       if(!detId.isValid())
       {
           Labels labels;
           labelsForToken(interestingDetIdCollections_[t], labels);
           edm::LogError("MissingInput")<<"no reason to skip detid from : (" << labels.module << ", "
                                                                             << labels.productInstance << ", "
                                                                             << labels.process << ")" << std::endl;
           continue;
       }
       collectedIds_.insert(detId->begin(),detId->end());
     }
 
 
   //screw it, cant think of a better solution, not the best but lets run over all the rec hits, remove the ones failing cleaning
   //and then merge in the collection not to be cleaned
   //mainly as I suspect its more efficient to find an object in the DetIdSet rather than the rec-hit in the rec-hit collecition
   //with only a det id
   //if its a CPU issues then revisit
   for(const auto& hit : *ESRecHits_) {
     if(hit.recoFlag()==1 || hit.recoFlag()==14 || (hit.recoFlag()<=10 && hit.recoFlag()>=5)){ //right we might need to erase it from the collection
       auto idIt = collectedIds_.find(hit.id());
       if(idIt!=collectedIds_.end()) collectedIds_.erase(idIt);
     }
   }
    
   
   for(const auto& token : interestingDetIdCollectionsNotToClean_) {
     e.getByToken(token,detId);    
     if(!detId.isValid()){ //meh might as well keep the warning
       Labels labels;
       labelsForToken(token, labels);
       edm::LogError("MissingInput")<<"no reason to skip detid from : (" << labels.module << ", "
                    << labels.productInstance << ", "
                    << labels.process << ")" << std::endl;
       continue;
     }
     collectedIds_.insert(detId->begin(),detId->end());
   }
   
 
   output->reserve( collectedIds_.size());
   EcalRecHitCollection::const_iterator it;
   for (it = ESRecHits_->begin(); it != ESRecHits_->end(); ++it) {
     if (collectedIds_.find(it->id())!=collectedIds_.end()){
       output->push_back(*it);
     }
   }
   collectedIds_.clear();
 
   e.put(output, OutputLabelES_);
 
 }