Print

---

Author: [log in to unmask]
Date: Mon Jan 23 11:07:45 2017
New Revision: 4687

Log:
smear ecal hit time in mc matches data

Modified:
    java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/HitTMCSmearDriver.java

Modified: java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/HitTMCSmearDriver.java
 =============================================================================
--- java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/HitTMCSmearDriver.java	(original)
+++ java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/HitTMCSmearDriver.java	Mon Jan 23 11:07:45 2017
@@ -30,9 +30,10 @@
     
    
     // Time resolution as derived for 2016 data
+    // Factor is reduced by 2.0 through MC testing to match data
     private static double calcSmear(double energy, double time){
         Random r = new Random();
-        double sigT = r.nextGaussian()*Math.sqrt(Math.pow(0.188/energy, 2) + Math.pow(0.152, 2));
+        double sigT = r.nextGaussian()*Math.sqrt(Math.pow(0.188/energy, 2) + Math.pow(0.152, 2))/2.0;
         return time + sigT;
     }
             
@@ -42,15 +43,14 @@
         if (event.hasCollection(CalorimeterHit.class, inputHitCollection)){
             //Get the hits in the event       
             List<CalorimeterHit> hits = event.get(CalorimeterHit.class, inputHitCollection);
-            
             for (CalorimeterHit iHit : hits){
                 double oldT = iHit.getTime();
                 double energy = iHit.getRawEnergy();          
                 double newT = calcSmear(energy, oldT);
-                
+
                 ((BaseCalorimeterHit) iHit).setTime(newT);
 
             }    
         } 
     }
-}
+}