#include <RecoHI/HiSpikeCleaner/src/HiSpikeCleaner.cc>

Inheritance diagram for HiSpikeCleaner:

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

	~HiSpikeCleaner ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob () override

virtual void	endJob () override

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

Private Attributes
double	etCut_

std::string	outputCollection_

edm::EDGetTokenT < EcalRecHitCollection >	rHInputProducerBToken_

edm::EDGetTokenT < EcalRecHitCollection >	rHInputProducerEToken_

edm::EDGetTokenT < reco::SuperClusterCollection >	sCInputProducerToken_

double	swissCutThr_

double	TimingCut_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 49 of file HiSpikeCleaner.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 71 of file HiSpikeCleaner.cc.

References etCut_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), outputCollection_, rHInputProducerBToken_, rHInputProducerEToken_, sCInputProducerToken_, AlCaHLTBitMon_QueryRunRegistry::string, swissCutThr_, and TimingCut_.

 {
    //register your products
 /* Examples
    produces<ExampleData2>();
 
    //if do put with a label
    produces<ExampleData2>("label");
 */
    //now do what ever other initialization is needed
    
   rHInputProducerBToken_  = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("recHitProducerBarrel"));
   rHInputProducerEToken_  = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("recHitProducerEndcap"));
   
   sCInputProducerToken_  = consumes<reco::SuperClusterCollection>(iConfig.getParameter<edm::InputTag>("originalSuperClusterProducer"));
   TimingCut_      = iConfig.getUntrackedParameter<double>  ("TimingCut",4.0);
   swissCutThr_      = iConfig.getUntrackedParameter<double>("swissCutThr",0.95);
   etCut_            = iConfig.getParameter<double>("etCut");
   
   outputCollection_ = iConfig.getParameter<std::string>("outputColl");
   produces<reco::SuperClusterCollection>(outputCollection_);  
 }

HiSpikeCleaner::~HiSpikeCleaner ( )

Definition at line 95 of file HiSpikeCleaner.cc.

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

Member Function Documentation

void HiSpikeCleaner::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 215 of file HiSpikeCleaner.cc.

216 {

217 }

void HiSpikeCleaner::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 221 of file HiSpikeCleaner.cc.

221 {

222 }

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 108 of file HiSpikeCleaner.cc.

References etCut_, edm::false, edm::EventSetup::get(), edm::Event::getByToken(), outputCollection_, edm::Event::put(), HI_PhotonSkim_cff::rechits, rHInputProducerBToken_, rHInputProducerEToken_, sCInputProducerToken_, fileCollector::seed, EcalTools::swissCross(), swissCutThr_, and TimingCut_.

 {
    using namespace edm;
 
    // Get raw SuperClusters from the event    
    Handle<reco::SuperClusterCollection> pRawSuperClusters;
    try { 
       iEvent.getByToken(sCInputProducerToken_, pRawSuperClusters);
    } catch ( cms::Exception& ex ) {
      edm::LogError("EgammaSCCorrectionMakerError") 
        << "Error! can't get the rawSuperClusters ";
    }    
    
    // Get the RecHits from the event
    Handle<EcalRecHitCollection> pRecHitsB;
    try { 
       iEvent.getByToken(rHInputProducerBToken_, pRecHitsB);
    } catch ( cms::Exception& ex ) {
      edm::LogError("EgammaSCCorrectionMakerError") 
        << "Error! can't get the RecHits ";
    }    
 
    // Get the RecHits from the event                                                                                                            
    Handle<EcalRecHitCollection> pRecHitsE;
    try {
      iEvent.getByToken(rHInputProducerEToken_, pRecHitsE);
    } catch ( cms::Exception& ex ) {
      edm::LogError("EgammaSCCorrectionMakerError")
        << "Error! can't get the RecHits ";
    }
 
    
    // get the channel status from the DB                                                                                                     
    //   edm::ESHandle<EcalChannelStatus> chStatus;
    //   iSetup.get<EcalChannelStatusRcd>().get(chStatus);
    
    edm::ESHandle<EcalSeverityLevelAlgo> ecalSevLvlAlgoHndl;
    iSetup.get<EcalSeverityLevelAlgoRcd>().get(ecalSevLvlAlgoHndl);
    
    
    // Create a pointer to the RecHits and raw SuperClusters
    const reco::SuperClusterCollection *rawClusters = pRawSuperClusters.product();
    
    
    EcalClusterLazyTools lazyTool(iEvent, iSetup, rHInputProducerBToken_,rHInputProducerEToken_);
 
    // Define a collection of corrected SuperClusters to put back into the event
    std::auto_ptr<reco::SuperClusterCollection> corrClusters(new reco::SuperClusterCollection);
    
    //  Loop over raw clusters and make corrected ones
    reco::SuperClusterCollection::const_iterator aClus;
    for(aClus = rawClusters->begin(); aClus != rawClusters->end(); aClus++)
       {
      double theEt = aClus->energy()/cosh( aClus->eta() ) ;
      //  std::cout << " et of SC = " << theEt << std::endl;
 
      if ( theEt < etCut_ )  continue;   // cut off low pT superclusters 
      
      bool flagS = true;
      float swissCrx(0);
      
      const reco::CaloClusterPtr seed = aClus->seed();
      DetId id = lazyTool.getMaximum(*seed).first;
      const EcalRecHitCollection & rechits = *pRecHitsB;
      EcalRecHitCollection::const_iterator it = rechits.find( id );
      
      if( it != rechits.end() ) {
         ecalSevLvlAlgoHndl->severityLevel(id, rechits);
         swissCrx = EcalTools::swissCross   (id, rechits, 0.,true);
         //      std::cout << "swissCross = " << swissCrx <<std::endl;
         // std::cout << " timing = " << it->time() << std::endl;
      }
      
      if ( fabs(it->time()) > TimingCut_ ) {
         flagS = false;
         //      std::cout << " timing = " << it->time() << std::endl;
         //   std::cout << " timing is bad........" << std::endl; 
      }
      if ( swissCrx > (float)swissCutThr_ ) {
         flagS = false ;     // swissCross cut
         //      std::cout << "swissCross = " << swissCrx <<std::endl;   
         //   std::cout << " removed by swiss cross cut" << std::endl;
      }
      // - kGood        --> good channel
      // - kProblematic --> problematic (e.g. noisy)
      // - kRecovered   --> recovered (e.g. an originally dead or saturated)
      // - kTime        --> the channel is out of time (e.g. spike)
      // - kWeird       --> weird (e.g. spike)
      // - kBad         --> bad, not suitable to be used in the reconstruction
      //   enum EcalSeverityLevel { kGood=0, kProblematic, kRecovered, kTime, kWeird, kBad };
         
      
      reco::SuperCluster newClus;
      if ( flagS == true)
         newClus=*aClus;
      else
         continue;
      corrClusters->push_back(newClus);
       }
    
    // Put collection of corrected SuperClusters into the event
    iEvent.put(corrClusters, outputCollection_);   
    
 }