Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // Gorad (Geant4 Open-source Radiation Analysis and Design) 27 // 28 // Author : Makoto Asai (SLAC National Accelerator Laboratory) 29 // 30 // Development of Gorad is funded by NASA Johnson Space Center (JSC) 31 // under the contract NNJ15HK11B. 32 // 33 // ******************************************************************** 34 // 35 // GRGeomImpBiasWorld.cc 36 // A parallel world class that defines the geometry 37 // of the geometry improtance biasing. 38 // 39 // History 40 // September 8th, 2020 : first implementation 41 // 42 // ******************************************************************** 43 44 #include "GRGeomImpBiasWorld.hh" 45 #include "GRDetectorConstruction.hh" 46 47 #include "G4Orb.hh" 48 #include "G4LogicalVolume.hh" 49 #include "G4VPhysicalVolume.hh" 50 #include "G4PVPlacement.hh" 51 52 #include "G4Region.hh" 53 #include "GRBiasingRegionInfo.hh" 54 55 #include "G4UIcommand.hh" 56 #include "G4VisAttributes.hh" 57 58 GRGeomImpBiasWorld::GRGeomImpBiasWorld(G4String& wName,GRDetectorConstruction* det) 59 : G4VUserParallelWorld(wName),detector(det) 60 { 61 stateNotifier = new GRGeomImpBiasWorldStateNotifier(this); 62 } 63 64 GRGeomImpBiasWorld::~GRGeomImpBiasWorld() 65 { 66 delete stateNotifier; 67 } 68 69 void GRGeomImpBiasWorld::Construct() 70 { 71 if(constructed) return; 72 73 constructed = true; 74 fWorld = GetWorld(); //physical volume of the world volume of parallel world 75 G4LogicalVolume* motherLog = fWorld->GetLogicalVolume(); 76 biasingRegion = new G4Region(fWorldName+"_Region"); 77 auto biasInfo = new GRBiasingRegionInfo(); 78 biasingRegion->SetUserInformation(biasInfo); 79 biasingRegion->AddRootLogicalVolume(motherLog); 80 biasingRegion->SetWorld(fWorld); 81 G4VisAttributes* wvisatt = new G4VisAttributes(G4Colour(.2,.2,0.)); 82 wvisatt->SetVisibility(false); 83 motherLog->SetVisAttributes(wvisatt); 84 85 G4double r0 = detector->geoImpP.radius; 86 if(r0 < 0.) 87 { 88 // radius of the outermost sphere is not specified, so seting it to the default 89 // value as the 80% of the world volume. 90 r0 = GRDetectorConstruction::GetWorldSize() * 0.8; 91 G4cout<<"############ Radius of the outermost biasing sphere is set to "<<r0<<" (mm)"<<G4endl; 92 } 93 94 G4int nL = detector->geoImpP.nLayer; 95 G4ThreeVector dp = (detector->geoImpP.posT - detector->geoImpP.pos0) / (nL-1); 96 G4double rt = detector->geoImpP.radiusT; 97 if(rt<0.) 98 { rt = r0/nL; } 99 G4double dr = (r0 - rt)/(nL-1); 100 101 for(G4int i=0;i<nL;i++) 102 { 103 G4double r = r0 - dr*i; 104 G4String vName = fWorldName + "_" + G4UIcommand::ConvertToString(i); 105 auto sph = new G4Orb(vName+"_solid",r); 106 auto lv = new G4LogicalVolume(sph,nullptr,vName+"_lv"); 107 G4ThreeVector pos = detector->geoImpP.pos0; 108 if(i!=0) pos = dp; 109 new G4PVPlacement(0,pos,lv,vName+"_pv",motherLog,false,i+1); 110 motherLog = lv; 111 G4VisAttributes* visatt = new G4VisAttributes(G4Colour(.7,.7,0.)); 112 visatt->SetVisibility(true); 113 lv->SetVisAttributes(visatt); 114 } 115 116 return; 117 } 118 119 void GRGeomImpBiasWorld::ConstructSD() 120 { ; } 121 122 void GRGeomImpBiasWorld::Update() 123 { 124 GRBiasingRegionInfo* biasInfo = static_cast<GRBiasingRegionInfo*>(biasingRegion->GetUserInformation()); 125 biasInfo->SetBiasingFactor(detector->geoImpP.factor); 126 biasInfo->SetProbability(detector->geoImpP.prob); 127 } 128 129 130