Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 /// \file ExGflash3ParallelWorld.cc 28 /// \brief Implementation of the ExGflash3Para 29 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 33 // User Classes 34 #include "ExGflash3ParallelWorld.hh" 35 36 #include "ExGflash3SensitiveDetector.hh" 37 38 // G4 Classes 39 #include "G4AutoDelete.hh" 40 #include "G4Box.hh" 41 #include "G4Colour.hh" 42 #include "G4LogicalVolume.hh" 43 #include "G4Material.hh" 44 #include "G4NistManager.hh" 45 #include "G4PVPlacement.hh" 46 #include "G4SDManager.hh" 47 #include "G4SystemOfUnits.hh" 48 #include "G4ThreeVector.hh" 49 #include "G4VPhysicalVolume.hh" 50 #include "G4VisAttributes.hh" 51 #include "globals.hh" 52 53 //....oooOO0OOooo........oooOO0OOooo........oo 54 55 ExGflash3ParallelWorld::ExGflash3ParallelWorld 56 : G4VUserParallelWorld(aWorldName) 57 { 58 G4cout << "ExGflash3ParallelWorld::Parralel 59 } 60 61 //....oooOO0OOooo........oooOO0OOooo........oo 62 63 ExGflash3ParallelWorld::~ExGflash3ParallelWorl 64 65 //....oooOO0OOooo........oooOO0OOooo........oo 66 67 void ExGflash3ParallelWorld::Construct() 68 { 69 // In parallel world material does not matte 70 G4Material* dummy = nullptr; 71 72 // Build parallel/ghost geometry: 73 auto ghostLogicalVolume = GetWorld()->GetLog 74 75 // Use part of the Ex1GflashDetectorConstruc 76 77 G4int nbOfCrystals = 10; // this are the cr 78 // cube of 10 x 10 79 // don't change it 80 // the readout in 81 // dimensions and 82 // in this version 83 // Simplified `CMS-like` PbWO4 crystal calor 84 G4double calo_xside = 31 * cm; 85 G4double calo_yside = 31 * cm; 86 G4double calo_zside = 24 * cm; 87 88 G4double crystalWidth = 3 * cm; 89 G4double crystalLength = 24 * cm; 90 91 calo_xside = (crystalWidth * nbOfCrystals) + 92 calo_yside = (crystalWidth * nbOfCrystals) + 93 calo_zside = crystalLength; 94 95 auto calo_box = new G4Box("CMS calorimeter", 96 calo_xside / 2., 97 calo_yside / 2., c 98 auto caloLog = new G4LogicalVolume(calo_box, 99 dummy, / 100 "calo log 101 nullptr, 102 nullptr, 103 nullptr); 104 105 G4double xpos = 0.0; 106 G4double ypos = 0.0; 107 G4double zpos = 100.0 * cm; 108 new G4PVPlacement(nullptr, G4ThreeVector(xpo 109 ghostLogicalVolume, false, 110 111 // Crystals 112 G4VSolid* crystal_box = new G4Box("Crystal", 113 crystalWid 114 // size 115 fCrystalLog = new G4LogicalVolume(crystal_bo 116 dummy, // 117 "CrystalLo 118 119 for (G4int i = 0; i < nbOfCrystals; i++) { 120 for (G4int j = 0; j < nbOfCrystals; j++) { 121 G4int n = i * 10 + j; 122 G4ThreeVector crystalPos((i * crystalWid 123 fCrystalPhys[n] = new G4PVPlacement(null 124 crys 125 fCry 126 "cry 127 calo 128 } 129 } 130 G4cout << "There are " << nbOfCrystals << " 131 << nbOfCrystals * nbOfCrystals << " c 132 G4cout << "They have width of " << crystalW 133 << crystalLength / cm << " cm. " << G 134 135 auto caloVisAtt = new G4VisAttributes(G4Colo 136 auto crystalVisAtt = new G4VisAttributes(G4C 137 caloLog->SetVisAttributes(caloVisAtt); 138 fCrystalLog->SetVisAttributes(crystalVisAtt) 139 } 140 141 //....oooOO0OOooo........oooOO0OOooo........oo 142 143 void ExGflash3ParallelWorld::ConstructSD() 144 { 145 // -- sensitive detectors: 146 G4SDManager* SDman = G4SDManager::GetSDMpoin 147 auto CaloSD = new ExGflash3SensitiveDetector 148 SDman->AddNewDetector(CaloSD); 149 fCrystalLog->SetSensitiveDetector(CaloSD); 150 } 151 152 //....oooOO0OOooo........oooOO0OOooo........oo 153