Implemented starting configuration.

[phys/u1casc.git] / u1casc-flux / sim.hpp

diff --git a/u1casc-flux/sim.hpp b/u1casc-flux/sim.hpp
index 7fd59aff106cd447785b33d106fd6b1fad5bc001..496a28db15a944f9c2c9c9192b7ab2be9eab3a9c 100644 (file)
--- a/u1casc-flux/sim.hpp
+++ b/u1casc-flux/sim.hpp
@@ -24,10 +24,11 @@ public:
   double WF(int nom, int denom, int iflav);
   int iPlaq(int i, int j);
   double I(int p);
   double WF(int nom, int denom, int iflav);
   int iPlaq(int i, int j);
   double I(int p);

-
+  void _newParas();
+  

 private:
 private:

+  void _resetConfig();

   void _makeSweep();
   void _makeSweep();

-  void _newParas();

   gsl_rng* rangsl;
   int lineSweep();
   int lineUpdate(int x, int vardir);
   gsl_rng* rangsl;
   int lineSweep();
   int lineUpdate(int x, int vardir);


@@ -92,7 +93,7 @@ sim::sim(o815 *_O815) : o815::sim( _O815,
 }
 
 void sim::_makeSweep() 
 }
 
 void sim::_makeSweep() 

-{  
+{ 

   lineSweep();
   
   for( int x=0; x<lsize4; x++) {
   lineSweep();
   
   for( int x=0; x<lsize4; x++) {


@@ -118,8 +119,10 @@ void sim::_newParas() {
   beta = (*O815->paraQ)["beta"];
   Mu[0] = (*O815->paraQ)["muone"];
   Mu[1] = (*O815->paraQ)["mutwo"];
   beta = (*O815->paraQ)["beta"];
   Mu[0] = (*O815->paraQ)["muone"];
   Mu[1] = (*O815->paraQ)["mutwo"];

+}

 
 

-  memset(confMem, 0, sizeof(int)*lsize4*(2*4+2*4+6));
+void sim::_resetConfig() {
+  memset(confMem, 0, confmemSize);

 }
 
 int sim::lineSweep()
 }
 
 int sim::lineSweep()