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 field/field02/include/F02DetectorCon 26 /// \file field/field02/include/F02DetectorConstruction.hh 27 /// \brief Definition of the F02DetectorConstr 27 /// \brief Definition of the F02DetectorConstruction class 28 // 28 // 29 // 29 // 30 // 30 // 31 // 31 // 32 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 33 //....oooOO0OOooo........oooOO0OOooo........oo 33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 34 34 35 #ifndef F02DetectorConstruction_h 35 #ifndef F02DetectorConstruction_h 36 #define F02DetectorConstruction_h 1 36 #define F02DetectorConstruction_h 1 37 37 38 #include "G4VUserDetectorConstruction.hh" 38 #include "G4VUserDetectorConstruction.hh" 39 #include "G4Cache.hh" 39 #include "G4Cache.hh" 40 #include "G4ThreeVector.hh" << 41 << 42 #include "CLHEP/Units/SystemOfUnits.h" << 43 40 44 class G4Box; 41 class G4Box; 45 class G4Tubs; 42 class G4Tubs; 46 class G4LogicalVolume; 43 class G4LogicalVolume; 47 class G4VPhysicalVolume; 44 class G4VPhysicalVolume; 48 45 49 class G4Material; 46 class G4Material; 50 class G4UniformMagField; 47 class G4UniformMagField; 51 48 52 class F02DetectorMessenger; 49 class F02DetectorMessenger; 53 class F02CalorimeterSD; 50 class F02CalorimeterSD; >> 51 class F02ElectricFieldSetup; 54 52 55 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 56 54 57 class F02DetectorConstruction : public G4VUser 55 class F02DetectorConstruction : public G4VUserDetectorConstruction 58 { 56 { 59 public: 57 public: 60 F02DetectorConstruction(); << 61 ~F02DetectorConstruction() override; << 62 58 63 void SetAbsorberMaterial(G4String); << 59 F02DetectorConstruction(); 64 void SetAbsorberThickness(G4double); << 60 virtual ~F02DetectorConstruction(); 65 void SetAbsorberRadius(G4double); << 66 << 67 void SetAbsorberZpos(G4double); << 68 61 69 void SetWorldMaterial(G4String); << 62 public: 70 void SetWorldSizeZ(G4double); << 71 void SetWorldSizeR(G4double); << 72 63 73 void SetFieldValue(G4ThreeVector value); << 64 void SetAbsorberMaterial (G4String); >> 65 void SetAbsorberThickness(G4double); >> 66 void SetAbsorberRadius(G4double); >> 67 >> 68 void SetAbsorberZpos(G4double); >> 69 >> 70 void SetWorldMaterial(G4String); >> 71 void SetWorldSizeZ(G4double); >> 72 void SetWorldSizeR(G4double); 74 73 75 G4VPhysicalVolume* Construct() override; << 74 virtual G4VPhysicalVolume* Construct(); 76 void ConstructSDandField() override; << 75 virtual void ConstructSDandField(); 77 76 78 void PrintCalorParameters(); << 77 public: 79 78 80 G4Material* GetWorldMaterial() { return fW << 79 void PrintCalorParameters(); 81 G4double GetWorldSizeZ() { return fWorldSi << 82 G4double GetWorldSizeR() { return fWorldSi << 83 80 84 G4double GetAbsorberZpos() { return fZAbso << 81 G4Material* GetWorldMaterial() {return fWorldMaterial;} 85 G4double GetZStartAbs() { return fZStartAb << 82 G4double GetWorldSizeZ() {return fWorldSizeZ;} 86 G4double GetZEndAbs() { return fZEndAbs; } << 83 G4double GetWorldSizeR() {return fWorldSizeR;} >> 84 >> 85 G4double GetAbsorberZpos() {return fZAbsorber;} >> 86 G4double GetZStartAbs() {return fZStartAbs;} >> 87 G4double GetZEndAbs() {return fZEndAbs;} >> 88 >> 89 G4Material* GetAbsorberMaterial() {return fAbsorberMaterial;} >> 90 G4double GetAbsorberThickness(){return fAbsorberThickness;} >> 91 G4double GetAbsorberRadius() {return fAbsorberRadius;} >> 92 >> 93 const G4VPhysicalVolume* GetPhysiWorld() {return fPhysiWorld;} >> 94 const G4VPhysicalVolume* GetAbsorber() {return fPhysiAbsorber;} >> 95 G4LogicalVolume* GetLogicalAbsorber() {return fLogicAbsorber;} 87 96 88 G4Material* GetAbsorberMaterial() { return << 97 private: 89 G4double GetAbsorberThickness() { return f << 90 G4double GetAbsorberRadius() { return fAbs << 91 98 92 const G4VPhysicalVolume* GetPhysiWorld() { << 99 F02DetectorMessenger* fDetectorMessenger; // pointer -> Messenger 93 const G4VPhysicalVolume* GetAbsorber() { r << 100 G4Cache<F02CalorimeterSD*> fCalorimeterSD; // pointer -> sensitive detector 94 G4LogicalVolume* GetLogicalAbsorber() { re << 101 G4Cache<F02ElectricFieldSetup*> fEmFieldSetup; >> 102 >> 103 G4Tubs* fSolidWorld; // pointer to the solid World >> 104 G4LogicalVolume* fLogicWorld; // pointer to the logical World >> 105 G4VPhysicalVolume* fPhysiWorld; // pointer to the physical World >> 106 >> 107 G4Tubs* fSolidAbsorber; // pointer to the solid Absorber >> 108 G4LogicalVolume* fLogicAbsorber; // pointer to the logical Absorber >> 109 G4VPhysicalVolume* fPhysiAbsorber; // pointer to the physical Absorber >> 110 >> 111 G4Material* fAbsorberMaterial; >> 112 G4double fAbsorberThickness; >> 113 G4double fAbsorberRadius; >> 114 G4bool fWorldChanged; >> 115 >> 116 G4double fZAbsorber; >> 117 G4double fZStartAbs, fZEndAbs; >> 118 >> 119 G4Material* fWorldMaterial; >> 120 G4double fWorldSizeR; >> 121 G4double fWorldSizeZ; 95 122 96 private: 123 private: 97 void DefineMaterials(); << 124 98 void ComputeCalorParameters(); << 125 void DefineMaterials(); 99 G4VPhysicalVolume* ConstructCalorimeter(); << 126 void ComputeCalorParameters(); 100 << 127 G4VPhysicalVolume* ConstructCalorimeter(); 101 F02DetectorMessenger* fDetectorMessenger = << 102 G4Cache<F02CalorimeterSD*> fCalorimeterSD << 103 << 104 G4Tubs* fSolidWorld = nullptr; // pointer << 105 G4LogicalVolume* fLogicWorld = nullptr; / << 106 G4VPhysicalVolume* fPhysiWorld = nullptr; << 107 << 108 G4Tubs* fSolidAbsorber = nullptr; // poin << 109 G4LogicalVolume* fLogicAbsorber = nullptr; << 110 G4VPhysicalVolume* fPhysiAbsorber = nullpt << 111 << 112 G4Material* fAbsorberMaterial = nullptr; << 113 G4double fAbsorberThickness = 4. * CLHEP:: << 114 G4double fAbsorberRadius = 10. * CLHEP::cm << 115 G4bool fWorldChanged; << 116 << 117 G4double fZAbsorber = 36. * CLHEP::cm; << 118 G4double fZStartAbs = 0.; << 119 G4double fZEndAbs = 0.; << 120 << 121 G4Material* fWorldMaterial = nullptr; << 122 G4double fWorldSizeR = 20. * CLHEP::cm; << 123 G4double fWorldSizeZ = 80. * CLHEP::cm; << 124 << 125 G4ThreeVector fFieldVector = {0., 100000.0 << 126 }; 128 }; 127 129 128 //....oooOO0OOooo........oooOO0OOooo........oo 130 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 129 131 130 inline void F02DetectorConstruction::ComputeCa 132 inline void F02DetectorConstruction::ComputeCalorParameters() 131 { 133 { 132 // Compute derived parameters of the calorim << 134 // Compute derived parameters of the calorimeter 133 135 134 fZStartAbs = fZAbsorber - 0.5 * fAbsorberThi << 136 fZStartAbs = fZAbsorber-0.5*fAbsorberThickness; 135 fZEndAbs = fZAbsorber + 0.5 * fAbsorberThick << 137 fZEndAbs = fZAbsorber+0.5*fAbsorberThickness; 136 } 138 } 137 139 138 #endif 140 #endif 139 141