#include <EcalClusterPUCleaningTools.h>

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

	EcalClusterPUCleaningTools (const edm::Event &ev, const edm::EventSetup &es, const edm::InputTag &redEBRecHits, const edm::InputTag &redEERecHits)

	~EcalClusterPUCleaningTools ()

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

void	getEERecHits (const edm::Event &ev, const edm::InputTag &redEERecHits)

void	getGeometry (const edm::EventSetup &es)

Private Attributes
const EcalRecHitCollection *	ebRecHits_

const EcalRecHitCollection *	eeRecHits_

const CaloGeometry *	geometry_

Detailed Description

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

Author: G. Franzoni - UMN

Version: $Id:

Definition at line 22 of file EcalClusterPUCleaningTools.h.

Constructor & Destructor Documentation

◆ EcalClusterPUCleaningTools()

EcalClusterPUCleaningTools::EcalClusterPUCleaningTools	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		const edm::InputTag &	redEBRecHits,
		const edm::InputTag &	redEERecHits
	)

Definition at line 14 of file EcalClusterPUCleaningTools.cc.

                                                                                         {
   getGeometry(es);
   getEBRecHits(ev, redEBRecHits);
   getEERecHits(ev, redEERecHits);
 }

References ev, getEBRecHits(), getEERecHits(), and getGeometry().

◆ ~EcalClusterPUCleaningTools()

EcalClusterPUCleaningTools::~EcalClusterPUCleaningTools ( )

Definition at line 23 of file EcalClusterPUCleaningTools.cc.

23 {}

Member Function Documentation

◆ CleanedSuperCluster()

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

Definition at line 43 of file EcalClusterPUCleaningTools.cc.

                                                                                        {
   //std::cout << "\nEcalClusterPUCleaningTools::CleanedSuperCluster called, this will give you back a cleaned supercluster" << std::endl;
  
   // 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 {
     std::cout << "The seed crystal of this SC is neither in EB nor in EE. This is a problem. Bailing out " << std::endl;
     assert(-1);
   }
  
   // return the cleaned supercluster SCluster, with covariances updated
   return reco::SuperCluster(ClusterE, math::XYZPoint(posX, posY, posZ), seed, thissc, Epreshower, phiWidth, etaWidth);
 }

References cms::cuda::assert(), SuperClusterShapeAlgo::Calculate_Covariances(), reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), gather_cfg::cout, ebRecHits_, DetId::Ecal, EcalBarrel, EcalEndcap, eeRecHits_, HCALHighEnergyHPDFilter_cfi::energy, SuperClusterShapeAlgo::etaWidth(), geometry_, CaloGeometry::getSubdetectorGeometry(), SuperClusterShapeAlgo::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.

◆ getEBRecHits()

void EcalClusterPUCleaningTools::getEBRecHits	(	const edm::Event &	ev,
		const edm::InputTag &	redEBRecHits
	)

private

Definition at line 31 of file EcalClusterPUCleaningTools.cc.

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

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

Referenced by EcalClusterPUCleaningTools().

◆ getEERecHits()

void EcalClusterPUCleaningTools::getEERecHits	(	const edm::Event &	ev,
		const edm::InputTag &	redEERecHits
	)

private

Definition at line 37 of file EcalClusterPUCleaningTools.cc.

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

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

Referenced by EcalClusterPUCleaningTools().

◆ getGeometry()

void EcalClusterPUCleaningTools::getGeometry ( const edm::EventSetup & es )

private

Definition at line 25 of file EcalClusterPUCleaningTools.cc.

                                                                     {
   edm::ESHandle<CaloGeometry> pGeometry;
   es.get<CaloGeometryRecord>().get(pGeometry);
   geometry_ = pGeometry.product();
 }