1 gaskelld 1.1 ; This is a CTP file
 2              
 3              begin parm beamandtargetinfo
 4                gen.xwid = 0.008868		;  beam width - one sigma (cm)  (89microns)
 5                gen.ywid = 0.004235		;  beam width - one sigma (cm)  (42microns)
 6                targ.fr_pattern = 2.		;  raster pattern: 1=square, 2=circular
 7                targ.fr1 = 0.015		;  horizontal size OR inner radius(2)
 8                targ.fr2 = 0.15		;  vertical size OR outer radius(2)
 9                targ.xoffset = 0.		;  target x-offset (cm): +x = beam right
10                targ.yoffset = 0.		;  target y-offset (cm): +y = up
11                targ.zoffset = 0.		;  target z-offset (cm): +z = downstream
12              				;   zreal = znominal + zoffset
13              end parm beamandtergetinfo
14              
15              ;These are offsets applied before the call to the single arm montecarlos.
16              begin parm spect_offset
17                spec.e.offset.x = 0.		;  x offset (cm)
18                spec.e.offset.y = 0.		;  y offset (cm)
19                spec.e.offset.z = 0.		;  z offset (cm)
20                spec.e.offset.xptar = 0.	;  xptar offset (mr)	!x(y)ptar is slope, so
21                spec.e.offset.yptar = 0.	;  yptar offset (mr)	!it's really unitless.
22 gaskelld 1.1   spec.p.offset.x = 0.		;  x offset (cm)
23                spec.p.offset.y = 0.		;  y offset (cm)
24                spec.p.offset.z = 0.		;  z offset (cm)
25                spec.p.offset.xptar = 2.62	;  xptar offset (mr)
26                spec.p.offset.yptar = 0.	;  yptar offset (mr)
27              end parm spect_offset
28              
29              begin parm simulate
30                hard_cuts = 0         ;  (ONE = TRUE) SPedge and Em.max are hard cuts(ntuple)
31                using_rad = 1		;  (ONE = TRUE)
32                use_expon = 0		;  (LEAVE AT 0)
33                one_tail = 0		;  0=all, 1=e, 2=e', 3=p, -3=all but p
34                intcor_mode = 1	;  (LEAVE AT 1)
35                spect_mode = 0	;  0=e+p arms, -1=p arm, -2=e arm only, 1=none
36                cuts.Em.min = 0.	;  (Em.min=Em.max=0.0 gives wide open cuts)
37                cuts.Em.max = 0.	; Must be wider than cuts in analysis(elastic or e,e'p)
38                using_Eloss = 1	;  (ONE = TRUE)
39                correct_Eloss = 1	;  ONE = correct reconstructed events for eloss.
40                correct_raster = 1	;  ONE = Reconstruct events using 'raster' matrix elements.
41                mc_smear = 1		;  ONE = target & hut mult scatt AND DC smearing.
42                deForest_flag = 0	;  0=sigcc1, 1=sigcc2, -1=sigcc1 ONSHELL
43 gaskelld 1.1   rad_flag = 0		;  (radiative option #1...see init.f)
44                extrad_flag = 2	;  (rad. option #2...see init.f)
45                lambda(1) = 0.0	;  if rad_flag.eq.4 then lambda(1) = {TF}
46                lambda(2) = 0.0	;  if rad_flag.eq.4 then lambda(2) = {TF}
47                lambda(3) = 0.0	;  if rad_flag.eq.4 then lambda(3) = {TF}
48                Nntu = 1		;  ONE = generate ntuples
49                using_Coulomb = 1	;  (ONE = TRUE)
50                dE_edge_test = 0.	;  (move around energy edges)
51                use_offshell_rad = 1	;  (ONE = TRUE)
52                Egamma_gen_max = 0. ;  Set >0 to hardwire the Egamma limits.
53              end parm simulate