From 1d568bc28fae0ba4c1cef92b8742ebbd6da623c8 Mon Sep 17 00:00:00 2001
From: morsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 8 Apr 2003 10:24:51 +0000
Subject: [PATCH] - Rapidity shift calculated in AliGenMC::Init() - Warning for
 depreciated method SetNclei. Prefer SetProjectile SetTarget.

---
 PYTHIA6/AliGenPythia.cxx | 23 +++++++++++++++++------
 PYTHIA6/AliGenPythia.h   |  1 +
 2 files changed, 18 insertions(+), 6 deletions(-)

diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx
index 838d9ebe374..8743ba83b7c 100644
--- a/PYTHIA6/AliGenPythia.cxx
+++ b/PYTHIA6/AliGenPythia.cxx
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.5  2003/04/04 08:12:56  morsch
+Boost for assymmetric systems added.
+
 Revision 1.4  2003/04/01 20:39:00  morsch
 Pycell configurable.
 
@@ -262,6 +265,7 @@ AliGenPythia::AliGenPythia()
   SetGammaPhiRange();
   SetGammaEtaRange();
   SetPtKick();
+  fSetNuclei = kFALSE;
 }
 
 AliGenPythia::AliGenPythia(Int_t npart)
@@ -304,6 +308,7 @@ AliGenPythia::AliGenPythia(Int_t npart)
     fCountMode = kCountAll;
     // Pycel
     SetPycellParameters();
+    fSetNuclei = kFALSE;
 }
 
 AliGenPythia::AliGenPythia(const AliGenPythia & Pythia)
@@ -463,7 +468,16 @@ void AliGenPythia::Init()
     fX2        = 0.;    
     fNev       = 0 ;
 //    
+//
+//
     AliGenMC::Init();
+//
+//
+//  
+    if (fSetNuclei) {
+	fDyBoost = 0;
+	Warning("Init","SetNuclei used. Use SetProjectile + SetTarget instead. fDyBoost has been reset to 0\n");
+    }
 }
 
 void AliGenPythia::Generate()
@@ -501,12 +515,7 @@ void AliGenPythia::Generate()
 	
 	fPythia->ImportParticles(fParticles,"All");
 
-	if (fZTarget != 0 && fAProjectile != 0) 
-	{
-	    Double_t dy    = - 0.5 * TMath::Log(Double_t(fZProjectile) * Double_t(fATarget) / 
-						(Double_t(fZTarget)    * Double_t(fAProjectile)));
-	    Boost(dy);
-	}
+	Boost();
 //
 //
 //
@@ -803,8 +812,10 @@ void AliGenPythia::AdjustWeights()
 void AliGenPythia::SetNuclei(Int_t a1, Int_t a2)
 {
 // Treat protons as inside nuclei with mass numbers a1 and a2  
+
     fAProjectile = a1;
     fATarget     = a2;
+    fSetNuclei   = kTRUE;
 }
 
 
diff --git a/PYTHIA6/AliGenPythia.h b/PYTHIA6/AliGenPythia.h
index 5acba60386f..4489dd180fb 100644
--- a/PYTHIA6/AliGenPythia.h
+++ b/PYTHIA6/AliGenPythia.h
@@ -173,6 +173,7 @@ class AliGenPythia : public AliGenMC
     Bool_t fFeedDownOpt;            // Option to set feed down from higher
                                     // quark families (e.g. b->c)
     Bool_t fFragmentation;          // Option to activate fragmentation by Pythia
+    Bool_t fSetNuclei;              // Flag indicating that SetNuclei has been called
     //
 
     CountMode_t fCountMode;         // Options for counting when the event will be finished.
-- 
2.31.1