Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 /// \file B1/src/DetectorConstruction.cc 28 /// \brief Implementation of the B1::DetectorConstruction class 29 30 #include "DetectorConstruction.hh" 31 32 #include "G4Box.hh" 33 #include "G4LogicalVolume.hh" 34 #include "G4NistManager.hh" 35 #include "G4PVPlacement.hh" 36 #include "G4RunManager.hh" 37 #include "G4SubtractionSolid.hh" 38 #include "G4SystemOfUnits.hh" 39 #include "G4Tubs.hh" 40 41 #include <cmath> 42 43 namespace VtkVis 44 { 45 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 48 G4VPhysicalVolume* DetectorConstruction::Construct() 49 { 50 // Get nist material manager 51 G4NistManager* nist = G4NistManager::Instance(); 52 53 // Envelope parameters 54 // 55 G4double env_sizeXY = 10 * m, env_sizeZ = 10 * m; 56 G4Material* env_mat = nist->FindOrBuildMaterial("G4_Galactic"); 57 58 // Option to switch on/off checking of volumes overlaps 59 // 60 G4bool checkOverlaps = true; 61 62 // 63 // World 64 // 65 G4double world_sizeXY = 1.2 * env_sizeXY; 66 G4double world_sizeZ = 1.2 * env_sizeZ; 67 G4Material* world_mat = nist->FindOrBuildMaterial("G4_AIR"); 68 69 auto solidWorld = 70 new G4Box("World", // its name 71 0.5 * world_sizeXY, 0.5 * world_sizeXY, 0.5 * world_sizeZ); // its size 72 73 auto logicWorld = new G4LogicalVolume(solidWorld, // its solid 74 world_mat, // its material 75 "World"); // its name 76 77 auto physWorld = new G4PVPlacement(nullptr, // no rotation 78 G4ThreeVector(), // at (0,0,0) 79 logicWorld, // its logical volume 80 "World", // its name 81 nullptr, // its mother volume 82 false, // no boolean operation 83 0, // copy number 84 checkOverlaps); // overlaps checking 85 86 // 87 // Envelope 88 // 89 auto solidEnv = new G4Box("Envelope", // its name 90 0.5 * env_sizeXY, 0.5 * env_sizeXY, 0.5 * env_sizeZ); // its size 91 92 auto logicEnv = new G4LogicalVolume(solidEnv, // its solid 93 env_mat, // its material 94 "Envelope"); // its name 95 96 new G4PVPlacement(nullptr, // no rotation 97 G4ThreeVector(), // at (0,0,0) 98 logicEnv, // its logical volume 99 "Envelope", // its name 100 logicWorld, // its mother volume 101 false, // no boolean operation 102 0, // copy number 103 checkOverlaps); // overlaps checking 104 105 auto dumpPos = G4ThreeVector(0, 0, 0); 106 auto dumpMat = nist->FindOrBuildMaterial("G4_CONCRETE"); 107 auto dumpShape1 = new G4Box("dumpShape1", 1 * m, 1 * m, 0.5 * m); 108 109 auto dumpCutoutPos = G4ThreeVector(0, 0, -5.1 * cm); 110 auto dumpCutoutShape = new G4Tubs("DumpCutoutShape", 0 * cm, 10 * cm, 90 / 2 * cm, 0, 2 * M_PI); 111 112 auto dumpShape2 = 113 new G4SubtractionSolid("dumpShape2", dumpShape1, dumpCutoutShape, 0, dumpCutoutPos); 114 115 auto dumpLogical = new G4LogicalVolume(dumpShape2, 116 dumpMat, // its material 117 "dump"); 118 // its name 119 new G4PVPlacement(nullptr, // no rotation 120 dumpPos, // at position 121 dumpLogical, // its logical volume 122 "dump", // its name 123 logicEnv, // its mother volume 124 false, // no boolean operation 125 0, // copy number 126 checkOverlaps); // overlaps checking 127 128 fScoringVolume = dumpLogical; 129 130 // 131 // always return the physical World 132 // 133 return physWorld; 134 } 135 136 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 137 138 } // namespace VtkVis 139