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 8.0.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
 26 /// \file electromagnetic/TestEm12/src/Detecto <<  23 // $Id: DetectorConstruction.cc,v 1.1 2005/07/22 11:08:48 maire Exp $
 27 /// \brief Implementation of the DetectorConst <<  24 // GEANT4 tag $Name: geant4-08-00-patch-01 $
 28 //                                             << 
 29 //                                                 25 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo     26 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     27 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    28 
 33 #include "DetectorConstruction.hh"                 29 #include "DetectorConstruction.hh"
 34                                                << 
 35 #include "DetectorMessenger.hh"                    30 #include "DetectorMessenger.hh"
 36                                                    31 
 37 #include "G4GeometryManager.hh"                << 
 38 #include "G4LogicalVolume.hh"                  << 
 39 #include "G4LogicalVolumeStore.hh"             << 
 40 #include "G4NistManager.hh"                        32 #include "G4NistManager.hh"
                                                   >>  33 #include "G4Sphere.hh"
                                                   >>  34 #include "G4LogicalVolume.hh"
                                                   >>  35 #include "G4VPhysicalVolume.hh"
 41 #include "G4PVPlacement.hh"                        36 #include "G4PVPlacement.hh"
 42 #include "G4PVReplica.hh"                          37 #include "G4PVReplica.hh"
 43 #include "G4PhysicalConstants.hh"              <<  38 #include "G4UniformMagField.hh"
                                                   >>  39 
                                                   >>  40 #include "G4GeometryManager.hh"
 44 #include "G4PhysicalVolumeStore.hh"                41 #include "G4PhysicalVolumeStore.hh"
 45 #include "G4RunManager.hh"                     <<  42 #include "G4LogicalVolumeStore.hh"
 46 #include "G4SolidStore.hh"                         43 #include "G4SolidStore.hh"
 47 #include "G4Sphere.hh"                         <<  44 
 48 #include "G4SystemOfUnits.hh"                  << 
 49 #include "G4UnitsTable.hh"                         45 #include "G4UnitsTable.hh"
 50 #include "G4VPhysicalVolume.hh"                << 
 51                                                    46 
 52 //....oooOO0OOooo........oooOO0OOooo........oo     47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 53                                                    48 
 54 DetectorConstruction::DetectorConstruction()       49 DetectorConstruction::DetectorConstruction()
 55 {                                                  50 {
 56   // default parameter values                      51   // default parameter values
 57   fAbsorRadius = 3 * cm;                       <<  52   absorRadius = 3*cm;
 58   fNbOfLayers = 1;                             <<  53   nbOfLayers = 1;
 59                                                <<  54   
                                                   >>  55   absorMaterial = 0;
                                                   >>  56   magField = 0;
                                                   >>  57   pAbsor   = 0;
                                                   >>  58   
 60   DefineMaterials();                               59   DefineMaterials();
 61   SetMaterial("G4_WATER");                         60   SetMaterial("G4_WATER");
 62                                                    61 
 63   // create commands for interactive definitio <<  62   // create commands for interactive definition of the detector  
 64   fDetectorMessenger = new DetectorMessenger(t <<  63   detectorMessenger = new DetectorMessenger(this);
 65 }                                                  64 }
 66                                                    65 
 67 //....oooOO0OOooo........oooOO0OOooo........oo     66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 68                                                    67 
 69 DetectorConstruction::~DetectorConstruction()      68 DetectorConstruction::~DetectorConstruction()
                                                   >>  69 { delete detectorMessenger;}
                                                   >>  70 
                                                   >>  71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  72 
                                                   >>  73 G4VPhysicalVolume* DetectorConstruction::Construct()
 70 {                                                  74 {
 71   delete fDetectorMessenger;                   <<  75   return ConstructVolumes();
 72 }                                                  76 }
 73                                                    77 
 74 //....oooOO0OOooo........oooOO0OOooo........oo     78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 75                                                    79 
 76 void DetectorConstruction::DefineMaterials()       80 void DetectorConstruction::DefineMaterials()
 77 {                                              <<  81 { 
 78   G4NistManager* man = G4NistManager::Instance     82   G4NistManager* man = G4NistManager::Instance();
 79                                                <<  83   
 80   man->FindOrBuildMaterial("G4_Al");           <<  84   G4bool isotopes = false;
 81   man->FindOrBuildMaterial("G4_Si");           <<  85   
 82   man->FindOrBuildMaterial("G4_Fe");           <<  86   man->FindOrBuildMaterial("G4_Al", isotopes);
 83   man->FindOrBuildMaterial("G4_Ge");           <<  87   man->FindOrBuildMaterial("G4_Si", isotopes);
 84   man->FindOrBuildMaterial("G4_Gd");           <<  88   man->FindOrBuildMaterial("G4_Fe", isotopes);
 85   man->FindOrBuildMaterial("G4_W");            <<  89   man->FindOrBuildMaterial("G4_Ge", isotopes);
 86   man->FindOrBuildMaterial("G4_Pb");           <<  90   man->FindOrBuildMaterial("G4_W" , isotopes);
 87                                                <<  91   man->FindOrBuildMaterial("G4_Pb", isotopes);
 88   man->FindOrBuildMaterial("G4_AIR");          <<  92   
 89   man->FindOrBuildMaterial("G4_WATER");        <<  93   man->FindOrBuildMaterial("G4_AIR"  , isotopes);
 90   man->FindOrBuildMaterial("G4_ALUMINUM_OXIDE" <<  94   man->FindOrBuildMaterial("G4_WATER", isotopes);
 91                                                <<  95 
 92   /// G4cout << *(G4Material::GetMaterialTable <<  96  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
 93 }                                                  97 }
 94                                                    98 
 95 //....oooOO0OOooo........oooOO0OOooo........oo     99 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 96                                                << 100   
 97 G4VPhysicalVolume* DetectorConstruction::Const << 101 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
 98 {                                              << 102 {
                                                   >> 103   G4GeometryManager::GetInstance()->OpenGeometry();
                                                   >> 104   G4PhysicalVolumeStore::GetInstance()->Clean();
                                                   >> 105   G4LogicalVolumeStore::GetInstance()->Clean();
                                                   >> 106   G4SolidStore::GetInstance()->Clean();
                                                   >> 107        
 99   // Absorber                                     108   // Absorber
100   //                                              109   //
101   G4Sphere* sAbsor = new G4Sphere("Absorber",  << 110   G4Sphere* 
102                                   0., fAbsorRa << 111   sAbsor = new G4Sphere("Absorber",     //name
103                                                << 112          0., absorRadius, 0., twopi, 0., pi); //size
104   fSpheres.push_back(sAbsor);                  << 113 
105                                                << 114   G4LogicalVolume*
106   G4LogicalVolume* lAbsor = new G4LogicalVolum << 115   lAbsor = new G4LogicalVolume(sAbsor,      //solid
107                                                << 116                    absorMaterial,   //material
108                                                << 117                   "Absorber");    //name
109   fLVolumes.push_back(lAbsor);                 << 118            
110                                                << 119   pAbsor = new G4PVPlacement(0,       //no rotation
111   fAbsor = new G4PVPlacement(0,  // no rotatio << 120                              G4ThreeVector(),   //at (0,0,0)
112                              G4ThreeVector(),  << 121                              lAbsor,      //logical volume
113                              lAbsor,  // logic << 122                             "Absorber",     //name
114                              "Absorber",  // n << 123                              0,       //mother  volume
115                              0,  // mother  vo << 124                              false,     //no boolean operation
116                              false,  // no boo << 125                              0);      //copy number
117                              0);  // copy numb << 
118                                                   126 
119   // Layers                                       127   // Layers
120   //                                              128   //
121   fLayerThickness = fAbsorRadius / fNbOfLayers << 129   layerThickness = absorRadius/nbOfLayers;
122                                                << 130                         
123   for (G4int i = 1; i <= fNbOfLayers; i++) {   << 131   for (G4int i=1; i<=nbOfLayers; i++) {        
124     G4Sphere* sLayer =                         << 132     G4Sphere*
125       new G4Sphere("Layer", (i - 1) * fLayerTh << 133     sLayer = new G4Sphere("Layer", (i-1)*layerThickness, i*layerThickness,
126                                                << 134                           0., twopi, 0., pi);
127     fSpheres.push_back(sLayer);                << 135      
128                                                << 136     G4LogicalVolume*                                  
129     G4LogicalVolume* lLayer = new G4LogicalVol << 137     lLayer = new G4LogicalVolume(sLayer,    //shape
130                                                << 138                                  absorMaterial,   //material
131                                                << 139                                  "Layer");    //name
132     fLVolumes.push_back(lLayer);               << 140          
133                                                << 141        new G4PVPlacement(0,     //no rotation
134     new G4PVPlacement(0,  // no rotation       << 142                                G4ThreeVector(),   //at (0,0,0)
135                       G4ThreeVector(),  // at  << 143                                lLayer,      //logical volume
136                       lLayer,  // logical volu << 144                                "Layer",     //name
137                       "Layer",  // name        << 145                                lAbsor,      //mother  volume
138                       lAbsor,  // mother  volu << 146                                false,     //no boolean operation
139                       false,  // no boolean op << 147                                i);      //copy number
140                       i);  // copy number      << 148                                  
141   }                                            << 149    }               
142                                                   150 
143   PrintParameters();                              151   PrintParameters();
144                                                << 152     
145   //                                           << 
146   // always return the root volume             << 
147   //                                              153   //
148   return fAbsor;                               << 154   //always return the root volume
                                                   >> 155   //  
                                                   >> 156   return pAbsor;
149 }                                                 157 }
150                                                   158 
151 //....oooOO0OOooo........oooOO0OOooo........oo    159 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
152                                                   160 
153 void DetectorConstruction::PrintParameters()      161 void DetectorConstruction::PrintParameters()
154 {                                                 162 {
155   G4cout << "\n-------------------------------    163   G4cout << "\n---------------------------------------------------------\n";
156   G4cout << "---> The Absorber is a sphere of  << 164   G4cout << "---> The Absorber is a sphere of " 
157          << " radius of " << fAbsorMaterial->G << 165          << G4BestUnit(absorRadius,"Length") << " radius of "
158          << " slices of " << G4BestUnit(fLayer << 166          << absorMaterial->GetName() << " divided in " << nbOfLayers 
159          << fAbsorMaterial << G4endl;          << 167    << " slices of " << G4BestUnit(layerThickness,"Length") << G4endl;
160   G4cout << "\n-------------------------------    168   G4cout << "\n---------------------------------------------------------\n";
161 }                                                 169 }
162                                                   170 
163 //....oooOO0OOooo........oooOO0OOooo........oo    171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
164                                                   172 
165 void DetectorConstruction::SetRadius(G4double     173 void DetectorConstruction::SetRadius(G4double value)
166 {                                                 174 {
167   // geometry was already constructed - scale  << 175   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 }                                                 176 }
183                                                   177 
184 //....oooOO0OOooo........oooOO0OOooo........oo    178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
185                                                   179 
186 void DetectorConstruction::SetMaterial(G4Strin    180 void DetectorConstruction::SetMaterial(G4String materialChoice)
187 {                                                 181 {
188   // search the material by its name           << 182   // search the material by its name   
189   G4Material* pttoMaterial = G4Material::GetMa << 183   G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);     
                                                   >> 184   if (pttoMaterial) absorMaterial = pttoMaterial;
                                                   >> 185 }
190                                                   186 
191   if (pttoMaterial && pttoMaterial != fAbsorMa << 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
192     fAbsorMaterial = pttoMaterial;             << 188 
193     G4RunManager::GetRunManager()->PhysicsHasB << 189 void DetectorConstruction::SetNbOfLayers(G4int value)
194                                                << 190 {
195     // geometry was already constructed - only << 191   nbOfLayers = value; 
196     if (fAbsor) {                              << 192 }
197       for (auto lv : fLVolumes) {              << 193  
198         lv->SetMaterial(fAbsorMaterial);       << 194 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
199       }                                        << 195 
                                                   >> 196 #include "G4FieldManager.hh"
                                                   >> 197 #include "G4TransportationManager.hh"
                                                   >> 198 
                                                   >> 199 void DetectorConstruction::SetMagField(G4double fieldValue)
                                                   >> 200 {
                                                   >> 201   //apply a global uniform magnetic field along Z axis
                                                   >> 202   G4FieldManager* fieldMgr 
                                                   >> 203    = G4TransportationManager::GetTransportationManager()->GetFieldManager();
                                                   >> 204     
                                                   >> 205   if (magField) delete magField;  //delete the existing magn field
                                                   >> 206   
                                                   >> 207   if (fieldValue!=0.)     // create a new one if non nul
                                                   >> 208     {
                                                   >> 209       magField = new G4UniformMagField(G4ThreeVector(0.,0.,fieldValue));        
                                                   >> 210       fieldMgr->SetDetectorField(magField);
                                                   >> 211       fieldMgr->CreateChordFinder(magField);
                                                   >> 212     }
                                                   >> 213    else
                                                   >> 214     {
                                                   >> 215       magField = 0;
                                                   >> 216       fieldMgr->SetDetectorField(magField);
200     }                                             217     }
201   }                                            << 
202 }                                                 218 }
203                                                   219 
204 //....oooOO0OOooo........oooOO0OOooo........oo    220 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
205                                                   221 
206 void DetectorConstruction::SetNbOfLayers(G4int << 222 #include "G4RunManager.hh" 
                                                   >> 223  
                                                   >> 224 void DetectorConstruction::UpdateGeometry()
207 {                                                 225 {
208   fNbOfLayers = value;                         << 226 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
209 }                                                 227 }
210                                                   228 
211 //....oooOO0OOooo........oooOO0OOooo........oo    229 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
212                                                   230