added docs

source/FinalWaveFunc.cpp

@@ -58,8 +58,7 @@ double finalWaveFunc::computeNorm(double rMax) {
 
 double DCAP::computeSquareWellDepth(shared_ptr<ParticleChannel>& channel, double energy, 
     double energyMin, double energyMax) {
-    double r = channel->getChannelRadius();
-    LogDerivDiff fobj(channel, energy, r);
+    LogDerivDiff fobj(channel, energy);
     TF1 * f = new TF1("f", fobj, energyMin-10, energyMax+10, 0);
     ROOT::Math::WrappedTF1 wf(*f);
     ROOT::Math::BrentRootFinder brf;

source/LogDerivDiff.cpp

@@ -23,8 +23,9 @@ using namespace ORM;
 using namespace DCAP;
 
 
-LogDerivDiff::LogDerivDiff(shared_ptr<ParticleChannel>& channel, double energy, double radius) 
-    : channel(channel), E(energy), r(radius) {
+LogDerivDiff::LogDerivDiff(shared_ptr<ParticleChannel>& channel, double energy) 
+    : channel(channel), E(energy) {
+    double r = channel->getChannelRadius();
     double W = channel->whittakerFunction(E);
     double dWdr = channel->whittakerFunctionDerivative(E);
     lhs = 1./(W/r) * (1./r*dWdr - W/(r*r));