hps-java/src/main/java/org/lcsim/hps/analysis/ecal
diff -u -r1.2 -r1.3
--- HPSEcalPlotsDriver.java 16 May 2011 19:31:28 -0000 1.2
+++ HPSEcalPlotsDriver.java 25 May 2011 20:55:52 -0000 1.3
@@ -28,7 +28,7 @@
* Diagnostic plots for HPS ECal.
*
* @author Jeremy McCormick <[log in to unmask]>
- * @version $Id: HPSEcalPlotsDriver.java,v 1.2 2011/05/16 19:31:28 jeremy Exp $
+ * @version $Id: HPSEcalPlotsDriver.java,v 1.3 2011/05/25 20:55:52 jeremy Exp $
*/
// TODO ZY plot can use h2d with small bins in Y to avoid weird binning effect
// TODO ZY of cluster positions (algo from Tim to find pos)
@@ -70,6 +70,7 @@
ICloud1D clusSeedEPlot;
ICloud1D clusSeedDistPlot;
ICloud2D leadClusAndPrimaryPlot;
+ IHistogram1D clusNHits;
//ICloud1D hit50MeVPlot;
//ICloud1D hit30MeVPlot;
@@ -218,6 +219,9 @@
clusSeedEPlot = aida.cloud1D(ecalClusterCollectionName + " : Cluster Seed Hit Distance from Beam Axis");
clusSeedEPlot.annotation().addItem("xAxisLabel", "Distance [mm]");
+ clusNHits = aida.histogram1D(ecalClusterCollectionName + " : Number of Hits per Cluster", 13, 0, 13);
+ clusNHits.annotation().addItem("xAxisLabel", "Number of Hits");
+
leadClusAndPrimaryPlot = aida.cloud2D(ecalClusterCollectionName + " : Lead Cluster E vs Highest Primary Particle E");
}
@@ -390,6 +394,8 @@
Set<CalorimeterHit> naughty = new HashSet<CalorimeterHit>();
for (Cluster clus : clusters)
{
+ clusNHits.fill(clus.getCalorimeterHits().size());
+
double e = clus.getEnergy();
clusEPlot.fill(e);
clusE += e;