Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/polarisation/Pol01/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/polarisation/Pol01/src/DetectorConstruction.cc (Version 11.3.0) and /examples/extended/polarisation/Pol01/src/DetectorConstruction.cc (Version 10.3)


  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 polarisation/Pol01/src/DetectorConst     26 /// \file polarisation/Pol01/src/DetectorConstruction.cc
 27 /// \brief Implementation of the DetectorConst     27 /// \brief Implementation of the DetectorConstruction class
 28 //                                                 28 //
                                                   >>  29 // $Id: DetectorConstruction.cc 98772 2016-08-09 14:25:31Z gcosmo $
 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 "G4Material.hh"
 37 #include "G4Box.hh"                                38 #include "G4Box.hh"
 38 #include "G4GeometryManager.hh"                << 
 39 #include "G4LogicalVolume.hh"                      39 #include "G4LogicalVolume.hh"
 40 #include "G4LogicalVolumeStore.hh"             << 
 41 #include "G4Material.hh"                       << 
 42 #include "G4NistManager.hh"                    << 
 43 #include "G4PVPlacement.hh"                        40 #include "G4PVPlacement.hh"
                                                   >>  41 
                                                   >>  42 #include "G4GeometryManager.hh"
 44 #include "G4PhysicalVolumeStore.hh"                43 #include "G4PhysicalVolumeStore.hh"
 45 #include "G4PolarizationManager.hh"            <<  44 #include "G4LogicalVolumeStore.hh"
 46 #include "G4RunManager.hh"                     << 
 47 #include "G4SolidStore.hh"                         45 #include "G4SolidStore.hh"
 48 #include "G4SystemOfUnits.hh"                  <<  46 
 49 #include "G4UnitsTable.hh"                         47 #include "G4UnitsTable.hh"
 50                                                    48 
                                                   >>  49 #include "G4PolarizationManager.hh"
                                                   >>  50 #include "G4NistManager.hh"
                                                   >>  51 #include "G4SystemOfUnits.hh"
                                                   >>  52 
 51 //....oooOO0OOooo........oooOO0OOooo........oo     53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 52                                                    54 
 53 DetectorConstruction::DetectorConstruction()       55 DetectorConstruction::DetectorConstruction()
 54   : G4VUserDetectorConstruction(), fWorld(0),  <<  56 : G4VUserDetectorConstruction(),
                                                   >>  57   fWorld(0), fBox(0), fTargetMaterial(0), fWorldMaterial(0)
 55 {                                                  58 {
 56   fBoxSizeXY = 50 * mm;                        <<  59   fBoxSizeXY = 50*mm;
 57   fBoxSizeZ = 5 * mm;                          <<  60   fBoxSizeZ = 5*mm;
 58   fWorldSize = 1. * m;                         <<  61   fWorldSize = 1.*m;
 59   ConstructMaterials();                        <<  62   SetTargetMaterial("G4_Fe");  
                                                   >>  63   SetWorldMaterial("G4_Galactic");  
 60   fMessenger = new DetectorMessenger(this);        64   fMessenger = new DetectorMessenger(this);
 61 }                                                  65 }
 62                                                    66 
 63 //....oooOO0OOooo........oooOO0OOooo........oo     67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 64                                                    68 
 65 DetectorConstruction::~DetectorConstruction()      69 DetectorConstruction::~DetectorConstruction()
 66 {                                              <<  70 { 
 67   delete fMessenger;                               71   delete fMessenger;
 68 }                                                  72 }
 69 //....oooOO0OOooo........oooOO0OOooo........oo << 
 70                                                << 
 71 void DetectorConstruction::ConstructMaterials( << 
 72 {                                              << 
 73   auto nistManager = G4NistManager::Instance() << 
 74                                                << 
 75   fTargetMaterial = nistManager->FindOrBuildMa << 
 76   if (fTargetMaterial == nullptr) {            << 
 77     G4cerr << "### ERROR - Material: <"        << 
 78            << "G4_Fe"                          << 
 79            << "> not found" << G4endl;         << 
 80   }                                            << 
 81                                                << 
 82   fWorldMaterial = nistManager->FindOrBuildMat << 
 83   if (fWorldMaterial == nullptr) {             << 
 84     G4cerr << "### ERROR -  Material: <"       << 
 85            << "G4_Galactic"                    << 
 86            << "> not found" << G4endl;         << 
 87   }                                            << 
 88 }                                              << 
 89                                                    73 
 90 //....oooOO0OOooo........oooOO0OOooo........oo     74 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 91                                                    75 
 92 G4LogicalVolume* lBox;                         << 
 93                                                << 
 94 G4VPhysicalVolume* DetectorConstruction::Const     76 G4VPhysicalVolume* DetectorConstruction::Construct()
 95 {                                                  77 {
 96   // G4PolarizationManager::GetInstance()->Cle <<  78   // Cleanup old geometry
                                                   >>  79   G4GeometryManager::GetInstance()->OpenGeometry();
                                                   >>  80   G4PhysicalVolumeStore::GetInstance()->Clean();
                                                   >>  81   G4LogicalVolumeStore::GetInstance()->Clean();
                                                   >>  82   G4SolidStore::GetInstance()->Clean();
                                                   >>  83   
                                                   >>  84   G4PolarizationManager::GetInstance()->Clean();
 97                                                    85 
 98   // World                                         86   // World
 99   //                                               87   //
100   G4Box* sWorld = new G4Box("World",  // name  <<  88   G4Box*
101                             fWorldSize / 2, fW <<  89   sWorld = new G4Box("World",                            //name
102                                                <<  90                    fWorldSize/2,fWorldSize/2,fWorldSize/2); //dimensions
103   G4LogicalVolume* lWorld = new G4LogicalVolum <<  91 
104                                                <<  92   G4LogicalVolume*                                                                 
105                                                <<  93   lWorld = new G4LogicalVolume(sWorld,                   //shape
106                                                <<  94                                fWorldMaterial,           //material
107   fWorld = new G4PVPlacement(0,  // no rotatio <<  95                               "World");                  //name
108                              G4ThreeVector(),  <<  96 
109                              lWorld,  // logic <<  97   fWorld = new G4PVPlacement(0,                          //no rotation
110                              "World",  // name <<  98                              G4ThreeVector(),            //at (0,0,0)
111                              0,  // mother vol <<  99                              lWorld,                     //logical volume
112                              false,  // no boo << 100                              "World",                    //name
113                              0);  // copy numb << 101                              0,                          //mother volume
114                                                << 102                              false,                      //no boolean operation
                                                   >> 103                              0);                         //copy number
                                                   >> 104                                                       
115   // Box                                          105   // Box
116   //                                           << 106   //                           
117   G4Box* sBox = new G4Box("Container",  // its << 107   G4Box*
118                           fBoxSizeXY / 2., fBo << 108   sBox = new G4Box("Container",                           //its name
119                                                << 109                    fBoxSizeXY/2.,fBoxSizeXY/2.,fBoxSizeZ/2.);//its dimensions
120   // G4LogicalVolume*                          << 110                    
121   lBox = new G4LogicalVolume(sBox,  // its sha << 111   G4LogicalVolume*
122                              fTargetMaterial,  << 112   lBox = new G4LogicalVolume(sBox,                        //its shape
123                              "theBox");  // it << 113                              fTargetMaterial,             //its material
124                                                << 114                              "theBox");                   //its name
125   fBox = new G4PVPlacement(0,  // no rotation  << 115 
126                            G4ThreeVector(),  / << 116   fBox = new G4PVPlacement(0,                             //no rotation
127                            lBox,  // its logic << 117                            G4ThreeVector(),               //at (0,0,0)
128                            fTargetMaterial->Ge << 118                            lBox,                          //its logical volume
129                            lWorld,  // its mot << 119                            fTargetMaterial->GetName(),    //its name
130                            false,  // no boole << 120                            lWorld,                        //its mother  volume
131                            0);  // copy number << 121                            false,                         //no boolean operation
                                                   >> 122                            0);                            //copy number
132                                                   123 
                                                   >> 124   // register logical Volume in PolarizationManager with zero polarization
                                                   >> 125   G4PolarizationManager * polMgr = G4PolarizationManager::GetInstance();
                                                   >> 126   polMgr->SetVolumePolarization(lBox,G4ThreeVector(0.,0.,0.));
                                                   >> 127                            
133   PrintParameters();                              128   PrintParameters();
134                                                << 129   
135   // always return the root volume             << 130   //always return the root volume
136   //                                              131   //
137   return fWorld;                                  132   return fWorld;
138 }                                                 133 }
139                                                   134 
140 //....oooOO0OOooo........oooOO0OOooo........oo    135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141                                                   136 
142 void DetectorConstruction::ConstructSDandField << 
143 {                                              << 
144   // register logical Volume in PolarizationMa << 
145   G4PolarizationManager* polMgr = G4Polarizati << 
146   polMgr->SetVolumePolarization(lBox, G4ThreeV << 
147 }                                              << 
148                                                << 
149 //....oooOO0OOooo........oooOO0OOooo........oo << 
150                                                << 
151 void DetectorConstruction::PrintParameters()      137 void DetectorConstruction::PrintParameters()
152 {                                                 138 {
153   G4cout << "\n The Box is " << G4BestUnit(fBo << 139   G4cout << "\n The Box is " << G4BestUnit(fBoxSizeXY,"Length")
154          << G4BestUnit(fBoxSizeXY, "Length") < << 140          << " x " << G4BestUnit(fBoxSizeXY,"Length")
155          << fTargetMaterial->GetName() << G4en << 141          << " x " << G4BestUnit(fBoxSizeZ,"Length")
                                                   >> 142          << " of " << fTargetMaterial->GetName() << G4endl;
156 }                                                 143 }
157                                                   144 
158 //....oooOO0OOooo........oooOO0OOooo........oo    145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
159                                                   146 
160 void DetectorConstruction::SetTargetMaterial(G    147 void DetectorConstruction::SetTargetMaterial(G4String materialChoice)
161 {                                                 148 {
162   // search the material by its name              149   // search the material by its name
163   G4Material* mat = G4NistManager::Instance()- << 150   G4Material* mat =
                                                   >> 151     G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
164   if (mat != fTargetMaterial) {                   152   if (mat != fTargetMaterial) {
165     if (mat) {                                 << 153     if(mat) {
166       fTargetMaterial = mat;                      154       fTargetMaterial = mat;
167       G4RunManager::GetRunManager()->PhysicsHa << 
168       UpdateGeometry();                           155       UpdateGeometry();
169     }                                          << 156     } else {
170     else {                                     << 157       G4cout << "### Warning!  Target material: <"
171       G4cout << "### Warning!  Target material << 158            << materialChoice << "> not found" << G4endl;  
172     }                                          << 159     }     
173   }                                               160   }
174 }                                                 161 }
175                                                   162 
176 //....oooOO0OOooo........oooOO0OOooo........oo    163 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
177                                                   164 
178 void DetectorConstruction::SetWorldMaterial(G4    165 void DetectorConstruction::SetWorldMaterial(G4String materialChoice)
179 {                                                 166 {
180   // search the material by its name              167   // search the material by its name
181   G4Material* mat = G4NistManager::Instance()- << 168   G4Material* mat =
                                                   >> 169     G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
182   if (mat != fWorldMaterial) {                    170   if (mat != fWorldMaterial) {
183     if (mat) {                                 << 171     if(mat) {
184       fWorldMaterial = mat;                       172       fWorldMaterial = mat;
185       G4RunManager::GetRunManager()->PhysicsHa << 
186       UpdateGeometry();                           173       UpdateGeometry();
187     }                                          << 174     } else {
188     else {                                     << 175       G4cout << "### Warning! World material: <"
189       G4cout << "### Warning! World material:  << 176            << materialChoice << "> not found" << G4endl;  
190     }                                          << 177     }     
191   }                                               178   }
192 }                                                 179 }
193                                                   180 
194 //....oooOO0OOooo........oooOO0OOooo........oo    181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
195                                                   182 
196 void DetectorConstruction::SetSizeXY(G4double     183 void DetectorConstruction::SetSizeXY(G4double value)
197 {                                                 184 {
198   fBoxSizeXY = value;                          << 185   fBoxSizeXY = value; 
199   if (fWorldSize < fBoxSizeXY) fWorldSize = 1. << 186   if (fWorldSize<fBoxSizeXY) fWorldSize = 1.2*fBoxSizeXY;
200   UpdateGeometry();                               187   UpdateGeometry();
201 }                                                 188 }
202                                                   189 
203 void DetectorConstruction::SetSizeZ(G4double v    190 void DetectorConstruction::SetSizeZ(G4double value)
204 {                                                 191 {
205   fBoxSizeZ = value;                           << 192   fBoxSizeZ = value; 
206   if (fWorldSize < fBoxSizeZ) fWorldSize = 1.2 << 193   if (fWorldSize<fBoxSizeZ) fWorldSize = 1.2*fBoxSizeZ;
207   UpdateGeometry();                               194   UpdateGeometry();
208 }                                                 195 }
209                                                   196 
210 //....oooOO0OOooo........oooOO0OOooo........oo    197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
211                                                   198 
212 #include "G4RunManager.hh"                        199 #include "G4RunManager.hh"
213                                                   200 
214 void DetectorConstruction::UpdateGeometry()       201 void DetectorConstruction::UpdateGeometry()
215 {                                                 202 {
216   // if (fWorld)                               << 203   if (fWorld) 
217   //   G4RunManager::GetRunManager()->DefineWo << 204     G4RunManager::GetRunManager()->DefineWorldVolume(Construct());
218   G4RunManager::GetRunManager()->GeometryHasBe << 
219 }                                                 205 }
220                                                   206 
221 //....oooOO0OOooo........oooOO0OOooo........oo    207 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
222                                                   208