l1tpf_calo::SimpleCaloLinker Class Reference

#include <CaloClusterer.h>

Inheritance diagram for l1tpf_calo::SimpleCaloLinker:

Public Member Functions
void	clear () override
void	run () override
	SimpleCaloLinker (const edm::ParameterSet &pset, const SingleCaloClusterer &ecal, const SingleCaloClusterer &hcal)
	~SimpleCaloLinker () override
Public Member Functions inherited from l1tpf_calo::SimpleCaloLinkerBase
void	clearBase ()
template<typename Corrector >
void	correct (const Corrector &corrector)
std::unique_ptr< l1t::PFClusterCollection >	fetch () const
std::unique_ptr< l1t::PFClusterCollection >	fetch (const edm::OrphanHandle< l1t::PFClusterCollection > &ecal, const edm::OrphanHandle< l1t::PFClusterCollection > &hcal) const
	SimpleCaloLinkerBase (const edm::ParameterSet &pset, const SingleCaloClusterer &ecal, const SingleCaloClusterer &hcal)
virtual	~SimpleCaloLinkerBase ()

Protected Attributes
PreClusterGrid	ecalToHCal_
Protected Attributes inherited from l1tpf_calo::SimpleCaloLinkerBase
IndexGrid	clusterIndex_
std::vector< CombinedCluster >	clusters_
const SingleCaloClusterer &	ecal_
std::vector< double >	etaBounds_
const Grid *	grid_
const SingleCaloClusterer &	hcal_
float	hoeCut_
std::vector< unsigned int >	maxClustersEtaPhi_
float	minHadronEt_
float	minHadronRawEt_
float	minPhotonEt_
bool	noEmInHGC_
std::vector< double >	phiBounds_

Detailed Description

Definition at line 298 of file CaloClusterer.h.

Constructor & Destructor Documentation

SimpleCaloLinker()

l1tpf_calo::SimpleCaloLinker::SimpleCaloLinker	(	const edm::ParameterSet &	pset,
		const SingleCaloClusterer &	ecal,
		const SingleCaloClusterer &	hcal
	)

Definition at line 538 of file CaloClusterer.cc.

    : SimpleCaloLinkerBase(pset, ecal, hcal), ecalToHCal_(*grid_) {}

~SimpleCaloLinker()

l1tpf_calo::SimpleCaloLinker::~SimpleCaloLinker ( )

override

Definition at line 543 of file CaloClusterer.cc.

{}

Member Function Documentation

clear()

void l1tpf_calo::SimpleCaloLinker::clear ( void  )

overridevirtual

Reimplemented from l1tpf_calo::SimpleCaloLinkerBase.

Definition at line 545 of file CaloClusterer.cc.

                                       {
  clearBase();
  ecalToHCal_.clear();
}

run()

void l1tpf_calo::SimpleCaloLinker::run ( )

overridevirtual

Implements l1tpf_calo::SimpleCaloLinkerBase.

Definition at line 550 of file CaloClusterer.cc.

References l1tpf_calo::GridData< T >::clear(), l1tpf_calo::CombinedCluster::clear(), l1tpf_calo::Cluster::constituents, SiPixelRawToDigiRegional_cfi::deltaPhi, l1tpf_calo::CombinedCluster::ecal_et, l1tpf_calo::CombinedCluster::ecal_eta, l1tpf_calo::CombinedCluster::ecal_phi, l1tnanotables_cff::et, l1tpf_calo::Cluster::et, l1tpf_calo::Cluster::eta, l1tpf_calo::CombinedCluster::hcal_et, mps_fire::i, l1tpf_calo::Cluster::phi, reco::reduceRange(), and compareTotals::tot.

                                     {
  unsigned int i, ncells = grid_->size();

  const EtGrid &hraw = hcal_.raw();
  const IndexGrid &ecals = ecal_.indexGrid();
  const IndexGrid &hcals = hcal_.indexGrid();

  // for each ECal cluster, get the corresponding HCal cluster and the sum of the neighbour HCal clusters
  ecalToHCal_.clear();
  for (i = 0; i < ncells; ++i) {
    if (ecals[i] >= 0) {
      if (hcals[i] >= 0) {
        ecalToHCal_[i].ptLocalMax = hcal_.cluster(i).et;
      } else {
        float tot = 0;
        for (int ineigh = 0; ineigh < 8; ++ineigh) {
          tot += hcal_.cluster(grid_->neighbour(i, ineigh)).et;
        }
        ecalToHCal_[i].ptOverNeighLocalMaxSum = tot ? ecal_.cluster(i).et / tot : 0;
      }
    }
  }

  clusterIndex_.fill(-1);
  clusters_.clear();
  CombinedCluster cluster;
  // promote HCal clusters to final clusters
  for (i = 0; i < ncells; ++i) {
    if (hcals[i] >= 0) {
      const Cluster &hcal = hcal_.cluster(i);
      cluster.clear();
      cluster.constituents.emplace_back(+i + 1, 1);
      if (ecalToHCal_[i].ptLocalMax > 0) {
        // direct linking is easy
        const Cluster &ecal = ecal_.cluster(i);
        if (ecal.et + hcal.et > minHadronRawEt_) {
          cluster.ecal_et = ecal.et;
          cluster.hcal_et = hcal.et;
          cluster.et = cluster.ecal_et + cluster.hcal_et;
          float wecal = cluster.ecal_et / cluster.et, whcal = 1.0 - wecal;
          cluster.eta = ecal.eta * wecal + hcal.eta * whcal;
          cluster.phi = ecal.phi * wecal + hcal.phi * whcal;
          cluster.ecal_eta = cluster.eta;
          cluster.ecal_phi = cluster.phi;
          // wrap around phi
          cluster.phi = reco::reduceRange(cluster.phi);
          cluster.constituents.emplace_back(-i - 1, 1);
        }
      } else {
        // sidewas linking is more annonying
        float myet = hcal.et;
        float etot = 0;
        float avg_eta = 0;
        float avg_phi = 0;
        for (int ineigh = 0; ineigh < 8; ++ineigh) {
          int ineighcell = grid_->neighbour(i, ineigh);
          if (ineighcell == -1)
            continue;  // skip dummy cells
          float fracet = myet * ecalToHCal_.neigh(i, ineigh).ptOverNeighLocalMaxSum;
          if (fracet == 0)
            continue;
          etot += fracet;
          avg_eta += fracet * (grid_->eta(ineighcell) - grid_->eta(i));
          avg_phi += fracet * deltaPhi(grid_->phi(ineighcell), grid_->phi(i));
          cluster.constituents.emplace_back(-i - 1, fracet / ecal_.cluster(ineighcell).et);
        }
        if (myet + etot > minHadronRawEt_) {
          cluster.hcal_et = hcal.et;
          cluster.ecal_et = etot;
          cluster.et = myet + etot;
          cluster.eta = hcal.eta + avg_eta / cluster.et;
          cluster.phi = hcal.phi + avg_phi / cluster.et;
          cluster.ecal_eta = cluster.eta;
          cluster.ecal_phi = cluster.phi;
          // wrap around phi
          cluster.phi = reco::reduceRange(cluster.phi);
        }
      }
      if (cluster.et > 0) {
        clusterIndex_[i] = clusters_.size();
        clusters_.push_back(cluster);
      }
    }
  }

  // promote Unlinked ECal clusters to final clusters
  for (i = 0; i < ncells; ++i) {
    if (ecals[i] >= 0 && ecalToHCal_[i].ptLocalMax == 0 && ecalToHCal_[i].ptOverNeighLocalMaxSum == 0) {
      cluster.clear();
      const Cluster &ecal = ecal_.cluster(i);
      cluster.ecal_et = ecal.et;
      cluster.hcal_et = hraw[i];
      cluster.et = cluster.ecal_et + cluster.hcal_et;
      cluster.eta = ecal.eta;
      cluster.phi = ecal.phi;
      cluster.constituents.emplace_back(-i - 1, 1);
      clusterIndex_[i] = clusters_.size();
      clusters_.push_back(cluster);
    }
  }
}

Member Data Documentation

ecalToHCal_

PreClusterGrid l1tpf_calo::SimpleCaloLinker::ecalToHCal_

protected

Definition at line 306 of file CaloClusterer.h.