From 918dee1c9c97f2d27382ff9ef33744bb01259362 Mon Sep 17 00:00:00 2001
From: "oliver.kirsebom@dal.ca" <oliver.kirsebom@dal.ca>
Date: Wed, 22 Apr 2020 08:32:35 -0300
Subject: [PATCH] added docs
---
source/FinalWaveFunc.cpp | 3 +--
source/LogDerivDiff.cpp | 5 +++--
2 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/source/FinalWaveFunc.cpp b/source/FinalWaveFunc.cpp
index 20573b6..c6bf35a 100644
--- a/source/FinalWaveFunc.cpp
+++ b/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;
diff --git a/source/LogDerivDiff.cpp b/source/LogDerivDiff.cpp
index e45fbcc..f5e0a3d 100644
--- a/source/LogDerivDiff.cpp
+++ b/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));
--
GitLab