Inheritance diagram for EcalPhiSymRecHitProducerBase:

Public Member Functions
	EcalPhiSymRecHitProducerBase (const edm::ParameterSet &pSet, edm::ConsumesCollector &&cc)

void	initializeConfigCache (edm::EventSetup const &setup, edm::ESGetToken< CaloGeometry, CaloGeometryRecord > const &geoToken, std::shared_ptr< ConfigCache > &cache) const

void	initializeJob ()

void	initializePhiSymCache (edm::EventSetup const &setup, edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > const &chStatusToken, ConfigCache const *config, std::shared_ptr< PhiSymCache > &cache) const

void	initializeStreamCache (ConfigCache const config, PhiSymCache cache) const

void	processEvent (edm::Event const &event, edm::EventSetup const &setup, ConfigCache const config, PhiSymCache cache) const

void	sumCache (PhiSymCache summaryc, PhiSymCache streamc) const

	~EcalPhiSymRecHitProducerBase ()

Protected Attributes
std::vector< double >	A_

std::vector< double >	B_

EcalRingCalibrationTools	calibRing_

edm::EDGetTokenT < EBRecHitCollection >	ebToken_

edm::EDGetTokenT < EBRecHitCollection >	eeToken_

float	etCutEB_

float	etCutEE_

std::vector< double >	eThresholdsEB_

float	eThresholdsEE_ [kNRingsEE]

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geoToken_

edm::ESGetToken < EcalLaserDbService, EcalLaserDbRecord >	laserDbToken_

std::vector< double >	misCalibRangeEB_

std::vector< double >	misCalibRangeEE_

std::vector< float >	misCalibStepsEB_

std::vector< float >	misCalibStepsEE_

int	nMisCalib_

int	nSumEtValues_

float	thrEEmod_

Static Protected Attributes
static const short	kNRingsEB = EcalRingCalibrationTools::N_RING_BARREL

static const short	kNRingsEE = EcalRingCalibrationTools::N_RING_ENDCAP

static const short	ringsInOneEE = kNRingsEE / 2

Detailed Description

Definition at line 69 of file EcalPhiSymRecHitProducers.cc.

Constructor & Destructor Documentation

EcalPhiSymRecHitProducerBase::EcalPhiSymRecHitProducerBase	(	const edm::ParameterSet &	pSet,
		edm::ConsumesCollector &&	cc
	)

explicit

Definition at line 121 of file EcalPhiSymRecHitProducers.cc.

     : geoToken_(cc.esConsumes()),
       laserDbToken_(cc.esConsumes()),
       ebToken_(cc.consumes<EBRecHitCollection>(pSet.getParameter<edm::InputTag>("barrelHitCollection"))),
       eeToken_(cc.consumes<EBRecHitCollection>(pSet.getParameter<edm::InputTag>("endcapHitCollection"))),
       etCutEB_(pSet.getParameter<double>("etCut_barrel")),
       eThresholdsEB_(pSet.getParameter<std::vector<double> >("eThresholds_barrel")),
       etCutEE_(pSet.getParameter<double>("etCut_endcap")),
       A_(pSet.getParameter<std::vector<double> >("A")),
       B_(pSet.getParameter<std::vector<double> >("B")),
       thrEEmod_(pSet.getParameter<double>("thrEEmod")),
       nMisCalib_(pSet.getParameter<int>("nMisCalib") / 2),
       nSumEtValues_(nMisCalib_ * 2 + 1),
       misCalibRangeEB_(pSet.getParameter<std::vector<double> >("misCalibRangeEB")),
       misCalibRangeEE_(pSet.getParameter<std::vector<double> >("misCalibRangeEE")) {}

EcalPhiSymRecHitProducerBase::~EcalPhiSymRecHitProducerBase ( )

inline

Definition at line 72 of file EcalPhiSymRecHitProducers.cc.

72 {};

Member Function Documentation

void EcalPhiSymRecHitProducerBase::initializeConfigCache	(	edm::EventSetup const &	setup,
		edm::ESGetToken< CaloGeometry, CaloGeometryRecord > const &	geoToken,
		std::shared_ptr< ConfigCache > &	cache
	)		const

Definition at line 313 of file EcalPhiSymRecHitProducers.cc.

References calibRing_, DetId::Ecal, EcalBarrel, EcalEndcap, PVValHelper::eta, etCutEB_, etCutEE_, eThresholdsEB_, eThresholdsEE_, relativeConstraints::geometry, edm::EventSetup::getData(), EcalRingCalibrationTools::getRingIndex(), gpuClustering::id, EEDetId::IX_MAX, kNRingsEB, relativeConstraints::ring, ringsInOneEE, and EcalRingCalibrationTools::setCaloGeometry().

Referenced by EcalPhiSymRecHitProducerLumi::globalBeginLuminosityBlock(), and EcalPhiSymRecHitProducerRun::globalBeginRun().

                                              {
   //---get the ecal geometry
   const auto* geometry = &setup.getData(geoToken);
   calibRing_.setCaloGeometry(geometry);
 
   const auto* barrelGeometry = geometry->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
   const auto* endcapGeometry = geometry->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
   cache->barrelDetIds = barrelGeometry->getValidDetIds(DetId::Ecal, EcalBarrel);
   cache->endcapDetIds = endcapGeometry->getValidDetIds(DetId::Ecal, EcalEndcap);
 
   for (auto& ebDetId : cache->barrelDetIds) {
     EBDetId id(ebDetId);
     int ring = calibRing_.getRingIndex(id);
     //---set etCut if first pass
     if (id.iphi() == 1) {
       auto cellGeometry = barrelGeometry->getGeometry(id);
       float eta = cellGeometry->getPosition().eta();
       cache->etCutsEB[ring] = eThresholdsEB_[ring] / cosh(eta) + etCutEB_;
     }
   }
   for (auto& eeDetId : cache->endcapDetIds) {
     EEDetId id(eeDetId);
     int ring = calibRing_.getRingIndex(id) - kNRingsEB;
     //---set eCutEE if first pass
     if (ring < ringsInOneEE && id.ix() == EEDetId::IX_MAX / 2) {
       auto cellGeometry = endcapGeometry->getGeometry(id);
       cache->etCutsEE[ring] = eThresholdsEE_[ring] / cosh(cellGeometry->getPosition().eta()) + etCutEE_;
       cache->etCutsEE[ring + ringsInOneEE] = cache->etCutsEE[ring];
     }
   }
 }

void EcalPhiSymRecHitProducerBase::initializeJob ( )

Definition at line 137 of file EcalPhiSymRecHitProducers.cc.

References A_, funct::abs(), B_, eThresholdsEE_, misCalibRangeEB_, misCalibRangeEE_, misCalibStepsEB_, misCalibStepsEE_, nMisCalib_, nSumEtValues_, ringsInOneEE, and thrEEmod_.

Referenced by EcalPhiSymRecHitProducerLumi::beginJob(), and EcalPhiSymRecHitProducerRun::beginJob().

                                                  {
   //---Compute the endcap thresholds using the provived parametric formula
   for (int iRing = 0; iRing < ringsInOneEE; ++iRing) {
     if (iRing < 30)
       eThresholdsEE_[iRing] = thrEEmod_ * (B_[0] + A_[0] * iRing) / 1000;
     else
       eThresholdsEE_[iRing] = thrEEmod_ * (B_[1] + A_[1] * iRing) / 1000;
     eThresholdsEE_[iRing + ringsInOneEE] = eThresholdsEE_[iRing];
   }
 
   //---misCalib value init (nMisCalib is half of the correct value!)
   float misCalibStepEB = std::abs(misCalibRangeEB_[1] - misCalibRangeEB_[0]) / (nMisCalib_ * 2);
   float misCalibStepEE = std::abs(misCalibRangeEE_[1] - misCalibRangeEE_[0]) / (nMisCalib_ * 2);
   misCalibStepsEB_.resize(nSumEtValues_);
   misCalibStepsEE_.resize(nSumEtValues_);
   for (int iMis = -nMisCalib_; iMis <= nMisCalib_; ++iMis) {
     //--- 0 -> 0; -i -> [1...n/2]; +i -> [n/2+1...n]
     int index = iMis > 0 ? iMis + nMisCalib_ : iMis == 0 ? 0 : iMis + nMisCalib_ + 1;
     misCalibStepsEB_[index] = iMis * misCalibStepEB;
     misCalibStepsEE_[index] = iMis * misCalibStepEE;
   }
 }

void EcalPhiSymRecHitProducerBase::initializePhiSymCache	(	edm::EventSetup const &	setup,
		edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > const &	chStatusToken,
		ConfigCache const *	config,
		std::shared_ptr< PhiSymCache > &	cache
	)		const

Definition at line 286 of file EcalPhiSymRecHitProducers.cc.

References ConfigCache::barrelDetIds, calibRing_, ConfigCache::endcapDetIds, edm::EventSetup::getData(), EcalRingCalibrationTools::getRingIndex(), gpuClustering::id, kNRingsEB, kNRingsEE, nSumEtValues_, relativeConstraints::ring, and EcalPhiSymRecHit::setEERing().

Referenced by EcalPhiSymRecHitProducerLumi::globalBeginLuminosityBlockSummary(), and EcalPhiSymRecHitProducerRun::globalBeginRunSummary().

                                              {
   cache->clear();
 
   //---get the channels status
   auto const& chStatus = setup.getData(chStatusToken);
 
   cache->recHitCollEB.resize(config->barrelDetIds.size());
   cache->recHitCollEE.resize(config->endcapDetIds.size());
   for (auto& ebDetId : config->barrelDetIds) {
     EBDetId id(ebDetId);
     cache->recHitCollEB.at(id.denseIndex()) =
         EcalPhiSymRecHit(ebDetId.rawId(), nSumEtValues_, chStatus[id].getStatusCode());
   }
   for (auto& eeDetId : config->endcapDetIds) {
     EEDetId id(eeDetId);
     int ring = calibRing_.getRingIndex(id) - kNRingsEB;
     cache->recHitCollEE.at(id.denseIndex()) =
         EcalPhiSymRecHit(eeDetId.rawId(), nSumEtValues_, chStatus[id].getStatusCode());
     cache->recHitCollEE.at(id.denseIndex())
         .setEERing(ring < kNRingsEE / 2 ? ring - kNRingsEE / 2 : ring - kNRingsEE / 2 + 1);
   }
 }

void EcalPhiSymRecHitProducerBase::initializeStreamCache	(	ConfigCache const *	config,
		PhiSymCache *	cache
	)		const

Definition at line 269 of file EcalPhiSymRecHitProducers.cc.

References ConfigCache::barrelDetIds, PhiSymCache::clear(), ConfigCache::endcapDetIds, gpuClustering::id, nSumEtValues_, PhiSymCache::recHitCollEB, and PhiSymCache::recHitCollEE.

Referenced by EcalPhiSymRecHitProducerLumi::streamBeginLuminosityBlock(), and EcalPhiSymRecHitProducerRun::streamBeginRun().

                                                                                                             {
   //---Initialize the per-stream RecHitCollection
   //   both collections are initialized to contain the total
   //   number of crystals, ordered accrodingly to the hashedIndex.
   cache->clear();
   cache->recHitCollEB.resize(config->barrelDetIds.size());
   cache->recHitCollEE.resize(config->endcapDetIds.size());
   for (auto& ebDetId : config->barrelDetIds) {
     EBDetId id(ebDetId);
     cache->recHitCollEB.at(id.denseIndex()) = EcalPhiSymRecHit(id.rawId(), nSumEtValues_);
   }
   for (auto& eeDetId : config->endcapDetIds) {
     EEDetId id(eeDetId);
     cache->recHitCollEE.at(id.denseIndex()) = EcalPhiSymRecHit(id.rawId(), nSumEtValues_);
   }
 }

void EcalPhiSymRecHitProducerBase::processEvent	(	edm::Event const &	event,
		edm::EventSetup const &	setup,
		ConfigCache const *	config,
		PhiSymCache *	cache
	)		const

Definition at line 160 of file EcalPhiSymRecHitProducers.cc.

References HLT_FULL_cff::barrelRecHits, calibRing_, EBDetId::denseIndex(), EEDetId::denseIndex(), ebToken_, DetId::Ecal, EcalBarrel, EcalEndcap, PhiSymCache::ecalLumiInfo, eeToken_, HLT_FULL_cff::endcapRecHits, relval_parameters_module::energy, PVValHelper::eta, ConfigCache::etCutsEB, ConfigCache::etCutsEE, eThresholdsEB_, eThresholdsEE_, relativeConstraints::geometry, geoToken_, edm::EventSetup::getData(), EcalRingCalibrationTools::getRingIndex(), kNRingsEB, EcalCondDBWriter_cfi::laser, laserDbToken_, misCalibRangeEB_, misCalibRangeEE_, misCalibStepsEB_, misCalibStepsEE_, nMisCalib_, nSumEtValues_, PhiSymCache::recHitCollEB, PhiSymCache::recHitCollEE, relativeConstraints::ring, edm::EventBase::time(), and edm::Timestamp::value().

Referenced by EcalPhiSymRecHitProducerLumi::accumulate(), and EcalPhiSymRecHitProducerRun::accumulate().

                                                                                 {
   uint64_t totHitsEB = 0;
   uint64_t totHitsEE = 0;
 
   //---get recHits collections
   auto barrelRecHits = event.get(ebToken_);
   auto endcapRecHits = event.get(eeToken_);
 
   //---get the laser corrections
   edm::Timestamp evtTimeStamp(event.time().value());
   auto const& laser = setup.getData(laserDbToken_);
 
   //---get the geometry
   auto const& geometry = setup.getData(geoToken_);
   auto barrelGeometry = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
   auto endcapGeometry = geometry.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
 
   //---EB---
   for (auto& recHit : barrelRecHits) {
     float energy = recHit.energy();
     EBDetId ebHit = EBDetId(recHit.id());
     int ring = calibRing_.getRingIndex(ebHit);
     //---if recHit energy is below thr even with the highest miscalib skip this recHit
     if (energy * misCalibRangeEB_[1] < eThresholdsEB_[ring])
       continue;
     float eta = barrelGeometry->getGeometry(ebHit)->getPosition().eta();
 
     //---compute et + miscalibration
     std::vector<float> etValues(nSumEtValues_, 0);
     //---one can do this in one for loop from -nMis to +nMis but in this way the
     //---program is faster
     //---NOTE: nMisCalib is half on the value set in the cfg python
     etValues[0] = recHit.energy() / cosh(eta);
     for (int iMis = -nMisCalib_; iMis < 0; ++iMis) {
       //--- 0 -> 0; -i -> [1...n/2]; +i -> [n/2+1...n]
       int index = iMis + nMisCalib_ + 1;
       etValues[index] = etValues[0] * (1 + misCalibStepsEB_[index]);
       //---set et to zero if out of range [e_thr, et_thr+1]
       if (etValues[index] * cosh(eta) < eThresholdsEB_[ring] || etValues[index] > configCache->etCutsEB[ring])
         etValues[index] = 0;
     }
     for (int iMis = 1; iMis <= nMisCalib_; ++iMis) {
       //--- 0 -> 0; -i -> [1...n/2]; +i -> [n/2+1...n]
       int index = iMis + nMisCalib_;
       etValues[index] = etValues[0] * (1 + misCalibStepsEB_[index]);
       //---set et to zero if out of range [e_thr, et_thr+1]
       if (etValues[index] * cosh(eta) < eThresholdsEB_[ring] || etValues[index] > configCache->etCutsEB[ring])
         etValues[index] = 0;
     }
     //---set et to zero if out of range [e_thr, et_thr+1]
     if (energy < eThresholdsEB_[ring] || etValues[0] > configCache->etCutsEB[ring])
       etValues[0] = 0;
     else
       ++totHitsEB;
     //---update the rechHit sumEt
     streamCache->recHitCollEB.at(ebHit.denseIndex())
         .addHit(etValues, laser.getLaserCorrection(recHit.id(), evtTimeStamp));
   }
 
   //---EE---
   for (auto& recHit : endcapRecHits) {
     EEDetId eeHit = EEDetId(recHit.id());
     int ring = calibRing_.getRingIndex(eeHit) - kNRingsEB;
     float energy = recHit.energy();
     //---if recHit energy is below thr even with the highest miscalib skip this recHit
     if (energy * misCalibRangeEE_[1] < eThresholdsEE_[ring])
       continue;
     float eta = endcapGeometry->getGeometry(eeHit)->getPosition().eta();
 
     //---compute et + miscalibration
     std::vector<float> etValues(nSumEtValues_, 0);
     //---one can do this in one for loop from -nMis to +nMis but in this way the
     //---program is faster
     //---NOTE: nMisCalib is half on the value set in the cfg python
     etValues[0] = recHit.energy() / cosh(eta);
     for (int iMis = -nMisCalib_; iMis < 0; ++iMis) {
       //--- 0 -> 0; -i -> [1...n/2]; +i -> [n/2+1...n]
       int index = iMis + nMisCalib_ + 1;
       etValues[index] = etValues[0] * (1 + misCalibStepsEE_[index]);
       //---set et to zero if out of range [e_thr, et_thr+1]
       if (etValues[index] * cosh(eta) < eThresholdsEE_[ring] || etValues[index] > configCache->etCutsEE[ring])
         etValues[index] = 0;
     }
     for (int iMis = 1; iMis <= nMisCalib_; ++iMis) {
       //--- 0 -> 0; -i -> [1...n/2]; +i -> [n/2+1...n]
       int index = iMis + nMisCalib_;
       etValues[index] = etValues[0] * (1 + misCalibStepsEE_[index]);
       //---set et to zero if out of range [e_thr, et_thr+1]
       if (etValues[index] * cosh(eta) < eThresholdsEE_[ring] || etValues[index] > configCache->etCutsEE[ring])
         etValues[index] = 0;
     }
     //---set et to zero if out of range [e_thr, et_thr+1]
     if (energy < eThresholdsEE_[ring] || etValues[0] > configCache->etCutsEE[ring])
       etValues[0] = 0;
     else
       ++totHitsEE;
     //---update the rechHit sumEt
     streamCache->recHitCollEE.at(eeHit.denseIndex())
         .addHit(etValues, laser.getLaserCorrection(recHit.id(), evtTimeStamp));
   }
 
   //---update the lumi info
   EcalPhiSymInfo thisEvent(totHitsEB, totHitsEE, 1, 0, 0, 0, 0);
   streamCache->ecalLumiInfo += thisEvent;
 }

void EcalPhiSymRecHitProducerBase::sumCache	(	PhiSymCache *	summaryc,
		PhiSymCache *	streamc
	)		const

Definition at line 348 of file EcalPhiSymRecHitProducers.cc.

References PhiSymCache::ecalLumiInfo, for(), mps_fire::i, PhiSymCache::recHitCollEB, and PhiSymCache::recHitCollEE.

Referenced by EcalPhiSymRecHitProducerLumi::streamEndLuminosityBlockSummary(), and EcalPhiSymRecHitProducerRun::streamEndRunSummary().

                                                                                              {
   //---The first argument is the summary cache that
   //   contains the lumi/run summary information.
   //   The stream partial sums are passed as second argument
   summaryc->ecalLumiInfo += streamc->ecalLumiInfo;
   for (unsigned int i = 0; i < summaryc->recHitCollEB.size(); ++i)
     summaryc->recHitCollEB[i] += streamc->recHitCollEB[i];
   for (unsigned int i = 0; i < summaryc->recHitCollEE.size(); ++i)
     summaryc->recHitCollEE[i] += streamc->recHitCollEE[i];
 }