EcalClusterPUCleaningTools Class Reference

#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	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		const edm::InputTag &	redEBRecHits,
		const edm::InputTag &	redEERecHits
	)

Definition at line 14 of file EcalClusterPUCleaningTools.cc.

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

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

EcalClusterPUCleaningTools::~EcalClusterPUCleaningTools

(

)

Definition at line 23 of file EcalClusterPUCleaningTools.cc.

{}

Member Function Documentation

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

Definition at line 43 of file EcalClusterPUCleaningTools.cc.

References 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(), SurveyInfoScenario_cff::seed, and mathSSE::sqrt().

                                                                                       {
  //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);
}

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

private

Definition at line 31 of file EcalClusterPUCleaningTools.cc.

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

Referenced by EcalClusterPUCleaningTools().

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

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

private

Definition at line 37 of file EcalClusterPUCleaningTools.cc.

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

Referenced by EcalClusterPUCleaningTools().

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

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

private

Definition at line 25 of file EcalClusterPUCleaningTools.cc.

References geometry_, edm::EventSetup::get(), and edm::ESHandle< T >::product().

Referenced by EcalClusterPUCleaningTools().

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

Member Data Documentation

const EcalRecHitCollection* EcalClusterPUCleaningTools::ebRecHits_

private

Definition at line 37 of file EcalClusterPUCleaningTools.h.

Referenced by CleanedSuperCluster(), and getEBRecHits().

const EcalRecHitCollection* EcalClusterPUCleaningTools::eeRecHits_

private

Definition at line 38 of file EcalClusterPUCleaningTools.h.

Referenced by CleanedSuperCluster(), and getEERecHits().

const CaloGeometry* EcalClusterPUCleaningTools::geometry_

private

Definition at line 36 of file EcalClusterPUCleaningTools.h.

Referenced by CleanedSuperCluster(), and getGeometry().