Author: [log in to unmask]
Date: Thu Jan 15 17:01:07 2015
New Revision: 1940
Log:
Add methods for finding hits that are shared or unshared in a cluster using contribution list comparison.
Modified:
java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ClusterUtilities.java
Modified: java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ClusterUtilities.java
=============================================================================
--- java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ClusterUtilities.java (original)
+++ java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ClusterUtilities.java Thu Jan 15 17:01:07 2015
@@ -429,4 +429,46 @@
}
return hits;
}
+
+ /**
+ * This method will determine which hits are likely shared with
+ * other clusters, by comparing a hit's energy to its contribution.
+ * If the energy is greater than the contribution, then its energy
+ * will assumed to be shared between clusters. This convention
+ * might not be true depending on the clustering algorithm, but
+ * it is probably the best that can be done without the complete
+ * list of clusters in the event.
+ * @param cluster The input cluster.
+ */
+ public static List<CalorimeterHit> findSharedHits(Cluster cluster) {
+ List<CalorimeterHit> sharedHits = new ArrayList<CalorimeterHit>();
+ for (int i = 0; i < cluster.getCalorimeterHits().size(); i++) {
+ if (cluster.getCalorimeterHits().get(i).getCorrectedEnergy() > cluster.getHitContributions()[i]) {
+ sharedHits.add(cluster.getCalorimeterHits().get(i));
+ }
+ }
+ return sharedHits;
+ }
+
+ /**
+ * This method will determine which hits are likely not shared
+ * other clusters, by comparing a hit's energy to its contribution.
+ * If the energy is equal to the contribution, then its energy
+ * will assumed to not be shared between clusters. This convention
+ * might not be true depending on the clustering algorithm, but
+ * it is probably the best that can be done without the complete list
+ * of clusters in the event.
+ * @param cluster The input cluster.
+ */
+ public static List<CalorimeterHit> findUnsharedHits(Cluster cluster) {
+ List<CalorimeterHit> uniqueHits = new ArrayList<CalorimeterHit>();
+ for (int i = 0; i < cluster.getCalorimeterHits().size(); i++) {
+ if (cluster.getCalorimeterHits().get(i).getCorrectedEnergy() == cluster.getHitContributions()[i]) {
+ uniqueHits.add(cluster.getCalorimeterHits().get(i));
+ }
+ }
+ return uniqueHits;
+ }
+
+
}
|