#include <SuperClusterShapeAlgo.h>

Public Member Functions
void	Calculate_Covariances (const reco::SuperCluster &passedCluster)

double	etaWidth ()

double	phiWidth ()

	SuperClusterShapeAlgo (const EcalRecHitCollection hits, const CaloSubdetectorGeometry geometry)

Private Attributes
double	etaWidth_

const CaloSubdetectorGeometry *	geometry_

double	phiWidth_

const EcalRecHitCollection *	recHits_

Detailed Description

Definition at line 27 of file SuperClusterShapeAlgo.h.

Constructor & Destructor Documentation

SuperClusterShapeAlgo::SuperClusterShapeAlgo	(	const EcalRecHitCollection *	hits,
		const CaloSubdetectorGeometry *	geometry
	)

Definition at line 14 of file SuperClusterShapeAlgo.cc.

15 : recHits_(hits), geometry_(geo) {}

SuperClusterShapeAlgo::geometry_

const CaloSubdetectorGeometry * geometry_

Definition: SuperClusterShapeAlgo.h:38

SuperClusterShapeAlgo::recHits_

const EcalRecHitCollection * recHits_

Definition: SuperClusterShapeAlgo.h:37

Member Function Documentation

void SuperClusterShapeAlgo::Calculate_Covariances ( const reco::SuperCluster & passedCluster )

Definition at line 17 of file SuperClusterShapeAlgo.cc.

References HLTTauDQMOffline_cfi::denominator, HLT_2018_cff::dEta, HLT_2018_cff::dPhi, edm::SortedCollection< T, SORT >::end(), reco::CaloCluster::energy(), PV3DBase< T, PVType, FrameType >::eta(), etaWidth_, edm::SortedCollection< T, SORT >::find(), geometry_, CaloSubdetectorGeometry::getGeometry(), reco::CaloCluster::hitsAndFractions(), PV3DBase< T, PVType, FrameType >::phi(), phiWidth_, Geom::pi(), reco::CaloCluster::position(), position, recHits_, mathSSE::sqrt(), and Geom::twoPi().

Referenced by EgammaSCEnergyCorrectionAlgo::applyCorrection(), HiEgammaSCEnergyCorrectionAlgo::applyCorrection(), and EcalClusterPUCleaningTools::CleanedSuperCluster().

                                                                                        {
   double numeratorEtaWidth = 0;
   double numeratorPhiWidth = 0;
 
   double scEnergy = passedCluster.energy();
   double denominator = scEnergy;
 
   double scEta = passedCluster.position().eta();
   double scPhi = passedCluster.position().phi();
 
   const std::vector<std::pair<DetId, float> >& detId = passedCluster.hitsAndFractions();
   // Loop over recHits associated with the given SuperCluster
   for (std::vector<std::pair<DetId, float> >::const_iterator hit = detId.begin(); hit != detId.end(); ++hit) {
     EcalRecHitCollection::const_iterator rHit = recHits_->find((*hit).first);
     //FIXME: THIS IS JUST A WORKAROUND A FIX SHOULD BE APPLIED
     if (rHit == recHits_->end()) {
       continue;
     }
     auto this_cell = geometry_->getGeometry(rHit->id());
     if (this_cell == nullptr) {
       //edm::LogInfo("SuperClusterShapeAlgo") << "pointer to the cell in Calculate_Covariances is NULL!";
       continue;
     }
     GlobalPoint position = this_cell->getPosition();
     //take into account energy fractions
     double energyHit = rHit->energy() * hit->second;
 
     //form differences
     double dPhi = position.phi() - scPhi;
     if (dPhi > +Geom::pi()) {
       dPhi = Geom::twoPi() - dPhi;
     }
     if (dPhi < -Geom::pi()) {
       dPhi = Geom::twoPi() + dPhi;
     }
 
     double dEta = position.eta() - scEta;
 
     if (energyHit > 0) {
       numeratorEtaWidth += energyHit * dEta * dEta;
       numeratorPhiWidth += energyHit * dPhi * dPhi;
     }
 
     etaWidth_ = sqrt(numeratorEtaWidth / denominator);
     phiWidth_ = sqrt(numeratorPhiWidth / denominator);
   }
 }