Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/hadronic/Hadr03/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/hadronic/Hadr03/src/DetectorConstruction.cc (Version 11.3.0) and /examples/extended/hadronic/Hadr03/src/DetectorConstruction.cc (Version 10.6.p3)


  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 DetectorConstruction.cc                  26 /// \file DetectorConstruction.cc
 27 /// \brief Implementation of the DetectorConst     27 /// \brief Implementation of the DetectorConstruction class
 28 //                                                 28 //
 29 //                                                 29 //
 30                                                    30 
 31 //....oooOO0OOooo........oooOO0OOooo........oo     31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32 //....oooOO0OOooo........oooOO0OOooo........oo     32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 33                                                    33 
 34 #include "DetectorConstruction.hh"                 34 #include "DetectorConstruction.hh"
 35                                                << 
 36 #include "DetectorMessenger.hh"                    35 #include "DetectorMessenger.hh"
 37                                                    36 
 38 #include "G4Box.hh"                            << 
 39 #include "G4GeometryManager.hh"                << 
 40 #include "G4LogicalVolume.hh"                  << 
 41 #include "G4LogicalVolumeStore.hh"             << 
 42 #include "G4Material.hh"                           37 #include "G4Material.hh"
 43 #include "G4NistManager.hh"                        38 #include "G4NistManager.hh"
                                                   >>  39 #include "G4Box.hh"
                                                   >>  40 #include "G4LogicalVolume.hh"
 44 #include "G4PVPlacement.hh"                        41 #include "G4PVPlacement.hh"
                                                   >>  42 
                                                   >>  43 #include "G4GeometryManager.hh"
 45 #include "G4PhysicalVolumeStore.hh"                44 #include "G4PhysicalVolumeStore.hh"
 46 #include "G4RunManager.hh"                     <<  45 #include "G4LogicalVolumeStore.hh"
 47 #include "G4SolidStore.hh"                         46 #include "G4SolidStore.hh"
 48 #include "G4SystemOfUnits.hh"                  <<  47 
 49 #include "G4UnitsTable.hh"                         48 #include "G4UnitsTable.hh"
                                                   >>  49 #include "G4SystemOfUnits.hh"
                                                   >>  50 
                                                   >>  51 #include "G4RunManager.hh"
 50                                                    52 
 51 //....oooOO0OOooo........oooOO0OOooo........oo     53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 52                                                    54 
 53 DetectorConstruction::DetectorConstruction()       55 DetectorConstruction::DetectorConstruction()
                                                   >>  56 :G4VUserDetectorConstruction(),
                                                   >>  57  fPBox(0), fLBox(0), fMaterial(0), fDetectorMessenger(0)
 54 {                                                  58 {
 55   fBoxSize = 10 * m;                           <<  59   fBoxSize = 10*m;
 56   DefineMaterials();                               60   DefineMaterials();
 57   SetMaterial("Molybdenum98");                 <<  61   SetMaterial("Molybdenum98");  
 58   fDetectorMessenger = new DetectorMessenger(t     62   fDetectorMessenger = new DetectorMessenger(this);
 59 }                                                  63 }
 60                                                    64 
 61 //....oooOO0OOooo........oooOO0OOooo........oo     65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 62                                                    66 
 63 DetectorConstruction::~DetectorConstruction()      67 DetectorConstruction::~DetectorConstruction()
 64 {                                              <<  68 { delete fDetectorMessenger;}
 65   delete fDetectorMessenger;                   << 
 66 }                                              << 
 67                                                    69 
 68 //....oooOO0OOooo........oooOO0OOooo........oo     70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 69                                                    71 
 70 G4VPhysicalVolume* DetectorConstruction::Const     72 G4VPhysicalVolume* DetectorConstruction::Construct()
 71 {                                                  73 {
 72   return ConstructVolumes();                       74   return ConstructVolumes();
 73 }                                                  75 }
 74                                                    76 
 75 //....oooOO0OOooo........oooOO0OOooo........oo     77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 76                                                    78 
 77 void DetectorConstruction::DefineMaterials()       79 void DetectorConstruction::DefineMaterials()
 78 {                                                  80 {
 79   // define a Material from isotopes           <<  81  // define a Material from isotopes
 80   //                                           <<  82  //
 81   MaterialWithSingleIsotope("Molybdenum98", "M <<  83  MaterialWithSingleIsotope("Molybdenum98", "Mo98",  10.28*g/cm3, 42, 98);
 82                                                <<  84  
 83   // NE213                                     <<  85  //NE213
 84   G4Element* H = new G4Element("Hydrogen", "H" <<  86  G4Element* H  = new G4Element("Hydrogen" ,"H" , 1.,  1.01*g/mole);
 85   G4Element* C = new G4Element("Carbon", "C",  <<  87  G4Element* C  = new G4Element("Carbon"   ,"C" , 6., 12.00*g/mole);
 86   G4Material* ne213 = new G4Material("NE213",  <<  88  G4Material* ne213 = 
 87   ne213->AddElement(H, 9.2 * perCent);         <<  89  new G4Material("NE213", 0.874*g/cm3, 2);
 88   ne213->AddElement(C, 90.8 * perCent);        <<  90  ne213->AddElement(H,    9.2*perCent);
 89                                                <<  91  ne213->AddElement(C,   90.8*perCent);
 90   G4Material* hydrogen = new G4Material("hydro <<  92  
 91   hydrogen->AddElement(H, 1);                  <<  93  G4Material* hydrogen = 
 92                                                <<  94  new G4Material("hydrogen", 1.0*g/cm3, 1);
 93   G4Material* carbon = new G4Material("carbon" <<  95  hydrogen->AddElement(H, 1);
 94   carbon->AddElement(C, 1);                    <<  96  
 95                                                <<  97  G4Material* carbon = 
 96   G4Material* plastic = new G4Material("plasti <<  98  new G4Material("carbon", 1.0*g/cm3, 1);
 97   plastic->AddElement(H, 1);                   <<  99  carbon->AddElement(C, 1);
 98   plastic->AddElement(C, 1);                   << 100  
 99                                                << 101  G4Material* plastic = 
100   // or use G4-NIST materials data base        << 102  new G4Material("plastic", 1.0*g/cm3, 2);
101   //                                           << 103  plastic->AddElement(H, 1);
102   G4NistManager* man = G4NistManager::Instance << 104  plastic->AddElement(C, 1);
103   man->FindOrBuildMaterial("G4_B");            << 105  
104                                                << 106  // or use G4-NIST materials data base
105   G4Element* O = man->FindOrBuildElement("O"); << 107  //
106   G4Element* Hf = man->FindOrBuildElement("Hf" << 108  G4NistManager* man = G4NistManager::Instance();
107                                                << 109  man->FindOrBuildMaterial("G4_B");
108   G4Material* HfO2 = new G4Material("HfO2", 9. << 110 
109   HfO2->AddElement(Hf, 1);                     << 111  G4Element*  O = man->FindOrBuildElement("O");
110   HfO2->AddElement(O, 2);                      << 112  G4Element* Hf = man->FindOrBuildElement("Hf");
111                                                << 113  
112   /// G4cout << *(G4Material::GetMaterialTable << 114  G4Material* HfO2 = new G4Material("HfO2", 9.68*g/cm3, 2);
                                                   >> 115  HfO2->AddElement(Hf, 1);
                                                   >> 116  HfO2->AddElement(O , 2);
                                                   >> 117  
                                                   >> 118  ///G4cout << *(G4Material::GetMaterialTable()) << G4endl;
113 }                                                 119 }
114                                                   120 
115 //....oooOO0OOooo........oooOO0OOooo........oo    121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
116                                                   122 
117 G4Material* DetectorConstruction::MaterialWith << 123 G4Material* DetectorConstruction::MaterialWithSingleIsotope( G4String name,
118                                                << 124                            G4String symbol, G4double density, G4int Z, G4int A)
119 {                                              << 125 {
120   // define a material from an isotope         << 126  // define a material from an isotope
121   //                                           << 127  //
122   G4int ncomponents;                           << 128  G4int ncomponents;
123   G4double abundance, massfraction;            << 129  G4double abundance, massfraction;
124                                                << 130 
125   G4Isotope* isotope = new G4Isotope(symbol, Z << 131  G4Isotope* isotope = new G4Isotope(symbol, Z, A);
126                                                << 132  
127   G4Element* element = new G4Element(name, sym << 133  G4Element* element  = new G4Element(name, symbol, ncomponents=1);
128   element->AddIsotope(isotope, abundance = 100 << 134  element->AddIsotope(isotope, abundance= 100.*perCent);
129                                                << 135  
130   G4Material* material = new G4Material(name,  << 136  G4Material* material = new G4Material(name, density, ncomponents=1);
131   material->AddElement(element, massfraction = << 137  material->AddElement(element, massfraction=100.*perCent);
132                                                   138 
133   return material;                             << 139  return material;
134 }                                                 140 }
135                                                   141 
136 //....oooOO0OOooo........oooOO0OOooo........oo    142 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
137                                                   143 
138 G4VPhysicalVolume* DetectorConstruction::Const    144 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
139 {                                                 145 {
140   // Cleanup old geometry                         146   // Cleanup old geometry
141   G4GeometryManager::GetInstance()->OpenGeomet    147   G4GeometryManager::GetInstance()->OpenGeometry();
142   G4PhysicalVolumeStore::GetInstance()->Clean(    148   G4PhysicalVolumeStore::GetInstance()->Clean();
143   G4LogicalVolumeStore::GetInstance()->Clean()    149   G4LogicalVolumeStore::GetInstance()->Clean();
144   G4SolidStore::GetInstance()->Clean();           150   G4SolidStore::GetInstance()->Clean();
145                                                   151 
146   G4Box* sBox = new G4Box("Container",  // its << 152   G4Box*
147                           fBoxSize / 2, fBoxSi << 153   sBox = new G4Box("Container",                         //its name
148                                                << 154                    fBoxSize/2,fBoxSize/2,fBoxSize/2);   //its dimensions
149   fLBox = new G4LogicalVolume(sBox,  // its sh << 155 
150                               fMaterial,  // i << 156   fLBox = new G4LogicalVolume(sBox,                     //its shape
151                               fMaterial->GetNa << 157                              fMaterial,                 //its material
152                                                << 158                              fMaterial->GetName());     //its name
153   fPBox = new G4PVPlacement(0,  // no rotation << 159 
154                             G4ThreeVector(),   << 160   fPBox = new G4PVPlacement(0,                          //no rotation
155                             fLBox,  // its log << 161                             G4ThreeVector(),            //at (0,0,0)
156                             fMaterial->GetName << 162                             fLBox,                      //its logical volume
157                             0,  // its mother  << 163                             fMaterial->GetName(),       //its name
158                             false,  // no bool << 164                             0,                          //its mother  volume
159                             0);  // copy numbe << 165                             false,                      //no boolean operation
160                                                << 166                             0);                         //copy number
                                                   >> 167                            
161   PrintParameters();                              168   PrintParameters();
162                                                << 169   
163   // always return the root volume             << 170   //always return the root volume
164   //                                              171   //
165   return fPBox;                                   172   return fPBox;
166 }                                                 173 }
167                                                   174 
168 //....oooOO0OOooo........oooOO0OOooo........oo    175 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
169                                                   176 
170 void DetectorConstruction::PrintParameters()      177 void DetectorConstruction::PrintParameters()
171 {                                                 178 {
172   G4cout << "\n The Box is " << G4BestUnit(fBo << 179   G4cout << "\n The Box is " << G4BestUnit(fBoxSize,"Length")
173          << "\n \n"                            << 180          << " of " << fMaterial->GetName() 
174          << fMaterial << G4endl;               << 181          << "\n \n" << fMaterial << G4endl;
175 }                                                 182 }
176                                                   183 
177 //....oooOO0OOooo........oooOO0OOooo........oo    184 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
178                                                   185 
179 void DetectorConstruction::SetMaterial(G4Strin    186 void DetectorConstruction::SetMaterial(G4String materialChoice)
180 {                                                 187 {
181   // search the material by its name              188   // search the material by its name
182   G4Material* pttoMaterial = G4NistManager::In << 189   G4Material* pttoMaterial = 
183                                                << 190      G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
184   if (pttoMaterial) {                          << 191   
185     if (fMaterial != pttoMaterial) {           << 192   if (pttoMaterial) { 
                                                   >> 193     if(fMaterial != pttoMaterial) {
186       fMaterial = pttoMaterial;                   194       fMaterial = pttoMaterial;
187       if (fLBox) {                             << 195       if(fLBox) { fLBox->SetMaterial(pttoMaterial); }
188         fLBox->SetMaterial(pttoMaterial);      << 
189       }                                        << 
190       G4RunManager::GetRunManager()->PhysicsHa    196       G4RunManager::GetRunManager()->PhysicsHasBeenModified();
191     }                                             197     }
192   }                                            << 198   } else {
193   else {                                       << 199     G4cout << "\n--> warning from DetectorConstruction::SetMaterial : "
194     G4cout << "\n--> warning from DetectorCons << 200            << materialChoice << " not found" << G4endl;
195            << " not found" << G4endl;          << 201   }              
196   }                                            << 
197 }                                                 202 }
198                                                   203 
199 //....oooOO0OOooo........oooOO0OOooo........oo    204 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
200                                                   205 
201 void DetectorConstruction::SetSize(G4double va    206 void DetectorConstruction::SetSize(G4double value)
202 {                                                 207 {
203   fBoxSize = value;                               208   fBoxSize = value;
204   G4RunManager::GetRunManager()->ReinitializeG    209   G4RunManager::GetRunManager()->ReinitializeGeometry();
205 }                                                 210 }
206                                                   211 
207 //....oooOO0OOooo........oooOO0OOooo........oo    212 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
208                                                   213