Fixed division by zero
authorpostrow <postrow@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 14 Aug 2009 08:43:23 +0000 (08:43 +0000)
committerpostrow <postrow@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 14 Aug 2009 08:43:23 +0000 (08:43 +0000)
EPOS/TEpos.cxx
EPOS/epos167/epos-tim-155.f

index 8f7c0bb..e993913 100644 (file)
@@ -105,11 +105,11 @@ void TEpos::Initialize() {
                IniModel(appli.model);
        ebin.nrebin = 1;
        ainit();
+       aseed(2);
 }
 
 void TEpos::GenerateEvent() {
-       cseed.seedj = gRandom->Rndm() * 1e10;
-       aseed(2);
+//     cseed.seedj = gRandom->Rndm() * 1e10;
        Int_t n = 1;
        evgen(n);
 }
index 29d6d7d..a439712 100644 (file)
@@ -426,7 +426,7 @@ c.....select flavor and z-value .....................................
       if(alfm.lt.rangen()*log(q2*z*(1.-z)/qcdlam)) goto 390
 
       if(ij(1).ne.j1.or.ij(2).eq.0)then
-        if(E.lt.sqrt(q2))goto 390
+        if(E.le.sqrt(q2))goto 390
         pzz=sqrt((E-sqrt(q2))*(E+sqrt(q2)))
         pt2=(E**2*(z*(1.-z)*q2-z*amm2(idflb)-(1.-z)*amm2(idfla))
      $       -.25*(amm2(idflb)-amm2(idfla)-q2)**2+q2*amm2(idfla))/pzz**2