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 ExGflash2DetectorConstruction.cc 28 /// \brief Implementation of the ExGflash2Dete 29 // 30 // Created by Joanna Weng 26.11.2004 31 32 //....oooOO0OOooo........oooOO0OOooo........oo 33 //....oooOO0OOooo........oooOO0OOooo........oo 34 35 // User Classes 36 #include "ExGflash2DetectorConstruction.hh" 37 38 #include "ExGflash2SensitiveDetector.hh" 39 40 // G4 Classes 41 #include "G4AutoDelete.hh" 42 #include "G4Box.hh" 43 #include "G4Colour.hh" 44 #include "G4LogicalVolume.hh" 45 #include "G4Material.hh" 46 #include "G4NistManager.hh" 47 #include "G4PVPlacement.hh" 48 #include "G4SDManager.hh" 49 #include "G4SystemOfUnits.hh" 50 #include "G4ThreeVector.hh" 51 #include "G4VPhysicalVolume.hh" 52 #include "G4VisAttributes.hh" 53 #include "globals.hh" 54 55 //....oooOO0OOooo........oooOO0OOooo........oo 56 57 ExGflash2DetectorConstruction::ExGflash2Detect 58 { 59 G4cout << "ExGflash2DetectorConstruction::De 60 } 61 62 //....oooOO0OOooo........oooOO0OOooo........oo 63 64 ExGflash2DetectorConstruction::~ExGflash2Detec 65 66 //....oooOO0OOooo........oooOO0OOooo........oo 67 68 G4VPhysicalVolume* ExGflash2DetectorConstructi 69 { 70 //--------- Definitions of Solids, Logical V 71 G4cout << "Defining the materials" << G4endl 72 // Get nist material manager 73 G4NistManager* nistManager = G4NistManager:: 74 // Build materials 75 G4Material* air = nistManager->FindOrBuildMa 76 G4Material* pbWO4 = nistManager->FindOrBuild 77 78 /******************************* 79 * The Experimental Hall * 80 *******************************/ 81 G4double experimentalHall_x = 1000. * cm; 82 G4double experimentalHall_y = 1000. * cm; 83 G4double experimentalHall_z = 1000. * cm; 84 85 G4VSolid* experimentalHall_box = new G4Box(" 86 e 87 e 88 e 89 90 auto experimentalHallLog = new G4LogicalVolu 91 92 93 94 G4VPhysicalVolume* experimentalHallPhys = 95 new G4PVPlacement(nullptr, 96 G4ThreeVector(), // at 97 "expHall", experimentalH 98 99 //------------------------------ 100 // Calorimeter segments 101 //------------------------------ 102 // Simplified `CMS-like` PbWO4 crystal calor 103 104 G4int nbOfCrystals = 10; // this are the cr 105 // cube of 10 x 10 106 // don't change it 107 // the readout in 108 // dimensions and 109 // in this version 110 // Simplified `CMS-like` PbWO4 crystal calor 111 G4double calo_xside = 31 * cm; 112 G4double calo_yside = 31 * cm; 113 G4double calo_zside = 24 * cm; 114 115 G4double crystalWidth = 3 * cm; 116 G4double crystalLength = 24 * cm; 117 118 calo_xside = (crystalWidth * nbOfCrystals) + 119 calo_yside = (crystalWidth * nbOfCrystals) + 120 calo_zside = crystalLength; 121 122 auto calo_box = new G4Box("CMS calorimeter", 123 calo_xside / 2., 124 calo_yside / 2., c 125 auto caloLog = new G4LogicalVolume(calo_box, 126 air, // 127 "calo log 128 nullptr, 129 nullptr, 130 nullptr); 131 132 G4double xpos = 0.0; 133 G4double ypos = 0.0; 134 G4double zpos = 100.0 * cm; 135 new G4PVPlacement(nullptr, G4ThreeVector(xpo 136 experimentalHallLog, false 137 138 // Crystals 139 G4VSolid* crystal_box = new G4Box("Crystal", 140 crystalWid 141 // size 142 fCrystalLog = new G4LogicalVolume(crystal_bo 143 pbWO4, // 144 "CrystalLo 145 146 for (G4int i = 0; i < nbOfCrystals; i++) { 147 for (G4int j = 0; j < nbOfCrystals; j++) { 148 G4int n = i * 10 + j; 149 G4ThreeVector crystalPos((i * crystalWid 150 fCrystalPhys[n] = new G4PVPlacement(null 151 crys 152 fCry 153 "cry 154 calo 155 } 156 } 157 G4cout << "There are " << nbOfCrystals << " 158 << nbOfCrystals * nbOfCrystals << " c 159 G4cout << "They have width of " << crystalW 160 << crystalLength / cm << " cm. The Ma 161 162 experimentalHallLog->SetVisAttributes(G4VisA 163 auto caloVisAtt = new G4VisAttributes(G4Colo 164 auto crystalVisAtt = new G4VisAttributes(G4C 165 caloLog->SetVisAttributes(caloVisAtt); 166 fCrystalLog->SetVisAttributes(crystalVisAtt) 167 168 // Parametrisation region defined in ExGflas 169 170 return experimentalHallPhys; 171 } 172 173 //....oooOO0OOooo........oooOO0OOooo........oo 174 175 void ExGflash2DetectorConstruction::ConstructS 176 { 177 // -- sensitive detectors: 178 G4SDManager* SDman = G4SDManager::GetSDMpoin 179 auto CaloSD = new ExGflash2SensitiveDetector 180 SDman->AddNewDetector(CaloSD); 181 fCrystalLog->SetSensitiveDetector(CaloSD); 182 183 // Fast simulation implemented in ExGflash2P 184 } 185 186 //....oooOO0OOooo........oooOO0OOooo........oo 187