simhadteam team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 10
committer: Daniel Cussol <cussol@xxxxxxxx>
branch nick: G4LPC
timestamp: Fri 2010-04-30 14:50:38 +0200
message:
  Ajout de materiaux "Eau_ioxx" permettant de
modified:
  src/G4LPCDetectorConstruction.cc
  src/G4LPCRunAction.cc


--
lp:simulationshadrontherapie/g4lpc
https://code.launchpad.net/~simhadteam/simulationshadrontherapie/G4LPC

Your team Simulateurs pour la Hadronthérapie is subscribed to branch lp:simulationshadrontherapie/g4lpc.
To unsubscribe from this branch go to https://code.launchpad.net/~simhadteam/simulationshadrontherapie/G4LPC/+edit-subscription

=== modified file 'src/G4LPCDetectorConstruction.cc'
--- src/G4LPCDetectorConstruction.cc	2010-04-08 11:30:31 +0000
+++ src/G4LPCDetectorConstruction.cc	2010-04-30 12:50:38 +0000
@@ -116,8 +116,42 @@
     new G4Material("Eau", density= 1.000*g/cm3, ncomponents=2);
   Eau->AddElement(H, natoms=2);
   Eau->AddElement(O, natoms=1);
+
+  G4Material* Eau_io75 = 
+    new G4Material("Eau_io75", density= 1.000*g/cm3, ncomponents=2);
+  Eau_io75->AddElement(H, natoms=2);
+  Eau_io75->AddElement(O, natoms=1);
 // overwrite computed meanExcitationEnergy with ICRU recommended value 
-  Eau->GetIonisation()->SetMeanExcitationEnergy(75.0*eV);
+  Eau_io75->GetIonisation()->SetMeanExcitationEnergy(75.0*eV);
+
+  G4Material* Eau_io70 = 
+    new G4Material("Eau_io70", density= 1.000*g/cm3, ncomponents=2);
+  Eau_io70->AddElement(H, natoms=2);
+  Eau_io70->AddElement(O, natoms=1);
+// overwrite computed meanExcitationEnergy
+  Eau_io70->GetIonisation()->SetMeanExcitationEnergy(70.0*eV);
+  //Eau_io70->GetIonisation()->SetMeanEnergyPerIonPair(70.0*eV);
+
+  G4Material* Eau_io65 = 
+    new G4Material("Eau_io65", density= 1.000*g/cm3, ncomponents=2);
+  Eau_io65->AddElement(H, natoms=2);
+  Eau_io65->AddElement(O, natoms=1);
+// overwrite computed meanExcitationEnergy
+  Eau_io65->GetIonisation()->SetMeanExcitationEnergy(65.0*eV);
+
+  G4Material* Eau_io60 = 
+    new G4Material("Eau_io60", density= 1.000*g/cm3, ncomponents=2);
+  Eau_io60->AddElement(H, natoms=2);
+  Eau_io60->AddElement(O, natoms=1);
+// overwrite computed meanExcitationEnergy
+  Eau_io60->GetIonisation()->SetMeanExcitationEnergy(60.0*eV);
+
+  G4Material* Eau_io80 = 
+    new G4Material("Eau_io80", density= 1.000*g/cm3, ncomponents=2);
+  Eau_io80->AddElement(H, natoms=2);
+  Eau_io80->AddElement(O, natoms=1);
+// overwrite computed meanExcitationEnergy
+  Eau_io80->GetIonisation()->SetMeanExcitationEnergy(80.0*eV);
 
 
  density = 1.290*mg/cm3;
@@ -1864,7 +1898,7 @@
 {
 if(IndexInRange(i))
  {
- G4Material *m= G4Material::GetMaterial(s); 
+ G4Material *m= G4Material::GetMaterial(s);
  if(!m)
   {
   G4NistManager* man = G4NistManager::Instance();

=== modified file 'src/G4LPCRunAction.cc'
--- src/G4LPCRunAction.cc	2009-10-21 09:27:24 +0000
+++ src/G4LPCRunAction.cc	2010-04-30 12:50:38 +0000
@@ -78,9 +78,14 @@
   cout << "Nombre de couches : " << det->GetNbLayers() << endl;
   for(Int_t i=0;i<det->GetNbLayers();i++)
    {
+   G4Material *mat=det->GetLayerMaterial(i);
    cout << "--> Couche " << i << " <-- "  << endl; 
-   cout << "  Materiau : " << det->GetLayerMaterial(i)->GetName() << endl;
-   cout << "  Epaisseur: " << det->GetLayerThickness(i)/mm << " mm." << endl;
+   cout << "  Materiau  : " << mat->GetName() << endl;
+   cout << "  Ionisation: " << mat->GetIonisation()
+                                ->GetMeanExcitationEnergy()/eV << " eV / " 
+                            << mat->GetIonisation()
+                                ->GetMeanEnergyPerIonPair()/eV << " eV." << endl;
+   cout << "  Epaisseur : " << det->GetLayerThickness(i)/mm << " mm." << endl;
    }
   }