plot 1
rung 3022
trig 10
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.
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'
COMM geom 'mcps' 'efla' 'ru11' 'mcp2' 'mcpa'  
view 1 60 40 0 10 10 4 4
stop


This is the Format Free Input File for GEANT_LEMSR.FOR:
=======================================================

how to use the data cards:

COMM            ---> comment line, no effect when reading

plot 0          ---> no graphics.
plot 1          ---> graphics on screen, and also information on detector hits.
plot 2          ---> graphics on screen and PS file.

rung 12         ---> define GEANT run number; will be appended to the NTP filename

trig 100000     ---> number of events to be processed; on PSICL1 10^5 takes
                     up to 3:30 min CPU time.

     1   2 3 4 5 6 7 8 9
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
                                        if ne 0, start position is at z=sets(2) 
                                        (here 74 cm upstream of MCP tube center)
                                         with beam momentum 
                                         sets(3).sets(4)sets(5),
                                         this is 2.50 MeV/c in this example
                             sets(6):   gaussian beam divergence (in degree)
                             sets(7):   if = 1, read initial values from file
                                        MC$INP:GEANT_LEMSR_INPUT.DAT
			     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%
                   		  and muon polarization=50% in sample.

ntva 'time' 'flag'        --->   see comments below


view 1 40 40 0 10 10 3 3  ---> in case of graphical output define view angle,
     (60)(40)(0)(10)(10)(2)(2) position and scaling of graph.
        default values         view(1): if == 1 then use the command lines
                                        parameters for graphical output;
                                        otherwise use default parameter
                               view(2): theta angle (degree)
                               view(3): phi angle 
                               view(4): psi angle
                               view(5): x0, x-position in HIGZ window
                               view(6): y0, y-position in HIGZ window
                               view(7): 10*x-scaling,  view(7)=3 <==> xscal=0.3
                               view(8): 10*y-scaling,  view(8)=3 <==> yscal=0.3
geom 'mcp2' 'efla'      ---> see comment below

stop   ---> stops reading this file and continues the program.

c------------------------------------------------------------------------------
c
c       flags for geometrical setup; it is possible to mount or remove the
c
c       'mcp2'  MCP2 plates
c       'mcpa'  MCP2 anode
c       'mcps'  MCP2 stainless steel vacuum tube
c       'samp'  sample  (this is at the moment the sample holder)
c       'cryo'  material of the cryostat without shielding
c       'crsh'  cryo shield
c       'efla'  the 100 CF flange at the end of the sample/MCP2 tube
c       'ru11'  use MCP setup of Run 11 (delay line anode)
c       'ru10'  use MCP setup of Run 10 and before (WSZ detector)
c               has no influence for cryo geometry              
c       'saal'  uses Al instead of Cu for sample holder plates
c
c       the volumes are activated when printing
c
c       geom 'mcp2' 'efla'  for example when data cards are requested
c
c---------------------------------------
c
c       flags to define NTuple variables for CWN
c
c       'init'  write initial parameter to NT
c       'code'  initial particle ID to NT
c       'endp'  end position and energy of e+
c       'time'  t-mcp and t-sci
c       'scip'  position when particle hits inner Sc.
c       'scie'  energy loss of particles in inner sc.
c       'scop'  position when particle hits outer Sc.
c       'scoe'  energy loss of particles in outer sc.
c       'ltrb'  energy losses in single detectors (inner and outer sc.)
c       'flag'  Vol. numbers and particle codes; needed to get the single
c               decay spectra
c       'nhit'  Number of Hits in detectors MCP, SCI and SCO
c       'posm'  Generated momentum of decay positron
c       'mcpe'  energy loss in MCP detector
c       'allv'  use all possible NT variables
c
c--------------------------------------------------------------------------------
