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 ExGflash1DetectorConstruction.cc 28 /// \brief Implementation of the ExGflash1Dete 29 // 30 // Created by Joanna Weng 26.11.2004 31 32 //....oooOO0OOooo........oooOO0OOooo........oo 33 //....oooOO0OOooo........oooOO0OOooo........oo 34 35 // User Classes 36 #include "ExGflash1DetectorConstruction.hh" 37 38 #include "ExGflash1SensitiveDetector.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 // fast simulation 56 #include "GFlashHitMaker.hh" 57 #include "GFlashHomoShowerParameterisation.hh" 58 #include "GFlashParticleBounds.hh" 59 #include "GFlashShowerModel.hh" 60 61 #include "G4FastSimulationManager.hh" 62 63 //....oooOO0OOooo........oooOO0OOooo........oo 64 65 ExGflash1DetectorConstruction::ExGflash1Detect 66 { 67 G4cout << "ExGflash1DetectorConstruction::De 68 } 69 70 //....oooOO0OOooo........oooOO0OOooo........oo 71 72 ExGflash1DetectorConstruction::~ExGflash1Detec 73 { 74 delete fFastShowerModel; 75 delete fParameterisation; 76 delete fParticleBounds; 77 delete fHitMaker; 78 } 79 80 //....oooOO0OOooo........oooOO0OOooo........oo 81 82 G4VPhysicalVolume* ExGflash1DetectorConstructi 83 { 84 //--------- Definitions of Solids, Logical V 85 G4cout << "Defining the materials" << G4endl 86 // Get nist material manager 87 G4NistManager* nistManager = G4NistManager:: 88 // Build materials 89 G4Material* air = nistManager->FindOrBuildMa 90 G4Material* pbWO4 = nistManager->FindOrBuild 91 92 /******************************* 93 * The Experimental Hall * 94 *******************************/ 95 G4double experimentalHall_x = 1000. * cm; 96 G4double experimentalHall_y = 1000. * cm; 97 G4double experimentalHall_z = 1000. * cm; 98 99 G4VSolid* experimentalHall_box = new G4Box(" 100 e 101 e 102 e 103 104 auto experimentalHallLog = new G4LogicalVolu 105 106 107 108 G4VPhysicalVolume* experimentalHallPhys = 109 new G4PVPlacement(nullptr, 110 G4ThreeVector(), // at 111 "expHall", experimentalH 112 113 //------------------------------ 114 // Calorimeter segments 115 //------------------------------ 116 // Simplified `CMS-like` PbWO4 crystal calor 117 118 G4int nbOfCrystals = 10; // this are the cr 119 // cube of 10 x 10 120 // don't change it 121 // the readout in 122 // dimensions and 123 // in this version 124 // Simplified `CMS-like` PbWO4 crystal calor 125 G4double calo_xside = 31 * cm; 126 G4double calo_yside = 31 * cm; 127 G4double calo_zside = 24 * cm; 128 129 G4double crystalWidth = 3 * cm; 130 G4double crystalLength = 24 * cm; 131 132 calo_xside = (crystalWidth * nbOfCrystals) + 133 calo_yside = (crystalWidth * nbOfCrystals) + 134 calo_zside = crystalLength; 135 136 auto calo_box = new G4Box("CMS calorimeter", 137 calo_xside / 2., 138 calo_yside / 2., c 139 auto caloLog = new G4LogicalVolume(calo_box, 140 air, // 141 "calo log 142 nullptr, 143 nullptr, 144 nullptr); 145 146 G4double xpos = 0.0; 147 G4double ypos = 0.0; 148 G4double zpos = 100.0 * cm; 149 new G4PVPlacement(nullptr, G4ThreeVector(xpo 150 experimentalHallLog, false 151 152 // Crystals 153 G4VSolid* crystal_box = new G4Box("Crystal", 154 crystalWid 155 // size 156 fCrystalLog = new G4LogicalVolume(crystal_bo 157 pbWO4, // 158 "CrystalLo 159 160 for (G4int i = 0; i < nbOfCrystals; i++) { 161 for (G4int j = 0; j < nbOfCrystals; j++) { 162 G4int n = i * 10 + j; 163 G4ThreeVector crystalPos((i * crystalWid 164 fCrystalPhys[n] = new G4PVPlacement(null 165 crys 166 fCry 167 "cry 168 calo 169 } 170 } 171 G4cout << "There are " << nbOfCrystals << " 172 << nbOfCrystals * nbOfCrystals << " c 173 G4cout << "They have width of " << crystalW 174 << crystalLength / cm << " cm. The Ma 175 176 experimentalHallLog->SetVisAttributes(G4VisA 177 auto caloVisAtt = new G4VisAttributes(G4Colo 178 auto crystalVisAtt = new G4VisAttributes(G4C 179 caloLog->SetVisAttributes(caloVisAtt); 180 fCrystalLog->SetVisAttributes(crystalVisAtt) 181 182 // define the fParameterisation region 183 fRegion = new G4Region("crystals"); 184 caloLog->SetRegion(fRegion); 185 fRegion->AddRootLogicalVolume(caloLog); 186 187 return experimentalHallPhys; 188 } 189 190 //....oooOO0OOooo........oooOO0OOooo........oo 191 192 void ExGflash1DetectorConstruction::ConstructS 193 { 194 // -- sensitive detectors: 195 G4SDManager* SDman = G4SDManager::GetSDMpoin 196 auto CaloSD = new ExGflash1SensitiveDetector 197 SDman->AddNewDetector(CaloSD); 198 fCrystalLog->SetSensitiveDetector(CaloSD); 199 200 // Get nist material manager 201 G4NistManager* nistManager = G4NistManager:: 202 G4Material* pbWO4 = nistManager->FindOrBuild 203 // -- fast simulation models: 204 // ***************************************** 205 // * Initializing shower modell 206 // ***************************************** 207 G4cout << "Creating shower parameterization 208 fFastShowerModel = new GFlashShowerModel("fF 209 fParameterisation = new GFlashHomoShowerPara 210 fFastShowerModel->SetParameterisation(*fPara 211 // Energy Cuts to kill particles: 212 fParticleBounds = new GFlashParticleBounds() 213 fFastShowerModel->SetParticleBounds(*fPartic 214 // Makes the EnergieSpots 215 fHitMaker = new GFlashHitMaker(); 216 fFastShowerModel->SetHitMaker(*fHitMaker); 217 G4cout << "end shower parameterization." << 218 // ***************************************** 219 } 220 221 //....oooOO0OOooo........oooOO0OOooo........oo 222