#include <EcalRecHitWorkerSimple.h>

Inheritance diagram for EcalRecHitWorkerSimple:

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

	EcalRecHitWorkerSimple (const edm::ParameterSet &, edm::ConsumesCollector &c)

bool	run (const edm::Event &evt, const EcalUncalibratedRecHit &uncalibRH, EcalRecHitCollection &result) override

void	set (const edm::EventSetup &es) override

	~EcalRecHitWorkerSimple () override

Public Member Functions inherited from EcalRecHitWorkerBaseClass
	EcalRecHitWorkerBaseClass (const edm::ParameterSet &)

	EcalRecHitWorkerBaseClass (const edm::ParameterSet &, edm::ConsumesCollector &c)

virtual	~EcalRecHitWorkerBaseClass ()

Protected Member Functions
uint32_t	setFlagBits (const std::vector< std::vector< uint32_t > > &map, const uint32_t &status)

Protected Attributes
edm::ESHandle< EcalADCToGeVConstant >	agc

edm::ESHandle< EcalChannelStatus >	chStatus

double	EBLaserMAX_

double	EBLaserMIN_

double	EELaserMAX_

double	EELaserMIN_

uint32_t	flagmask_

edm::ESHandle< EcalIntercalibConstants >	ical

edm::ESHandle< EcalTimeCalibConstants >	itime

bool	killDeadChannels_

edm::ESHandle< EcalLaserDbService >	laser

bool	laserCorrection_

edm::ESHandle< EcalTimeOffsetConstant >	offtime

EcalRecHitSimpleAlgo *	rechitMaker_

bool	skipTimeCalib_

std::vector< int >	v_chstatus_

std::vector< std::vector< uint32_t > >	v_DB_reco_flags_

Additional Inherited Members
Public Types inherited from EcalRecHitWorkerBaseClass
enum	RecoveryFlags { EB_single, EE_single, EB_VFE, EE_VFE, EB_FE, EE_FE }

Detailed Description

Definition at line 23 of file EcalRecHitWorkerSimple.h.

Constructor & Destructor Documentation

◆ EcalRecHitWorkerSimple() [1/2]

EcalRecHitWorkerSimple::EcalRecHitWorkerSimple	(	const edm::ParameterSet &	ps,
		edm::ConsumesCollector &	c
	)

Definition at line 15 of file EcalRecHitWorkerSimple.cc.

     : EcalRecHitWorkerBaseClass(ps, c) {
   rechitMaker_ = new EcalRecHitSimpleAlgo();
   v_chstatus_ = StringToEnumValue<EcalChannelStatusCode::Code>(
       ps.getParameter<std::vector<std::string> >("ChannelStatusToBeExcluded"));
   killDeadChannels_ = ps.getParameter<bool>("killDeadChannels");
   laserCorrection_ = ps.getParameter<bool>("laserCorrection");
   EBLaserMIN_ = ps.getParameter<double>("EBLaserMIN");
   EELaserMIN_ = ps.getParameter<double>("EELaserMIN");
   EBLaserMAX_ = ps.getParameter<double>("EBLaserMAX");
   EELaserMAX_ = ps.getParameter<double>("EELaserMAX");
  
   skipTimeCalib_ = ps.getParameter<bool>("skipTimeCalib");
  
   // Traslate string representation of flagsMapDBReco into enum values
   const edm::ParameterSet& p = ps.getParameter<edm::ParameterSet>("flagsMapDBReco");
   std::vector<std::string> recoflagbitsStrings = p.getParameterNames();
   v_DB_reco_flags_.resize(32);
  
   for (unsigned int i = 0; i != recoflagbitsStrings.size(); ++i) {
     EcalRecHit::Flags recoflagbit = (EcalRecHit::Flags)StringToEnumValue<EcalRecHit::Flags>(recoflagbitsStrings[i]);
     std::vector<std::string> dbstatus_s = p.getParameter<std::vector<std::string> >(recoflagbitsStrings[i]);
     std::vector<uint32_t> dbstatuses;
     for (unsigned int j = 0; j != dbstatus_s.size(); ++j) {
       EcalChannelStatusCode::Code dbstatus =
           (EcalChannelStatusCode::Code)StringToEnumValue<EcalChannelStatusCode::Code>(dbstatus_s[j]);
       dbstatuses.push_back(dbstatus);
     }
  
     v_DB_reco_flags_[recoflagbit] = dbstatuses;
   }
  
   flagmask_ = 0;
   flagmask_ |= 0x1 << EcalRecHit::kNeighboursRecovered;
   flagmask_ |= 0x1 << EcalRecHit::kTowerRecovered;
   flagmask_ |= 0x1 << EcalRecHit::kDead;
   flagmask_ |= 0x1 << EcalRecHit::kKilled;
   flagmask_ |= 0x1 << EcalRecHit::kTPSaturated;
   flagmask_ |= 0x1 << EcalRecHit::kL1SpikeFlag;
 }

References EBLaserMAX_, EBLaserMIN_, EELaserMAX_, EELaserMIN_, flagmask_, edm::ParameterSet::getParameter(), mps_fire::i, dqmiolumiharvest::j, EcalRecHit::kDead, killDeadChannels_, EcalRecHit::kKilled, EcalRecHit::kL1SpikeFlag, EcalRecHit::kNeighboursRecovered, EcalRecHit::kTowerRecovered, EcalRecHit::kTPSaturated, laserCorrection_, AlCaHLTBitMon_ParallelJobs::p, rechitMaker_, skipTimeCalib_, v_chstatus_, and v_DB_reco_flags_.

◆ EcalRecHitWorkerSimple() [2/2]

EcalRecHitWorkerSimple::EcalRecHitWorkerSimple ( const edm::ParameterSet & )

◆ ~EcalRecHitWorkerSimple()

EcalRecHitWorkerSimple::~EcalRecHitWorkerSimple ( )

override

Definition at line 157 of file EcalRecHitWorkerSimple.cc.

157 { delete rechitMaker_; }

References rechitMaker_.

Member Function Documentation

◆ run()

bool EcalRecHitWorkerSimple::run	(	const edm::Event &	evt,
		const EcalUncalibratedRecHit &	uncalibRH,
		EcalRecHitCollection &	result
	)

overridevirtual

Implements EcalRecHitWorkerBaseClass.

Definition at line 70 of file EcalRecHitWorkerSimple.cc.

                                                                {
   DetId detid = uncalibRH.id();
  
   EcalChannelStatusMap::const_iterator chit = chStatus->find(detid);
   EcalChannelStatusCode::Code dbstatus = chit->getStatusCode();
  
   // check for channels to be excluded from reconstruction
   if (!v_chstatus_.empty()) {
     std::vector<int>::const_iterator res = std::find(v_chstatus_.begin(), v_chstatus_.end(), dbstatus);
     if (res != v_chstatus_.end())
       return false;
   }
  
   uint32_t flagBits = setFlagBits(v_DB_reco_flags_, dbstatus);
  
   float offsetTime = 0;  // the global time phase
   const EcalIntercalibConstantMap& icalMap = ical->getMap();
   if (detid.subdetId() == EcalEndcap) {
     rechitMaker_->setADCToGeVConstant(float(agc->getEEValue()));
     if (!skipTimeCalib_)
       offsetTime = offtime->getEEValue();
   } else {
     rechitMaker_->setADCToGeVConstant(float(agc->getEBValue()));
     if (!skipTimeCalib_)
       offsetTime = offtime->getEBValue();
   }
  
   // first intercalibration constants
   EcalIntercalibConstantMap::const_iterator icalit = icalMap.find(detid);
   EcalIntercalibConstant icalconst = 1;
   if (icalit != icalMap.end()) {
     icalconst = (*icalit);
   } else {
     edm::LogError("EcalRecHitError") << "No intercalib const found for xtal " << detid.rawId()
                                      << "! something wrong with EcalIntercalibConstants in your DB? ";
   }
  
   // get laser coefficient
   float lasercalib = 1.;
   if (laserCorrection_)
     lasercalib = laser->getLaserCorrection(detid, evt.time());
  
   // get time calibration coefficient
   EcalTimeCalibConstant itimeconst = 0;
  
   if (!skipTimeCalib_) {
     const EcalTimeCalibConstantMap& itimeMap = itime->getMap();
     EcalTimeCalibConstantMap::const_iterator itime = itimeMap.find(detid);
  
     if (itime != itimeMap.end()) {
       itimeconst = (*itime);
     } else {
       edm::LogError("EcalRecHitError") << "No time calib const found for xtal " << detid.rawId()
                                        << "! something wrong with EcalTimeCalibConstants in your DB? ";
     }
   }
  
   // make the rechit and put in the output collection, unless recovery has to take care of it
   if (!(flagmask_ & flagBits) || !killDeadChannels_) {
     EcalRecHit myrechit(rechitMaker_->makeRecHit(uncalibRH,
                                                  icalconst * lasercalib,
                                                  (itimeconst + offsetTime),
                                                  /*recoflags_ 0*/
                                                  flagBits));
  
     if (detid.subdetId() == EcalBarrel && (lasercalib < EBLaserMIN_ || lasercalib > EBLaserMAX_))
       myrechit.setFlag(EcalRecHit::kPoorCalib);
     if (detid.subdetId() == EcalEndcap && (lasercalib < EELaserMIN_ || lasercalib > EELaserMAX_))
       myrechit.setFlag(EcalRecHit::kPoorCalib);
     result.push_back(myrechit);
   }
  
   return true;
 }

References agc, chStatus, EBLaserMAX_, EcalBarrel, EcalEndcap, EELaserMAX_, EcalCondObjectContainer< T >::end(), spr::find(), EcalCondObjectContainer< T >::find(), flagmask_, EcalADCToGeVConstant::getEBValue(), EcalTimeOffsetConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), EcalTimeOffsetConstant::getEEValue(), EcalLaserDbService::getLaserCorrection(), ical, EcalUncalibratedRecHit::id(), itime, killDeadChannels_, EcalRecHit::kPoorCalib, laser, laserCorrection_, EcalRecHitSimpleAlgo::makeRecHit(), offtime, DetId::rawId(), rechitMaker_, mps_fire::result, EcalRecHitSimpleAlgo::setADCToGeVConstant(), EcalRecHit::setFlag(), setFlagBits(), skipTimeCalib_, DetId::subdetId(), edm::EventBase::time(), v_chstatus_, and v_DB_reco_flags_.

◆ set()

void EcalRecHitWorkerSimple::set ( const edm::EventSetup & es )

overridevirtual

Implements EcalRecHitWorkerBaseClass.

Definition at line 56 of file EcalRecHitWorkerSimple.cc.

                                                         {
   es.get<EcalIntercalibConstantsRcd>().get(ical);
  
   if (!skipTimeCalib_) {
     es.get<EcalTimeCalibConstantsRcd>().get(itime);
     es.get<EcalTimeOffsetConstantRcd>().get(offtime);
   }
  
   es.get<EcalADCToGeVConstantRcd>().get(agc);
   es.get<EcalChannelStatusRcd>().get(chStatus);
   if (laserCorrection_)
     es.get<EcalLaserDbRecord>().get(laser);
 }

References agc, chStatus, edm::EventSetup::get(), get, ical, itime, laser, laserCorrection_, offtime, and skipTimeCalib_.

◆ setFlagBits()

uint32_t EcalRecHitWorkerSimple::setFlagBits	(	const std::vector< std::vector< uint32_t > > &	map,
		const uint32_t &	status
	)

protected

Definition at line 148 of file EcalRecHitWorkerSimple.cc.

                                                                                                                {
   for (unsigned int i = 0; i != map.size(); ++i) {
     if (std::find(map[i].begin(), map[i].end(), status) != map[i].end())
       return 0x1 << i;
   }
  
   return 0;
 }

References begin, end, spr::find(), mps_fire::i, genParticles_cff::map, mps_update::status, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by run().

Member Data Documentation

◆ agc

edm::ESHandle<EcalADCToGeVConstant> EcalRecHitWorkerSimple::agc

protected

Definition at line 41 of file EcalRecHitWorkerSimple.h.

Referenced by run(), and set().

◆ chStatus

edm::ESHandle<EcalChannelStatus> EcalRecHitWorkerSimple::chStatus

protected

Definition at line 42 of file EcalRecHitWorkerSimple.h.

Referenced by run(), and set().

◆ EBLaserMAX_

double EcalRecHitWorkerSimple::EBLaserMAX_

protected

Definition at line 35 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), and run().

◆ EBLaserMIN_

double EcalRecHitWorkerSimple::EBLaserMIN_

protected

Definition at line 33 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple().

◆ EELaserMAX_

double EcalRecHitWorkerSimple::EELaserMAX_

protected

Definition at line 36 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), and run().

◆ EELaserMIN_

double EcalRecHitWorkerSimple::EELaserMIN_

protected

Definition at line 34 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple().

◆ flagmask_

uint32_t EcalRecHitWorkerSimple::flagmask_

protected

Definition at line 51 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), and run().

◆ ical

edm::ESHandle<EcalIntercalibConstants> EcalRecHitWorkerSimple::ical

protected

Definition at line 38 of file EcalRecHitWorkerSimple.h.

Referenced by run(), and set().

◆ itime

edm::ESHandle<EcalTimeCalibConstants> EcalRecHitWorkerSimple::itime

protected

Definition at line 39 of file EcalRecHitWorkerSimple.h.

Referenced by run(), and set().

◆ killDeadChannels_

bool EcalRecHitWorkerSimple::killDeadChannels_

protected

Definition at line 53 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), and run().

◆ laser

edm::ESHandle<EcalLaserDbService> EcalRecHitWorkerSimple::laser

protected

Definition at line 44 of file EcalRecHitWorkerSimple.h.

Referenced by run(), and set().

◆ laserCorrection_

bool EcalRecHitWorkerSimple::laserCorrection_

protected

Definition at line 54 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), run(), and set().

◆ offtime

edm::ESHandle<EcalTimeOffsetConstant> EcalRecHitWorkerSimple::offtime

protected

Definition at line 40 of file EcalRecHitWorkerSimple.h.

Referenced by run(), and set().

◆ rechitMaker_

EcalRecHitSimpleAlgo* EcalRecHitWorkerSimple::rechitMaker_

protected

Definition at line 57 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), run(), and ~EcalRecHitWorkerSimple().

◆ skipTimeCalib_

bool EcalRecHitWorkerSimple::skipTimeCalib_

protected

Definition at line 55 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), run(), and set().

◆ v_chstatus_

std::vector<int> EcalRecHitWorkerSimple::v_chstatus_

protected

Definition at line 43 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), and run().

◆ v_DB_reco_flags_

std::vector<std::vector<uint32_t> > EcalRecHitWorkerSimple::v_DB_reco_flags_

protected

Definition at line 47 of file EcalRecHitWorkerSimple.h.

Referenced by EcalRecHitWorkerSimple(), and run().

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalRecHitWorkerSimple() [1/2]

◆ EcalRecHitWorkerSimple() [2/2]

◆ ~EcalRecHitWorkerSimple()

Member Function Documentation

◆ run()

◆ set()

◆ setFlagBits()

Member Data Documentation

◆ agc

◆ chStatus

◆ EBLaserMAX_

◆ EBLaserMIN_

◆ EELaserMAX_

◆ EELaserMIN_

◆ flagmask_

◆ ical

◆ itime

◆ killDeadChannels_

◆ laser

◆ laserCorrection_

◆ offtime

◆ rechitMaker_

◆ skipTimeCalib_

◆ v_chstatus_

◆ v_DB_reco_flags_