Changeset 1340 for trunk/source/processes/electromagnetic/lowenergy/src/G4DNAExcitation.cc

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/electromagnetic/lowenergy/src/G4DNAExcitation.cc (modified) (4 diffs)

trunk/source/processes/electromagnetic/lowenergy/src/G4DNAExcitation.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4DNAExcitation.cc,v 1.4 2010/03/27 11:32:41 sincerti Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4DNAExcitation.cc,v 1.7 2010/10/08 08:53:17 sincerti Exp $
+// GEANT4 tag $Name: emlowen-V09-03-54 $
 
 #include "G4DNAExcitation.hh"
@@ -57 +57 @@
        &p == G4Electron::Electron() 
     || &p == G4Proton::ProtonDefinition()
+    || &p == instance->GetIon("hydrogen")
     || &p == instance->GetIon("alpha++")
     || &p == instance->GetIon("alpha+")
@@ -77 +78 @@
     {
 
-      // First model
-
+      // Emfietzoglou model
+/*
       if(!Model()) SetModel(new G4DNAEmfietzoglouExcitationModel);
       Model()->SetLowEnergyLimit(8.23*eV);
       Model()->SetHighEnergyLimit(10*MeV);
+*/
+      // Born model
 
-      // Alternative model
-/*
       if(!Model()) SetModel(new G4DNABornExcitationModel);
       Model()->SetLowEnergyLimit(9*eV);
       Model()->SetHighEnergyLimit(1*MeV);
-*/
+
       AddEmModel(1, Model());   
     }
@@ -106 +107 @@
     }
 
+    if(name == "hydrogen")
+    {
+      if(!Model()) SetModel(new G4DNAMillerGreenExcitationModel);
+      Model()->SetLowEnergyLimit(10*eV);
+      Model()->SetHighEnergyLimit(500*keV);
+   
+      AddEmModel(1, Model());   
+    }
+
+
     if( name == "alpha" || name == "alpha+" || name == "helium" )
     {
       if(!Model()) SetModel(new G4DNAMillerGreenExcitationModel);
       Model()->SetLowEnergyLimit(1*keV);
-      Model()->SetHighEnergyLimit(10*MeV);
+      Model()->SetHighEnergyLimit(400*MeV);
 
       AddEmModel(1, Model());