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

Updated geant_lemsr.input.

Browse Source
This commit is contained in:
prokscha 2007-02-09 08:06:01 +00:00
parent 9494ea7b29
commit bea4687d35
1 changed files with 17 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
geant3/inp/geant_lemsr.input
View File

@ -1,17 +1,22 @@
plot 1 plot 0
rung 3022 rung 3138
trig 10 trig 1000000
rndm 9871 59367 rndm 9871 59367
COMM 1 2 3 4 5 6 7 8 9 COMM 1 2 3 4 5 6 7 8 9
sets 5 -74 1 9 0 0 0 0 0 COMM sets 500 -74 1 9 0 0 0 0 0
bfie 0. 0. 0. 0. 0. sets 2 0 0 0 0 0 0 0 0
bfie 0. 0. 0. 100. 100.
spar ' e+' ' e-' ' gam' spar ' e+' ' e-' ' gam'
cuts 0.0001 0.0001 0.01 0.01 0.00001 cuts 0.0001 0.0001 0.01 0.01 0.00001
COMM CutGam CutEle CutNeu CutHad CutMuo COMM CutGam CutEle CutNeu CutHad CutMuo
ntva 'init' 'time' 'flag' 'scie' 'scoe' 'ltrb' ntva 'init' 'time' 'flag' 'scie' 'scoe' 'ltrb'
geom 'mcps' 'samp' 'efla' 'cryo' 'crsh' 'saal' 'ru11' COMM geom 'mcps' 'samp' 'efla' 'cryo' 'crsh'
geom 'mcps' 'samp' 'efla' 'cryo' 'crsh' 'saal'
COMM geom 'mcps' 'efla' 'ru11' 'mcp2' 'mcpa' COMM geom 'mcps' 'efla' 'ru11' 'mcp2' 'mcpa'
COMM geom 'mcps'
view 1 60 40 0 10 10 4 4 view 1 60 40 0 10 10 4 4
bpar 0.8 0. 0.5 0.001
COMM dx0 dy0 r*R_MCP2 RelaxRate
stop stop
@ -36,7 +41,7 @@ sets 5 -74 2 5 0 1 0 2 5 ---> first number : particle ID, 5 = mu+ (2 = e+)
ID = 500; user defined mu+ ID = 500; user defined mu+
with anistropic decay prob. with anistropic decay prob.
second number: if == 0, michel-distributed at second number: if == 0, michel-distributed at
MCP2 position MCP2 position ( use particle id 2, e+!!!)
if ne 0, start position is at z=sets(2) if ne 0, start position is at z=sets(2)
(here 74 cm upstream of MCP tube center) (here 74 cm upstream of MCP tube center)
with beam momentum with beam momentum
@ -48,8 +53,9 @@ sets 5 -74 2 5 0 1 0 2 5 ---> first number : particle ID, 5 = mu+ (2 = e+)
sets(8): gaussian dp/p distribution, sets(8): gaussian dp/p distribution,
sigma_p = sets(8).sets(9)% sigma_p = sets(8).sets(9)%
spar ' e+' ' e-' ' gam' --> select secondary particles to be tracked. spar ' e+' ' e-' ' gam' --> select secondary particles to be tracked.
bpar 0.25 -0.1 0.55 0.001 ---> beam offset: x0+0.25cm, y0-0.1cm, r=0.55*R_MCP2 (2cm)
BFIE 0. 0. 100. 90. 50. ---> Sample B-field (Gauss), Bx, By, Bz, and muon polarization=90% exponential relaxation 1MHz = 0.001/ns
BFIE 0. 0. 100. 90. 50. ---> Sample B-field (Gauss), Bx, By, Bz, and muon polarization=90%
and muon polarization=50% in sample. and muon polarization=50% in sample.
ntva 'time' 'flag' ---> see comments below ntva 'time' 'flag' ---> see comments below
@ -98,6 +104,7 @@ c
c 'init' write initial parameter to NT c 'init' write initial parameter to NT
c 'code' initial particle ID to NT c 'code' initial particle ID to NT
c 'endp' end position and energy of e+ c 'endp' end position and energy of e+
c 'endm' end position and energy of muon (particle code 500)
c 'time' t-mcp and t-sci c 'time' t-mcp and t-sci
c 'scip' position when particle hits inner Sc. c 'scip' position when particle hits inner Sc.
c 'scie' energy loss of particles in inner sc. c 'scie' energy loss of particles in inner sc.