diff --git a/source/processes/electromagnetic/dna/models/src/G4DNAPTBElasticModel.cc b/source/processes/electromagnetic/dna/models/src/G4DNAPTBElasticModel.cc
index f8ab025a499..297039a1529 100644
--- a/source/processes/electromagnetic/dna/models/src/G4DNAPTBElasticModel.cc
+++ b/source/processes/electromagnetic/dna/models/src/G4DNAPTBElasticModel.cc
@@ -368,7 +368,7 @@ void G4DNAPTBElasticModel::SampleSecondaries(std::vector<G4DynamicParticle*>* /*
   }
 
   G4double electronEnergy0 = aDynamicElectron->GetKineticEnergy();
-  const std::size_t& materialID = couple->GetIndex();
+  const std::size_t& materialID = couple->GetMaterial()->GetIndex();
   auto p = aDynamicElectron->GetParticleDefinition();
 
   // set killBelowEnergy value for material
diff --git a/source/processes/electromagnetic/dna/models/src/G4DNAPTBExcitationModel.cc b/source/processes/electromagnetic/dna/models/src/G4DNAPTBExcitationModel.cc
index 33e6e4c4c49..63360cd4f2d 100644
--- a/source/processes/electromagnetic/dna/models/src/G4DNAPTBExcitationModel.cc
+++ b/source/processes/electromagnetic/dna/models/src/G4DNAPTBExcitationModel.cc
@@ -255,7 +255,7 @@ void G4DNAPTBExcitationModel::SampleSecondaries(std::vector<G4DynamicParticle*>*
                                                 const G4DynamicParticle* aDynamicParticle,
                                                 G4double /*tmin*/, G4double /*tmax*/)
 {
-  const std::size_t& materialID = (std::size_t)couple->GetIndex();
+  const std::size_t& materialID = (std::size_t)couple->GetMaterial()->GetIndex();
 
   // Get the incident particle kinetic energy
   G4double k = aDynamicParticle->GetKineticEnergy();
diff --git a/source/processes/electromagnetic/dna/utils/src/G4DNAPTBIonisationStructure.cc b/source/processes/electromagnetic/dna/utils/src/G4DNAPTBIonisationStructure.cc
index ce961ba61ae..69338098957 100644
--- a/source/processes/electromagnetic/dna/utils/src/G4DNAPTBIonisationStructure.cc
+++ b/source/processes/electromagnetic/dna/utils/src/G4DNAPTBIonisationStructure.cc
@@ -52,6 +52,7 @@ G4DNAPTBIonisationStructure::G4DNAPTBIonisationStructure()
     energyConstant[index].push_back(17.07 * eV);
     energyConstant[index].push_back(21.00 * eV);
     energyConstant[index].push_back(41.72 * eV);
+    energyConstant[index].push_back(409.9*eV);
   }
 
   // MPietrzak