Division by zero fixed
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 18 Oct 2000 19:11:27 +0000 (19:11 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 18 Oct 2000 19:11:27 +0000 (19:11 +0000)
EVGEN/AliGenParam.cxx
EVGEN/AliGenPythia.cxx

index 14639ca..a6fa964 100644 (file)
@@ -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
index 9bd8a6c..50fbb57 100644 (file)
@@ -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;
+        }
     }
 
 //