--- trunk/src/brains/SimInfo.cpp	2006/07/03 13:18:43	998
+++ trunk/src/brains/SimInfo.cpp	2006/09/21 18:25:17	1045
@@ -53,6 +53,7 @@
 #include "brains/SimInfo.hpp"
 #include "math/Vector3.hpp"
 #include "primitives/Molecule.hpp"
+#include "primitives/StuntDouble.hpp"
 #include "UseTheForce/fCutoffPolicy.h"
 #include "UseTheForce/DarkSide/fElectrostaticSummationMethod.h"
 #include "UseTheForce/DarkSide/fElectrostaticScreeningMethod.h"
@@ -599,6 +600,8 @@ namespace oopse {
     //setup fortran force field
     /** @deprecate */    
     int isError = 0;
+    
+    setupCutoff();
     
     setupElectrostaticSummationMethod( isError );
     setupSwitchingFunction();
@@ -610,9 +613,6 @@ namespace oopse {
       painCave.isFatal = 1;
       simError();
     }
-  
-    
-    setupCutoff();
 
     calcNdf();
     calcNdfRaw();
@@ -1073,9 +1073,8 @@ namespace oopse {
     int sm = UNDAMPED;
     RealType alphaVal;
     RealType dielectric;
-
+    
     errorOut = isError;
-    alphaVal = simParams_->getDampingAlpha();
     dielectric = simParams_->getDielectric();
 
     if (simParams_->haveElectrostaticSummationMethod()) {
@@ -1121,10 +1120,16 @@ namespace oopse {
 	if (myScreen == "DAMPED") {
 	  sm = DAMPED;
 	  if (!simParams_->haveDampingAlpha()) {
-	    //throw error
+	    // first set a cutoff dependent alpha value
+	    // we assume alpha depends linearly with rcut from 0 to 20.5 ang
+	    alphaVal = 0.5125 - rcut_* 0.025;
+	    // for values rcut > 20.5, alpha is zero
+	    if (alphaVal < 0) alphaVal = 0;
+
+	    // throw warning
 	    sprintf( painCave.errMsg,
 		     "SimInfo warning: dampingAlpha was not specified in the input file.\n"
-                     "\tA default value of %f (1/ang) will be used.\n", alphaVal);
+                     "\tA default value of %f (1/ang) will be used for the cutoff of\n\t%f (ang).\n", alphaVal, rcut_);
 	    painCave.isFatal = 0;
 	    simError();
 	  }
@@ -1436,6 +1441,24 @@ namespace oopse {
       return angularMomentum;
    }
    
-   
+  StuntDouble* SimInfo::getIOIndexToIntegrableObject(int index) {
+    return IOIndexToIntegrableObject.at(index);
+  }
+  
+  void SimInfo::setIOIndexToIntegrableObject(const std::vector<StuntDouble*>& v) {
+    IOIndexToIntegrableObject= v;
+  }
+
+/*
+   void SimInfo::setStuntDoubleFromGlobalIndex(std::vector<StuntDouble*> v) {
+      assert( v.size() == nAtoms_ + nRigidBodies_);
+      sdByGlobalIndex_ = v;
+    }
+
+    StuntDouble* SimInfo::getStuntDoubleFromGlobalIndex(int index) {
+      //assert(index < nAtoms_ + nRigidBodies_);
+      return sdByGlobalIndex_.at(index);
+    }   
+*/   
 }//end namespace oopse