ecal::MyClusterWeight Class Reference

Computes the weight (distance) between two Ecal clusters. More...

#include <MyClusterWeight.h>

Public Types
using	ClusterType = recon::WorkingCluster<ldmx::EcalHit>

Public Member Functions
double	operator() (const ClusterType &a, const ClusterType &b)
	Compute the weight between two clusters.

Detailed Description

Computes the weight (distance) between two Ecal clusters.

The weight is used by TemplatedClusterFinder to decide which clusters to merge. Smaller weights indicate clusters that should be merged first. The weight is based on both transverse and longitudinal separation.

Definition at line 23 of file MyClusterWeight.h.

Member Typedef Documentation

ClusterType

using ecal::MyClusterWeight::ClusterType = recon::WorkingCluster<ldmx::EcalHit>

Definition at line 25 of file MyClusterWeight.h.

Member Function Documentation

operator()()

double ecal::MyClusterWeight::operator() ( const ClusterType & a, const ClusterType & b )

inline

Compute the weight between two clusters.

Parameters

a	First cluster
b	Second cluster

Returns

Weight value (smaller = should merge first)

Definition at line 34 of file MyClusterWeight.h.

                                                                {
    // Moliere radius of detector, roughly. In mm
    double rmol = 10.00;
    // Characteristic cluster longitudinal variable in mm
    double dzchar = 100.0;
 
    double a_e = a.energy();
    double a_x = a.centroidX();
    double a_y = a.centroidY();
    double a_z = a.centroidZ();
 
    double b_e = b.energy();
    double b_x = b.centroidX();
    double b_y = b.centroidY();
    double b_z = b.centroidZ();
 
    double dijz;
    if (a_e >= b_e) {
      dijz = b_z - a_z;
    } else {
      dijz = a_z - b_z;
    }
 
    double dij_t = std::sqrt(std::pow(a_x - b_x, 2) + std::pow(a_y - b_y, 2));
 
    double weight_t = std::exp(std::pow(dij_t / rmol, 2)) - 1;
    double weight_z = std::exp(std::abs(dijz) / dzchar) - 1;
 
    // Return the highest of the two weights
    return std::max(weight_t, weight_z);
  }

References recon::WorkingCluster< HitType >::centroidX(), recon::WorkingCluster< HitType >::centroidY(), recon::WorkingCluster< HitType >::centroidZ(), and recon::WorkingCluster< HitType >::energy().

---

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

• Ecal/include/Ecal/MyClusterWeight.h