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

bool	hasAbilityToProduceInRuns () 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 35 of file ReducedRecHitCollectionProducer.h.

Constructor & Destructor Documentation

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

explicit

ctor

Definition at line 25 of file ReducedRecHitCollectionProducer.cc.

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

 {
 
   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_) ;
   
 }

ReducedRecHitCollectionProducer::~ReducedRecHitCollectionProducer ( )

override

Definition at line 47 of file ReducedRecHitCollectionProducer.cc.

48 {}

Member Function Documentation

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

override

producer

Definition at line 53 of file ReducedRecHitCollectionProducer.cc.

References popcon2dropbox::copy(), edm::SortedCollection< T, SORT >::end(), spr::find(), edm::SortedCollection< T, SORT >::find(), edm::Event::getByToken(), cuy::ii, interestingDetIdCollections_, edm::HandleBase::isValid(), tablePrinter::labels, eostools::move(), edm::Event::put(), recHitsToken_, reducedHitsCollection_, lumiQTWidget::t, and tier0::unique().

 {
    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_ );
 }