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