Print

Print


Author: [log in to unmask]
Date: Mon Sep 28 10:54:33 2015
New Revision: 3728

Log:
fix bugs

Modified:
    java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnDriver.java
    java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnSSPDriver.java

Modified: java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnDriver.java
 =============================================================================
--- java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnDriver.java	(original)
+++ java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnDriver.java	Mon Sep 28 10:54:33 2015
@@ -34,7 +34,7 @@
  */
 public class TriggerTurnOnDriver extends Driver {
 
-    private static Logger logger = LogUtil.create(TriggerTurnOnDriver.class, new BasicLogFormatter(), Level.FINE);
+    private static Logger logger = LogUtil.create(TriggerTurnOnDriver.class, new BasicLogFormatter(), Level.INFO);
     private final String ecalClusterCollectionName = "EcalClustersCorr";
     IPlotter plotter;
     IPlotter plotter2;
@@ -163,22 +163,25 @@
                 if( Math.abs(clusterPosIdy) > 5 ) 
                     throw new RuntimeException("invalid crystal position " + clusterPosIdy);
                 int half = clusterPosIdy > 0 ? 0 : 1;
-                clusterEOne_Random_thetaY[half][Math.abs(clusterPosIdy)-1].fill(clusterEMax.getEnergy());
-                
-            }
-        }
+                int ypos = Math.abs(clusterPosIdy)-1;
+                clusterEOne_Random_thetaY[half][ypos].fill(clusterEMax.getEnergy());
+            }
+        }
+
+
 
         // fill numerator
         if (triggerDecisions.passed(TriggerType.SINGLES1_SIM)) {
             logger.fine("Eureka. They both fired.");
             if(clusterEMax != null) {
                 clusterE_RandomSingles1.fill(clusterEMax.getEnergy());
-                if(clusters.size() == 1) 
+                if(clusters.size() == 1) {
                     clusterEOne_RandomSingles1.fill(clusterEMax.getEnergy());
                     int clusterPosIdy = clusterEMax.getCalorimeterHits().get(0).getIdentifierFieldValue("iy");
                     int half = clusterPosIdy > 0 ? 0 : 1;
-                    clusterEOne_RandomSingles1_thetaY[half][Math.abs(clusterPosIdy)-1].fill(clusterEMax.getEnergy());
-
+                    int ypos = Math.abs(clusterPosIdy)-1;
+                    clusterEOne_RandomSingles1_thetaY[half][ypos].fill(clusterEMax.getEnergy());
+                }
             }
         }
 
@@ -196,14 +199,13 @@
         int r = 0;
         for(int i=0;i<2;++i) {
             for(int y=0;y<5;++y) {
-                clusterEOne_RandomSingles1_thetaY_trigEff[i][y] = aida.histogramFactory().divide("trigEffEone", clusterEOne_Random_thetaY[i][y], clusterEOne_RandomSingles1_thetaY[i][y]);
-                plotter333.region(r);
+                clusterEOne_RandomSingles1_thetaY_trigEff[i][y] = aida.histogramFactory().divide("trigEffEone_" + (i==0?"top":"bottom") + "_" + y, clusterEOne_RandomSingles1_thetaY[i][y], clusterEOne_Random_thetaY[i][y]);
+                plotter333.region(r).plot(clusterEOne_RandomSingles1_thetaY_trigEff[i][y]);
                 r++;
             }
         }
         logger.info("entries in clusterE_RandomSingles1_trigEff: " + Integer.toString(clusterE_RandomSingles1_trigEff.allEntries()));
         plotter.region(2).plot(clusterE_RandomSingles1_trigEff);
-        plotter2.region(2).plot(clusterEOne_RandomSingles1_trigEff);
         plotter2.region(2).plot(clusterEOne_RandomSingles1_trigEff);
         
     }

Modified: java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnSSPDriver.java
 =============================================================================
--- java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnSSPDriver.java	(original)
+++ java/trunk/analysis/src/main/java/org/hps/analysis/trigger/TriggerTurnOnSSPDriver.java	Mon Sep 28 10:54:33 2015
@@ -34,7 +34,7 @@
  */
 public class TriggerTurnOnSSPDriver extends Driver {
     
-    private static Logger logger = LogUtil.create(TriggerTurnOnSSPDriver.class, new BasicLogFormatter(), Level.FINE);
+    private static Logger logger = LogUtil.create(TriggerTurnOnSSPDriver.class, new BasicLogFormatter(), Level.INFO);
     private final String ecalClusterCollectionName = "EcalClustersCorr";
     IPlotter plotter;
     IPlotter plotter2;
@@ -174,12 +174,12 @@
             logger.fine("Eureka. They both fired.");
             if(clusterEMax != null) {
                 clusterE_RandomSingles1.fill(clusterEMax.getEnergy());
-                if(clusters.size() == 1) 
+                if(clusters.size() == 1) {
                     clusterEOne_RandomSingles1.fill(clusterEMax.getEnergy());
                     int clusterPosIdy = clusterEMax.getYIndex();
                     int half = clusterPosIdy > 0 ? 0 : 1;
                     clusterEOne_RandomSingles1_thetaY[half][Math.abs(clusterPosIdy)-1].fill(clusterEMax.getEnergy());
-
+                }
             }
         }
 
@@ -192,19 +192,18 @@
     protected void endOfData() {
         logger.info("Processed " + nEventsProcessed);
         logger.info("nResSingles1 " + nResultSingles1 + " nSimSingles1 " + nSimSingles1);
-        clusterE_RandomSingles1_trigEff = aida.histogramFactory().divide("trigEff", clusterE_RandomSingles1, clusterE_Random);
-        clusterEOne_RandomSingles1_trigEff = aida.histogramFactory().divide("trigEffEone", clusterEOne_RandomSingles1, clusterEOne_Random);
+        clusterE_RandomSingles1_trigEff = aida.histogramFactory().divide("SSP_trigEff", clusterE_RandomSingles1, clusterE_Random);
+        clusterEOne_RandomSingles1_trigEff = aida.histogramFactory().divide("SSP_trigEffEone", clusterEOne_RandomSingles1, clusterEOne_Random);
         int r = 0;
         for(int i=0;i<2;++i) {
             for(int y=0;y<5;++y) {
-                clusterEOne_RandomSingles1_thetaY_trigEff[i][y] = aida.histogramFactory().divide("trigEffEone", clusterEOne_Random_thetaY[i][y], clusterEOne_RandomSingles1_thetaY[i][y]);
-                plotter333.region(r);
+                clusterEOne_RandomSingles1_thetaY_trigEff[i][y] = aida.histogramFactory().divide("SSP_trigEffEone_" + (i==0?"top":"bottom") + "_" + y, clusterEOne_RandomSingles1_thetaY[i][y], clusterEOne_Random_thetaY[i][y]);
+                plotter333.region(r).plot(clusterEOne_RandomSingles1_thetaY_trigEff[i][y]);
                 r++;
             }
         }
         logger.info("entries in clusterE_RandomSingles1_trigEff: " + Integer.toString(clusterE_RandomSingles1_trigEff.allEntries()));
         plotter.region(2).plot(clusterE_RandomSingles1_trigEff);
-        plotter2.region(2).plot(clusterEOne_RandomSingles1_trigEff);
         plotter2.region(2).plot(clusterEOne_RandomSingles1_trigEff);
         
     }