#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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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_);
 }