reco::PFCluster Class Reference

Particle flow cluster, see clustering algorithm in PFClusterAlgo. More...

#include <PFCluster.h>

Inheritance diagram for reco::PFCluster:

Public Types
typedef std::vector< std::pair < CaloClusterPtr::key_type, edm::Ptr< PFCluster > > >	EEtoPSAssociation
typedef ROOT::Math::PositionVector3D < ROOT::Math::CylindricalEta3D < double > >	REPPoint
Public Types inherited from reco::CaloCluster
enum	AlgoId {   island = 0, hybrid = 1, fixedMatrix = 2, dynamicHybrid = 3,   multi5x5 = 4, particleFlow = 5, hgcal_em = 6, hgcal_had = 7,   hgcal_mixed = 8, hfnose = 9, undefined = 1000 }
typedef AlgoId	AlgoID
enum	HCalFlags { badHcalMarker = 1 }
enum	SCFlags { cleanOnly = 0, common = 100, uncleanOnly = 200 }

Public Member Functions
void	addRecHitFraction (const reco::PFRecHitFraction &frac)
	add a given fraction of the rechit More...
void	calculatePositionREP ()
	computes posrep_ once and for all More...
double	charge () const
	dummy charge More...
double	depth () const
	cluster depth More...
double	energy () const
	cluster energy More...
PFLayer::Layer	layer () const
	cluster layer, see PFLayer.h in this directory More...
PFCluster &	operator= (const PFCluster &)
	PFCluster ()
	PFCluster (PFLayer::Layer layer, double energy, double x, double y, double z)
	constructor More...
const REPPoint &	positionREP () const
	cluster position: rho, eta, phi More...
template<typename pruner >
void	pruneUsing (pruner prune)
double	pt () const
	transverse momentum, massless approximation More...
const std::vector < reco::PFRecHitFraction > &	recHitFractions () const
	vector of rechit fractions More...
void	reset ()
	resets clusters parameters More...
void	resetHitsAndFractions ()
	reset only hits and fractions More...
void	setDepth (double depth)
void	setLayer (PFLayer::Layer layer)
	set layer More...
void	setTime (float time, float timeError=0)
void	setTimeError (float timeError)
double	theta () const
	angle More...
float	time () const
float	timeError () const
math::XYZPoint const &	vertex () const
	dummy vertex access More...
double	vx () const
double	vy () const
double	vz () const
Public Member Functions inherited from reco::CaloCluster
void	addHitAndFraction (DetId id, float fraction)
AlgoId	algo () const
	algorithm identifier More...
AlgoID	algoID () const
	CaloCluster ()
	default constructor. Sets energy and position to zero More...
	CaloCluster (AlgoID algoID)
	constructor with algoId, to be used in all child classes More...
	CaloCluster (double energy, const math::XYZPoint &position, const CaloID &caloID)
	CaloCluster (double energy, const math::XYZPoint &position)
	constructor from values More...
	CaloCluster (double energy, const math::XYZPoint &position, const CaloID &caloID, const AlgoID &algoID, uint32_t flags=0)
	CaloCluster (double energy, const math::XYZPoint &position, const CaloID &caloID, const std::vector< std::pair< DetId, float > > &usedHitsAndFractions, const AlgoId algoId, const DetId seedId=DetId(0), uint32_t flags=0)
	CaloCluster (double energy, const math::XYZPoint &position, float chi2, const std::vector< DetId > &usedHits, const AlgoId algoId, uint32_t flags=0)
	temporary compatibility constructor More...
const CaloID &	caloID () const
double	correctedEnergy () const
float	correctedEnergyUncertainty () const
double	energy () const
	cluster energy More...
double	eta () const
	pseudorapidity of cluster centroid More...
uint32_t	flags () const
const std::vector< std::pair < DetId, float > > &	hitsAndFractions () const
bool	isInClean () const
bool	isInUnclean () const
bool	operator< (const CaloCluster &rhs) const
	comparison < operator More...
bool	operator<= (const CaloCluster &rhs) const
	comparison <= operator More...
bool	operator== (const CaloCluster &rhs) const
	comparison == operator More...
bool	operator> (const CaloCluster &rhs) const
	comparison > operator More...
bool	operator>= (const CaloCluster &rhs) const
	comparison >= operator More...
double	phi () const
	azimuthal angle of cluster centroid More...
const math::XYZPoint &	position () const
	cluster centroid position More...
std::string	printHitAndFraction (unsigned i) const
	print hitAndFraction More...
void	reset ()
	resets the CaloCluster (position, energy, hitsAndFractions) More...
DetId	seed () const
	return DetId of seed More...
void	setAlgoId (const AlgoId &id)
void	setCaloId (const CaloID &id)
void	setCorrectedEnergy (double cenergy)
void	setCorrectedEnergyUncertainty (float energyerr)
void	setEnergy (double energy)
void	setFlags (uint32_t flags)
void	setPosition (const math::XYZPoint &p)
void	setSeed (const DetId &id)
size_t	size () const
	size in number of hits (e.g. in crystals for ECAL) More...
double	x () const
	x coordinate of cluster centroid More...
double	y () const
	y coordinate of cluster centroid More...
double	z () const
	z coordinate of cluster centroid More...
virtual	~CaloCluster ()
	destructor More...

Static Public Member Functions
static double	getDepthCorrection (double energy, bool isBelowPS=false, bool isHadron=false)

Private Attributes
double	depth_
PFLayer::Layer	layer_
	transient layer More...
REPPoint	posrep_
	cluster position: rho, eta, phi (transient) More...
std::vector < reco::PFRecHitFraction >	rechits_
	vector of rechit fractions (transient) More...
float	time_
	Michalis: add timing and depth information. More...
float	timeError_

Static Private Attributes
static const constexpr double	depthCorA_ = 0.89
	depth corrections More...
static const constexpr double	depthCorAp_ = 0.89
static const constexpr double	depthCorB_ = 7.3
static const constexpr double	depthCorBp_ = 4.0
static const math::XYZPoint	dummyVtx_

Additional Inherited Members
Protected Attributes inherited from reco::CaloCluster
AlgoID	algoID_
CaloID	caloID_
	bitmask for detector information More...
double	correctedEnergy_
float	correctedEnergyUncertainty_
double	energy_
	cluster energy More...
uint32_t	flags_
std::vector< std::pair< DetId, float > >	hitsAndFractions_
math::XYZPoint	position_
	cluster centroid position More...
DetId	seedId_
	DetId of seed. More...
Static Protected Attributes inherited from reco::CaloCluster
static const uint32_t	flagsMask_ = 0x0FFFFFFF
static const uint32_t	flagsOffset_ = 28

Detailed Description

Particle flow cluster, see clustering algorithm in PFClusterAlgo.

A particle flow cluster is defined by its energy and position, which are calculated from a vector of PFRecHitFraction. This calculation is performed in PFClusterAlgo.

Author

Colin Bernet

Date

July 2006

Definition at line 42 of file PFCluster.h.

Member Typedef Documentation

typedef std::vector<std::pair<CaloClusterPtr::key_type, edm::Ptr<PFCluster> > > reco::PFCluster::EEtoPSAssociation

Definition at line 44 of file PFCluster.h.

typedef ROOT::Math::PositionVector3D<ROOT::Math::CylindricalEta3D<double> > reco::PFCluster::REPPoint

Definition at line 48 of file PFCluster.h.

Constructor & Destructor Documentation

reco::PFCluster::PFCluster ( )

inline

Definition at line 50 of file PFCluster.h.

: CaloCluster(CaloCluster::particleFlow), time_(-99.0), layer_(PFLayer::NONE) {}

PFCluster::PFCluster	(	PFLayer::Layer	layer,
		double	energy,
		double	x,
		double	y,
		double	z
	)

constructor

Definition at line 10 of file PFCluster.cc.

    : CaloCluster(energy, math::XYZPoint(x, y, z), PFLayer::toCaloID(layer), CaloCluster::particleFlow),
      posrep_(position_.Rho(), position_.Eta(), position_.Phi()),
      time_(-99.),
      depth_(0.),
      layer_(layer) {}

Member Function Documentation

void PFCluster::addRecHitFraction

(

const reco::PFRecHitFraction &

frac

)

add a given fraction of the rechit

Definition at line 33 of file PFCluster.cc.

References reco::CaloCluster::addHitAndFraction(), reco::PFRecHitFraction::fraction(), reco::PFRecHitFraction::recHitRef(), and rechits_.

Referenced by PFMultiDepthClusterizer::absorbCluster(), Basic2DClusterForEachSeed::buildClusters(), GenericSimClusterMapper::buildClusters(), PFClusterFromHGCalTrackster::buildClusters(), RealisticSimClusterMapper::buildClusters(), Basic2DGenericTopoClusterizer::buildTopoCluster(), Basic2DGenericPFlowClusterizer::seedPFClustersFromTopo(), and PFlow2DClusterizerWithTime::seedPFClustersFromTopo().

                                                                  {
  rechits_.push_back(frac);

  addHitAndFraction(frac.recHitRef()->detId(), frac.fraction());
}

void reco::PFCluster::calculatePositionREP ( )

inline

computes posrep_ once and for all

Definition at line 95 of file PFCluster.h.

References reco::CaloCluster::position_, and posrep_.

Referenced by Basic2DClusterForEachSeed::buildClusters(), Cluster3DPCACalculator::calculateAndSetPositionActual(), ECAL2DPositionCalcWithDepthCorr::calculateAndSetPositionActual(), Basic2DGenericPFlowPositionCalc::calculateAndSetPositionActual(), PFElecTkProducer::isSharingEcalEnergyWithEgSC(), and ConvBremPFTrackFinder::runConvBremFinder().

{ posrep_.SetCoordinates(position_.Rho(), position_.Eta(), position_.Phi()); }

double reco::PFCluster::charge ( ) const

inline

dummy charge

some classes to make this fit into a template footprint for RecoPFClusterRefCandidate so we can make jets and MET out of PFClusters.

Definition at line 107 of file PFCluster.h.

{ return 0; }

double reco::PFCluster::depth ( ) const

inline

cluster depth

Definition at line 82 of file PFCluster.h.

References depth_.

Referenced by PFMultiDepthClusterizer::link(), and setDepth().

{ return depth_; }

double reco::PFCluster::energy ( ) const

inline

cluster energy

Definition at line 74 of file PFCluster.h.

References reco::CaloCluster::energy_.

Referenced by HGCalHitCalibration::analyze(), FWPFClusterRPZProxyBuilder::build(), FWPFClusterRPZUtils::buildRhoPhiClusterLineSet(), FWPFClusterRPZUtils::buildRhoZClusterLineSet(), ECAL2DPositionCalcWithDepthCorr::calculateAndSetPositionActual(), PFMultiDepthClusterizer::calculateShowerShapes(), PFClusterEMEnergyCorrector::correctEnergies(), PFPhotonTranslator::createBasicCluster(), PFEnergyCalibration::energyEm(), PFlow2DClusterizerWithTime::growPFClusters(), PFElecTkProducer::isSharingEcalEnergyWithEgSC(), PFMultiDepthClusterizer::link(), reco::operator<<(), reco::HGCalMultiCluster::operator>(), pt(), Jet.Jet::rawEnergy(), FWPFBlockProxyBuilder::setupClusterElement(), FWPFClusterRPZProxyBuilder::sharedBuild(), and FWPFClusterLegoProxyBuilder::sharedBuild().

{ return energy_; }

double PFCluster::getDepthCorrection ( double  energy, bool  isBelowPS = false, bool  isHadron = false  )

static

Definition at line 39 of file PFCluster.cc.

References depthCorA_, depthCorAp_, depthCorB_, depthCorBp_, and log.

Referenced by l1tpf::propagateToCalo(), and LowPtGsfElectronSeedProducer::propagateTrackToCalo().

                                                                                 {
  double corrA = depthCorA_;
  double corrB = depthCorB_;
  if (isBelowPS) {
    corrA = depthCorAp_;
    corrB = depthCorBp_;
  }
  return isHadron ? corrA : corrA * (corrB + log(energy));
}

PFLayer::Layer PFCluster::layer

(

)

const

cluster layer, see PFLayer.h in this directory

Definition at line 56 of file PFCluster.cc.

References reco::CaloCluster::caloID(), PFLayer::fromCaloID(), layer_, and PFLayer::NONE.

Referenced by EgammaIsoESDetIdCollectionProducer::addDetIds(), FWPFEcalClusterRPZProxyBuilder::build(), FWPFEcalClusterLegoProxyBuilder::build(), FWPFHcalClusterLegoProxyBuilder::build(), PFClusterEMEnergyCorrector::correctEnergies(), geometryXMLparser.DTAlignable::index(), geometryXMLparser.CSCAlignable::index(), reco::operator<<(), PFAlgo::reconstructCluster(), KDTreeLinkerPSEcal::searchLinks(), setLayer(), LinkByRecHit::testECALAndPSByRecHit(), and LinkByRecHit::testTrackAndClusterByRecHit().

                                    {
  // cout<<"calling PFCluster::layer "<<caloID()<<" "<<PFLayer::fromCaloID( caloID() )<<endl;
  if (layer_ != PFLayer::NONE)
    return layer_;
  return PFLayer::fromCaloID(caloID());
}

PFCluster & PFCluster::operator=

(

const PFCluster &

other

)

Definition at line 63 of file PFCluster.cc.

References reco::CaloCluster::energy_, reco::CaloCluster::position_, posrep_, and rechits_.

Referenced by reco::PFSuperCluster::operator=().

                                                      {
  CaloCluster::operator=(other);
  rechits_ = other.rechits_;
  energy_ = other.energy_;
  position_ = other.position_;
  posrep_ = other.posrep_;

  return *this;
}

const REPPoint& reco::PFCluster::positionREP ( ) const

inline

cluster position: rho, eta, phi

Definition at line 92 of file PFCluster.h.

References posrep_.

Referenced by PFMultiDepthClusterizer::calculateShowerShapes(), PFMultiDepthClusterizer::link(), reco::operator<<(), KDTreeLinkerPSEcal::searchLinks(), and LinkByRecHit::testTrackAndClusterByRecHit().

{ return posrep_; }

template<typename pruner >

void reco::PFCluster::pruneUsing ( pruner  prune )

inline

Definition at line 123 of file PFCluster.h.

References first, reco::CaloCluster::hitsAndFractions_, mps_fire::i, eostools::move(), and rechits_.

                                  {
      // remove_if+erase algo applied to both vectors...
      auto iter = std::find_if_not(rechits_.begin(), rechits_.end(), prune);
      if (iter == rechits_.end())
        return;
      auto first = iter - rechits_.begin();
      for (auto i = first; ++i < int(rechits_.size());) {
        if (prune(rechits_[i])) {
          rechits_[first] = std::move(rechits_[i]);
          hitsAndFractions_[first] = std::move(hitsAndFractions_[i]);
          ++first;
        }
      }
      rechits_.erase(rechits_.begin() + first, rechits_.end());
      hitsAndFractions_.erase(hitsAndFractions_.begin() + first, hitsAndFractions_.end());
    }

double reco::PFCluster::pt ( ) const

inline

transverse momentum, massless approximation

Definition at line 110 of file PFCluster.h.

References energy(), reco::CaloCluster::position_, and funct::sin().

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), PFClusterEMEnergyCorrector::correctEnergies(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

{ return (energy() * sin(position_.theta())); }

const std::vector<reco::PFRecHitFraction>& reco::PFCluster::recHitFractions ( ) const

inline

vector of rechit fractions

Definition at line 65 of file PFCluster.h.

References rechits_.

Referenced by PFMultiDepthClusterizer::absorbCluster(), Cluster3DPCACalculator::calculateAndSetPositionActual(), ECAL2DPositionCalcWithDepthCorr::calculateAndSetPositionActual(), Basic2DGenericPFlowPositionCalc::calculateAndSetPositionActual(), PFMultiDepthClusterizer::calculateShowerShapes(), PFlow2DClusterizerWithTime::clusterTimeResolution(), PFlow2DClusterizerWithTime::clusterTimeResolutionFromSeed(), EcalBarrelClusterFastTimer::getTimeForECALPFCluster(), Basic2DGenericPFlowClusterizer::growPFClusters(), PFlow2DClusterizerWithTime::growPFClusters(), reco::operator<<(), Basic2DGenericPFlowClusterizer::seedPFClustersFromTopo(), PFlow2DClusterizerWithTime::seedPFClustersFromTopo(), PFAlgo::setHcalDepthInfo(), LinkByRecHit::testECALAndPSByRecHit(), and LinkByRecHit::testTrackAndClusterByRecHit().

{ return rechits_; }

void PFCluster::reset

(

void

)

resets clusters parameters

Definition at line 17 of file PFCluster.cc.

References reco::CaloCluster::energy_, layer_, PFLayer::NONE, reco::CaloCluster::position_, posrep_, rechits_, reco::CaloCluster::reset(), and time_.

Referenced by Basic2DClusterForEachSeed::buildClusters(), and reco::PFSuperCluster::reset().

                      {
  energy_ = 0;
  position_ *= 0;
  posrep_ *= 0;
  time_ = -99.;
  layer_ = PFLayer::NONE;
  rechits_.clear();

  CaloCluster::reset();
}

void PFCluster::resetHitsAndFractions

(

)

reset only hits and fractions

Definition at line 28 of file PFCluster.cc.

References reco::CaloCluster::hitsAndFractions_, and rechits_.

                                      {
  rechits_.clear();
  hitsAndFractions_.clear();
}

void reco::PFCluster::setDepth ( double  depth )

inline

Definition at line 89 of file PFCluster.h.

References depth(), and depth_.

Referenced by Basic2DClusterForEachSeed::buildClusters(), and Basic2DGenericPFlowPositionCalc::calculateAndSetPositionActual().

{ depth_ = depth; }

void PFCluster::setLayer

(

PFLayer::Layer

layer

)

set layer

Definition at line 49 of file PFCluster.cc.

References reco::CaloCluster::caloID_, layer(), layer_, and PFLayer::toCaloID().

Referenced by Basic2DClusterForEachSeed::buildClusters(), Cluster3DPCACalculator::calculateAndSetPositionActual(), ECAL2DPositionCalcWithDepthCorr::calculateAndSetPositionActual(), Basic2DGenericPFlowPositionCalc::calculateAndSetPositionActual(), and reco::HGCalMultiCluster::HGCalMultiCluster().

                                           {
  // cout<<"calling PFCluster::setLayer "<<layer<<endl;
  layer_ = layer;
  caloID_ = PFLayer::toCaloID(layer);
  // cout<<"done "<<caloID_<<endl;
}

void reco::PFCluster::setTime ( float  time, float  timeError = 0  )

inline

Definition at line 84 of file PFCluster.h.

References time(), time_, timeError(), and timeError_.

Referenced by Basic2DClusterForEachSeed::buildClusters(), RealisticSimClusterMapper::buildClusters(), Cluster3DPCACalculator::calculateAndSetPositionActual(), ECAL2DPositionCalcWithDepthCorr::calculateAndSetPositionActual(), Basic2DGenericPFlowPositionCalc::calculateAndSetPositionActual(), PFlow2DClusterizerWithTime::clusterTimeResolution(), PFlow2DClusterizerWithTime::clusterTimeResolutionFromSeed(), and MultiClustersFromTrackstersProducer::produce().

                                                  {
      time_ = time;
      timeError_ = timeError;
    }

void reco::PFCluster::setTimeError ( float  timeError )

inline

Definition at line 88 of file PFCluster.h.

References timeError(), and timeError_.

Referenced by ECAL2DPositionCalcWithDepthCorr::calculateAndSetPositionActual(), Basic2DGenericPFlowPositionCalc::calculateAndSetPositionActual(), PFlow2DClusterizerWithTime::clusterTimeResolution(), and PFlow2DClusterizerWithTime::clusterTimeResolutionFromSeed().

{ timeError_ = timeError; }

double reco::PFCluster::theta ( ) const

inline

angle

Definition at line 113 of file PFCluster.h.

References reco::CaloCluster::position_.

Referenced by Tau.Tau::zImpact().

{ return position_.theta(); }

float reco::PFCluster::time ( ) const

inline

Returns

cluster time

Definition at line 77 of file PFCluster.h.

References time_.

Referenced by PFlow2DClusterizerWithTime::dist2Time(), and setTime().

{ return time_; }

float reco::PFCluster::timeError ( ) const

inline

Returns

the timing uncertainty

Definition at line 79 of file PFCluster.h.

References timeError_.

Referenced by setTime(), and setTimeError().

{ return timeError_; }

math::XYZPoint const& reco::PFCluster::vertex ( ) const

inline

dummy vertex access

Definition at line 116 of file PFCluster.h.

References dummyVtx_.

Referenced by Tau.Tau::dxy(), vx(), vy(), and vz().

{ return dummyVtx_; }

double reco::PFCluster::vx ( ) const

inline

Definition at line 117 of file PFCluster.h.

References vertex().

{ return vertex().x(); }

double reco::PFCluster::vy ( ) const

inline

Definition at line 118 of file PFCluster.h.

References vertex().

{ return vertex().y(); }

double reco::PFCluster::vz ( ) const

inline

Definition at line 119 of file PFCluster.h.

References vertex().

{ return vertex().z(); }

Member Data Documentation

double reco::PFCluster::depth_

private

Definition at line 150 of file PFCluster.h.

Referenced by depth(), and setDepth().

const constexpr double reco::PFCluster::depthCorA_ = 0.89

staticprivate

depth corrections

Definition at line 156 of file PFCluster.h.

Referenced by getDepthCorrection().

const constexpr double reco::PFCluster::depthCorAp_ = 0.89

staticprivate

Definition at line 158 of file PFCluster.h.

Referenced by getDepthCorrection().

const constexpr double reco::PFCluster::depthCorB_ = 7.3

staticprivate

Definition at line 157 of file PFCluster.h.

Referenced by getDepthCorrection().

const constexpr double reco::PFCluster::depthCorBp_ = 4.0

staticprivate

Definition at line 159 of file PFCluster.h.

Referenced by getDepthCorrection().

const math::XYZPoint PFCluster::dummyVtx_

staticprivate

Definition at line 161 of file PFCluster.h.

Referenced by vertex().

PFLayer::Layer reco::PFCluster::layer_

private

transient layer

Definition at line 153 of file PFCluster.h.

Referenced by layer(), reset(), and setLayer().

REPPoint reco::PFCluster::posrep_

private

cluster position: rho, eta, phi (transient)

Definition at line 146 of file PFCluster.h.

Referenced by calculatePositionREP(), operator=(), positionREP(), and reset().

std::vector<reco::PFRecHitFraction> reco::PFCluster::rechits_

private

vector of rechit fractions (transient)

Definition at line 143 of file PFCluster.h.

Referenced by addRecHitFraction(), operator=(), pruneUsing(), recHitFractions(), reset(), and resetHitsAndFractions().

float reco::PFCluster::time_

private

Michalis: add timing and depth information.

Definition at line 149 of file PFCluster.h.

Referenced by reset(), setTime(), and time().

float reco::PFCluster::timeError_

private

Definition at line 149 of file PFCluster.h.

Referenced by setTime(), setTimeError(), and timeError().