plot 0 rung 3138 trig 1000000 rndm 9871 59367 COMM 1 2 3 4 5 6 7 8 9 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' 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 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 ( 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 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. 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 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 '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. 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--------------------------------------------------------------------------------