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 9.6.p4)


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