#include <EcalClusterPUCleaningTools.h>

Public Member Functions
reco::SuperCluster	CleanedSuperCluster (float xi, const reco::SuperCluster &cluster, const edm::Event &ev, const edm::EventSetup &es)

	EcalClusterPUCleaningTools (edm::ConsumesCollector &cc, const edm::InputTag &redEBRecHits, const edm::InputTag &redEERecHits)

	~EcalClusterPUCleaningTools ()

Private Member Functions
void	getEBRecHits (const edm::Event &ev)

void	getEERecHits (const edm::Event &ev)

Private Attributes
const EcalRecHitCollection *	ebRecHits_

const EcalRecHitCollection *	eeRecHits_

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geometryToken_

const edm::EDGetTokenT< EcalRecHitCollection >	pEBRecHitsToken_

const edm::EDGetTokenT< EcalRecHitCollection >	pEERecHitsToken_

Detailed Description

tool to clean reco::Supercluster from effects of multiple interactions ( PU )

Author: G. Franzoni - UMN

Version: $Id:

Definition at line 24 of file EcalClusterPUCleaningTools.h.

Constructor & Destructor Documentation

◆ EcalClusterPUCleaningTools()

EcalClusterPUCleaningTools::EcalClusterPUCleaningTools	(	edm::ConsumesCollector &	cc,
		const edm::InputTag &	redEBRecHits,
		const edm::InputTag &	redEERecHits
	)

Definition at line 13 of file EcalClusterPUCleaningTools.cc.

     : pEBRecHitsToken_(cc.consumes<EcalRecHitCollection>(redEBRecHits)),
       pEERecHitsToken_(cc.consumes<EcalRecHitCollection>(redEERecHits)),
       geometryToken_(cc.esConsumes()) {}

◆ ~EcalClusterPUCleaningTools()

EcalClusterPUCleaningTools::~EcalClusterPUCleaningTools ( )

Definition at line 20 of file EcalClusterPUCleaningTools.cc.

20 {}

Member Function Documentation

◆ CleanedSuperCluster()

reco::SuperCluster EcalClusterPUCleaningTools::CleanedSuperCluster	(	float	xi,
		const reco::SuperCluster &	cluster,
		const edm::Event &	ev,
		const edm::EventSetup &	es
	)

Definition at line 34 of file EcalClusterPUCleaningTools.cc.

References cms::cuda::assert(), SuperClusterShapeAlgo::Calculate_Covariances(), reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), ebRecHits_, DetId::Ecal, EcalBarrel, EcalEndcap, eeRecHits_, HCALHighEnergyHPDFilter_cfi::energy, SuperClusterShapeAlgo::etaWidth(), photons_cff::etaWidth, makeMEIFBenchmarkPlots::ev, geometryToken_, edm::EventSetup::getData(), getEBRecHits(), getEERecHits(), SuperClusterShapeAlgo::phiWidth(), photons_cff::phiWidth, RecoTauValidation_cfi::posX, RecoTauValidation_cfi::posY, reco::SuperCluster::preshowerEnergy(), edm::PtrVector< T >::push_back(), reco::SuperCluster::seed(), fileCollector::seed, mathSSE::sqrt(), and protons_cff::xi.

                                                                                             {
   // get the geometry
   const auto &geometry = es.getData(geometryToken_);
 
   // get the RecHits
   getEBRecHits(ev);
   getEERecHits(ev);
 
   // seed basic cluster of initial SC: this will remain in the cleaned SC, by construction
   const reco::CaloClusterPtr &seed = scluster.seed();
 
   // this should be replaced by the 5x5 around the seed; a good approx of how E_seed is defined
   float seedBCEnergy = (scluster.seed())->energy();
   float eSeed = 0.35;  // standard eSeed in EB ; see CMS IN-2010/008
   int numBcRemoved = 0;
 
   double ClusterE = 0;  //Sum of cluster energies for supercluster.
   //Holders for position of this supercluster.
   double posX = 0;
   double posY = 0;
   double posZ = 0;
 
   reco::CaloClusterPtrVector thissc;
 
   // looping on basic clusters within the Supercluster
   for (reco::CaloCluster_iterator bcIt = scluster.clustersBegin(); bcIt != scluster.clustersEnd(); bcIt++) {
     // E_seed is an Et  selection on 5x1 dominoes (around which BC's are built), see CMS IN-2010/008
     // here such selection is implemented on the full BC around it
     if ((*bcIt)->energy() >
         sqrt(eSeed * eSeed + xi * xi * seedBCEnergy * seedBCEnergy / cosh((*bcIt)->eta()) / cosh((*bcIt)->eta()))) {
       ;
     }  // the sum only of the BC's that pass the Esee selection
     else {
       numBcRemoved++;
       continue;
     }  // count how many basic cluster get removed by the cleaning
 
     // if BC passes dynamic selection, include it in the 'cleaned' supercluster
     thissc.push_back((*bcIt));
     // cumulate energy and position of the cleaned supercluster
     ClusterE += (*bcIt)->energy();
     posX += (*bcIt)->energy() * (*bcIt)->position().X();
     posY += (*bcIt)->energy() * (*bcIt)->position().Y();
     posZ += (*bcIt)->energy() * (*bcIt)->position().Z();
 
   }  // loop on basic clusters of the original SC
 
   posX /= ClusterE;
   posY /= ClusterE;
   posZ /= ClusterE;
 
   // make temporary 'cleaned' supercluster in oder to compute the covariances
   double Epreshower = scluster.preshowerEnergy();
   double phiWidth = 0.;
   double etaWidth = 0.;
   reco::SuperCluster suCltmp(ClusterE, math::XYZPoint(posX, posY, posZ), seed, thissc, Epreshower, phiWidth, etaWidth);
 
   // construct cluster shape to compute ieta and iphi covariances of the SC
   const CaloSubdetectorGeometry *geometry_p = nullptr;
   if (seed->seed().det() == DetId::Ecal && seed->seed().subdetId() == EcalBarrel) {
     geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
     SuperClusterShapeAlgo SCShape(ebRecHits_, geometry_p);
     SCShape.Calculate_Covariances(suCltmp);
     phiWidth = SCShape.phiWidth();
     etaWidth = SCShape.etaWidth();
   } else if (seed->seed().det() == DetId::Ecal && seed->seed().subdetId() == EcalEndcap) {
     geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
     SuperClusterShapeAlgo SCShape(eeRecHits_, geometry_p);
     SCShape.Calculate_Covariances(suCltmp);
     phiWidth = SCShape.phiWidth();
     etaWidth = SCShape.etaWidth();
   } else {
     edm::LogError("SeedError")
         << "The seed crystal of this SC is neither in EB nor in EE. This is a problem. Bailing out";
     assert(-1);
   }
 
   // return the cleaned supercluster SCluster, with covariances updated
   return reco::SuperCluster(ClusterE, math::XYZPoint(posX, posY, posZ), seed, thissc, Epreshower, phiWidth, etaWidth);
 }

◆ getEBRecHits()

void EcalClusterPUCleaningTools::getEBRecHits ( const edm::Event & ev )

private

Definition at line 22 of file EcalClusterPUCleaningTools.cc.

References ebRecHits_, makeMEIFBenchmarkPlots::ev, pEBRecHitsToken_, and edm::Handle< T >::product().

Referenced by CleanedSuperCluster().

                                                                 {
   edm::Handle<EcalRecHitCollection> pEBRecHits;
   ev.getByToken(pEBRecHitsToken_, pEBRecHits);
   ebRecHits_ = pEBRecHits.product();
 }

◆ getEERecHits()

void EcalClusterPUCleaningTools::getEERecHits ( const edm::Event & ev )

private

Definition at line 28 of file EcalClusterPUCleaningTools.cc.

References eeRecHits_, makeMEIFBenchmarkPlots::ev, pEERecHitsToken_, and edm::Handle< T >::product().

Referenced by CleanedSuperCluster().

                                                                 {
   edm::Handle<EcalRecHitCollection> pEERecHits;
   ev.getByToken(pEERecHitsToken_, pEERecHits);
   eeRecHits_ = pEERecHits.product();
 }