Division by zero fixed

diff --git a/EVGEN/AliGenParam.cxx b/EVGEN/AliGenParam.cxx
index 14639ca..a6fa964 100644 (file)
--- a/EVGEN/AliGenParam.cxx
+++ b/EVGEN/AliGenParam.cxx
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.23  2000/10/02 21:28:06  fca
+Removal of useless dependecies via forward declarations
+
 Revision 1.22  2000/09/12 14:14:55  morsch
 Call fDecayer->ForceDecay() at the beginning of Generate().
 
@@ -340,6 +343,7 @@ void AliGenParam::Generate()
              wgtch=fChildWeight*fPtParaFunc(& ptd, &dummy);
          }
          xmt=sqrt(pt*pt+am*am);
+      if (TMath::Abs(ty)==1) ty=0;
          pl=xmt*ty/sqrt(1.-ty*ty);
          theta=TMath::ATan2(pt,pl);
 // Cut on theta

diff --git a/EVGEN/AliGenPythia.cxx b/EVGEN/AliGenPythia.cxx
index 9bd8a6c..50fbb57 100644 (file)
--- a/EVGEN/AliGenPythia.cxx
+++ b/EVGEN/AliGenPythia.cxx
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.24  2000/09/18 10:41:35  morsch
+Add possibility to use nuclear structure functions from PDF library V8.
+
 Revision 1.23  2000/09/14 14:05:40  morsch
 dito
 
@@ -420,11 +423,16 @@ Bool_t AliGenPythia::KinematicSelection(TParticle *particle)
 
 //
 // rapidity cut
-    Float_t y = 0.5*TMath::Log((e+pz)/(e-pz));
-    if (y > fYMax || y < fYMin)
-    {
+    if (e==pz) {
+      return kFALSE;
+    }
+    else {
+      Float_t y = 0.5*TMath::Log((e+pz)/(e-pz));
+      if (y > fYMax || y < fYMin)
+        {
 //     printf("\n failed y cut %f %f %f \n",y,fYMin,fYMax);
-       return kFALSE;
+          return kFALSE;
+        }
     }
 
 //