| 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 RE06/src/RE06ParallelWorld.cc 26 /// \file RE06/src/RE06ParallelWorld.cc
27 /// \brief Implementation of the RE06ParallelW 27 /// \brief Implementation of the RE06ParallelWorld class
28 // 28 //
29 // << 29 // $Id: RE06ParallelWorld.cc 83381 2014-08-21 14:18:13Z gcosmo $
>> 30 //
30 31
31 #include "RE06ParallelWorld.hh" 32 #include "RE06ParallelWorld.hh"
32 33
>> 34 #include "G4Tubs.hh"
33 #include "G4LogicalVolume.hh" 35 #include "G4LogicalVolume.hh"
34 #include "G4MultiFunctionalDetector.hh" <<
35 #include "G4PSEnergyDeposit.hh" <<
36 #include "G4PSMinKinEAtGeneration.hh" <<
37 #include "G4PSNofSecondary.hh" <<
38 #include "G4PSNofStep.hh" <<
39 #include "G4PSTrackLength.hh" <<
40 #include "G4PVPlacement.hh" 36 #include "G4PVPlacement.hh"
41 #include "G4PVReplica.hh" 37 #include "G4PVReplica.hh"
>> 38
42 #include "G4SDManager.hh" 39 #include "G4SDManager.hh"
43 #include "G4SDParticleFilter.hh" << 40 #include "G4MultiFunctionalDetector.hh"
44 #include "G4SystemOfUnits.hh" <<
45 #include "G4Tubs.hh" <<
46 #include "G4VPrimitiveScorer.hh" 41 #include "G4VPrimitiveScorer.hh"
>> 42 #include "G4PSEnergyDeposit.hh"
>> 43 #include "G4PSNofSecondary.hh"
>> 44 #include "G4PSTrackLength.hh"
>> 45 #include "G4PSNofStep.hh"
>> 46 #include "G4PSMinKinEAtGeneration.hh"
47 #include "G4VSDFilter.hh" 47 #include "G4VSDFilter.hh"
>> 48 #include "G4SDParticleFilter.hh"
48 #include "G4ios.hh" 49 #include "G4ios.hh"
>> 50 #include "G4SystemOfUnits.hh"
49 51
50 //....oooOO0OOooo........oooOO0OOooo........oo 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
51 53
52 G4ThreadLocal G4bool RE06ParallelWorld::fSDCon 54 G4ThreadLocal G4bool RE06ParallelWorld::fSDConstructed = false;
53 55
54 RE06ParallelWorld::RE06ParallelWorld(G4String 56 RE06ParallelWorld::RE06ParallelWorld(G4String worldName)
55 : G4VUserParallelWorld(worldName), << 57 :G4VUserParallelWorld(worldName),
56 fConstructed(false), << 58 fConstructed(false),
57 fSerial(false), << 59 fSerial(false),
58 fTotalThickness(2.0 * m), << 60 fTotalThickness(2.0*m),
59 fNumberOfLayers(20) << 61 fNumberOfLayers(20)
60 { 62 {
61 for (size_t i = 0; i < 3; i++) { << 63 for(size_t i=0;i<3;i++)
>> 64 {
62 fCalorLogical[i] = 0; 65 fCalorLogical[i] = 0;
63 fLayerLogical[i] = 0; 66 fLayerLogical[i] = 0;
64 fCalorPhysical[i] = 0; 67 fCalorPhysical[i] = 0;
65 fLayerPhysical[i] = 0; 68 fLayerPhysical[i] = 0;
66 } 69 }
67 fCalName[0] = "Calor-AP"; 70 fCalName[0] = "Calor-AP";
68 fCalName[1] = "Calor-BP"; 71 fCalName[1] = "Calor-BP";
69 fCalName[2] = "Calor-CP"; 72 fCalName[2] = "Calor-CP";
70 } 73 }
71 74
72 //....oooOO0OOooo........oooOO0OOooo........oo 75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73 76
74 RE06ParallelWorld::~RE06ParallelWorld() 77 RE06ParallelWorld::~RE06ParallelWorld()
75 { << 78 {;}
76 ; <<
77 } <<
78 79
79 //....oooOO0OOooo........oooOO0OOooo........oo 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
80 81
81 void RE06ParallelWorld::Construct() 82 void RE06ParallelWorld::Construct()
82 { 83 {
83 if (!fConstructed) { << 84 if(!fConstructed)
>> 85 {
84 fConstructed = true; 86 fConstructed = true;
85 SetupGeometry(); 87 SetupGeometry();
86 } 88 }
87 } 89 }
88 90
89 //....oooOO0OOooo........oooOO0OOooo........oo 91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
90 92
91 void RE06ParallelWorld::ConstructSD() 93 void RE06ParallelWorld::ConstructSD()
92 { 94 {
93 if (!fSDConstructed) { << 95 if(!fSDConstructed)
>> 96 {
94 fSDConstructed = true; 97 fSDConstructed = true;
95 SetupDetectors(); 98 SetupDetectors();
96 } 99 }
97 } 100 }
98 101
99 //....oooOO0OOooo........oooOO0OOooo........oo 102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
100 void RE06ParallelWorld::SetupGeometry() 103 void RE06ParallelWorld::SetupGeometry()
101 { 104 {
102 // << 105 //
103 // World 106 // World
104 // 107 //
105 G4VPhysicalVolume* ghostWorld = GetWorld(); 108 G4VPhysicalVolume* ghostWorld = GetWorld();
106 G4LogicalVolume* worldLogical = ghostWorld-> 109 G4LogicalVolume* worldLogical = ghostWorld->GetLogicalVolume();
107 << 110
108 // << 111 //
109 // Calorimeter 112 // Calorimeter
110 // << 113 //
111 G4VSolid* calorSolid = new G4Tubs("Calor", 0 << 114 G4VSolid* calorSolid
>> 115 = new G4Tubs("Calor",0.0,0.5*m,fTotalThickness/2.,0.0,360.*deg);
112 G4int i; 116 G4int i;
113 for (i = 0; i < 3; i++) { << 117 for(i=0;i<3;i++)
114 fCalorLogical[i] = new G4LogicalVolume(cal << 118 {
115 if (fSerial) { << 119 fCalorLogical[i] = new G4LogicalVolume(calorSolid,0,fCalName[i]);
116 fCalorPhysical[i] = << 120 if(fSerial)
117 new G4PVPlacement(0, G4ThreeVector(0., << 121 {
118 fCalorLogical[i], fC << 122 fCalorPhysical[i] = new G4PVPlacement(0,
>> 123 G4ThreeVector(0.,0.,G4double(i-1)*fTotalThickness),
>> 124 fCalorLogical[i],fCalName[i],worldLogical,false,i);
119 } 125 }
120 else { << 126 else
121 fCalorPhysical[i] = new G4PVPlacement(0, << 127 {
122 fC << 128 fCalorPhysical[i] = new G4PVPlacement(0,
>> 129 G4ThreeVector(0.,G4double(i-1)*m,0.),
>> 130 fCalorLogical[i],fCalName[i],worldLogical,false,i);
123 } 131 }
124 } 132 }
125 << 133
126 // << 134 //
127 // Layers --- as absorbers 135 // Layers --- as absorbers
128 // 136 //
129 G4VSolid* layerSolid = new G4Tubs("Layer", 0 << 137 G4VSolid* layerSolid
130 for (i = 0; i < 3; i++) { << 138 = new G4Tubs("Layer",0.0,0.5*m,fTotalThickness/2.,0.0,360.*deg);
131 fLayerLogical[i] = new G4LogicalVolume(lay << 139 for(i=0;i<3;i++)
132 fLayerPhysical[i] = new G4PVReplica(fCalNa << 140 {
133 kRho, << 141 fLayerLogical[i]
>> 142 = new G4LogicalVolume(layerSolid,0,fCalName[i]+"_LayerLog");
>> 143 fLayerPhysical[i]
>> 144 = new G4PVReplica(fCalName[i]+"_Layer",fLayerLogical[i],fCalorLogical[i],
>> 145 kRho,fNumberOfLayers,0.5*m/fNumberOfLayers);
134 } 146 }
>> 147
135 } 148 }
136 149
137 //....oooOO0OOooo........oooOO0OOooo........oo 150 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
138 151
139 void RE06ParallelWorld::SetupDetectors() 152 void RE06ParallelWorld::SetupDetectors()
140 { 153 {
141 G4SDManager::GetSDMpointer()->SetVerboseLeve 154 G4SDManager::GetSDMpointer()->SetVerboseLevel(1);
142 G4String filterName, particleName; 155 G4String filterName, particleName;
143 156
144 G4SDParticleFilter* gammaFilter = << 157 G4SDParticleFilter* gammaFilter
145 new G4SDParticleFilter(filterName = "gamma << 158 = new G4SDParticleFilter(filterName="gammaFilter",particleName="gamma");
146 G4SDParticleFilter* electronFilter = << 159 G4SDParticleFilter* electronFilter
147 new G4SDParticleFilter(filterName = "elect << 160 = new G4SDParticleFilter(filterName="electronFilter",particleName="e-");
148 G4SDParticleFilter* positronFilter = << 161 G4SDParticleFilter* positronFilter
149 new G4SDParticleFilter(filterName = "posit << 162 = new G4SDParticleFilter(filterName="positronFilter",particleName="e+");
150 G4SDParticleFilter* epFilter = new G4SDParti << 163 G4SDParticleFilter* epFilter
151 epFilter->add(particleName = "e-"); << 164 = new G4SDParticleFilter(filterName="epFilter");
152 epFilter->add(particleName = "e+"); << 165 epFilter->add(particleName="e-");
153 << 166 epFilter->add(particleName="e+");
154 for (G4int i = 0; i < 3; i++) { << 167
155 G4String detName = fCalName[i] + "_para"; << 168 for(G4int i=0;i<3;i++)
>> 169 {
>> 170 G4String detName = fCalName[i]+"_para";
156 G4MultiFunctionalDetector* det = new G4Mul 171 G4MultiFunctionalDetector* det = new G4MultiFunctionalDetector(detName);
157 172
158 G4VPrimitiveScorer* primitive; 173 G4VPrimitiveScorer* primitive;
159 primitive = new G4PSEnergyDeposit("eDep"); 174 primitive = new G4PSEnergyDeposit("eDep");
160 det->RegisterPrimitive(primitive); 175 det->RegisterPrimitive(primitive);
161 primitive = new G4PSNofSecondary("nGamma") 176 primitive = new G4PSNofSecondary("nGamma");
162 primitive->SetFilter(gammaFilter); 177 primitive->SetFilter(gammaFilter);
163 det->RegisterPrimitive(primitive); 178 det->RegisterPrimitive(primitive);
164 primitive = new G4PSNofSecondary("nElectro 179 primitive = new G4PSNofSecondary("nElectron");
165 primitive->SetFilter(electronFilter); 180 primitive->SetFilter(electronFilter);
166 det->RegisterPrimitive(primitive); 181 det->RegisterPrimitive(primitive);
167 primitive = new G4PSNofSecondary("nPositro 182 primitive = new G4PSNofSecondary("nPositron");
168 primitive->SetFilter(positronFilter); 183 primitive->SetFilter(positronFilter);
169 det->RegisterPrimitive(primitive); 184 det->RegisterPrimitive(primitive);
170 primitive = new G4PSTrackLength("trackLeng 185 primitive = new G4PSTrackLength("trackLength");
171 primitive->SetFilter(epFilter); 186 primitive->SetFilter(epFilter);
172 det->RegisterPrimitive(primitive); 187 det->RegisterPrimitive(primitive);
173 primitive = new G4PSNofStep("nStep"); 188 primitive = new G4PSNofStep("nStep");
174 primitive->SetFilter(epFilter); 189 primitive->SetFilter(epFilter);
175 det->RegisterPrimitive(primitive); 190 det->RegisterPrimitive(primitive);
176 191
177 G4SDManager::GetSDMpointer()->AddNewDetect 192 G4SDManager::GetSDMpointer()->AddNewDetector(det);
178 SetSensitiveDetector(fLayerLogical[i], det 193 SetSensitiveDetector(fLayerLogical[i], det);
179 } 194 }
180 G4SDManager::GetSDMpointer()->SetVerboseLeve 195 G4SDManager::GetSDMpointer()->SetVerboseLevel(0);
>> 196
181 } 197 }
182 198
183 //....oooOO0OOooo........oooOO0OOooo........oo 199 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
184 200
185 void RE06ParallelWorld::SetSerialGeometry(G4bo 201 void RE06ParallelWorld::SetSerialGeometry(G4bool serial)
186 { 202 {
187 if (fSerial == serial) return; << 203 if(fSerial==serial) return;
188 fSerial = serial; << 204 fSerial=serial;
189 if (!fConstructed) return; << 205 if(!fConstructed) return;
190 for (G4int i = 0; i < 3; i++) { << 206 for(G4int i=0;i<3;i++)
191 if (fSerial) { << 207 {
192 fCalorPhysical[i]->SetTranslation(G4Thre << 208 if(fSerial)
>> 209 {
>> 210 fCalorPhysical[i]
>> 211 ->SetTranslation(G4ThreeVector(0.,0.,G4double(i-1)*2.*m));
193 } 212 }
194 else { << 213 else
195 fCalorPhysical[i]->SetTranslation(G4Thre << 214 {
>> 215 fCalorPhysical[i]
>> 216 ->SetTranslation(G4ThreeVector(0.,G4double(i-1)*m,0.));
196 } 217 }
197 } 218 }
198 } 219 }
199 220
200 //....oooOO0OOooo........oooOO0OOooo........oo 221 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
201 222