minor changes

main.cpp

@@ -33,17 +33,17 @@ int main(int argc, char* argv[]) {
     Particle p{1, 1, Parity::P, 0.5_J};   // proton
     Particle B11{5, 11, Parity::P, 1.5_J}; // B11 g.s. 3/2+
 
-    double beamEnergy = 2000;
+    double beamEnergy = 2000; //keV
+    double r0 = 1.36; //reduced channel radius in fm
     
     vector<double> excEnergy{9640, 10844, 11820, 12710};
     vector<Spin> spin{3_J, 1_J, 2_J, 1_J};
     vector<OrbitalAngularMomentum> lOrb{2_L, 2_L, 2_L, 1_L};
-    vector<double> S{0.3, 0.10, 0.15, 1}; //spectroscopic factors relative to 12.71
-    for (auto& ss : S) ss *= 0.86;
+    vector<double> S{0.30, 0.23, 0.11, 0.86}; //spectroscopic factors relative to 12.71
 
 
     cout << Form("\nDC cross section @ %.2f MeV\n", beamEnergy/1e3) << endl;
-    cout << " Ex (MeV)   l    S    cs (ub)" << endl;
+    cout << " Ex (MeV)   li   lf    S    cs (ub)" << endl;
     cout << "--------------------------------" << endl;
 
     for (size_t i=0; i<excEnergy.size(); i++) {
@@ -51,19 +51,22 @@ int main(int argc, char* argv[]) {
         vector<shared_ptr<ParticleChannel>> channels;
         for (size_t l=0; l<6; l++) {
             auto c = make_shared<ParticleChannel>(p, B11, spin[i], OrbitalAngularMomentum(l));
-            c->setR0(1.36);
+            c->setR0(r0);
             channels.push_back(c);
         }
 
         auto& final = channels[lOrb[i].getValue()];
 
         double sigmaTot = 0;
+        string li = "";
         for (size_t l=0; l<6; l++) {
             double sigma = directCaptureCrossSec(channels[l], final, beamEnergy, excEnergy[i], S[i]);
+            if (sigma > 0) li += Form("%d,", int(l));
             sigmaTot += sigma;
         } 
 
-        cout << Form("  %.2f     %.0f   %.2f   %.2f", excEnergy[i]/1e3, lOrb[i].getValue(), S[i], sigmaTot) << endl;
+        li.pop_back();
+        cout << Form("  %.2f     %s   %.0f   %.2f   %.2f", excEnergy[i]/1e3, li.c_str(), lOrb[i].getValue(), S[i], sigmaTot) << endl;
     }
     cout << "--------------------------------\n" << endl;