#include <EventFilter/EcalRawToRecHitProducer/src/EcalRawToRecHitProducer.cc>

Inheritance diagram for EcalRawToRecHitProducer:

Public Types
typedef edm::LazyGetter < EcalRecHit >	EcalRecHitLazyGetter

typedef edm::RefGetter < EcalRecHit >	EcalRecHitRefGetter

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

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

Public Member Functions
	EcalRawToRecHitProducer (const edm::ParameterSet &)

	~EcalRawToRecHitProducer ()

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 Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
EcalCleaningAlgo *	cleaningAlgo_

std::string	EBrechitCollection_

std::string	EErechitCollection_

edm::InputTag	lsourceTag_

std::string	rechitCollection_

edm::InputTag	sourceTag_

bool	splitOutput_

Additional Inherited Members
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (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

Description: <one line="" class="" summary>="">

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

Definition at line 55 of file EcalRawToRecHitProducer.h.

Member Typedef Documentation

typedef edm::LazyGetter<EcalRecHit> EcalRawToRecHitProducer::EcalRecHitLazyGetter

Definition at line 58 of file EcalRawToRecHitProducer.h.

typedef edm::RefGetter<EcalRecHit> EcalRawToRecHitProducer::EcalRecHitRefGetter

Definition at line 59 of file EcalRawToRecHitProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 4 of file EcalRawToRecHitProducer.cc.

References cleaningAlgo_, EBrechitCollection_, EErechitCollection_, edm::ParameterSet::empty(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), LogDebug, lsourceTag_, rechitCollection_, sourceTag_, and splitOutput_.

 {
   lsourceTag_=iConfig.getParameter<edm::InputTag>("lazyGetterTag");
   sourceTag_=iConfig.getParameter<edm::InputTag>("sourceTag");
 
   splitOutput_=iConfig.getParameter<bool>("splitOutput");
   if (splitOutput_){
     EBrechitCollection_=iConfig.getParameter<std::string>("EBrechitCollection");
     EErechitCollection_=iConfig.getParameter<std::string>("EErechitCollection");
     produces<EBRecHitCollection>(EBrechitCollection_);
     produces<EERecHitCollection>(EErechitCollection_);
     LogDebug("EcalRawToRecHit|Producer")<<"ready to create rechits from lazy getter: "<<lsourceTag_
                     <<"\n using region ref from: "<<sourceTag_
                     <<"\n splitting in two collections"
                     <<"\n EB instance: "<<EBrechitCollection_
                     <<"\n EE instance: "<<EErechitCollection_;
   }
   else{
     rechitCollection_=iConfig.getParameter<std::string>("rechitCollection");
     produces<EcalRecHitCollection>(rechitCollection_);
     LogDebug("EcalRawToRecHit|Producer")<<"ready to create rechits from lazy getter: "<<lsourceTag_
                     <<"\n using region ref from: "<<sourceTag_
                     <<"\n not splitting the output collection.";
   }
 
   cleaningAlgo_=0;
   if (iConfig.exists("cleaningConfig")){
     const edm::ParameterSet & cleaning=iConfig.getParameter<edm::ParameterSet>("cleaningConfig");
     if (!cleaning.empty())
       cleaningAlgo_ = new EcalCleaningAlgo(cleaning);
   }
 }

EcalRawToRecHitProducer::~EcalRawToRecHitProducer ( )

Definition at line 38 of file EcalRawToRecHitProducer.cc.

References cleaningAlgo_.

 {
   if (cleaningAlgo_) delete cleaningAlgo_;
 }

Member Function Documentation

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

privatevirtual

Implements edm::EDProducer.

Definition at line 50 of file EcalRawToRecHitProducer.cc.

References cleaningAlgo_, cond::rpcobgas::detid, EBrechitCollection_, EErechitCollection_, edm::Event::getByLabel(), MyWatcher::lap(), LogDebug, lsourceTag_, edm::Event::put(), rechitCollection_, EcalCleaningAlgo::setFlags(), sourceTag_, splitOutput_, and DetId::subdetId().

Referenced by python.JSONExport.JsonExport::export(), and python.HTMLExport.HTMLExport::export().

 {
   using namespace edm;
 
   MyWatcher watcher("Producer");
   LogDebug("EcalRawToRecHit|Producer")<<watcher.lap();
 
   //retrieve a lazygetter
   edm::Handle<EcalRecHitLazyGetter> lgetter;
   iEvent.getByLabel(lsourceTag_, lgetter);
   LogDebug("EcalRawToRecHit|Producer")<<"lazy getter retreived from: "<<lsourceTag_<<(lgetter.failedToGet()?" not valid ":"valid")
                       <<watcher.lap();
   
   //retrieve a refgetter
   edm::Handle<EcalRecHitRefGetter> rgetter;
   iEvent.getByLabel(sourceTag_ ,rgetter);
   LogDebug("EcalRawToRecHit|Producer")<<"ref getter retreived from: "<<sourceTag_<<(rgetter.failedToGet()?" not valid ":"valid")
                       <<watcher.lap();
 
  
   if (splitOutput_){
     //prepare the output collection
     std::auto_ptr<EBRecHitCollection> EBrechits( new EBRecHitCollection );
     std::auto_ptr<EERecHitCollection> EErechits( new EERecHitCollection );
     //loop the refgetter
     unsigned int iR=0;
     EcalRecHitRefGetter::const_iterator iRegion=rgetter->begin();
     EcalRecHitRefGetter::const_iterator iRegionEnd=rgetter->end();
     for (;iRegion!=iRegionEnd;++iRegion){
       LogDebug("EcalRawToRecHit|Producer")<<"looping over refgetter region: "<<iR<<watcher.lap();
       std::vector<EcalRecHit>::const_iterator iRecHit=lgetter->begin_record()+iRegion->start();
       std::vector<EcalRecHit>::const_iterator iRecHitEnd =lgetter->begin_record()+iRegion->finish();
       for (;iRecHit!=iRecHitEnd;iRecHit++){
     DetId detid =iRecHit->id();
     //split barrel and endcap
     int EcalNum=detid.subdetId(); //1 stands for Barrel, 2 for endcaps
     LogDebug("EcalRawToRecHit|Producer")<<"subdetId is: "<<EcalNum;
     if (EcalNum==1) EBrechits->push_back(*iRecHit);
     else if (EcalNum==2) EErechits->push_back(*iRecHit);
     else {
       edm::LogError("IncorrectRecHit")<<" a subdetid is not recognized. recHit on :"<< iRecHit->id().rawId() 
                            <<" is lost.";
     }//subdetid
       }//loop over things in region
       LogDebug("EcalRawToRecHit|Producer")<<"looping over refgetter region: "<<iR++<<" done"
                       <<watcher.lap();
     }//loop over regions
 
     LogDebug("EcalRawToRecHit|Producer")<<EBrechits->size()<<" EB recHits to be put with instance: "<<EBrechitCollection_
                     <<"\n"<<EErechits->size()<<" EE recHits to be put with instance: "<<EErechitCollection_
                     << watcher.lap();
     
     // cleaning of anomalous signals, aka spikes
     // only doable once we have a "global" collection of hits
     if (cleaningAlgo_){
       EBrechits->sort();
       EErechits->sort();
       cleaningAlgo_->setFlags(*EBrechits);
       cleaningAlgo_->setFlags(*EErechits);
     }
 
     iEvent.put(EBrechits, EBrechitCollection_);
     iEvent.put(EErechits, EErechitCollection_);
     LogDebug("EcalRawToRecHit|Producer")<<"collections uploaded."
                     << watcher.lap();
   }
   else{
     //prepare the output collection
     std::auto_ptr< EcalRecHitCollection > rechits( new EcalRecHitCollection);
     //loop the refgetter
     unsigned int iR=0;
     EcalRecHitRefGetter::const_iterator iRegion=rgetter->begin();
     EcalRecHitRefGetter::const_iterator iRegionEnd=rgetter->end();
     for (;iRegion!=iRegionEnd;++iRegion){
       LogDebug("EcalRawToRecHit|Producer")<<"looping over refgetter region: "<<iR<<watcher.lap();
       std::vector<EcalRecHit>::const_iterator iRecHit=lgetter->begin_record()+iRegion->start();
       std::vector<EcalRecHit>::const_iterator iRecHitEnd=lgetter->begin_record()+iRegion->finish();
       for (;iRecHit!=iRecHitEnd;iRecHit++){
     LogDebug("EcalRawToRecHit|Producer")<<"dereferencing rechit ref.";
     DetId detid =iRecHit->id();
     int EcalNum=detid.subdetId(); //1 stands for Barrel, 2 for endcaps
     LogDebug("EcalRawToRecHit|Producer")<<"subdetId is: "<<EcalNum;
     rechits->push_back(*iRecHit);
       }//loop over things in region
       LogDebug("EcalRawToRecHit|Producer")<<"looping over refgetter region: "<<iR++<<" done"<<watcher.lap();
     }//loop over regions
 
     if (cleaningAlgo_){
       rechits->sort();
       cleaningAlgo_->setFlags(*rechits);
     }
     LogDebug("EcalRawToRecHit|Producer")<<rechits->size()<<" rechits to be put."<< watcher.lap();
     iEvent.put(rechits,rechitCollection_);
     LogDebug("EcalRawToRecHit|Producer")<<"collections uploaded."
                     << watcher.lap();
   }
 
 }