diff --git a/scripts/checkPulseHam_diff.py b/scripts/checkPulseHam_diff.py
index b2ae650..d639f11 100644
--- a/scripts/checkPulseHam_diff.py
+++ b/scripts/checkPulseHam_diff.py
@@ -6,8 +6,8 @@
 
 macros = '/home/hep/matzeni/CCE/macros/CCEScan/'
 
-filelumi = TFile('../CCE/data/db/lumi.root')
-file_v_depl = TFile('../CCE/data/db/vdepl.root')
+filelumi = TFile('../data/db/lumi.root')
+file_v_depl = TFile('../data/db/vdepl.root')
 fileHam = TFile('/disk/data12/lhcb/STAgeingData/data/ST/Aging/Hamburg.root')
 fileData = TFile('/disk/data12/lhcb/STAgeingData/data/ST/Aging/CCEScan.root')
 
@@ -71,12 +71,15 @@
     #Data
     vD_time =[]
     vData=[]
+    vData_err=[]
 
     for ifill in fills:
 
 
         #  vData =np.concatenate((vData, np.array(fileData.Get('TT/TTaU/'+isect+'/'+str(int(ifill))+'/v_volt_val7'))),axis=0)    
         vData = np.append( np.array(fileData.Get('TT/TTaU/'+isect+'/'+ifill+'/v_volt_val7'))[0],vData)
+        vData_err = np.append( np.array(fileData.Get('TT/TTaU/'+isect+'/'+ifill+'/v_volt_val7'))[1],vData_err)
+        
         print(df['Time'][df.Fill == int(ifill)].values[0])
         vD_time = np.append( df['Time'][df.Fill == int(ifill)].values[0],vD_time)
 
@@ -89,9 +92,22 @@
     #####################################
 
     vH_inter = np.interp(vD_time, vH_time, vH_V)
+    print 'vH_time: ',vH_time
+    print 'vD_time: ',vD_time
     print 'vH_V: ',vH_V
+    print 'vData: ',vData
+    print 'vData_err: ',vData_err
     print 'vH_inter: ',vH_inter
-    print 'Difference" ',vData-vH_inter
+    vDiff = np.absolute(vData-vH_inter)
+    print 'Difference: ',vDiff
+
+    for i in xrange(len(vD_time)):
+        print(i)
+        if (vDiff[i] >= 3*vData_err[i] ):
+            print 'WARNING! Check Sector: {} for Fill: {} '.format(isect,df['Fill'][df.Time==vD_time[i]])
+
+
+
     plt.plot(vH_time,vH_V, 'r--',vD_time,vH_inter, 'bs',vD_time,vData, 'gs')
     plt.show()
     raw_input('>')