Macro momres mass=0.14
| kmikhail@jlab.org
| exe ~/.pawlogon.kumac
| updated by elton for Hall D
h/del *; v/del *
|
h/fil 1 momres.hbook
|
for/fil 66 'momres_m'//[mass]//'.ps'
* meta 66 -100212
*meta 66 -100112
meta 66 -111
|
opt date
option nbox
Set date 3



ZONE 1 2

vector/create m(1) r [mass]
*vector/create m(1) r 0.140
mess 'mass=' $eval(m(1))
set gsiz 0.35
set ygti 1.0
title 'Hall D CDC at 90 degrees; sigx=150mu, scat, st, 23 layers, B=1.74T'

clr
* opt liny
opt logy
igset txci 4

* opt linx;null 0.001 4.9 1e-6 0.01
opt linx;null 0.001 4.9 1e-6 3.
|opt linx;null 0.003 4.5 1e-6 3.
atit 'Distance (m)' 'Thickness/X?0!'
set hcol 2
set hwid 5
|B-field x-2,y-3,z-4
*h/pl 2 s
*h/pl 3 s
h/pl 4 s
set hcol 4
h/pl 1 s
|igset txci 2;itx  2. 0.02 'Bxdl=2.37 Txm' |5deg
|igset txci 2;itx  2. 0.02 'Bxdl=2.75 Txm' |10deg
|igset txci 2;itx  2. 0.02 'Bxdl=1.97 Txm' |15deg
|igset txci 2;itx  2. 0.02 'Bxdl=1.49 Txm' |20deg
|igset txci 2;itx  2. 0.02 'Bxdl=1.17 Txm' |25deg
* igset txci 2;itx  2. 0.02 'Bxdl=0.91 Txm' |30deg
|igset txci 2;itx  2. 0.02 'Bxdl=0.67 Txm' |35deg
|igset txci 2;itx  2. 0.02 'Bxdl=0.45 Txm' |40deg
|igset txci 2;itx  2. 0.02 'Bxdl=0.27 Txm' |45deg

igset txci 4
* opt linx;null 0.001 .7 1e-6 0.01
opt linx;null 0.001 .7 1e-6 3.
atit 'Distance (m)' 'Thickness/X?0!'
set hcol 4
h/pl 1 s
*for Hall D there is no field free region
* igset txci 2;itx  .25 0.02 'No Field region'

WAIT
|
OPT ndat
opt liny; opt linx
*opt logy
set xwin 2.3; set xlab 1.5
size 20 20
zon 2 2
set hwid 5
opt utit
red=2;gre=3;blu=4
|null 0.9 5.1 0. 0.01
PS=15|7
|PS=1
* max 10 0.02|0.006|
max 10 0.04 |0.006|
set hcol [blu]; h/pl 10 
set hcol [gre]; h/pl 11 s
set hcol [red]; h/pl 12 s
set hcol [blu]; h/pl 10 s
atit 'P (GeV/c)' '[D]P/P'|*100 ("Y#)'
|v/cr ppar(3) r 3*1;His/fit 10(1:100) pfit.f sqo 3 ppar
v/cr ppar(3) r 3*0.001;His/fit 10([PS]:100) pfit.f sqo 2 ppar
**|set fcol 6
**|FUN1 99 'cdrp.f' 100 0. 5. s
**|set fcol 1
 
opt liny; opt linx
|null 0.9 5.1 0. 0.006
* max 20 0.005
max 20 0.01
I1=6
set hcol [blu]; h/pl 20([I1]:) 
set hcol [gre]; h/pl 21 s
set hcol [red]; h/pl 22 s
set hcol [blu]; h/pl 20 s
atit 'P (GeV/c)' '[D][Q] (radians)'
v/cr tpar(2) r 2*0.01;His/fit 20([I1]:100) tfit.f sqo 2 tpar

|null 0.9 5.1 0. 0.0003
* max 30 500
max 30 1000
set hcol [blu]; h/pl 30([I1]:) 
set hcol [gre]; h/pl 31 s
set hcol [red]; h/pl 32 s
set hcol [blu]; h/pl 30 s
atit 'P (GeV/c)' '[D]X ([m]m)'
v/cr xpar(2) r 2*0.01;His/fit 30([I1]:100) tfit.f sqo 2 xpar

null 0 10 0 10 AB
igset txci [red];itx  -2. -2 'red   - mult.'
igset txci [gre];itx  2. -2 'green - pos.'
igset txci [blu];itx  6. -2 'blue  - sum'
igset txci 1;itx 6 -1 'mass=' $eval(m(1))

igset txci 1;
xs1=xpar(1);xs2=xpar(2)
itx  -2. 2.5 '[D]X=sqrt([s]^2!&?1!+([s]?2!/P[b])^2!)'
itx  0. 1. [s]?1!=[xs1]([m]m) 
itx  0. 0. [s]?2!=[xs2]([m]m.GeV) 
v/pri ppar;v/pri tpar;v/pri xpar;
ps1=ppar(1);ps2=ppar(2)
itx -2. 10.5 '[D]P/P=sqrt([s]^2!&?1!P^2!+([s]?2!/[b])^2!)'
itx  0. 9. [s]?1!=[ps1](GeV^-1!)
itx  0. 8. [s]?2!=[ps2] 
ts1=tpar(1);ts2=tpar(2)
itx -2. 6.5 '[D][Q]=sqrt([s]^2!&?1!+([s]?2!/P[b])^2!)'
itx  0. 5. [s]?1!=[ts1](rad) 
itx  0. 4. [s]?2!=[ts2](rad.GeV)
*
*
close 66