#include <Calibration/EcalAlCaRecoProducers/src/ReducedRecHitCollectionProducer.cc>

Inheritance diagram for ReducedRecHitCollectionProducer:

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

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

	~ReducedRecHitCollectionProducer () override

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

edm::EDGetTokenT< EcalRecHitCollection >	recHitsToken_

std::string	reducedHitsCollection_

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

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

Definition at line 33 of file ReducedRecHitCollectionProducer.h.

Constructor & Destructor Documentation

◆ ReducedRecHitCollectionProducer()

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

explicit

ctor

Definition at line 21 of file ReducedRecHitCollectionProducer.cc.

                                                                                                {
   recHitsToken_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("recHitsLabel"));
  
   interestingDetIdCollections_ =
       edm::vector_transform(iConfig.getParameter<std::vector<edm::InputTag>>("interestingDetIdCollections"),
                             [this](edm::InputTag const& tag) { return consumes<DetIdCollection>(tag); });
  
   reducedHitsCollection_ = iConfig.getParameter<std::string>("reducedHitsCollection");
  
   //register your products
   produces<EcalRecHitCollection>(reducedHitsCollection_);
 }

References edm::ParameterSet::getParameter(), interestingDetIdCollections_, recHitsToken_, reducedHitsCollection_, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and edm::vector_transform().

◆ ~ReducedRecHitCollectionProducer()

ReducedRecHitCollectionProducer::~ReducedRecHitCollectionProducer ( )

override

Definition at line 34 of file ReducedRecHitCollectionProducer.cc.

34 {}

Member Function Documentation

◆ produce()

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

override

producer

Definition at line 37 of file ReducedRecHitCollectionProducer.cc.

                                                                                            {
   using namespace edm;
   using namespace std;
  
   if (interestingDetIdCollections_.empty()) {
     edm::LogError("ReducedRecHitCollectionProducer") << "VInputTag collections empty";
     return;
   }
  
   Handle<DetIdCollection> detIds;
   iEvent.getByToken(interestingDetIdCollections_[0], detIds);
   std::vector<DetId> xtalsToStore((*detIds).size());
   std::copy((*detIds).begin(), (*detIds).end(), xtalsToStore.begin());
  
   //Merging DetIds from different collections
   for (unsigned int t = 1; t < interestingDetIdCollections_.size(); ++t) {
     Handle<DetIdCollection> detId;
     iEvent.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;
     }
  
     for (unsigned int ii = 0; ii < (*detId).size(); ii++) {
       if (std::find(xtalsToStore.begin(), xtalsToStore.end(), (*detId)[ii]) == xtalsToStore.end())
         xtalsToStore.push_back((*detId)[ii]);
     }
   }
  
   Handle<EcalRecHitCollection> recHitsHandle;
   iEvent.getByToken(recHitsToken_, recHitsHandle);
   if (!recHitsHandle.isValid()) {
     edm::LogError("ReducedRecHitCollectionProducer") << "RecHit collection not found";
     return;
   }
  
   //Create empty output collections
   auto miniRecHitCollection = std::make_unique<EcalRecHitCollection>();
  
   for (unsigned int iCry = 0; iCry < xtalsToStore.size(); iCry++) {
     EcalRecHitCollection::const_iterator iRecHit = recHitsHandle->find(xtalsToStore[iCry]);
     if ((iRecHit != recHitsHandle->end()) &&
         (miniRecHitCollection->find(xtalsToStore[iCry]) == miniRecHitCollection->end()))
       miniRecHitCollection->push_back(*iRecHit);
   }
  
   std::sort(xtalsToStore.begin(), xtalsToStore.end());
   std::unique(xtalsToStore.begin(), xtalsToStore.end());
  
   //   std::cout << "New Collection " << reducedHitsCollection_ << " size is " << miniRecHitCollection->size() << " original is " << recHitsHandle->size() << std::endl;
   iEvent.put(std::move(miniRecHitCollection), reducedHitsCollection_);
 }