# Simplest TOPAS example. # A box inside another box with a different camera view (see # last block of lines) # A simple unidirectional, linear and monoenergetic beam #====== #Materials #====== #Se utilizan solo los materiales ya implementados en TOPAS #======== #Geometry #======== #****************************** #**** World ****** #****************************** b:Ge/World/Invisible = "True" d:Ge/World/HLX = 2. m d:Ge/World/HLY = 2. m d:Ge/World/HLZ = 2. m s:Ge/World/Material = "Vacuum" #****************************** #**** Target ****** #****************************** s:Ge/MyBox/Type = "TsBox" dc:Ge/MyBox/HLX = 1 mm dc:Ge/MyBox/HLY = 1 mm dc:Ge/MyBox/HLZ = 1 um sc:Ge/MyBox/Material = "G4_Au" d:Ge/MyBox/TransX = 0. m d:Ge/MyBox/TransY = 0. m d:Ge/MyBox/TransZ = 0. m d:Ge/MyBox/RotX = 0. deg d:Ge/MyBox/RotY = 0. deg d:Ge/MyBox/RotZ = 0. deg s:Ge/MyBox/Parent = "World" #****************************** #**** Cylindrical detector ****** #****************************** s:Ge/Detector/Type = "TsCylinder" d:Ge/Detector/RMin = 5.0 mm d:Ge/Detector/RMax = 5.1 mm d:Ge/Detector/SPhi = 5.0 deg d:Ge/Detector/DPhi = 350. deg d:Ge/Detector/HL = 1 mm s:Ge/Detector/Material = "Vacuum" s:Ge/Detector/Parent = "World" ic:Ge/Detector/PhiBins = 70 dc:Ge/Detector/RotX = 90 deg dc:Ge/Detector/RotY = 0 deg dc:Ge/Detector/RotZ = -90 deg #====================== # Source component # Ubicación de la fuente #====================== s:Ge/BeamPosition/Parent = "World" s:Ge/BeamPosition/Type = "Group" d:Ge/BeamPosition/TransX = 0. m d:Ge/BeamPosition/TransY = 0. m d:Ge/BeamPosition/TransZ = -6 mm d:Ge/BeamPosition/RotX = 0. deg d:Ge/BeamPosition/RotY = 0. deg d:Ge/BeamPosition/RotZ = 0. deg #====================== # Primary source (beam) #====================== #s:So/MySource/Type = "Beam" s:So/MySource/Type = "Isotropic" s:So/MySource/Component = "BeamPosition" s:So/MySource/BeamParticle = "alpha" dc:So/MySource/BeamEnergy = 5.5 MeV u:So/MySource/BeamEnergySpread = 0.5 s:So/MySource/BeamPositionDistribution = "None" #s:So/MySource/BeamPositionCutoffShape = "Ellipse" #d:So/MySource/BeamPositionCutoffX = 0.1 mm # X extent of position (if Flat or Gaussian) #d:So/MySource/BeamPositionCutoffY = 0.1 mm #s:So/MySource/BeamAngularDistribution = "Flat" #d:So/MySource/BeamAngularCutoffX = 1. deg #d:So/MySource/BeamAngularCutoffY = 1. deg i:Ts/ShowHistoryCountAtInterval = 1 ic:So/MySource/NumberOfHistoriesInRun = 1000 #====================== # Physics List #====================== sv:Ph/Default/Modules = 1 "g4em-standard_opt4" #sv:Ph/Default/Modules = 1 "Transportation_Only" #==================== #Graphical part #==================== s:Gr/ViewA/Type = "OpenGL" Ts/UseQt = "True" i:Gr/ViewA/WindowSizeX = 1024 i:Gr/ViewA/WindowSizeY = 768 #b:Gr/ViewA/IncludeAxes = "True" #s:Gr/ViewA/AxesComponent = "World" #Se puede ajustar al volumen requerido #d:Gr/ViewA/AxesSize = 0.1 m # #***** Camera view ******* d:Gr/ViewA/Theta = -90 deg #d:Gr/ViewA/Phi = 150 deg #s:Gr/ViewA/Projection = "Perspective" #d:Gr/ViewA/PerspectiveAngle = 30 deg u:Gr/ViewA/Zoom = 1.3 #b:Ts/PauseBeforeQuit = "True" i:Ts/TrackingVerbosity = 0 i:Ts/NumberOfThreads = 1