Inserting TMath.h where required by the new version of ROOT
[u/mrichter/AliRoot.git] / PYTHIA6 / AliPythiaRndm.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id$ */
17
18 //-----------------------------------------------------------------------------
19 //   Class: AliPythiaRndm
20 //   Responsibilities: Interface to Root random number generator 
21 //                     from Fortran (re-implements FINCTION PYR from PYTHIA)
22 //                     Very similar to AliHijingRndm
23 //   Collaborators: AliPythia and AliGenPythia classes
24 //   Example:
25 //
26 //   root> AliPythia::Instance();
27 //   root> AliPythiaRndm::SetPythiaRandom(new TRandom3()); 
28 //   root> AliPythiaRndm::GetPythiaRandom()->SetSeed(0);
29 //   root> cout<<"Seed "<< AliPythiaRndm::GetPythiaRandom()->GetSeed() <<endl;
30 //
31 //-----------------------------------------------------------------------------
32
33 #include <TMath.h>
34 #include <TRandom.h>
35
36 #include "AliPythiaRndm.h"
37
38 TRandom * AliPythiaRndm::fgPythiaRandom=0;
39
40 ClassImp(AliPythiaRndm)
41
42
43 //_______________________________________________________________________
44 void AliPythiaRndm::SetPythiaRandom(TRandom *ran) {
45   //
46   // Sets the pointer to an existing random numbers generator
47   //
48   if(ran) fgPythiaRandom=ran;
49   else fgPythiaRandom=gRandom;
50 }
51
52 //_______________________________________________________________________
53 TRandom * AliPythiaRndm::GetPythiaRandom() {
54   //
55   // Retrieves the pointer to the random numbers generator
56   //
57   return fgPythiaRandom;
58 }
59
60 //_______________________________________________________________________
61 #define pyr        pyr_
62 #define pygauss    pygauss_
63 #define pyrset     pyrset_
64 #define pyrget     pyrget_
65
66 extern "C" {
67     Double_t pyr(Int_t*) 
68     {
69         // Wrapper to FINCTION PYR from PYTHIA
70         // Uses static method to retrieve the pointer to the (C++) generator
71         Double_t r;
72         do r=AliPythiaRndm::GetPythiaRandom()->Rndm();
73         while(0 >= r || r >= 1);
74         return r;
75     }
76     
77     Double_t pygauss(Double_t x0, Double_t sig)
78     {
79         Double_t s = 2.;
80         Double_t v1 = 0.;
81         Double_t v2 = 0.;
82         
83         while (s > 1.) {
84             v1 = 2. * pyr(0) - 1.;
85             v2 = 2. * pyr(0) - 1.;
86             s = v1 * v1 + v2 * v2;
87         }
88         return v1 * TMath::Sqrt(-2. * TMath::Log(s) / s) * sig + x0;
89     }
90
91     void pyrset(Int_t*,Int_t*) {}
92     void pyrget(Int_t*,Int_t*) {}
93 }