#include <RecoLocalCalo/EcalDeadChannelRecoveryProducers/src/EcalDeadChannelRecoveryProducers.cc>

Inheritance diagram for EcalDeadChannelRecoveryProducers:

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

	~EcalDeadChannelRecoveryProducers ()

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	beginJob ()

virtual void	endJob ()

virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
std::vector< EBDetId >	ChannelsDeadID

bool	CorrectDeadCells_

std::string	CorrectionMethod_

std::string	DeadChannelFileName_

std::string	hitCollection_

std::string	hitProducer_

std::string	reducedHitCollection_

double	Sum8GeVThreshold_

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)

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

Date:: 2009/12/14 22:23:53

Revision:: 1.2

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

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

Definition at line 30 of file EcalDeadChannelRecoveryProducers.h.

Constructor & Destructor Documentation

EcalDeadChannelRecoveryProducers::EcalDeadChannelRecoveryProducers ( const edm::ParameterSet & ps )

explicit

Definition at line 57 of file EcalDeadChannelRecoveryProducers.cc.

References edm::ParameterSet::getParameter().

 {
 
    //now do what ever other initialization is needed
   CorrectDeadCells_     = ps.getParameter<bool>("CorrectDeadCells");
   CorrectionMethod_     = ps.getParameter<std::string>("CorrectionMethod");
   hitProducer_          = ps.getParameter<std::string>("hitProducer");
   hitCollection_        = ps.getParameter<std::string>("hitCollection");
   reducedHitCollection_ = ps.getParameter<std::string>("reducedHitCollection");
   DeadChannelFileName_ = ps.getParameter<std::string>("DeadChannelsFile");
   Sum8GeVThreshold_= ps.getParameter<double>("Sum8GeVThreshold");
 
    produces< EcalRecHitCollection >(reducedHitCollection_);
 }

EcalDeadChannelRecoveryProducers::~EcalDeadChannelRecoveryProducers ( )

Definition at line 73 of file EcalDeadChannelRecoveryProducers.cc.

 {
  
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void EcalDeadChannelRecoveryProducers::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 143 of file EcalDeadChannelRecoveryProducers.cc.

References gather_cfg::cout.

 {
     FILE* DeadCha;
     printf("Dead Channels FILE: %s\n",DeadChannelFileName_.c_str());
     DeadCha = fopen(DeadChannelFileName_.c_str(),"r");
 
     int fileStatus=0;
     int ieta=-10000;
     int iphi=-10000;
     while(fileStatus != EOF) {
     fileStatus = fscanf(DeadCha,"%d %d\n",&ieta,&iphi);
     //    std::cout<<" ieta "<<ieta<<" iphi "<<iphi<<std::endl;
     if(ieta==-10000||iphi==-10000){std::cout << "Problem reading Dead Channels file "<<std::endl;break;}
     EBDetId cell(ieta,iphi);
     ChannelsDeadID.push_back(cell);
     } //end while       
     fclose(DeadCha);
     //  std::cout<<" Read "<<ChannelsDeadID.size()<<" dead channels "<<std::endl;
   
 }

void EcalDeadChannelRecoveryProducers::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 166 of file EcalDeadChannelRecoveryProducers.cc.

166 {

167 }

void EcalDeadChannelRecoveryProducers::produce	(	edm::Event &	evt,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 88 of file EcalDeadChannelRecoveryProducers.cc.

References gather_cfg::cout, edm::EventSetup::get(), edm::Event::getByLabel(), and edm::Event::put().

 {
   using namespace edm;
   
   edm::ESHandle<CaloTopology> theCaloTopology;
   iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
   
   // get the hit collection from the event:
   edm::Handle<EcalRecHitCollection> rhcHandle;
   evt.getByLabel(hitProducer_, hitCollection_, rhcHandle);
   if (!(rhcHandle.isValid())) 
     {
       std::cout << "could not get a handle on the EcalRecHitCollection!" << std::endl;
       return;
     }
   const EcalRecHitCollection* hit_collection = rhcHandle.product();
   
   // create an auto_ptr to a EcalRecHitCollection, copy the RecHits into it and put it in the Event:
   std::auto_ptr< EcalRecHitCollection > redCollection(new EcalRecHitCollection);
   
   EcalDeadChannelRecoveryAlgos *DeadChannelCorrector = new EcalDeadChannelRecoveryAlgos(theCaloTopology.product());
   
   //Dead Cells are read from a text file
   std::vector<EBDetId>::const_iterator DeadCell;
 
   //
   //This should work only if we REMOVE the DC RecHit from the reduced RecHit collection 
   //
   for(EcalRecHitCollection::const_iterator it = hit_collection->begin(); it != hit_collection->end(); ++it) {     
     std::vector<EBDetId>::const_iterator CheckDead = ChannelsDeadID.begin();
     bool OverADeadRecHit=false;
     while(CheckDead<ChannelsDeadID.end()){
       if(it->detid()==*CheckDead){OverADeadRecHit=true;break;}
       CheckDead++;
     }
     if(!OverADeadRecHit)redCollection->push_back( *it );
   }
   for(DeadCell=ChannelsDeadID.begin();DeadCell<ChannelsDeadID.end();DeadCell++){
     EcalRecHit NewRecHit = DeadChannelCorrector->Correct(*DeadCell,hit_collection,CorrectionMethod_,Sum8GeVThreshold_);
     redCollection->push_back( NewRecHit );
   }
 
 
 
   
   delete DeadChannelCorrector ;
   
   evt.put(redCollection, reducedHitCollection_);
   
 }