Print

Print


Commit in java/trunk/users/src/main/java/org/lcsim/hps/users/phansson on MAIN
TwoTrackAnlysis.java+21-388 -> 89
Add region count

java/trunk/users/src/main/java/org/lcsim/hps/users/phansson
TwoTrackAnlysis.java 88 -> 89
--- java/trunk/users/src/main/java/org/lcsim/hps/users/phansson/TwoTrackAnlysis.java	2013-12-13 21:11:11 UTC (rev 88)
+++ java/trunk/users/src/main/java/org/lcsim/hps/users/phansson/TwoTrackAnlysis.java	2013-12-13 21:11:45 UTC (rev 89)
@@ -68,6 +68,7 @@
     private int totalTwoTrackEvents=0;
     private int totalMCEvents=0;
     private int totalTwoTrackMCEvents=0;
+	private int nEventsWithGoodRegionCluster=0;
     private boolean hideFrame = false;
     private String outputPlotFileName;
     private String ecalClusterCollectionName = "EcalClusters";
@@ -321,6 +322,21 @@
             if(_debug) {
                 System.out.println(this.getClass().getSimpleName() + ": found " + clusters.size() + " ecal clusters " + event.getEventNumber());
             }
+
+			boolean goodRegion = false;
+			for(HPSEcalCluster c : clusters) {
+                int iy = c.getSeedHit().getIdentifierFieldValue("iy");
+                int ix = c.getSeedHit().getIdentifierFieldValue("ix");
+                double E = c.getEnergy();
+				int evtnr = event.getEventNumber();
+                int clsize = c.getSize();
+				if( iy>0 && ix>0 && E>0.6) {
+					//printWriter.format("%d %5.5f %5d %5d %5d ",evtnr,E,ix,iy,clsize);
+					goodRegion = true;
+				}
+			}
+			if( goodRegion ) nEventsWithGoodRegionCluster++;
+
         }
         
         Hep3Vector vtxPosMC = null;
@@ -384,9 +400,11 @@
     public void endOfData() {
         
         System.out.println(this.getClass().getSimpleName() + ": Total Number of Events = "+this.totalEvents);
-        System.out.println(this.getClass().getSimpleName() + ": Total Number of Two-Track events Processed = "+this.totalTwoTrackEvents);
-        System.out.println(this.getClass().getSimpleName() + ": Total Number of MCEvents = "+this.totalMCEvents);
-        System.out.println(this.getClass().getSimpleName() + ": Total Number of Two-Track MC events Processed = "+this.totalTwoTrackMCEvents);
+        System.out.println(this.getClass().getSimpleName() + ": Total Number of Events filling tuple = "+this.totalTwoTrackEvents);
+        System.out.println(this.getClass().getSimpleName() + ": Total Number of events with MC collection = "+this.totalMCEvents);
+        System.out.println(this.getClass().getSimpleName() + ": Total Number of events with e+e- ID'd = "+this.totalTwoTrackMCEvents);
+        System.out.println(this.getClass().getSimpleName() + ": Total Number of events with cl in good region = "+this.nEventsWithGoodRegionCluster);
+		
         
         if (!"".equals(outputPlotFileName)) {
             try {
SVNspam 0.1