Minor corrections
[u/mrichter/AliRoot.git] / ZDC / AliGenZDC.cxx
index bcbaa6ef81b2ac0fc48ca202672e26a6d7eeee6f..13ff20efdb42ec00e0a19f763f69f42c301ee854 100644 (file)
 
 /*
 $Log$
+Revision 1.10  2001/09/26 16:00:47  coppedis
+Minor change
+
+Revision 1.8  2001/04/20 10:10:25  coppedis
+Minor changes
+
+Revision 1.7  2001/03/15 16:13:28  coppedis
+Code review
+
 Revision 1.6  2000/11/30 17:16:14  coppedis
 Changes suggested by fca
 
@@ -135,12 +144,6 @@ void AliGenZDC::Generate()
      fP[i] = pLab[i];
   }
   
-  if(fDebugOpt == 1){
-    printf("\n\n               Initial momentum :      \n");
-    for(i=0; i<=2; i++){
-       printf("        p[%d] = %f\n",i,pLab[i]);
-    }
-  }
   
   // Beam divergence and crossing angle
   if(fBeamCrossAngle!=0.) {
@@ -171,9 +174,7 @@ void AliGenZDC::Generate()
     TVector3 b(fP[0]/fP0, fP[1]/fP0, fP[2]/fP0);
     TLorentzVector pFermi(dddp[0], dddp[1], dddp[2], dddp0);
 
-    
     pFermi.Boost(b);
-
     for(i=0; i<=2; i++){
        fBoostP[i] = pFermi[i];
        fP[i] = pFermi[i];
@@ -241,7 +242,7 @@ void AliGenZDC::ExtractFermi(Int_t id, Double_t *fPp, Double_t *fProbintp,
   Float_t xx = gRandom->Rndm();
   assert ( id==kProton || id==kNeutron );
   if(id==kProton){
-    for(i=0; i<=200; i++){
+    for(i=1; i<=200; i++){
        if((xx>=fProbintp[i-1]) && (xx<fProbintp[i])) break;
        }
   }