Vertex smearing possible.

diff --git a/EVGEN/AliGenFixed.cxx b/EVGEN/AliGenFixed.cxx
index f359cbb917b62ef9af0777e6f078c519ec80f4d3..ad6f5b0bef96a7e1d66a8d9bfed2c2cbb476edab 100644 (file)
--- a/EVGEN/AliGenFixed.cxx
+++ b/EVGEN/AliGenFixed.cxx
@@ -60,29 +60,19 @@ void AliGenFixed::Generate()
   //
   Float_t polar[3]= {0,0,0};
   if(!fExplicit) {
-    fP[0] = fPMin*TMath::Cos(fPhiMin)*TMath::Sin(fThetaMin);
-    fP[1] = fPMin*TMath::Sin(fPhiMin)*TMath::Sin(fThetaMin);
-    fP[2] = fPMin*TMath::Cos(fThetaMin);
+    fP[0] = fPMin * TMath::Cos(fPhiMin) * TMath::Sin(fThetaMin);
+    fP[1] = fPMin * TMath::Sin(fPhiMin) * TMath::Sin(fThetaMin);
+    fP[2] = fPMin * TMath::Cos(fThetaMin);
   }
-  Int_t i, nt;
+  Int_t i, j, nt;
   //
   Float_t o[3];
-  o[0] = fOrigin[0];
-  o[1] = fOrigin[1];
-  o[2] = fOrigin[2];
-
-  printf("Origin %f %f %f \n", o[0], o[1], o[2]);
-  
+  if(fVertexSmear == kPerEvent) {
+      Vertex();
+      for (j = 0;j < 3; j++) o[j] = fVertex[j];
+  }
   
-  for(i=0;i<fNpart;i++) 
-    PushTrack(fTrackIt,-1,fIpart,fP,o ,polar,0,kPPrimary,nt);
+  for(i = 0; i < fNpart; i++) 
+    PushTrack(fTrackIt, -1, fIpart, fP, o , polar, 0, kPPrimary, nt);
 }
   
-//_____________________________________________________________________________
-void AliGenFixed::SetSigma(Float_t /*sx*/, Float_t /*sy*/, Float_t /*sz*/)
-{
-  //
-  // Set the interaction point sigma
-  //
-  printf("Vertex smearing not implemented for fixed generator\n");
-}

diff --git a/EVGEN/AliGenFixed.h b/EVGEN/AliGenFixed.h
index f0443cdd9fb9d7fb77287737592bd80e5619cab3..75a33479c3561954a2c431fee3f212ab53b03858 100644 (file)
--- a/EVGEN/AliGenFixed.h
+++ b/EVGEN/AliGenFixed.h
@@ -20,7 +20,6 @@ class AliGenFixed : public AliGenerator
   virtual ~AliGenFixed() {}
   virtual void Generate();
   virtual void Init() {}
-  virtual void SetSigma(Float_t sx, Float_t sy, Float_t sz);
   virtual void SetMomentum(Float_t pmom) {fPMin=pmom; fPMax=pmom; SetBit(kMomentumRange);}
   virtual void SetPhi(Float_t phi) {fPhiMin=phi*TMath::Pi()/180; fPhiMax=phi*TMath::Pi()/180; SetBit(kPhiRange);}
   virtual void SetTheta(Float_t theta) {fThetaMin=theta*TMath::Pi()/180; fThetaMax=theta*TMath::Pi()/180; SetBit(kThetaRange);}