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)
11                doing_kaon = 1		;  (ONE = TRUE)
12                which_kaon = 0		;  (0=p-->K+lamda0,1=p-->K+sigma0,2=n-->K+sigma)
13                                              ;  (+10 for bound hypernucleus in final state)
14                doing_pion = 0		;  (ONE = TRUE)
15                doing_decay = 0		;  1=decay ON, 0=decay OFF.
16                ctau = 371.3			;  decay length (cm). 780.4 for pi, 371.3 for K
17                extra_dbase_file='extra_KAON' ;  additional input parameters.
18              end parm experiment
19              
20              begin parm kinematics_main
21              ; kinematics for hydrogen november
22 gaskelld 1.1   Ebeam = 3244.8                 ;  (MeV)
23                dEbeam = 0.05                 ;  beam energy variation (%)
24                electron_arm = 1              ;  1=hms,2=sos,3=hrsr,4=hrsl
25                hadron_arm = 2                ;  1=hms,2=sos,3=hrsr,4=hrsl
26                spec.e.P = 1571.272              ;  e arm central momentum (MeV/c)
27                spec.e.theta = 14.92          ;  e arm angle setting (degrees)
28                spec.p.P = 1292.58              ;  p arm central momentum (MeV/c)
29                spec.p.theta = 13.399           ;  p arm angle setting (degrees)
30              end parm kinematics_main
31              
32              begin parm target
33                targ.A = 1.			;  target A
34                targ.Z = 1.			;  target Z
35                targ.mass_amu = 1.007276	;  target mass in amu
36                targ.mrec_amu = 0.		;  recoil mass in amu (eep=A-1 system,pion=A-2)
37                targ.rho = 0.07332		;  target density (g/cm^3)
38                targ.thick = 293.28		;  target thick (mg/cm^2) 4.0cm long
39                targ.angle = 0.		;  target angle (for solid target) (degrees)
40                targ.abundancy = 100.		;  target purity (%)
41                targ.can = 2			;  1=beer can (fpi), 2=pudding can (nucpi)
42              end parm target
43 gaskelld 1.1 
44              begin parm debug		;  (ONES give helpful debug info)
45                debug(1) = 0			;  turns on output from brem.f
46                debug(2) = 0			;  into/outa subs
47                debug(3) = 0			;  spit out values (init. and main loop).
48                debug(4) = 0			;  mostly comp_ev, gen_rad diagnostics.
49                debug(5) = 0			;  a bit of everything
50              end parm debug
51              
52              begin parm e_arm_accept
53                SPedge.e.delta.min = -14.0    ;  delta min (%) (SPECTROMETER ACCEPTANCE!)
54                SPedge.e.delta.max =  14.0    ;  delta max (%)
55                SPedge.e.yptar.min = -60.0    ;  yptar min = (mrad)
56                SPedge.e.yptar.max =  60.0    ;  yptar max = (mrad)
57                SPedge.e.xptar.min = -100.0   ;  xptar min = (mrad)
58                SPedge.e.xptar.max =  100.0   ;  xptar max = (mrad)
59              end parm e_arm_accept
60              
61              begin parm p_arm_accept
62                SPedge.p.delta.min = -22.0    ;  delta min (%) (SPECTROMETER ACCEPTANCE!)
63                SPedge.p.delta.max =  22.0    ;  delta max (%)
64 gaskelld 1.1   SPedge.p.yptar.min = -90.0    ;  yptar min = (mrad)
65                SPedge.p.yptar.max =  90.0    ;  yptar max = (mrad)
66                SPedge.p.xptar.min = -50.0    ;  xptar min = (mrad)
67                SPedge.p.xptar.max =  50.0    ;  xptar max = (mrad)
68              end parm p_arm_accept