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