Access to common block ALLDLT containing primary ionisation electron information.

diff --git a/TFluka/Falldlt.h b/TFluka/Falldlt.h
index 71c90fcc6815c734f04971c46ec0473f48d7fb40..dc4b89a66fd322f414d8aa56aa31bba9514a1103 100644 (file)
--- a/TFluka/Falldlt.h
+++ b/TFluka/Falldlt.h
@@ -51,7 +51,7 @@ extern "C" {
 //      SAVE / ALLDLT /
     const Int_t mxalld = 5000;
     typedef struct {
-       Double_t talldl[kmxalld];
+       Double_t talldl[mxalld];
        Double_t tallmn;
        Int_t    nalldl;
        Int_t    lalldl;

diff --git a/TFluka/TFluka.cxx b/TFluka/TFluka.cxx
index 30c8e0080bd342edf23da9f79b266bab63f96410..30dea6944cd08bef4b17beae0a56a0ecea35009b 100644 (file)
--- a/TFluka/TFluka.cxx
+++ b/TFluka/TFluka.cxx
@@ -48,6 +48,7 @@
 #include "Fstepsz.h"       //(STEPSZ) fluka common
 #include "Fopphst.h"       //(OPPHST) fluka common
 #include "Fltclcm.h"       //(LTCLCM) fluka common
+#include "Falldlt.h"       //(ALLDLT) fluka common
 
 #include "TVirtualMC.h"
 #include "TMCProcess.h"
@@ -2162,3 +2163,22 @@ void TFluka::AddParticlesToPdgDataBase() const
                      0,6,"Ion",kion+20030);
 }
 
+  //
+  // Info about primary ionization electrons
+  Int_t TFluka::GetNPrimaryElectrons()
+{
+    // Get number of primary electrons
+    return ALLDLT.nalldl;
+}
+
+Double_t GetPrimaryElectronKineticEnergy(Int_t i)
+{
+    // Returns kinetic energy of primary electron i
+    if (i >= 0 && i < ALLDLT.nalldl) {
+       return ALLDLT.talldl[i];
+    } else {
+       Warning("GetPrimaryElectronKineticEnergy", 
+               "Primary electron index out of range %d %d \n", 
+               i, ALLDLT.nalldl);
+    }
+}

diff --git a/TFluka/TFluka.h b/TFluka/TFluka.h
index d5aed8ce075dfb54c99a8b4bdab45c32bec2f928..e2d9a27bdc610c577db0779d102e7bf3ac592328 100644 (file)
--- a/TFluka/TFluka.h
+++ b/TFluka/TFluka.h
@@ -386,6 +386,11 @@ class TFluka : public TVirtualMC {
   private:
   void PrintHeader();
   void AddParticlesToPdgDataBase() const;
+  //
+  // Info about primary ionization electrons
+  Int_t    GetNPrimaryElectrons();
+  Double_t GetPrimaryElectronKineticEnergy(Int_t i);
+  //
   TFluka(const TFluka &mc): TVirtualMC(mc) {;}
   TFluka & operator=(const TFluka &) {return (*this);}