patch for random seed
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 24 Aug 2012 19:56:26 +0000 (19:56 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 24 Aug 2012 19:56:26 +0000 (19:56 +0000)
<http://savannah.cern.ch/task/?32889>

 Fiorella Fionda

TEvtGen/EvtGenModels/EvtPHOTOS.cxx
TEvtGen/PHOTOS/phocin.F
TEvtGen/PHOTOS/phoini.F

index 807eba3..bf9b1a4 100644 (file)
@@ -28,6 +28,7 @@
 #include "EvtGenBase/EvtPDL.hh"
 #include "EvtGenModels/EvtPHOTOS.hh"
 #include "EvtGenBase/EvtReport.hh"
+#include "TRandom.h"
 #include <stdlib.h>
 
 extern "C" void begevtgenstorex_(int *,int *,int *,int *,
@@ -46,7 +47,7 @@ extern "C" void heplst_(int *);
 
 extern "C" void photos_(int *);
 
-extern "C" void phoini_();
+extern "C" void phoini_(int *, int *);
 
 
 EvtPHOTOS::EvtPHOTOS(std::string photontype){
@@ -71,7 +72,9 @@ void EvtPHOTOS::doRadCorr( EvtParticle *p){
 
   if (first) {
     first=0;
-    phoini_();
+    int iseed1=(int)gRandom->Integer(31327);
+    int iseed2=(int)gRandom->Integer(30080);
+    phoini_(&iseed1,&iseed2);
   }
 
   double mpho=EvtPDL::getMeanMass(GAMM);
index 5ea1789..d0385d4 100644 (file)
@@ -1,4 +1,4 @@
-      SUBROUTINE PHOCIN
+      SUBROUTINE PHOCIN(ISEEDF,ISEEDS)
 C.----------------------------------------------------------------------
 C.
 C.    PHOTOS:   PHOton Common INitialisation
@@ -15,6 +15,7 @@ C.                Z. Was                          Last Update: 10/08/93
 C.
 C.----------------------------------------------------------------------
       IMPLICIT NONE
+      INTEGER ISEEDF, ISEEDS
       INTEGER NMXHEP
       PARAMETER (NMXHEP=10000)
       LOGICAL QEDRAD
@@ -60,8 +61,11 @@ C--   Define some constants
       TWOPI=6.28318530717958648D0
 C--
 C--   Default seeds Marsaglia and Zaman random number generator
-      ISEED(1)=1802
-      ISEED(2)=9373
+C      ISEED(1)=1802
+C      ISEED(2)=9373
+C      updated on 23.08.2012 to introduce proper random number sequence
+       ISEED(1)=ISEEDF
+       ISEED(2)=ISEEDS
 C--
 C--   Iitialization for extra options
 C--   (1)
index 6a52602..fef1fb8 100644 (file)
@@ -86,7 +86,8 @@ C.   PHOPRO   3)               4      var. for photon rep. (in branch)
 C.   PHOCMS   2)               3      parameters of boost to branch CMS
 C.   PHNUM    4)               1      event number from outside         
 C.----------------------------------------------------------------------
-      SUBROUTINE PHOINI
+      SUBROUTINE PHOINI(ISEEDF,ISEEDS)
+
 C.----------------------------------------------------------------------
 C.
 C.    PHOTOS:   PHOton radiation in decays INItialisation
@@ -104,6 +105,7 @@ C.                                                Last Update: 12/04/90
 C.
 C.----------------------------------------------------------------------
       IMPLICIT NONE
+      INTEGER ISEEDF, ISEEDS
       INTEGER INIT
       SAVE INIT
       DATA INIT/ 0/
@@ -113,7 +115,7 @@ C--   Return if already initialized...
       INIT=1
 C--
 C--   Preset parameters in PHOTOS commons
-      CALL PHOCIN
+      CALL PHOCIN(ISEEDF,ISEEDS)
 C--
 C--   Print info
       CALL PHOINF