Introducing Id
[u/mrichter/AliRoot.git] / ZDC / AliGenZDC.cxx
index bcbaa6e..96b9f05 100644 (file)
 
 /*
 $Log$
+Revision 1.11  2001/11/11 17:27:50  hristov
+Minor corrections
+
+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
 
@@ -45,12 +57,12 @@ Introduction of the Copyright and cvs Log
 #include <TLorentzVector.h>
 #include <TVector3.h>
 #include "TDatabasePDG.h"
+#include "TMCProcess.h"
 
 #include "AliGenZDC.h"
 #include "AliConst.h"
 #include "AliPDG.h"
 #include "AliRun.h"
-#include "AliMCProcess.h"
  
 ClassImp(AliGenZDC)
  
@@ -135,12 +147,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 +177,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 +245,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;
        }
   }