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