* twicking details in random generator

Examples/Example4/run.mac

-/run/beamOn 10
+/run/beamOn 1000

NPLib/Physics/NPReaction.cxx

@@ -87,7 +87,6 @@ ClassImp(Reaction)
     fshoot3=true;
     fshoot4=true;
     fUseExInGeant4=true;
-    RandGen=gRandom;
 
     fLabCrossSection=false; // flag if the provided cross-section is in the lab or not
 
@@ -156,7 +155,6 @@ Reaction::Reaction(string reaction){
 
   fshoot3=true;
   fshoot4=true;
-  RandGen=gRandom;
 
   fLabCrossSection=false;
 
@@ -357,7 +355,7 @@ double  Reaction::EnergyLabFromThetaLab(double ThetaLab){
  
   if(B>D) {
     ThetaLabMax = asin(sqrt(D/B));
-    if(RandGen->Rndm()<0.5) sign=-1;
+    if(gRandom->Rndm()<0.5) sign=-1;
   }
   if(ThetaLab>ThetaLabMax) return -1;
 

NPSimulation/Process/BeamReaction.cc

@@ -320,7 +320,7 @@ void NPS::BeamReaction::DoIt(const G4FastTrack& fastTrack,
         // Angles
         // Shoot and Set a Random ThetaCM
         m_Reaction.ShootRandomThetaCM();
-        double phi = RandFlat::shoot() * 2. * pi;
+        double phi = G4RandFlat::shoot() * 2. * pi;
 
         //////////////////////////////////////////////////
         /////  Momentum and angles from  kinematics  /////
@@ -453,8 +453,8 @@ void NPS::BeamReaction::DoIt(const G4FastTrack& fastTrack,
         // Shoot and Set a Random ThetaCM
         //m_QFS.ShootRandomThetaCM();
         //m_QFS.ShootRandomPhiCM();
-        double theta = RandFlat::shoot() *  pi;
-        double phi = RandFlat::shoot() * 2. * pi - pi; //rand in [-pi,pi]
+        double theta = G4RandFlat::shoot() *  pi;
+        double phi = G4RandFlat::shoot() * 2. * pi - pi; //rand in [-pi,pi]
 
         m_QFS.SetThetaCM(theta);
         m_QFS.SetPhiCM(phi);

NPSimulation/Simulation.cc

@@ -67,6 +67,7 @@ int main(int argc, char** argv){
     
     // initialize the state of the root and geant4 random generator
     if(OptionManager->GetRandomSeed()>0){
+      std::cout << " Seeds for random generators set to: " << OptionManager->GetRandomSeed() << std::endl;
       gRandom->SetSeed(OptionManager->GetRandomSeed()); 
       CLHEP::HepRandom::setTheSeed(OptionManager->GetRandomSeed(),3);
     }