diff --git a/geant3/inp/geant_lemsr.input b/geant3/inp/geant_lemsr.input
index 5385139..6f963e2 100644
--- a/geant3/inp/geant_lemsr.input
+++ b/geant3/inp/geant_lemsr.input
@@ -1,17 +1,22 @@
-plot 1
-rung 3022
-trig 10
+plot 0
+rung 3138
+trig 1000000
 rndm 9871 59367
 COMM  1   2  3 4 5 6 7 8 9
-sets 5  -74  1 9 0 0 0 0 0
-bfie 0. 0. 0. 0. 0.
+COMM sets 500 -74 1 9 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'   
 cuts 0.0001  0.0001  0.01   0.01 0.00001
 COMM CutGam CutEle CutNeu CutHad CutMuo
-ntva 'init' 'time' 'flag' 'scie' 'scoe' 'ltrb' 
-geom 'mcps' 'samp'  'efla'  'cryo'  'crsh' 'saal' 'ru11'
+ntva 'init' 'time' 'flag' 'scie' 'scoe' 'ltrb'
+COMM geom 'mcps' 'samp'  'efla'  'cryo'  'crsh'
+geom 'mcps' 'samp'  'efla'  'cryo'  'crsh' 'saal'
 COMM geom 'mcps' 'efla' 'ru11' 'mcp2' 'mcpa'  
+COMM geom 'mcps'
 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
 
 
@@ -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+
                                                 with anistropic decay prob.
                              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) 
                                         (here 74 cm upstream of MCP tube center)
                                          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,
 					sigma_p = sets(8).sets(9)%
 spar '  e+' '  e-' ' gam'       --> select secondary particles to be tracked.
-
-BFIE 0. 0. 100.  90. 50.   --->   Sample B-field (Gauss), Bx, By, Bz, and muon polarization=90%
+bpar 0.25 -0.1 0.55 0.001 ---> beam offset: x0+0.25cm, y0-0.1cm, r=0.55*R_MCP2 (2cm)
+                                            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.
 
 ntva 'time' 'flag'        --->   see comments below
@@ -98,6 +104,7 @@ c
 c       'init'  write initial parameter to NT
 c       'code'  initial particle ID to NT
 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       'scip'  position when particle hits inner Sc.
 c       'scie'  energy loss of particles in inner sc.