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              
 7              begin parm experiment
 8                ngen = 1000			;  POS: # of successes; NEG: # of tries
 9                EXPER.charge =  1.0		;  total charge (mC)
10                doing_phsp = 0		;  (ONE = TRUE)	- If all of the doing_* are
11                doing_kaon = 0		;  (ONE = TRUE)   false, then doing (e,e'p).
12                doing_pion = 0		;  (ONE = TRUE)
13                doing_decay = 0		;  1=decay ON, 0=decay OFF.
14                ctau = 780.4			;  decay length (cm)
15                extra_dbase_file='extra_hydrogen' ;  additional input parameters.
16              end parm experiment
17              
18              begin parm kinematics_main
19                Ebeam = 3056.			;  (MeV)
20                dEbeam = 0.05			;  beam energy variation (%)
21                electron_arm = 1              ;  1=hms,2=sos,3=hrsr,4=hrsl
22 gaskelld 1.1   hadron_arm = 2                ;  1=hms,2=sos,3=hrsr,4=hrsl
23                spec.e.P = 1300.0		;  e arm central momentum (MeV/c)
24                spec.e.theta = 54.0		;  e arm angle setting (degrees)
25                spec.p.P = 2520.0		;  p arm central momentum (MeV/c)
26                spec.p.theta = 24.80		;  p arm angle setting (degrees)
27              end parm kinematics_main
28              
29              begin parm target
30                targ.A = 56.0			;  target A
31                targ.Z = 26.0			;  target Z
32                targ.mass_amu = 55.847	;  target mass in amu
33                targ.mrec_amu = 54.847	;  recoil mass in amu (eep=A-1 system,pion=A-2)
34                targ.rho = 7.87		;  target density (g/cm^3)
35                targ.thick = 212.9		;  target thick (mg/cm^2)
36                targ.angle = 0.		;  target angle (for solid target) (degrees)
37                targ.abundancy = 100.		;  target purity (%)
38                targ.can = 1                  ;  1=beer can (fpi), 2=pudding can (nucpi)
39              end parm target
40              
41              begin parm debug		;  (ONES give helpful debug info)
42                debug(1) = 0			;  turns on output from brem.f
43 gaskelld 1.1   debug(2) = 0			;  into/outa subs.
44                debug(3) = 0			;  spit out values (init. and main loop).
45                debug(4) = 0			;  mostly comp_ev, gen_rad diagnostics.
46                debug(5) = 0			;  a bit of everything.
47              end parm debug
48              
49              begin parm e_arm_accept
50                SPedge.e.delta.min = -10.0	;  delta min (SPECTROMETER ACCEPTANCE!)
51                SPedge.e.delta.max =  10.0	;  delta max
52                SPedge.e.yptar.min = -50.0	;  yptar min = {TF} / 1000 (mrad)
53                SPedge.e.yptar.max =  50.0	;  yptar max = {TF} / 1000
54                SPedge.e.xptar.min = -100.0	;  xptar min = {TF} / 1000 (mrad)
55                SPedge.e.xptar.max =  100.0	;  xptar max = {TF} / 1000
56              end parm e_arm_accept
57              
58              begin parm p_arm_accept
59                SPedge.p.delta.min = -15.0	;  delta min (SPECTROMETER ACCEPTANCE!)
60                SPedge.p.delta.max =  15.0	;  delta max
61                SPedge.p.yptar.min = -90.0	;  yptar min = {TF} / 1000 (mrad)
62                SPedge.p.yptar.max =  90.0	;  yptar max = {TF} / 1000
63                SPedge.p.xptar.min = -50.0	;  xptar min = {TF} / 1000 (mrad)
64 gaskelld 1.1   SPedge.p.xptar.max =  50.0	;  xptar max = {TF} / 1000
65              end parm p_arm_accept