Modified absorber correction. Added function FieldCorrection() to account

[u/mrichter/AliRoot.git] / EVGEN / AliPythia.cxx

diff --git a/EVGEN/AliPythia.cxx b/EVGEN/AliPythia.cxx
index 9b50d8f88fb3c86f78feaaccccac6a02d1bf3824..cf3c951dd5329c2feec6520e8ec6ba2de2aee65f 100644 (file)
--- a/EVGEN/AliPythia.cxx
+++ b/EVGEN/AliPythia.cxx
@@ -15,6 +15,12 @@
 
 /*
 $Log$
+Revision 1.23  2002/05/06 07:17:29  morsch
+Pyr gives random number r in interval 0 < r < 1.
+
+Revision 1.22  2002/04/26 10:28:48  morsch
+Option kPyBeautyPbMNR added (N. Carrer).
+
 Revision 1.21  2002/03/25 14:46:16  morsch
 Case  kPyD0PbMNR added (N. Carrer).
 
@@ -182,7 +188,21 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
        SetPARP(82,3.47);   // set value pT_0  for turn-off of the cross section of                  
                             // multiple interaction at a reference energy = 14000 GeV
        SetPARP(89,14000.); // reference energy for the above parameter
-       SetPARP(90,0.174);  // set exponent for energy dependence of pT_0 
+       SetPARP(90,0.174);  // set exponent for energy dependence of pT_0
+    case kPyMbNonDiffr:
+// Minimum Bias pp-Collisions
+//
+//   
+//      select Pythia min. bias model
+       SetMSEL(0);
+       SetMSUB(95,1);      // low pt production
+       SetMSTP(81,1);      // multiple interactions switched on
+       SetMSTP(82,3);      // model with varying impact param. & a single Gaussian
+       SetPARP(82,3.47);   // set value pT_0  for turn-off of the cross section of                  
+                            // multiple interaction at a reference energy = 14000 GeV
+       SetPARP(89,14000.); // reference energy for the above parameter
+       SetPARP(90,0.174);  // set exponent for energy dependence of pT_0
+ 
        break;
     case kPyJets:
        SetMSEL(1);
@@ -378,7 +398,12 @@ void  AliPythia::SetDecayTable()
 #endif
 
 extern "C" {
-  Double_t pyr(Int_t*) {return sRandom->Rndm();}
+  Double_t pyr(Int_t*) 
+{
+      Float_t r;
+      do r=sRandom->Rndm(); while(0 >= r || r >= 1);
+      return r;
+}
   void pyrset(Int_t*,Int_t*) {}
   void pyrget(Int_t*,Int_t*) {}
 }