LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Adapted LEM mac files to Geant4 10.3

Browse Source
This commit is contained in:
prokscha
2017-03-18 16:17:41 +01:00
parent b613c0dde1
commit ff4cccadd1
5 changed files with 2855 additions and 26 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
run/LEM/templates/Template_15kV_muons.mac
View File

@ -414,7 +414,8 @@
# trigger foil is 52-55mm upstream of TrigV center
# Carbon Foil (default HALF-thickness 0.000005147 mm, see below => CFoil thick = 10.3 nm).
# USE THE NAME CFoil or coulombCFoil, otherwise musrMuFormation won't work!
/musr/command construct box CFoil 60 60 0.0000025 G4_GRAPHITE 0 0 -550 log_World norot dead 621
#/musr/command construct box CFoil 60 60 0.0000025 G4_GRAPHITE 0 0 -550 log_World norot dead 621
/musr/command construct box CFoil 60 60 0.000002 G4_GRAPHITE 0 0 -550 log_World norot dead 621
# Notes: NIST tables use G4_GRAPHITE with 1.7 g/cm3 and 78 eV ioniz. energy.
# An area density of 2.20 ug/cm2 implies a CF thickn. = (2.20*1.e-6/1.70)*cm = 1.294e-5 mm - Total thickness
@ -689,13 +690,13 @@
/musr/command globalfield Lens1EField 251.5 0.0 -1678.0 fromfile 2DE L1_Erz.map log_L1EField 9.0
# LENS 2. Best L3 settings: 7.0 kV for 10 keV muons
/musr/command globalfield Lens3EField 0. 0. -863.0 fromfile 2DE L3_Erz.map log_L3EField 9.5
/musr/command globalfield Lens3EField 0. 0. -863.0 fromfile 2DE L3_Erz.map log_L3EField 10.5
# Spin rotator
# Electric field should be 370mm plates
# Magnetic field should be 260mm poles as measured srlem_ext.map
/musr/command globalfield Ele_field 0. 0. -1284.5 fromfile 3DE E_Sep_P37.5cm.map log_SREField 2.57
/musr/command globalfield Mag_field 0. 0. -1284.51 fromfile 3DB srlem_ext.map log_SRBField -0.006923
/musr/command globalfield Ele_field 0. 0. -1284.5 fromfile 3DE E_Sep_P37.5cm.map log_SREField 2.717
/musr/command globalfield Mag_field 0. 0. -1284.51 fromfile 3DB srlem_ext.map log_SRBField -0.00692
# Electrical Field areas in the Trigger Detector
# En = Electrical Field n: TnFieldMgr (n = 1-3)
@ -721,10 +722,10 @@
### Electric field at RING ANODE - from 3DE field map
# To create an arbitrary configuration, switch on all fields and set different potentials.
/musr/command globalfield RngAnT_field 0. 0. -143.00 fromfile 3DE EM_3D_ext_gridf.map log_RA_T 0.0
/musr/command globalfield RngAnR_field 0. 0. -143.02 fromfile 3DE EM_3D_ext_gridf.map log_RA_R 0.0
/musr/command globalfield RngAnB_field 0. 0. -143.04 fromfile 3DE EM_3D_ext_gridf.map log_RA_B 0.0
/musr/command globalfield RngAnL_field 0. 0. -143.06 fromfile 3DE EM_3D_ext_gridf.map log_RA_L 0.0
/musr/command globalfield RngAnT_field 0. 0. -143.00 fromfile 3DE EM_3D_ext_gridf.map log_RA_T 11.9
/musr/command globalfield RngAnR_field 0. 0. -143.02 fromfile 3DE EM_3D_ext_gridf.map log_RA_R 11.9
/musr/command globalfield RngAnB_field 0. 0. -143.04 fromfile 3DE EM_3D_ext_gridf.map log_RA_B 11.9
/musr/command globalfield RngAnL_field 0. 0. -143.06 fromfile 3DE EM_3D_ext_gridf.map log_RA_L 11.9
### Magnetic field at the sample
/musr/command construct tubs SBField 0 0.01 0.005 0 360 G4_Galactic 0.0 0.0 15.0 log_World norot dead 709 nofield
@ -750,7 +751,7 @@
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
# Geant 4.10
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
@ -777,7 +778,8 @@
#/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
#/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#/musr/command process addDiscreteProcess mu+ G4CoulombScattering
#/musr/command process addProcess mu+ G4StepLimiter -1 -1 5
#/musr/command process addProcess mu+ G4UserSpecialCuts -1 -1 3
#/musr/command process addProcess mu+ G4StepLimiter -1 -1 4
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
@ -810,7 +812,11 @@
# Set particle energy cuts on particular volumes (in eV)
#/musr/command SetUserLimits log_World ustepMax(mm) utrakMax(mm) utimeMax(ns) uekinMin(MeV) urangMin(mm)
/musr/command SetUserLimits log_World -1 -1 -1 1e-7 -1
#/musr/command SetUserLimits log_World -1 -1 -1 1e-7 -1
# Geant4 10.2
#------------
/process/em/lowestMuHadEnergy 1e-6 MeV
# Store ALL the events in a ROOT tree or just the interesting ones? (default is true)
#*/musr/command storeOnlyEventsWithHits false
@ -1000,7 +1006,7 @@
#*/hits/inactivate /musr/ScintSD
# Only for code debugging!
#/tracking/verbose 1
/tracking/verbose 0
# BEAM ON
#/run/beamOn 1000000