Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm12/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/TestEm12/src/DetectorConstruction.cc (Version 11.3.0) and /examples/extended/electromagnetic/TestEm12/src/DetectorConstruction.cc (Version 11.2)


  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/TestEm12/src/Detecto     26 /// \file electromagnetic/TestEm12/src/DetectorConstruction.cc
 27 /// \brief Implementation of the DetectorConst     27 /// \brief Implementation of the DetectorConstruction class
 28 //                                                 28 //
 29 //                                                 29 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo     30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    32 
 33 #include "DetectorConstruction.hh"                 33 #include "DetectorConstruction.hh"
 34                                                << 
 35 #include "DetectorMessenger.hh"                    34 #include "DetectorMessenger.hh"
 36                                                    35 
 37 #include "G4GeometryManager.hh"                << 
 38 #include "G4LogicalVolume.hh"                  << 
 39 #include "G4LogicalVolumeStore.hh"             << 
 40 #include "G4NistManager.hh"                        36 #include "G4NistManager.hh"
                                                   >>  37 #include "G4Sphere.hh"
                                                   >>  38 #include "G4LogicalVolume.hh"
                                                   >>  39 #include "G4VPhysicalVolume.hh"
 41 #include "G4PVPlacement.hh"                        40 #include "G4PVPlacement.hh"
 42 #include "G4PVReplica.hh"                          41 #include "G4PVReplica.hh"
 43 #include "G4PhysicalConstants.hh"              <<  42 
                                                   >>  43 #include "G4GeometryManager.hh"
 44 #include "G4PhysicalVolumeStore.hh"                44 #include "G4PhysicalVolumeStore.hh"
 45 #include "G4RunManager.hh"                     <<  45 #include "G4LogicalVolumeStore.hh"
 46 #include "G4SolidStore.hh"                         46 #include "G4SolidStore.hh"
 47 #include "G4Sphere.hh"                         <<  47 #include "G4RunManager.hh"
 48 #include "G4SystemOfUnits.hh"                  <<  48 
 49 #include "G4UnitsTable.hh"                         49 #include "G4UnitsTable.hh"
 50 #include "G4VPhysicalVolume.hh"                <<  50 #include "G4PhysicalConstants.hh"
                                                   >>  51 #include "G4SystemOfUnits.hh"
 51                                                    52 
 52 //....oooOO0OOooo........oooOO0OOooo........oo     53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 53                                                    54 
 54 DetectorConstruction::DetectorConstruction()       55 DetectorConstruction::DetectorConstruction()
 55 {                                                  56 {
 56   // default parameter values                      57   // default parameter values
 57   fAbsorRadius = 3 * cm;                       <<  58   fAbsorRadius = 3*cm;
 58   fNbOfLayers = 1;                                 59   fNbOfLayers = 1;
 59                                                <<  60   
 60   DefineMaterials();                               61   DefineMaterials();
 61   SetMaterial("G4_WATER");                         62   SetMaterial("G4_WATER");
 62                                                    63 
 63   // create commands for interactive definitio <<  64   // create commands for interactive definition of the detector  
 64   fDetectorMessenger = new DetectorMessenger(t     65   fDetectorMessenger = new DetectorMessenger(this);
 65 }                                                  66 }
 66                                                    67 
 67 //....oooOO0OOooo........oooOO0OOooo........oo     68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 68                                                    69 
 69 DetectorConstruction::~DetectorConstruction()      70 DetectorConstruction::~DetectorConstruction()
 70 {                                              <<  71 { delete fDetectorMessenger;}
 71   delete fDetectorMessenger;                   << 
 72 }                                              << 
 73                                                    72 
 74 //....oooOO0OOooo........oooOO0OOooo........oo     73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 75                                                    74 
 76 void DetectorConstruction::DefineMaterials()       75 void DetectorConstruction::DefineMaterials()
 77 {                                              <<  76 { 
 78   G4NistManager* man = G4NistManager::Instance     77   G4NistManager* man = G4NistManager::Instance();
 79                                                <<  78   
 80   man->FindOrBuildMaterial("G4_Al");               79   man->FindOrBuildMaterial("G4_Al");
 81   man->FindOrBuildMaterial("G4_Si");               80   man->FindOrBuildMaterial("G4_Si");
 82   man->FindOrBuildMaterial("G4_Fe");               81   man->FindOrBuildMaterial("G4_Fe");
 83   man->FindOrBuildMaterial("G4_Ge");               82   man->FindOrBuildMaterial("G4_Ge");
 84   man->FindOrBuildMaterial("G4_Gd");               83   man->FindOrBuildMaterial("G4_Gd");
 85   man->FindOrBuildMaterial("G4_W");                84   man->FindOrBuildMaterial("G4_W");
 86   man->FindOrBuildMaterial("G4_Pb");               85   man->FindOrBuildMaterial("G4_Pb");
 87                                                <<  86   
 88   man->FindOrBuildMaterial("G4_AIR");              87   man->FindOrBuildMaterial("G4_AIR");
 89   man->FindOrBuildMaterial("G4_WATER");            88   man->FindOrBuildMaterial("G4_WATER");
 90   man->FindOrBuildMaterial("G4_ALUMINUM_OXIDE"     89   man->FindOrBuildMaterial("G4_ALUMINUM_OXIDE");
 91                                                <<  90   
 92   /// G4cout << *(G4Material::GetMaterialTable <<  91   ///G4cout << *(G4Material::GetMaterialTable()) << G4endl;
 93 }                                                  92 }
 94                                                    93 
 95 //....oooOO0OOooo........oooOO0OOooo........oo     94 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 96                                                <<  95   
 97 G4VPhysicalVolume* DetectorConstruction::Const     96 G4VPhysicalVolume* DetectorConstruction::Construct()
 98 {                                                  97 {
 99   // Absorber                                      98   // Absorber
100   //                                               99   //
101   G4Sphere* sAbsor = new G4Sphere("Absorber",  << 100   G4Sphere* 
102                                   0., fAbsorRa << 101   sAbsor = new G4Sphere("Absorber",                           //name
                                                   >> 102                      0., fAbsorRadius, 0., twopi, 0., pi);    //size
103                                                   103 
104   fSpheres.push_back(sAbsor);                     104   fSpheres.push_back(sAbsor);
105                                                   105 
106   G4LogicalVolume* lAbsor = new G4LogicalVolum << 106   G4LogicalVolume*
107                                                << 107   lAbsor = new G4LogicalVolume(sAbsor,                        //solid
108                                                << 108                                      fAbsorMaterial,          //material
                                                   >> 109                                     "Absorber");              //name
109   fLVolumes.push_back(lAbsor);                    110   fLVolumes.push_back(lAbsor);
110                                                << 111                                    
111   fAbsor = new G4PVPlacement(0,  // no rotatio << 112   fAbsor = new G4PVPlacement(0,                         //no rotation
112                              G4ThreeVector(),  << 113                              G4ThreeVector(),           //at (0,0,0)
113                              lAbsor,  // logic << 114                              lAbsor,                    //logical volume
114                              "Absorber",  // n << 115                             "Absorber",                 //name
115                              0,  // mother  vo << 116                              0,                         //mother  volume
116                              false,  // no boo << 117                              false,                     //no boolean operation
117                              0);  // copy numb << 118                              0);                        //copy number
118                                                   119 
119   // Layers                                       120   // Layers
120   //                                              121   //
121   fLayerThickness = fAbsorRadius / fNbOfLayers << 122   fLayerThickness = fAbsorRadius/fNbOfLayers;
122                                                << 123                         
123   for (G4int i = 1; i <= fNbOfLayers; i++) {   << 124   for (G4int i=1; i<=fNbOfLayers; i++) {                           
124     G4Sphere* sLayer =                         << 125     G4Sphere*
125       new G4Sphere("Layer", (i - 1) * fLayerTh << 126     sLayer = new G4Sphere("Layer", (i-1)*fLayerThickness, i*fLayerThickness,
                                                   >> 127                           0., twopi, 0., pi);
126                                                   128 
127     fSpheres.push_back(sLayer);                   129     fSpheres.push_back(sLayer);
128                                                << 130                  
129     G4LogicalVolume* lLayer = new G4LogicalVol << 131     G4LogicalVolume* 
130                                                << 132     lLayer = new G4LogicalVolume(sLayer,                //shape
131                                                << 133                                  fAbsorMaterial,        //material
                                                   >> 134                                  "Layer");              //name
132     fLVolumes.push_back(lLayer);                  135     fLVolumes.push_back(lLayer);
133                                                << 136                                  
134     new G4PVPlacement(0,  // no rotation       << 137              new G4PVPlacement(0,                      //no rotation
135                       G4ThreeVector(),  // at  << 138                                G4ThreeVector(),        //at (0,0,0)
136                       lLayer,  // logical volu << 139                                lLayer,                 //logical volume
137                       "Layer",  // name        << 140                                "Layer",                //name
138                       lAbsor,  // mother  volu << 141                                lAbsor,                 //mother  volume
139                       false,  // no boolean op << 142                                false,                  //no boolean operation
140                       i);  // copy number      << 143                                i);                     //copy number
141   }                                            << 144                                                    
                                                   >> 145    }                                              
142                                                   146 
143   PrintParameters();                              147   PrintParameters();
144                                                << 148     
145   //                                           << 
146   // always return the root volume             << 
147   //                                              149   //
                                                   >> 150   //always return the root volume
                                                   >> 151   //  
148   return fAbsor;                                  152   return fAbsor;
149 }                                                 153 }
150                                                   154 
151 //....oooOO0OOooo........oooOO0OOooo........oo    155 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
152                                                   156 
153 void DetectorConstruction::PrintParameters()      157 void DetectorConstruction::PrintParameters()
154 {                                                 158 {
155   G4cout << "\n-------------------------------    159   G4cout << "\n---------------------------------------------------------\n";
156   G4cout << "---> The Absorber is a sphere of  << 160   G4cout << "---> The Absorber is a sphere of " 
157          << " radius of " << fAbsorMaterial->G << 161          << G4BestUnit(fAbsorRadius,"Length") << " radius of "
158          << " slices of " << G4BestUnit(fLayer << 162          << fAbsorMaterial->GetName() << " divided in " << fNbOfLayers 
159          << fAbsorMaterial << G4endl;          << 163          << " slices of " << G4BestUnit(fLayerThickness,"Length")
                                                   >> 164          << "\n \n" << fAbsorMaterial << G4endl;
160   G4cout << "\n-------------------------------    165   G4cout << "\n---------------------------------------------------------\n";
161 }                                                 166 }
162                                                   167 
163 //....oooOO0OOooo........oooOO0OOooo........oo    168 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
164                                                   169 
165 void DetectorConstruction::SetRadius(G4double     170 void DetectorConstruction::SetRadius(G4double value)
166 {                                                 171 {
167   // geometry was already constructed - scale     172   // geometry was already constructed - scale radii of all spheres
168   if (fAbsor) {                                << 173   if(fAbsor) {
169     G4double scale = value / fAbsorRadius;     << 174     G4double scale = value/fAbsorRadius;
170     for (auto solid : fSpheres) {              << 175     for (auto solid : fSpheres) { 
171       if (scale > 1.0) {                       << 176       if(scale > 1.0) {
172         solid->SetOuterRadius(solid->GetOuterR << 177         solid->SetOuterRadius(solid->GetOuterRadius()*scale); 
173         solid->SetInnerRadius(solid->GetInnerR << 178         solid->SetInnerRadius(solid->GetInnerRadius()*scale);
174       }                                        << 179       } else { 
175       else {                                   << 180         solid->SetInnerRadius(solid->GetInnerRadius()*scale);
176         solid->SetInnerRadius(solid->GetInnerR << 181         solid->SetOuterRadius(solid->GetOuterRadius()*scale); 
177         solid->SetOuterRadius(solid->GetOuterR << 
178       }                                           182       }
179     }                                             183     }
180   }                                               184   }
181   fAbsorRadius = value;                           185   fAbsorRadius = value;
182 }                                                 186 }
183                                                   187 
184 //....oooOO0OOooo........oooOO0OOooo........oo    188 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
185                                                   189 
186 void DetectorConstruction::SetMaterial(G4Strin    190 void DetectorConstruction::SetMaterial(G4String materialChoice)
187 {                                                 191 {
188   // search the material by its name           << 192   // search the material by its name   
189   G4Material* pttoMaterial = G4Material::GetMa << 193   G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);     
190                                                   194 
191   if (pttoMaterial && pttoMaterial != fAbsorMa    195   if (pttoMaterial && pttoMaterial != fAbsorMaterial) {
192     fAbsorMaterial = pttoMaterial;                196     fAbsorMaterial = pttoMaterial;
193     G4RunManager::GetRunManager()->PhysicsHasB    197     G4RunManager::GetRunManager()->PhysicsHasBeenModified();
194                                                   198 
195     // geometry was already constructed - only    199     // geometry was already constructed - only change material
196     if (fAbsor) {                              << 200     if(fAbsor) {
197       for (auto lv : fLVolumes) {              << 201       for (auto lv : fLVolumes) { lv->SetMaterial(fAbsorMaterial); }
198         lv->SetMaterial(fAbsorMaterial);       << 
199       }                                        << 
200     }                                             202     }
201   }                                               203   }
202 }                                                 204 }
203                                                   205 
204 //....oooOO0OOooo........oooOO0OOooo........oo    206 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
205                                                   207 
206 void DetectorConstruction::SetNbOfLayers(G4int    208 void DetectorConstruction::SetNbOfLayers(G4int value)
207 {                                                 209 {
208   fNbOfLayers = value;                            210   fNbOfLayers = value;
209 }                                                 211 }
210                                                << 212  
211 //....oooOO0OOooo........oooOO0OOooo........oo    213 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
212                                                   214