Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm5/src/DetectorConstruction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/electromagnetic/TestEm5/src/DetectorConstruction.cc (Version 11.3.0) and /examples/extended/electromagnetic/TestEm5/src/DetectorConstruction.cc (Version 9.3.p2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 /// \file electromagnetic/TestEm5/src/Detector <<  26 // $Id: DetectorConstruction.cc,v 1.15 2009/01/22 17:41:43 vnivanch Exp $
 27 /// \brief Implementation of the DetectorConst <<  27 // GEANT4 tag $Name: geant4-09-03-patch-02 $
 28 //                                             << 
 29 //                                                 28 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo     29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    31 
 33 #include "DetectorConstruction.hh"                 32 #include "DetectorConstruction.hh"
 34                                                << 
 35 #include "DetectorMessenger.hh"                    33 #include "DetectorMessenger.hh"
 36                                                    34 
 37 #include "G4AutoDelete.hh"                     <<  35 #include "G4Material.hh"
 38 #include "G4Box.hh"                                36 #include "G4Box.hh"
 39 #include "G4GeometryManager.hh"                << 
 40 #include "G4GlobalMagFieldMessenger.hh"        << 
 41 #include "G4LogicalVolume.hh"                      37 #include "G4LogicalVolume.hh"
 42 #include "G4LogicalVolumeStore.hh"             << 
 43 #include "G4Material.hh"                       << 
 44 #include "G4NistManager.hh"                    << 
 45 #include "G4PVPlacement.hh"                        38 #include "G4PVPlacement.hh"
 46 #include "G4PhysicalConstants.hh"              <<  39 #include "G4UniformMagField.hh"
                                                   >>  40 
                                                   >>  41 #include "G4GeometryManager.hh"
 47 #include "G4PhysicalVolumeStore.hh"                42 #include "G4PhysicalVolumeStore.hh"
 48 #include "G4RunManager.hh"                     <<  43 #include "G4LogicalVolumeStore.hh"
 49 #include "G4SolidStore.hh"                         44 #include "G4SolidStore.hh"
 50 #include "G4SystemOfUnits.hh"                  <<  45 
 51 #include "G4UniformMagField.hh"                << 
 52 #include "G4UnitsTable.hh"                         46 #include "G4UnitsTable.hh"
                                                   >>  47 #include "G4NistManager.hh"
                                                   >>  48 #include "G4RunManager.hh"
 53                                                    49 
 54 //....oooOO0OOooo........oooOO0OOooo........oo     50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 55                                                    51 
 56 DetectorConstruction::DetectorConstruction()       52 DetectorConstruction::DetectorConstruction()
                                                   >>  53 :AbsorberMaterial(0),WorldMaterial(0),defaultWorld(true),
                                                   >>  54  solidWorld(0),logicWorld(0),physiWorld(0),
                                                   >>  55  solidAbsorber(0),logicAbsorber(0),physiAbsorber(0),
                                                   >>  56  magField(0)
 57 {                                                  57 {
 58   // default parameter values of the calorimet     58   // default parameter values of the calorimeter
 59   fAbsorberThickness = 1. * cm;                <<  59   AbsorberThickness = 1.*cm;
 60   fAbsorberSizeYZ = 2. * cm;                   <<  60   AbsorberSizeYZ    = 2.*cm;
 61   fXposAbs = 0. * cm;                          <<  61   XposAbs           = 0.*cm;
 62   ComputeGeomParameters();                     <<  62   ComputeCalorParameters();
 63                                                <<  63   
 64   // materials                                 <<  64   // materials  
 65   DefineMaterials();                               65   DefineMaterials();
 66   SetWorldMaterial("G4_Galactic");             <<  66   SetWorldMaterial   ("Galactic");
 67   SetAbsorberMaterial("G4_Si");                <<  67   SetAbsorberMaterial("Silicon");
 68                                                <<  68  
 69   // create commands for interactive definitio <<  69   // create commands for interactive definition of the calorimeter  
 70   fDetectorMessenger = new DetectorMessenger(t <<  70   detectorMessenger = new DetectorMessenger(this);
 71 }                                                  71 }
 72                                                    72 
 73 //....oooOO0OOooo........oooOO0OOooo........oo     73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 74                                                    74 
 75 DetectorConstruction::~DetectorConstruction()      75 DetectorConstruction::~DetectorConstruction()
 76 {                                              <<  76 { 
 77   delete fDetectorMessenger;                   <<  77   delete detectorMessenger;
 78 }                                                  78 }
 79                                                    79 
 80 //....oooOO0OOooo........oooOO0OOooo........oo     80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 81                                                    81 
 82 void DetectorConstruction::DefineMaterials()   <<  82 G4VPhysicalVolume* DetectorConstruction::Construct()
 83 {                                                  83 {
 84   // This function illustrates the possible wa <<  84   return ConstructCalorimeter();
                                                   >>  85 }
 85                                                    86 
 86   G4String symbol;  // a=mass of a mole;       <<  87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 87   G4double a, z, density;  // z=mean number of <<  88 
                                                   >>  89 void DetectorConstruction::DefineMaterials()
                                                   >>  90 { 
                                                   >>  91   //This function illustrates the possible ways to define materials
                                                   >>  92  
                                                   >>  93   G4String symbol;             //a=mass of a mole;
                                                   >>  94   G4double a, z, density;      //z=mean number of protons;  
 88                                                    95 
 89   G4int ncomponents, natoms;                       96   G4int ncomponents, natoms;
 90   G4double fractionmass;                           97   G4double fractionmass;
 91   G4double temperature, pressure;                  98   G4double temperature, pressure;
 92                                                <<  99   
 93   //                                              100   //
 94   // define Elements                              101   // define Elements
 95   //                                              102   //
 96                                                   103 
 97   G4Element* H = new G4Element("Hydrogen", sym << 104   G4Element* H  = new G4Element("Hydrogen",symbol="H",  z= 1, a=   1.01*g/mole);
 98   G4Element* C = new G4Element("Carbon", symbo << 105   G4Element* C  = new G4Element("Carbon",  symbol="C",  z= 6, a=  12.01*g/mole);
 99   G4Element* N = new G4Element("Nitrogen", sym << 106   G4Element* N  = new G4Element("Nitrogen",symbol="N",  z= 7, a=  14.01*g/mole);
100   G4Element* O = new G4Element("Oxygen", symbo << 107   G4Element* O  = new G4Element("Oxygen",  symbol="O",  z= 8, a=  16.00*g/mole);
101   G4Element* Na = new G4Element("Sodium", symb << 108   G4Element* Na = new G4Element("Sodium",  symbol="Na", z=11, a=  22.99*g/mole);
102   G4Element* Ar = new G4Element("Argon", symbo << 109   G4Element* Ar = new G4Element("Argon",   symbol="Ar", z=18, a=  39.95*g/mole);
103   G4Element* I = new G4Element("Iodine", symbo << 110   G4Element* I  = new G4Element("Iodine",  symbol="I" , z=53, a= 126.90*g/mole);
104   G4Element* Xe = new G4Element("Xenon", symbo << 111   G4Element* Xe = new G4Element("Xenon",   symbol="Xe", z=54, a= 131.29*g/mole);
105                                                   112 
106   //                                              113   //
107   // define simple materials                      114   // define simple materials
108   //                                              115   //
109                                                   116 
110   new G4Material("H2Liq", z = 1, a = 1.01 * g  << 117   new G4Material("H2Liq"    , z= 1, a= 1.01*g/mole, density= 70.8*mg/cm3);
111   new G4Material("Beryllium", z = 4, a = 9.01  << 118   new G4Material("Beryllium", z= 4, a= 9.01*g/mole, density= 1.848*g/cm3);
112   new G4Material("Aluminium", z = 13, a = 26.9 << 119   new G4Material("Aluminium", z=13, a=26.98*g/mole, density= 2.700*g/cm3);
113   new G4Material("Silicon", z = 14, a = 28.09  << 120   new G4Material("Silicon"  , z=14, a=28.09*g/mole, density= 2.330*g/cm3);
114                                                << 121 
115   G4Material* lAr = new G4Material("liquidArgo << 122   G4Material* lAr = 
116   lAr->AddElement(Ar, natoms = 1);             << 123     new G4Material("liquidArgon", density= 1.390*g/cm3, ncomponents=1);
117                                                << 124   lAr->AddElement(Ar, natoms=1);
118   new G4Material("Iron", z = 26, a = 55.85 * g << 125 
119   new G4Material("Copper", z = 29, a = 63.55 * << 126   new G4Material("Iron",     z=26, a= 55.85*g/mole, density= 7.870*g/cm3);
120   new G4Material("Germanium", z = 32, a = 72.6 << 127   new G4Material("Copper",   z=29, a= 63.55*g/mole, density= 8.960*g/cm3);
121   new G4Material("Silver", z = 47, a = 107.87  << 128   new G4Material("Germanium",z=32, a= 72.61*g/mole, density= 5.323*g/cm3);
122   new G4Material("Tungsten", z = 74, a = 183.8 << 129   new G4Material("Silver",   z=47, a=107.87*g/mole, density= 10.50*g/cm3);
123   new G4Material("Gold", z = 79, a = 196.97 *  << 130   new G4Material("Tungsten", z=74, a=183.85*g/mole, density= 19.30*g/cm3);
124   new G4Material("Lead", z = 82, a = 207.19 *  << 131   new G4Material("Gold",     z=79, a=196.97*g/mole, density= 19.32*g/cm3);
                                                   >> 132   new G4Material("Lead",     z=82, a=207.19*g/mole, density= 11.35*g/cm3);
125                                                   133 
126   //                                              134   //
127   // define a material from elements.   case 1    135   // define a material from elements.   case 1: chemical molecule
128   //                                              136   //
129                                                   137 
130   G4Material* H2O = new G4Material("Water", de << 138   G4Material* H2O = new G4Material("Water", density= 1.000*g/cm3, ncomponents=2);
131   H2O->AddElement(H, natoms = 2);              << 139   H2O->AddElement(H, natoms=2);
132   H2O->AddElement(O, natoms = 1);              << 140   H2O->AddElement(O, natoms=1);
133   H2O->GetIonisation()->SetMeanExcitationEnerg << 141   H2O->GetIonisation()->SetMeanExcitationEnergy(75*eV);
134                                                << 142 
135   G4Material* CH = new G4Material("Plastic", d << 143   G4Material* CH = new G4Material("Plastic", density= 1.04*g/cm3, ncomponents=2);
136   CH->AddElement(C, natoms = 1);               << 144   CH->AddElement(C, natoms=1);
137   CH->AddElement(H, natoms = 1);               << 145   CH->AddElement(H, natoms=1);
138                                                << 146 
139   G4Material* NaI = new G4Material("NaI", dens << 147   G4Material* NaI = new G4Material("NaI", density= 3.67*g/cm3, ncomponents=2);
140   NaI->AddElement(Na, natoms = 1);             << 148   NaI->AddElement(Na, natoms=1);
141   NaI->AddElement(I, natoms = 1);              << 149   NaI->AddElement(I , natoms=1);
142   NaI->GetIonisation()->SetMeanExcitationEnerg << 150   NaI->GetIonisation()->SetMeanExcitationEnergy(452*eV);
143                                                   151 
144   //                                              152   //
145   // define a material from elements.   case 2    153   // define a material from elements.   case 2: mixture by fractional mass
146   //                                              154   //
147                                                   155 
148   G4Material* Air = new G4Material("Air", dens << 156   G4Material* Air = new G4Material("Air", density= 1.290*mg/cm3, ncomponents=2);
149   Air->AddElement(N, fractionmass = 0.7);      << 157   Air->AddElement(N, fractionmass=0.7);
150   Air->AddElement(O, fractionmass = 0.3);      << 158   Air->AddElement(O, fractionmass=0.3);
151                                                << 159 
152   G4Material* Air20 = new G4Material("Air20",  << 160   G4Material* Air20 = 
153                                      kStateGas << 161     new G4Material("Air20", density= 1.205*mg/cm3, ncomponents=2,
154   Air20->AddElement(N, fractionmass = 0.7);    << 162        kStateGas, 293.*kelvin, 1.*atmosphere);
155   Air20->AddElement(O, fractionmass = 0.3);    << 163   Air20->AddElement(N, fractionmass=0.7);
156                                                << 164   Air20->AddElement(O, fractionmass=0.3);
157   // Graphite                                  << 165 
158   //                                           << 166   //Graphite
159   G4Material* Graphite = new G4Material("Graph << 167   //
160   Graphite->AddElement(C, fractionmass = 1.);  << 168   G4Material* Graphite = 
161                                                << 169     new G4Material("Graphite", density= 1.7*g/cm3, ncomponents=1);
162   // Havar                                     << 170   Graphite->AddElement(C, fractionmass=1.);
163   //                                           << 171 
164   G4Element* Cr = new G4Element("Chrome", "Cr" << 172   //Havar
165   G4Element* Fe = new G4Element("Iron", "Fe",  << 173   //
166   G4Element* Co = new G4Element("Cobalt", "Co" << 174   G4Element* Cr = new G4Element("Chrome", "Cr", z=25, a=  51.996*g/mole);
167   G4Element* Ni = new G4Element("Nickel", "Ni" << 175   G4Element* Fe = new G4Element("Iron"  , "Fe", z=26, a=  55.845*g/mole);
168   G4Element* W = new G4Element("Tungsten", "W" << 176   G4Element* Co = new G4Element("Cobalt", "Co", z=27, a=  58.933*g/mole);
169                                                << 177   G4Element* Ni = new G4Element("Nickel", "Ni", z=28, a=  58.693*g/mole);
170   G4Material* Havar = new G4Material("Havar",  << 178   G4Element* W  = new G4Element("Tungsten","W", z=74, a= 183.850*g/mole);
171   Havar->AddElement(Cr, fractionmass = 0.1785) << 179 
172   Havar->AddElement(Fe, fractionmass = 0.1822) << 180   G4Material* Havar = 
173   Havar->AddElement(Co, fractionmass = 0.4452) << 181     new G4Material("Havar", density= 8.3*g/cm3, ncomponents=5);
174   Havar->AddElement(Ni, fractionmass = 0.1310) << 182   Havar->AddElement(Cr, fractionmass=0.1785);
175   Havar->AddElement(W, fractionmass = 0.0631); << 183   Havar->AddElement(Fe, fractionmass=0.1822);
                                                   >> 184   Havar->AddElement(Co, fractionmass=0.4452);
                                                   >> 185   Havar->AddElement(Ni, fractionmass=0.1310);
                                                   >> 186   Havar->AddElement(W , fractionmass=0.0631);
176                                                   187 
177   //                                              188   //
178   // examples of gas                              189   // examples of gas
179   //                                           << 190   //  
180   new G4Material("ArgonGas", z = 18, a = 39.94 << 191   new G4Material("ArgonGas", z=18, a=39.948*g/mole, density= 1.782*mg/cm3,
181                  273.15 * kelvin, 1 * atmosphe << 192      kStateGas, 273.15*kelvin, 1*atmosphere);
182                                                << 193          
183   new G4Material("XenonGas", z = 54, a = 131.2 << 194   new G4Material("XenonGas", z=54, a=131.29*g/mole, density= 5.458*mg/cm3,
184                  293.15 * kelvin, 1 * atmosphe << 195      kStateGas, 293.15*kelvin, 1*atmosphere);
185                                                << 196          
186   G4Material* CO2 = new G4Material("CarbonicGa << 197   G4Material* CO2 =
187   CO2->AddElement(C, natoms = 1);              << 198     new G4Material("CarbonicGas", density= 1.977*mg/cm3, ncomponents=2);
188   CO2->AddElement(O, natoms = 2);              << 199   CO2->AddElement(C, natoms=1);
                                                   >> 200   CO2->AddElement(O, natoms=2);
                                                   >> 201 
                                                   >> 202   G4Material* ArCO2 =
                                                   >> 203     new G4Material("ArgonCO2",   density= 1.8223*mg/cm3, ncomponents=2);
                                                   >> 204   ArCO2->AddElement (Ar,  fractionmass=0.7844);
                                                   >> 205   ArCO2->AddMaterial(CO2, fractionmass=0.2156);
189                                                   206 
190   G4Material* ArCO2 = new G4Material("ArgonCO2 << 207   //another way to define mixture of gas per volume
191   ArCO2->AddElement(Ar, fractionmass = 0.7844) << 
192   ArCO2->AddMaterial(CO2, fractionmass = 0.215 << 
193                                                << 
194   // another way to define mixture of gas per  << 
195   G4Material* NewArCO2 =                          208   G4Material* NewArCO2 =
196     new G4Material("NewArgonCO2", density = 1. << 209     new G4Material("NewArgonCO2", density= 1.8223*mg/cm3, ncomponents=3);
197   NewArCO2->AddElement(Ar, natoms = 8);        << 210   NewArCO2->AddElement (Ar, natoms=8);
198   NewArCO2->AddElement(C, natoms = 2);         << 211   NewArCO2->AddElement (C,  natoms=2);
199   NewArCO2->AddElement(O, natoms = 4);         << 212   NewArCO2->AddElement (O,  natoms=4);
200                                                << 213 
201   G4Material* ArCH4 = new G4Material("ArgonCH4 << 214   G4Material* ArCH4 = 
202   ArCH4->AddElement(Ar, natoms = 93);          << 215     new G4Material("ArgonCH4",    density= 1.709*mg/cm3,  ncomponents=3);
203   ArCH4->AddElement(C, natoms = 7);            << 216   ArCH4->AddElement (Ar, natoms=93);
204   ArCH4->AddElement(H, natoms = 28);           << 217   ArCH4->AddElement (C,  natoms=7);
205                                                << 218   ArCH4->AddElement (H,  natoms=28);
206   G4Material* XeCH = new G4Material("XenonMeth << 219 
207                                     ncomponent << 220   G4Material* XeCH = 
208   XeCH->AddElement(Xe, natoms = 875);          << 221     new G4Material("XenonMethanePropane", density= 4.9196*mg/cm3, ncomponents=3,
209   XeCH->AddElement(C, natoms = 225);           << 222        kStateGas, 293.15*kelvin, 1*atmosphere);
210   XeCH->AddElement(H, natoms = 700);           << 223   XeCH->AddElement (Xe, natoms=875);
211                                                << 224   XeCH->AddElement (C,  natoms=225);
212   G4Material* steam = new G4Material("WaterSte << 225   XeCH->AddElement (H,  natoms=700);
213   steam->AddMaterial(H2O, fractionmass = 1.);  << 226 
214   steam->GetIonisation()->SetMeanExcitationEne << 227   G4Material* steam = 
215                                                << 228     new G4Material("WaterSteam", density= 1.0*mg/cm3, ncomponents=1);
216   G4Material* rock1 = new G4Material("Standard << 229   steam->AddMaterial(H2O, fractionmass=1.);
217   rock1->AddElement(Na, 1);                    << 230   steam->GetIonisation()->SetMeanExcitationEnergy(71.6*eV);  
218                                                   231 
219   //                                              232   //
220   // example of vacuum                            233   // example of vacuum
221   //                                              234   //
222   density = universe_mean_density;  // from Ph << 235 
223   pressure = 3.e-18 * pascal;                  << 236   density     = universe_mean_density;    //from PhysicalConstants.h
224   temperature = 2.73 * kelvin;                 << 237   pressure    = 3.e-18*pascal;
225   new G4Material("Galactic", z = 1, a = 1.01 * << 238   temperature = 2.73*kelvin;
                                                   >> 239   new G4Material("Galactic", z=1, a=1.01*g/mole,density,
                                                   >> 240      kStateGas,temperature,pressure);
226 }                                                 241 }
227                                                   242 
228 //....oooOO0OOooo........oooOO0OOooo........oo    243 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
229                                                   244 
230 void DetectorConstruction::ComputeGeomParamete << 245 void DetectorConstruction::ComputeCalorParameters()
231 {                                                 246 {
232   // Compute derived parameters of the calorim    247   // Compute derived parameters of the calorimeter
233   fXstartAbs = fXposAbs - 0.5 * fAbsorberThick << 248   xstartAbs = XposAbs-0.5*AbsorberThickness; 
234   fXendAbs = fXposAbs + 0.5 * fAbsorberThickne << 249   xendAbs   = XposAbs+0.5*AbsorberThickness;
235                                                << 250      
236   G4double xmax = std::max(std::abs(fXstartAbs << 251   if (defaultWorld) {
237   fWorldSizeX = 2.4 * xmax;                    << 252      WorldSizeX = 1.5*AbsorberThickness; WorldSizeYZ= 1.2*AbsorberSizeYZ;
238   fWorldSizeYZ = 1.2 * fAbsorberSizeYZ;        << 253   }   
239   if (nullptr != fPhysiWorld) {                << 
240     ChangeGeometry();                          << 
241   }                                            << 
242 }                                                 254 }
243                                                   255 
244 //....oooOO0OOooo........oooOO0OOooo........oo    256 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
245                                                << 257   
246 G4VPhysicalVolume* DetectorConstruction::Const << 258 G4VPhysicalVolume* DetectorConstruction::ConstructCalorimeter()
247 {                                              << 259 { 
248   if (nullptr != fPhysiWorld) {                << 260   // Cleanup old geometry
249     return fPhysiWorld;                        << 261   //
250   }                                            << 262   G4GeometryManager::GetInstance()->OpenGeometry();
                                                   >> 263   G4PhysicalVolumeStore::GetInstance()->Clean();
                                                   >> 264   G4LogicalVolumeStore::GetInstance()->Clean();
                                                   >> 265   G4SolidStore::GetInstance()->Clean();
                                                   >> 266   
                                                   >> 267   // complete the Calor parameters definition 
                                                   >> 268   ComputeCalorParameters();
                                                   >> 269         
251   // World                                        270   // World
252   //                                              271   //
253   fSolidWorld = new G4Box("World",  // its nam << 272   solidWorld = new G4Box("World",       //its name
254                           fWorldSizeX / 2, fWo << 273                    WorldSizeX/2,WorldSizeYZ/2,WorldSizeYZ/2);   //its size
255                                                << 274                          
256   fLogicWorld = new G4LogicalVolume(fSolidWorl << 275   logicWorld = new G4LogicalVolume(solidWorld,    //its solid
257                                     fWorldMate << 276                                    WorldMaterial, //its material
258                                     "World");  << 277                                    "World");    //its name
259                                                << 278                                    
260   fPhysiWorld = new G4PVPlacement(0,  // no ro << 279   physiWorld = new G4PVPlacement(0,     //no rotation
261                                   G4ThreeVecto << 280            G4ThreeVector(), //at (0,0,0)
262                                   fLogicWorld, << 281                                  logicWorld,    //its logical volume
263                                   "World",  // << 282                                  "World",   //its name
264                                   0,  // its m << 283                                  0,     //its mother  volume
265                                   false,  // n << 284                                  false,     //no boolean operation
266                                   0);  // copy << 285                                  0);      //copy number
267                                                << 286                                  
268   // Absorber                                     287   // Absorber
269   //                                           << 288   // 
270   fSolidAbsorber =                             << 289   solidAbsorber = new G4Box("Absorber", 
271     new G4Box("Absorber", fAbsorberThickness / << 290                       AbsorberThickness/2,AbsorberSizeYZ/2,AbsorberSizeYZ/2); 
272                                                << 291                           
273   fLogicAbsorber = new G4LogicalVolume(fSolidA << 292   logicAbsorber = new G4LogicalVolume(solidAbsorber,    //its solid
274                                        fAbsorb << 293                                 AbsorberMaterial, //its material
275                                        "Absorb << 294                              "Absorber");       //its name
276                                                << 295                               
277   fPhysiAbsorber = new G4PVPlacement(0,  // no << 296   physiAbsorber = new G4PVPlacement(0,       //no rotation
278                                      G4ThreeVe << 297             G4ThreeVector(XposAbs,0.,0.),    //its position
279                                      fLogicAbs << 298                                 logicAbsorber,     //its logical volume
280                                      "Absorber << 299         "Absorber",         //its name
281                                      fLogicWor << 300                                 logicWorld,        //its mother
282                                      false,  / << 301                                 false,             //no boulean operat
283                                      0);  // c << 302                                 0);                //copy number
284                                                << 303                                         
285   PrintGeomParameters();                       << 304   PrintCalorParameters();         
286                                                << 305   
287   // always return the physical World          << 306   //always return the physical World
288   //                                           << 307   //
289   return fPhysiWorld;                          << 308   return physiWorld;
290 }                                                 309 }
291                                                   310 
292 //....oooOO0OOooo........oooOO0OOooo........oo    311 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
293                                                   312 
294 void DetectorConstruction::PrintGeomParameters << 313 void DetectorConstruction::PrintCalorParameters()
295 {                                              << 314 {
296   G4cout << "\n" << fWorldMaterial << G4endl;  << 315   G4cout << "\n" << WorldMaterial    << G4endl;
297   G4cout << "\n" << fAbsorberMaterial << G4end << 316   G4cout << "\n" << AbsorberMaterial << G4endl;
298                                                << 317     
299   G4cout << "\n The  WORLD   is made of " << G << 318   G4cout << "\n The  WORLD   is made of "  << G4BestUnit(WorldSizeX,"Length")
300          << fWorldMaterial->GetName();         << 319          << " of " << WorldMaterial->GetName();
301   G4cout << ". The transverse size (YZ) of the << 320   G4cout << ". The transverse size (YZ) of the world is " 
302          << G4endl;                            << 321          << G4BestUnit(WorldSizeYZ,"Length") << G4endl;
303   G4cout << " The ABSORBER is made of " << G4B << 322   G4cout << " The ABSORBER is made of " 
304          << fAbsorberMaterial->GetName();      << 323          <<G4BestUnit(AbsorberThickness,"Length")
305   G4cout << ". The transverse size (YZ) is " < << 324    << " of " << AbsorberMaterial->GetName();
306   G4cout << " X position of the middle of the  << 325   G4cout << ". The transverse size (YZ) is " 
                                                   >> 326          << G4BestUnit(AbsorberSizeYZ,"Length") << G4endl;
                                                   >> 327   G4cout << " X position of the middle of the absorber "
                                                   >> 328          << G4BestUnit(XposAbs,"Length");
307   G4cout << G4endl;                               329   G4cout << G4endl;
308 }                                                 330 }
309                                                   331 
310 //....oooOO0OOooo........oooOO0OOooo........oo    332 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
311                                                   333 
312 void DetectorConstruction::SetAbsorberMaterial << 334 void DetectorConstruction::SetAbsorberMaterial(G4String materialChoice)
313 {                                                 335 {
314   // search the material by its name              336   // search the material by its name
315   G4Material* pttoMaterial = G4NistManager::In << 337   G4Material* pttoMaterial =
                                                   >> 338     G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
316                                                   339 
317   if (pttoMaterial && fAbsorberMaterial != ptt << 340   if (pttoMaterial && AbsorberMaterial != pttoMaterial) {
318     fAbsorberMaterial = pttoMaterial;          << 341     AbsorberMaterial = pttoMaterial;                  
319     if (fLogicAbsorber) {                      << 342     if(logicAbsorber) logicAbsorber->SetMaterial(AbsorberMaterial);
320       fLogicAbsorber->SetMaterial(fAbsorberMat << 
321     }                                          << 
322     G4RunManager::GetRunManager()->PhysicsHasB    343     G4RunManager::GetRunManager()->PhysicsHasBeenModified();
323   }                                               344   }
324 }                                                 345 }
325                                                   346 
326 //....oooOO0OOooo........oooOO0OOooo........oo    347 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
327                                                   348 
328 void DetectorConstruction::SetWorldMaterial(co << 349 void DetectorConstruction::SetWorldMaterial(G4String materialChoice)
329 {                                                 350 {
330   // search the material by its name              351   // search the material by its name
331   G4Material* pttoMaterial = G4NistManager::In << 352   G4Material* pttoMaterial =
                                                   >> 353     G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
332                                                   354 
333   if (pttoMaterial && fWorldMaterial != pttoMa << 355   if (pttoMaterial && WorldMaterial != pttoMaterial) {
334     fWorldMaterial = pttoMaterial;             << 356     WorldMaterial = pttoMaterial;                  
335     if (fLogicWorld) {                         << 357     if(logicWorld) logicWorld->SetMaterial(WorldMaterial);
336       fLogicWorld->SetMaterial(fWorldMaterial) << 
337     }                                          << 
338     G4RunManager::GetRunManager()->PhysicsHasB    358     G4RunManager::GetRunManager()->PhysicsHasBeenModified();
339   }                                               359   }
340 }                                                 360 }
341                                                << 361     
342 //....oooOO0OOooo........oooOO0OOooo........oo    362 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
343                                                   363 
344 void DetectorConstruction::SetAbsorberThicknes    364 void DetectorConstruction::SetAbsorberThickness(G4double val)
345 {                                                 365 {
346   fAbsorberThickness = val;                    << 366   AbsorberThickness = val;
347   ComputeGeomParameters();                     << 367   G4RunManager::GetRunManager()->GeometryHasBeenModified();
348 }                                              << 368 }  
349                                                   369 
350 //....oooOO0OOooo........oooOO0OOooo........oo    370 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
351                                                   371 
352 void DetectorConstruction::SetAbsorberSizeYZ(G    372 void DetectorConstruction::SetAbsorberSizeYZ(G4double val)
353 {                                                 373 {
354   fAbsorberSizeYZ = val;                       << 374   AbsorberSizeYZ = val;
355   ComputeGeomParameters();                     << 375   G4RunManager::GetRunManager()->GeometryHasBeenModified();
356 }                                              << 376 }  
357                                                   377 
358 //....oooOO0OOooo........oooOO0OOooo........oo    378 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
359                                                   379 
360 void DetectorConstruction::SetWorldSizeX(G4dou    380 void DetectorConstruction::SetWorldSizeX(G4double val)
361 {                                                 381 {
362   fWorldSizeX = val;                           << 382   WorldSizeX = val;
363   ComputeGeomParameters();                     << 383   defaultWorld = false;
364 }                                              << 384   G4RunManager::GetRunManager()->GeometryHasBeenModified();
                                                   >> 385 }  
365                                                   386 
366 //....oooOO0OOooo........oooOO0OOooo........oo    387 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
367                                                   388 
368 void DetectorConstruction::SetWorldSizeYZ(G4do    389 void DetectorConstruction::SetWorldSizeYZ(G4double val)
369 {                                                 390 {
370   fWorldSizeYZ = val;                          << 391   WorldSizeYZ = val;
371   ComputeGeomParameters();                     << 392   defaultWorld = false;
372 }                                              << 393   G4RunManager::GetRunManager()->GeometryHasBeenModified();
                                                   >> 394 }  
373                                                   395 
374 //....oooOO0OOooo........oooOO0OOooo........oo    396 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
375                                                   397 
376 void DetectorConstruction::SetAbsorberXpos(G4d    398 void DetectorConstruction::SetAbsorberXpos(G4double val)
377 {                                                 399 {
378   fXposAbs = val;                              << 400   XposAbs  = val;
379   ComputeGeomParameters();                     << 401   G4RunManager::GetRunManager()->GeometryHasBeenModified();
380 }                                              << 402 }  
381                                                   403 
382 //....oooOO0OOooo........oooOO0OOooo........oo    404 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
383                                                   405 
384 void DetectorConstruction::ConstructSDandField << 406 #include "G4FieldManager.hh"
                                                   >> 407 #include "G4TransportationManager.hh"
                                                   >> 408 
                                                   >> 409 void DetectorConstruction::SetMagField(G4double fieldValue)
385 {                                                 410 {
386   if (fFieldMessenger.Get() == 0) {            << 411   //apply a global uniform magnetic field along Z axis
387     // Create global magnetic field messenger. << 412   G4FieldManager* fieldMgr 
388     // Uniform magnetic field is then created  << 413    = G4TransportationManager::GetTransportationManager()->GetFieldManager();
389     // the field value is not zero.            << 414     
390     G4ThreeVector fieldValue = G4ThreeVector() << 415   if(magField) delete magField;   //delete the existing magn field
391     G4GlobalMagFieldMessenger* msg = new G4Glo << 416   
392     // msg->SetVerboseLevel(1);                << 417   if(fieldValue!=0.)      // create a new one if non nul
393     G4AutoDelete::Register(msg);               << 418   { magField = new G4UniformMagField(G4ThreeVector(0.,0.,fieldValue));        
394     fFieldMessenger.Put(msg);                  << 419     fieldMgr->SetDetectorField(magField);
                                                   >> 420     fieldMgr->CreateChordFinder(magField);
                                                   >> 421   } else {
                                                   >> 422     magField = NULL;
                                                   >> 423     fieldMgr->SetDetectorField(magField);
395   }                                               424   }
396 }                                                 425 }
397                                                   426 
398 //....oooOO0OOooo........oooOO0OOooo........oo << 427 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
399                                                   428 
400 void DetectorConstruction::ChangeGeometry()    << 429 void DetectorConstruction::UpdateGeometry()
401 {                                                 430 {
402   fSolidWorld->SetXHalfLength(fWorldSizeX * 0. << 431   G4RunManager::GetRunManager()->PhysicsHasBeenModified();
403   fSolidWorld->SetYHalfLength(fWorldSizeYZ * 0 << 432   G4RunManager::GetRunManager()->DefineWorldVolume(ConstructCalorimeter());
404   fSolidWorld->SetZHalfLength(fWorldSizeYZ * 0 << 
405                                                << 
406   fSolidAbsorber->SetXHalfLength(fAbsorberThic << 
407   fSolidAbsorber->SetYHalfLength(fAbsorberSize << 
408   fSolidAbsorber->SetZHalfLength(fAbsorberSize << 
409 }                                                 433 }
410                                                   434 
411 //....oooOO0OOooo........oooOO0OOooo........oo    435 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 436 
412                                                   437