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
|