#include <EcalRecHitProducer.h>

Inheritance diagram for EcalRecHitProducer:

Public Member Functions
	EcalRecHitProducer (const edm::ParameterSet &ps)

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

	~EcalRecHitProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
std::unique_ptr< EcalCleaningAlgo >	cleaningAlgo_

edm::EDGetTokenT< std::set< EBDetId > >	ebDetIdToBeRecoveredToken_

edm::EDGetTokenT< std::set< EcalTrigTowerDetId > >	ebFEToBeRecoveredToken_

std::string	ebRechitCollection_

edm::EDGetTokenT< EBUncalibratedRecHitCollection >	ebUncalibRecHitToken_

edm::EDGetTokenT< std::set< EEDetId > >	eeDetIdToBeRecoveredToken_

edm::EDGetTokenT< std::set< EcalScDetId > >	eeFEToBeRecoveredToken_

std::string	eeRechitCollection_

edm::EDGetTokenT< EEUncalibratedRecHitCollection >	eeUncalibRecHitToken_

bool	killDeadChannels_

bool	recoverEBFE_

bool	recoverEBIsolatedChannels_

bool	recoverEBVFE_

bool	recoverEEFE_

bool	recoverEEIsolatedChannels_

bool	recoverEEVFE_

std::unique_ptr< EcalRecHitWorkerBaseClass >	worker_

std::unique_ptr< EcalRecHitWorkerBaseClass >	workerRecover_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

produce ECAL rechits from uncalibrated rechits

Author: Shahram Rahatlou, University of Rome & INFN, March 2006

Definition at line 26 of file EcalRecHitProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 of file EcalRecHitProducer.cc.

References EnergyCorrector::c, cleaningAlgo_, ebDetIdToBeRecoveredToken_, ebFEToBeRecoveredToken_, ebRechitCollection_, ebUncalibRecHitToken_, eeDetIdToBeRecoveredToken_, eeFEToBeRecoveredToken_, eeRechitCollection_, eeUncalibRecHitToken_, reco::get(), edm::ParameterSet::getParameter(), killDeadChannels_, recoverEBFE_, recoverEBIsolatedChannels_, recoverEBVFE_, recoverEEFE_, recoverEEIsolatedChannels_, recoverEEVFE_, AlCaHLTBitMon_QueryRunRegistry::string, worker_, workerRecover_, and ~EcalRecHitProducer().

 {
        
         ebRechitCollection_        = ps.getParameter<std::string>("EBrechitCollection");
         eeRechitCollection_        = ps.getParameter<std::string>("EErechitCollection");
 
         recoverEBIsolatedChannels_   = ps.getParameter<bool>("recoverEBIsolatedChannels");
         recoverEEIsolatedChannels_   = ps.getParameter<bool>("recoverEEIsolatedChannels");
         recoverEBVFE_                = ps.getParameter<bool>("recoverEBVFE");
         recoverEEVFE_                = ps.getParameter<bool>("recoverEEVFE");
         recoverEBFE_                 = ps.getParameter<bool>("recoverEBFE");
         recoverEEFE_                 = ps.getParameter<bool>("recoverEEFE");
         killDeadChannels_            = ps.getParameter<bool>("killDeadChannels");
 
    
         produces< EBRecHitCollection >(ebRechitCollection_);
         produces< EERecHitCollection >(eeRechitCollection_);
 
     
     ebUncalibRecHitToken_ = 
       consumes<EBUncalibratedRecHitCollection>( ps.getParameter<edm::InputTag>("EBuncalibRecHitCollection"));
     
     eeUncalibRecHitToken_ = 
       consumes<EEUncalibratedRecHitCollection>( ps.getParameter<edm::InputTag>("EEuncalibRecHitCollection"));
 
         ebDetIdToBeRecoveredToken_ = 
       consumes<std::set<EBDetId>>(ps.getParameter<edm::InputTag>("ebDetIdToBeRecovered"));  
     
     eeDetIdToBeRecoveredToken_=        
       consumes<std::set<EEDetId>>(ps.getParameter<edm::InputTag>("eeDetIdToBeRecovered"));
 
     ebFEToBeRecoveredToken_ = consumes<std::set<EcalTrigTowerDetId>>(ps.getParameter<edm::InputTag>("ebFEToBeRecovered"));
 
         eeFEToBeRecoveredToken_= consumes<std::set<EcalScDetId>>( ps.getParameter<edm::InputTag>("eeFEToBeRecovered"))   ;
 
         std::string componentType = ps.getParameter<std::string>("algo");
     edm::ConsumesCollector c{consumesCollector()};
         worker_ = std::unique_ptr<EcalRecHitWorkerBaseClass>{EcalRecHitWorkerFactory::get()->create(componentType, ps, c)};
 
         // to recover problematic channels
         componentType = ps.getParameter<std::string>("algoRecover");
         workerRecover_ = std::unique_ptr<EcalRecHitWorkerBaseClass>{EcalRecHitWorkerFactory::get()->create(componentType, ps, c)};
 
     edm::ParameterSet cleaningPs = 
       ps.getParameter<edm::ParameterSet>("cleaningConfig");
     cleaningAlgo_ = std::make_unique<EcalCleaningAlgo>(cleaningPs);
 }

EcalRecHitProducer::~EcalRecHitProducer ( )

overridedefault

Referenced by EcalRecHitProducer().

Member Function Documentation

void EcalRecHitProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 269 of file EcalRecHitProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, edm::ParameterSetDescription::reserve(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                     {
   edm::ParameterSetDescription desc;
   desc.add<bool>("recoverEEVFE",false);
   desc.add<std::string>("EErechitCollection","EcalRecHitsEE");
   desc.add<bool>("recoverEBIsolatedChannels",false);
   desc.add<bool>("recoverEBVFE",false);
   desc.add<bool>("laserCorrection",true);
   desc.add<double>("EBLaserMIN",0.5);
   desc.add<bool>("killDeadChannels",true);
   {
     std::vector<int> temp1;
     temp1.reserve(3);
     temp1.push_back(14);
     temp1.push_back(78);
     temp1.push_back(142);
     desc.add<std::vector<int> >("dbStatusToBeExcludedEB",temp1);
   }
   desc.add<edm::InputTag>("EEuncalibRecHitCollection",edm::InputTag("ecalMultiFitUncalibRecHit","EcalUncalibRecHitsEE"));
   {
     std::vector<int> temp1;
     temp1.reserve(3);
     temp1.push_back(14);
     temp1.push_back(78);
     temp1.push_back(142);
     desc.add<std::vector<int> >("dbStatusToBeExcludedEE",temp1);
   }
   desc.add<double>("EELaserMIN",0.5);
   desc.add<edm::InputTag>("ebFEToBeRecovered",edm::InputTag("ecalDetIdToBeRecovered","ebFE"));
   {
     edm::ParameterSetDescription psd0;
     psd0.add<double>("e6e2thresh",0.04);
     psd0.add<double>("tightenCrack_e6e2_double",3);
     psd0.add<double>("e4e1Threshold_endcap",0.3);
     psd0.add<double>("tightenCrack_e4e1_single",3);
     psd0.add<double>("tightenCrack_e1_double",2);
     psd0.add<double>("cThreshold_barrel",4);
     psd0.add<double>("e4e1Threshold_barrel",0.08);
     psd0.add<double>("tightenCrack_e1_single",2);
     psd0.add<double>("e4e1_b_barrel",-0.024);
     psd0.add<double>("e4e1_a_barrel",0.04);
     psd0.add<double>("ignoreOutOfTimeThresh",1000000000.0);
     psd0.add<double>("cThreshold_endcap",15);
     psd0.add<double>("e4e1_b_endcap",-0.0125);
     psd0.add<double>("e4e1_a_endcap",0.02);
     psd0.add<double>("cThreshold_double",10);
     desc.add<edm::ParameterSetDescription>("cleaningConfig",psd0);
   }
   desc.add<double>("logWarningEtThreshold_EE_FE",50);
   desc.add<edm::InputTag>("eeDetIdToBeRecovered",edm::InputTag("ecalDetIdToBeRecovered","eeDetId"));
   desc.add<bool>("recoverEBFE",true);
   desc.add<edm::InputTag>("eeFEToBeRecovered",edm::InputTag("ecalDetIdToBeRecovered","eeFE"));
   desc.add<edm::InputTag>("ebDetIdToBeRecovered",edm::InputTag("ecalDetIdToBeRecovered","ebDetId"));
   desc.add<double>("singleChannelRecoveryThreshold",8);
   {
     std::vector<std::string> temp1;
     temp1.reserve(9);
     temp1.push_back("kNoisy");
     temp1.push_back("kNNoisy");
     temp1.push_back("kFixedG6");
     temp1.push_back("kFixedG1");
     temp1.push_back("kFixedG0");
     temp1.push_back("kNonRespondingIsolated");
     temp1.push_back("kDeadVFE");
     temp1.push_back("kDeadFE");
     temp1.push_back("kNoDataNoTP");
     desc.add<std::vector<std::string> >("ChannelStatusToBeExcluded",temp1);
   }
   desc.add<std::string>("EBrechitCollection","EcalRecHitsEB");
   desc.add<edm::InputTag>("triggerPrimitiveDigiCollection",edm::InputTag("ecalDigis","EcalTriggerPrimitives"));
   desc.add<bool>("recoverEEFE",true);
   desc.add<std::string>("singleChannelRecoveryMethod","NeuralNetworks");
   desc.add<double>("EBLaserMAX",3.0);
   {
     edm::ParameterSetDescription psd0;
     {
       std::vector<std::string> temp2;
       temp2.reserve(4);
       temp2.push_back("kOk");
       temp2.push_back("kDAC");
       temp2.push_back("kNoLaser");
       temp2.push_back("kNoisy");
       psd0.add<std::vector<std::string> >("kGood",temp2);
     }
     {
       std::vector<std::string> temp2;
       temp2.reserve(3);
       temp2.push_back("kFixedG0");
       temp2.push_back("kNonRespondingIsolated");
       temp2.push_back("kDeadVFE");
       psd0.add<std::vector<std::string> >("kNeighboursRecovered",temp2);
     }
     {
       std::vector<std::string> temp2;
       temp2.reserve(1);
       temp2.push_back("kNoDataNoTP");
       psd0.add<std::vector<std::string> >("kDead",temp2);
     }
     {
       std::vector<std::string> temp2;
       temp2.reserve(3);
       temp2.push_back("kNNoisy");
       temp2.push_back("kFixedG6");
       temp2.push_back("kFixedG1");
       psd0.add<std::vector<std::string> >("kNoisy",temp2);
     }
     {
       std::vector<std::string> temp2;
       temp2.reserve(1);
       temp2.push_back("kDeadFE");
       psd0.add<std::vector<std::string> >("kTowerRecovered",temp2);
     }
     desc.add<edm::ParameterSetDescription>("flagsMapDBReco",psd0);
   }
   desc.add<edm::InputTag>("EBuncalibRecHitCollection",edm::InputTag("ecalMultiFitUncalibRecHit","EcalUncalibRecHitsEB"));
   desc.add<std::string>("algoRecover","EcalRecHitWorkerRecover");
   desc.add<std::string>("algo","EcalRecHitWorkerSimple");
   desc.add<double>("EELaserMAX",8.0);
   desc.add<double>("logWarningEtThreshold_EB_FE",50);
   desc.add<bool>("recoverEEIsolatedChannels",false);
   desc.add<bool>("skipTimeCalib",false);
   descriptions.add("ecalRecHit",desc);
 }

void EcalRecHitProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 81 of file EcalRecHitProducer.cc.

References edm::SortedCollection< T, SORT >::begin(), cleaningAlgo_, EcalRecHitWorkerBaseClass::EB_FE, EcalRecHitWorkerBaseClass::EB_single, EcalRecHitWorkerBaseClass::EB_VFE, ebDetIdToBeRecoveredToken_, ebFEToBeRecoveredToken_, ebRechitCollection_, ebUncalibRecHitToken_, EcalRecHitWorkerBaseClass::EE_FE, EcalRecHitWorkerBaseClass::EE_single, EcalRecHitWorkerBaseClass::EE_VFE, eeDetIdToBeRecoveredToken_, eeFEToBeRecoveredToken_, eeRechitCollection_, eeUncalibRecHitToken_, EcalCondObjectContainer< T >::end(), edm::SortedCollection< T, SORT >::end(), EBDetId::ETAPHIMODE, EcalCondObjectContainer< T >::find(), edm::EventSetup::get(), edm::Event::getByToken(), EcalChannelStatusCode::getStatusCode(), EcalChannelStatusCode::kDeadVFE, killDeadChannels_, LogDebug, eostools::move(), edm::Handle< T >::product(), edm::Event::put(), recoverEBFE_, recoverEBIsolatedChannels_, recoverEBVFE_, recoverEEFE_, recoverEEIsolatedChannels_, recoverEEVFE_, edm::SortedCollection< T, SORT >::size(), EEDetId::validDetId(), worker_, and workerRecover_.

 {
         using namespace edm;
 
         Handle< EBUncalibratedRecHitCollection > pEBUncalibRecHits;
         Handle< EEUncalibratedRecHitCollection > pEEUncalibRecHits;
 
         const EBUncalibratedRecHitCollection*  ebUncalibRecHits = nullptr;
         const EEUncalibratedRecHitCollection*  eeUncalibRecHits = nullptr; 
 
         // get the barrel uncalib rechit collection
        
     evt.getByToken( ebUncalibRecHitToken_, pEBUncalibRecHits);
     ebUncalibRecHits = pEBUncalibRecHits.product();
     LogDebug("EcalRecHitDebug") << "total # EB uncalibrated rechits: " << ebUncalibRecHits->size();
         
 
        
     evt.getByToken( eeUncalibRecHitToken_, pEEUncalibRecHits);
     eeUncalibRecHits = pEEUncalibRecHits.product(); // get a ptr to the product
     LogDebug("EcalRecHitDebug") << "total # EE uncalibrated rechits: " << eeUncalibRecHits->size();
        
         // collection of rechits to put in the event
         auto ebRecHits = std::make_unique<EBRecHitCollection>();
         auto eeRecHits = std::make_unique<EERecHitCollection>();
 
         worker_->set(es);
 
         if ( recoverEBIsolatedChannels_ || recoverEEIsolatedChannels_
                 || recoverEBFE_ || recoverEEFE_
                 || recoverEBVFE_ || recoverEEVFE_
                 || killDeadChannels_ ) {
                 workerRecover_->set(es);
         }
 
         if (ebUncalibRecHits)
         {
                 // loop over uncalibrated rechits to make calibrated ones
                 for(EBUncalibratedRecHitCollection::const_iterator it  = ebUncalibRecHits->begin(); it != ebUncalibRecHits->end(); ++it) {
                         worker_->run(evt, *it, *ebRecHits);
                 }
         }
 
         if (eeUncalibRecHits)
         {
                 // loop over uncalibrated rechits to make calibrated ones
                 for(EEUncalibratedRecHitCollection::const_iterator it  = eeUncalibRecHits->begin(); it != eeUncalibRecHits->end(); ++it) {
                         worker_->run(evt, *it, *eeRecHits);
                 }
         }
 
         // sort collections before attempting recovery, to avoid insertion of double recHits
         ebRecHits->sort();
         eeRecHits->sort();
         
         if ( recoverEBIsolatedChannels_ || recoverEBFE_ || killDeadChannels_ )
         {
                 edm::Handle< std::set<EBDetId> > pEBDetId;
                 const std::set<EBDetId> * detIds = nullptr;
         evt.getByToken( ebDetIdToBeRecoveredToken_, pEBDetId);
         detIds = pEBDetId.product();
                        
 
                 if ( detIds ) {
                         edm::ESHandle<EcalChannelStatus> chStatus;
                         es.get<EcalChannelStatusRcd>().get(chStatus);
                         for( std::set<EBDetId>::const_iterator it = detIds->begin(); it != detIds->end(); ++it ) {
                                 // get channel status map to treat dead VFE separately
                                 EcalChannelStatusMap::const_iterator chit = chStatus->find( *it );
                                 EcalChannelStatusCode chStatusCode;
                                 if ( chit != chStatus->end() ) {
                                         chStatusCode = *chit;
                                 } else {
                                         edm::LogError("EcalRecHitProducerError") << "No channel status found for xtal "
                                                 << (*it).rawId()
                                                 << "! something wrong with EcalChannelStatus in your DB? ";
                                 }
                                 EcalUncalibratedRecHit urh;
                                 if ( chStatusCode.getStatusCode()  == EcalChannelStatusCode::kDeadVFE ) { // dead VFE (from DB info)
                                         // uses the EcalUncalibratedRecHit to pass the DetId info
                                         urh = EcalUncalibratedRecHit( *it, 0, 0, 0, 0, EcalRecHitWorkerBaseClass::EB_VFE );
                                         if ( recoverEBVFE_ || killDeadChannels_ ) workerRecover_->run( evt, urh, *ebRecHits );
                                 } else {
                                         // uses the EcalUncalibratedRecHit to pass the DetId info
                                         urh = EcalUncalibratedRecHit( *it, 0, 0, 0, 0, EcalRecHitWorkerBaseClass::EB_single );
                                         if ( recoverEBIsolatedChannels_ || killDeadChannels_ ) workerRecover_->run( evt, urh, *ebRecHits );
                                 }
                                 
                         }
                 }
         }
 
         if ( recoverEEIsolatedChannels_ || recoverEEVFE_ || killDeadChannels_ )
         {
                 edm::Handle< std::set<EEDetId> > pEEDetId;
                 const std::set<EEDetId> * detIds = nullptr;
                
         evt.getByToken( eeDetIdToBeRecoveredToken_, pEEDetId);
         detIds = pEEDetId.product();
             
                 if ( detIds ) {
                         edm::ESHandle<EcalChannelStatus> chStatus;
                         es.get<EcalChannelStatusRcd>().get(chStatus);
                         for( std::set<EEDetId>::const_iterator it = detIds->begin(); it != detIds->end(); ++it ) {
                                 // get channel status map to treat dead VFE separately
                                 EcalChannelStatusMap::const_iterator chit = chStatus->find( *it );
                                 EcalChannelStatusCode chStatusCode;
                                 if ( chit != chStatus->end() ) {
                                         chStatusCode = *chit;
                                 } else {
                                         edm::LogError("EcalRecHitProducerError") << "No channel status found for xtal "
                                                 << (*it).rawId()
                                                 << "! something wrong with EcalChannelStatus in your DB? ";
                                 }
                                 EcalUncalibratedRecHit urh;
                                 if ( chStatusCode.getStatusCode()  == EcalChannelStatusCode::kDeadVFE) { // dead VFE (from DB info)
                                         // uses the EcalUncalibratedRecHit to pass the DetId info
                                         urh = EcalUncalibratedRecHit( *it, 0, 0, 0, 0, EcalRecHitWorkerBaseClass::EE_VFE );
                                         if ( recoverEEVFE_ || killDeadChannels_ ) workerRecover_->run( evt, urh, *eeRecHits );
                                 } else {
                                         // uses the EcalUncalibratedRecHit to pass the DetId info
                                         urh = EcalUncalibratedRecHit( *it, 0, 0, 0, 0, EcalRecHitWorkerBaseClass::EE_single );
                                         if ( recoverEEIsolatedChannels_ || killDeadChannels_ ) workerRecover_->run( evt, urh, *eeRecHits );
                                 }
                         }
                 }
         }
 
         if ( recoverEBFE_ || killDeadChannels_ )
         {
                 edm::Handle< std::set<EcalTrigTowerDetId> > pEBFEId;
                 const std::set<EcalTrigTowerDetId> * ttIds = nullptr;
                
         evt.getByToken( ebFEToBeRecoveredToken_, pEBFEId);
         ttIds = pEBFEId.product();
                 
                 if ( ttIds ) {
                         for( std::set<EcalTrigTowerDetId>::const_iterator it = ttIds->begin(); it != ttIds->end(); ++it ) {
                                 // uses the EcalUncalibratedRecHit to pass the DetId info
                                 int ieta = (((*it).ietaAbs()-1)*5+1)*(*it).zside(); // from EcalTrigTowerConstituentsMap
                                 int iphi = (((*it).iphi()-1)*5+11)%360;             // from EcalTrigTowerConstituentsMap
                                 if( iphi <= 0 ) iphi += 360;                        // from EcalTrigTowerConstituentsMap
                                 EcalUncalibratedRecHit urh( EBDetId(ieta, iphi, EBDetId::ETAPHIMODE), 0, 0, 0, 0, EcalRecHitWorkerBaseClass::EB_FE );
                                 workerRecover_->run( evt, urh, *ebRecHits );
                         }
                 }
         }
 
         if ( recoverEEFE_ || killDeadChannels_ )
         {
                 edm::Handle< std::set<EcalScDetId> > pEEFEId;
                 const std::set<EcalScDetId> * scIds = nullptr;
             
         evt.getByToken( eeFEToBeRecoveredToken_, pEEFEId);
         scIds = pEEFEId.product();
     
             
                 if ( scIds ) {
                         for( std::set<EcalScDetId>::const_iterator it = scIds->begin(); it != scIds->end(); ++it ) {
                                 // uses the EcalUncalibratedRecHit to pass the DetId info
                                 if (EEDetId::validDetId( ((*it).ix()-1)*5+1, ((*it).iy()-1)*5+1, (*it).zside() )) {
                                         EcalUncalibratedRecHit urh( EEDetId( ((*it).ix()-1)*5+1, ((*it).iy()-1)*5+1, (*it).zside() ), 0, 0, 0, 0, EcalRecHitWorkerBaseClass::EE_FE );
                                         workerRecover_->run( evt, urh, *eeRecHits );
                                 }
                         }
                 }
         }
 
     // without re-sorting, find (used below in cleaning) will lead
         // to undefined results
     ebRecHits->sort();
         eeRecHits->sort();
     
     // apply spike cleaning
     if (cleaningAlgo_){
       cleaningAlgo_->setFlags(*ebRecHits);
       cleaningAlgo_->setFlags(*eeRecHits);
     }
 
 
         // put the collection of recunstructed hits in the event   
         LogInfo("EcalRecHitInfo") << "total # EB calibrated rechits: " << ebRecHits->size();
         LogInfo("EcalRecHitInfo") << "total # EE calibrated rechits: " << eeRecHits->size();
 
         evt.put(std::move(ebRecHits), ebRechitCollection_);
         evt.put(std::move(eeRecHits), eeRechitCollection_);
 }