ecal::WorkingCluster Class Reference

Public Member Functions
	WorkingCluster (const ldmx::EcalHit *eh, const ldmx::EcalGeometry &geom)
void	add (const ldmx::EcalHit *eh, const ldmx::EcalGeometry &geom)
void	add (const WorkingCluster &wc)
const TLorentzVector &	centroid () const
std::vector< const ldmx::EcalHit * >	getHits () const
bool	empty () const
void	clear ()

Private Attributes
std::vector< const ldmx::EcalHit * >	hits_
TLorentzVector	centroid_

Detailed Description

Definition at line 16 of file WorkingCluster.h.

Constructor & Destructor Documentation

WorkingCluster()

ecal::WorkingCluster::WorkingCluster	(	const ldmx::EcalHit *	eh,
		const ldmx::EcalGeometry &	geom
	)

Definition at line 11 of file WorkingCluster.cxx.

                                                            {
  add(eh, hex);
}

~WorkingCluster()

ecal::WorkingCluster::~WorkingCluster ( )

inline

Definition at line 20 of file WorkingCluster.h.

{};

Member Function Documentation

add() [1/2]

void ecal::WorkingCluster::add	(	const ldmx::EcalHit *	eh,
		const ldmx::EcalGeometry &	geom
	)

The ID number is implicitly converted to EcalID

Definition at line 16 of file WorkingCluster.cxx.

                                                      {
  double hitE = eh->getEnergy();
 
  auto [hitX, hitY, hitZ] = hex.getPosition(eh->getID());
 
  double newE = hitE + centroid_.E();
  double newCentroidX = (centroid_.Px() * centroid_.E() + hitE * hitX) / newE;
  double newCentroidY = (centroid_.Py() * centroid_.E() + hitE * hitY) / newE;
  double newCentroidZ = (centroid_.Pz() * centroid_.E() + hitE * hitZ) / newE;
 
  centroid_.SetPxPyPzE(newCentroidX, newCentroidY, newCentroidZ, newE);
 
  hits_.push_back(eh);
}

References ldmx::CalorimeterHit::getEnergy(), and ldmx::CalorimeterHit::getID().

add() [2/2]

void ecal::WorkingCluster::add

(

const WorkingCluster &

wc

)

Definition at line 33 of file WorkingCluster.cxx.

                                                 {
  double clusterE = wc.centroid().E();
  double centroidX = wc.centroid().Px();
  double centroidY = wc.centroid().Py();
  double centroidZ = wc.centroid().Pz();
 
  double newE = clusterE + centroid_.E();
  double newCentroidX =
      (centroid_.Px() * centroid_.E() + clusterE * centroidX) / newE;
  double newCentroidY =
      (centroid_.Py() * centroid_.E() + clusterE * centroidY) / newE;
  double newCentroidZ =
      (centroid_.Pz() * centroid_.E() + clusterE * centroidZ) / newE;
 
  centroid_.SetPxPyPzE(newCentroidX, newCentroidY, newCentroidZ, newE);
 
  std::vector<const ldmx::EcalHit*> clusterHits = wc.getHits();
 
  for (size_t i = 0; i < clusterHits.size(); i++) {
    hits_.push_back(clusterHits[i]);
  }
}

centroid()

const TLorentzVector & ecal::WorkingCluster::centroid ( ) const

inline

Definition at line 26 of file WorkingCluster.h.

{ return centroid_; }

clear()

void ecal::WorkingCluster::clear ( )

inline

Definition at line 32 of file WorkingCluster.h.

{ hits_.clear(); }

empty()

bool ecal::WorkingCluster::empty ( ) const

inline

Definition at line 30 of file WorkingCluster.h.

{ return hits_.empty(); }

getHits()

std::vector< const ldmx::EcalHit * > ecal::WorkingCluster::getHits ( ) const

inline

Definition at line 28 of file WorkingCluster.h.

{ return hits_; }

Member Data Documentation

centroid_

TLorentzVector ecal::WorkingCluster::centroid_

private

Definition at line 36 of file WorkingCluster.h.

hits_

std::vector<const ldmx::EcalHit*> ecal::WorkingCluster::hits_

private

Definition at line 35 of file WorkingCluster.h.

---

The documentation for this class was generated from the following files:

• Ecal/include/Ecal/WorkingCluster.h
• Ecal/src/Ecal/WorkingCluster.cxx