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