Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/fanoCavity2/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/medical/fanoCavity2/src/DetectorConstruction.cc (Version 11.3.0) and /examples/extended/medical/fanoCavity2/src/DetectorConstruction.cc (Version 9.6.p3)


  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 medical/fanoCavity2/src/DetectorCons     26 /// \file medical/fanoCavity2/src/DetectorConstruction.cc
 27 /// \brief Implementation of the DetectorConst     27 /// \brief Implementation of the DetectorConstruction class
 28 //                                                 28 //
                                                   >>  29 // $Id$
                                                   >>  30 
 29 //                                                 31 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo     32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    34 
 33 #include "DetectorConstruction.hh"                 35 #include "DetectorConstruction.hh"
 34                                                << 
 35 #include "DetectorMessenger.hh"                    36 #include "DetectorMessenger.hh"
 36                                                    37 
 37 #include "G4GeometryManager.hh"                << 
 38 #include "G4LogicalVolume.hh"                  << 
 39 #include "G4LogicalVolumeStore.hh"             << 
 40 #include "G4Material.hh"                           38 #include "G4Material.hh"
 41 #include "G4NistManager.hh"                    <<  39 #include "G4Tubs.hh"
                                                   >>  40 #include "G4LogicalVolume.hh"
                                                   >>  41 #include "G4VPhysicalVolume.hh"
 42 #include "G4PVPlacement.hh"                        42 #include "G4PVPlacement.hh"
 43 #include "G4PhysicalConstants.hh"              <<  43 
                                                   >>  44 #include "G4GeometryManager.hh"
 44 #include "G4PhysicalVolumeStore.hh"                45 #include "G4PhysicalVolumeStore.hh"
                                                   >>  46 #include "G4LogicalVolumeStore.hh"
 45 #include "G4SolidStore.hh"                         47 #include "G4SolidStore.hh"
 46 #include "G4SystemOfUnits.hh"                  <<  48 
 47 #include "G4Tubs.hh"                           << 
 48 #include "G4UnitsTable.hh"                         49 #include "G4UnitsTable.hh"
 49 #include "G4VPhysicalVolume.hh"                <<  50 #include "G4PhysicalConstants.hh"
                                                   >>  51 #include "G4SystemOfUnits.hh"
 50                                                    52 
 51 //....oooOO0OOooo........oooOO0OOooo........oo     53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 52                                                    54 
 53 DetectorConstruction::DetectorConstruction() : <<  55 DetectorConstruction::DetectorConstruction()
                                                   >>  56 :fWall(0), fCavity(0)
 54 {                                                  57 {
 55   // default parameter values                      58   // default parameter values
 56   fWallThickness = 5 * cm;                     <<  59   fWallThickness   = 5*cm;  
 57   fCavityThickness = 2 * mm;                   <<  60   fCavityThickness = 2*mm;
 58   fCavityRadius = 10 * m;                      <<  61   fWorldRadius     = 10*m;    
 59                                                <<  62   
 60   DefineMaterials();                               63   DefineMaterials();
 61   SetWallMaterial("G4_WATER");                 <<  64   SetWallMaterial("Water");
 62   SetCavityMaterial("g4Water_gas");            <<  65   
 63                                                <<  66   // create commands for interactive definition of the detector  
 64   // create commands for interactive definitio << 
 65   fDetectorMessenger = new DetectorMessenger(t     67   fDetectorMessenger = new DetectorMessenger(this);
 66 }                                                  68 }
 67                                                    69 
 68 //....oooOO0OOooo........oooOO0OOooo........oo     70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 69                                                    71 
 70 DetectorConstruction::~DetectorConstruction()      72 DetectorConstruction::~DetectorConstruction()
                                                   >>  73 { delete fDetectorMessenger;}
                                                   >>  74 
                                                   >>  75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  76 
                                                   >>  77 G4VPhysicalVolume* DetectorConstruction::Construct()
 71 {                                                  78 {
 72   delete fDetectorMessenger;                   <<  79   return ConstructVolumes();
 73 }                                                  80 }
 74                                                    81 
 75 //....oooOO0OOooo........oooOO0OOooo........oo     82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 76                                                    83 
 77 void DetectorConstruction::DefineMaterials()       84 void DetectorConstruction::DefineMaterials()
 78 {                                              <<  85 { 
 79   G4double z, a;                               <<  86   //
 80                                                <<  87   // define Elements
 81   G4Element* H = new G4Element("Hydrogen", "H" <<  88   //
 82   G4Element* O = new G4Element("Oxygen", "O",  <<  89   G4double z,a;
                                                   >>  90   
                                                   >>  91   G4Element* H  = new G4Element("Hydrogen" ,"H" , z= 1., a=   1.01*g/mole);
                                                   >>  92   G4Element* O  = new G4Element("Oxygen"   ,"O" , z= 8., a=  16.00*g/mole);
 83                                                    93 
 84   G4Material* H2O = new G4Material("Water", 1. <<  94   //
                                                   >>  95   // define materials
                                                   >>  96   //  
                                                   >>  97   G4Material* H2O = 
                                                   >>  98   new G4Material("Water", 1.0*g/cm3, 2);
 85   H2O->AddElement(H, 2);                           99   H2O->AddElement(H, 2);
 86   H2O->AddElement(O, 1);                          100   H2O->AddElement(O, 1);
 87   H2O->GetIonisation()->SetMeanExcitationEnerg << 101   H2O->GetIonisation()->SetMeanExcitationEnergy(78.0*eV);
 88                                                << 102   
 89   G4Material* gas = new G4Material("Water_gas" << 103   G4Material* gas = 
                                                   >> 104   new G4Material("Water_gas", 1.0*mg/cm3, 2);
 90   gas->AddElement(H, 2);                          105   gas->AddElement(H, 2);
 91   gas->AddElement(O, 1);                          106   gas->AddElement(O, 1);
 92   gas->GetIonisation()->SetMeanExcitationEnerg << 107   gas->GetIonisation()->SetMeanExcitationEnergy(78.0*eV);
 93                                                << 108   
 94   new G4Material("Graphite", 6, 12.01 * g / mo << 109   new G4Material("Graphite",     6, 12.01*g/mole, 2.265*g/cm3);
 95   new G4Material("Graphite_gas", 6, 12.01 * g  << 110   new G4Material("Graphite_gas", 6, 12.01*g/mole, 2.265*mg/cm3);  
 96                                                << 111   
 97   new G4Material("Aluminium", 13, 26.98 * g /  << 112   new G4Material("Aluminium",     13, 26.98*g/mole, 2.700*g/cm3);
 98   new G4Material("Aluminium_gas", 13, 26.98 *  << 113   new G4Material("Aluminium_gas", 13, 26.98*g/mole, 2.700*mg/cm3);  
 99                                                << 114           
100   // alternatively, use G4 data base           << 115  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
101   //                                           << 
102   G4NistManager* nist = G4NistManager::Instanc << 
103                                                << 
104   nist->FindOrBuildMaterial("G4_WATER");       << 
105   nist->BuildMaterialWithNewDensity("g4Water_g << 
106                                                << 
107   // printout                                  << 
108   G4cout << *(G4Material::GetMaterialTable())  << 
109 }                                                 116 }
110                                                   117 
111 //....oooOO0OOooo........oooOO0OOooo........oo    118 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
112                                                << 119   
113 G4VPhysicalVolume* DetectorConstruction::Const << 120 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
114 {                                              << 121 {
115   if (fWall) {                                 << 122   G4GeometryManager::GetInstance()->OpenGeometry();
116     return fWall;                              << 123   G4PhysicalVolumeStore::GetInstance()->Clean();
117   }                                            << 124   G4LogicalVolumeStore::GetInstance()->Clean();
118                                                << 125   G4SolidStore::GetInstance()->Clean();
119   // Chamber                                   << 126                    
                                                   >> 127   // Wall
120   //                                              128   //
121   fTotalThickness = fCavityThickness + 2 * fWa << 129   fWorldThickness = fCavityThickness + 2*fWallThickness;
122   fWallRadius = fCavityRadius;                 << 130   
123                                                << 131   G4Tubs* 
124   G4Tubs* sChamber = new G4Tubs("Chamber",  // << 132   sWall = new G4Tubs("Wall",                                         //name
125                                 0., fWallRadiu << 133                       0.,fWorldRadius,0.5*fWorldThickness,0.,twopi); //size
126                                                << 134 
127   G4LogicalVolume* lChamber = new G4LogicalVol << 135   G4LogicalVolume*
128                                                << 136   lWall = new G4LogicalVolume(sWall,                   //solid
129                                                << 137                               fWallMaterial,           //material
130                                                << 138                               "Wall");                 //name
131   fWall = new G4PVPlacement(0,  // no rotation << 139                                    
132                             G4ThreeVector(),   << 140   fWall = new G4PVPlacement(0,                         //no rotation
133                             lChamber,  // logi << 141                             G4ThreeVector(),           //at (0,0,0)
134                             "Wall",  // name   << 142                             lWall,                     //logical volume
135                             0,  // mother  vol << 143                             "Wall",                    //name
136                             false,  // no bool << 144                             0,                         //mother  volume
137                             0);  // copy numbe << 145                             false,                     //no boolean operation
                                                   >> 146                             0);                        //copy number
138                                                   147 
139   // Cavity                                       148   // Cavity
140   //                                           << 149   //                             
141   G4Tubs* sCavity = new G4Tubs("Cavity", 0., f << 150   G4Tubs*
142                                                << 151   sCavity = new G4Tubs("Cavity",        
143   G4LogicalVolume* lCavity = new G4LogicalVolu << 152                        0.,fWorldRadius,0.5*fCavityThickness,0.,twopi);
144                                                << 153                  
145                                                << 154   G4LogicalVolume*                                                                 
146                                                << 155   lCavity = new G4LogicalVolume(sCavity,               //shape
147   fCavity = new G4PVPlacement(0,  // no rotati << 156                                 fCavityMaterial,       //material
148                               G4ThreeVector(), << 157                                 "Cavity");             //name
149                               lCavity,  // log << 158                                 
150                               "Cavity",  // na << 159   fCavity = new G4PVPlacement(0,                       //no rotation
151                               lChamber,  // mo << 160                              G4ThreeVector(),          //at (0,0,0)
152                               false,  // no bo << 161                              lCavity,                  //logical volume
153                               1);  // copy num << 162                              "Cavity",                 //name
154                                                << 163                              lWall,                    //mother  volume
                                                   >> 164                              false,                    //no boolean operation
                                                   >> 165                              1);                       //copy number
                                                   >> 166                                 
155   PrintParameters();                              167   PrintParameters();
156                                                << 168     
157   //                                           << 
158   // always return the root volume             << 
159   //                                              169   //
                                                   >> 170   //always return the root volume
                                                   >> 171   //  
160   return fWall;                                   172   return fWall;
161 }                                                 173 }
162                                                   174 
163 //....oooOO0OOooo........oooOO0OOooo........oo    175 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
164                                                   176 
165 void DetectorConstruction::PrintParameters()      177 void DetectorConstruction::PrintParameters()
166 {                                                 178 {
167   G4cout << "\n-------------------------------    179   G4cout << "\n---------------------------------------------------------\n";
168   G4cout << "---> The Wall is " << G4BestUnit( << 180   G4cout << "---> The Wall is " << G4BestUnit(fWallThickness,"Length")
169          << fWallMaterial->GetName() << " ( "  << 181          << " of " << fWallMaterial->GetName() << " ( " 
170          << G4BestUnit(fWallMaterial->GetDensi << 182          << G4BestUnit(fWallMaterial->GetDensity(),"Volumic Mass") << " )\n";
171   G4cout << "     The Cavity is " << G4BestUni << 183   G4cout << "     The Cavity is " << G4BestUnit(fCavityThickness,"Length")
172          << fCavityMaterial->GetName() << " (  << 184          << " of " << fCavityMaterial->GetName() << " ( " 
173          << G4BestUnit(fCavityMaterial->GetDen << 185          << G4BestUnit(fCavityMaterial->GetDensity(),"Volumic Mass") << " )";                  
174   G4cout << "\n-------------------------------    186   G4cout << "\n---------------------------------------------------------\n";
175   G4cout << G4endl;                               187   G4cout << G4endl;
176 }                                                 188 }
177                                                << 189  
178 //....oooOO0OOooo........oooOO0OOooo........oo    190 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
179                                                   191 
180 void DetectorConstruction::SetWallThickness(G4 << 192 void DetectorConstruction::SetWallMaterial(G4String materialChoice)
181 {                                                 193 {
182   fWallThickness = value;                      << 194   // search the material by its name   
                                                   >> 195   G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);     
                                                   >> 196   if (pttoMaterial) fWallMaterial = pttoMaterial;
                                                   >> 197   
                                                   >> 198   pttoMaterial = G4Material::GetMaterial(materialChoice + "_gas");
                                                   >> 199   if (pttoMaterial) fCavityMaterial = pttoMaterial;       
183 }                                                 200 }
184                                                   201 
185 //....oooOO0OOooo........oooOO0OOooo........oo    202 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
186                                                   203 
187 void DetectorConstruction::SetWallMaterial(con << 204 void DetectorConstruction::SetWallThickness(G4double value)
188 {                                                 205 {
189   // search the material by its name           << 206   fWallThickness = value;
190   G4Material* pttoMaterial = G4Material::GetMa << 
191   if (pttoMaterial) fWallMaterial = pttoMateri << 
192 }                                                 207 }
193                                                   208 
194 //....oooOO0OOooo........oooOO0OOooo........oo    209 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
195                                                   210 
196 void DetectorConstruction::SetCavityThickness(    211 void DetectorConstruction::SetCavityThickness(G4double value)
197 {                                                 212 {
198   fCavityThickness = value;                       213   fCavityThickness = value;
                                                   >> 214   
199 }                                                 215 }
200                                                << 
201 //....oooOO0OOooo........oooOO0OOooo........oo    216 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
202                                                   217 
203 void DetectorConstruction::SetCavityRadius(G4d << 218 void DetectorConstruction::SetWorldRadius(G4double value)
204 {                                                 219 {
205   fCavityRadius = value;                       << 220   fWorldRadius = value;
206 }                                                 221 }
207                                                << 222   
208 //....oooOO0OOooo........oooOO0OOooo........oo    223 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
209                                                   224 
210 void DetectorConstruction::SetCavityMaterial(c << 225 #include "G4RunManager.hh" 
                                                   >> 226  
                                                   >> 227 void DetectorConstruction::UpdateGeometry()
211 {                                                 228 {
212   // search the material by its name           << 229 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
213   G4Material* pttoMaterial = G4Material::GetMa << 
214   if (pttoMaterial) fCavityMaterial = pttoMate << 
215 }                                                 230 }
216                                                   231 
217 //....oooOO0OOooo........oooOO0OOooo........oo    232 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
218                                                   233