1 gaskelld 1.1 ; This is a CTP file
2
3 begin parm experiment
4 ngen = 10000 ; 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 end parm experiment
14
15 begin parm kinematics_main
16 Ebeam = 5750. ; (MeV)
17 dEbeam = 0.05 ; beam energy variation (%)
18 electron_arm = 2 ; 1=hms,2=sos,3=hrsr,4=hrsl
19 hadron_arm = 1 ; 1=hms,2=sos,3=hrsr,4=hrsl
20 spec.e.P = 1410. ; e arm central momentum (MeV/c)
21 spec.e.theta = 57.4 ; e arm angle setting (degrees)
22 gaskelld 1.1 spec.p.P = 3039. ; p arm central momentum (MeV/c)
23 spec.p.theta = 13.4 ; p arm angle setting (degrees)
24 end parm kinematics_main
25
26 begin parm target
27 targ.A = 1. ; target A
28 targ.Z = 1. ; target Z
29 targ.mass_amu = 1.007276 ; target mass in amu
30 targ.mrec_amu = 0. ; recoil mass in amu (eep=A-1 system,pion=A-2)
31 targ.rho = 0.07332 ; target density (g/cm^3)
32 targ.thick = 295.179 ; target thick (mg/cm^2)
33 targ.angle = 0. ; target angle (for solid target) (degrees)
34 targ.abundancy = 100. ; target purity (%)
35 targ.can = 2 ; 1=beer can (fpi), 2=pudding can (nucpi)
36 end parm target
37
38 begin parm debug ; (ONES give helpful debug info)
39 debug(1) = 0 ; turns on output from brem.f
40 debug(2) = 0 ; into/outa subs.
41 debug(3) = 0 ; spit out values (init. and main loop).
42 debug(4) = 0 ; mostly comp_ev, gen_rad diagnostics.
43 gaskelld 1.1 debug(5) = 0 ; a bit of everything.
44 end parm debug
45
46 begin parm e_arm_accept
47 SPedge.e.delta.min = -18.0 ; delta min (SPECTROMETER ACCEPTANCE!)
48 SPedge.e.delta.max = 18.0 ; delta max
49 SPedge.e.yptar.min = -90.0 ; .yptar.min = {TF} / 1000 (mrad)
50 SPedge.e.yptar.max = 90.0 ; .yptar.max = {TF} / 1000
51 SPedge.e.xptar.min = -50.0 ; .xptar.min = {TF} / 1000 (mrad)
52 SPedge.e.xptar.max = 50.0 ; .xptar.max = {TF} / 1000
53 end parm e_arm_accept
54
55 begin parm p_arm_accept
56 SPedge.p.delta.min = -12.0 ; delta min (SPECTROMETER ACCEPTANCE!)
57 SPedge.p.delta.max = 12.0 ; delta max
58 SPedge.p.yptar.min = -50.0 ; .yptar.min = {TF} / 1000 (mrad)
59 SPedge.p.yptar.max = 50.0 ; .yptar.max = {TF} / 1000
60 SPedge.p.xptar.min = -100.0 ; .xptar.min = {TF} / 1000 (mrad)
61 SPedge.p.xptar.max = 100.0 ; .xptar.max = {TF} / 1000
62 end parm p_arm_accept
63
64 gaskelld 1.1 ; This is a CTP file
65
66 begin parm beamandtargetinfo
67 gen.xwid = 0.008868 ; beam width - one sigma (cm) (89microns)
68 gen.ywid = 0.004235 ; beam width - one sigma (cm) (42microns)
69 targ.fr_pattern = 2. ; raster pattern: 1=square, 2=circular
70 targ.fr1 = 0.015 ; horizontal size OR inner radius(2)
71 targ.fr2 = 0.15 ; vertical size OR outer radius(2)
72 targ.xoffset = 0. ; target x-offset (cm): +x = beam right
73 targ.yoffset = 0. ; target y-offset (cm): +y = up
74 targ.zoffset = 0. ; target z-offset (cm): +z = downstream
75 ; zreal = znominal + zoffset
76 end parm beamandtergetinfo
77
78 ;These are offsets applied before the call to the single arm montecarlos.
79 begin parm spect_offset
80 spec.e.offset.x = 0. ; x offset (cm)
81 spec.e.offset.y = 0. ; y offset (cm)
82 spec.e.offset.z = 0. ; z offset (cm)
83 spec.e.offset.xptar = 0. ; xptar offset (mr) !x(y)ptar is slope, so
84 spec.e.offset.yptar = 0. ; yptar offset (mr) !it's really unitless.
85 gaskelld 1.1 spec.p.offset.x = 0. ; x offset (cm)
86 spec.p.offset.y = 0. ; y offset (cm)
87 spec.p.offset.z = 0. ; z offset (cm)
88 spec.p.offset.xptar = 2.62 ; xptar offset (mr)
89 spec.p.offset.yptar = 0. ; yptar offset (mr)
90 end parm spect_offset
91
92 begin parm simulate
93 hard_cuts = 0 ; (ONE = TRUE) SPedge and Em.max are hard cuts(ntuple)
94 using_rad = 1 ; (ONE = TRUE)
95 use_expon = 0 ; (LEAVE AT 0)
96 one_tail = 0 ; 0=all, 1=e, 2=e', 3=p, -3=all but p
97 intcor_mode = 1 ; (LEAVE AT 1)
98 spect_mode = 0 ; 0=e+p arms, -1=p arm, -2=e arm only, 1=none
99 cuts.Em.min = 0. ; (Em.min=Em.max=0.0 gives wide open cuts)
100 cuts.Em.max = 0. ; Must be wider than cuts in analysis(elastic or e,e'p)
101 using_Eloss = 1 ; (ONE = TRUE)
102 correct_Eloss = 1 ; ONE = correct reconstructed events for eloss.
103 correct_raster = 1 ; ONE = Reconstruct events using 'raster' matrix elements.
104 mc_smear = 1 ; ONE = target & hut mult scatt AND DC smearing.
105 deForest_flag = 0 ; 0=sigcc1, 1=sigcc2, -1=sigcc1 ONSHELL
106 gaskelld 1.1 rad_flag = 0 ; (radiative option #1...see init.f)
107 extrad_flag = 2 ; (rad. option #2...see init.f)
108 lambda(1) = 0.0 ; if rad_flag.eq.4 then lambda(1) = {TF}
109 lambda(2) = 0.0 ; if rad_flag.eq.4 then lambda(2) = {TF}
110 lambda(3) = 0.0 ; if rad_flag.eq.4 then lambda(3) = {TF}
111 Nntu = 1 ; ONE = generate ntuples
112 using_Coulomb = 1 ; (ONE = TRUE)
113 dE_edge_test = 0. ; (move around energy edges)
114 use_offshell_rad = 1 ; (ONE = TRUE)
115 Egamma_gen_max = 0. ; Set >0 to hardwire the Egamma limits.
116 end parm simulate
|