hcal::MyClusterWeight Class Reference

Public Member Functions
double	operator() (const WorkingCluster &a, const WorkingCluster &b)

Detailed Description

Definition at line 11 of file MyClusterWeight.h.

Member Function Documentation

operator()()

double hcal::MyClusterWeight::operator() ( const WorkingCluster &  a, const WorkingCluster &  b  )

inline

Definition at line 13 of file MyClusterWeight.h.

                               {  // returns weighting function, where smallest
                                  // weights will be combined first
 
    double rmol = 10.00;    // Moliere radius of detector, roughly. In mm TODO
    double dzchar = 100.0;  // lateral shower development in mm TODO
 
    double aE = a.centroid().E();
    double aX = a.centroid().Px();
    double aY = a.centroid().Py();
    double aZ = a.centroid().Pz();
 
    double bE = b.centroid().E();
    double bX = b.centroid().Px();
    double bY = b.centroid().Py();
    double bZ = b.centroid().Pz();
 
    double dijz;
    if (aE >= bE) {
      // differences in Z
      dijz = bZ - aZ;
    } else {
      dijz = aZ - bZ;
    }
 
    // Transverse Difference
    double dijT = pow(pow(aX - bX, 2) + pow(aY - bY, 2), 0.5);
    // Trans --> massive
    double weightT = exp(pow(dijT / rmol, 2)) - 1;
    // Long
    double weightZ = (exp(abs(dijz) / dzchar) - 1);
 
    // Return the highest of the two weights
    if (weightT <= weightZ) {
      return weightZ;
    } else {
      return weightT;
    }
  }

---

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

• Hcal/include/Hcal/MyClusterWeight.h