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 <TRandom.h>
  34
  35 #include "AliPythiaRndm.h"
  36
  37 TRandom * AliPythiaRndm::fgPythiaRandom=0;
  38
  39 ClassImp(AliPythiaRndm)
  40
  41
  42 //_______________________________________________________________________
  43 void AliPythiaRndm::SetPythiaRandom(TRandom *ran) {
  44   //
  45   // Sets the pointer to an existing random numbers generator
  46   //
  47   if(ran) fgPythiaRandom=ran;
  48   else fgPythiaRandom=gRandom;
  49 }
  50
  51 //_______________________________________________________________________
  52 TRandom * AliPythiaRndm::GetPythiaRandom() {
  53   //
  54   // Retrieves the pointer to the random numbers generator
  55   //
  56   return fgPythiaRandom;
  57 }
  58
  59 //_______________________________________________________________________
  60 #define pyr        pyr_
  61 #define pygauss    pygauss_
  62 #define pyrset     pyrset_
  63 #define pyrget     pyrget_
  64
  65 extern "C" {
  66     Double_t pyr(Int_t*)
  67     {
  68         // Wrapper to FINCTION PYR from PYTHIA
  69         // Uses static method to retrieve the pointer to the (C++) generator
  70         Double_t r;
  71         do r=AliPythiaRndm::GetPythiaRandom()->Rndm();
  72         while(0 >= r || r >= 1);
  73         return r;
  74     }
  75
  76     Double_t pygauss(Double_t x0, Double_t sig)
  77     {
  78         Double_t s = 2.;
  79         Double_t v1 = 0.;
  80         Double_t v2 = 0.;
  81
  82         while (s > 1.) {
  83             v1 = 2. * pyr(0) - 1.;
  84             v2 = 2. * pyr(0) - 1.;
  85             s = v1 * v1 + v2 * v2;
  86         }
  87         return v1 * TMath::Sqrt(-2. * TMath::Log(s) / s) * sig + x0;
  88     }
  89
  90     void pyrset(Int_t*,Int_t*) {}
  91     void pyrget(Int_t*,Int_t*) {}
  92 }