Geant4 Cross Reference |
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.hh 26 /// \file DetectorConstruction.hh 27 /// \brief Definition of the DetectorConstruct 27 /// \brief Definition of the DetectorConstruction class 28 28 29 //....oooOO0OOooo........oooOO0OOooo........oo 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 30 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 31 32 #ifndef DetectorConstruction_h 32 #ifndef DetectorConstruction_h 33 #define DetectorConstruction_h 1 33 #define DetectorConstruction_h 1 34 34 35 #include "CLHEP/Units/SystemOfUnits.h" << 36 << 37 #include "G4VUserDetectorConstruction.hh" 35 #include "G4VUserDetectorConstruction.hh" 38 #include "globals.hh" 36 #include "globals.hh" 39 37 >> 38 #include "CLHEP/Units/SystemOfUnits.h" >> 39 40 class G4Box; 40 class G4Box; 41 class G4LogicalVolume; 41 class G4LogicalVolume; 42 class G4VPhysicalVolume; 42 class G4VPhysicalVolume; 43 class G4Material; 43 class G4Material; 44 class DetectorMessenger; 44 class DetectorMessenger; 45 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 47 48 class DetectorConstruction : public G4VUserDet 48 class DetectorConstruction : public G4VUserDetectorConstruction 49 { 49 { 50 public: 50 public: 51 DetectorConstruction(); << 51 52 ~DetectorConstruction() override; << 52 DetectorConstruction(); >> 53 ~DetectorConstruction() override; 53 54 54 public: 55 public: 55 G4VPhysicalVolume* Construct() override; 56 G4VPhysicalVolume* Construct() override; 56 57 57 void SetAbsorberMaterial(const G4String&); << 58 void SetAbsorberMaterial (const G4String&); 58 void SetAbsorberThickness(G4double); 59 void SetAbsorberThickness(G4double); 59 60 60 void SetGapMaterial(const G4String&); << 61 void SetGapMaterial (const G4String&); 61 void SetGapThickness(G4double); 62 void SetGapThickness(G4double); 62 63 63 void SetCalorSizeYZ(G4double); 64 void SetCalorSizeYZ(G4double); 64 void SetNbOfLayers(G4int); << 65 void SetNbOfLayers (G4int); 65 66 66 void PrintCalorParameters(); 67 void PrintCalorParameters(); 67 68 68 G4double GetWorldSizeX() { return fWorldSi << 69 G4double GetWorldSizeX() { return fWorldSizeX; } 69 G4double GetWorldSizeYZ() { return fWorldS 70 G4double GetWorldSizeYZ() { return fWorldSizeYZ; } 70 71 71 G4double GetCalorThickness() { return fCal 72 G4double GetCalorThickness() { return fCalorThickness; } 72 G4double GetCalorSizeYZ() { return fCalorS << 73 G4double GetCalorSizeYZ() { return fCalorSizeYZ; } 73 74 74 G4int GetNbOfLayers() { return fNbOfLayers 75 G4int GetNbOfLayers() { return fNbOfLayers; } 75 76 76 G4Material* GetAbsorberMaterial() { return << 77 G4Material* GetAbsorberMaterial() { return fAbsorberMaterial; } 77 G4double GetAbsorberThickness() { return f << 78 G4double GetAbsorberThickness() { return fAbsorberThickness; } 78 79 79 G4Material* GetGapMaterial() { return fGap << 80 G4Material* GetGapMaterial() { return fGapMaterial; } 80 G4double GetGapThickness() { return fGapTh << 81 G4double GetGapThickness() { return fGapThickness; } 81 82 82 const G4VPhysicalVolume* GetphysiWorld() { 83 const G4VPhysicalVolume* GetphysiWorld() { return fPhysiWorld; } 83 const G4VPhysicalVolume* GetAbsorber() { r << 84 const G4VPhysicalVolume* GetAbsorber() { return fPhysiAbsorber; } 84 const G4VPhysicalVolume* GetGap() { return << 85 const G4VPhysicalVolume* GetGap() { return fPhysiGap; } 85 86 86 private: 87 private: 87 void DefineMaterials(); 88 void DefineMaterials(); 88 void ComputeCalorParameters(); 89 void ComputeCalorParameters(); 89 G4VPhysicalVolume* ConstructCalorimeter(); 90 G4VPhysicalVolume* ConstructCalorimeter(); 90 91 91 G4Material* fAbsorberMaterial = nullptr; 92 G4Material* fAbsorberMaterial = nullptr; 92 G4Material* fGapMaterial = nullptr; 93 G4Material* fGapMaterial = nullptr; 93 G4Material* fDefaultMaterial = nullptr; 94 G4Material* fDefaultMaterial = nullptr; 94 95 95 G4int fNbOfLayers = 10; << 96 G4int fNbOfLayers = 10; 96 G4double fAbsorberThickness = 10. * CLHEP: << 97 G4double fAbsorberThickness = 10.*CLHEP::mm; 97 G4double fGapThickness = 5. * CLHEP::mm; << 98 G4double fGapThickness = 5.*CLHEP::mm; 98 G4double fCalorSizeYZ = 10. * CLHEP::cm; << 99 G4double fCalorSizeYZ = 10.*CLHEP::cm; 99 << 100 100 G4double fCalorThickness = 0.; // will be << 101 G4double fCalorThickness = 0.; // will be computed 101 G4double fLayerThickness = 0.; // will be << 102 G4double fLayerThickness= 0.; // will be computed 102 G4double fWorldSizeYZ = 0.; // will be co << 103 G4double fWorldSizeYZ = 0.; // will be computed; 103 G4double fWorldSizeX = 0.; // will be com << 104 G4double fWorldSizeX = 0.; // will be computed; 104 << 105 105 G4Box* fSolidWorld = nullptr; // pointer << 106 G4Box* fSolidWorld = nullptr; //pointer to the solid World 106 G4LogicalVolume* fLogicWorld = nullptr; / << 107 G4LogicalVolume* fLogicWorld = nullptr; //pointer to the logical World 107 G4VPhysicalVolume* fPhysiWorld = nullptr; << 108 G4VPhysicalVolume* fPhysiWorld = nullptr; //pointer to the physical World 108 << 109 109 G4Box* fSolidCalor = nullptr; // pointer << 110 G4Box* fSolidCalor = nullptr; //pointer to the solid Calor 110 G4LogicalVolume* fLogicCalor = nullptr; / << 111 G4LogicalVolume* fLogicCalor = nullptr; //pointer to the logical Calor 111 G4VPhysicalVolume* fPhysiCalor = nullptr; << 112 G4VPhysicalVolume* fPhysiCalor = nullptr; //pointer to the physical Calor 112 << 113 113 G4Box* fSolidLayer = nullptr; // pointer << 114 G4Box* fSolidLayer = nullptr; //pointer to the solid Layer 114 G4LogicalVolume* fLogicLayer = nullptr; / << 115 G4LogicalVolume* fLogicLayer = nullptr; //pointer to the logical Layer 115 G4VPhysicalVolume* fPhysiLayer = nullptr; << 116 G4VPhysicalVolume* fPhysiLayer = nullptr; //pointer to the physical Layer 116 << 117 117 G4Box* fSolidAbsorber = nullptr; // point << 118 G4Box* fSolidAbsorber = nullptr; //pointer to the solid Absorber 118 G4LogicalVolume* fLogicAbsorber = nullptr; << 119 G4LogicalVolume* fLogicAbsorber = nullptr; //pointer to the logical Absorber 119 G4VPhysicalVolume* fPhysiAbsorber = nullpt << 120 G4VPhysicalVolume* fPhysiAbsorber = nullptr; //pointer to the physical Absorber 120 << 121 121 G4Box* fSolidGap = nullptr; // pointer to << 122 G4Box* fSolidGap = nullptr; //pointer to the solid Gap 122 G4LogicalVolume* fLogicGap = nullptr; // << 123 G4LogicalVolume* fLogicGap = nullptr; //pointer to the logical Gap 123 G4VPhysicalVolume* fPhysiGap = nullptr; / << 124 G4VPhysicalVolume* fPhysiGap = nullptr; //pointer to the physical Gap 124 125 125 DetectorMessenger* fDetectorMessenger = nu << 126 DetectorMessenger* fDetectorMessenger = nullptr; //pointer to the Messenger 126 }; 127 }; 127 128 128 //....oooOO0OOooo........oooOO0OOooo........oo 129 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 129 130 130 inline void DetectorConstruction::ComputeCalor 131 inline void DetectorConstruction::ComputeCalorParameters() 131 { 132 { 132 // Compute derived parameters of the calorim 133 // Compute derived parameters of the calorimeter 133 fLayerThickness = fAbsorberThickness + fGapT << 134 fLayerThickness = fAbsorberThickness + fGapThickness; 134 fCalorThickness = fNbOfLayers * fLayerThickn << 135 fCalorThickness = fNbOfLayers*fLayerThickness; 135 136 136 fWorldSizeX = 1.2 * fCalorThickness; << 137 fWorldSizeX = 1.2*fCalorThickness; fWorldSizeYZ = 1.2*fCalorSizeYZ; 137 fWorldSizeYZ = 1.2 * fCalorSizeYZ; << 138 } 138 } 139 139 140 //....oooOO0OOooo........oooOO0OOooo........oo 140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 141 141 142 #endif 142 #endif >> 143 143 144