PFECALSuperClusterAlgo Class Reference

\ Algorithm for box particle flow super clustering in the ECAL More...

#include <PFECALSuperClusterAlgo.h>

Classes
class	CalibratedPFCluster

Public Types
typedef std::shared_ptr< CalibratedPFCluster >	CalibratedClusterPtr
typedef std::vector< CalibratedClusterPtr >	CalibratedClusterPtrVector
enum	clustering_type { kBOX =1, kMustache =2 }
enum	energy_weight { kRaw, kCalibratedNoPS, kCalibratedTotal }

Public Member Functions
std::unique_ptr< reco::SuperClusterCollection > &	getEBOutputSCCollection ()
std::unique_ptr< reco::SuperClusterCollection > &	getEEOutputSCCollection ()
void	loadAndSortPFClusters (const edm::Event &evt)
	PFECALSuperClusterAlgo ()
	constructor More...
void	run ()
void	setClusteringType (clustering_type thetype)
void	setCrackCorrections (bool applyCrackCorrections)
void	setDropUnseedable (const bool d)
void	setEnergyWeighting (energy_weight thetype)
void	setEtawidthSuperClusterBarrel (double etawidth)
void	setEtawidthSuperClusterEndcap (double etawidth)
void	setIsOOTCollection (bool isOOTCollection)
void	setMajorityFraction (const double f)
void	setPFClusterCalibration (const std::shared_ptr< PFEnergyCalibration > &)
void	setPhiwidthSuperClusterBarrel (double phiwidth)
void	setPhiwidthSuperClusterEndcap (double phiwidth)
void	setSatelliteMerging (const bool doit)
void	setSatelliteThreshold (const double t)
void	setThreshPFClusterBarrel (double thresh)
void	setThreshPFClusterEndcap (double thresh)
void	setThreshPFClusterSeedBarrel (double thresh)
void	setThreshPFClusterSeedEndcap (double thresh)
void	setThreshSuperClusterEt (double thresh)
void	setTokens (const edm::ParameterSet &, edm::ConsumesCollector &&)
void	setUseDynamicDPhi (bool useit)
void	setUseETForSeeding (bool useET)
void	setUseRegression (bool useRegression)
void	setVerbosityLevel (bool verbose)
void	update (const edm::EventSetup &)

Private Member Functions
void	buildAllSuperClusters (CalibratedClusterPtrVector &, double seedthresh)
void	buildSuperCluster (CalibratedClusterPtr &, CalibratedClusterPtrVector &)

Private Attributes
CalibratedClusterPtrVector	_clustersEB
CalibratedClusterPtrVector	_clustersEE
clustering_type	_clustype
energy_weight	_eweight
std::shared_ptr< PFEnergyCalibration >	_pfEnergyCalibration
bool	applyCrackCorrections_
const EcalRecHitCollection *	barrelRecHits_
const reco::BeamSpot *	beamSpot_
const ESChannelStatus *	channelStatus_
bool	doSatelliteClusterMerge_
bool	dropUnseedable_
const reco::PFCluster::EEtoPSAssociation *	EEtoPS_
const EcalRecHitCollection *	endcapRecHits_
double	etawidthSuperClusterBarrel_
double	etawidthSuperClusterEndcap_
double	fractionForMajority_
edm::EDGetTokenT< EcalRecHitCollection >	inputTagBarrelRecHits_
edm::EDGetTokenT< reco::BeamSpot >	inputTagBeamSpot_
edm::EDGetTokenT< EcalRecHitCollection >	inputTagEndcapRecHits_
edm::EDGetTokenT< edm::View< reco::PFCluster > >	inputTagPFClusters_
edm::EDGetTokenT< reco::PFCluster::EEtoPSAssociation >	inputTagPFClustersES_
bool	isOOTCollection_
double	phiwidthSuperClusterBarrel_
double	phiwidthSuperClusterEndcap_
std::unique_ptr< SCEnergyCorrectorSemiParm >	regr_
double	satelliteThreshold_
std::unique_ptr< reco::SuperClusterCollection >	superClustersEB_
std::unique_ptr< reco::SuperClusterCollection >	superClustersEE_
bool	threshIsET_
double	threshPFClusterBarrel_
double	threshPFClusterEndcap_
double	threshPFClusterSeedBarrel_
double	threshPFClusterSeedEndcap_
double	threshSuperClusterEt_
bool	useDynamicDPhi_
bool	useRegression_
bool	verbose_

Detailed Description

\ Algorithm for box particle flow super clustering in the ECAL

Original Author: Nicolas Chanon Additional Authors (Mustache): Y. Gershtein, R. Patel, L. Gray

Date

July 2012

Definition at line 43 of file PFECALSuperClusterAlgo.h.

Member Typedef Documentation

typedef std::shared_ptr<CalibratedPFCluster> PFECALSuperClusterAlgo::CalibratedClusterPtr

Definition at line 63 of file PFECALSuperClusterAlgo.h.

typedef std::vector<CalibratedClusterPtr> PFECALSuperClusterAlgo::CalibratedClusterPtrVector

Definition at line 64 of file PFECALSuperClusterAlgo.h.

Member Enumeration Documentation

enum PFECALSuperClusterAlgo::clustering_type

Enumerator
kBOX
kMustache

Definition at line 45 of file PFECALSuperClusterAlgo.h.

{kBOX=1, kMustache=2};

enum PFECALSuperClusterAlgo::energy_weight

Enumerator
kRaw
kCalibratedNoPS
kCalibratedTotal

Definition at line 46 of file PFECALSuperClusterAlgo.h.

{kRaw, kCalibratedNoPS, kCalibratedTotal};

Constructor & Destructor Documentation

PFECALSuperClusterAlgo::PFECALSuperClusterAlgo

(

)

constructor

Definition at line 122 of file PFECALSuperClusterAlgo.cc.

References setPFClusterCalibration().

: beamSpot_(nullptr) { }

Member Function Documentation

void PFECALSuperClusterAlgo::buildAllSuperClusters ( CalibratedClusterPtrVector &  , double  seedthresh  )

private

Definition at line 264 of file PFECALSuperClusterAlgo.cc.

References buildSuperCluster(), fastPrimaryVertexProducer_cfi::clusters, and threshIsET_.

Referenced by run().

                                 {
  auto seedable = std::bind(isSeed, _1, seedthresh, threshIsET_);
  // make sure only seeds appear at the front of the list of clusters
  std::stable_partition(clusters.begin(),clusters.end(),seedable);
  // in each iteration we are working on a list that is already sorted
  // in the cluster energy and remains so through each iteration
  // NB: since clusters is sorted in loadClusters any_of has O(1)
  //     timing until you run out of seeds!  
  while( std::any_of(clusters.cbegin(), clusters.cend(), seedable) ) {    
    buildSuperCluster(clusters.front(),clusters);
  }
}

void PFECALSuperClusterAlgo::buildSuperCluster ( CalibratedClusterPtr &  , CalibratedClusterPtrVector &    )

private

Definition at line 278 of file PFECALSuperClusterAlgo.cc.

References _clustype, _eweight, _pfEnergyCalibration, reco::SuperCluster::addCluster(), reco::CaloCluster::addHitAndFraction(), reco::SuperCluster::addPreshowerCluster(), applyCrackCorrections_, barrelRecHits_, beamSpot_, channelStatus_, TauDecayModes::dec, doSatelliteClusterMerge_, dropUnseedable_, PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, EEtoPS_, endcapRecHits_, reco::PFCluster::energy(), reco::CaloCluster::energy(), etawidthSuperClusterBarrel_, etawidthSuperClusterEndcap_, Exception, ESCondObjectContainer< T >::find(), edm::SortedCollection< T, SORT >::find(), fractionForMajority_, edm::Ptr< T >::get(), ESCondObjectContainer< T >::getMap(), PFLayer::HGCAL, mps_fire::i, isOOTCollection_, kCalibratedNoPS, kCalibratedTotal, EcalRecHit::kOutOfTime, kRaw, reco::PFCluster::layer(), PFClusterWidthAlgo::pflowEtaWidth(), PFClusterWidthAlgo::pflowPhiWidth(), phiwidthSuperClusterBarrel_, phiwidthSuperClusterEndcap_, reco::BeamSpot::position(), reco::CaloCluster::position(), RecoTauValidation_cfi::posX, RecoTauValidation_cfi::posY, reco::SuperCluster::preshowerClustersBegin(), reco::SuperCluster::preshowerClustersEnd(), PFLayer::PS1, PFLayer::PS2, ptFast(), reco::SuperCluster::rawEnergy(), reco::PFCluster::recHitFractions(), regr_, satelliteThreshold_, reco::SuperCluster::seed(), reco::CaloCluster::setCorrectedEnergy(), reco::SuperCluster::setEtaWidth(), reco::SuperCluster::setPhiWidth(), reco::SuperCluster::setPreshowerEnergy(), reco::SuperCluster::setPreshowerEnergyPlane1(), reco::SuperCluster::setPreshowerEnergyPlane2(), reco::SuperCluster::setSeed(), sortByKey(), superClustersEB_, superClustersEE_, threshSuperClusterEt_, useDynamicDPhi_, useRegression_, and verbose_.

Referenced by buildAllSuperClusters().

{
  double etawidthSuperCluster = 0.0;
  double phiwidthSuperCluster = 0.0;
  bool isEE = false;
  switch( seed->the_ptr()->layer() )
  {
  case PFLayer::ECAL_BARREL:
    phiwidthSuperCluster = phiwidthSuperClusterBarrel_;
    etawidthSuperCluster = etawidthSuperClusterBarrel_;
    edm::LogInfo("PFClustering") << "Building SC number "  
                 << superClustersEB_->size() + 1
                 << " in the ECAL barrel!";
    break;
  case PFLayer::HGCAL:  
  case PFLayer::ECAL_ENDCAP:  
  
    phiwidthSuperCluster = phiwidthSuperClusterEndcap_; 
    etawidthSuperCluster = etawidthSuperClusterEndcap_;
    edm::LogInfo("PFClustering") << "Building SC number "  
                 << superClustersEE_->size() + 1
                 << " in the ECAL endcap!" << std::endl;
    isEE = true;
    break;
  default:
    break;
  }
  auto isClusteredWithSeed = std::bind(isClustered, _1, seed,_clustype,useDynamicDPhi_, etawidthSuperCluster, phiwidthSuperCluster);
  auto matchesSeedByRecHit = std::bind(isLinkedByRecHit, _1, seed,satelliteThreshold_, fractionForMajority_,0.1,0.2);
  
  // this function shuffles the list of clusters into a list
  // where all clustered sub-clusters are at the front 
  // and returns a pointer to the first unclustered cluster.
  // The relative ordering of clusters is preserved 
  // (i.e. both resulting sub-lists are sorted by energy).
  auto not_clustered = std::stable_partition(clusters.begin(),clusters.end(),
                         isClusteredWithSeed);
  // satellite cluster merging
  // it was found that large clusters can split!
  if( doSatelliteClusterMerge_ ) {    
    not_clustered = std::stable_partition(not_clustered,clusters.end(),
                      matchesSeedByRecHit);
  }

  if(verbose_) {
    edm::LogInfo("PFClustering") << "Dumping cluster detail";
    edm::LogVerbatim("PFClustering")
      << "\tPassed seed: e = " << seed->energy_nocalib() 
      << " eta = " << seed->eta() << " phi = " << seed->phi() 
      << std::endl;  
    for( auto clus = clusters.cbegin(); clus != not_clustered; ++clus ) {
      edm::LogVerbatim("PFClustering") 
    << "\t\tClustered cluster: e = " << (*clus)->energy_nocalib() 
    << " eta = " << (*clus)->eta() << " phi = " << (*clus)->phi() 
    << std::endl;
    }
    for( auto clus = not_clustered; clus != clusters.end(); ++clus ) {
      edm::LogVerbatim("PFClustering") 
    << "\tNon-Clustered cluster: e = " << (*clus)->energy_nocalib() 
    << " eta = " << (*clus)->eta() << " phi = " << (*clus)->phi() 
    << std::endl;
    }    
  }

  if (not_clustered == clusters.begin()) {
    if(dropUnseedable_){
      clusters.erase(clusters.begin());
      return;
    }
    else {
      throw cms::Exception("PFECALSuperClusterAlgo::buildSuperCluster")
        << "Cluster is not seedable!" << std::endl 
        << "\tNon-Clustered cluster: e = " << (*not_clustered)->energy_nocalib()
        << " eta = " << (*not_clustered)->eta() << " phi = " << (*not_clustered)->phi()
        << std::endl;
    }
  }

  // move the clustered clusters out of available cluster list
  // and into a temporary vector for building the SC  
  CalibratedClusterPtrVector clustered(clusters.begin(),not_clustered);
  clusters.erase(clusters.begin(),not_clustered);    
  // need the vector of raw pointers for a PF width class
  std::vector<const reco::PFCluster*> bare_ptrs;
  // calculate necessary parameters and build the SC
  double posX(0), posY(0), posZ(0),
    corrSCEnergy(0), corrPS1Energy(0), corrPS2Energy(0), 
    ePS1(0), ePS2(0), energyweight(0), energyweighttot(0);
  std::vector<double> ps1_energies, ps2_energies;
  int condP1(1), condP2(1);
  for( auto& clus : clustered ) {    
    ePS1 = ePS2 = 0;
    energyweight = clus->energy_nocalib();
    bare_ptrs.push_back(clus->the_ptr().get());
    // update EE calibrated super cluster energies
    if( isEE ) {
      ePS1 = ePS2 = 0;
      condP1 = condP2 = 1;
      ps1_energies.clear();
      ps2_energies.clear();
      auto ee_key_val = 
    std::make_pair(clus->the_ptr().key(),edm::Ptr<reco::PFCluster>());
      const auto clustops = std::equal_range(EEtoPS_->begin(),
                         EEtoPS_->end(),
                         ee_key_val,
                         sortByKey);      
      for( auto i_ps = clustops.first; i_ps != clustops.second; ++i_ps) {
    edm::Ptr<reco::PFCluster> psclus(i_ps->second);

    auto const& recH_Frac = psclus->recHitFractions();

    switch( psclus->layer() ) {
    case PFLayer::PS1:
      ps1_energies.push_back(psclus->energy());
      for (auto const& recH : recH_Frac){
            ESDetId strip1 = recH.recHitRef()->detId();
        if(strip1 != ESDetId(0)){
          ESChannelStatusMap::const_iterator status_p1 = channelStatus_->getMap().find(strip1);
          // getStatusCode() == 1 => dead channel
          //apply correction if all recHits in dead region
          if(status_p1->getStatusCode() == 0) condP1 = 0; //active  
        }
      }
      break;
    case PFLayer::PS2:
      ps2_energies.push_back(psclus->energy());
      for (auto const& recH : recH_Frac){
            ESDetId strip2 = recH.recHitRef()->detId();
        if(strip2 != ESDetId(0)) {
          ESChannelStatusMap::const_iterator status_p2 = channelStatus_->getMap().find(strip2);
          if(status_p2->getStatusCode() == 0) condP2 = 0;
        }
      }   
      break;
    default:
      break;
    }
      }    
      if(condP1 == 1) ePS1 = -1.;
      if(condP2 == 1) ePS2 = -1.;  
      _pfEnergyCalibration->energyEm(*(clus->the_ptr()),
                     ps1_energies,ps2_energies,
                     ePS1,ePS2,
                     applyCrackCorrections_);
    }

    if(ePS1 == -1.) ePS1 = 0;    
    if(ePS2 == -1.) ePS2 = 0;    

    switch( _eweight ) {
    case kRaw: // energyweight is initialized to raw cluster energy
      break;
    case kCalibratedNoPS:
      energyweight = clus->energy() - ePS1 - ePS2;
      break;
    case kCalibratedTotal:
      energyweight = clus->energy();
      break;
    default:
      break;
    }    
    const math::XYZPoint& cluspos = clus->the_ptr()->position();
    posX += energyweight * cluspos.X();
    posY += energyweight * cluspos.Y();
    posZ += energyweight * cluspos.Z();    

    energyweighttot += energyweight;
    corrSCEnergy    += clus->energy();   
    corrPS1Energy   += ePS1;
    corrPS2Energy   += ePS2;
  }
  posX /= energyweighttot;
  posY /= energyweighttot;
  posZ /= energyweighttot;    
  
  // now build the supercluster
  reco::SuperCluster new_sc(corrSCEnergy,math::XYZPoint(posX,posY,posZ));   
  new_sc.setCorrectedEnergy(corrSCEnergy);
  new_sc.setSeed(clustered.front()->the_ptr());
  new_sc.setPreshowerEnergy(corrPS1Energy+corrPS2Energy);
  new_sc.setPreshowerEnergyPlane1(corrPS1Energy);
  new_sc.setPreshowerEnergyPlane2(corrPS2Energy);
  for( const auto& clus : clustered ) {
    new_sc.addCluster(clus->the_ptr());
    
    auto& hits_and_fractions = clus->the_ptr()->hitsAndFractions();
    for( auto& hit_and_fraction : hits_and_fractions ) {
      new_sc.addHitAndFraction(hit_and_fraction.first,hit_and_fraction.second);
    }
    if( isEE ) {
      auto ee_key_val = 
    std::make_pair(clus->the_ptr().key(),edm::Ptr<reco::PFCluster>());
      const auto clustops = std::equal_range(EEtoPS_->begin(),
                         EEtoPS_->end(),
                         ee_key_val,
                         sortByKey);
      // EE rechits should be uniquely matched to sets of pre-shower
      // clusters at this point, so we throw an exception if otherwise
      // now wrapped in EDM debug flags
      for( auto i_ps = clustops.first; i_ps != clustops.second; ++i_ps) {
    edm::Ptr<reco::PFCluster> psclus(i_ps->second);
#ifdef EDM_ML_DEBUG
    auto found_pscluster = std::find_if(new_sc.preshowerClustersBegin(),
                     new_sc.preshowerClustersEnd(),
                     [&i_ps](const auto& i) { return i.key() == i_ps->first;});
    if( found_pscluster == new_sc.preshowerClustersEnd() ) {
#endif    
      new_sc.addPreshowerCluster(psclus);         
#ifdef EDM_ML_DEBUG
    } else {
      throw cms::Exception("PFECALSuperClusterAlgo::buildSuperCluster")
        << "Found a PS cluster matched to more than one EE cluster!" 
        << std::endl << std::hex << psclus.get() << " == " 
        << found_pscluster->get() << std::dec << std::endl;
    }
#endif
      }
    }
  }  
  
  // calculate linearly weighted cluster widths
  PFClusterWidthAlgo pfwidth(bare_ptrs);
  new_sc.setEtaWidth(pfwidth.pflowEtaWidth());
  new_sc.setPhiWidth(pfwidth.pflowPhiWidth());
  
  // cache the value of the raw energy  
  new_sc.rawEnergy();

  //apply regression energy corrections
  if( useRegression_ ) {  
    regr_->modifyObject(new_sc);
  }
  
  // save the super cluster to the appropriate list (if it passes the final
  // Et threshold)
  //Note that Et is computed here with respect to the beamspot position
  //in order to be consistent with the cut applied in the
  //ElectronSeedProducer
  double scEtBS = 
    ptFast(new_sc.energy(),new_sc.position(),beamSpot_->position());

  if ( scEtBS > threshSuperClusterEt_ ) {
    switch( seed->the_ptr()->layer() ) {
    case PFLayer::ECAL_BARREL:
      if(isOOTCollection_) {
    DetId seedId = new_sc.seed()->seed();
    EcalRecHitCollection::const_iterator seedRecHit = barrelRecHits_->find(seedId);
    if (!seedRecHit->checkFlag(EcalRecHit::kOutOfTime)) break;
      }
      superClustersEB_->push_back(new_sc);
      break;
    case PFLayer::HGCAL:
    case PFLayer::ECAL_ENDCAP:    
      if(isOOTCollection_) {
    DetId seedId = new_sc.seed()->seed();
    EcalRecHitCollection::const_iterator seedRecHit = endcapRecHits_->find(seedId);
    if (!seedRecHit->checkFlag(EcalRecHit::kOutOfTime)) break;
      }
      superClustersEE_->push_back(new_sc);    
      break;
    default:
      break;
    }
  }
}

std::unique_ptr<reco::SuperClusterCollection>& PFECALSuperClusterAlgo::getEBOutputSCCollection ( )

inline

Definition at line 110 of file PFECALSuperClusterAlgo.h.

References superClustersEB_.

{ return superClustersEB_; }

std::unique_ptr<reco::SuperClusterCollection>& PFECALSuperClusterAlgo::getEEOutputSCCollection ( )

inline

Definition at line 112 of file PFECALSuperClusterAlgo.h.

References loadAndSortPFClusters(), run(), and superClustersEE_.

{ return superClustersEE_; }

void PFECALSuperClusterAlgo::loadAndSortPFClusters

(

const edm::Event &

evt

)

Definition at line 172 of file PFECALSuperClusterAlgo.cc.

References _clustersEB, _clustersEE, barrelRecHits_, beamSpot_, fastPrimaryVertexProducer_cfi::clusters, PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, EEtoPS_, endcapRecHits_, Exception, edm::Event::getByToken(), PFLayer::HGCAL, mps_fire::i, inputTagBarrelRecHits_, inputTagBeamSpot_, inputTagEndcapRecHits_, inputTagPFClusters_, inputTagPFClustersES_, isOOTCollection_, edm::HandleBase::isValid(), LogDebug, edm::Handle< T >::product(), regr_, jetUpdater_cfi::sort, superClustersEB_, superClustersEE_, threshPFClusterBarrel_, threshPFClusterEndcap_, and useRegression_.

Referenced by getEEOutputSCCollection(), and update().

                                              { 
  
  //load input collections
  //Load the pfcluster collections
  edm::Handle<edm::View<reco::PFCluster> > pfclustersHandle;
  iEvent.getByToken( inputTagPFClusters_, pfclustersHandle );  

  edm::Handle<reco::PFCluster::EEtoPSAssociation > psAssociationHandle;
  iEvent.getByToken( inputTagPFClustersES_,  psAssociationHandle);  
  
  const PFClusterView& clusters = *pfclustersHandle.product();
  const reco::PFCluster::EEtoPSAssociation& psclusters = *psAssociationHandle.product();
  
  //load BeamSpot
  edm::Handle<reco::BeamSpot> bsHandle;
  iEvent.getByToken( inputTagBeamSpot_, bsHandle);
  beamSpot_ = bsHandle.product();
  
  //initialize regression for this event
  if (useRegression_) {
    regr_->setEvent(iEvent);
  }  
  
  // reset the system for running
  superClustersEB_ = std::make_unique<reco::SuperClusterCollection>();
  _clustersEB.clear();
  superClustersEE_ = std::make_unique<reco::SuperClusterCollection>();
  _clustersEE.clear();  
  EEtoPS_ = &psclusters;

  //Select PF clusters available for the clustering
  for ( size_t i = 0; i < clusters.size(); ++i ){
    auto cluster = clusters.ptrAt(i);
    LogDebug("PFClustering")
      << "Loading PFCluster i="<<cluster.key()
      <<" energy="<<cluster->energy()<<std::endl;
        
    // protection for sim clusters
    if( cluster->caloID().detectors() == 0 && 
    cluster->hitsAndFractions().empty() ) continue;

    CalibratedClusterPtr calib_cluster(new CalibratedPFCluster(cluster));
    switch( cluster->layer() ) {
    case PFLayer::ECAL_BARREL:
      if( calib_cluster->energy() > threshPFClusterBarrel_ ) {
    _clustersEB.push_back(calib_cluster);   
      }
      break;
    case PFLayer::HGCAL:
    case PFLayer::ECAL_ENDCAP:   
      if( calib_cluster->energy() > threshPFClusterEndcap_ ) {
    _clustersEE.push_back(calib_cluster);
      }
      break;      
    default:
      break;
    }
  }
  // sort full cluster collections by their calibrated energy
  // this will put all the seeds first by construction
  std::sort(_clustersEB.begin(), _clustersEB.end(), greaterByEt);
  std::sort(_clustersEE.begin(), _clustersEE.end(), greaterByEt);

  // set recHit collections for OOT photons
  if (isOOTCollection_)
  {
    edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
    iEvent.getByToken(inputTagBarrelRecHits_, barrelRecHitsHandle);
    if (!barrelRecHitsHandle.isValid()) {
      throw cms::Exception("PFECALSuperClusterAlgo") 
    << "If you use OOT photons, need to specify proper barrel rec hit collection";
    }
    barrelRecHits_ = barrelRecHitsHandle.product();

    edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
    iEvent.getByToken(inputTagEndcapRecHits_, endcapRecHitsHandle);
    if (!endcapRecHitsHandle.isValid()) {
      throw cms::Exception("PFECALSuperClusterAlgo") 
    << "If you use OOT photons, need to specify proper endcap rec hit collection";
    }
    endcapRecHits_ = endcapRecHitsHandle.product();
  }
}

void PFECALSuperClusterAlgo::run

(

)

Definition at line 256 of file PFECALSuperClusterAlgo.cc.

References _clustersEB, _clustersEE, buildAllSuperClusters(), threshPFClusterSeedBarrel_, and threshPFClusterSeedEndcap_.

Referenced by getEEOutputSCCollection().

                                 {  
  // clusterize the EB
  buildAllSuperClusters(_clustersEB, threshPFClusterSeedBarrel_);
  // clusterize the EE
  buildAllSuperClusters(_clustersEE, threshPFClusterSeedEndcap_);
}

void PFECALSuperClusterAlgo::setClusteringType ( clustering_type  thetype )

inline

Definition at line 72 of file PFECALSuperClusterAlgo.h.

References _clustype.

{ _clustype = thetype; } 

void PFECALSuperClusterAlgo::setCrackCorrections ( bool  applyCrackCorrections )

inline

Definition at line 103 of file PFECALSuperClusterAlgo.h.

References particleFlowSuperClusterECAL_cfi::applyCrackCorrections, applyCrackCorrections_, setTokens(), and update().

{ applyCrackCorrections_ = applyCrackCorrections;}

void PFECALSuperClusterAlgo::setDropUnseedable ( const bool  d )

inline

Definition at line 99 of file PFECALSuperClusterAlgo.h.

References edmIntegrityCheck::d, and dropUnseedable_.

{ dropUnseedable_ = d; }

void PFECALSuperClusterAlgo::setEnergyWeighting ( energy_weight  thetype )

inline

Definition at line 74 of file PFECALSuperClusterAlgo.h.

References _eweight.

{ _eweight = thetype; } 

void PFECALSuperClusterAlgo::setEtawidthSuperClusterBarrel ( double  etawidth )

inline

Definition at line 90 of file PFECALSuperClusterAlgo.h.

References etawidthSuperClusterBarrel_.

{ etawidthSuperClusterBarrel_ = etawidth;}

void PFECALSuperClusterAlgo::setEtawidthSuperClusterEndcap ( double  etawidth )

inline

Definition at line 92 of file PFECALSuperClusterAlgo.h.

References etawidthSuperClusterEndcap_, and setPFClusterCalibration().

{ etawidthSuperClusterEndcap_ = etawidth;}

void PFECALSuperClusterAlgo::setIsOOTCollection ( bool  isOOTCollection )

inline

Definition at line 101 of file PFECALSuperClusterAlgo.h.

References particleFlowSuperClusterOOTECAL_cff::isOOTCollection, and isOOTCollection_.

{ isOOTCollection_ = isOOTCollection; }

void PFECALSuperClusterAlgo::setMajorityFraction ( const double  f )

inline

Definition at line 98 of file PFECALSuperClusterAlgo.h.

References f, and fractionForMajority_.

{ fractionForMajority_ = f; }

void PFECALSuperClusterAlgo::setPFClusterCalibration

(

const std::shared_ptr< PFEnergyCalibration > &

calib

)

Definition at line 125 of file PFECALSuperClusterAlgo.cc.

References _pfEnergyCalibration, and calib.

Referenced by PFECALSuperClusterAlgo(), and setEtawidthSuperClusterEndcap().

                                                                         {
  _pfEnergyCalibration = calib;
}

void PFECALSuperClusterAlgo::setPhiwidthSuperClusterBarrel ( double  phiwidth )

inline

Definition at line 89 of file PFECALSuperClusterAlgo.h.

References phiwidthSuperClusterBarrel_.

{ phiwidthSuperClusterBarrel_ = phiwidth;}

void PFECALSuperClusterAlgo::setPhiwidthSuperClusterEndcap ( double  phiwidth )

inline

Definition at line 91 of file PFECALSuperClusterAlgo.h.

References phiwidthSuperClusterEndcap_.

{ phiwidthSuperClusterEndcap_ = phiwidth;}

void PFECALSuperClusterAlgo::setSatelliteMerging ( const bool  doit )

inline

Definition at line 96 of file PFECALSuperClusterAlgo.h.

References doSatelliteClusterMerge_.

{ doSatelliteClusterMerge_ = doit; }

void PFECALSuperClusterAlgo::setSatelliteThreshold ( const double  t )

inline

Definition at line 97 of file PFECALSuperClusterAlgo.h.

References satelliteThreshold_, and protons_cff::t.

{ satelliteThreshold_ = t; }

void PFECALSuperClusterAlgo::setThreshPFClusterBarrel ( double  thresh )

inline

Definition at line 85 of file PFECALSuperClusterAlgo.h.

References scrapingFilter_cfi::thresh, and threshPFClusterBarrel_.

{ threshPFClusterBarrel_ = thresh;}

void PFECALSuperClusterAlgo::setThreshPFClusterEndcap ( double  thresh )

inline

Definition at line 87 of file PFECALSuperClusterAlgo.h.

References scrapingFilter_cfi::thresh, and threshPFClusterEndcap_.

{ threshPFClusterEndcap_ = thresh;}

void PFECALSuperClusterAlgo::setThreshPFClusterSeedBarrel ( double  thresh )

inline

Definition at line 84 of file PFECALSuperClusterAlgo.h.

References scrapingFilter_cfi::thresh, and threshPFClusterSeedBarrel_.

{ threshPFClusterSeedBarrel_ = thresh;}

void PFECALSuperClusterAlgo::setThreshPFClusterSeedEndcap ( double  thresh )

inline

Definition at line 86 of file PFECALSuperClusterAlgo.h.

References scrapingFilter_cfi::thresh, and threshPFClusterSeedEndcap_.

{ threshPFClusterSeedEndcap_ = thresh;}

void PFECALSuperClusterAlgo::setThreshSuperClusterEt ( double  thresh )

inline

Definition at line 82 of file PFECALSuperClusterAlgo.h.

References scrapingFilter_cfi::thresh, and threshSuperClusterEt_.

{ threshSuperClusterEt_ = thresh; }

void PFECALSuperClusterAlgo::setTokens	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	cc
	)

Definition at line 129 of file PFECALSuperClusterAlgo.cc.

References edm::ParameterSet::getParameter(), inputTagBarrelRecHits_, inputTagBeamSpot_, inputTagEndcapRecHits_, inputTagPFClusters_, inputTagPFClustersES_, isOOTCollection_, regr_, and useRegression_.

Referenced by setCrackCorrections().

                                                                                                {
  
  inputTagPFClusters_ = 
    cc.consumes<edm::View<reco::PFCluster> >(iConfig.getParameter<edm::InputTag>("PFClusters"));
  inputTagPFClustersES_ = 
    cc.consumes<reco::PFCluster::EEtoPSAssociation>(iConfig.getParameter<edm::InputTag>("ESAssociation"));  
  inputTagBeamSpot_ = 
    cc.consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("BeamSpot"));    
    
  if (useRegression_) {
    const edm::ParameterSet &regconf = iConfig.getParameter<edm::ParameterSet>("regressionConfig");
    
    regr_.reset(new SCEnergyCorrectorSemiParm());
    regr_->setTokens(regconf, cc);  
  }

  if (isOOTCollection_) { // OOT photons only
    inputTagBarrelRecHits_ =
      cc.consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("barrelRecHits"));
    inputTagEndcapRecHits_ =
      cc.consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("endcapRecHits"));
  }  
}

void PFECALSuperClusterAlgo::setUseDynamicDPhi ( bool  useit )

inline

Definition at line 78 of file PFECALSuperClusterAlgo.h.

References useDynamicDPhi_.

{ useDynamicDPhi_ = useit; } 

void PFECALSuperClusterAlgo::setUseETForSeeding ( bool  useET )

inline

Definition at line 76 of file PFECALSuperClusterAlgo.h.

References threshIsET_.

{ threshIsET_ = useET; } 

void PFECALSuperClusterAlgo::setUseRegression ( bool  useRegression )

inline

Definition at line 80 of file PFECALSuperClusterAlgo.h.

References particleFlowSuperClusterECAL_cfi::useRegression, and useRegression_.

{ useRegression_ = useRegression; }

void PFECALSuperClusterAlgo::setVerbosityLevel ( bool  verbose )

inline

Definition at line 70 of file PFECALSuperClusterAlgo.h.

References verbose, and verbose_.

{ verbose_ = verbose;}

void PFECALSuperClusterAlgo::update

(

const edm::EventSetup &

setup

)

Definition at line 153 of file PFECALSuperClusterAlgo.cc.

References _pfEnergyCalibration, channelStatus_, edm::EventSetup::get(), loadAndSortPFClusters(), edm::ESHandle< T >::product(), regr_, and useRegression_.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), setCrackCorrections(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                              {
 
  if (useRegression_) {
    regr_->setEventSetup(setup);
  }
  
  edm::ESHandle<ESEEIntercalibConstants> esEEInterCalibHandle_;
  setup.get<ESEEIntercalibConstantsRcd>().get(esEEInterCalibHandle_);
  _pfEnergyCalibration->initAlphaGamma_ESplanes_fromDB(esEEInterCalibHandle_.product());

  edm::ESHandle<ESChannelStatus> esChannelStatusHandle_;
  setup.get<ESChannelStatusRcd>().get(esChannelStatusHandle_);
  channelStatus_ = esChannelStatusHandle_.product();

}

Member Data Documentation

CalibratedClusterPtrVector PFECALSuperClusterAlgo::_clustersEB

private

Definition at line 127 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and run().

CalibratedClusterPtrVector PFECALSuperClusterAlgo::_clustersEE

private

Definition at line 128 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and run().

clustering_type PFECALSuperClusterAlgo::_clustype

private

Definition at line 133 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setClusteringType().

energy_weight PFECALSuperClusterAlgo::_eweight

private

Definition at line 134 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setEnergyWeighting().

std::shared_ptr<PFEnergyCalibration> PFECALSuperClusterAlgo::_pfEnergyCalibration

private

Definition at line 132 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), setPFClusterCalibration(), and update().

bool PFECALSuperClusterAlgo::applyCrackCorrections_

private

Definition at line 164 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setCrackCorrections().

const EcalRecHitCollection* PFECALSuperClusterAlgo::barrelRecHits_

private

Definition at line 171 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and loadAndSortPFClusters().

const reco::BeamSpot* PFECALSuperClusterAlgo::beamSpot_

private

Definition at line 124 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and loadAndSortPFClusters().

const ESChannelStatus* PFECALSuperClusterAlgo::channelStatus_

private

Definition at line 125 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and update().

bool PFECALSuperClusterAlgo::doSatelliteClusterMerge_

private

Definition at line 158 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setSatelliteMerging().

bool PFECALSuperClusterAlgo::dropUnseedable_

private

Definition at line 160 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setDropUnseedable().

const reco::PFCluster::EEtoPSAssociation* PFECALSuperClusterAlgo::EEtoPS_

private

Definition at line 131 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and loadAndSortPFClusters().

const EcalRecHitCollection* PFECALSuperClusterAlgo::endcapRecHits_

private

Definition at line 172 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and loadAndSortPFClusters().

double PFECALSuperClusterAlgo::etawidthSuperClusterBarrel_

private

Definition at line 154 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setEtawidthSuperClusterBarrel().

double PFECALSuperClusterAlgo::etawidthSuperClusterEndcap_

private

Definition at line 156 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setEtawidthSuperClusterEndcap().

double PFECALSuperClusterAlgo::fractionForMajority_

private

Definition at line 159 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setMajorityFraction().

edm::EDGetTokenT<EcalRecHitCollection> PFECALSuperClusterAlgo::inputTagBarrelRecHits_

private

Definition at line 169 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and setTokens().

edm::EDGetTokenT<reco::BeamSpot> PFECALSuperClusterAlgo::inputTagBeamSpot_

private

Definition at line 122 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and setTokens().

edm::EDGetTokenT<EcalRecHitCollection> PFECALSuperClusterAlgo::inputTagEndcapRecHits_

private

Definition at line 170 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and setTokens().

edm::EDGetTokenT<edm::View<reco::PFCluster> > PFECALSuperClusterAlgo::inputTagPFClusters_

private

Definition at line 120 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and setTokens().

edm::EDGetTokenT<reco::PFCluster::EEtoPSAssociation> PFECALSuperClusterAlgo::inputTagPFClustersES_

private

Definition at line 121 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and setTokens().

bool PFECALSuperClusterAlgo::isOOTCollection_

private

Definition at line 168 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), loadAndSortPFClusters(), setIsOOTCollection(), and setTokens().

double PFECALSuperClusterAlgo::phiwidthSuperClusterBarrel_

private

Definition at line 153 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setPhiwidthSuperClusterBarrel().

double PFECALSuperClusterAlgo::phiwidthSuperClusterEndcap_

private

Definition at line 155 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setPhiwidthSuperClusterEndcap().

std::unique_ptr<SCEnergyCorrectorSemiParm> PFECALSuperClusterAlgo::regr_

private

Definition at line 144 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), loadAndSortPFClusters(), setTokens(), and update().

double PFECALSuperClusterAlgo::satelliteThreshold_

private

Definition at line 159 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setSatelliteThreshold().

std::unique_ptr<reco::SuperClusterCollection> PFECALSuperClusterAlgo::superClustersEB_

private

Definition at line 129 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), getEBOutputSCCollection(), and loadAndSortPFClusters().

std::unique_ptr<reco::SuperClusterCollection> PFECALSuperClusterAlgo::superClustersEE_

private

Definition at line 130 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), getEEOutputSCCollection(), and loadAndSortPFClusters().

bool PFECALSuperClusterAlgo::threshIsET_

private

Definition at line 165 of file PFECALSuperClusterAlgo.h.

Referenced by buildAllSuperClusters(), and setUseETForSeeding().

double PFECALSuperClusterAlgo::threshPFClusterBarrel_

private

Definition at line 149 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and setThreshPFClusterBarrel().

double PFECALSuperClusterAlgo::threshPFClusterEndcap_

private

Definition at line 151 of file PFECALSuperClusterAlgo.h.

Referenced by loadAndSortPFClusters(), and setThreshPFClusterEndcap().

double PFECALSuperClusterAlgo::threshPFClusterSeedBarrel_

private

Definition at line 148 of file PFECALSuperClusterAlgo.h.

Referenced by run(), and setThreshPFClusterSeedBarrel().

double PFECALSuperClusterAlgo::threshPFClusterSeedEndcap_

private

Definition at line 150 of file PFECALSuperClusterAlgo.h.

Referenced by run(), and setThreshPFClusterSeedEndcap().

double PFECALSuperClusterAlgo::threshSuperClusterEt_

private

Definition at line 146 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setThreshSuperClusterEt().

bool PFECALSuperClusterAlgo::useDynamicDPhi_

private

Definition at line 162 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setUseDynamicDPhi().

bool PFECALSuperClusterAlgo::useRegression_

private

Definition at line 143 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), loadAndSortPFClusters(), setTokens(), setUseRegression(), and update().

bool PFECALSuperClusterAlgo::verbose_

private

Definition at line 140 of file PFECALSuperClusterAlgo.h.

Referenced by buildSuperCluster(), and setVerbosityLevel().