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

Inheritance diagram for ReducedRecHitCollectionProducer:

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

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

	~ReducedRecHitCollectionProducer ()

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
std::vector< edm::InputTag >	interestingDetIdCollections_

edm::InputTag	recHitsLabel_

std::string	reducedHitsCollection_

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)

static void	prevalidate (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 ReducedRecHitCollectionProducer.h.

Constructor & Destructor Documentation

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

explicit

ctor

Definition at line 23 of file ReducedRecHitCollectionProducer.cc.

References edm::ParameterSet::getParameter(), interestingDetIdCollections_, recHitsLabel_, and reducedHitsCollection_.

 {
 
   recHitsLabel_ = iConfig.getParameter< edm::InputTag > ("recHitsLabel");
   interestingDetIdCollections_ = iConfig.getParameter< std::vector<edm::InputTag> > ("interestingDetIdCollections");
 
   reducedHitsCollection_ = iConfig.getParameter<std::string>("reducedHitsCollection");
   
    //register your products
   produces< EcalRecHitCollection > (reducedHitsCollection_) ;
   
 }

ReducedRecHitCollectionProducer::~ReducedRecHitCollectionProducer ( )

Definition at line 37 of file ReducedRecHitCollectionProducer.cc.

38 {}

Member Function Documentation

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

virtual

producer

Implements edm::EDProducer.

Definition at line 43 of file ReducedRecHitCollectionProducer.cc.

References filterCSVwithJSON::copy, spr::find(), edm::Event::getByLabel(), interestingDetIdCollections_, edm::HandleBase::isValid(), edm::Event::put(), recHitsLabel_, reducedHitsCollection_, python.multivaluedict::sort(), and lumiQTWidget::t.

 {
    using namespace edm;
    using namespace std;
 
    if (interestingDetIdCollections_.size() < 1)
      {
        edm::LogError("ReducedRecHitCollectionProducer") << "VInputTag collections empty" ;
        return;
      }
    
 
    Handle< DetIdCollection > detIds;
    iEvent.getByLabel(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.getByLabel(interestingDetIdCollections_[t],detId);
        if( !detId.isValid() ){
      edm::LogError("MissingInput")<<"no reason to skip detid from :"<<interestingDetIdCollections_[t];
      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.getByLabel(recHitsLabel_,recHitsHandle);
    if( !recHitsHandle.isValid() ) 
      {
        edm::LogError("ReducedRecHitCollectionProducer") << "RecHit collection not found";
        return;
      }
    
    //Create empty output collections
    std::auto_ptr< EcalRecHitCollection > miniRecHitCollection (new 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( miniRecHitCollection,reducedHitsCollection_ );
 }