EcalChannelKiller Class Reference

#include <RecoCaloTools/EcalChannelKiller/src/EcalChannelKiller.cc>

Inheritance diagram for EcalChannelKiller:

Public Member Functions
	EcalChannelKiller (const edm::ParameterSet &)
	~EcalChannelKiller ()
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
std::string	DeadChannelFileName_
std::string	hitCollection_
std::string	hitProducer_
std::string	reducedHitCollection_

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

Date:

2010/08/06 20:24:49

Revision:

1.3

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

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

Definition at line 30 of file EcalChannelKiller.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 48 of file EcalChannelKiller.cc.

References edm::ParameterSet::getParameter().

{

   hitProducer_          = ps.getParameter<std::string>("hitProducer");
   hitCollection_        = ps.getParameter<std::string>("hitCollection");
   reducedHitCollection_ = ps.getParameter<std::string>("reducedHitCollection");
   DeadChannelFileName_  = ps.getParameter<std::string>("DeadChannelsFile");

   produces< EcalRecHitCollection >(reducedHitCollection_);


}

EcalChannelKiller::~EcalChannelKiller

(

)

Definition at line 62 of file EcalChannelKiller.cc.

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

}

Member Function Documentation

void EcalChannelKiller::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 133 of file EcalChannelKiller.cc.

{

  //Open the DeadChannel file, read it.
  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 EcalChannelKiller::endJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 160 of file EcalChannelKiller.cc.

                          {
}

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

privatevirtual

Implements edm::EDProducer.

Definition at line 77 of file EcalChannelKiller.cc.

References edm::Event::getByLabel(), and edm::Event::put().

{
   using namespace edm;

   // get the hit collection from the event:
   edm::Handle<EcalRecHitCollection> rhcHandle;
   iEvent.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();
 
   int nRed = 0;
   
   // create an auto_ptr to a EcalRecHitCollection, copy the RecHits into it and put in the Event:
   std::auto_ptr< EcalRecHitCollection > redCollection(new EcalRecHitCollection);
   
   
   for(EcalRecHitCollection::const_iterator it = hit_collection->begin(); it != hit_collection->end(); ++it) {
     
     double NewEnergy =it->energy();
     bool ItIsDead=false;
     //Dead Cells are read from text files
     std::vector<EBDetId>::const_iterator DeadCell;
     for(DeadCell=ChannelsDeadID.begin();DeadCell<ChannelsDeadID.end();DeadCell++){
       if(it->detid()==*DeadCell){
     ItIsDead=true;
     NewEnergy =0.;
     nRed++;
     
       }
     }//End looping on vector of Dead Cells

     // Make a new RecHit
     //
     // TODO what will be the it->time() for D.C. ?
     // Could we use it for "correction" identification?
     //
     if(!ItIsDead){
       EcalRecHit NewHit(it->id(),NewEnergy,it->time());
       redCollection->push_back( NewHit );
     }
   }
   //   std::cout << "total # hits: " << nTot << "  #hits with E = " << 0 << " GeV : " << nRed << std::endl;
   
   iEvent.put(redCollection, reducedHitCollection_);
   
}

Member Data Documentation

std::vector<EBDetId> EcalChannelKiller::ChannelsDeadID

private

Definition at line 45 of file EcalChannelKiller.h.

std::string EcalChannelKiller::DeadChannelFileName_

private

Definition at line 44 of file EcalChannelKiller.h.

std::string EcalChannelKiller::hitCollection_

private

Definition at line 42 of file EcalChannelKiller.h.

std::string EcalChannelKiller::hitProducer_

private

Definition at line 41 of file EcalChannelKiller.h.

std::string EcalChannelKiller::reducedHitCollection_

private

Definition at line 43 of file EcalChannelKiller.h.