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