*Adding support for beam excitation energy

 - optional token, keeping retrocompatibility
 - allow to use "beam" at rest as source for gamma cascade

Inputs/EventGenerator/24Nadp.reaction

@@ -2,6 +2,7 @@
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Beam
   Particle= 24Na
+  ExcitationEnergy= 0
   Energy= 192
   SigmaEnergy= 1
   SigmaThetaX= 0.01 

NPLib/Physics/NPBeam.cxx

@@ -82,6 +82,7 @@ Beam::Beam(){
 Beam::Beam(string isotope){
   SetUp(isotope);
   fEnergy = 0;
+  fExcitationEnergy = 0; 
   fSigmaEnergy = -1 ;
   fMeanX = 0 ;
   fMeanY = 0 ;
@@ -183,6 +184,12 @@ void Beam::ReadConfigurationFile(string Path){
         if(fVerboseLevel==1) cout << "Beam Particle: " << GetName() << endl;
       }
       
+      else if (DataBuffer == "ExcitationEnergy=") {
+        ReactionFile >> DataBuffer;
+        fExcitationEnergy = atof(DataBuffer.c_str()) * MeV;
+        if(fVerboseLevel==1) cout << "Excitation Energy: " << fExcitationEnergy << " MeV" << endl;
+      }
+      
       else if (DataBuffer == "Energy=") {
         check_Energy = true ;
         ReactionFile >> DataBuffer;

NPLib/Physics/NPBeam.h

@@ -57,6 +57,7 @@ namespace NPL{
   private:
     //Nucleus* fBeamNucleus;
     double fEnergy;
+    double fExcitationEnergy;
     double fSigmaEnergy;
     double fMeanX;
     double fMeanY;
@@ -94,6 +95,7 @@ namespace NPL{
     // Get
     // Nucleus*  GetNucleus     () const {return fBeamNucleus;}
     double    GetEnergy      () const {return fEnergy;}
+    double    GetExcitationEnergy() const {return fExcitationEnergy;}
     double    GetSigmaEnergy () const {return fSigmaEnergy;}
     double    GetMeanX       () const {return fMeanX;}
     double    GetMeanY       () const {return fMeanY;}

NPSimulation/src/EventGeneratorBeam.cc

@@ -75,7 +75,7 @@ void EventGeneratorBeam::GenerateEvent(G4Event* anEvent){
 
   if( anEvent->GetEventID()==0){
     // Define the particle to be shoot
-    m_particle = G4ParticleTable::GetParticleTable()->GetIon(m_Beam->GetZ(), m_Beam->GetA() , 0.);
+    m_particle = G4ParticleTable::GetParticleTable()->GetIon(m_Beam->GetZ(), m_Beam->GetA() ,m_Beam->GetExcitationEnergy());
   }
   
   ///////////////////////////////////////////////////////////////////////