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