X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=EMCAL%2Fjet_finder_ua1.F;h=4cdc1cdeda9ad8a096bf2b2d9aee2765a3d6c655;hp=87cbc404eed8dc05ddf51668c5a926b0b7de462c;hb=1a434afd893be48093e65e493c4b5d413864ffb7;hpb=ae70b41643709ab6463356830d24148d1d0f7f8d

diff --git a/EMCAL/jet_finder_ua1.F b/EMCAL/jet_finder_ua1.F
index 87cbc404eed..4cdc1cdeda9 100644
--- a/EMCAL/jet_finder_ua1.F
+++ b/EMCAL/jet_finder_ua1.F
@@ -18,6 +18,7 @@ c#include "calorSize.inc"
       real C_2PI 
       real etaCellSize 
       real phiCellSize 
+      real arg
       INTEGER NMAX,JMAX
       parameter(NMAX=30000,JMAX=100) ! 10-oct-2201
       integer ncell, ierror, mode, ncell_tot
@@ -105,8 +106,9 @@ c*-sum up total energy of all cells
 c*-Watch out!!! For mode=1, it can jump back to here for next iteration!!!
  999  continue
 c*-kill cells (flag=2) with Et below ET_MIN after background subtraction
-      call vzero(flag,ncell)
+cfca      call vzero(flag,ncell)
       do i=1, ncell
+         flag(i)=0
          if(etc(i)-et_ave .le. et_min) flag(i)=2
       enddo
       njet = 0
@@ -225,8 +227,18 @@ c*-sum up unused cells within required distance of given eta/phi
 ! 5-oct-2001 by PAI
             
 c*-reject cluster below minimum Ej_min
+c* protection (am)
             etas=eta+etas/ets
-            if(ets*cosh(etas/ets).lt.ej_min) then
+            arg = 0.
+            if (arg .ne. 0.) then
+               if (abs(etas/ets) .lt. 23.719) then
+                  arg = ets * cosh(etas/ets)
+               else
+                  arg = 1.e10
+               endif
+            endif
+            
+            if(arg .lt. ej_min) then
                do k=1,ncell
                   if(flag(k).le.0) flag(k)=0
                enddo